Genomic structural variation (SV), a common hallmark of cancer, has important predictive and therapeutic implications. However, accurately detecting SV using high-throughput sequencing data remains challenging, especially for 'targeted' resequencing efforts. This is critically important in the clinical setting where targeted resequencing is frequently being applied to rapidly assess clinically actionable mutations in tumor biopsies in a cost-effective manner. We present BreaKmer, a novel approach that uses a 'kmer' strategy to assemble misaligned sequence reads for predicting insertions, deletions, inversions, tandem duplications and translocations at base-pair resolution in targeted resequencing data. Variants are predicted by realigning an assembled consensus sequence created from sequence reads that were abnormally aligned to the reference genome. Using targeted resequencing data from tumor specimens with orthogonally validated SV, non-tumor samples and whole-genome sequencing data, BreaKmer had a 97.4% overall sensitivity for known events and predicted 17 positively validated, novel variants. Relative to four publically available algorithms, BreaKmer detected SV with increased sensitivity and limited calls in non-tumor samples, key features for variant analysis of tumor specimens in both the clinical and research settings.

Nearly all prostate cancer deaths are from metastatic castration-resistant prostate cancer (mCRPC), but there have been few whole-genome sequencing (WGS) studies of this disease state. We performed linked-read WGS on 23 mCRPC biopsy specimens and analyzed cell-free DNA sequencing data from 86 patients with mCRPC. In addition to frequent rearrangements affecting known prostate cancer genes, we observed complex rearrangements of the AR locus in most cases. Unexpectedly, these rearrangements include highly recurrent tandem duplications involving an upstream enhancer of AR in 70%-87% of cases compared with <2% of primary prostate cancers. A subset of cases displayed AR or MYC enhancer duplication in the context of a genome-wide tandem duplicator phenotype associated with CDK12 inactivation. Our findings highlight the complex genomic structure of mCRPC, nominate alterations that may inform prostate cancer treatment, and suggest that additional recurrent events in the non-coding mCRPC genome remain to be discovered.

Macroautophagy (hereafter referred to as autophagy) is an evolutionarily conserved catabolic process necessary for normal recycling of cellular constituents and for appropriate response to cellular stress. Although several genes belonging to the core molecular machinery involved in autophagosome formation have been discovered, relatively little is known about the nature of signaling networks controlling autophagy upon intracellular or extracellular stimuli. We discovered ATG8-like proteins (MAP1LC3B, GABARAP and GABARAPL1) as novel interactors of MAPK15/ERK8, a MAP kinase involved in cell proliferation and transformation. Based on the role of these proteins in the autophagic process, we demonstrated that MAPK15 is indeed localized to autophagic compartments and increased, in a kinase-dependent fashion, ATG8-like proteins lipidation, autophagosome formation and SQSTM1 degradation, while decreasing LC3B inhibitory phosphorylation. Interestingly, we also identified a conserved LC3-interacting region (LIR) in MAPK15 responsible for its interaction with ATG8-like proteins, for its localization to autophagic structures and, consequently, for stimulation of the formation of these compartments. Furthermore, we reveal that MAPK15 activity was induced in response to serum and amino-acid starvation and that this stimulus, in turn, required endogenous MAPK15 expression to induce the autophagic process. Altogether, these results suggested a new function for MAPK15 as a regulator of autophagy, acting through interaction with ATG8 family proteins. Also, based on the key role of this process in several human diseases, these results supported the use of this MAP kinase as a potential novel therapeutic target.

Although a few cancer genes are mutated in a high proportion of tumours of a given type (>20%), most are mutated at intermediate frequencies (2-20%). To explore the feasibility of creating a comprehensive catalogue of cancer genes, we analysed somatic point mutations in exome sequences from 4,742 human cancers and their matched normal-tissue samples across 21 cancer types. We found that large-scale genomic analysis can identify nearly all known cancer genes in these tumour types. Our analysis also identified 33 genes that were not previously known to be significantly mutated in cancer, including genes related to proliferation, apoptosis, genome stability, chromatin regulation, immune evasion, RNA processing and protein homeostasis. Down-sampling analysis indicates that larger sample sizes will reveal many more genes mutated at clinically important frequencies. We estimate that near-saturation may be achieved with 600-5,000 samples per tumour type, depending on background mutation frequency. The results may help to guide the next stage of cancer genomics.

Cancer cells that express oncogenic alleles of RAS typically require sustained expression of the mutant allele for survival, but the molecular basis of this oncogene dependency remains incompletely understood. To identify genes that can functionally substitute for oncogenic RAS, we systematically expressed 15,294 open reading frames in a human KRAS-dependent colon cancer cell line engineered to express an inducible KRAS-specific shRNA. We found 147 genes that promoted survival upon KRAS suppression. In particular, the transcriptional coactivator YAP1 rescued cell viability in KRAS-dependent cells upon suppression of KRAS and was required for KRAS-induced cell transformation. Acquired resistance to Kras suppression in a Kras-driven murine lung cancer model also involved increased YAP1 signaling. KRAS and YAP1 converge on the transcription factor FOS and activate a transcriptional program involved in regulating the epithelial-mesenchymal transition (EMT). Together, these findings implicate transcriptional regulation of EMT by YAP1 as a significant component of oncogenic RAS signaling.

We describe the landscape of somatic genomic alterations of 66 chromophobe renal cell carcinomas (ChRCCs) on the basis of multidimensional and comprehensive characterization, including mtDNA and whole-genome sequencing. The result is consistent that ChRCC originates from the distal nephron compared with other kidney cancers with more proximal origins. Combined mtDNA and gene expression analysis implicates changes in mitochondrial function as a component of the disease biology, while suggesting alternative roles for mtDNA mutations in cancers relying on oxidative phosphorylation. Genomic rearrangements lead to recurrent structural breakpoints within TERT promoter region, which correlates with highly elevated TERT expression and manifestation of kataegis, representing a mechanism of TERT upregulation in cancer distinct from previously observed amplifications and point mutations.

The analysis of exonic DNA from prostate cancers has identified recurrently mutated genes, but the spectrum of genome-wide alterations has not been profiled extensively in this disease. We sequenced the genomes of 57 prostate tumors and matched normal tissues to characterize somatic alterations and to study how they accumulate during oncogenesis and progression. By modeling the genesis of genomic rearrangements, we identified abundant DNA translocations and deletions that arise in a highly interdependent manner. This phenomenon, which we term "chromoplexy," frequently accounts for the dysregulation of prostate cancer genes and appears to disrupt multiple cancer genes coordinately. Our modeling suggests that chromoplexy may induce considerable genomic derangement over relatively few events in prostate cancer and other neoplasms, supporting a model of punctuated cancer evolution. By characterizing the clonal hierarchy of genomic lesions in prostate tumors, we charted a path of oncogenic events along which chromoplexy may drive prostate carcinogenesis.

We describe a computational method that infers tumor purity and malignant cell ploidy directly from analysis of somatic DNA alterations. The method, named ABSOLUTE, can detect subclonal heterogeneity and somatic homozygosity, and it can calculate statistical sensitivity for detection of specific aberrations. We used ABSOLUTE to analyze exome sequencing data from 214 ovarian carcinoma tumor-normal pairs. This analysis identified both pervasive subclonal somatic point-mutations and a small subset of predominantly clonal and homozygous mutations, which were overrepresented in the tumor suppressor genes TP53 and NF1 and in a candidate tumor suppressor gene CDK12. We also used ABSOLUTE to infer absolute allelic copy-number profiles from 3,155 diverse cancer specimens, revealing that genome-doubling events are common in human cancer, likely occur in cells that are already aneuploid, and influence pathways of tumor progression (for example, with recessive inactivation of NF1 being less common after genome doubling). ABSOLUTE will facilitate the design of clinical sequencing studies and studies of cancer genome evolution and intra-tumor heterogeneity.

The Mixed-Lineage Leukemia (MLL) protein is a histone methyltransferase that is mutated in clinically and biologically distinctive subsets of acute leukemia. MLL normally associates with a cohort of highly conserved cofactors to form a macromolecular complex that includes menin, a product of the MEN1 tumor suppressor gene, which is mutated in heritable and sporadic endocrine tumors. We demonstrate here that oncogenic MLL fusion proteins retain an ability to stably associate with menin through a high-affinity, amino-terminal, conserved binding motif and that this interaction is required for the initiation of MLL-mediated leukemogenesis. Furthermore, menin is essential for maintenance of MLL-associated but not other oncogene induced myeloid transformation. Acute genetic ablation of menin reverses aberrant Hox gene expression mediated by MLL-menin promoter-associated complexes, and specifically abrogates the differentiation arrest and oncogenic properties of MLL-transformed leukemic blasts. These results demonstrate that a human oncoprotein is critically dependent on direct physical interaction with a tumor suppressor protein for its oncogenic activity, validate a potential target for molecular therapy, and suggest central roles for menin in altered epigenetic functions underlying the pathogenesis of hematopoietic cancers.

Neuroblastoma, an embryonal tumour of the peripheral sympathetic nervous system, accounts for approximately 15% of all deaths due to childhood cancer. High-risk neuroblastomas are rapidly progressive; even with intensive myeloablative chemotherapy, relapse is common and almost uniformly fatal. Here we report the detection of previously unknown mutations in the ALK gene, which encodes a receptor tyrosine kinase, in 8% of primary neuroblastomas. Five non-synonymous sequence variations were identified in the kinase domain of ALK, of which three were somatic and two were germ line. The most frequent mutation, F1174L, was also identified in three different neuroblastoma cell lines. ALK complementary DNAs encoding the F1174L and R1275Q variants, but not the wild-type ALK cDNA, transformed interleukin-3-dependent murine haematopoietic Ba/F3 cells to cytokine-independent growth. Ba/F3 cells expressing these mutations were sensitive to the small-molecule inhibitor of ALK, TAE684 (ref. 4). Furthermore, two human neuroblastoma cell lines harbouring the F1174L mutation were also sensitive to the inhibitor. Cytotoxicity was associated with increased amounts of apoptosis as measured by TdT-mediated dUTP nick end labelling (TUNEL). Short hairpin RNA (shRNA)-mediated knockdown of ALK expression in neuroblastoma cell lines with the F1174L mutation also resulted in apoptosis and impaired cell proliferation. Thus, activating alleles of the ALK receptor tyrosine kinase are present in primary neuroblastoma tumours and in established neuroblastoma cell lines, and confer sensitivity to ALK inhibition with small molecules, providing a molecular rationale for targeted therapy of this disease.

Multiple members of the let-7 family of miRNAs are often repressed in human cancers, thereby promoting oncogenesis by derepressing targets such as HMGA2, K-Ras and c-Myc. However, the mechanism by which let-7 miRNAs are coordinately repressed is unclear. The RNA-binding proteins LIN28 and LIN28B block let-7 precursors from being processed to mature miRNAs, suggesting that their overexpression might promote malignancy through repression of let-7. Here we show that LIN28 and LIN28B are overexpressed in primary human tumors and human cancer cell lines (overall frequency approximately 15%), and that overexpression is linked to repression of let-7 family miRNAs and derepression of let-7 targets. LIN28 and LIN28b facilitate cellular transformation in vitro, and overexpression is associated with advanced disease across multiple tumor types. Our work provides a mechanism for the coordinate repression of let-7 miRNAs observed in a subset of human cancers, and associates activation of LIN28 and LIN28B with poor clinical prognosis.

Soft-tissue sarcomas, which result in approximately 10,700 diagnoses and 3,800 deaths per year in the United States, show remarkable histologic diversity, with more than 50 recognized subtypes. However, knowledge of their genomic alterations is limited. We describe an integrative analysis of DNA sequence, copy number and mRNA expression in 207 samples encompassing seven major subtypes. Frequently mutated genes included TP53 (17% of pleomorphic liposarcomas), NF1 (10.5% of myxofibrosarcomas and 8% of pleomorphic liposarcomas) and PIK3CA (18% of myxoid/round-cell liposarcomas, or MRCs). PIK3CA mutations in MRCs were associated with Akt activation and poor clinical outcomes. In myxofibrosarcomas and pleomorphic liposarcomas, we found both point mutations and genomic deletions affecting the tumor suppressor NF1. Finally, we found that short hairpin RNA (shRNA)-based knockdown of several genes amplified in dedifferentiated liposarcoma, including CDK4 and YEATS4, decreased cell proliferation. Our study yields a detailed map of molecular alterations across diverse sarcoma subtypes and suggests potential subtype-specific targets for therapy.

Aneuploidy, whole chromosome or chromosome arm imbalance, is a near-universal characteristic of human cancers. In 10,522 cancer genomes from The Cancer Genome Atlas, aneuploidy was correlated with TP53 mutation, somatic mutation rate, and expression of proliferation genes. Aneuploidy was anti-correlated with expression of immune signaling genes, due to decreased leukocyte infiltrates in high-aneuploidy samples. Chromosome arm-level alterations show cancer-specific patterns, including loss of chromosome arm 3p in squamous cancers. We applied genome engineering to delete 3p in lung cells, causing decreased proliferation rescued in part by chromosome 3 duplication. This study defines genomic and phenotypic correlates of cancer aneuploidy and provides an experimental approach to study chromosome arm aneuploidy.

Genetic alterations in signaling pathways that control cell-cycle progression, apoptosis, and cell growth are common hallmarks of cancer, but the extent, mechanisms, and co-occurrence of alterations in these pathways differ between individual tumors and tumor types. Using mutations, copy-number changes, mRNA expression, gene fusions and DNA methylation in 9,125 tumors profiled by The Cancer Genome Atlas (TCGA), we analyzed the mechanisms and patterns of somatic alterations in ten canonical pathways: cell cycle, Hippo, Myc, Notch, Nrf2, PI-3-Kinase/Akt, RTK-RAS, TGFβ signaling, p53 and β-catenin/Wnt. We charted the detailed landscape of pathway alterations in 33 cancer types, stratified into 64 subtypes, and identified patterns of co-occurrence and mutual exclusivity. Eighty-nine percent of tumors had at least one driver alteration in these pathways, and 57% percent of tumors had at least one alteration potentially targetable by currently available drugs. Thirty percent of tumors had multiple targetable alterations, indicating opportunities for combination therapy.

Whole-exome sequencing of cell-free DNA (cfDNA) could enable comprehensive profiling of tumors from blood but the genome-wide concordance between cfDNA and tumor biopsies is uncertain. Here we report ichorCNA, software that quantifies tumor content in cfDNA from 0.1× coverage whole-genome sequencing data without prior knowledge of tumor mutations. We apply ichorCNA to 1439 blood samples from 520 patients with metastatic prostate or breast cancers. In the earliest tested sample for each patient, 34% of patients have ≥10% tumor-derived cfDNA, sufficient for standard coverage whole-exome sequencing. Using whole-exome sequencing, we validate the concordance of clonal somatic mutations (88%), copy number alterations (80%), mutational signatures, and neoantigens between cfDNA and matched tumor biopsies from 41 patients with ≥10% cfDNA tumor content. In summary, we provide methods to identify patients eligible for comprehensive cfDNA profiling, revealing its applicability to many patients, and demonstrate high concordance of cfDNA and metastatic tumor whole-exome sequencing.

We previously piloted the concept of a Connectivity Map (CMap), whereby genes, drugs, and disease states are connected by virtue of common gene-expression signatures. Here, we report more than a 1,000-fold scale-up of the CMap as part of the NIH LINCS Consortium, made possible by a new, low-cost, high-throughput reduced representation expression profiling method that we term L1000. We show that L1000 is highly reproducible, comparable to RNA sequencing, and suitable for computational inference of the expression levels of 81% of non-measured transcripts. We further show that the expanded CMap can be used to discover mechanism of action of small molecules, functionally annotate genetic variants of disease genes, and inform clinical trials. The 1.3 million L1000 profiles described here, as well as tools for their analysis, are available at https://clue.io.

Evidence suggests a possible role of Fusobacterium nucleatum in colorectal carcinogenesis, especially in right-sided proximal colorectum. Considering a change in bowel contents and microbiome from proximal to distal colorectal segments, we hypothesized that the proportion of colorectal carcinoma enriched with F. nucleatum might gradually increase along the bowel subsites from rectum to cecum.

Genome sequencing has uncovered a new mutational phenomenon in cancer and congenital disorders called chromothripsis. Chromothripsis is characterized by extensive genomic rearrangements and an oscillating pattern of DNA copy number levels, all curiously restricted to one or a few chromosomes. The mechanism for chromothripsis is unknown, but we previously proposed that it could occur through the physical isolation of chromosomes in aberrant nuclear structures called micronuclei. Here, using a combination of live cell imaging and single-cell genome sequencing, we demonstrate that micronucleus formation can indeed generate a spectrum of genomic rearrangements, some of which recapitulate all known features of chromothripsis. These events are restricted to the mis-segregated chromosome and occur within one cell division. We demonstrate that the mechanism for chromothripsis can involve the fragmentation and subsequent reassembly of a single chromatid from a micronucleus. Collectively, these experiments establish a new mutational process of which chromothripsis is one extreme outcome.

Genetic and chemical perturbations impact diverse cellular phenotypes, including multiple indicators of cell health. These readouts reveal toxicity and antitumorigenic effects relevant to drug discovery and personalized medicine. We developed two customized microscopy assays, one using four targeted reagents and the other three targeted reagents, to collectively measure 70 specific cell health phenotypes including proliferation, apoptosis, reactive oxygen species, DNA damage, and cell cycle stage. We then tested an approach to predict multiple cell health phenotypes using Cell Painting, an inexpensive and scalable image-based morphology assay. In matched CRISPR perturbations of three cancer cell lines, we collected both Cell Painting and cell health data. We found that simple machine learning algorithms can predict many cell health readouts directly from Cell Painting images, at less than half the cost. We hypothesized that these models can be applied to accurately predict cell health assay outcomes for any future or existing Cell Painting dataset. For Cell Painting images from a set of 1500+ compound perturbations across multiple doses, we validated predictions by orthogonal assay readouts. We provide a web app to browse predictions: http://broad.io/cell-health-app. Our approach can be used to add cell health annotations to Cell Painting datasets.

Cancer evolution poses a central obstacle to cure, as resistant clones expand under therapeutic selection pressures. Genome sequencing of relapsed disease can nominate genomic alterations conferring resistance but sample collection lags behind, limiting therapeutic innovation. Genome-wide screens offer a complementary approach to chart the compendium of escape genotypes, anticipating clinical resistance. We report genome-wide open reading frame (ORF) resistance screens for first- and third-generation epidermal growth factor receptor (EGFR) inhibitors and a MEK inhibitor. Using serial sampling, dose gradients, and mathematical modeling, we generate genotype-fitness maps across therapeutic contexts and identify alterations that escape therapy. Our data expose varying dose-fitness relationship across genotypes, ranging from complete dose invariance to paradoxical dose dependency where fitness increases in higher doses. We predict fitness with combination therapy and compare these estimates to genome-wide fitness maps of drug combinations, identifying genotypes where combination therapy results in unexpected inferior effectiveness. These data are applied to nominate combination optimization strategies to forestall resistant disease.