Cancer primarily develops because of somatic alterations in the genome. Advances in sequencing have enabled large-scale sequencing studies across many tumor types, emphasizing the discovery of ...alterations in protein-coding genes. However, the protein-coding exome comprises less than 2% of the human genome. Here we analyze the complete genome sequences of 863 human tumors from The Cancer Genome Atlas and other sources to systematically identify noncoding regions that are recurrently mutated in cancer. We use new frequency- and sequence-based approaches to comprehensively scan the genome for noncoding mutations with potential regulatory impact. These methods identify recurrent mutations in regulatory elements upstream of PLEKHS1, WDR74 and SDHD, as well as previously identified mutations in the TERT promoter. SDHD promoter mutations are frequent in melanoma and are associated with reduced gene expression and poor prognosis. The non-protein-coding cancer genome remains widely unexplored, and our findings represent a step toward targeting the entire genome for clinical purposes.
Tumors with mismatch repair deficiency (MMR-d) are characterized by sequence alterations in microsatellites and can accumulate thousands of mutations. This high mutational burden renders tumors ...immunogenic and sensitive to programmed cell death-1 (PD-1) immune checkpoint inhibitors. Yet, despite their tumor immunogenicity, patients with MMR-deficient tumors experience highly variable responses, and roughly half are refractory to treatment. We present experimental and clinical evidence showing that the degree of microsatellite instability (MSI) and resultant mutational load, in part, underlies the variable response to PD-1 blockade immunotherapy in MMR-d human and mouse tumors. The extent of response is particularly associated with the accumulation of insertion-deletion (indel) mutational load. This study provides a rationale for the genome-wide characterization of MSI intensity and mutational load to better profile responses to anti-PD-1 immunotherapy across MMR-deficient human cancers.
BRCA1 and BRCA2 are involved in homologous recombination (HR) DNA repair and are germ-line cancer pre-disposition genes that result in a syndrome of hereditary breast and ovarian cancer (HBOC). ...Whether germ-line or somatic alterations in these genes or other members of the HR pathway and if mono- or bi-allelic alterations of HR-related genes have a phenotypic impact on other cancers remains to be fully elucidated. Here, we perform a pan-cancer analysis of The Cancer Genome Atlas (TCGA) data set and observe that bi-allelic pathogenic alterations in homologous recombination (HR) DNA repair-related genes are prevalent across many malignancies. These bi-allelic alterations often associate with genomic features of HR deficiency. Further, in ovarian, breast and prostate cancers, bi-allelic alterations are mutually exclusive of each other. The combination of these two properties facilitates reclassification of variants of unknown significance affecting DNA repair genes, and may help personalize HR directed therapies in the clinic.Germline mutations in homologous recombination (HR) DNA repair genes are linked to breast and ovarian cancer. Here, the authors show that mutually exclusive bi-allelic inactivation of HR genes are present in other cancer types and associated with genomic features of HR deficiency, expanding the potential use of HR-directed therapies.
Tumor-infiltrating immune cells have been linked to prognosis and response to immunotherapy; however, the levels of distinct immune cell subsets and the signals that draw them into a tumor, such as ...the expression of antigen presenting machinery genes, remain poorly characterized. Here, we employ a gene expression-based computational method to profile the infiltration levels of 24 immune cell populations in 19 cancer types.
We compare cancer types using an immune infiltration score and a T cell infiltration score and find that clear cell renal cell carcinoma (ccRCC) is among the highest for both scores. Using immune infiltration profiles as well as transcriptomic and proteomic datasets, we characterize three groups of ccRCC tumors: T cell enriched, heterogeneously infiltrated, and non-infiltrated. We observe that the immunogenicity of ccRCC tumors cannot be explained by mutation load or neo-antigen load, but is highly correlated with MHC class I antigen presenting machinery expression (APM). We explore the prognostic value of distinct T cell subsets and show in two cohorts that Th17 cells and CD8
T/Treg ratio are associated with improved survival, whereas Th2 cells and Tregs are associated with negative outcomes. Investigation of the association of immune infiltration patterns with the subclonal architecture of tumors shows that both APM and T cell levels are negatively associated with subclone number.
Our analysis sheds light on the immune infiltration patterns of 19 human cancers and unravels mRNA signatures with prognostic utility and immunotherapeutic biomarker potential in ccRCC.
Only a fraction of patients with cancer respond to immune checkpoint blockade (ICB) treatment, but current decision-making procedures have limited accuracy. In this study, we developed a machine ...learning model to predict ICB response by integrating genomic, molecular, demographic and clinical data from a comprehensively curated cohort (MSK-IMPACT) with 1,479 patients treated with ICB across 16 different cancer types. In a retrospective analysis, the model achieved high sensitivity and specificity in predicting clinical response to immunotherapy and predicted both overall survival and progression-free survival in the test data across different cancer types. Our model significantly outperformed predictions based on tumor mutational burden, which was recently approved by the U.S. Food and Drug Administration for this purpose
. Additionally, the model provides quantitative assessments of the model features that are most salient for the predictions. We anticipate that this approach will substantially improve clinical decision-making in immunotherapy and inform future interventions.
Tumors with somatic mutations in the proofreading exonuclease domain of DNA polymerase epsilon (POLE-exo*) exhibit a novel mutator phenotype, with markedly elevated TCT→TAT and TCG→TTG mutations and ...overall mutation frequencies often exceeding 100 mutations/Mb. Here, we identify POLE-exo* tumors in numerous cancers and classify them into two groups, A and B, according to their mutational properties. Group A mutants are found only in POLE, whereas Group B mutants are found in POLE and POLD1 and appear to be nonfunctional. In Group A, cell-free polymerase assays confirm that mutations in the exonuclease domain result in high mutation frequencies with a preference for C→A mutation. We describe the patterns of amino acid substitutions caused by POLE-exo* and compare them to other tumor types. The nucleotide preference of POLE-exo* leads to increased frequencies of recurrent nonsense mutations in key tumor suppressors such as TP53, ATM, and PIK3R1. We further demonstrate that strand-specific mutation patterns arise from some of these POLE-exo* mutants during genome duplication. This is the first direct proof of leading strand-specific replication by human POLE, which has only been demonstrated in yeast so far. Taken together, the extremely high mutation frequency and strand specificity of mutations provide a unique identifier of eukaryotic origins of replication.
We sought to define the prevalence and co-occurrence of actionable genomic alterations in patients with high-grade bladder cancer to serve as a platform for therapeutic drug discovery.
An integrative ...analysis of 97 high-grade bladder tumors was conducted to identify actionable drug targets, which are defined as genomic alterations that have been clinically validated in another cancer type (eg, BRAF mutation) or alterations for which a selective inhibitor of the target or pathway is under clinical investigation. DNA copy number alterations (CNAs) were defined by using array comparative genomic hybridization. Mutation profiling was performed by using both mass spectroscopy-based genotyping and Sanger sequencing.
Sixty-one percent of tumors harbored potentially actionable genomic alterations. A core pathway analysis of the integrated data set revealed a nonoverlapping pattern of mutations in the RTK-RAS-RAF and phosphoinositide 3-kinase/AKT/mammalian target of rapamycin pathways and regulators of G1-S cell cycle progression. Unsupervised clustering of CNAs defined two distinct classes of bladder tumors that differed in the degree of their CNA burden. Integration of mutation and copy number analyses revealed that mutations in TP53 and RB1 were significantly more common in tumors with a high CNA burden (P < .001 and P < .003, respectively).
High-grade bladder cancer possesses substantial genomic heterogeneity. The majority of tumors harbor potentially tractable genomic alterations that may predict for response to target-selective agents. Given the genomic diversity of bladder cancers, optimal development of target-specific agents will require pretreatment genomic characterization.
Genome wide association studies (GWAS) of human diseases have generally identified many loci associated with risk with relatively small effect sizes. The omnigenic model attempts to explain this ...observation by suggesting that diseases can be thought of as networks, where genes with direct involvement in disease-relevant biological pathways are named 'core genes', while peripheral genes influence disease risk via their interactions or regulatory effects on core genes. Here, we demonstrate a method for identifying candidate core genes solely from genes in or near disease-associated SNPs (GWAS hits) in conjunction with protein-protein interaction network data. Applied to 1,381 GWAS studies from 5 ancestries, we identify a total of 1,865 candidate core genes in 343 GWAS studies. Our analysis identifies several well-known disease-related genes that are not identified by GWAS, including BRCA1 in Breast Cancer, Amyloid Precursor Protein (APP) in Alzheimer's Disease, INS in A1C measurement and Type 2 Diabetes, and PCSK9 in LDL cholesterol, amongst others. Notably candidate core genes are preferentially enriched for disease relevance over GWAS hits and are enriched for both Clinvar pathogenic variants and known drug targets-consistent with the predictions of the omnigenic model. We subsequently use parent term annotations provided by the GWAS catalog, to merge related GWAS studies and identify candidate core genes in over-arching disease processes such as cancer-where we identify 109 candidate core genes.
PBRM1 is the second most commonly mutated gene after VHL in clear cell renal cell carcinoma (ccRCC). However, the biological consequences of PBRM1 mutations for kidney tumorigenesis are unknown. ...Here, we find that kidney-specific deletion of Vhl and Pbrm1, but not either gene alone, results in bilateral, multifocal, transplantable clear cell kidney cancers. PBRM1 loss amplified the transcriptional outputs of HIF1 and STAT3 incurred by Vhl deficiency. Analysis of mouse and human ccRCC revealed convergence on mTOR activation, representing the third driver event after genetic inactivation of VHL and PBRM1. Our study reports a physiological preclinical ccRCC mouse model that recapitulates somatic mutations in human ccRCC and provides mechanistic and therapeutic insights into PBRM1 mutated subtypes of human ccRCC.
Display omitted
•PBRM1 is a bona fide tumor suppressor in the pathogenesis of ccRCC•PBRM1 prevents self-perpetuating amplification of HIF1/STAT3 signaling in Vhl−/− cell•Loss of Vhl and Pbrm1 in mouse kidney results in multifocal, transplantable ccRCC•In ccRCC, mTORC1 activation is the third driver event after loss of VHL and PBRM1
Nargund et al. present a three-step process in the pathogenesis of mouse and human clear cell kidney cancer. After the loss of VHL, the loss of SWI/SNF tumor suppressor protein PBRM1/BAF180 further activates HIF1/STAT3 signaling in mouse kidney and positions mTORC1 activation as the preferred third driver event.
Many long noncoding RNA (lncRNA) genes exist near cancer-associated loci, yet evidence connecting lncRNA functions to recurrent genetic alterations in cancer are lacking. Here, we report that DINO, ...the lncRNA transcribed from the cancer-associated DINO/CDKN1A locus, suppresses tumor formation independent of p21, the protein encoded at the locus. Loss of one or two alleles of Dino impairs p53 signaling and apoptosis, resulting in a haplo-insufficient tumor suppressor phenotype in genetically defined mouse models of tumorigenesis. A discrete region of the DINO/CDKN1A locus is recurrently hypermethylated in human cancers, silencing DINO but not CDKN1A, the gene encoding p21. Hypermethylation silences DINO, impairs p53 signaling pathway in trans, and is mutually exclusive with TP53 alterations, indicating that DINO and TP53 comprise a common tumor suppressor module. Therefore, DINO encodes a lncRNA essential for tumor suppression that is recurrently silenced in human cancers as a mechanism to escape p53-dependent tumor suppression.
Display omitted
•DINO hypermethylation in human cancers is mutually exclusive with TP53 mutations•Dino is a haplo-insufficient tumor suppressor encoded in the Dino/Cdkn1a locus•Dino’s tumor suppressor mechanism requires p53•Hypermethylation silences DINO in human cancers to impair p53 signaling
Tumorigenesis requires escape from the p53 tumor suppressor protein. How 50% of cancers develop despite an intact p53 pathway remains largely unknown. Marney et al. report here that loss of the lncRNA DINO provides a route to tumorigenesis in p53-intact cells, and DINO is recurrently silenced in human cancers.