Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout ageing.

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mitochondrial genome in different tissues throughout ageing. We used ultrasensitive duplex sequencing to profile ~2.5 million mitochondrial genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloguing ~1.2 million mitochondrial somatic and ultralow-frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the light strand origin of replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared with primates with a surfeit of reactive oxygen species-associated G > T/C > A mutations, and that somatic mutations in protein-coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that 're-align' mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

Error-corrected duplex sequencing enables direct detection and quantification of mutations in human TK6 cells with strong inter-laboratory consistency.

Error-corrected duplex sequencing (DS) enables direct quantification of low-frequency mutations and offers tremendous potential for chemical mutagenicity assessment. We investigated the utility of DS to quantify induced mutation frequency (MF) and spectrum in human lymphoblastoid TK6 cells exposed to a prototypical DNA alkylating agent, N-ethyl-N-nitrosourea (ENU). Furthermore, we explored appropriate experimental parameters for this application, and assessed inter-laboratory reproducibility. In two independent experiments in two laboratories, TK6 cells were exposed to ENU (25-200 µM) and DNA was sequenced 48, 72, and 96 h post-exposure. A DS mutagenicity panel targeting twenty 2.4-kb regions distributed across the genome was used to sample diverse, genome-representative sequence contexts. A significant increase in MF that was unaffected by time was observed in both laboratories. Concentration-response in the MF from the two laboratories was strongly positively correlated (r = 0.97). C:G>T:A, T:A>C:G, T:A>A:T, and T:A>G:C mutations increased in consistent, concentration-dependent manners in both laboratories, with high proportions of C:G>T:A at all time points. The consistent results across the three time points suggest that 48 h may be sufficient for mutation analysis post-exposure. The target sites responded similarly between the two laboratories and revealed a higher average MF in intergenic regions. These results, demonstrating remarkable reproducibility across time and laboratory for both MF and spectrum, support the high value of DS for characterizing chemical mutagenicity in both research and regulatory evaluation.

Next Generation Sequencing Workshop at the Royal Society of Medicine (London, May 2022): how genomics is on the path to modernizing genetic toxicology.

The use of error-corrected Next Generation Sequencing (ecNG) to determine mutagenicity has been a subject of growing interest and potentially a disruptive technology that could supplement, and in time, replace current testing paradigms in preclinical safety assessment. Considering this, a Next Generation Sequencing Workshop was held at the Royal Society of Medicine in London in May 2022, supported by the United Kingdom Environmental Mutagen Society (UKEMS) and TwinStrand Biosciences (WA, USA), to discuss progress and future applications of this technology. In this meeting report, the invited speakers provide an overview of the Workshop topics covered and identify future directions for research. In the area of somatic mutagenesis, several speakers reviewed recent progress made with correlating ecNGS to classic in vivo transgenic rodent mutation assays as well as exploring the use of this technology directly in humans and animals, and in complex organoid models. Additionally, ecNGS has been used for detecting off-target effects of gene editing tools and emerging data suggest ecNGS potential to measure clonal expansion of cells carrying mutations in cancer driver genes as an early marker of carcinogenic potential and for direct human biomonitoring. As such, the workshop demonstrated the importance of raising awareness and support for advancing the science of ecNGS for mutagenesis, gene editing, and carcinogenesis research. Furthermore, the potential of this new technology to contribute to advances in drug and product development and improve safety assessment was extensively explored.

Adopting Duplex Sequencing™ Technology for Genetic Toxicity Testing: A Proof-of-Concept Mutagenesis Experiment with N-Ethyl-N-Nitrosourea (ENU)-Exposed Rats.

Duplex sequencing (DuplexSeq) is an error-corrected next-generation sequencing (ecNGS) method in which molecular barcodes informatically link PCR-copies back to their source DNA strands, enabling computational removal of errors by comparing grouped strand sequencing reads. The resulting background of less than one artifactual mutation per 10 7 nucleotides allows for direct detection of somatic mutations. TwinStrand Biosciences, Inc. has developed a DuplexSeq-based mutagenesis assay to sample the rat genome, which can be applied to genetic toxicity testing. To evaluate this assay for early detection of mutagenesis, a time-course study was conducted using male Hsd:Sprague Dawley SD rats (3 per group) administered a single dose of 40 mg/kg N-ethyl-N-nitrosourea (ENU) via gavage, with mutation frequency (MF) and spectrum analyzed in stomach, bone marrow, blood, and liver tissues at 3 h, 24 h, 7 d, and 28 d post-exposure. Significant increases in MF were observed in ENU-exposed rats as early as 24 h for stomach (site of contact) and bone marrow (a highly proliferative tissue) and at 7 d for liver and blood. The canonical, mutational signature of ENU was established by 7 d post-exposure in all four tissues. Interlaboratory analysis of a subset of samples from different tissues and time points demonstrated remarkable reproducibility for both MF and spectrum. These results demonstrate that MF and spectrum can be evaluated successfully by directly sequencing targeted regions of DNA obtained from various tissues, a considerable advancement compared to currently used in vivo gene mutation assays. HIGHLIGHTS: DuplexSeq is an ultra-accurate NGS technology that directly quantifies mutationsENU-dependent mutagenesis was detected 24 h post-exposure in proliferative tissuesMultiple tissues exhibited the canonical ENU mutation spectrum 7 d after exposureResults obtained with DuplexSeq were highly concordant between laboratoriesThe Rat-50 Mutagenesis Assay is promising for applications in genetic toxicology.

Quantification of measurable residual disease using duplex sequencing in adults with acute myeloid leukemia.

The presence of measurable residual disease (MRD) is strongly associated with treatment outcomes in acute myeloid leukemia (AML). Despite the correlation with clinical outcomes, MRD assessment has yet to be standardized or routinely incorporated into clinical trials. Discrepancies have been observed between different techniques for MRD assessment and there remains a need to compare centralized, high-quality multiparametric flow cytometry (MFC) and ultrasensitive next-generation sequencing (NGS) in AML patients with diverse mutational profiles. In 62 patients with AML, aged 18-60, in first complete remission after intensive induction therapy on the randomized phase 3 SWOG-S0106 clinical trial, MRD detection by MFC was compared with a 29 gene panel utilizing duplex sequencing (DS), an NGS method that generates double-stranded consensus sequences to reduce false positive errors. Using DS, detection of a persistent mutation utilizing defined criteria was seen in 22 (35%) patients and was strongly associated with higher rates of relapse (68% vs 13% at year 5; HR, 8.8; 95% CI, 3.2-24.5; P<0.001) and decreased survival (32% vs 82% at year 5; HR, 5.6; 95% CI, 2.3-13.8; P<0.001). MRD as defined by DS strongly outperformed MFC, which was observed in 10 (16%) patients and marginally associated with higher rates of relapse (50% vs 30% at year 5; HR, 2.4; 95% CI, 0.9-6.7; P=0.087) and decreased survival (40% vs 68% at year 5; HR, 2.5; 95% CI, 1.0-6.3; P=0.059). Furthermore, the prognostic significance of DS MRD status at the time of remission was similar on both randomized arms of the trial, predicting S0106 clinical trial outcomes. These findings suggest that DS is a powerful tool that could be used in patient management and for early treatment assessment in clinical trials.

Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout aging.

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here, we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mt-genome in different tissues throughout aging. We used ultra-sensitive Duplex Sequencing to profile ~2.5 million mt-genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloging ~1.2 million mitochondrial somatic and ultra low frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the Light Strand Origin of Replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared to primates with a surfeit of reactive oxygen species-associated G>T/C>A mutations, and that somatic mutations in protein coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that "re-align" mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

Genetic toxicity testing using human in vitro organotypic airway cultures: Assessing DNA damage with the CometChip and mutagenesis by Duplex Sequencing.

The organotypic human air-liquid-interface (ALI) airway tissue model has been used as an in vitro cell culture system for evaluating the toxicity of inhaled substances. ALI airway cultures are highly differentiated, which has made it challenging to evaluate genetic toxicology endpoints. In the current study, we assayed DNA damage with the high-throughput CometChip assay and quantified mutagenesis with Duplex Sequencing, an error-corrected next-generation sequencing method capable of detecting a single mutation per 107 base pairs. Fully differentiated human ALI airway cultures were treated from the basolateral side with 6.25 to 100 µg/mL ethyl methanesulfonate (EMS) over a period of 28 days. CometChip assays were conducted after 3 and 28 days of treatment, and Duplex Sequencing after 28 days of treatment. Treating the airway cultures with EMS resulted in time- and concentration-dependent increases in DNA damage and a concentration-dependent increase in mutant frequency. The mutations observed in the EMS-treated cultures were predominantly C → T transitions and exhibited a unique trinucleotide signature relative to the negative control. Measurement of physiological endpoints indicated that the EMS treatments had no effect on anti-p63-positive basal cell frequency, but produced concentration-responsive increases in cytotoxicity and perturbations in cell morphology, along with concentration-responsive decreases in culture viability, goblet cell and anti-Ki67-positive proliferating cell frequency, cilia beating frequency, and mucin secretion. The results indicate that a unified 28-day study can be used to measure several important safety endpoints in physiologically relevant human in vitro ALI airway cultures, including DNA damage, mutagenicity, and tissue-specific general toxicity.

Direct quantification of in vivo mutagenesis and carcinogenesis using duplex sequencing.

The ability to accurately measure mutations is critical for basic research and identifying potential drug and chemical carcinogens. Current methods for in vivo quantification of mutagenesis are limited because they rely on transgenic rodent systems that are low-throughput, expensive, prolonged, and do not fully represent other species such as humans. Next-generation sequencing (NGS) is a conceptually attractive alternative for detecting mutations in the DNA of any organism; however, the limit of resolution for standard NGS is poor. Technical error rates (∼1 × 10-3) of NGS obscure the true abundance of somatic mutations, which can exist at per-nucleotide frequencies ≤1 × 10-7 Using duplex sequencing, an extremely accurate error-corrected NGS (ecNGS) technology, we were able to detect mutations induced by three carcinogens in five tissues of two strains of mice within 31 d following exposure. We observed a strong correlation between mutation induction measured by duplex sequencing and the gold-standard transgenic rodent mutation assay. We identified exposure-specific mutation spectra of each compound through trinucleotide patterns of base substitution. We observed variation in mutation susceptibility by genomic region, as well as by DNA strand. We also identified a primordial marker of carcinogenesis in a cancer-predisposed strain of mice, as evidenced by clonal expansions of cells carrying an activated oncogene, less than a month after carcinogen exposure. These findings demonstrate that ecNGS is a powerful method for sensitively detecting and characterizing mutagenesis and the early clonal evolutionary hallmarks of carcinogenesis. Duplex sequencing can be broadly applied to basic mutational research, regulatory safety testing, and emerging clinical applications.

Ultra-accurate Duplex Sequencing for the assessment of pretreatment ABL1 kinase domain mutations in Ph+ ALL.

Mutations of ABL1 are the dominant mechanism of relapse in Philadelphia chromosome-positive acute lymphoblastic leukemia (Ph + ALL). We performed highly accurate Duplex Sequencing of exons 4-10 of ABL1 on bone marrow or peripheral blood samples from 63 adult patients with previously untreated Ph + ALL who received induction with intensive chemotherapy plus a BCR-ABL1 TKI. We identified ABL1 mutations prior to BCR-ABL1 TKI exposure in 78% of patients. However, these mutations were generally present at extremely low levels (median variant allelic frequency 0.008% [range, 0.004%-3.71%] and did not clonally expand and lead to relapse in any patient, even when the pretreatment mutation was known to confer resistance to the TKI received. In relapse samples harboring a TKI-resistant ABL1 mutation, the corresponding mutation could not be detected pretreatment, despite validated sequencing sensitivity of Duplex Sequencing down to 0.005%. In samples under the selective pressure of ongoing TKI therapy, we detected low-level, emerging resistance mutations up to 5 months prior to relapse. These findings suggest that pretreatment ABL1 mutation assessment should not guide upfront TKI selection in Ph + ALL, although serial testing while on TKI therapy may allow for early detection of clinically actionable resistant clones.

Ultra-Sensitive TP53 Sequencing for Cancer Detection Reveals Progressive Clonal Selection in Normal Tissue over a Century of Human Lifespan.

High-accuracy next-generation DNA sequencing promises a paradigm shift in early cancer detection by enabling the identification of mutant cancer molecules in minimally invasive body fluid samples. We demonstrate 80% sensitivity for ovarian cancer detection using ultra-accurate Duplex Sequencing to identify TP53 mutations in uterine lavage. However, in addition to tumor DNA, we also detect low-frequency TP53 mutations in nearly all lavages from women with and without cancer. These mutations increase with age and share the selection traits of clonal TP53 mutations commonly found in human tumors. We show that low-frequency TP53 mutations exist in multiple healthy tissues, from newborn to centenarian, and progressively increase in abundance and pathogenicity with older age across tissue types. Our results illustrate that subclonal cancer evolutionary processes are a ubiquitous part of normal human aging, and great care must be taken to distinguish tumor-derived from age-associated mutations in high-sensitivity clinical cancer diagnostics.

Mutational spectra of aflatoxin B1 in vivo establish biomarkers of exposure for human hepatocellular carcinoma.

Aflatoxin B1 (AFB1) and/or hepatitis B and C viruses are risk factors for human hepatocellular carcinoma (HCC). Available evidence supports the interpretation that formation of AFB1-DNA adducts in hepatocytes seeds a population of mutations, mainly G:C→T:A, and viral processes synergize to accelerate tumorigenesis, perhaps via inflammation. Responding to a need for early-onset evidence predicting disease development, highly accurate duplex sequencing was used to monitor acquisition of high-resolution mutational spectra (HRMS) during the process of hepatocarcinogenesis. Four-day-old male mice were treated with AFB1 using a regimen that induced HCC within 72 wk. For analysis, livers were separated into tumor and adjacent cellular fractions. HRMS of cells surrounding the tumors revealed predominantly G:C→T:A mutations characteristic of AFB1 exposure. Importantly, 25% of all mutations were G→T in one trinucleotide context (CGC; the underlined G is the position of the mutation), which is also a hotspot mutation in human liver tumors whose incidence correlates with AFB1 exposure. The technology proved sufficiently sensitive that the same distinctive spectrum was detected as early as 10 wk after dosing, well before evidence of neoplasia. Additionally, analysis of tumor tissue revealed a more complex pattern than observed in surrounding hepatocytes; tumor HRMS were a composite of the 10-wk spectrum and a more heterogeneous set of mutations that emerged during tumor outgrowth. We propose that the 10-wk HRMS reflects a short-term mutational response to AFB1, and, as such, is an early detection metric for AFB1-induced liver cancer in this mouse model that will be a useful tool to reconstruct the molecular etiology of human hepatocarcinogenesis.