Evolutionary dynamics and significance of multiple subclonal mutations in cancer

Single-strand specific nuclease enhances accuracy of error-corrected sequencing and improves rare mutation-detection sensitivity

Error-corrected sequences (ECSs) that utilize double-stranded DNA sequences are useful in detecting mutagen-induced mutations. However, relatively higher frequencies of G:C > T:A (1 × 10⁻⁷ bp) and G:C > C:G (2 × 10⁻⁷ bp) errors decrease the accuracy of detection of rare G:C mutations (approximately 10⁻⁷ bp). Oxidized guanines in single-strand (SS) overhangs generated after shearing could serve as the source of these errors. To remove these errors, we first computationally discarded up to 20 read bases corresponding to the ends of the DNA fragments. Error frequencies decreased proportionately with trimming length; however, the results indicated that they were not sufficiently removed. To efficiently remove SS overhangs, we evaluated three mechanistically distinct SS-specific nucleases (S1 Nuclease, mung bean nuclease, and RecJf exonuclease) and found that they were more efficient than computational trimming. Consequently, we established Jade-Seq™, an ECS protocol with S1 Nuclease treatment, which reduced G:C > T:A and G:C > C:G errors to 0.50 × 10⁻⁷ bp and 0.12 × 10⁻⁷ bp, respectively. This was probably because S1 Nuclease removed SS regions, such as gaps and nicks, depending on its wide substrate specificity. Subsequently, we evaluated the mutation-detection sensitivity of Jade-Seq™ using DNA samples from TA100 cells exposed to 3-methylcholanthrene and 7,12-dimethylbenz[a]anthracene, which contained the rare G:C > T:A mutation (i.e., 2 × 10⁻⁷ bp). Fold changes of G:C > T:A compared to the vehicle control were 1.2- and 1.3-times higher than those of samples without S1 Nuclease treatment, respectively. These findings indicate the potential of Jade-Seq™ for detecting rare mutations and determining the mutagenicity of environmental mutagens.

How Should Cancer Models Be Constructed?

Choosing and optimizing treatment strategies for cancer requires capturing its complex dynamics sufficiently well for understanding but without being overwhelmed. Mathematical models are essential to achieve this understanding, and we discuss the challenge of choosing the right level of complexity to address the full range of tumor complexity from growth, the generation of tumor heterogeneity, and interactions within tumors and with treatments and the tumor microenvironment. We discuss the differences between conceptual and descriptive models, and compare the use of predator-prey models, evolutionary game theory, and dynamic precision medicine approaches in the face of uncertainty about mechanisms and parameter values. Although there is of course no one-size-fits-all approach, we conclude that broad and flexible thinking about cancer, based on combined modeling approaches, will play a key role in finding creative and improved treatments.

Molecular dynamics approach to identification of new OGG1 cancer-associated somatic variants with impaired activity

DNA of living cells is always exposed to damaging factors. To counteract the consequences of DNA lesions, cells have evolved several DNA repair systems, among which base excision repair is one of the most important systems. Many currently used antitumor drugs act by damaging DNA, and DNA repair often interferes with chemotherapy and radiotherapy in cancer cells. Tumors are usually extremely genetically heterogeneous, often bearing mutations in DNA repair genes. Thus, knowledge of the functionality of cancer-related variants of proteins involved in DNA damage response and repair is of great interest for personalization of cancer therapy. Although computational methods to predict the variant functionality have attracted much attention, at present, they are mostly based on sequence conservation and make little use of modern capabilities in computational analysis of 3D protein structures. We have used molecular dynamics (MD) to model the structures of 20 clinically observed variants of a DNA repair enzyme, 8-oxoguanine DNA glycosylase. In parallel, we have experimentally characterized the activity, thermostability, and DNA binding in a subset of these mutant proteins. Among the analyzed variants of 8-oxoguanine DNA glycosylase, three (I145M, G202C, and V267M) were significantly functionally impaired and were successfully predicted by MD. Alone or in combination with sequence-based methods, MD may be an important functional prediction tool for cancer-related protein variants of unknown significance.

Rate volatility and asymmetric segregation diversify mutation burden in cells with mutator alleles

Mutations that compromise mismatch repair (MMR) or DNA polymerase ε or δ exonuclease domains produce mutator phenotypes capable of fueling cancer evolution. Here, we investigate how combined defects in these pathways expands genetic heterogeneity in cells of the budding yeast, Saccharomyces cerevisiae, using a single-cell resolution approach that tallies all mutations arising from individual divisions. The distribution of replication errors present in mother cells after the initial S-phase was broader than expected for a single uniform mutation rate across all cell divisions, consistent with volatility of the mutator phenotype. The number of mismatches that then segregated to the mother and daughter cells co-varied, suggesting that each division is governed by a different underlying genome-wide mutation rate. The distribution of mutations that individual cells inherit after the second S-phase is further broadened by the sequential actions of semiconservative replication and mitotic segregation of chromosomes. Modeling suggests that this asymmetric segregation may diversify mutation burden in mutator-driven tumors.

Dowsett et al use a single-cell resolution approach to analyse the distribution of mutations across several divisions in yeast diploid strains mutated in mismatch repair and polymerase delta proofreading. They find that the underlying mutation rate varies from one division to another, and that new mutations segregate unequally between sister chromatids at each division, expanding genetic heterogeneity in the population.

Detection of genome-wide low-frequency mutations with Paired-End and Complementary Consensus Sequencing (PECC-Seq) revealed end-repair-derived artifacts as residual errors

To improve the accuracy and the cost-efficiency of next-generation sequencing in ultralow-frequency mutation detection, we developed the Paired-End and Complementary Consensus Sequencing (PECC-Seq), a PCR-free duplex consensus sequencing approach. PECC-Seq employed shear points as endogenous barcodes to identify consensus sequences from the overlap in the shortened, complementary DNA strand-derived paired-end reads for sequencing error correction. With the high accuracy of PECC-Seq, we identified the characteristic base substitution errors introduced by the end-repair process of mechanical fragmentation-based library preparations, which were prominent at the terminal 7 bp of the library fragments in the 5'-NpCpA-3' and 5'-NpCpT-3' trinucleotide context. As demonstrated at the human genome scale (TK6 cells), after removing these potential end-repair artifacts from the terminal 7 bp, PECC-Seq could reduce the sequencing error frequency to mid-10^-7 with a relatively low sequencing depth. For TA base pairs, the background error rate could be suppressed to mid-10^-8. In mutagen-treated (6 μg/mL methyl methanesulfonate or 12 μg/mL N-nitroso-N-ethylurea) TK6, increases in mutagen treatment-related mutant frequencies could be detected, indicating the potential of PECC-Seq in detecting genome-wide ultra-rare mutations. In addition, our finding on the patterns of end-repair artifacts may provide new insights into further reducing technical errors not only for PECC-Seq, but also for other next-generation sequencing techniques.

Evidence for immortality and autonomy in animal cancer models is often not provided, which causes confusion on key issues of cancer biology

Modern research into carcinogenesis has undergone three phases. Surgeons and pathologists started the first phase roughly 250 years ago, establishing morphological traits of tumors for pathologic diagnosis, and setting immortality and autonomy as indispensable criteria for neoplasms. A century ago, medical doctors, biologists and chemists started to enhance "experimental cancer research" by establishing many animal models of chemical-induced carcinogenesis for studies of cellular mechanisms. In this second phase, the two-hit theory and stepwise carcinogenesis of "initiation-promotion" or "initiation-promotion-progression" were established, with an illustrious finding that outgrowths induced in animals depend on the inducers, and thus are not authentically neoplastic, until late stages. The last 40 years are the third incarnation, molecular biologists have gradually dominated the carcinogenesis research fraternity and have established numerous genetically-modified animal models of carcinogenesis. However, evidence has not been provided for immortality and autonomy of the lesions from most of these models. Probably, many lesions had already been collected from animals for analyses of molecular mechanisms of "cancer" before the lesions became autonomous. We herein review the monumental work of many predecessors to reinforce that evidence for immortality and autonomy is essential for confirming a neoplastic nature. We extrapolate that immortality and autonomy are established early during sporadic human carcinogenesis, unlike the late establishment in most animal models. It is imperative to resume many forerunners' work by determining the genetic bases for initiation, promotion and progression, the genetic bases for immortality and autonomy, and which animal models are, in fact, good for identifying such genetic bases.

Comparative Study of NGS Platform Ion Torrent Personal Genome Machine and Therascreen Rotor-Gene Q for the Detection of Somatic Variants in Cancer

Molecular profiling of a tumor allows the opportunity to design specific therapies which are able to interact only with cancer cells characterized by the accumulation of several genomic aberrations. This study investigates the usefulness of next-generation sequencing (NGS) and mutation-specific analysis methods for the detection of target genes for current therapies in non-small-cell lung cancer (NSCLC), metastatic colorectal cancer (mCRC), and melanoma patients. We focused our attention on EGFR, BRAF, KRAS, and BRAF genes for NSCLC, melanoma, and mCRC samples, respectively. Our study demonstrated that in about 2% of analyzed cases, the two techniques did not show the same or overlapping results. Two patients affected by mCRC resulted in wild-type (WT) for BRAF and two cases with NSCLC were WT for EGFR according to PGM analysis. In contrast, these samples were mutated for the evaluated genes using the therascreen test on Rotor-Gene Q. In conclusion, our experience suggests that it would be appropriate to confirm the WT status of the genes of interest with a more sensitive analysis method to avoid the presence of a small neoplastic clone and drive the clinician to correct patient monitoring.

Mechanistic understanding of cellular responses to genomic stress

Within the past half century we have learned of multiple pathways for repairing damaged DNA, based upon the intrinsic redundancy of information in its complementary double strands. Mechanistic details of these pathways have provided insights into environmental and endogenous threats to genomic stability. Studies on bacterial responses to ultraviolet light led to the discovery of excision repair, as well as the inducible SOS response to DNA damage. Similar responses in eukaryotes promote upregulation of error-prone translesion DNA polymerases. Recent advances in this burgeoning field include duplex DNA sequencing to provide strikingly accurate profiling of mutational signatures, analyses of gene expression patterns in single cells, CRISPR/Cas9 to generate changes at precise genomic positions, novel roles for RNA in gene expression and DNA repair, phase-separated aqueous environments for specialized cellular transactions, and DNA lesions as epigenetic signals for gene expression. The Environmental Mutagenesis and Genomics Society (EMGS), through the broad range of expertise in its membership, stands at the crossroad of basic understanding of mechanisms for genomic maintenance and the field of genetic toxicology, with the need for regulation of exposures to toxic substances. Our future challenges include devising strategies and technologies to identify individuals who are susceptible to specific genomic stresses, along with basic research on the underlying mechanisms of cellular stress responses that promote disease-causing mutations. As the science moves forward it should also be a responsibility for the EMGS to expand its outreach programs for the enlightenment and benefit of all humans and the biosphere. Environ. Mol. Mutagen. 61:25-33, 2020. © 2019 Wiley Periodicals, Inc.

Treatment resistance in urothelial carcinoma: an evolutionary perspective

The emergence of treatment-resistant clones is a critical barrier to cure in patients with urothelial carcinoma. Setting the stage for the evolution of resistance, urothelial carcinoma is characterized by extensive mutational heterogeneity, which is detectable even in patients with early stage disease. Chemotherapy and immunotherapy both act as selective pressures that shape the evolutionary trajectory of urothelial carcinoma throughout the course of the disease. A detailed understanding of the dynamics of evolutionary drivers is required for the rational development of curative therapies. Herein, we describe the molecular basis of the clonal evolution of urothelial carcinomas and the use of genomic approaches to predict treatment responses. We discuss various mechanisms of resistance to chemotherapy with a focus on the mutagenic effects of the DNA dC->dU-editing enzymes APOBEC3 family of proteins. We also review the evolutionary mechanisms underlying resistance to immunotherapy, such as the loss of clonal tumour neoantigens. By dissecting treatment resistance through an evolutionary lens, the field will advance towards true precision medicine for urothelial carcinoma.

Potential roles of reactive oxygen species derived from chemical substances involved in cancer development in the female reproductive system

Reactive oxygen species (ROS) are major sources of cellular oxidative stress. Specifically, cancer cells harbor genetic alterations that promote a continuous and elevated production of ROS. While such oxidative stress conditions could be harmful to normal cells, they facilitate cancer cell growth in multiple ways by causing DNA damage and genomic instability, and ultimately by reprogramming cancer cell metabolism. This review provides up to date findings regarding the roles of ROS generation induced by diverse biological molecules and chemicals in representative women’s cancer. Specifically, we describe the cellular signaling pathways that regulate direct or indirect interactions between ROS homeostasis and metabolism within female genital cancer cells.

Multiclonal tumor origin: Evidence and implications

An accurate understanding of the clonal origins of tumors is critical for designing effective strategies to treat or prevent cancer and for guiding the field of cancer risk assessment. The intent of this review is to summarize evidence of multiclonal tumor origin and, thereby, contest the commonly held assumption of monoclonal tumor origin. This review describes relevant studies of X chromosome inactivation, analyses of tumor heterogeneity using other markers, single cell sequencing, and lineage tracing studies in aggregation chimeras and engineered rodent models. Methods for investigating tumor clonality have an inherent bias against detecting multiclonality. Despite this, multiclonality has been observed within all tumor stages and within 53 different types of tumors. For myeloid tumors, monoclonal tumor origin may be the predominant path to cancer and a monoclonal tumor origin cannot be ruled out for a fraction of other cancer types. Nevertheless, a large body of evidence supports the conclusion that most cancers are multiclonal in origin. Cooperation between different cell types and between clones of cells carrying different genetic and/or epigenetic lesions is discussed, along with how polyclonal tumor origin can be integrated with current perspectives on the genesis of tumors. In order to develop biologically sound and useful approaches to cancer risk assessment and precision medicine, mathematical models of carcinogenesis are needed, which incorporate multiclonal tumor origin and the contributions of spontaneous mutations in conjunction with the selective advantages conferred by particular mutations and combinations of mutations. Adherence to the idea that a growth must develop from a single progenitor cell to be considered neoplastic has outlived its usefulness. Moving forward, explicit examination of tumor clonality, using advanced tools, like lineage tracing models, will provide a strong foundation for future advances in clinical oncology and better training for the next generation of oncologists and pathologists.