Barrett's esophagus is thought to progress to esophageal adenocarcinoma (EAC) through a stepwise progression with loss of CDKN2A followed by TP53 inactivation and aneuploidy. Here we present whole-exome sequencing from 25 pairs of EAC and Barrett's esophagus and from 5 patients whose Barrett's esophagus and tumor were extensively sampled. Our analysis showed that oncogene amplification typically occurred as a late event and that TP53 mutations often occurred early in Barrett's esophagus progression, including in non-dysplastic epithelium. Reanalysis of additional EAC exome data showed that the majority (62.5%) of EACs emerged following genome doubling and that tumors with genomic doubling had different patterns of genomic alterations, with more frequent oncogenic amplification and less frequent inactivation of tumor suppressors, including CDKN2A. These data suggest that many EACs emerge not through the gradual accumulation of tumor-suppressor alterations but rather through a more direct path whereby a TP53-mutant cell undergoes genome doubling, followed by the acquisition of oncogenic amplifications.

Tumor mutational burden correlates with response to immune checkpoint blockade in multiple solid tumors, although in microsatellite-stable tumors this association is of uncertain clinical utility. Here we uniformly analyzed whole-exome sequencing (WES) of 249 tumors and matched normal tissue from patients with clinically annotated outcomes to immune checkpoint therapy, including radiographic response, across multiple cancer types to examine additional tumor genomic features that contribute to selective response. Our analyses identified genomic correlates of response beyond mutational burden, including somatic events in individual driver genes, certain global mutational signatures, and specific HLA-restricted neoantigens. However, these features were often interrelated, highlighting the complexity of identifying genetic driver events that generate an immunoresponsive tumor environment. This study lays a path forward in analyzing large clinical cohorts in an integrated and multifaceted manner to enhance the ability to discover clinically meaningful predictive features of response to immune checkpoint blockade.

Current genomics methods are designed to handle tens to thousands of samples but will need to scale to millions to match the pace of data and hypothesis generation in biomedical science. Here, we show that high efficiency at low cost can be achieved by leveraging general-purpose libraries for computing using graphics processing units (GPUs), such as PyTorch and TensorFlow. We demonstrate > 200-fold decreases in runtime and ~ 5-10-fold reductions in cost relative to CPUs. We anticipate that the accessibility of these libraries will lead to a widespread adoption of GPUs in computational genomics.

Nonalcoholic steatohepatitis (NASH) is the most common chronic liver disease globally and a leading cause for liver transplantation in the US. Its pathogenesis remains imprecisely defined. We combined two high-resolution modalities to tissue samples from NASH clinical trials, machine learning (ML)-based quantification of histological features and transcriptomics, to identify genes that are associated with disease progression and clinical events. A histopathology-driven 5-gene expression signature predicted disease progression and clinical events in patients with NASH with F3 (pre-cirrhotic) and F4 (cirrhotic) fibrosis. Notably, the Notch signaling pathway and genes implicated in liver-related diseases were enriched in this expression signature. In a validation cohort where pharmacologic intervention improved disease histology, multiple Notch signaling components were suppressed.

Craniopharyngiomas are epithelial tumors that typically arise in the suprasellar region of the brain. Patients experience substantial clinical sequelae from both extension of the tumors and therapeutic interventions that damage the optic chiasm, the pituitary stalk and the hypothalamic area. Using whole-exome sequencing, we identified mutations in CTNNB1 (β-catenin) in nearly all adamantinomatous craniopharyngiomas examined (11/12, 92%) and recurrent mutations in BRAF (resulting in p.Val600Glu) in all papillary craniopharyngiomas (3/3, 100%). Targeted genotyping revealed BRAF p.Val600Glu in 95% of papillary craniopharyngiomas (36 of 39 tumors) and mutation of CTNNB1 in 96% of adamantinomatous craniopharyngiomas (51 of 53 tumors). The CTNNB1 and BRAF mutations were clonal in each tumor subtype, and we detected no other recurrent mutations or genomic aberrations in either subtype. Adamantinomatous and papillary craniopharyngiomas harbor mutations that are mutually exclusive and clonal. These findings have important implications for the diagnosis and treatment of these neoplasms.

Resistance to the Bruton's tyrosine kinase (BTK) inhibitor ibrutinib has been attributed solely to mutations in BTK and related pathway molecules. Using whole-exome and deep-targeted sequencing, we dissect evolution of ibrutinib resistance in serial samples from five chronic lymphocytic leukaemia patients. In two patients, we detect BTK-C481S mutation or multiple PLCG2 mutations. The other three patients exhibit an expansion of clones harbouring del(8p) with additional driver mutations (EP300, MLL2 and EIF2A), with one patient developing trans-differentiation into CD19-negative histiocytic sarcoma. Using droplet-microfluidic technology and growth kinetic analyses, we demonstrate the presence of ibrutinib-resistant subclones and estimate subclone size before treatment initiation. Haploinsufficiency of TRAIL-R, a consequence of del(8p), results in TRAIL insensitivity, which may contribute to ibrutinib resistance. These findings demonstrate that the ibrutinib therapy favours selection and expansion of rare subclones already present before ibrutinib treatment, and provide insight into the heterogeneity of genetic changes associated with ibrutinib resistance.

In Rspondin-based 3D cultures, Lgr5 stem cells from multiple organs form ever-expanding epithelial organoids that retain their tissue identity. We report the establishment of tumor organoid cultures from 20 consecutive colorectal carcinoma (CRC) patients. For most, organoids were also generated from adjacent normal tissue. Organoids closely recapitulate several properties of the original tumor. The spectrum of genetic changes within the "living biobank" agrees well with previous large-scale mutational analyses of CRC. Gene expression analysis indicates that the major CRC molecular subtypes are represented. Tumor organoids are amenable to high-throughput drug screens allowing detection of gene-drug associations. As an example, a single organoid culture was exquisitely sensitive to Wnt secretion (porcupine) inhibitors and carried a mutation in the negative Wnt feedback regulator RNF43, rather than in APC. Organoid technology may fill the gap between cancer genetics and patient trials, complement cell-line- and xenograft-based drug studies, and allow personalized therapy design. PAPERCLIP.

Obesity is linked to increased incidence of endometrioid endometrial cancer (EEC) and complex atypical hyperplasia (CAH). We here explore pattern and sequence of molecular alterations characterizing endometrial carcinogenesis in general and related to body mass index (BMI), to improve diagnostic stratification and treatment strategies. We performed molecular characterization of 729 prospectively collected EEC and CAH. Candidate biomarkers were identified in frozen samples by whole-exome and Sanger sequencing, oligonucleotide gene expression and Reverse Phase Protein Arrays (investigation cohort) and further explored in formalin fixed tissues by immunohistochemistry and Fluorescent in Situ Hybridization (validation cohort). We here demonstrate that PIK3CA mutations, PTEN loss, PI3K and KRAS activation are early events in endometrial carcinogenesis. Molecular changes related to KRAS activation and inflammation are more common in obese CAH patients, suggesting different prevention and systemic treatment strategies in obese and non-obese patients. We also found that oncoprotein Stathmin might improve preoperative diagnostic distinction between premalignant and malignant endometrial lesions.

Cancer genomes contain large numbers of somatic mutations but few of these mutations drive tumor development. Current approaches either identify driver genes on the basis of mutational recurrence or approximate the functional consequences of nonsynonymous mutations by using bioinformatic scores. Passenger mutations are enriched in characteristic nucleotide contexts, whereas driver mutations occur in functional positions, which are not necessarily surrounded by a particular nucleotide context. We observed that mutations in contexts that deviate from the characteristic contexts around passenger mutations provide a signal in favor of driver genes. We therefore developed a method that combines this feature with the signals traditionally used for driver-gene identification. We applied our method to whole-exome sequencing data from 11,873 tumor-normal pairs and identified 460 driver genes that clustered into 21 cancer-related pathways. Our study provides a resource of driver genes across 28 tumor types with additional driver genes identified according to mutations in unusual nucleotide contexts.

Knowledge of the genomic landscape of chronic lymphocytic leukemia (CLL) grows increasingly detailed, providing challenges in contextualizing the accumulated information. To define the underlying networks, we here perform a multi-platform molecular characterization. We identify major subgroups characterized by genomic instability (GI) or activation of epithelial-mesenchymal-transition (EMT)-like programs, which subdivide into non-inflammatory and inflammatory subtypes. GI CLL exhibit disruption of genome integrity, DNA-damage response and are associated with mutagenesis mediated through activation-induced cytidine deaminase or defective mismatch repair. TP53 wild-type and mutated/deleted cases constitute a transcriptionally uniform entity in GI CLL and show similarly poor progression-free survival at relapse. EMT-like CLL exhibit high genomic stability, reduced benefit from the addition of rituximab and EMT-like differentiation is inhibited by induction of DNA damage. This work extends the perspective on CLL biology and risk categories in TP53 wild-type CLL. Furthermore, molecular targets identified within each subgroup provide opportunities for new treatment approaches.

Pancreatic cancer currently holds the position of third deadliest cancer in the United States and the 5-year survival rate is among the lowest for major cancers at just 12%. Thus, continued research efforts to better understand the clinical and molecular underpinnings of pancreatic cancer are critical to developing both early detection methodologies as well as improved therapeutic options. This study introduces Pancreatic Cancer Action Network's (PanCAN's) SPARK, a cloud-based data and analytics platform that integrates patient health data from the PanCAN's research initiatives and aims to accelerate pancreatic cancer research by making real-world patient health data and analysis tools easier to access and use.

Malignant pleural mesothelioma (MPM) is an aggressive neoplasm associated with asbestos exposure. Although previous studies based on candidate gene approaches have identified important common somatic mutations in MPM, these studies have focused on small sets of genes and have provided a limited view of the genetic alterations underlying this disease. Here, we performed whole-exome sequencing on DNA from 22 MPMs and matched blood samples, and identified 517 somatic mutations across 490 mutated genes. Integrative analysis of mutations and somatic copy-number alterations revealed frequent genetic alterations in BAP1, NF2, CDKN2A, and CUL1. Our study presents the first unbiased view of the genomic basis of MPM.

Angiocentric gliomas are pediatric low-grade gliomas (PLGGs) without known recurrent genetic drivers. We performed genomic analysis of new and published data from 249 PLGGs, including 19 angiocentric gliomas. We identified MYB-QKI fusions as a specific and single candidate driver event in angiocentric gliomas. In vitro and in vivo functional studies show that MYB-QKI rearrangements promote tumorigenesis through three mechanisms: MYB activation by truncation, enhancer translocation driving aberrant MYB-QKI expression and hemizygous loss of the tumor suppressor QKI. To our knowledge, this represents the first example of a single driver rearrangement simultaneously transforming cells via three genetic and epigenetic mechanisms in a tumor.

The estrogen receptor alpha (ERα) is highly expressed in both endometrial and breast cancers, and represents the most prevalent therapeutic target in breast cancer. However, anti-estrogen therapy has not been shown to be effective in endometrial cancer. Recently it has been shown that hormone-binding domain alterations of ERα in breast cancer contribute to acquired resistance to anti-estrogen therapy. In analyses of genomic data from The Cancer Genome Atlas (TCGA), we observe that endometrial carcinomas manifest recurrent ESR1 gene amplifications that truncate the hormone-binding domain encoding region of ESR1 and are associated with reduced mRNA expression of exons encoding the hormone-binding domain. These findings support a role for hormone-binding alterations of ERα in primary endometrial cancer, with potentially important therapeutic implications.

Small cell lung carcinoma (SCLC) is a highly lethal, smoking-associated cancer with few known targetable genetic alterations. Using genome sequencing, we characterized the somatic evolution of a genetically engineered mouse model (GEMM) of SCLC initiated by loss of Trp53 and Rb1. We identified alterations in DNA copy number and complex genomic rearrangements and demonstrated a low somatic point mutation frequency in the absence of tobacco mutagens. Alterations targeting the tumor suppressor Pten occurred in the majority of murine SCLC studied, and engineered Pten deletion accelerated murine SCLC and abrogated loss of Chr19 in Trp53; Rb1; Pten compound mutant tumors. Finally, we found evidence for polyclonal and sequential metastatic spread of murine SCLC by comparative sequencing of families of related primary tumors and metastases. We propose a temporal model of SCLC tumorigenesis with implications for human SCLC therapeutics and the nature of cancer-genome evolution in GEMMs.

We performed harmonized molecular and clinical analysis on 1,048 melanomas and discovered markedly different global genomic properties among subtypes (BRAF, (N)RAS, NF1, triple wild-type (TWT)), subtype-specific preferences for secondary driver genes and active mutational processes previously unreported in melanoma. Secondary driver genes significantly enriched in specific subtypes reflected preferential dysregulation of additional pathways, such as induction of transforming growth factor-β signaling in BRAF melanomas and inactivation of the SWI/SNF complex in (N)RAS melanomas, and select co-mutation patterns coordinated selective response to immune checkpoint blockade. We also defined the mutational landscape of TWT melanomas and revealed enrichment of DNA-repair-defect signatures in this subtype, which were associated with transcriptional downregulation of key DNA-repair genes, and may revive previously discarded or currently unconsidered therapeutic modalities for genomically stratified melanoma patient subsets. Broadly, harmonized meta-analysis of melanoma whole exomes revealed distinct molecular drivers that may point to multiple opportunities for biological and therapeutic investigation.

Recent studies have detailed the genomic landscape of primary endometrial cancers, but the evolution of these cancers into metastases has not been characterized. We performed whole-exome sequencing of 98 tumor biopsies including complex atypical hyperplasias, primary tumors and paired abdominopelvic metastases to survey the evolutionary landscape of endometrial cancer. We expanded and reanalyzed The Cancer Genome Atlas (TCGA) data, identifying new recurrent alterations in primary tumors, including mutations in the estrogen receptor cofactor gene NRIP1 in 12% of patients. We found that likely driver events were present in both primary and metastatic tissue samples, with notable exceptions such as ARID1A mutations. Phylogenetic analyses indicated that the sampled metastases typically arose from a common ancestral subclone that was not detected in the primary tumor biopsy. These data demonstrate extensive genetic heterogeneity in endometrial cancers and relative homogeneity across metastatic sites.

Hemangioblastomas consist of 10-20% neoplastic "stromal" cells within a vascular tumor cell mass of reactive pericytes, endothelium and lymphocytes. Familial cases of central nervous system hemangioblastoma uniformly result from mutations in the Von Hippel-Lindau (VHL) gene. In contrast, inactivation of VHL has been previously observed in only a minority of sporadic hemangioblastomas, suggesting an alternative genetic etiology. We performed deep-coverage DNA sequencing on 32 sporadic hemangioblastomas (whole exome discovery cohort n = 10, validation n = 22), followed by analysis of clonality, copy number alteration, and somatic mutation. We identified somatic mutation, loss of heterozygosity and/or deletion of VHL in 8 of 10 discovery cohort tumors. VHL inactivating events were ultimately detected in 78% (25/32) of cases. No other gene was significantly mutated. Overall, deep-coverage sequence analysis techniques uncovered VHL alterations within the neoplastic fraction of these tumors at higher frequencies than previously reported. Our findings support the central role of VHL inactivation in the molecular pathogenesis of both familial and sporadic hemangioblastomas.

Comparison of sequencing data from a tumor sample with data from a matched germline control is a key step for accurate detection of somatic mutations. Detection sensitivity for somatic variants is greatly reduced when the matched normal sample is contaminated with tumor cells. To overcome this limitation, we developed deTiN, a method that estimates the tumor-in-normal (TiN) contamination level and, in cases affected by contamination, improves sensitivity by reclassifying initially discarded variants as somatic.

To further our understanding of the somatic genetic basis of uveal melanoma, we sequenced the protein-coding regions of 52 primary tumors and 3 liver metastases together with paired normal DNA. Known recurrent mutations were identified in GNAQ, GNA11, BAP1, EIF1AX, and SF3B1. The role of mutated EIF1AX was tested using loss of function approaches including viability and translational efficiency assays. Knockdown of both wild type and mutant EIF1AX was lethal to uveal melanoma cells. We probed the function of N-terminal tail EIF1AX mutations by performing RNA sequencing of polysome-associated transcripts in cells expressing endogenous wild type or mutant EIF1AX. Ribosome occupancy of the global translational apparatus was sensitive to suppression of wild type but not mutant EIF1AX. Together, these studies suggest that cells expressing mutant EIF1AX may exhibit aberrant translational regulation, which may provide clonal selective advantage in the subset of uveal melanoma that harbors this mutation.