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.

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 γδ 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.

BACKGROUND: Current 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. METHODS: We 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. RESULTS: We 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. CONCLUSIONS: DNA 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.

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.

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

BACKGROUND: 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. RESULTS: RNA sequencing was performed for 51 patient tumor tissues extracted from patients undergoing surgical resection, and expression was associated with overall survival time from diagnosis. Our analysis uncovered 323 transcripts whose expression correlates with survival time in our pancreatic patient cohort. This genomic signature was validated in an independent RNA-seq dataset of 68 additional patients from the International Cancer Genome Consortium. We demonstrate that this transcriptional profile is largely independent of markers of cellular division and present a 19-transcript predictive model built from a subset of the 323 transcripts that can distinguish patients with differing survival times across both the training and validation patient cohorts. We present evidence that a subset of the survival-associated transcripts is associated with resistance to gemcitabine treatment in vitro, and reveal that reduced expression of one of the survival-associated transcripts, Angiopoietin-like 4, impairs growth of a gemcitabine-resistant pancreatic cancer cell line. CONCLUSIONS: Gene expression patterns in pancreatic adenocarcinoma tumors can distinguish patients with differing survival outcomes after undergoing surgical resection, and the survival difference could be associated with the intrinsic gemcitabine sensitivity of primary patient tumors. Thus, these transcriptional differences may impact patient care by distinguishing patients who would benefit from a non-gemcitabine based therapy.

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.

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.

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.