A glycolytic metabolite bypasses "two-hit" tumor suppression by BRCA2.

Knudson's "two-hit" paradigm posits that carcinogenesis requires inactivation of both copies of an autosomal tumor suppressor gene. Here, we report that the glycolytic metabolite methylglyoxal (MGO) transiently bypasses Knudson's paradigm by inactivating the breast cancer suppressor protein BRCA2 to elicit a cancer-associated, mutational single-base substitution (SBS) signature in nonmalignant mammary cells or patient-derived organoids. Germline monoallelic BRCA2 mutations predispose to these changes. An analogous SBS signature, again without biallelic BRCA2 inactivation, accompanies MGO accumulation and DNA damage in Kras-driven, Brca2-mutant murine pancreatic cancers and human breast cancers. MGO triggers BRCA2 proteolysis, temporarily disabling BRCA2's tumor suppressive functions in DNA repair and replication, causing functional haploinsufficiency. Intermittent MGO exposure incites episodic SBS mutations without permanent BRCA2 inactivation. Thus, a metabolic mechanism wherein MGO-induced BRCA2 haploinsufficiency transiently bypasses Knudson's two-hit requirement could link glycolysis activation by oncogenes, metabolic disorders, or dietary challenges to mutational signatures implicated in cancer evolution.

Mutational signature assignment heterogeneity is widespread and can be addressed by ensemble approaches.

Single-base substitution (SBS) mutational signatures have become standard practice in cancer genomics. In lieu of de novo signature extraction, reference signature assignment allows users to estimate the activities of pre-established SBS signatures within individual malignancies. Several tools have been developed for this purpose, each with differing methodologies. However, due to a lack of standardization, there may be inter-tool variability in signature assignment. We deeply characterized three assignment strategies and five SBS signature assignment tools. We observed that assignment strategy choice can significantly influence results and interpretations. Despite varying recommendations by tools, Refit performed best by reducing overfitting and maximizing reconstruction of the original mutational spectra. Even after uniform application of Refit, tools varied remarkably in signature assignments both qualitatively (Jaccard index = 0.38-0.83) and quantitatively (Kendall tau-b = 0.18-0.76). This phenomenon was exacerbated for 'flat' signatures such as the homologous recombination deficiency signature SBS3. An ensemble approach (EnsembleFit), which leverages output from all five tools, increased SBS3 assignment accuracy in BRCA1/2-deficient breast carcinomas. After generating synthetic mutational profiles for thousands of pan-cancer tumors, EnsembleFit reduced signature activity assignment error 15.9-24.7% on average using Catalogue of Somatic Mutations In Cancer and non-standard reference signature sets. We have also released the EnsembleFit web portal (https://www.ensemblefit.pittlabgenomics.com) for users to generate or download ensemble-based SBS signature assignments using any strategy and combination of tools. Overall, we show that signature assignment heterogeneity across tools and strategies is non-negligible and propose a viable, ensemble solution.

Immune pathway upregulation and lower genomic instability distinguish EBV-positive nodal T/NK-cell lymphoma from ENKTL and PTCL-NOS.

Primary Epstein-Barr virus (EBV)-positive nodal T/NK-cell lymphoma (PTCL-EBV) is a poorly understood disease which shows features resembling extranodal NK/T-cell lymphoma (ENKTL) and is currently not recognized as a distinct entity but categorized as a variant of primary T-cell lymphoma not otherwise specified (PTCL-NOS). Herein, we analyzed copynumber aberrations (n=77) with a focus on global measures of genomic instability and homologous recombination deficiency and performed gene expression (n=84) and EBV miRNA expression (n=24) profiling as well as targeted mutational analysis (n=16) to further characterize PTCL-EBV in relation to ENKTL and PTCL-NOS. Multivariate analysis revealed that patients with PTCL-EBV had a significantly worse outcome compared to patients with PTCL-NOS (P=0.002) but not to those with ENKTL. Remarkably, PTCL-EBV exhibited significantly lower genomic instability and homologous recombination deficiency scores compared to ENKTL and PTCL-NOS. Gene set enrichment analysis revealed that many immune-related pathways, interferon α/γ response, and IL6_JAK_STAT3 signaling were significantly upregulated in PTCLEBV and correlated with lower genomic instability scores. We also identified that NFκB-associated genes, BIRC3, NFKB1 (P50) and CD27, and their proteins are upregulated in PTCL-EBV. Most PTCL-EBV demonstrated a type 2 EBV latency pattern and, strikingly, exhibited downregulated expression of most EBV miRNA compared to ENKTL and their target genes were also enriched in immune-related pathways. PTCL-EBV also showed frequent mutations of TET2, PIK3CD and STAT3, and are characterized by microsatellite stability. Overall, poor outcome, low genomic instability, upregulation of immune pathways and downregulation of EBV miRNA are distinctive features of PTCL-EBV. Our data support the concept that PTCL-EBV could be considered as a distinct entity, provide novel insights into the pathogenesis of the disease and offer potential new therapeutic targets for this tumor.

Germline variants and somatic mutation signatures of breast cancer across populations of African and European ancestry in the US and Nigeria.

Somatic mutation signatures may represent footprints of genetic and environmental exposures that cause different cancer. Few studies have comprehensively examined their association with germline variants, and none in an indigenous African population. SomaticSignatures was employed to extract mutation signatures based on whole-genome or whole-exome sequencing data from female patients with breast cancer (TCGA, training set, n = 1,011; Nigerian samples, validation set, n = 170), and to estimate contributions of signatures in each sample. Association between somatic signatures and common single nucleotide polymorphisms (SNPs) or rare deleterious variants were examined using linear regression. Nine stable signatures were inferred, and four signatures (APOBEC C>T, APOBEC C>G, aging and homologous recombination deficiency) were highly similar to known COSMIC signatures and explained the majority (60-85%) of signature contributions. There were significant heritable components associated with APOBEC C>T signature (h2 = 0.575, p = 0.010) and the combined APOBEC signatures (h2 = 0.432, p = 0.042). In TCGA dataset, seven common SNPs within or near GNB5 were significantly associated with an increased proportion (beta = 0.33, 95% CI = 0.21-0.45) of APOBEC signature contribution at genome-wide significance, while rare germline mutations in MTCL1 was also significantly associated with a higher contribution of this signature (p = 6.1 × 10-6 ). This is the first study to identify associations between germline variants and mutational patterns in breast cancer across diverse populations and geography. The findings provide evidence to substantiate causal links between germline genetic risk variants and carcinogenesis.

Identification of novel molecular regulators of tumor necrosis factor-related apoptosis-inducing ligand (TRAIL)-induced apoptosis in breast cancer cells by RNAi screening.

INTRODUCTION: Tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) binds to its receptors, TRAIL-receptor 1 (TRAIL-R1) and TRAIL-receptor 2 (TRAIL-R2), leading to apoptosis by activation of caspase-8 and the downstream executioner caspases, caspase-3 and caspase-7 (caspase-3/7). Triple-negative breast cancer (TNBC) cell lines with a mesenchymal phenotype are sensitive to TRAIL, whereas other breast cancer cell lines are resistant. The underlying mechanisms that control TRAIL sensitivity in breast cancer cells are not well understood. Here, we performed small interfering RNA (siRNA) screens to identify molecular regulators of the TRAIL pathway in breast cancer cells. METHODS: We conducted siRNA screens of the human kinome (691 genes), phosphatome (320 genes), and about 300 additional genes in the mesenchymal TNBC cell line MB231. Forty-eight hours after transfection of siRNA, parallel screens measuring caspase-8 activity, caspase-3/7 activity, or cell viability were conducted in the absence or presence of TRAIL for each siRNA, relative to a negative control siRNA (siNeg). A subset of genes was screened in cell lines representing epithelial TNBC (MB468), HER2-amplified breast cancer (SKBR3), and estrogen receptor-positive breast cancer (T47D). Selected putative negative regulators of the TRAIL pathway were studied by using small-molecule inhibitors. RESULTS: The primary screens in MB231 identified 150 genes, including 83 kinases, 4 phosphatases, and 63 nonkinases, as potential negative regulators of TRAIL. The identified genes are involved in many critical cell processes, including apoptosis, growth factor-receptor signaling, cell-cycle regulation, transcriptional regulation, and DNA repair. Gene-network analysis identified four genes (PDPK1, IKBKB, SRC, and BCL2L1) that formed key nodes within the interaction network of negative regulators. A secondary screen of a subset of the genes identified in additional cell lines representing different breast cancer subtypes and sensitivities to TRAIL validated and extended these findings. Further, we confirmed that small-molecule inhibition of SRC or BCL2L1, in combination with TRAIL, sensitizes breast cancer cells to TRAIL-induced apoptosis, including cell lines resistant to TRAIL-induced cytotoxicity. CONCLUSIONS: These data identify novel molecular regulators of TRAIL-induced apoptosis in breast cancer cells and suggest strategies for the enhanced application of TRAIL as a therapy for breast cancer.

Patterns of Oncogene Coexpression at Single-Cell Resolution Influence Survival in Lymphoma.

Cancers often overexpress multiple clinically relevant oncogenes, but it is not known if combinations of oncogenes in cellular subpopulations within a cancer influence clinical outcomes. Using quantitative multispectral imaging of the prognostically relevant oncogenes MYC, BCL2, and BCL6 in diffuse large B-cell lymphoma (DLBCL), we show that the percentage of cells with a unique combination MYC+BCL2+BCL6- (M+2+6-) consistently predicts survival across four independent cohorts (n = 449), an effect not observed with other combinations including M+2+6+. We show that the M+2+6- percentage can be mathematically derived from quantitative measurements of the individual oncogenes and correlates with survival in IHC (n = 316) and gene expression (n = 2,521) datasets. Comparative bulk/single-cell transcriptomic analyses of DLBCL samples and MYC/BCL2/BCL6-transformed primary B cells identify molecular features, including cyclin D2 and PI3K/AKT as candidate regulators of M+2+6- unfavorable biology. Similar analyses evaluating oncogenic combinations at single-cell resolution in other cancers may facilitate an understanding of cancer evolution and therapy resistance. SIGNIFICANCE: Using single-cell-resolved multiplexed imaging, we show that selected subpopulations of cells expressing specific combinations of oncogenes influence clinical outcomes in lymphoma. We describe a probabilistic metric for the estimation of cellular oncogenic coexpression from IHC or bulk transcriptomes, with possible implications for prognostication and therapeutic target discovery in cancer. This article is highlighted in the In This Issue feature, p. 1027.

Systems-wide RNAi analysis of CASP8AP2/FLASH shows transcriptional deregulation of the replication-dependent histone genes and extensive effects on the transcriptome of colorectal cancer cells.

BACKGROUND: Colorectal carcinomas (CRC) carry massive genetic and transcriptional alterations that influence multiple cellular pathways. The study of proteins whose loss-of-function (LOF) alters the growth of CRC cells can be used to further understand the cellular processes cancer cells depend upon for survival. RESULTS: A small-scale RNAi screen of ~400 genes conducted in SW480 CRC cells identified several candidate genes as required for the viability of CRC cells, most prominently CASP8AP2/FLASH. To understand the function of this gene in maintaining the viability of CRC cells in an unbiased manner, we generated gene specific expression profiles following RNAi. Silencing of CASP8AP2/FLASH resulted in altered expression of over 2500 genes enriched for genes associated with cellular growth and proliferation. Loss of CASP8AP2/FLASH function was significantly associated with altered transcription of the genes encoding the replication-dependent histone proteins as a result of the expression of the non-canonical polyA variants of these transcripts. Silencing of CASP8AP2/FLASH also mediated enrichment of changes in the expression of targets of the NFκB and MYC transcription factors. These findings were confirmed by whole transcriptome analysis of CASP8AP2/FLASH silenced cells at multiple time points. Finally, we identified and validated that CASP8AP2/FLASH LOF increases the expression of neurofilament heavy polypeptide (NEFH), a protein recently linked to regulation of the AKT1/ß-catenin pathway. CONCLUSIONS: We have used unbiased RNAi based approaches to identify and characterize the function of CASP8AP2/FLASH, a protein not previously reported as required for cell survival. This study further defines the role CASP8AP2/FLASH plays in the regulating expression of the replication-dependent histones and shows that its LOF results in broad and reproducible effects on the transcriptome of colorectal cancer cells including the induction of expression of the recently described tumor suppressor gene NEFH.

Pan-cancer pervasive upregulation of 3' UTR splicing drives tumourigenesis.

Most mammalian genes generate messenger RNAs with variable untranslated regions (UTRs) that are important post-transcriptional regulators. In cancer, shortening at 3' UTR ends via alternative polyadenylation can activate oncogenes. However, internal 3' UTR splicing remains poorly understood as splicing studies have traditionally focused on protein-coding alterations. Here we systematically map the pan-cancer landscape of 3' UTR splicing and present this in SpUR ( http://www.cbrc.kaust.edu.sa/spur/home/ ). 3' UTR splicing is widespread, upregulated in cancers, correlated with poor prognosis and more prevalent in oncogenes. We show that antisense oligonucleotide-mediated inhibition of 3' UTR splicing efficiently reduces oncogene expression and impedes tumour progression. Notably, CTNNB1 3' UTR splicing is the most consistently dysregulated event across cancers. We validate its upregulation in hepatocellular carcinoma and colon adenocarcinoma, and show that the spliced 3' UTR variant is the predominant contributor to its oncogenic functions. Overall, our study highlights the importance of 3' UTR splicing in cancer and may launch new avenues for RNA-based anti-cancer therapeutics.

Identification of the receptor tyrosine kinase AXL in breast cancer as a target for the human miR-34a microRNA.

The identification of molecular features that contribute to the progression of breast cancer can provide valuable insight into the pathogenesis of this disease. Deregulated microRNA expression represents one type of molecular event that has been associated with many different human cancers. In order to identify a miRNA/mRNA regulatory interaction that is biologically relevant to the triple-negative breast cancer genotype/phenotype, we initially conducted a miRNA profiling experiment to detect differentially expressed miRNAs in cell line models representing triple-negative (MDA-MB-231), ER(+) (MCF7), and HER-2 overexpressed (SK-BR-3) histotypes. We identified human miR-34a expression as being >3-fold down (from its median expression value across all cell lines) in MDA-MB-231 cells, and identified AXL as a putative mRNA target using multiple miRNA/target prediction algorithms. The miR-34a/AXL interaction was functionally characterized through ectopic overexpression experiments with a miR-34a mimic in two independent triple-negative breast cancer cell lines. In reporter assays, miR-34a binds to its putative target site within the AXL 3'UTR to inhibit luciferase expression. We also observed degradation of AXL mRNA and decreased AXL protein levels, as well as cell signaling effects on AKT phosphorylation and phenotypic effects on cell migration. Finally, we present an inverse correlative trend in miR-34a and AXL expression for both cell line and patient tumor samples.

Highly recurrent CBS epimutations in gastric cancer CpG island methylator phenotypes and inflammation.

BACKGROUND: CIMP (CpG island methylator phenotype) is an epigenetic molecular subtype, observed in multiple malignancies and associated with the epigenetic silencing of tumor suppressors. Currently, for most cancers including gastric cancer (GC), mechanisms underlying CIMP remain poorly understood. We sought to discover molecular contributors to CIMP in GC, by performing global DNA methylation, gene expression, and proteomics profiling across 14 gastric cell lines, followed by similar integrative analysis in 50 GC cell lines and 467 primary GCs. RESULTS: We identify the cystathionine beta-synthase enzyme (CBS) as a highly recurrent target of epigenetic silencing in CIMP GC. Likewise, we show that CBS epimutations are significantly associated with CIMP in various other cancers, occurring even in premalignant gastroesophageal conditions and longitudinally linked to clinical persistence. Of note, CRISPR deletion of CBS in normal gastric epithelial cells induces widespread DNA methylation changes that overlap with primary GC CIMP patterns. Reflecting its metabolic role as a gatekeeper interlinking the methionine and homocysteine cycles, CBS loss in vitro also causes reductions in the anti-inflammatory gasotransmitter hydrogen sulfide (H2S), with concomitant increase in NF-κB activity. In a murine genetic model of CBS deficiency, preliminary data indicate upregulated immune-mediated transcriptional signatures in the stomach. CONCLUSIONS: Our results implicate CBS as a bi-faceted modifier of aberrant DNA methylation and inflammation in GC and highlights H2S donors as a potential new therapy for CBS-silenced lesions.

Whole-genome analysis of Nigerian patients with breast cancer reveals ethnic-driven somatic evolution and distinct genomic subtypes.

Black women across the African diaspora experience more aggressive breast cancer with higher mortality rates than white women of European ancestry. Although inter-ethnic germline variation is known, differential somatic evolution has not been investigated in detail. Analysis of deep whole genomes of 97 breast cancers, with RNA-seq in a subset, from women in Nigeria in comparison with The Cancer Genome Atlas (n = 76) reveal a higher rate of genomic instability and increased intra-tumoral heterogeneity as well as a unique genomic subtype defined by early clonal GATA3 mutations with a 10.5-year younger age at diagnosis. We also find non-coding mutations in bona fide drivers (ZNF217 and SYPL1) and a previously unreported INDEL signature strongly associated with African ancestry proportion, underscoring the need to expand inclusion of diverse populations in biomedical research. Finally, we demonstrate that characterizing tumors for homologous recombination deficiency has significant clinical relevance in stratifying patients for potentially life-saving therapies.

Quantitative imaging of RAD51 expression as a marker of platinum resistance in ovarian cancer.

Early relapse after platinum chemotherapy in epithelial ovarian cancer (EOC) portends poor survival. A-priori identification of platinum resistance is therefore crucial to improve on standard first-line carboplatin-paclitaxel treatment. The DNA repair pathway homologous recombination (HR) repairs platinum-induced damage, and the HR recombinase RAD51 is overexpressed in cancer. We therefore designed a REMARK-compliant study of pre-treatment RAD51 expression in EOC, using fluorescent quantitative immunohistochemistry (qIHC) to overcome challenges in quantitation of protein expression in situ. In a discovery cohort (n = 284), RAD51-High tumours had shorter progression-free and overall survival compared to RAD51-Low cases in univariate and multivariate analyses. The association of RAD51 with relapse/survival was validated in a carboplatin monotherapy SCOTROC4 clinical trial cohort (n = 264) and was predominantly noted in HR-proficient cancers (Myriad HRDscore < 42). Interestingly, overexpression of RAD51 modified expression of immune-regulatory pathways in vitro, while RAD51-High tumours showed exclusion of cytotoxic T cells in situ. Our findings highlight RAD51 expression as a determinant of platinum resistance and suggest possible roles for therapy to overcome immune exclusion in RAD51-High EOC. The qIHC approach is generalizable to other proteins with a continuum instead of discrete/bimodal expression.

Author Correction: Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors.

The original version of this Article omitted a declaration from the competing interests statement, which should have included the following: 'K.P.W. is President of Tempus Lab, Inc., Chicago, IL, USA'. This has now been corrected in both the PDF and HTML versions of the Article.

Author Correction: Characterization of Nigerian breast cancer reveals prevalent homologous recombination deficiency and aggressive molecular features.

The original version of this Article contained an error in the author affiliations. The affiliation of Kevin P. White with Tempus Labs, Inc. Chicago, IL, USA was inadvertently omitted. This has now been corrected in both the PDF and HTML versions of the Article.

Genetic Ancestry May Influence the Evolutionary Trajectory of Cancers.

In this issue of Cancer Cell, Yuan et al. compared the somatic alterations harbored by tumors from European and African ancestry individuals. They determined that the latter group has a propensity for aberrations that are consistent with genomic instability, potentially lending insight to the genomic basis of cancer health disparities.

Author Correction: Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors.

The original version of this Article contained an error in the author affiliations. The affiliation of Kevin P. White with Tempus Labs, Inc., Chicago, IL, USA was inadvertently omitted.This has now been corrected in both the PDF and HTML versions of the Article.

Characterization of Nigerian breast cancer reveals prevalent homologous recombination deficiency and aggressive molecular features.

Racial/ethnic disparities in breast cancer mortality continue to widen but genomic studies rarely interrogate breast cancer in diverse populations. Through genome, exome, and RNA sequencing, we examined the molecular features of breast cancers using 194 patients from Nigeria and 1037 patients from The Cancer Genome Atlas (TCGA). Relative to Black and White cohorts in TCGA, Nigerian HR + /HER2 - tumors are characterized by increased homologous recombination deficiency signature, pervasive TP53 mutations, and greater structural variation-indicating aggressive biology. GATA3 mutations are also more frequent in Nigerians regardless of subtype. Higher proportions of APOBEC-mediated substitutions strongly associate with PIK3CA and CDH1 mutations, which are underrepresented in Nigerians and Blacks. PLK2, KDM6A, and B2M are also identified as previously unreported significantly mutated genes in breast cancer. This dataset provides novel insights into potential molecular mechanisms underlying outcome disparities and lay a foundation for deployment of precision therapeutics in underserved populations.

Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors.

Homozygous deletions are rare in cancers and often target tumour suppressor genes. Here, we build a compendium of 2218 primary tumours across 12 human cancer types and systematically screen for homozygous deletions, aiming to identify rare tumour suppressors. Our analysis defines 96 genomic regions recurrently targeted by homozygous deletions. These recurrent homozygous deletions occur either over tumour suppressors or over fragile sites, regions of increased genomic instability. We construct a statistical model that separates fragile sites from regions showing signatures of positive selection for homozygous deletions and identify candidate tumour suppressors within those regions. We find 16 established tumour suppressors and propose 27 candidate tumour suppressors. Several of these genes (including MGMT, RAD17, and USP44) show prior evidence of a tumour suppressive function. Other candidate tumour suppressors, such as MAFTRR, KIAA1551, and IGF2BP2, are novel. Our study demonstrates how rare tumour suppressors can be identified through copy number meta-analysis.

Comparison of Breast Cancer Molecular Features and Survival by African and European Ancestry in The Cancer Genome Atlas.

Importance: African Americans have the highest breast cancer mortality rate. Although racial difference in the distribution of intrinsic subtypes of breast cancer is known, it is unclear if there are other inherent genomic differences that contribute to the survival disparities. Objectives: To investigate racial differences in breast cancer molecular features and survival and to estimate the heritability of breast cancer subtypes. Design, Setting, and Participants: Among a convenience cohort of patients with invasive breast cancer, breast tumor and matched normal tissue sample data (as of September 18, 2015) were obtained from The Cancer Genome Atlas. Main Outcomes and Measures: Breast cancer­free interval, tumor molecular features, and genetic variants. Results: Participants were 930 patients with breast cancer, including 154 black patients of African ancestry (mean [SD] age at diagnosis, 55.66 [13.01] years; 98.1% [n = 151] female) and 776 white patients of European ancestry (mean [SD] age at diagnosis, 59.51 [13.11] years; 99.0% [n = 768] female). Compared with white patients, black patients had a worse breast cancer-free interval (hazard ratio, HR=1.67; 95% CI, 1.02-2.74; P = .043). They had a higher likelihood of basal-like (odds ratio, 3.80; 95% CI, 2.46-5.87; P < .001) and human epidermal growth factor receptor 2 (ERBB2 [formerly HER2])­enriched (odds ratio, 2.22; 95% CI, 1.10-4.47; P = .027) breast cancer subtypes, with the Luminal A subtype as the reference. Blacks had more TP53 mutations and fewer PIK3CA mutations than whites. While most molecular differences were eliminated after adjusting for intrinsic subtype, the study found 16 DNA methylation probes, 4 DNA copy number segments, 1 protein, and 142 genes that were differentially expressed, with the gene-based signature having an excellent capacity for distinguishing breast tumors from black vs white patients (cross-validation C index, 0.878). Using germline genotypes, the heritability of breast cancer subtypes (basal vs nonbasal) was estimated to be 0.436 (P = 1.5 × 10−14). The estrogen receptor­positive polygenic risk score built from 89 known susceptibility variants was higher in blacks than in whites (difference, 0.24; P = 2.3 × 10−5), while the estrogen receptor­negative polygenic risk score was much higher in blacks than in whites (difference, 0.48; P = 2.8 × 10−11). Conclusions and Relevance: On the molecular level, after adjusting for intrinsic subtype frequency differences, this study found a modest number of genomic differences but a significant clinical survival outcome difference between blacks and whites in The Cancer Genome Atlas data set. Moreover, more than 40% of breast cancer subtype frequency differences could be explained by genetic variants. These data could form the basis for the development of molecular targeted therapies to improve clinical outcomes for the specific subtypes of breast cancers that disproportionately affect black women. Findings also indicate that personalized risk assessment and optimal treatment could reduce deaths from aggressive breast cancers for black women.