Different Facets of Copy Number Changes: Permanent, Transient, and Adaptive

CONNECTING THE DOTS: EPIGENETIC REGULATION OF EXTRACHROMOSOMAL AND INHERITED DNA AMPLIFICATIONS

DNA amplification has intrigued scientists for decades. Since its discovery, significant progress has been made in understanding the mechanisms promoting DNA amplification and their associated function(s). While DNA copy gains were once thought to be regulated purely by stochastic processes, recent findings have revealed the important role of epigenetic modifications in driving these amplifications and their integration into the genome. Furthermore, advances in genomic technology have enabled detailed characterization of these genomic events in terms of size, structure, formation, and regulation. This review highlights how our understanding of DNA amplifications has evolved over time, tracing its trajectory from initial discovery to the contemporary landscape. We describe how recent discoveries have started to uncover how these genomic events occur by controlled biological processes rather than stochastic mechanisms, presenting opportunities for therapeutic modulation.

NF-κB signaling directs a program of transient amplifications at innate immune response genes

The cellular response to pathogens involves an intricate response directed by key innate immune signaling pathways which is characterized by cell-to-cell heterogeneity. How this heterogeneity is established and regulated remains unclear. We describe a program of transient site-specific gains (TSSG) producing extrachromosomal DNA (ecDNA) of immune-related genes in response to innate immune signaling. Activation of NF-κB drives TSSG of the interferon receptor gene cluster through inducible recruitment of the transcription factor RelA and the pre-replication complex member MCM2 to an epigenetically regulated TSSG control element. Targeted recruitment of RelA or p300 are sufficient to induce TSSG formation. RelA and MCM2 specify a program of TSSG for at least six and as many as 179 regions enriched in innate immune response genes. Identification of this program reveals regulated production of ecDNA as a mechanism of heterogeneity in the host response.

Analysis of clinical Candida parapsilosis isolates reveals copy number variation in key fluconazole resistance genes

We used whole-genome sequencing to analyze a collection of 35 fluconazole-resistant and 7 susceptible Candida parapsilosis isolates together with coverage analysis and GWAS techniques to identify new mechanisms of fluconazole resistance. Phylogenetic analysis shows that although the collection is diverse, two persistent clinical lineages were identified. We identified copy number variation (CNV) of two genes, ERG11 and CDR1B, in resistant isolates. Two strains have a CNV at the ERG11 locus; the entire ORF is amplified in one, and only the promoter region is amplified in the other. We show that the annotated telomeric gene CDR1B is actually an artifactual in silico fusion of two highly similar neighboring CDR genes due to an assembly error in the C. parapsilosis CDC317 reference genome. We report highly variable copy numbers of the CDR1B region across the collection. Several strains have increased the expansion of the two genes into a tandem array of new chimeric genes. Other strains have experienced a deletion between the two genes creating a single gene with a reciprocal chimerism. We find translocations, duplications, and gene conversion across the CDR gene family in the C. parapsilosis species complex, showing that it is a highly dynamic family.

Analysis of clinical Candida parapsilosis isolates reveals copy number variation in key fluconazole resistance genes

We used whole-genome sequencing to analyse a collection of 35 fluconazole resistant and 7 susceptible Candida parapsilosis isolates together with coverage analysis and GWAS techniques to identify new mechanisms of fluconazole resistance. Phylogenetic analysis shows that although the collection is diverse, two probable outbreak groups were identified. We identified copy number variation of two genes, ERG11 and CDR1B, in resistant isolates. Two strains have a CNV at the ERG11 locus; the entire ORF is amplified in one, and only the promoter region is amplified in the other. We show the annotated telomeric gene CDR1B is actually an artefactual in silico fusion of two highly similar neighbouring CDR genes due to an assembly error in the C. parapsilosis CDC317 reference genome. We report highly variable copy numbers of the CDR1B region across the collection. Several strains have increased expansion of the two genes into a tandem array of new chimeric genes. Other strains have experienced a deletion between the two genes creating a single gene with a reciprocal chimerism. We find translocations, duplications, and gene conversion across the CDR gene family in the C. parapsilosis species complex, showing that it is a highly dynamic family.

Assessment of MYC and TERT copy number variations in lung cancer using digital PCR

OBJECTIVE

Lung cancer is the second most frequent cancer type and the most common cause of cancer-related deaths worldwide. Alteration of gene copy numbers are associated with lung cancer and the determination of copy number variations (CNV) is appropriate for the discrimination between tumor and non-tumor tissue in lung cancer. As telomerase reverse transcriptase (TERT) and v-myc avian myelocytomatosis viral oncogene homolog (MYC) play a role in lung cancer the aims of this study were the verification of our recent results analyzing MYC CNV in tumor and non-tumor tissue of lung cancer patients using an independent study group and the assessment of TERT CNV as an additional marker.

RESULTS

TERT and MYC status was analyzed using digital PCR (dPCR) in tumor and adjacent non-tumor tissue samples of 114 lung cancer patients. The difference between tumor and non-tumor samples were statistically significant (p < 0.0001) for TERT and MYC. Using a predefined specificity of 99% a sensitivity of 41% and 51% was observed for TERT and MYC, respectively. For the combination of TERT and MYC the overall sensitivity increased to 60% at 99% specificity. We demonstrated that a combination of markers increases the performance in comparison to individual markers. Additionally, the determination of CNV using dPCR might be an appropriate tool in precision medicine.

Epigenetic balance ensures mechanistic control of MLL amplification and rearrangement

MLL/KMT2A amplifications and translocations are prevalent in infant, adult, and therapy-induced leukemia. However, the molecular contributor(s) to these alterations are unclear. Here, we demonstrate that histone H3 lysine 9 mono- and di-methylation (H3K9me1/2) balance at the MLL/KMT2A locus regulates these amplifications and rearrangements. This balance is controlled by the crosstalk between lysine demethylase KDM3B and methyltransferase G9a/EHMT2. KDM3B depletion increases H3K9me1/2 levels and reduces CTCF occupancy at the MLL/KMT2A locus, in turn promoting amplification and rearrangements. Depleting CTCF is also sufficient to generate these focal alterations. Furthermore, the chemotherapy doxorubicin (Dox), which associates with therapy-induced leukemia and promotes MLL/KMT2A amplifications and rearrangements, suppresses KDM3B and CTCF protein levels. KDM3B and CTCF overexpression rescues Dox-induced MLL/KMT2A alterations. G9a inhibition in human cells or mice also suppresses MLL/KMT2A events accompanying Dox treatment. Therefore, MLL/KMT2A amplifications and rearrangements are controlled by epigenetic regulators that are tractable drug targets, which has clinical implications.

Best Practices in Microbial Experimental Evolution: Using Reporters and Long-Read Sequencing to Identify Copy Number Variation in Experimental Evolution

Copy number variants (CNVs), comprising gene amplifications and deletions, are a pervasive class of heritable variation. CNVs play a key role in rapid adaptation in both natural, and experimental, evolution. However, despite the advent of new DNA sequencing technologies, detection and quantification of CNVs in heterogeneous populations has remained challenging. Here, we summarize recent advances in the use of CNV reporters that provide a facile means of quantifying de novo CNVs at a specific locus in the genome, and nanopore sequencing, for resolving the often complex structures of CNVs. We provide guidance for the engineering and analysis of CNV reporters and practical guidelines for single-cell analysis of CNVs using flow cytometry. We summarize recent advances in nanopore sequencing, discuss the utility of this technology, and provide guidance for the bioinformatic analysis of these data to define the molecular structure of CNVs. The combination of reporter systems for tracking and isolating CNV lineages and long-read DNA sequencing for characterizing CNV structures enables unprecedented resolution of the mechanisms by which CNVs are generated and their evolutionary dynamics.

Transition from Transient DNA Rereplication to Inherited Gene Amplification Following Prolonged Environmental Stress

Cells require the ability to adapt to changing environmental conditions, however, it is unclear how these changes elicit stable permanent changes in genomes. We demonstrate that, in response to environmental metal exposure, the metallothionein (MT) locus undergoes DNA rereplication generating transient site-specific gene amplifications (TSSGs). Chronic metal exposure allows transition from MT TSSG to inherited MT gene amplification through homologous recombination within and outside of the MT locus. DNA rereplication of the MT locus is suppressed by H3K27me3 and EZH2. Long-term ablation of EZH2 activity eventually leads to integration and inheritance of MT gene amplifications without the selective pressure of metal exposure. The rereplication and inheritance of MT gene amplification is an evolutionarily conserved response to environmental metal from yeast to human. Our results describe a new paradigm for adaptation to environmental stress where targeted, transient DNA rereplication precedes stable inherited gene amplification.

A pan-cancer analysis of ABI3BP: a potential biomarker for prognosis and immunoinfiltration

Objective

ABI Family Member 3 Binding Protein (ABI3BP) is an extracellular matrix protein that affects the carcinogenesis of lung and esophageal cancer. However, the relevance of ABI3BP in different forms of cancer is uncertain.

Methods

ABI3BP expression was interpreted using the Cancer Genome Atlas (TCGA) database, the Genotype Tissue Expression Atlas (GTEx) database, the Human Protein Atlas (HPA) database, the Cancer Cell Line Encyclopedia (CCLE) database, and immunohistochemistry. The R programming language was used to analyze the association between ABI3BP expression and patient prognosis and evaluate the relationship between ABI3BP and the immune characteristics of tumors. Using the GDSC and CTRP databases, a drug sensitivity analysis of ABI3BP was conducted.

Results

ABI3BP mRNA expression was shown by differential analysis to be down-regulated in 16 tumor types relative to normal tissues, corresponding with its protein expression level as determined by immunohistochemistry. Abnormal expression of ABI3BP accurately predicts the prognosis of patients with renal chromophobe carcinoma (KICH), mesothelioma (MESO), and pancreatic adenocarcinoma (PAAD). Meanwhile, aberrant expression of ABI3BP was associated with immune checkpoints, TMB, MSI, tumor purity, HRD, LOH, and drug sensitivity. A correlation between ABI3BP expression and the amount of infiltration of several immune-related cells in pan-cancer was determined by Immune Score, Stromal Score, and Estimated Score.

Conclusion

Our results show that ABI3BP might be employed as a molecular biomarker to predict prognosis, treatment susceptibility, and immunological response in patients with pan-cancer.

DNA repair as a shared hallmark in cancer and ageing

Increasing evidence demonstrates that DNA damage and genome instability play a crucial role in ageing. Mammalian cells have developed a wide range of complex and well‐orchestrated DNA repair pathways to respond to and resolve many different types of DNA lesions that occur from exogenous and endogenous sources. Defects in these repair pathways lead to accelerated or premature ageing syndromes and increase the likelihood of cancer development. Understanding the fundamental mechanisms of DNA repair will help develop novel strategies to treat ageing‐related diseases. Here, we revisit the processes involved in DNA damage repair and how these can contribute to diseases, including ageing and cancer. We also review recent mechanistic insights into DNA repair and discuss how these insights are being used to develop novel therapeutic strategies for treating human disease. We discuss the use of PARP inhibitors in the clinic for the treatment of breast and ovarian cancer and the challenges associated with acquired drug resistance. Finally, we discuss how DNA repair pathway‐targeted therapeutics are moving beyond PARP inhibition in the search for ever more innovative and efficacious cancer therapies.

The population genetics of adaptation through copy‐number variation in a fungal plant pathogen

Microbial pathogens can adapt rapidly to changing environments such as the application of pesticides or host resistance. Copy number variations (CNVs) are a major source of adaptive genetic variation for recent adaptation. Here, we analyse how a major fungal pathogen of barley, Rhynchosporium commune, has adapted to the host environment and fungicide applications. We screen the genomes of 125 isolates sampled across a worldwide set of populations and identify a total of 7,879 gene duplications and 116 gene deletions. Most gene duplications result from segmental chromosomal duplications. Although CNVs are generally under negative selection, we find that genes affected by CNVs are enriched in functions related to host exploitation (i.e., effectors and cell-wall-degrading enzymes). We perform genome-wide association studies (GWAS) and identify a large segmental duplication of CYP51A that has contributed to the emergence of azole resistance and a duplication encompassing an effector gene affecting virulence. We show that the adaptive CNVs were probably created by recently active transposable element families. Moreover, we find that specific transposable element families are important drivers of recent gene CNV. Finally, we use a genome-wide single nucleotide polymorphism data set to replicate the GWAS and contrast it with the CNV-focused analysis. Together, our findings show how extensive segmental duplications create the raw material for recent adaptation in global populations of a fungal pathogen.

Natural variation in the consequences of gene overexpression and its implications for evolutionary trajectories

Copy number variation through gene or chromosome amplification provides a route for rapid phenotypic variation and supports the long-term evolution of gene functions. Although the evolutionary importance of copy-number variation is known, little is understood about how genetic background influences its tolerance. Here, we measured fitness costs of over 4000 overexpressed genes in 15 Saccharomyces cerevisiae strains representing different lineages, to explore natural variation in tolerating gene overexpression (OE). Strain-specific effects dominated the fitness costs of gene OE. We report global differences in the consequences of gene OE, independent of the amplified gene, as well as gene-specific effects that were dependent on the genetic background. Natural variation in the response to gene OE could be explained by several models, including strain-specific physiological differences, resource limitations, and regulatory sensitivities. This work provides new insight on how genetic background influences tolerance to gene amplification and the evolutionary trajectories accessible to different backgrounds.

The cancer puzzle: Welcome to organicism

During the fifty years since President Nixon declared the "War on Cancer", those inside and outside the cancer community have witnessed the systematic moving of the goalposts attitude to accommodate evidence into an inadequate theory, that is, the Somatic Mutation Theory (SMT). This sorry state promoted a renewable yearly promise that at the end of the next 10-year period the promises uttered in 1971 would become reality. Each failure triggered calls to do more of the same research under the same theory, routinely using more and more sophisticated technology. Meanwhile, in the last few years, an unambiguous general consensus has emerged acknowledging that this overall long, intensive effort has failed, and that it is likely that the solution to the cancer problem resides elsewhere, namely, in alternative theoretical principles of biology. In this essay we concentrate, first, on the big picture, from the philosophical stance (reductionism versus organicism) to the need to adopt rigorous theories. From this novel perspective we conceptualize cancer as a disease of tissue organization akin to development gone awry. Finally, having identified both a promising stance and a useful theory, i.e., the tissue organization field theory (TOFT), we call for abandoning the SMT and for adopting the more promising TOFT.

Signatures of Discriminative Copy Number Aberrations in 31 Cancer Subtypes

Copy number aberrations (CNA) are one of the most important classes of genomic mutations related to oncogenetic effects. In the past three decades, a vast amount of CNA data has been generated by molecular-cytogenetic and genome sequencing based methods. While this data has been instrumental in the identification of cancer-related genes and promoted research into the relation between CNA and histo-pathologically defined cancer types, the heterogeneity of source data and derived CNV profiles pose great challenges for data integration and comparative analysis. Furthermore, a majority of existing studies have been focused on the association of CNA to pre-selected "driver" genes with limited application to rare drivers and other genomic elements. In this study, we developed a bioinformatics pipeline to integrate a collection of 44,988 high-quality CNA profiles of high diversity. Using a hybrid model of neural networks and attention algorithm, we generated the CNA signatures of 31 cancer subtypes, depicting the uniqueness of their respective CNA landscapes. Finally, we constructed a multi-label classifier to identify the cancer type and the organ of origin from copy number profiling data. The investigation of the signatures suggested common patterns, not only of physiologically related cancer types but also of clinico-pathologically distant cancer types such as different cancers originating from the neural crest. Further experiments of classification models confirmed the effectiveness of the signatures in distinguishing different cancer types and demonstrated their potential in tumor classification.

Copy Number Variations of CEP63, FOSL2 and PAQR6 Serve as Novel Signatures for the Prognosis of Bladder Cancer

Background

Finding effective prognostic signatures is of great urgency due to the high risk of recurrence and progression of bladder cancer (BC). Although a lot of genetic alterations are involved in the carcinogenesis, none of them were referred in the current risk group stratifications. In this study, we aimed to find significant copy number variations (CNVs) to predict prognosis for BC patients.

Methods

CNVs with high aberration frequencies in BC were explored by array-based comparative genomic hybridization in 65 tumor samples. Candidates were validated in independent groups of BC tumor samples (n=219) and urine samples (n=123). 3D digital PCR was applied for detecting accurate gene copy numbers in BC urine. In order to explore the prognostic value of candidate CNVs, all enrolled patients were followed up for the disease-free survival (DFS). Cox proportional hazards regression analysis was performed to find the independent prognostic factors for DFS.

Results

CNVs of CEP63, FOSL2 and PAQR6 with high aberration frequencies (67.7%, 56.9% and 60.0%, respectively) were found in BC tumors. Copy numbers of CEP63, FOSL2 and PAQR6 were gained in 219 tumor samples. CNVs of CEP63 and FOSL2 were correlated with advanced tumor stage and high grade. Retrospective analysis (median follow-up time: 69 months) revealed that CNVs of CEP63 and FOSL2 were independent prognostic factors for DFS of non-muscle-invasive bladder cancer (NMIBC) patients, while CNVs of FOSL2 and PAQR6 were independent prognostic factors for DFS of muscle-invasive bladder cancer (MIBC) patients. Models for predicting DFS were constructed based on CNVs of three genes. Patients with high prognostic indexes tended to have poor DFS. Prognostic index can also help to identify those with worse outcomes among high risk NMIBC patients. Copy number gains of CEP63 and FOSL2 in urine were found to be significantly correlated with poor DFS of NMIBC patients.

Conclusions

CNVs of CEP63, FOSL2 and PAQR6 were capable of predicting DFS and may serve as promising signatures for prognosis of BC.

Identification and utilization of copy number information for correcting Hi-C contact map of cancer cell lines

Background

Hi-C and its variant techniques have been developed to capture the spatial organization of chromatin. Normalization of Hi-C contact map is essential for accurate modeling and interpretation of high-throughput chromatin conformation capture (3C) experiments. Hi-C correction tools were originally developed to normalize systematic biases of karyotypically normal cell lines. However, a vast majority of available Hi-C datasets are derived from cancer cell lines that carry multi-level DNA copy number variations (CNVs). CNV regions display over- or under-representation of interaction frequencies compared to CN-neutral regions. Therefore, it is necessary to remove CNV-driven bias from chromatin interaction data of cancer cell lines to generate a euploid-equivalent contact map.

Results

We developed the HiCNAtra framework to compute high-resolution CNV profiles from Hi-C or 3C-seq data of cancer cell lines and to correct chromatin contact maps from systematic biases including CNV-associated bias. First, we introduce a novel ‘entire-fragment’ counting method for better estimation of the read depth (RD) signal from Hi-C reads that recapitulates the whole-genome sequencing (WGS)-derived coverage signal. Second, HiCNAtra employs a multimodal-based hierarchical CNV calling approach, which outperformed OneD and HiNT tools, to accurately identify CNVs of cancer cell lines. Third, incorporating CNV information with other systematic biases, HiCNAtra simultaneously estimates the contribution of each bias and explicitly corrects the interaction matrix using Poisson regression. HiCNAtra normalization abolishes CNV-induced artifacts from the contact map generating a heatmap with homogeneous signal. When benchmarked against OneD, CAIC, and ICE methods using MCF7 cancer cell line, HiCNAtra-corrected heatmap achieves the least 1D signal variation without deforming the inherent chromatin interaction signal. Additionally, HiCNAtra-corrected contact frequencies have minimum correlations with each of the systematic bias sources compared to OneD’s explicit method. Visual inspection of CNV profiles and contact maps of cancer cell lines reveals that HiCNAtra is the most robust Hi-C correction tool for ameliorating CNV-induced bias.

Conclusions

HiCNAtra is a Hi-C-based computational tool that provides an analytical and visualization framework for DNA copy number profiling and chromatin contact map correction of karyotypically abnormal cell lines. HiCNAtra is an open-source software implemented in MATLAB and is available at https://github.com/AISKhalil/HiCNAtra.

Integrative radiogenomics for virtual biopsy and treatment monitoring in ovarian cancer

BACKGROUND

Ovarian cancer survival rates have not changed in the last 20 years. The majority of cases are High-grade serous ovarian carcinomas (HGSOCs), which are typically diagnosed at an advanced stage with multiple metastatic lesions. Taking biopsies of all sites of disease is infeasible, which challenges the implementation of stratification tools based on molecular profiling.

MAIN BODY

In this review, we describe how these challenges might be overcome by integrating quantitative features extracted from medical imaging with the analysis of paired genomic profiles, a combined approach called radiogenomics, to generate virtual biopsies. Radiomic studies have been used to model different imaging phenotypes, and some radiomic signatures have been associated with paired molecular profiles to monitor spatiotemporal changes in the heterogeneity of tumours. We describe different strategies to integrate radiogenomic information in a global and local manner, the latter by targeted sampling of tumour habitats, defined as regions with distinct radiomic phenotypes.

CONCLUSION

Linking radiomics and biological correlates in a targeted manner could potentially improve the clinical management of ovarian cancer. Radiogenomic signatures could be used to monitor tumours during the course of therapy, offering additional information for clinical decision making. In summary, radiogenomics may pave the way to virtual biopsies and treatment monitoring tools for integrative tumour analysis.

A systematic literature review assessing if genetic biomarkers are predictors for platinum-based chemotherapy response in ovarian cancer patients

BACKGROUND

Ovarian cancer is the deadliest of gynecologic malignancies with the 5-year overall survival rate remaining at approximately 30%, a rate that has not improved over the last three decades. Standard of care for epithelial ovarian cancer patients consists of a platinum compound with a taxane given intravenously following debulking surgery; however, 80% of cases relapse within 2 years of diagnosis. This review sought to identify key underlying biomarkers related to platinum resistance in ovarian cancer to establish possible prognostic biomarkers of chemoresponse.

METHODS

A systematic literature review was conducted across three databases PubMed, EMBASE and SCOPUS to summarise the evidence for prognostic biomarkers in platinum-resistant ovarian cancer patients.

RESULTS

Forty-eight human studies were used in the review encompassing 6719 participants in retrospective and prospective study designs. A total of 68 biomarkers were reported that were significantly correlated with chemoresponse and/or survival reporting a p value less than or equal to 0.05.

CONCLUSION

This review accentuates the pleiotropic phenotypic complexities related to the response to platinum therapy in ovarian cancer. A one-size-fits-all approach may be ineffective in a large portion of patients, emphasising the need for a whole system-based approach and personalised treatment strategies. Identifying key biomarkers to aid clinical decision-making is the first essential step in developing and appropriating therapies for at-risk patients, reducing toxicity and improving quality of life.

Multi-omic signatures identify pan-cancer classes of tumors beyond tissue of origin

Despite recent advances in treatment, cancer continues to be one of the most lethal human maladies. One of the challenges of cancer treatment is the diversity among similar tumors that exhibit different clinical outcomes. Most of this variability comes from wide-spread molecular alterations that can be summarized by omic integration. Here, we have identified eight novel tumor groups (C1-8) via omic integration, characterized by unique cancer signatures and clinical characteristics. C3 had the best clinical outcomes, while C2 and C5 had poorest. C1, C7, and C8 were upregulated for cellular and mitochondrial translation, and relatively low proliferation. C6 and C4 were also downregulated for cellular and mitochondrial translation, and had high proliferation rates. C4 was represented by copy losses on chromosome 6, and had the highest number of metastatic samples. C8 was characterized by copy losses on chromosome 11, having also the lowest lymphocytic infiltration rate. C6 had the lowest natural killer infiltration rate and was represented by copy gains of genes in chromosome 11. C7 was represented by copy gains on chromosome 6, and had the highest upregulation in mitochondrial translation. We believe that, since molecularly alike tumors could respond similarly to treatment, our results could inform therapeutic action.

Hierarchical discovery of large-scale and focal copy number alterations in low-coverage cancer genomes

Background

Detection of DNA copy number alterations (CNAs) is critical to understand genetic diversity, genome evolution and pathological conditions such as cancer. Cancer genomes are plagued with widespread multi-level structural aberrations of chromosomes that pose challenges to discover CNAs of different length scales, and distinct biological origins and functions. Although several computational tools are available to identify CNAs using read depth (RD) signal, they fail to distinguish between large-scale and focal alterations due to inaccurate modeling of the RD signal of cancer genomes. Additionally, RD signal is affected by overdispersion-driven biases at low coverage, which significantly inflate false detection of CNA regions.

Results

We have developed CNAtra framework to hierarchically discover and classify ‘large-scale’ and ‘focal’ copy number gain/loss from a single whole-genome sequencing (WGS) sample. CNAtra first utilizes a multimodal-based distribution to estimate the copy number (CN) reference from the complex RD profile of the cancer genome. We implemented Savitzky-Golay smoothing filter and Modified Varri segmentation to capture the change points of the RD signal. We then developed a CN state-driven merging algorithm to identify the large segments with distinct copy numbers. Next, we identified focal alterations in each large segment using coverage-based thresholding to mitigate the adverse effects of signal variations. Using cancer cell lines and patient datasets, we confirmed CNAtra’s ability to detect and distinguish the segmental aneuploidies and focal alterations. We used realistic simulated data for benchmarking the performance of CNAtra against other single-sample detection tools, where we artificially introduced CNAs in the original cancer profiles. We found that CNAtra is superior in terms of precision, recall and f-measure. CNAtra shows the highest sensitivity of 93 and 97% for detecting large-scale and focal alterations respectively. Visual inspection of CNAs revealed that CNAtra is the most robust detection tool for low-coverage cancer data.

Conclusions

CNAtra is a single-sample CNA detection tool that provides an analytical and visualization framework for CNA profiling without relying on any reference control. It can detect chromosome-level segmental aneuploidies and high-confidence focal alterations, even from low-coverage data. CNAtra is an open-source software implemented in MATLAB^®. It is freely available at https://github.com/AISKhalil/CNAtra.

Over a century of cancer research: Inconvenient truths and promising leads

Despite over a century of intensive efforts, the great gains promised by the War on Cancer nearly 50 years ago have not materialized. Since 1999, we have analyzed the lack of progress in explaining and "curing" cancer by examining the merits of the premises that determine how cancer is understood and treated. Our ongoing critical analyses have aimed at clarifying the sources of misunderstandings at the root of the cancer puzzle while providing a plausible and comprehensive biomedical perspective as well as a new theory of carcinogenesis that is compatible with evolutionary theory. In this essay, we explain how this new theory, the tissue organization field theory (TOFT), can help chart a path to progress for cancer researchers by explaining features of cancer that remain unexplainable from the perspective of the still hegemonic somatic mutation theory (SMT) and its variants. Of equal significance, the premises underlying the TOFT offer new perspectives on basic biological phenomena.

Immuno-genomic landscape of osteosarcoma

Limited clinical activity has been seen in osteosarcoma (OS) patients treated with immune checkpoint inhibitors (ICI). To gain insights into the immunogenic potential of these tumors, we conducted whole genome, RNA, and T-cell receptor sequencing, immunohistochemistry and reverse phase protein array profiling (RPPA) on OS specimens from 48 pediatric and adult patients with primary, relapsed, and metastatic OS. Median immune infiltrate level was lower than in other tumor types where ICI are effective, with concomitant low T-cell receptor clonalities. Neoantigen expression in OS was lacking and significantly associated with high levels of nonsense-mediated decay (NMD). Samples with low immune infiltrate had higher number of deleted genes while those with high immune infiltrate expressed higher levels of adaptive resistance pathways. PARP2 expression levels were significantly negatively associated with the immune infiltrate. Together, these data reveal multiple immunosuppressive features of OS and suggest immunotherapeutic opportunities in OS patients.

The efficacy of immune checkpoint inhibitors (ICI) in osteosarcoma has been limited. Here, the authors investigate the immunogenomic landscape of osteosarcoma, and integrated analyses highlight features related to a suppressed immune microenvironment.

Histone Lysine Methylation Dynamics Control EGFR DNA Copy-Number Amplification

Acquired chromosomal DNA copy gains are a feature of many tumors; however, the mechanisms that underpin oncogene amplification are poorly understood. Recent studies have begun to uncover the importance of epigenetic states and histone lysine methyltransferases (KMT) and demethylases (KDM) in regulating transient site-specific DNA copy-number gains (TSSG). In this study, we reveal a critical interplay between a myriad of lysine methyltransferases and demethylases in modulating H3K4/9/27 methylation balance to control extrachromosomal amplification of the EGFR oncogene. This study further establishes that cellular signals (hypoxia and EGF) are able to directly promote EGFR amplification through modulation of the enzymes controlling EGFR copy gains. Moreover, we demonstrate that chemical inhibitors targeting specific KMTs and KDMs are able to promote or block extrachromosomal EGFR amplification, which identifies potential therapeutic strategies for controlling EGFR copy-number heterogeneity in cancer, and, in turn, drug response. SIGNIFICANCE: This study identifies a network of epigenetic factors and cellular signals that directly control EGFR DNA amplification. We demonstrate that chemical inhibitors targeting enzymes controlling this amplification can be used to rheostat EGFR copy number, which uncovers therapeutic opportunities for controlling EGFR DNA amplification heterogeneity and the associated drug response.This article is highlighted in the In This Issue feature, p. 161.

DNA Copy Number Variations as Markers of Mutagenic Impact

DNA copy number variation (CNV) occurs due to deletion or duplication of DNA segments resulting in a different number of copies of a specific DNA-stretch on homologous chromosomes. Implications of CNVs in evolution and development of different diseases have been demonstrated although contribution of environmental factors, such as mutagens, in the origin of CNVs, is poorly understood. In this review, we summarize current knowledge about mutagen-induced CNVs in human, animal and plant cells. Differences in CNV frequencies induced by radiation and chemical mutagens, distribution of CNVs in the genome, as well as adaptive effects in plants, are discussed. Currently available information concerning impact of mutagens in induction of CNVs in germ cells is presented. Moreover, the potential of CNVs as a new endpoint in mutagenicity test-systems is discussed.

MEK1/2 inhibitor withdrawal reverses acquired resistance driven by BRAFV600E amplification whereas KRASG13D amplification promotes EMT-chemoresistance

Acquired resistance to MEK1/2 inhibitors (MEKi) arises through amplification of BRAFV600E or KRASG13D to reinstate ERK1/2 signalling. Here we show that BRAFV600E amplification and MEKi resistance are reversible following drug withdrawal. Cells with BRAFV600E amplification are addicted to MEKi to maintain a precise level of ERK1/2 signalling that is optimal for cell proliferation and survival, and tumour growth in vivo. Robust ERK1/2 activation following MEKi withdrawal drives a p57KIP2-dependent G1 cell cycle arrest and senescence or expression of NOXA and cell death, selecting against those cells with amplified BRAFV600E. p57KIP2 expression is required for loss of BRAFV600E amplification and reversal of MEKi resistance. Thus, BRAFV600E amplification confers a selective disadvantage during drug withdrawal, validating intermittent dosing to forestall resistance. In contrast, resistance driven by KRASG13D amplification is not reversible; rather ERK1/2 hyperactivation drives ZEB1-dependent epithelial-to-mesenchymal transition and chemoresistance, arguing strongly against the use of drug holidays in cases of KRASG13D amplification.

Cross-talk between Lysine-Modifying Enzymes Controls Site-Specific DNA Amplifications

Acquired chromosomal DNA amplifications are features of many tumors. Although overexpression and stabilization of the histone H3 lysine 9/36 (H3K9/36) tri-demethylase KDM4A generates transient site-specific copy number gains (TSSGs), additional mechanisms directly controlling site-specific DNA copy gains are not well defined. In this study, we uncover a collection of H3K4-modifying chromatin regulators that function with H3K9 and H3K36 regulators to orchestrate TSSGs. Specifically, the H3K4 tri-demethylase KDM5A and specific COMPASS/KMT2 H3K4 methyltransferases modulate different TSSG loci through H3K4 methylation states and KDM4A recruitment. Furthermore, a distinct chromatin modifier network, MLL1-KDM4B-KDM5B, controls copy number regulation at a specific genomic locus in a KDM4A-independent manner. These pathways comprise an epigenetic addressing system for defining site-specific DNA rereplication and amplifications.

Implementing precision cancer medicine in the genomic era

The utilization of genomic data to direct treatment for cancer patients represents the central tenet in precision oncology, in which a patient is matched to a specific drug or therapy based on the genetic drivers detected in his or her tumor rather than the tumor's histologic classification. The expected but not always realized outcomes of molecularly matched therapies include increased response rates, more durable responses, deeper responses, and decreased number of therapy-related side effects. In this review, we will discuss different facets of utilizing genomic data to direct the increasingly complex care of cancer patients. We discuss the enlarging compendium of actionable genomic alterations and the development of novel molecular diagnostic assays for clinical application. Finally, we present an overview of the growing number of genomics-driven clinical trials and conclude with a discussion of future challenges in the implementation of precision oncology.

Genome plasticity in response to stress in Tetrahymena thermophila: selective and reversible chromosome amplification and paralogous expansion of metallothionein genes

Extreme stress situations can induce genetic variations including genome reorganization. In ciliates like Tetrahymena thermophila, the approximately 45-fold ploidy of the somatic macronucleus may enable adaptive responses that depend on genome plasticity. To identify potential genome-level adaptations related to metal toxicity, we isolated three Tetrahymena thermophila strains after an extended adaptation period to extreme metal concentrations (Cd²⁺ , Cu²⁺ or Pb²⁺ ). In the Cd-adapted strain, we found a approximately five-fold copy number increase of three genes located in the same macronuclear chromosome, including two metallothionein genes, MTT1 and MTT3. The apparent amplification of this macronuclear chromosome was reversible and reproducible, depending on the presence of environmental metal. We also analysed three knockout (KO) and/or knockdown (KD) strains for MTT1 and/or MTT5. In the MTT5KD strain, we found at least two new genes arising from paralogous expansion of MTT1, which encode truncated variants of MTT1. The expansion can be explained by a model based on somatic recombination between MTT1 genes on pairs of macronuclear chromosomes. At least two of the new paralogs are transcribed and upregulated in response to Cd²⁺ . Altogether, we have thus identified two distinct mechanisms, both involving genomic plasticity in the polyploid macronucleus that may represent adaptive responses to metal-related stress.

Overlooked: Extrachromosomal DNA and Their Possible Impact on Whole Genome Sequencing

Extrachromosomal (ec) DNA in eukaryotic cells has been known for decades. The structures described range from linear double stranded (ds) DNA to circular dsDNA, distinct from mitochondrial (mt) DNA. The sizes of circular forms are described from some hundred base pairs (bp) up to more than 150 kbp. The number of molecules per cell ranges from several hundred to a thousand. Semi-quantitative determinations of circular dsDNA show proportions as high as several percentages of the total DNA per cell. These ecDNA fractions harbor sequences that are known to be present in chromosomal DNA (chrDNA) too. Sequencing projects on, for example the human genome, have to take into account the ecDNA sequences which are simultaneously ascertained; corrections cannot be performed retrospectively. Concerning the results of sequencings derived from extracted whole DNA: if the ecDNA fractions contained therein are not taken into account, erroneous conclusions at the chromosomal level may result.

Site-Specific Genomic Alterations in a Well-Differentiated Pancreatic Neuroendocrine Tumor With High-Grade Progression

The major categories of pancreatic neuroendocrine tumor (PanNET) are well-differentiated NET and poorly differentiated neuroendocrine carcinoma. Sequencing of these tumors has identified multiple important genes in the pathogenesis of PanNETs, such as DAXX/ATRX, MEN1, TP53, RB, and mTOR pathway genes. We identified a case of well-differentiated PanNET with high-grade progression with simultaneous low- and high-grade histologic regions containing variable genomic profiles. We performed tumor microdissection and analyzed both regions using a 409-gene comprehensive cancer panel using next-generation sequencing in addition to immunohistochemical and morphologic studies. The low-grade region showed a change in the DAXX gene as a copy number variant (CNV) deletion. The high-grade region showed CNV deletion changes in the DAXX gene as well as the MEN1 gene. We observed additional mutational changes in the PTEN gene and SMAD4 gene in the high-grade region. Our data support that high-grade progression in PanNETs may be the result of the progressive accumulation of genetic changes (CNVs and point mutational changes) within the body of the tumor. Next generation sequencing may provide pathologists and clinicians with ancillary information to accurately characterize and treat these tumors.

Glyphosate Resistance and EPSPS Gene Duplication: Convergent Evolution in Multiple Plant Species

One of the increasingly widespread mechanisms of resistance to the herbicide glyphosate is copy number variation (CNV) of the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene. EPSPS gene duplication has been reported in 8 weed species, ranging from 3 to 5 extra copies to more than 150 extra copies. In the case of Palmer amaranth (Amaranthus palmeri), a section of >300 kb containing EPSPS and many other genes has been replicated and inserted at new loci throughout the genome, resulting in significant increase in total genome size. The replicated sequence contains several classes of mobile genetic elements including helitrons, raising the intriguing possibility of extra-chromosomal replication of the EPSPS-containing sequence. In kochia (Kochia scoparia), from 3 to more than 10 extra EPSPS copies are arranged as a tandem gene duplication at one locus. In the remaining 6 weed species that exhibit EPSPS gene duplication, little is known about the underlying mechanisms of gene duplication or their entire sequence. There is mounting evidence that adaptive gene amplification is an important mode of evolution in the face of intense human-mediated selection pressure. The convergent evolution of CNVs for glyphosate resistance in weeds, through at least 2 different mechanisms, may be indicative of a more general importance for this mechanism of adaptation in plants. CNVs warrant further investigation across plant functional genomics for adaptation to biotic and abiotic stresses, particularly for adaptive evolution on rapid time scales.

SNP co-association and network analyses identify E2F3, KDM5A and BACH2 as key regulators of the bovine milk fatty acid profile

The fatty acid (FA) profile has a considerable impact on the nutritional and technological quality of milk and dairy products. The molecular mechanism underlying the regulation of fat metabolism in bovine mammary gland have been not completely elucidated. We conducted genome-wide association studies (GWAS) across 65 milk FAs and fat percentage in 1,152 Brown Swiss cows. In total, we identified 175 significant single nucleotide polymorphism (SNPs) spanning all chromosomes. Pathway analyses revealed that 12:0 was associated with the greatest number of overrepresented categories/pathways (e.g. mitogen-activated protein kinase (MAPK) activity and protein phosphorylation), suggesting that it might play an important biological role in controlling milk fat composition. An Associated Weight Matrix approach based on SNP co-associations predicted a network of 791 genes related to the milk FA profile, which were involved in several connected molecular pathways (e.g., MAPK, lipid metabolism and hormone signalling) and undetectable through standard GWAS. Analysis of transcription factors and their putative target genes within the network identified BACH2, E2F3 and KDM5A as key regulators of milk FA metabolism. These findings contribute to increasing knowledge of FA metabolism and mammary gland functionality in dairy cows and may be useful in developing targeted breeding practices to improve milk quality.

The Influence of Polyploidy on the Evolution of Yeast Grown in a Sub-Optimal Carbon Source

Polyploidization events have occurred during the evolution of many fungi, plant, and animal species and are thought to contribute to speciation and tumorigenesis, however little is known about how ploidy level contributes to adaptation at the molecular level. Here we integrate whole genome sequencing, RNA expression analysis, and relative fitness of ∼100 evolved clones at three ploidy levels. Independent haploid, diploid, and tetraploid populations were grown in a low carbon environment for 250 generations. We demonstrate that the key adaptive mutation in the evolved clones is predicted by a gene expression signature of just five genes. All of the adaptive mutations identified encompass a narrow set of genes, however the tetraploid clones gain a broader spectrum of adaptive mutations than haploid or diploid clones. While many of the adaptive mutations occur in genes that encode proteins with known roles in glucose sensing and transport, we discover mutations in genes with no canonical role in carbon utilization (IPT1 and MOT3), as well as identify novel dominant mutations in glucose signal transducers thought to only accumulate recessive mutations in carbon limited environments (MTH1 and RGT1). We conclude that polyploid cells explore more genotypic and phenotypic space than lower ploidy cells. Our study provides strong evidence for the beneficial role of polyploidization events that occur during the evolution of many species and during tumorigenesis.

A Gene Signature for Selecting Benefit from Hypoxia Modification of Radiotherapy for High-Risk Bladder Cancer Patients

Purpose: Hypoxia modification improves overall survival in muscle-invasive bladder cancer patients who undergo radiotherapy. There is evidence that hypoxic tumors benefit most from hypoxia modification. The study aimed to identify or derive a hypoxia gene signature that predicts benefit from hypoxia-modifying treatment in bladder cancer.Experimental Design: Published hypoxia signatures were tested and a new one derived by analyzing bladder cancer transcriptomic data from public databases. Tumor samples were available from the BCON phase III randomized trial of radiotherapy alone or with carbogen and nicotinamide (CON). Gene expression data were generated for 151 tumors using Affymetrix Human 1.0 Exon ST arrays and used for independent validation.Results: A 24-gene signature was derived, which was prognostic in four of six independent surgical cohorts (n = 679; meta HR, 2.32; 95% CI, 1.73-3.12; P < 0.0001). The signature was also prognostic in BCON patients receiving radiotherapy alone (n = 75; HR for local relapse-free survival, 2.37; 95% CI, 1.26-4.47; P = 0.0076). The signature predicted benefit from CON (n = 76; HR, 0.47; 95% CI, 0.26-0.86; P = 0.015). Prognostic significance (P = 0.017) and predictive significance (P = 0.058) remained after adjusting for clinicopathologic variables. A test for interaction between hypoxia status and treatment arms was significant (P = 0.0094).Conclusions: A 24-gene hypoxia signature has strong and independent prognostic and predictive value for muscle-invasive bladder cancer patients. The signature can aid identification of patients likely to benefit from the addition of carbogen and nicotinamide to radiotherapy. Clin Cancer Res; 23(16); 4761-8. ©2017 AACR.

Extrachromosomal oncogene amplification drives tumour evolution and genetic heterogeneity

Human cells have twenty-three pairs of chromosomes. In cancer, however, genes can be amplified in chromosomes or in circular extrachromosomal DNA (ecDNA), although the frequency and functional importance of ecDNA are not understood. We performed whole-genome sequencing, structural modelling and cytogenetic analyses of 17 different cancer types, including analysis of the structure and function of chromosomes during metaphase of 2,572 dividing cells, and developed a software package called ECdetect to conduct unbiased, integrated ecDNA detection and analysis. Here we show that ecDNA was found in nearly half of human cancers; its frequency varied by tumour type, but it was almost never found in normal cells. Driver oncogenes were amplified most commonly in ecDNA, thereby increasing transcript level. Mathematical modelling predicted that ecDNA amplification would increase oncogene copy number and intratumoural heterogeneity more effectively than chromosomal amplification. We validated these predictions by quantitative analyses of cancer samples. The results presented here suggest that ecDNA contributes to accelerated evolution in cancer.

Regulation of Transient Site-specific Copy Gain by MicroRNA

Intra-tumor copy number heterogeneity is commonly observed in cancer; however, the molecular mechanisms that contribute to heterogeneity remain poorly understood. Up-regulation of the histone demethylase KDM4A promotes transient site-specific copy gain (TSSG) in cells; therefore, uncovering how KDM4A levels are controlled is important for understanding the regulation of copy number heterogeneity. Here, we demonstrate that KDM4A is regulated by hsa-mir-23a-3p, hsa-mir-23b-3p, and hsa-mir-137. Altering expression of these microRNAs (miRNAs) regulates KDM4A-dependent TSSG. miRNA inhibition promoted copy gains and increased expression of the drug-resistant oncogene CKS1B, which was further substantiated in primary breast tumors. Consistent with increased CKS1B expression, miRNA inhibition reduced breast cancer cell sensitivity to cisplatin. Our data identify these miRNAs as regulators of TSSG and copy gains of a drug resistance gene.