Nicholas S. Davis

Transposase mediated construction of RNA-seq libraries

RNA-seq has been widely adopted as a gene-expression measurement tool due to the detail, resolution, and sensitivity of transcript characterization that the technique provides. Here we present two transposon-based methods that efficiently construct high-quality RNA-seq libraries. We first describe a method that creates RNA-seq libraries for Illumina sequencing from double-stranded cDNA with only two enzymatic reactions. We generated high-quality RNA-seq libraries from as little as 10 pg of mRNA (∼1 ng of total RNA) with this approach. We also present a strand-specific RNA-seq library construction protocol that combines transposon-based library construction with uracil DNA glycosylase and endonuclease VIII to specifically degrade the second strand constructed during cDNA synthesis. The directional RNA-seq libraries maintain the same quality as the nondirectional libraries, while showing a high degree of strand specificity, such that 99.5% of reads map to the expected genomic strand. Each transposon-based library construction method performed well when compared with standard RNA-seq library construction methods with regard to complexity of the libraries, correlation between biological replicates, and the percentage of reads that align to the genome as well as exons. Our results show that high-quality RNA-seq libraries can be constructed efficiently and in an automatable fashion using transposition technology.

Recurrent read-through fusion transcripts in breast cancer

Read-through fusion transcripts that result from the splicing of two adjacent genes in the same coding orientation are a recently discovered type of chimeric RNA. We sought to determine if read-through fusion transcripts exist in breast cancer. We performed paired-end RNA-seq of 168 breast samples, including 28 breast cancer cell lines, 42 triple negative breast cancer primary tumors, 42 estrogen receptor positive (ER+) breast cancer primary tumors, and 56 non-malignant breast tissue samples. We analyzed the sequencing data to identify breast cancer associated read-through fusion transcripts. We discovered two recurrent read-through fusion transcripts that were identified in breast cancer cell lines, confirmed across breast cancer primary tumors, and were not detected in normal tissues (SCNN1A-TNFRSF1A and CTSD-IFITM10). Both fusion transcripts use canonical splice sites to join the last splice donor of the 5′ gene to the first splice acceptor of the 3′ gene, creating an in-frame fusion transcript. Western blots indicated that the fusion transcripts are translated into fusion proteins in breast cancer cells. Custom small interfering RNAs targeting the CTSD-IFITM10 fusion junction reduced expression of the fusion transcript and reduced breast cancer cell proliferation. Read-through fusion transcripts between adjacent genes with different biochemical functions represent a new type of recurrent molecular defect in breast cancer that warrant further investigation as potential biomarkers and therapeutic targets. Both breast cancer associated fusion transcripts identified in this study involve membrane proteins (SCNN1A-TNFRSF1A and CTSD-IFITM10), which raises the possibility that they could be breast cancer-specific cell surface markers.

Genetic and Functional Drivers of Diffuse Large B Cell Lymphoma

Diffuse large B cell lymphoma (DLBCL) is the most common form of blood cancer and is characterized by a striking degree of genetic and clinical heterogeneity. This heterogeneity poses a major barrier to understanding the genetic basis of the disease and its response to therapy. Here, we performed an integrative analysis of whole-exome sequencing and transcriptome sequencing in a cohort of 1,001 DLBCL patients to comprehensively define the landscape of 150 genetic drivers of the disease. We characterized the functional impact of these genes using an unbiased CRISPR screen of DLBCL cell lines to define oncogenes that promote cell growth. A prognostic model comprising these genetic alterations outperformed current established methods: cell of origin, the International Prognostic Index comprising clinical variables, and dual MYC and BCL2 expression. These results comprehensively define the genetic drivers and their functional roles in DLBCL to identify new therapeutic opportunities in the disease.

The genetic basis of hepatosplenic T-cell lymphoma

Hepatosplenic T-cell lymphoma (HSTL) is a rare and lethal lymphoma; the genetic drivers of this disease are unknown. Through whole-exome sequencing of 68 HSTLs, we define recurrently mutated driver genes and copy-number alterations in the disease. Chromatin-modifying genes, including SETD2, INO80, and ARID1B, were commonly mutated in HSTL, affecting 62% of cases. HSTLs manifest frequent mutations in STAT5B (31%), STAT3 (9%), and PIK3CD (9%), for which there currently exist potential targeted therapies. In addition, we noted less frequent events in EZH2, KRAS, and TP53SETD2 was the most frequently silenced gene in HSTL. We experimentally demonstrated that SETD2 acts as a tumor suppressor gene. In addition, we found that mutations in STAT5B and PIK3CD activate critical signaling pathways important to cell survival in HSTL. Our work thus defines the genetic landscape of HSTL and implicates gene mutations linked to HSTL pathogenesis and potential treatment targets.Significance: We report the first systematic application of whole-exome sequencing to define the genetic basis of HSTL, a rare but lethal disease. Our work defines SETD2 as a tumor suppressor gene in HSTL and implicates genes including INO80 and PIK3CD in the disease. Cancer Discov; 7(4); 369-79. ©2017 AACR.See related commentary by Yoshida and Weinstock, p. 352This article is highlighted in the In This Issue feature, p. 339.

Enteropathy-associated T cell lymphoma subtypes are characterized by loss of function of SETD2

Enteropathy-associated T cell lymphoma (EATL) is a lethal, and the most common, neoplastic complication of celiac disease. Here, we defined the genetic landscape of EATL through whole-exome sequencing of 69 EATL tumors. SETD2 was the most frequently silenced gene in EATL (32% of cases). The JAK-STAT pathway was the most frequently mutated pathway, with frequent mutations in STAT5B as well as JAK1, JAK3, STAT3, and SOCS1. We also identified mutations in KRAS, TP53, and TERT. Type I EATL and type II EATL (monomorphic epitheliotropic intestinal T cell lymphoma) had highly overlapping genetic alterations indicating shared mechanisms underlying their pathogenesis. We modeled the effects of SETD2 loss in vivo by developing a T cell–specific knockout mouse. These mice manifested an expansion of &ggr;&dgr; T cells, indicating novel roles for SETD2 in T cell development and lymphomagenesis. Our data render the most comprehensive genetic portrait yet of this uncommon but lethal disease and may inform future classification schemes.

Genome-wide DNA methylation measurements in prostate tissues uncovers novel prostate cancer diagnostic biomarkers and transcription factor binding patterns

BackgroundCurrent diagnostic tools for prostate cancer lack specificity and sensitivity for detecting very early lesions. DNA methylation is a stable genomic modification that is detectable in peripheral patient fluids such as urine and blood plasma that could serve as a non-invasive diagnostic biomarker for prostate cancer.MethodsWe measured genome-wide DNA methylation patterns in 73 clinically annotated fresh-frozen prostate cancers and 63 benign-adjacent prostate tissues using the Illumina Infinium HumanMethylation450 BeadChip array. We overlaid the most significantly differentially methylated sites in the genome with transcription factor binding sites measured by the Encyclopedia of DNA Elements consortium. We used logistic regression and receiver operating characteristic curves to assess the performance of candidate diagnostic models.ResultsWe identified methylation patterns that have a high predictive power for distinguishing malignant prostate tissue from benign-adjacent prostate tissue, and these methylation signatures were validated using data from The Cancer Genome Atlas Project. Furthermore, by overlaying ENCODE transcription factor binding data, we observed an enrichment of enhancer of zeste homolog 2 binding in gene regulatory regions with higher DNA methylation in malignant prostate tissues.ConclusionsDNA methylation patterns are greatly altered in prostate cancer tissue in comparison to benign-adjacent tissue. We have discovered patterns of DNA methylation marks that can distinguish prostate cancers with high specificity and sensitivity in multiple patient tissue cohorts, and we have identified transcription factors binding in these differentially methylated regions that may play important roles in prostate cancer development.

RNA sequencing of pancreatic adenocarcinoma tumors yields novel expression patterns associated with long‐term survival and reveals a role for ANGPTL4

Pancreatic adenocarcinoma patients have low survival rates due to late‐stage diagnosis and high rates of cancer recurrence even after surgical resection. It is important to understand the molecular characteristics associated with survival differences in pancreatic adenocarcinoma tumors that may inform patient care.

GNA13 loss in germinal center B cells leads to impaired apoptosis and promotes lymphoma in vivo

GNA13 is the most frequently mutated gene in germinal center (GC)-derived B-cell lymphomas, including nearly a quarter of Burkitt lymphoma and GC-derived diffuse large B-cell lymphoma. These mutations occur in a pattern consistent with loss of function. We have modeled the GNA13-deficient state exclusively in GC B cells by crossing the Gna13 conditional knockout mouse strain with the GC-specific AID-Cre transgenic strain. AID-Cre(+) GNA13-deficient mice demonstrate disordered GC architecture and dark zone/light zone distribution in vivo, and demonstrate altered migration behavior, decreased levels of filamentous actin, and attenuated RhoA activity in vitro. We also found that GNA13-deficient mice have increased numbers of GC B cells that display impaired caspase-mediated cell death and increased frequency of somatic hypermutation in the immunoglobulin VH locus. Lastly, GNA13 deficiency, combined with conditional MYC transgene expression in mouse GC B cells, promotes lymphomagenesis. Thus, GNA13 loss is associated with GC B-cell persistence, in which impaired apoptosis and ongoing somatic hypermutation may lead to an increased risk of lymphoma development.

Post-mortem molecular profiling of three psychiatric disorders reveals widespread dysregulation of cell-type associated transcripts and refined disease-related transcription changes

Background Psychiatric disorders are multigenic diseases with complex etiology contributing significantly to human morbidity and mortality. Although clinically distinct, several disorders share many symptoms suggesting common underlying molecular changes exist that may implicate important regulators of pathogenesis and new therapeutic targets. Results We compared molecular signatures across brain regions and disorders in the transcriptomes of postmortem human brain samples. We performed RNA sequencing on tissue from the anterior cingulate cortex, dorsolateral prefrontal cortex, and nucleus accumbens from three groups of 24 patients each diagnosed with schizophrenia, bipolar disorder, or major depressive disorder, and from 24 control subjects, and validated the results in an independent cohort. The most significant disease differences were in the anterior cingulate cortex of schizophrenia samples compared to controls. Transcriptional changes were assessed in an independent cohort, revealing the transcription factor EGR1 as significantly down regulated in both cohorts and as a potential regulator of broader transcription changes observed in schizophrenia patients. Additionally, broad down regulation of genes specific to neurons and concordant up regulation of genes specific to astrocytes was observed in SZ and BPD patients relative to controls. We also assessed the biochemical consequences of gene expression changes with untargeted metabolomic profiling and identified disruption of GABA levels in schizophrenia patients. Conclusions We provide a comprehensive post-mortem transcriptome profile of three psychiatric disorders across three brain regions. We highlight a high-confidence set of independently validated genes differentially expressed between schizophrenia and control patients in the anterior cingulate cortex and integrate transcriptional changes with untargeted metabolite profiling.