Faculty Bibliography

Single-cell RNA sequencing (sc-RNASeq) is a recently developed technique used to evaluate the transcriptome of individual cells. As opposed to conventional RNASeq in which entire populations are sequenced in bulk, sc-RNASeq can be beneficial when trying to better understand gene expression patterns in markedly heterogeneous populations of cells or when trying to identify transcriptional signatures of rare cells that may be underrepresented when using conventional bulk RNASeq. In this method, we describe the generation and analysis of cDNA libraries from single patient-derived glioblastoma cells using the C1 Fluidigm system. The protocol details the use of the C1 integrated fluidics circuit (IFC) for capturing, imaging and lysing cells; performing reverse transcription; and generating cDNA libraries that are ready for sequencing and analysis.

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.

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.

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.

Aflatoxin B1 (AFB1) is a mutagen and IARC (International Agency for Research on Cancer) Group 1 carcinogen that causes hepatocellular carcinoma (HCC). Here, we present the first whole-genome data on the mutational signatures of AFB1 exposure from a total of >40,000 mutations in four experimental systems: two different human cell lines, in liver tumors in wild-type mice, and in mice that carried a hepatitis B surface antigen transgene-this to model the multiplicative effects of aflatoxin exposure and hepatitis B in causing HCC. AFB1 mutational signatures from all four experimental systems were remarkably similar. We integrated the experimental mutational signatures with data from newly sequenced HCCs from Qidong County, China, a region of well-studied aflatoxin exposure. This indicated that COSMIC mutational signature 24, previously hypothesized to stem from aflatoxin exposure, indeed likely represents AFB1 exposure, possibly combined with other exposures. Among published somatic mutation data, we found evidence of AFB1 exposure in 0.7% of HCCs treated in North America, 1% of HCCs from Japan, but 16% of HCCs from Hong Kong. Thus, aflatoxin exposure apparently remains a substantial public health issue in some areas. This aspect of our study exemplifies the promise of future widespread resequencing of tumor genomes in providing new insights into the contribution of mutagenic exposures to cancer incidence.

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.

Plakophilin-2 (PKP2) is a component of the desmosome and known for its role in cell-cell adhesion. Mutations in human PKP2 associate with a life-threatening arrhythmogenic cardiomyopathy, often of right ventricular predominance. Here, we use a range of state-of-the-art methods and a cardiomyocyte-specific, tamoxifen-activated, PKP2 knockout mouse to demonstrate that in addition to its role in cell adhesion, PKP2 is necessary to maintain transcription of genes that control intracellular calcium cycling. Lack of PKP2 reduces expression of Ryr2 (coding for Ryanodine Receptor 2), Ank2 (coding for Ankyrin-B), Cacna1c (coding for CaV1.2) and Trdn (coding for triadin), and protein levels of calsequestrin-2 (Casq2). These factors combined lead to disruption of intracellular calcium homeostasis and isoproterenol-induced arrhythmias that are prevented by flecainide treatment. We propose a previously unrecognized arrhythmogenic mechanism related to PKP2 expression and suggest that mutations in PKP2 in humans may cause life-threatening arrhythmias even in the absence of structural disease.It is believed that mutations in desmosomal adhesion complex protein plakophilin 2 (PKP2) cause arrhythmia due to loss of cell-cell communication. Here the authors show that PKP2 controls the expression of proteins involved in calcium cycling in adult mouse hearts, and that lack of PKP2 can cause arrhythmia in a structurally normal heart.