Faculty Bibliography

We report here an update to the draft genome sequence of Kluyvera intestini sp. nov. strain GT-16, generated using MinION long-read sequencing technology. The complete genome sequence of the human-derived strain GT-16 measured 5,768,848 bp. An improved high-quality complete genome sequence provides insights into the mobility potential of resistance genes in this species.

Exposure to aristolochic acid (AA, IARC Group 1 carcinogen) present in some traditional herbal medicines leads to aristolochic acid nephropathy (AAN), often complicated by development of multiple urothelial carcinomas of sequential onset. We used genome-scale mutational signature analysis of multiple urinary tract tumors of AAN cases from a unique patient group to determine the relationships of the patients' late-onset second cancers to the AA exposure as well as to the first cancers. Aristolactam-DNA adduct-positive AAN patients (n=4) who developed cancer within 8 years following the initial exposure to AA were chosen for analysis of their first cancers (upper tract urothelial carcinomas, UTUC) and second cancers of delayed onset (1-9 years after first-cancer diagnosis, involving the bladder or ureteral meatus). All patients had received a kidney transplant before developing the second cancers and had a functional renal graft prior to prophylactic nephroureterectomy. Genomic DNAs were isolated from FFPE sections of the renal cortex and from the upper and lower tract tumors of each patient using laser capture micro-dissection or macrodissection of the tumor areas. Low-coverage (~15x) exome 100-bp paired-end sequencing was performed using Illumina HiSeq2500. Customized variant calling identified somatic variants absent in non-tumor tissues. Non-negative matrix factorization was applied to extract mutational signatures in the tumor tissues. In all cases, we established the mutational signature of AA (the COSMIC signature 22) in the first UTUC as well as second cancers involving the bladder or lower ureter (meatus). Additionally, the first and second cancers harbored considerable overlaps in exposure-specific (A>T) somatic mutations. This finding provides evidence that the delayed onset of bladder urothelial carcinomas in AAN patients is likely due to distal seeding of cancer cells originating from the primary UTUC tumors. Our first-of-its-kind study addresses the risk as well as mechanistic factors leading to the second, lateonset bladder urothelial carcinomas following kidney transplantation and primary UTUC development. Our results underline the importance of long-term bladder follow-up in high-risk populations with established or suspected AA exposure

Atlantic tomcod in the Hudson River Estuary bioaccumulate high hepatic burdens of environmental toxicants. Previously, we demonstrated that Hudson River tomcod developed resistance to TCDD and PCB toxicity probably through strong natural selection during their early life-stages for a variant of the Aryl Hydrocarbon Receptor2 (AHR2). Here, we evaluated the genomic consequences of the resistant genotype by comparing global gene expression in larval tomcod from the Hudson River with expression in larvae from a nearby sensitive population (Shinnecock Bay). We developed an annotated draft tomcod genome to explore the effects of multigenerational exposure to toxicants and a functionally impaired AHR2 on the transcriptome. We used the tomcod genome as a reference in RNA-Seq to compare global gene expression in tomcod larvae from the Hudson River and Shinnecock Bay after experimental exposure of larvae to graded doses of TCDD. We found dramatic differences between offspring from the two populations in the number of genes that were differentially expressed at all doses (0.01, 0.1, and 1 ppb) and even in the vehicle controls. At the two lowest TCDD doses, 250 and 1,141 genes were differentially expressed in Shinnecock Bay larvae compared with 14 and 12, respectively, in Hudson River larvae. At the highest dose (1.0 ppb), 934 genes were differentially expressed in Shinnecock Bay larvae and 173 in Hudson River larvae, but only 28 (16%) of affected genes were shared among both populations. Given the large difference between the two populations in the number and identity of differentially expressed genes, it is likely that the AHR2 pathway interacts directly or indirectly with many genes beyond those known in the AHR2 battery and that other regulatory systems may also respond to TCDD exposure. The effects of chronic multi-generational exposure to environmental toxicants on the genome of Hudson River tomcod are much greater than previously expected.

Treating KRAS-mutant lung adenocarcinoma (LUAD) remains a major challenge in cancer treatment given the difficulties associated with directly inhibiting the KRAS oncoprotein. One approach to addressing this challenge is to define mutations that frequently co-occur with those in KRAS, which themselves may lead to therapeutic vulnerabilities in tumors. Approximately 20% of KRAS-mutant LUAD tumors carry loss-of-function mutations in the KEAP1 gene encoding Kelch-like ECH-associated protein 1 (refs. 2, 3, 4), a negative regulator of nuclear factor erythroid 2-like 2 (NFE2L2; hereafter NRF2), which is the master transcriptional regulator of the endogenous antioxidant response. The high frequency of mutations in KEAP1 suggests an important role for the oxidative stress response in lung tumorigenesis. Using a CRISPR-Cas9-based approach in a mouse model of KRAS-driven LUAD, we examined the effects of Keap1 loss in lung cancer progression. We show that loss of Keap1 hyperactivates NRF2 and promotes KRAS-driven LUAD in mice. Through a combination of CRISPR-Cas9-based genetic screening and metabolomic analyses, we show that Keap1- or Nrf2-mutant cancers are dependent on increased glutaminolysis, and this property can be therapeutically exploited through the pharmacological inhibition of glutaminase. Finally, we provide a rationale for stratification of human patients with lung cancer harboring KRAS/KEAP1- or KRAS/NRF2-mutant lung tumors as likely to respond to glutaminase inhibition.

OBJECTIVE: Atherosclerosis is exaggerated in African American (AA) systemic lupus erythematosus (SLE) patients, with doubled cardiovascular disease (CVD) risk compared to White patients. The extent to which common Apolipoprotein L1 (APOL1) risk alleles (RA) contribute to this trend is unknown. This retrospective cohort study assessed prevalent atherosclerotic disease across APOL1 genotypes in AA SLE patients. METHODS: One hundred thirteen AA SLE subjects were APOL1-genotyped and stratified as having: zero risk alleles, one risk allele, or two risk alleles. Chart review assessed CVD manifestations including abdominal aortic aneurysm, angina, carotid artery disease, coronary artery disease, myocardial infarction, peripheral vascular disease, stroke, and vascular calcifications. Associations between the genotypes and a composite endpoint defined as one or more CVD manifestations were calculated using logistic regression. Symptomatic atherosclerotic disease, excluding incidental vascular calcifications, was also assessed. RESULTS: The 0-risk-allele, 1-risk-allele and 2-risk-allele groups, respectively, comprised 34%, 53%, and 13% of the cohort. Respectively, 13.2%, 41.7%, and 60.0% of the 0-risk allele, 1-risk-allele, and 2-risk-allele groups met the composite endpoint of atherosclerotic CVD (p = 0.001). Adjusting for risk factors-including smoking, ESRD, BMI >25 and hypertension-we observed an association between carrying one or more RA and atherosclerotic CVD (OR = 7.1; p = 0.002). For symptomatic disease, the OR was 3.5 (p = 0.02). In a time-to-event analysis, the proportion of subjects free from the composite primary endpoint, symptomatic atherosclerotic CVD, was higher in the 0-risk-allele group compared to the 1-risk-allele and 2-risk-allele groups (chi2 = 6.5; p = 0.04). CONCLUSIONS: Taken together, the APOL1 RAs associate with prevalent atherosclerotic CVD in this cohort of AA SLE patients, perhaps reflecting a potentiating effect of SLE on APOL1-related cardiovascular phenotypes.

Increased intestinal permeability and translocation of gut bacteria trigger various polyaetiological diseases associated with chronic inflammation and underlie a variety of poorly treatable pathologies. Previous studies have established a primary role of the microbiota composition and intestinal permeability in such pathologies. Using a rat model, we examined the effects of exposure to a bacteriophage cocktail on intestinal permeability and relative abundance of taxonomic units in the gut bacterial community. There was an increase in markers of impaired gut permeability, such as the lactulose/mannitol ratio, plasma endotoxin concentrations, and serum levels of inflammation-related cytokines, following the bacteriophage challenge. We observed significant differences in the alpha diversity of faecal bacterial species and found that richness and diversity index values increased following the bacteriophage challenge. There was a reduction in the abundance of Blautia, Catenibacterium, Lactobacillus, and Faecalibacterium species and an increase in Butyrivibrio, Oscillospira and Ruminococcus after bacteriophage administration. These findings provide novel insights into the role of bacteriophages as potentially pathogenic for mammals and their possible implication in the development of diseases associated with increased intestinal permeability.