Perspective on Oncogenic Processes at the End of the Beginning of Cancer Genomics.

The Cancer Genome Atlas (TCGA) has catalyzed systematic characterization of diverse genomic alterations underlying human cancers. At this historic junction marking the completion of genomic characterization of over 11,000 tumors from 33 cancer types, we present our current understanding of the molecular processes governing oncogenesis. We illustrate our insights into cancer through synthesis of the findings of the TCGA PanCancer Atlas project on three facets of oncogenesis: (1) somatic driver mutations, germline pathogenic variants, and their interactions in the tumor; (2) the influence of the tumor genome and epigenome on transcriptome and proteome; and (3) the relationship between tumor and the microenvironment, including implications for drugs targeting driver events and immunotherapies. These results will anchor future characterization of rare and common tumor types, primary and relapsed tumors, and cancers across ancestry groups and will guide the deployment of clinical genomic sequencing.

Comprehensive Characterization of Cancer Driver Genes and Mutations.

Identifying molecular cancer drivers is critical for precision oncology. Multiple advanced algorithms to identify drivers now exist, but systematic attempts to combine and optimize them on large datasets are few. We report a PanCancer and PanSoftware analysis spanning 9,423 tumor exomes (comprising all 33 of The Cancer Genome Atlas projects) and using 26 computational tools to catalog driver genes and mutations. We identify 299 driver genes with implications regarding their anatomical sites and cancer/cell types. Sequence- and structure-based analyses identified >3,400 putative missense driver mutations supported by multiple lines of evidence. Experimental validation confirmed 60%-85% of predicted mutations as likely drivers. We found that >300 MSI tumors are associated with high PD-1/PD-L1, and 57% of tumors analyzed harbor putative clinically actionable events. Our study represents the most comprehensive discovery of cancer genes and mutations to date and will serve as a blueprint for future biological and clinical endeavors.

Comprehensive Molecular Characterization of Muscle-Invasive Bladder Cancer.

We report a comprehensive analysis of 412 muscle-invasive bladder cancers characterized by multiple TCGA analytical platforms. Fifty-eight genes were significantly mutated, and the overall mutational load was associated with APOBEC-signature mutagenesis. Clustering by mutation signature identified a high-mutation subset with 75% 5-year survival. mRNA expression clustering refined prior clustering analyses and identified a poor-survival "neuronal" subtype in which the majority of tumors lacked small cell or neuroendocrine histology. Clustering by mRNA, long non-coding RNA (lncRNA), and miRNA expression converged to identify subsets with differential epithelial-mesenchymal transition status, carcinoma in situ scores, histologic features, and survival. Our analyses identified 5 expression subtypes that may stratify response to different treatments.

Mutational Strand Asymmetries in Cancer Genomes Reveal Mechanisms of DNA Damage and Repair.

Mutational processes constantly shape the somatic genome, leading to immunity, aging, cancer, and other diseases. When cancer is the outcome, we are afforded a glimpse into these processes by the clonal expansion of the malignant cell. Here, we characterize a less explored layer of the mutational landscape of cancer: mutational asymmetries between the two DNA strands. Analyzing whole-genome sequences of 590 tumors from 14 different cancer types, we reveal widespread asymmetries across mutagenic processes, with transcriptional ("T-class") asymmetry dominating UV-, smoking-, and liver-cancer-associated mutations and replicative ("R-class") asymmetry dominating POLE-, APOBEC-, and MSI-associated mutations. We report a striking phenomenon of transcription-coupled damage (TCD) on the non-transcribed DNA strand and provide evidence that APOBEC mutagenesis occurs on the lagging-strand template during DNA replication. As more genomes are sequenced, studying and classifying their asymmetries will illuminate the underlying biological mechanisms of DNA damage and repair.

KDM4A lysine demethylase induces site-specific copy gain and rereplication of regions amplified in tumors.

Acquired chromosomal instability and copy number alterations are hallmarks of cancer. Enzymes capable of promoting site-specific copy number changes have yet to be identified. Here, we demonstrate that H3K9/36me3 lysine demethylase KDM4A/JMJD2A overexpression leads to localized copy gain of 1q12, 1q21, and Xq13.1 without global chromosome instability. KDM4A-amplified tumors have increased copy gains for these same regions. 1q12h copy gain occurs within a single cell cycle, requires S phase, and is not stable but is regenerated each cell division. Sites with increased copy number are rereplicated and have increased KDM4A, MCM, and DNA polymerase occupancy. Suv39h1/KMT1A or HP1γ overexpression suppresses the copy gain, whereas H3K9/K36 methylation interference promotes gain. Our results demonstrate that overexpression of a chromatin modifier results in site-specific copy gains. This begins to establish how copy number changes could originate during tumorigenesis and demonstrates that transient overexpression of specific chromatin modulators could promote these events.

The repertoire of mutational signatures in human cancer.

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature1. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Analyses of non-coding somatic drivers in 2,658 cancer whole genomes.

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

Next-generation characterization of the Cancer Cell Line Encyclopedia.

Large panels of comprehensively characterized human cancer models, including the Cancer Cell Line Encyclopedia (CCLE), have provided a rigorous framework with which to study genetic variants, candidate targets, and small-molecule and biological therapeutics and to identify new marker-driven cancer dependencies. To improve our understanding of the molecular features that contribute to cancer phenotypes, including drug responses, here we have expanded the characterizations of cancer cell lines to include genetic, RNA splicing, DNA methylation, histone H3 modification, microRNA expression and reverse-phase protein array data for 1,072 cell lines from individuals of various lineages and ethnicities. Integration of these data with functional characterizations such as drug-sensitivity, short hairpin RNA knockdown and CRISPR-Cas9 knockout data reveals potential targets for cancer drugs and associated biomarkers. Together, this dataset and an accompanying public data portal provide a resource for the acceleration of cancer research using model cancer cell lines.

High-throughput genetic clustering of type 2 diabetes loci reveals heterogeneous mechanistic pathways of metabolic disease.

AIMS/HYPOTHESIS: Type 2 diabetes is highly polygenic and influenced by multiple biological pathways. Rapid expansion in the number of type 2 diabetes loci can be leveraged to identify such pathways. METHODS: We developed a high-throughput pipeline to enable clustering of type 2 diabetes loci based on variant-trait associations. Our pipeline extracted summary statistics from genome-wide association studies (GWAS) for type 2 diabetes and related traits to generate a matrix of 323 variants × 64 trait associations and applied Bayesian non-negative matrix factorisation (bNMF) to identify genetic components of type 2 diabetes. Epigenomic enrichment analysis was performed in 28 cell types and single pancreatic cells. We generated cluster-specific polygenic scores and performed regression analysis in an independent cohort (N=25,419) to assess for clinical relevance. RESULTS: We identified ten clusters of genetic loci, recapturing the five from our prior analysis as well as novel clusters related to beta cell dysfunction, pronounced insulin secretion, and levels of alkaline phosphatase, lipoprotein A and sex hormone-binding globulin. Four clusters related to mechanisms of insulin deficiency, five to insulin resistance and one had an unclear mechanism. The clusters displayed tissue-specific epigenomic enrichment, notably with the two beta cell clusters differentially enriched in functional and stressed pancreatic beta cell states. Additionally, cluster-specific polygenic scores were differentially associated with patient clinical characteristics and outcomes. The pipeline was applied to coronary artery disease and chronic kidney disease, identifying multiple overlapping clusters with type 2 diabetes. CONCLUSIONS/INTERPRETATION: Our approach stratifies type 2 diabetes loci into physiologically interpretable genetic clusters associated with distinct tissues and clinical outcomes. The pipeline allows for efficient updating as additional GWAS become available and can be readily applied to other conditions, facilitating clinical translation of GWAS findings. Software to perform this clustering pipeline is freely available.

Recurrent and functional regulatory mutations in breast cancer.

Genomic analysis of tumours has led to the identification of hundreds of cancer genes on the basis of the presence of mutations in protein-coding regions. By contrast, much less is known about cancer-causing mutations in non-coding regions. Here we perform deep sequencing in 360 primary breast cancers and develop computational methods to identify significantly mutated promoters. Clear signals are found in the promoters of three genes. FOXA1, a known driver of hormone-receptor positive breast cancer, harbours a mutational hotspot in its promoter leading to overexpression through increased E2F binding. RMRP and NEAT1, two non-coding RNA genes, carry mutations that affect protein binding to their promoters and alter expression levels. Our study shows that promoter regions harbour recurrent mutations in cancer with functional consequences and that the mutations occur at similar frequencies as in coding regions. Power analyses indicate that more such regions remain to be discovered through deep sequencing of adequately sized cohorts of patients.

Hypoxia drives transient site-specific copy gain and drug-resistant gene expression.

Copy number heterogeneity is a prominent feature within tumors. The molecular basis for this heterogeneity remains poorly characterized. Here, we demonstrate that hypoxia induces transient site-specific copy gains (TSSGs) in primary, nontransformed, and transformed human cells. Hypoxia-driven copy gains are not dependent on HIF1α or HIF2α; however, they are dependent on the KDM4A histone demethylase and are blocked by inhibition of KDM4A with a small molecule or the natural metabolite succinate. Furthermore, this response is conserved at a syntenic region in zebrafish cells. Regions with site-specific copy gain are also enriched for amplifications in hypoxic primary tumors. These tumors exhibited amplification and overexpression of the drug resistance gene CKS1B, which we recapitulated in hypoxic breast cancer cells. Our results demonstrate that hypoxia provides a biological stimulus to create transient site-specific copy alterations that could result in heterogeneity within tumors and cell populations. These findings have major implications in our understanding of copy number heterogeneity and the emergence of drug resistance genes in cancer.

Quantification of somatic mutation flow across individual cell division events by lineage sequencing.

Mutation data reveal the dynamic equilibrium between DNA damage and repair processes in cells and are indispensable to the understanding of age-related diseases, tumor evolution, and the acquisition of drug resistance. However, available genome-wide methods have a limited ability to resolve rare somatic variants and the relationships between these variants. Here, we present lineage sequencing, a new genome sequencing approach that enables somatic event reconstruction by providing quality somatic mutation call sets with resolution as high as the single-cell level in subject lineages. Lineage sequencing entails sampling single cells from a population and sequencing subclonal sample sets derived from these cells such that knowledge of relationships among the cells can be used to jointly call variants across the sample set. This approach integrates data from multiple sequence libraries to support each variant and precisely assigns mutations to lineage segments. We applied lineage sequencing to a human colon cancer cell line with a DNA polymerase epsilon (POLE) proofreading deficiency (HT115) and a human retinal epithelial cell line immortalized by constitutive telomerase expression (RPE1). Cells were cultured under continuous observation to link observed single-cell phenotypes with single-cell mutation data. The high sensitivity, specificity, and resolution of the data provide a unique opportunity for quantitative analysis of variation in mutation rate, spectrum, and correlations among variants. Our data show that mutations arrive with nonuniform probability across sublineages and that DNA lesion dynamics may cause strong correlations between certain mutations.

Subependymal giant cell astrocytomas are characterized by mTORC1 hyperactivation, a very low somatic mutation rate, and a unique gene expression profile.

Subependymal giant-cell astrocytomas (SEGAs) are slow-growing brain tumors that are a hallmark feature seen in 5-10% of patients with Tuberous Sclerosis Complex (TSC). Though histologically benign, they can cause serious neurologic symptoms, leading to death if untreated. SEGAs consistently show biallelic loss of TSC1 or TSC2. Herein, we aimed to define other somatic events beyond TSC1/TSC2 loss and identify potential transcriptional drivers that contribute to SEGA formation. Paired tumor-normal whole-exome sequencing was performed on 21 resected SEGAs from 20 TSC patients. Pathogenic variants in TSC1/TSC2 were identified in 19/21 (90%) SEGAs. Copy neutral loss of heterozygosity (size range: 2.2-46 Mb) was seen in 76% (16/21) of SEGAs (44% chr9q and 56% chr16p). An average of 1.4 other somatic variants (range 0-7) per tumor were identified, unlikely of pathogenic significance. Whole transcriptome RNA-sequencing analyses revealed 190 common differentially expressed genes in SEGA (n = 16, 13 from a prior study) in pairwise comparison to each of: low grade diffuse gliomas (n = 530) and glioblastoma (n = 171) from The Cancer Genome Atlas (TCGA) consortium, ganglioglioma (n = 10), TSC cortical tubers (n = 15), and multiple normal tissues. Among these, homeobox transcription factors (TFs) HMX3, HMX2, VAX1, SIX3; and TFs IRF6 and EOMES were all expressed >12-fold higher in SEGAs (FDR/q-value < 0.05). Immunohistochemistry supported the specificity of IRF6, VAX1, SIX3 for SEGAs in comparison to other tumor entities and normal brain. We conclude that SEGAs have an extremely low somatic mutation rate, suggesting that TSC1/TSC2 loss is sufficient to drive tumor growth. The unique and highly expressed SEGA-specific TFs likely reflect the neuroepithelial cell of origin, and may also contribute to the transcriptional and epigenetic state that enables SEGA growth following two-hit loss of TSC1 or TSC2 and mTORC1 activation.

Genomic analyses of PMBL reveal new drivers and mechanisms of sensitivity to PD-1 blockade.

Primary mediastinal large B-cell lymphomas (PMBLs) are aggressive tumors that typically present as large mediastinal masses in young women. PMBLs share clinical, transcriptional, and molecular features with classical Hodgkin lymphoma (cHL), including constitutive activation of nuclear factor κB (NF-κB), JAK/STAT signaling, and programmed cell death protein 1 (PD-1)-mediated immune evasion. The demonstrated efficacy of PD-1 blockade in relapsed/refractory PMBLs led to recent approval by the US Food and Drug Administration and underscored the importance of characterizing targetable genetic vulnerabilities in this disease. Here, we report a comprehensive analysis of recurrent genetic alterations -somatic mutations, somatic copy number alterations, and structural variants-in a cohort of 37 newly diagnosed PMBLs. We identified a median of 9 genetic drivers per PMBL, including known and newly identified components of the JAK/STAT and NF-κB signaling pathways and frequent B2M alterations that limit major histocompatibility complex class I expression, as in cHL. PMBL also exhibited frequent, newly identified driver mutations in ZNF217 and an additional epigenetic modifier, EZH2. The majority of these alterations were clonal, which supports their role as early drivers. In PMBL, we identified several previously uncharacterized molecular features that may increase sensitivity to PD-1 blockade, including high tumor mutational burden, microsatellite instability, and an apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) mutational signature. The shared genetic features between PMBL and cHL provide a framework for analyzing the mechanism of action of PD-1 blockade in these related lymphoid malignancies.

A model combining clinical and genomic factors to predict response to PD-1/PD-L1 blockade in advanced urothelial carcinoma.

BACKGROUND: In metastatic urothelial carcinoma (mUC), predictive biomarkers that correlate with response to immune checkpoint inhibitors (ICIs) are lacking. Here, we interrogated genomic and clinical features associated with response to ICIs in mUC. METHODS: Sixty two mUC patients treated with ICI who had targeted tumour sequencing were studied. We examined associations between candidate biomarkers and clinical benefit (CB, any objective reduction in tumour size) versus no clinical benefit (NCB, no change or objective increase in tumour size). Both univariable and multivariable analyses for associations were conducted. A comparator cohort of 39 mUC patients treated with taxanes was analysed by using the same methodology. RESULTS: Nine clinical and seven genomic factors correlated with clinical outcomes in univariable analysis in the ICI cohort. Among the 16 factors, neutrophil-to-lymphocyte ratio (NLR) ≥5 (OR = 0.12, 95% CI, 0.01-1.15), visceral metastasis (OR = 0.05, 95% CI, 0.01-0.43) and single-nucleotide variant (SNV) count < 10 (OR = 0.04, 95% CI, 0.006-0.27) were identified as independent predictors of NCB to ICI in multivariable analysis (c-statistic = 0.90). None of the 16 variables were associated with clinical benefit in the taxane cohort. CONCLUSIONS: This three-factor model includes genomic (SNV count >9) and clinical (NLR <5, lack of visceral metastasis) variables predictive for benefit to ICI but not taxane therapy for mUC. External validation of these hypothesis-generating results is warranted to enable use in routine clinical care.

Mutational heterogeneity in cancer and the search for new cancer-associated genes.

Major international projects are underway that are aimed at creating a comprehensive catalogue of all the genes responsible for the initiation and progression of cancer. These studies involve the sequencing of matched tumour-normal samples followed by mathematical analysis to identify those genes in which mutations occur more frequently than expected by random chance. Here we describe a fundamental problem with cancer genome studies: as the sample size increases, the list of putatively significant genes produced by current analytical methods burgeons into the hundreds. The list includes many implausible genes (such as those encoding olfactory receptors and the muscle protein titin), suggesting extensive false-positive findings that overshadow true driver events. We show that this problem stems largely from mutational heterogeneity and provide a novel analytical methodology, MutSigCV, for resolving the problem. We apply MutSigCV to exome sequences from 3,083 tumour-normal pairs and discover extraordinary variation in mutation frequency and spectrum within cancer types, which sheds light on mutational processes and disease aetiology, and in mutation frequency across the genome, which is strongly correlated with DNA replication timing and also with transcriptional activity. By incorporating mutational heterogeneity into the analyses, MutSigCV is able to eliminate most of the apparent artefactual findings and enable the identification of genes truly associated with cancer.