DNA copy number motifs are strong and independent predictors of survival in breast cancer

Copy number alterations: a catastrophic orchestration of the breast cancer genome

Breast cancer (BCa) is a prevalent malignancy that predominantly affects women around the world. Somatic copy number alterations (CNAs) are tumor-specific amplifications or deletions of DNA segments that often drive BCa development and therapy resistance. Hence, the complex patterns of CNAs complement BCa classification systems. In addition, understanding the precise contributions of CNAs is essential for tailoring personalized treatment approaches. This review highlights how tumor evolution drives the acquisition of CNAs, which in turn shape the genomic landscapes of BCas. It also discusses advanced methodologies for identifying recurrent CNAs, studying CNAs in BCa and their clinical impact.

Copy number signatures and CCNE1 amplification reveal the involvement of replication stress in high-grade endometrial tumors oncogenesis

PURPOSE

Managing high-grade endometrial cancer in Martinique poses significant challenges. The diversity of copy number alterations in high-grade endometrial tumors, often associated with a TP53 mutation, is a key factor complicating treatment. Due to the high incidence of high-grade tumors with poor prognosis, our study aimed to characterize the molecular signature of these tumors within a cohort of 25 high-grade endometrial cases.

METHODS

We conducted a comprehensive pangenomic analysis to categorize the copy number alterations involved in these tumors. Whole-Exome Sequencing (WES) and Homologous Recombination (HR) analysis were performed. The alterations obtained from the WES were classified into various signatures using the Copy Number Signatures tool available in COSMIC.

RESULTS

We identified several signatures that correlated with tumor stage and disctinct prognoses. These signatures all seem to be linked to replication stress, with CCNE1 amplification identified as the primary driver of oncogenesis in over 70% of tumors analyzed.

CONCLUSION

The identification of CCNE1 amplification, which is currently being explored as a therapeutic target in clinical trials, suggests new treatment strategies for high-grade endometrial cancer. This finding holds particular significance for Martinique, where access to care is challenging.

Txnrd1 as a prognosticator for recurrence, metastasis and response to neoadjuvant chemotherapy and radiotherapy in breast cancer patients

Thioredoxin reductase 1 (Txnrd1) is known to have prognostic significance in a subset of breast cancer patients. Despite the pivotal role of Txnrd1 in regulating several cellular and physiological processes in cancer progression and metastasis, its clinical significance is largely unrecognized. Here, we undertook a retrospective comprehensive meta-analysis of 13,322 breast cancer patients from 43 independent cohorts to assess prognostic and predictive roles of Txnrd1. We observed that Txnrd1 has a positive correlation with tumor grade and size and it is over-expressed in higher-grade and larger tumors. Further, hormone receptor-negative and HER2-positive tumors exhibit elevated Txnrd1 gene expression. Patients with elevated Txnrd1 expression exhibit significant hazards for shorter disease-specific and overall survival. While Txnrd1 has a positive correlation with tumor recurrence and metastasis, it has a negative correlation with time to recurrence and metastasis. Txnrd1High patients exhibit 2.5 years early recurrence and 1.3 years early metastasis as compared to Txnrd1Low cohort. Interestingly, patients with high Txnrd1 gene expression exhibit a pathologic complete response (pCR) to neoadjuvant chemotherapy, but they experience early recurrence after radiotherapy. Txnrd1High MDA-MB-231 cells exhibit significant ROS generation and reduced viability after doxorubicin treatment compared to Txnrd1Low MCF7 cells. Corroborating with findings from meta-analysis, Txnrd1 depletion leads to decreased survival, enhanced sensitivity to radiation induced killing, poor scratch-wound healing, and reduced invasion potential in MDA-MB-231 cells. Thus, Txnrd1 appears to be a potential predictor of recurrence, metastasis and therapy response in breast cancer patients.

Pan-Cancer Analysis of Copy-Number Features Identifies Recurrent Signatures and a Homologous Recombination Deficiency Biomarker to Predict Poly (ADP-Ribose) Polymerase Inhibitor Response

PURPOSE

Copy-number (CN) features reveal the molecular state of cancers and may have predictive and prognostic value in the treatment of cancer. We sought to apply published CN analysis methods to a large pan-cancer data set and characterize ubiquitous CN signatures across tumor types, including potential utility for treatment selection.

METHODS

We analyzed the landscape of CN features in 260,333 pan-cancer samples. We examined the association of 10 signatures with genomic alterations and clinical characteristics and trained a machine learning classifier using CN and insertion and deletion features to detect homologous recombination deficiency signature (HRDsig) positivity. Clinical outcomes were assessed using a real-world clinicogenomic database (CGDB) of comprehensive genomic profiling linked to deidentified, electronic health record-derived clinical data.

RESULTS

CN signatures were prevalent across cancer types and associated with diverse processes including focal tandem duplications, seismic amplifications, genome-wide loss of heterozygosity (gLOH), and HRD. Our novel HRDsig outperformed gLOH in predicting BRCAness and effectively distinguished biallelic BRCA and homologous recombination-repair wild-type (HRRwt) samples pan-tumor, demonstrating high sensitivity to detect biallelic BRCA in ovarian (93%) and other HRD-associated cancers (80%-87%). Pan-tumor prevalence of HRDsig was 6.4%. HRRwt cases represented a significant fraction of the HRDsig-positive cohort, likely reflecting a population with nongenomic mechanisms of HRD. In ovarian and prostate CGDBs, HRDsig identified more patients than gLOH and had predictive value for poly (ADP-ribose) polymerase inhibitor (PARPi) benefit.

CONCLUSION

Tumor CN profiles are informative, revealing diverse processes active in cancer. We describe the landscape of 10 CN signatures in a large pan-cancer cohort, including two associated with HRD. We trained a machine learning-based HRDsig that robustly identified BRCAness and associated with biallelic BRCA pan-tumor, and was predictive of PARPi benefit in real-world ovarian and prostate data sets.

Comparative Genomics Provides Etiologic and Biological Insight into Melanoma Subtypes

Melanoma is a cancer of melanocytes, with multiple subtypes based on body site location. Cutaneous melanoma is associated with skin exposed to ultraviolet radiation; uveal melanoma occurs in the eyes; mucosal melanoma occurs in internal mucous membranes; and acral melanoma occurs on the palms, soles, and nail beds. Here, we present the largest whole-genome sequencing study of melanoma to date, with 570 tumors profiled, as well as methylation and RNA sequencing for subsets of tumors. Uveal melanoma is genomically distinct from other melanoma subtypes, harboring the lowest tumor mutation burden and with significantly mutated genes in the G-protein signaling pathway. Most cutaneous, acral, and mucosal melanomas share alterations in components of the MAPK, PI3K, p53, p16, and telomere pathways. However, the mechanism by which these pathways are activated or inactivated varies between melanoma subtypes. Additionally, we identify potential novel germline predisposition genes for some of the less common melanoma subtypes.

SIGNIFICANCE

This is the largest whole-genome analysis of melanoma to date, comprehensively comparing the genomics of the four major melanoma subtypes. This study highlights both similarities and differences between the subtypes, providing insights into the etiology and biology of melanoma. This article is highlighted in the In This Issue feature, p. 2711.

Integrated Tissue and Blood miRNA Expression Profiles Identify Novel Biomarkers for Accurate Non-Invasive Diagnosis of Breast Cancer: Preliminary Results and Future Clinical Implications

We aimed to identify miRNAs that were closely related to breast cancer (BRCA). By integrating several methods including significance analysis of microarrays, fold change, Pearson’s correlation analysis, t test, and receiver operating characteristic analysis, we developed a decision-tree-based scoring algorithm, called Optimized Scoring Mechanism for Primary Synergy MicroRNAs (O-PSM). Five synergy miRNAs (hsa-miR-139-5p, hsa-miR-331-3p, hsa-miR-342-5p, hsa-miR-486-5p, and hsa-miR-654-3p) were identified using O-PSM, which were used to distinguish normal samples from pathological ones, and showed good results in blood data and in multiple sets of tissue data. These five miRNAs showed accurate categorization efficiency in BRCA typing and staging and had better categorization efficiency than experimentally verified miRNAs. In the Protein-Protein Interaction (PPI) network, the target genes of hsa-miR-342-5p have the most regulatory relationships, which regulate carcinogenesis proliferation and metastasis by regulating Glycosaminoglycan biosynthesis and the Rap1 signaling pathway. Moreover, hsa-miR-342-5p showed potential clinical application in survival analysis. We also used O-PSM to generate an R package uploaded on github (SuFei-lab/OPSM accessed on 22 October 2021). We believe that miRNAs included in O-PSM could have clinical implications for diagnosis, prognostic stratification and treatment of BRCA, proposing potential significant biomarkers that could be utilized to design personalized treatment plans in BRCA patients in the future.

Prognostic Value of Copy Number Alteration Burden in Early-Stage Breast Cancer and the Construction of an 11-Gene Copy Number Alteration Model

The increasing burden of breast cancer has prompted a wide range of researchers to search for new prognostic markers. Considering that tumor mutation burden (TMB) is low and copy number alteration burden (CNAB) is high in breast cancer, we built a CNAB-based model using a public database and validated it with a Chinese population. We collected formalin-fixed, paraffin-embedded (FFPE) tissue samples from 31 breast cancer patients who were treated between 2010 and 2014 at the National Cancer Center (CICAMS). METABRIC and TCGA data were downloaded via cBioPortal. In total, 2295 patients with early-stage breast cancer were enrolled in the study, including 1427 in the METABRIC cohort, 837 in the TCGA cohort, and 31 in the CICAMS cohort. Based on the ROC curve, we consider 2.2 CNA/MBp as the threshold for the CNAB-high and CNAB-low groupings. In both the TCGA cohort and the CICAMS cohort, CNAB-high had a worse prognosis than CNAB-low. We further simplified this model by establishing a prognostic nomogram for early breast cancer patients by 11 core genes, and this nomogram was highly effective in both the TCGA cohort and the CICAMS cohort. We hope that this model will subsequently help clinicians with prognostic assessments.

Whole genome copy number analyses reveal a highly aberrant genome in TP53 mutant lung adenocarcinoma tumors

BACKGROUND

Genetic alterations are common in non-small cell lung cancer (NSCLC), and DNA mutations and translocations are targets for therapy. Copy number aberrations occur frequently in NSCLC tumors and may influence gene expression and further alter signaling pathways. In this study we aimed to characterize the genomic architecture of NSCLC tumors and to identify genomic differences between tumors stratified by histology and mutation status. Furthermore, we sought to integrate DNA copy number data with mRNA expression to find genes with expression putatively regulated by copy number aberrations and the oncogenic pathways associated with these affected genes.

METHODS

Copy number data were obtained from 190 resected early-stage NSCLC tumors and gene expression data were available from 113 of the adenocarcinomas. Clinical and histopathological data were known, and EGFR-, KRAS- and TP53 mutation status was determined. Allele-specific copy number profiles were calculated using ASCAT, and regional copy number aberration were subsequently obtained and analyzed jointly with the gene expression data.

RESULTS

The NSCLC tumors tissue displayed overall complex DNA copy number profiles with numerous recurrent aberrations. Despite histological differences, tissue samples from squamous cell carcinomas and adenocarcinomas had remarkably similar copy number patterns. The TP53-mutated lung adenocarcinomas displayed a highly aberrant genome, with significantly altered copy number profiles including gains, losses and focal complex events. The EGFR-mutant lung adenocarcinomas had specific arm-wise aberrations particularly at chromosome7p and 9q. A large number of genes displayed correlation between copy number and expression level, and the PI(3)K-mTOR pathway was highly enriched for such genes.

CONCLUSIONS

The genomic architecture in NSCLC tumors is complex, and particularly TP53-mutated lung adenocarcinomas displayed highly aberrant copy number profiles. We suggest to always include TP53-mutation status when studying copy number aberrations in NSCLC tumors. Copy number may further impact gene expression and alter cellular signaling pathways.

The somatic molecular evolution of cancer: Mutation, selection, and epistasis

Cancer progression has been attributed to somatic changes in single-nucleotide variants, copy-number aberrations, loss of heterozygosity, chromosomal instability, epistatic interactions, and the tumor microenvironment. It is not entirely clear which of these changes are essential and which are ancillary to cancer. The dynamic nature of cancer evolution in a patient can be illuminated using several concepts and tools from classical evolutionary biology. Neutral mutation rates in cancer cells are calculable from genomic data such as synonymous mutations, and selective pressures are calculable from rates of fixation occurring beyond the expectation by neutral mutation and drift. However, these cancer effect sizes of mutations are complicated by epistatic interactions that can determine the likely sequence of gene mutations. In turn, longitudinal phylogenetic analyses of somatic cancer progression offer an opportunity to identify key moments in cancer evolution, relating the timing of driver mutations to corresponding landmarks in the clinical timeline. These analyses reveal temporal aspects of genetic and phenotypic change during tumorigenesis and across clinical timescales. Using a related framework, clonal deconvolution, physical locations of clones, and their phylogenetic relations can be used to infer tumor migration histories. Additionally, genetic interactions with the tumor microenvironment can be analyzed with longstanding approaches applied to organismal genotype-by-environment interactions. Fitness landscapes for cancer evolution relating to genotype, phenotype, and environment could enable more accurate, personalized therapeutic strategies. An understanding of the trajectories underlying the evolution of neoplasms, primary, and metastatic tumors promises fundamental advances toward accurate and personalized predictions of therapeutic response.

Patient-Specific Network for Personalized Breast Cancer Therapy with Multi-Omics Data

The development of new computational approaches that are able to design the correct personalized drugs is the crucial therapeutic issue in cancer research. However, tumor heterogeneity is the main obstacle to developing patient-specific single drugs or combinations of drugs that already exist in clinics. In this study, we developed a computational approach that integrates copy number alteration, gene expression, and a protein interaction network of 73 basal breast cancer samples. 2509 prognostic genes harboring a copy number alteration were identified using survival analysis, and a protein–protein interaction network considering the direct interactions was created. Each patient was described by a specific combination of seven altered hub proteins that fully characterize the 73 basal breast cancer patients. We suggested the optimal combination therapy for each patient considering drug–protein interactions. Our approach is able to confirm well-known cancer related genes and suggest novel potential drug target genes. In conclusion, we presented a new computational approach in breast cancer to deal with the intra-tumor heterogeneity towards personalized cancer therapy.

Akyerli C, K. Yüksel Ş, E. Danyeli A, Tihan T, Sezerman OU, Yakıcıer MC, Pamir MN, Özduman K. Mutations and Copy Number Alterations in IDH Wild-Type Glioblastomas Are Shaped by Different Oncogenic Mechanisms

Little is known about the mutational processes that shape the genetic landscape of gliomas. Numerous mutational processes leave marks on the genome in the form of mutations, copy number alterations, rearrangements or their combinations. To explore gliomagenesis, we hypothesized that gliomas with different underlying oncogenic mechanisms would have differences in the burden of various forms of these genomic alterations. This was an analysis on adult diffuse gliomas, but IDH-mutant gliomas as well as diffuse midline gliomas H3-K27M were excluded to search for the possible presence of new entities among the very heterogenous group of IDH-WT glioblastomas. The cohort was divided into two molecular subsets: (1) Molecularly-defined GBM (mGBM) as those that carried molecular features of glioblastomas (including TERT promoter mutations, 7/10 pattern, or EGFR-amplification), and (2) those who did not (others). Whole exome sequencing was performed for 37 primary tumors and matched blood samples as well as 8 recurrences. Single nucleotide variations (SNV), short insertion or deletions (indels) and copy number alterations (CNA) were quantified using 5 quantitative metrics (SNV burden, indel burden, copy number alteration frequency-wGII, chromosomal arm event ratio-CAER, copy number amplitude) as well as 4 parameters that explored underlying oncogenic mechanisms (chromothripsis, double minutes, microsatellite instability and mutational signatures). Findings were validated in the TCGA pan-glioma cohort. mGBM and "Others" differed significantly in their SNV (only in the TCGA cohort) and CNA metrics but not indel burden. SNV burden increased with increasing age at diagnosis and at recurrences and was driven by mismatch repair deficiency. On the contrary, indel and CNA metrics remained stable over increasing age at diagnosis and with recurrences. Copy number alteration frequency (wGII) correlated significantly with chromothripsis while CAER and CN amplitude correlated significantly with the presence of double minutes, suggesting separate underlying mechanisms for different forms of CNA.