Direct quantification of in vivo mutagenesis and carcinogenesis using duplex sequencing

Estimation of Acceptable Daily Intake Values based on Modeling and In Vivo Mutagenicity of NDSRIs of Fluoxetine, Duloxetine and Atomoxetine

Nitrosamine drug substance related impurities or NDSRIs can be formed if an active pharmaceutical ingredient (API) has an intrinsic secondary amine that can undergo nitrosation. This is a concern as 1) nitrosamines are potentially highly potent carcinogens, 2) secondary amines in API are common, and 3) NDSRIs that might form from such secondary amines will be of unknown carcinogenic potency. Approaches for evaluating NDSRIs include read across, quantum mechanical modeling of reactivity, in vitro mutation data, and transgenic in vivo mutation data. These approaches were used here to assess NDSRIs that could potentially form from the drugs fluoxetine, duloxetine and atomoxetine. Based on a read across informed by modeling of physicochemical properties and mechanistic activation from quantum mechanical modeling, NDSRIs of fluoxetine, duloxetine, and atomoxetine were 10-100-fold less potent compared with highly potent nitrosamines such as NDMA or NDEA. While the NDSRIs were all confirmed to be mutagenic in vitro (Ames assay) and in vivo (TGR) studies, the latter data indicated that the potency of the mutation response was > 4400 ng/day for all compounds- an order of magnitude higher than published regulatory limits for these NDSRIs. The approaches described herein can be used qualitatively to better categorize NDSRIs with respect to potency and inform whether they are in the ICH M7R2-designated Cohort of Concern.

Human cancer genomes harbor the mutational signature of tobacco-specific nitrosamines NNN and NNK

Tobacco usage is linked to multiple cancer types and accounts for a quarter of all cancer-related deaths. Tobacco smoke contains various carcinogenic compounds, including polycyclic aromatic hydrocarbons (PAH), though the mutagenic potential of many tobacco-related chemicals remains largely unexplored. In particular, the highly carcinogenic tobacco-specific nitrosamines NNN and NNK form pre-mutagenic pyridyloxobutyl (POB) DNA adducts. In the study presented here, we identified genome-scale POB-induced mutational signatures in cell lines and rat tumors, while also investigating their role in human cancer. These signatures are characterized by T>N and C>T mutations forming from specific POB adducts damaging dT and dC residues. Analysis of 2,780 cancer genomes uncovered POB signatures in ∼180 tumors; from cancer types distinct from the ones linked to smoking-related signatures SBS4 and SBS92. This suggests that, unlike PAH compounds, the POB pathway may contribute uniquely to the mutational landscapes of certain hematological malignancies and cancers of the kidney, breast, prostate and pancreas.

Severity of effect considerations regarding the use of mutation as a toxicological endpoint for risk assessment: A report from the 8th International Workshop on Genotoxicity Testing (IWGT)

Exposure levels without appreciable human health risk may be determined by dividing a point of departure on a dose-response curve (e.g., benchmark dose) by a composite adjustment factor (AF). An "effect severity" AF (ESAF) is employed in some regulatory contexts. An ESAF of 10 may be incorporated in the derivation of a health-based guidance value (HBGV) when a "severe" toxicological endpoint, such as teratogenicity, irreversible reproductive effects, neurotoxicity, or cancer was observed in the reference study. Although mutation data have been used historically for hazard identification, this endpoint is suitable for quantitative dose-response modeling and risk assessment. As part of the 8th International Workshops on Genotoxicity Testing, a sub-group of the Quantitative Analysis Work Group (WG) explored how the concept of effect severity could be applied to mutation. To approach this question, the WG reviewed the prevailing regulatory guidance on how an ESAF is incorporated into risk assessments, evaluated current knowledge of associations between germline or somatic mutation and severe disease risk, and mined available data on the fraction of human germline mutations expected to cause severe disease. Based on this review and given that mutations are irreversible and some cause severe human disease, in regulatory settings where an ESAF is used, a majority of the WG recommends applying an ESAF value between 2 and 10 when deriving a HBGV from mutation data. This recommendation may need to be revisited in the future if direct measurement of disease-causing mutations by error-corrected next generation sequencing clarifies selection of ESAF values.

Origins of Second Malignancies in Children and Mutational Footprint of Chemotherapy in Normal Tissues

Pediatric cancers are rare diseases, and children without known germline predisposing conditions who develop a second malignancy during developmental ages are extremely rare. We present four such clinical cases and, through whole-genome and error-correcting ultra-deep duplex sequencing of tumor and normal samples, we explored the origin of the second malignancy in four children, uncovering different routes of development. The exposure to cytotoxic therapies was linked to the emergence of a secondary acute myeloid leukemia. A common somatic mutation acquired early during embryonic development was the driver of two solid malignancies in another child. In two cases, the two tumors developed from completely independent clones diverging during embryogenesis. Importantly, we demonstrate that platinum-based therapies contributed at least one order of magnitude more mutations per day of exposure than aging to normal tissues in these children.

SIGNIFICANCE

Using whole-genome and error-correcting ultra-deep duplex sequencing, we uncover different origins for second neoplasms in four children. We also uncover the presence of platinum-related mutations across 10 normal tissues of exposed individuals, highlighting the impact that the use of cytotoxic therapies may have on cancer survivors. See related commentary by Pacyna and Nangalia, p. 900. This article is featured in Selected Articles from This Issue, p. 897.

Effects of urban-induced mutations on ecology, evolution and health

Increasing evidence suggests that urbanization is associated with higher mutation rates, which can affect the health and evolution of organisms that inhabit cities. Elevated pollution levels in urban areas can induce DNA damage, leading to de novo mutations. Studies on mutations induced by urban pollution are most prevalent in humans and microorganisms, whereas studies of non-human eukaryotes are rare, even though increased mutation rates have the potential to affect organisms and their populations in contemporary time. Our Perspective explores how higher mutation rates in urban environments could impact the fitness, ecology and evolution of populations. Most mutations will be neutral or deleterious, and higher mutation rates associated with elevated pollution in urban populations can increase the risk of cancer in humans and potentially other species. We highlight the potential for urban-driven increased deleterious mutational loads in some organisms, which could lead to a decline in population growth of a wide diversity of organisms. Although beneficial mutations are expected to be rare, we argue that higher mutation rates in urban areas could influence adaptive evolution, especially in organisms with short generation times. Finally, we explore avenues for future research to better understand the effects of urban-induced mutations on the fitness, ecology and evolution of city-dwelling organisms.

ICH S1 prospective evaluation study and weight of evidence assessments: commentary from industry representatives

Industry representatives on the ICH S1B(R1) Expert Working Group (EWG) worked closely with colleagues from the Drug Regulatory Authorities to develop an addendum to the ICH S1B guideline on carcinogenicity studies that allows for a weight-of-evidence (WoE) carcinogenicity assessment in some cases, rather than conducting a 2-year rat carcinogenicity study. A subgroup of the EWG composed of regulators have published in this issue a detailed analysis of the Prospective Evaluation Study (PES) conducted under the auspices of the ICH S1B(R1) EWG. Based on the experience gained through the Prospective Evaluation Study (PES) process, industry members of the EWG have prepared the following commentary to aid sponsors in assessing the standard WoE factors, considering how novel investigative approaches may be used to support a WoE assessment, and preparing appropriate documentation of the WoE assessment for presentation to regulatory authorities. The commentary also reviews some of the implementation challenges sponsors must consider in developing a carcinogenicity assessment strategy. Finally, case examples drawn from previously marketed products are provided as a supplement to this commentary to provide additional examples of how WoE criteria may be applied. The information and opinions expressed in this commentary are aimed at increasing the quality of WoE assessments to ensure the successful implementation of this approach.

Liver-on-chip model and application in predictive genotoxicity and mutagenicity of drugs

Currently, there is no test system, whether in vitro or in vivo, capable of examining all endpoints required for genotoxicity evaluation used in pre-clinical drug safety assessment. The objective of this study was to develop a model which could assess all the required endpoints and possesses robust human metabolic activity, that could be used in a streamlined, animal-free manner. Liver-on-chip (LOC) models have intrinsic human metabolic activity that mimics the in vivo environment, making it a preferred test system. For our assay, the LOC was assembled using primary human hepatocytes or HepaRG cells, in a MPS-T12 plate, maintained under microfluidic flow conditions using the PhysioMimix® Microphysiological System (MPS), and co-cultured with human lymphoblastoid (TK6) cells in transwells. This system allows for interaction between two compartments and for the analysis of three different genotoxic endpoints, i.e. DNA strand breaks (comet assay) in hepatocytes, chromosome loss or damage (micronucleus assay) and mutation (Duplex Sequencing) in TK6 cells. Both compartments were treated at 0, 24 and 45 h with two direct genotoxicants: methyl methanesulfonate (MMS) and ethyl methanesulfonate (EMS), and two genotoxicants requiring metabolic activation: benzo[a]pyrene (B[a]P) and cyclophosphamide (CP). Assessment of cytochrome activity, RNA expression, albumin, urea and lactate dehydrogenase production, demonstrated functional metabolic capacities. Genotoxicity responses were observed for all endpoints with MMS and EMS. Increases in the micronucleus and mutations (MF) frequencies were also observed with CP, and %Tail DNA with B[a]P, indicating the metabolic competency of the test system. CP did not exhibit an increase in the %Tail DNA, which is in line with in vivo data. However, B[a]P did not exhibit an increase in the % micronucleus and MF, which might require an optimization of the test system. In conclusion, this proof-of-principle experiment suggests that LOC-MPS technology is a promising tool for in vitro hazard identification genotoxicants.

Principles of digital sequencing using unique molecular identifiers

Massively parallel sequencing technologies have long been used in both basic research and clinical routine. The recent introduction of digital sequencing has made previously challenging applications possible by significantly improving sensitivity and specificity to now allow detection of rare sequence variants, even at single molecule level. Digital sequencing utilizes unique molecular identifiers (UMIs) to minimize sequencing-induced errors and quantification biases. Here, we discuss the principles of UMIs and how they are used in digital sequencing. We outline the properties of different UMI types and the consequences of various UMI approaches in relation to experimental protocols and bioinformatics. Finally, we describe how digital sequencing can be applied in specific research fields, focusing on cancer management where it can be used in screening of asymptomatic individuals, diagnosis, treatment prediction, prognostication, monitoring treatment efficacy and early detection of treatment resistance as well as relapse.

Subclonal Cancer Driver Mutations Are Prevalent in the Unresected Peritumoral Edema of Adult Diffuse Gliomas

Adult diffuse gliomas commonly recur regardless of therapy. As recurrence typically arises from the peritumoral edema adjacent to the resected bulk tumor, the profiling of somatic mutations from infiltrative malignant cells within this critical, unresected region could provide important insights into residual disease. A key obstacle has been the inability to distinguish between next-generation sequencing (NGS) noise and the true but weak signal from tumor cells hidden among the noncancerous brain tissue of the peritumoral edema. Here, we developed and validated True2 sequencing to reduce NGS-associated errors to <1 false positive/100 kb panel positions while detecting 97.6% of somatic mutations with an allele frequency ≥0.1%. True2 was then used to study the tumor and peritumoral edema of 22 adult diffuse gliomas including glioblastoma, astrocytoma, oligodendroglioma, and NF1-related low-grade neuroglioma. The tumor and peritumoral edema displayed a similar mutation burden, indicating that surgery debulks these cancers physically but not molecularly. Moreover, variants in the peritumoral edema included unique cancer driver mutations absent in the bulk tumor. Finally, analysis of multiple samples from each patient revealed multiple subclones with unique mutations in the same gene in 17 of 22 patients, supporting the occurrence of convergent evolution in response to patient-specific selective pressures in the tumor microenvironment that may form the molecular foundation of recurrent disease. Collectively, True2 enables the detection of ultralow frequency mutations during molecular analyses of adult diffuse gliomas, which is necessary to understand cancer evolution, recurrence, and individual response to therapy.

SIGNIFICANCE

True2 is a next-generation sequencing workflow that facilitates unbiased discovery of somatic mutations across the full range of variant allele frequencies, which could help identify residual disease vulnerabilities for targeted adjuvant therapies.

Frequencies and spectra of aflatoxin B1-induced mutations in liver genomes of NEIL1-deficient mice as revealed by duplex sequencing

Increased risk for the development of hepatocellular carcinoma (HCC) is driven by a number of etiological factors including hepatitis viral infection and dietary exposures to foods contaminated with aflatoxin-producing molds. Intracellular metabolic activation of aflatoxin B1 (AFB1) to a reactive epoxide generates highly mutagenic AFB1-Fapy-dG adducts. Previously, we demonstrated that repair of AFB1-Fapy-dG adducts can be initiated by the DNA glycosylase NEIL1 and that male Neil1-/- mice were significantly more susceptible to AFB1-induced HCC relative to wild-type mice. To investigate the mechanisms underlying this enhanced carcinogenesis, WT and Neil1-/- mice were challenged with a single, 4 mg/kg dose of AFB1 and frequencies and spectra of mutations were analyzed in liver DNAs 2.5 months post-injection using duplex sequencing. The analyses of DNAs from AFB1-challenged mice revealed highly elevated mutation frequencies in the nuclear genomes of both males and females, but not the mitochondrial genomes. In both WT and Neil1-/- mice, mutation spectra were highly similar to the AFB1-specific COSMIC signature SBS24. Relative to wild-type, the NEIL1 deficiency increased AFB1-induced mutagenesis with concomitant elevated HCCs in male Neil1-/- mice. Our data establish a critical role of NEIL1 in limiting AFB1-induced mutagenesis and ultimately carcinogenesis.

5-Chloro-2'-deoxycytidine Induces a Distinctive High-Resolution Mutational Spectrum of Transition Mutations In Vivo

The biomarker 5-chlorocytosine (5ClC) appears in the DNA of inflamed tissues. Replication of a site-specific 5ClC in a viral DNA genome results in C → T mutations, which is consistent with 5ClC acting as a thymine mimic in vivo. Direct damage of nucleic acids by immune-cell-derived hypochlorous acid is one mechanism by which 5ClC could appear in the genome. A second, nonmutually exclusive mechanism involves damage of cytosine nucleosides or nucleotides in the DNA precursor pool, with subsequent utilization of the 5ClC deoxynucleotide triphosphate as a precursor for DNA synthesis. The present work characterized the mutagenic properties of 5ClC in the nucleotide pool by exposing cells to the nucleoside 5-chloro-2'-deoxycytidine (5CldC). In both Escherichia coli and mouse embryonic fibroblasts (MEFs), 5CldC in the growth media was potently mutagenic, indicating that 5CldC enters cells and likely is erroneously incorporated into the genome from the nucleotide pool. High-resolution sequencing of DNA from MEFs derived from the gptΔ C57BL/6J mouse allowed qualitative and quantitative characterization of 5CldC-induced mutations; CG → TA transitions in 5'-GC(Y)-3' contexts (Y = a pyrimidine) were dominant, while TA → CG transitions appeared at a much lower frequency. The high-resolution mutational spectrum of 5CldC revealed a notable similarity to the Catalogue of Somatic Mutations in Cancer mutational signatures SBS84 and SBS42, which appear in human lymphoid tumors and in occupationally induced cholangiocarcinomas, respectively. SBS84 is associated with the expression of activation-induced cytidine deaminase (AID), a cytosine deaminase associated with inflammation, as well as immunoglobulin gene diversification during antibody maturation. The similarity between the spectra of AID activation and 5CldC could be coincidental; however, the administration of 5CldC did induce some AID expression in MEFs, which have no inherent expression of its gene. In summary, this work shows that 5CldC induces a distinct pattern of mutations in cells. Moreover, that pattern resembles human mutational signatures induced by inflammatory processes, such as those triggered in certain malignancies.

Application of duplex sequencing to evaluate mutagenicity of aristolochic acid and methapyrilene in Fisher 344 rats

Duplex sequencing (DS) is an error-corrected next-generation sequencing (NGS) method that can overcome notorious high error rate from the process of NGS and detect ultralow-frequency mutations. In this study, we evaluated the mutagenicity of aristolochic acid, a known genotoxic carcinogen, and methapyrilene, a known nongenotoxic carcinogen using DS. Four male Fisher 344 rats were treated with aristolochic acid, methapyrilene, or the vehicle control for 6 weeks, liver tissues were collected one day after the treatment, and the DNA was isolated for analysis. The mutation frequency for the aristolochic acid-treated group was significantly increased over the vehicle control (44-fold), whereas no significant difference in the mutation frequency was observed between the methapyrilene-treated and the control groups. The primary type of mutation induced by aristolochic acid was A:T > T:A transversion, which occurred frequently at ApT sites, whereas the major type of mutation in the control and methapyrilene-treated groups was G:C > A:T transition, which occurred frequently at CpG sites. These findings are consistent with previously published data obtained with other in vivo mutation assays. Thus, our results suggest that the DS mutation assay is a promising technology for assessing mutagenicity of chemicals in vivo.

MtDNA deletions and aging

Aging is the major risk factor in most of the leading causes of mortality worldwide, yet its fundamental causes mostly remain unclear. One of the clear hallmarks of aging is mitochondrial dysfunction. Mitochondria are best known for their roles in cellular energy generation, but they are also critical biosynthetic and signaling organelles. They also undergo multiple changes with organismal age, including increased genetic errors in their independent, circular genome. A key group of studies looking at mice with increased mtDNA mutations showed that premature aging phenotypes correlated with increased deletions but not point mutations. This generated an interest in mitochondrial deletions as a potential fundamental cause of aging. However, subsequent studies in different models have yielded diverse results. This review summarizes the research on mitochondrial deletions in various organisms to understand their possible roles in causing aging while identifying the key complications in quantifying deletions across all models.

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 and discrepancies have been observed between different techniques for MRD assessment. In 62 patients with AML, aged 18-60 years, in first complete remission after intensive induction therapy on the randomized phase III SWOG-S0106 clinical trial (clinicaltrials gov. Identifier: NCT00085709), MRD detection by centralized, high-quality multiparametric flow cytometry was compared with a 29-gene panel utilizing duplex sequencing (DS), an ultrasensitive next-generation sequencing method that generates double-stranded consensus sequences to reduce false positive errors. MRD as defined by DS was observed in 22 (35%) patients and was strongly associated with higher rates of relapse (68% vs. 13%; hazard ratio [HR] =8.8; 95% confidence interval [CI]: 3.2-24.5; P<0.001) and decreased survival (32% vs. 82%; HR=5.6; 95% CI: 2.3-13.8; P<0.001) at 5 years. DS MRD strongly outperformed multiparametric flow cytometry MRD, which was observed in ten (16%) patients and marginally associated with higher rates of relapse (50% vs. 30%; HR=2.4; 95% CI: 0.9-6.7; P=0.087) and decreased survival (40% vs. 68%; HR=2.5; 95% CI: 1.0-6.3; P=0.059) at 5 years. Furthermore, the prognostic significance of DS MRD status at the time of remission for subsequent relapse was similar on both randomized arms of the trial. These findings suggest that next-generation sequencing-based AML MRD testing is a powerful tool that could be developed for use in patient management and for early anti-leukemic treatment assessment in clinical trials.

Genetic profiling of rat gliomas and cardiac schwannomas from life-time radiofrequency radiation exposure study using a targeted next-generation sequencing gene panel

The cancer hazard associated with lifetime exposure to radiofrequency radiation (RFR) was examined in Sprague Dawley (SD) rats at the Ramazzini Institute (RI), Italy. There were increased incidences of gliomas and cardiac schwannomas. The translational relevance of these rare rat tumors for human disease is poorly understood. We examined the genetic alterations in RFR-derived rat tumors through molecular characterization of important cancer genes relevant for human gliomagenesis. A targeted next-generation sequencing (NGS) panel was designed for rats based on the top 23 orthologous human glioma-related genes. Single-nucleotide variants (SNVs) and small insertion and deletions (indels) were characterized in the rat gliomas and cardiac schwannomas. Translational relevance of these genetic alterations in rat tumors to human disease was determined through comparison with the Catalogue of Somatic Mutations in Cancer (COSMIC) database. These data suggest that rat gliomas resulting from life-time exposure to RFR histologically resemble low grade human gliomas but surprisingly no mutations were detected in rat gliomas that had homology to the human IDH1 p.R132 or IDH2 p.R172 suggesting that rat gliomas are primarily wild-type for IDH hotspot mutations implicated in human gliomas. The rat gliomas appear to share some genetic alterations with IDH1 wildtype human gliomas and rat cardiac schwannomas also harbor mutations in some of the queried cancer genes. These data demonstrate that targeted NGS panels based on tumor specific orthologous human cancer driver genes are an important tool to examine the translational relevance of rodent tumors resulting from chronic/life-time rodent bioassays.

CRISPR/Cas9 as a Mutagenic Factor

The discovery of the CRISPR/Cas9 microbial adaptive immune system has revolutionized the field of genetics, by greatly enhancing the capacity for genome editing. CRISPR/Cas9-based editing starts with DNA breaks (or other lesions) predominantly at target sites and, unfortunately, at off-target genome sites. DNA repair systems differing in accuracy participate in establishing desired genetic changes but also introduce unwanted mutations, that may lead to hereditary, oncological, and other diseases. New approaches to alleviate the risks associated with genome editing include attenuating the off-target activity of editing complex through the use of modified forms of Cas9 nuclease and single guide RNA (sgRNA), improving delivery methods for sgRNA/Cas9 complex, and directing DNA lesions caused by the sgRNA/Cas9 to non-mutagenic repair pathways. Here, we have described CRISPR/Cas9 as a new powerful mutagenic factor, discussed its mutagenic properties, and reviewed factors influencing the mutagenic activity of CRISPR/Cas9.

Interpretation of in vitro concentration-response data for risk assessment and regulatory decision-making: Report from the 2022 IWGT quantitative analysis expert working group meeting

Quantitative risk assessments of chemicals are routinely performed using in vivo data from rodents; however, there is growing recognition that non-animal approaches can be human-relevant alternatives. There is an urgent need to build confidence in non-animal alternatives given the international support to reduce the use of animals in toxicity testing where possible. In order for scientists and risk assessors to prepare for this paradigm shift in toxicity assessment, standardization and consensus on in vitro testing strategies and data interpretation will need to be established. To address this issue, an Expert Working Group (EWG) of the 8th International Workshop on Genotoxicity Testing (IWGT) evaluated the utility of quantitative in vitro genotoxicity concentration-response data for risk assessment. The EWG first evaluated available in vitro methodologies and then examined the variability and maximal response of in vitro tests to estimate biologically relevant values for the critical effect sizes considered adverse or unacceptable. Next, the EWG reviewed the approaches and computational models employed to provide human-relevant dose context to in vitro data. Lastly, the EWG evaluated risk assessment applications for which in vitro data are ready for use and applications where further work is required. The EWG concluded that in vitro genotoxicity concentration-response data can be interpreted in a risk assessment context. However, prior to routine use in regulatory settings, further research will be required to address the remaining uncertainties and limitations.

Clonal selection of hematopoietic stem cells after gene therapy for sickle cell disease

Gene therapy (GT) provides a potentially curative treatment option for patients with sickle cell disease (SCD); however, the occurrence of myeloid malignancies in GT clinical trials has prompted concern, with several postulated mechanisms. Here, we used whole-genome sequencing to track hematopoietic stem cells (HSCs) from six patients with SCD at pre- and post-GT time points to map the somatic mutation and clonal landscape of gene-modified and unmodified HSCs. Pre-GT, phylogenetic trees were highly polyclonal and mutation burdens per cell were elevated in some, but not all, patients. Post-GT, no clonal expansions were identified among gene-modified or unmodified cells; however, an increased frequency of potential driver mutations associated with myeloid neoplasms or clonal hematopoiesis (DNMT3A- and EZH2-mutated clones in particular) was observed in both genetically modified and unmodified cells, suggesting positive selection of mutant clones during GT. This work sheds light on HSC clonal dynamics and the mutational landscape after GT in SCD, highlighting the enhanced fitness of some HSCs harboring pre-existing driver mutations. Future studies should define the long-term fate of mutant clones, including any contribution to expansions associated with myeloid neoplasms.

Genome-wide direct quantification of in vivo mutagenesis using high-accuracy paired-end and complementary consensus sequencing

Error-corrected next-generation sequencing (ecNGS) is an emerging technology for accurately measuring somatic mutations. Here, we report paired-end and complementary consensus sequencing (PECC-Seq), a high-accuracy ecNGS approach for genome-wide somatic mutation detection. We characterize a novel 2-aminoimidazolone lesion besides 7,8-dihydro-8-oxoguanine and the resulting end-repair artifacts originating from NGS library preparation that obscure the sequencing accuracy of NGS. We modify library preparation protocol for the enzymatic removal of end-repair artifacts and improve the accuracy of our previously developed duplex consensus sequencing method. Optimized PECC-Seq shows an error rate of <5 × 10-8 with consensus bases compressed from approximately 25 Gb of raw sequencing data, enabling the accurate detection of low-abundance somatic mutations. We apply PECC-Seq to the quantification of in vivo mutagenesis. Compared with the classic gpt gene mutation assay using gpt delta transgenic mice, PECC-Seq exhibits high sensitivity in quantitatively measuring dose-dependent mutagenesis induced by Aristolochic acid I (AAI). Moreover, PECC-Seq specifically characterizes the distinct genome-wide mutational signatures of AAI, Benzo[a]pyrene, N-Nitroso-N-ethylurea and N-nitrosodiethylamine and reveals the mutational signature of Quinoline in common mouse models. Overall, our findings demonstrate that high-accuracy PECC-Seq is a promising tool for genome-wide somatic mutagenesis quantification and for in vivo mutagenicity testing.

Adopting duplex sequencing technology for genetic toxicity testing: A proof-of-concept mutagenesis experiment with N-ethyl-N-nitrosourea (ENU)-exposed rats

Duplex sequencing (DS) is an error-corrected next-generation sequencing method in which molecular barcodes informatically link PCR-copies back to their source DNA strands, enabling computational removal of errors in consensus sequences. The resulting background of less than one artifactual mutation per 107 nucleotides allows for direct detection of somatic mutations. TwinStrand Biosciences, Inc. has developed a DS-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.

Whole-genome high-fidelity sequencing: A novel approach to detecting and characterization of mutagenicity in vivo

Direct DNA sequencing can be used for characterizing mutagenicity in simple and complex biological models. Recently we described a method of whole-genome sequencing for detecting mutations in simple models of cultured bacteria, mammalian cells, and nematode. In the current proof-of-concept study, we expand and improve our method for evaluating a more complex mammalian biological model in outbred mice. We detail the method by applying it to a small set of animals treated with a mutagen with known mutagenicity profiles, N-ethyl-N-nitrosourea (ENU), for consistency with the known data. Whole-genome high-fidelity sequencing (HiFi Sequencing) showed frequencies and spectra of background mutations in tissues of untreated mice that were consistent with normal ageing and characterized by spontaneous or enzymatic deamination of 5-methylcytosine. In mice treated with a single 40 mg/kg dose of ENU, the frequency of mutations in the genomic DNA of solid tissues increased up to 7-fold, with the greatest increase observed in the spleen and the smallest increase in the liver. The most common mutations detected in ENU-treated mice were T > A transitions and T > C transversions, consistent with the types of mutations caused by alkylating agents. The data suggest that HiFi Sequencing may be useful for characterizing mutagenicity of novel compounds in various biological models.

Comparison of the transgenic rodent mutation assay, error corrected next generation duplex sequencing, and the alkaline comet assay to detect dose-related mutations following exposure to N-nitrosodiethylamine

N-Nitrosodiethylamine (NDEA), a well-studied N-nitrosamine, was tested in rats to compare the dose-response relationship of three genotoxicity endpoints. Mutant / mutation frequencies were determined using the transgenic rodent (TGR) gene mutation assay and error corrected next generation sequencing (ecNGS) (i.e., duplex sequencing (DS)), and genetic damage was detected by the alkaline comet assay. Big Blue® (cII Locus) animals (n = 6 per dose group) were administered doses of 0.001, 0.01, 0.1, 1, 3 mg/kg/day NDEA by oral gavage. Samples were collected for cII mutation and DS analyses following 28-days of exposure and 3 days recovery. In a separate study, male Sprague-Dawley (SD) rats (n = 6 per dose group) were administered the same doses by oral gavage for two consecutive days and then samples collected for the alkaline comet assay. A dose-related increase in mutant / mutation frequencies of the liver but not duodenum was observed using the TGR assay and DS with DS resulting in a slightly more sensitive response, with a lower benchmark dose (BMD). In addition, a dose-related increase in percent tail DNA was observed in the liver using the alkaline comet assay. Therefore, DS and comet assays showed good utility for hazard identification and dose-response analysis of a representative N-nitrosamine comparable to the TGR gene mutation assay.

Genotoxicity assessment: opportunities, challenges and perspectives for quantitative evaluations of dose-response data

Genotoxicity data are mainly interpreted in a qualitative way, which typically results in a binary classification of chemical entities. For more than a decade, there has been a discussion about the need for a paradigm shift in this regard. Here, we review current opportunities, challenges and perspectives for a more quantitative approach to genotoxicity assessment. Currently discussed opportunities mainly include the determination of a reference point (e.g., a benchmark dose) from genetic toxicity dose-response data, followed by calculation of a margin of exposure (MOE) or derivation of a health-based guidance value (HBGV). In addition to new opportunities, major challenges emerge with the quantitative interpretation of genotoxicity data. These are mainly rooted in the limited capability of standard in vivo genotoxicity testing methods to detect different types of genetic damage in multiple target tissues and the unknown quantitative relationships between measurable genotoxic effects and the probability of experiencing an adverse health outcome. In addition, with respect to DNA-reactive mutagens, the question arises whether the widely accepted assumption of a non-threshold dose-response relationship is at all compatible with the derivation of a HBGV. Therefore, at present, any quantitative genotoxicity assessment approach remains to be evaluated case-by-case. The quantitative interpretation of in vivo genotoxicity data for prioritization purposes, e.g., in connection with the MOE approach, could be seen as a promising opportunity for routine application. However, additional research is needed to assess whether it is possible to define a genotoxicity-derived MOE that can be considered indicative of a low level of concern. To further advance quantitative genotoxicity assessment, priority should be given to the development of new experimental methods to provide a deeper mechanistic understanding and a more comprehensive basis for the analysis of dose-response relationships.

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.

Structural underpinnings of mutation rate variations in the human genome

Single nucleotide mutation rates have critical implications for human evolution and genetic diseases. Importantly, the rates vary substantially across the genome and the principles underlying such variations remain poorly understood. A recent model explained much of this variation by considering higher-order nucleotide interactions in the 7-mer sequence context around mutated nucleotides. This model's success implicates a connection between DNA shape and mutation rates. DNA shape, i.e. structural properties like helical twist and tilt, is known to capture interactions between nucleotides within a local context. Thus, we hypothesized that changes in DNA shape features at and around mutated positions can explain mutation rate variations in the human genome. Indeed, DNA shape-based models of mutation rates showed similar or improved performance over current nucleotide sequence-based models. These models accurately characterized mutation hotspots in the human genome and revealed the shape features whose interactions underlie mutation rate variations. DNA shape also impacts mutation rates within putative functional regions like transcription factor binding sites where we find a strong association between DNA shape and position-specific mutation rates. This work demonstrates the structural underpinnings of nucleotide mutations in the human genome and lays the groundwork for future models of genetic variations to incorporate DNA shape.

Error-corrected next generation sequencing - Promises and challenges for genotoxicity and cancer risk assessment

Error-corrected Next Generation Sequencing (ecNGS) is rapidly emerging as a valuable, highly sensitive and accurate method for detecting and characterizing mutations in any cell type, tissue or organism from which DNA can be isolated. Recent mutagenicity and carcinogenicity studies have used ecNGS to quantify drug-/chemical-induced mutations and mutational spectra associated with cancer risk. ecNGS has potential applications in genotoxicity assessment as a new readout for traditional models, for mutagenesis studies in 3D organotypic cultures, and for detecting off-target effects of gene editing tools. Additionally, early data suggest that ecNGS can measure clonal expansion of mutations as a mechanism-agnostic early marker of carcinogenic potential and can evaluate mutational load directly in human biomonitoring studies. In this review, we discuss promising applications, challenges, limitations, and key data initiatives needed to enable regulatory testing and adoption of ecNGS - including for advancing safety assessment, augmenting weight-of-evidence for mutagenicity and carcinogenicity mechanisms, identifying early biomarkers of cancer risk, and managing human health risk from chemical exposures.

Next-generation sequencing methodologies to detect low-frequency mutations: "Catch me if you can"

Mutations, the irreversible changes in an organism's DNA sequence, are present in tissues at a variant allele frequency (VAF) ranging from ∼10-8 per bp for a founder mutation to ∼10-3 for a histologically normal tissue sample containing several independent clones - compared to 1%- 50% for a heterozygous tumor mutation or a polymorphism. The rarity of these events poses a challenge for accurate clinical diagnosis and prognosis, toxicology, and discovering new disease etiologies. Standard Next-Generation Sequencing (NGS) technologies report VAFs as low as 0.5% per nt, but reliably observing rarer precursor events requires additional sophistication to measure ultralow-frequency mutations. We detail the challenge; define terms used to characterize the results, which vary between laboratories and sometimes conflict between biologists and bioinformaticists; and describe recent innovations to improve standard NGS methodologies including: single-strand consensus sequence methods such as Safe-SeqS and SiMSen-Seq; tandem-strand consensus sequence methods such as o2n-Seq and SMM-Seq; and ultrasensitive parent-strand consensus sequence methods such as DuplexSeq, PacBio HiFi, SinoDuplex, OPUSeq, EcoSeq, BotSeqS, Hawk-Seq, NanoSeq, SaferSeq, and CODEC. Practical applications are also noted. Several methods quantify VAF down to 10-5 at a nt and mutation frequency (MF) in a target region down to 10-7 per nt. By expanding to > 1 Mb of sites never observed twice, thus forgoing VAF, other methods quantify MF < 10-9 per nt or < 15 errors per haploid genome. Clonal expansion cannot be directly distinguished from independent mutations by sequencing, so it is essential for a paper to report whether its MF counted only different mutations - the minimum independent-mutation frequency MFminI - or all mutations observed including recurrences - the larger maximum independent-mutation frequency MFmaxI which may reflect clonal expansion. Ultrasensitive methods reveal that, without their use, even mutations with VAF 0.5-1% are usually spurious.

International Workshops on Genotoxicity Testing (IWGT): Origins, achievements and ambitions

Over the past thirty years, the International Workshops on Genotoxicity Testing (IWGT) became one of the leading groups in the field of regulatory genotoxicology, not only due to the diversity of participants with respect to geography and professional affiliation, but also due to the unique setup of recurring IWGT meetings every four years. The hallmarks of the IWGT process have been diligent initial planning approaches of the working groups, collection of data so as to stimulate data-driven discussions and debate, and striving to reach consensus recommendations. The scientific quality of the Working Groups (WGs) has been exceptional due to the selection of highly regarded experts on each topic. As a result, the IWGT working group reports have become important documents. The deliberations and publications have provided guidance on test systems and testing protocols that have influenced the development or revision of test guidelines of the Organisation for Economic Co-operation and Development (OECD), guidance by the International Council for Harmonisation (ICH), and strategic testing or data analysis approaches in general. This article summarizes the history of the IWGT, identifies some of its major achievements, and provides an outlook for the future.

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.

svCapture: efficient and specific detection of very low frequency structural variant junctions by error-minimized capture sequencing

Error-corrected sequencing of genomic targets enriched by probe-based capture has become a standard approach for detecting single-nucleotide variants (SNVs) and small insertion/deletions (indels) present at very low variant allele frequencies. Less attention has been given to comparable strategies for rare structural variant (SV) junctions, where different error mechanisms must be addressed. Working from samples with known SV properties, we demonstrate that duplex sequencing (DuplexSeq), which demands confirmation of variants on both strands of a source DNA molecule, eliminates false SV junctions arising from chimeric PCR. DuplexSeq could not address frequent intermolecular ligation artifacts that arise during Y-adapter addition prior to strand denaturation without requiring multiple source molecules. In contrast, tagmentation libraries coupled with data filtering based on strand family size greatly reduced both artifact classes and enabled efficient and specific detection of single-molecule SV junctions. The throughput of SV capture sequencing (svCapture) and base-level accuracy of DuplexSeq provided detailed views of the microhomology profile and limited occurrence of de novo SNVs near the junctions of hundreds of newly created SVs, suggesting end joining as a possible formation mechanism. The open source svCapture pipeline enables rare SV detection as a routine addition to SNVs/indels in properly prepared capture sequencing libraries.

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.

Lung adenocarcinoma promotion by air pollutants

A complete understanding of how exposure to environmental substances promotes cancer formation is lacking. More than 70 years ago, tumorigenesis was proposed to occur in a two-step process: an initiating step that induces mutations in healthy cells, followed by a promoter step that triggers cancer development1. Here we propose that environmental particulate matter measuring ≤2.5 μm (PM2.5), known to be associated with lung cancer risk, promotes lung cancer by acting on cells that harbour pre-existing oncogenic mutations in healthy lung tissue. Focusing on EGFR-driven lung cancer, which is more common in never-smokers or light smokers, we found a significant association between PM2.5 levels and the incidence of lung cancer for 32,957 EGFR-driven lung cancer cases in four within-country cohorts. Functional mouse models revealed that air pollutants cause an influx of macrophages into the lung and release of interleukin-1β. This process results in a progenitor-like cell state within EGFR mutant lung alveolar type II epithelial cells that fuels tumorigenesis. Ultradeep mutational profiling of histologically normal lung tissue from 295 individuals across 3 clinical cohorts revealed oncogenic EGFR and KRAS driver mutations in 18% and 53% of healthy tissue samples, respectively. These findings collectively support a tumour-promoting role for  PM2.5 air pollutants  and provide impetus for public health policy initiatives to address air pollution to reduce disease burden.

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.

Strategies for Assessing Acceptable Intakes for Novel N-Nitrosamines Derived from Active Pharmaceutical Ingredients

The detection of N-nitrosamines, derived from solvents and reagents and, on occasion, the active pharmaceutical ingredient (API) at higher than acceptable levels in drug products, has led regulators to request a detailed review for their presence in all medicinal products. In the absence of rodent carcinogenicity data for novel N-nitrosamines derived from amine-containing APIs, a conservative class limit of 18 ng/day (based on the most carcinogenic N-nitrosamines) or the derivation of acceptable intakes (AIs) using structurally related surrogates with robust rodent carcinogenicity data is recommended. The guidance has implications for the pharmaceutical industry given the vast number of marketed amine-containing drugs. In this perspective, the rate-limiting step in N-nitrosamine carcinogenicity, involving cytochrome P450-mediated α-carbon hydroxylation to yield DNA-reactive diazonium or carbonium ion intermediates, is discussed with reference to the selection of read-across analogs to derive AIs. Risk-mitigation strategies for managing putative N-nitrosamines in the preclinical discovery setting are also presented.

Liposome-Mediated Gene Transfer in Differentiated HepaRG™ Cells: Expression of Liver Specific Functions and Application to the Cytochrome P450 2D6 Expression

The goal of this study was to establish a procedure for gene delivery mediated by cationic liposomes in quiescent differentiated HepaRG™ human hepatoma cells. We first identified several cationic lipids promoting efficient gene transfer with low toxicity in actively dividing HepG2, HuH7, BC2 and progenitor HepaRG™ human hepatoma cells. The lipophosphoramidate Syn1-based nanovector, which allowed the highest transfection efficiencies of progenitor HepaRG™ cells, was next used to transfect differentiated HepaRG™ cells. Lipofection of these cells using Syn1-based liposome was poorly efficient most likely because the differentiated HepaRG™ cells are highly quiescent. Thus, we engineered the differentiated HepaRG™ Mitogenic medium supplement (ADD1001) that triggered robust proliferation of differentiated cells. Importantly, we characterized the phenotypical changes occurring during proliferation of differentiated HepaRG™ cells and demonstrated that mitogenic stimulation induced a partial and transient decrease in the expression levels of some liver specific functions followed by a fast recovery of the full differentiation status upon removal of the mitogens. Taking advantage of the proliferation of HepaRG™ cells, we defined lipofection conditions using Syn1-based liposomes allowing transient expression of the cytochrome P450 2D6, a phase I enzyme poorly expressed in HepaRG cells, which opens new means for drug metabolism studies in HepaRG™ cells.

Resistance of mitochondrial DNA to cadmium and Aflatoxin B1 damage-induced germline mutation accumulation in C. elegans

Mitochondrial DNA (mtDNA) is prone to mutation in aging and over evolutionary time, yet the processes that regulate the accumulation of de novo mtDNA mutations and modulate mtDNA heteroplasmy are not fully elucidated. Mitochondria lack certain DNA repair processes, which could contribute to polymerase error-induced mutations and increase susceptibility to chemical-induced mtDNA mutagenesis. We conducted error-corrected, ultra-sensitive Duplex Sequencing to investigate the effects of two known nuclear genome mutagens, cadmium and Aflatoxin B1, on germline mtDNA mutagenesis in Caenorhabditis elegans. Detection of thousands of mtDNA mutations revealed pervasive heteroplasmy in C. elegans and that mtDNA mutagenesis is dominated by C:G → A:T mutations generally attributed to oxidative damage. However, there was no effect of either exposure on mtDNA mutation frequency, spectrum, or trinucleotide context signature despite a significant increase in nuclear mutation rate after aflatoxin B1 exposure. Mitophagy-deficient mutants pink-1 and dct-1 accumulated significantly higher levels of mtDNA damage compared to wild-type C. elegans after exposures. However, there were only small differences in mtDNA mutation frequency, spectrum, or trinucleotide context signature compared to wild-type after 3050 generations, across all treatments. These findings suggest mitochondria harbor additional previously uncharacterized mechanisms that regulate mtDNA mutational processes across generations.

Duplex sequencing identifies genomic features that determine susceptibility to benzo(a)pyrene-induced in vivo mutations

Exposure to environmental mutagens increases the risk of cancer and genetic disorders. We used Duplex Sequencing (DS), a high-accuracy error-corrected sequencing technology, to analyze mutation induction across twenty 2.4 kb intergenic and genic targets in the bone marrow of MutaMouse males exposed to benzo(a)pyrene (BaP), a widespread environmental pollutant. DS revealed a linear dose-related induction of mutations across all targets with low intra-group variability. Heterochromatic and intergenic regions exhibited the highest mutation frequencies (MF). C:G > A:T transversions at CCA, CCC and GCC trinucleotides were enriched in BaP-exposed mice consistent with the known etiology of BaP mutagenesis. However, GC-content had no effect on mutation susceptibility. A positive correlation was observed between DS and the “gold-standard” transgenic rodent gene mutation assay. Overall, we demonstrate that DS is a promising approach to study in vivo mutagenesis and yields critical insight into the genomic features governing mutation susceptibility, spectrum, and variability across the genome.

OPUSeq simplifies detection of low-frequency DNA variants and uncovers fragmentase-associated artifacts

Detection of low-frequency DNA variants (below 1%) is becoming increasingly important in biomedical research and clinical practice, but is challenging to do with standard sequencing approaches due to high error rates. The use of double-stranded unique molecular identifiers (dsUMIs) allows correction of errors by comparing reads arising from the same original DNA duplex. However, the implementation of such approaches is still challenging. Here, we present a novel method, one-pot dsUMI sequencing (OPUSeq), which allows incorporation of dsUMIs in the same reaction as the library PCR. This obviates the need for adapter pre-synthesis or additional enzymatic steps. OPUSeq can be incorporated into standard DNA library preparation approaches and coupled with hybridization target capture. We demonstrate successful error correction and detection of variants down to allele frequency of 0.01%. Using OPUSeq, we also show that the use of enzymatic fragmentation can lead to the appearance of spurious double-stranded variants, interfering with detection of variant fractions below 0.1%.

Non-canonical genomic driver mutations of urethane carcinogenesis

The carcinogen urethane induces pulmonary tumors in mice initiated by an incredibly specific Q61L/R oncogenic mutation in the proto-oncogene Kras. Previous Whole-Exome Sequencing of urethane-induced tumors revealed a bias towards A➙T/G and G➙A substitutions. Subsequent ultra-sensitive Maximum-Depth Sequencing of Kras shortly after urethane exposure suggest a further refinement to CA➙CT/G substitutions. As C182AA➙C182T/GA substitutions in Kras result in Q61L/R mutations, the extreme bias of urethane towards these genomic driver mutations can be ascribed to the specificity of the carcinogen for CA➙CT/G substitutions. However, we previously found that changing rare codons to common in the Kras gene to increase protein expression shifted mutations in urethane-induced tumors away from Kras, or when detected in Kras, to G12D mutations that are usually rarely detected in such tumors. Moreover, the loss of p53 partially reversed this effect, generating tumors with either Q61L/R or G12D oncogenic Kras mutations, or no Kras mutations, presumably due to other genomic driver mutations. Determining the origin of these G12D and other unknown non-canonical genomic driver mutations would provide critical insight into the extreme bias of carcinogens for specific genomic driver mutations. We thus compared the types of Single Nucleotide Variations detected by previously performed Maximum-Depth Sequencing immediately after urethane exposure to the mutation signatures derived from Whole Exome Sequencing of urethane-induced tumors. This identified two types of non-canonical mutations. First, a V637E oncogenic mutation in the proto-oncogene Braf that conforms to the mutation signature of urethane, suggesting that the mutational bias of the carcinogen may account for this non-canonical mutation, similar to that for canonical Q61L/R mutations in Kras. Second, G12D and Q61H mutations in Kras that did not fit this mutation signature, and instead shared similarity with Single Nucleotide Variations detected by Maximum-Depth Sequencing from normal cells, suggesting that perhaps these mutations were pre-existing. We thus posit that when canonical Kras mutations are selected against that the carcinogen may instead promote the expansion of pre-existing genomic driver mutations, although admittedly we cannot rule out other mechanisms. Interrogating the mutation signatures of human lung cancers similarly identified KRAS genomic driver mutations that failed to match the mutation signature of the tumor. Thus, we also speculate that the selection for non-canonical genomic driver mutations during urethane carcinogenesis may reflect the process by which discordance between genomic driver mutations and mutational signatures arises in human cancers.

Preparation of Duplex Sequencing Libraries for Archival Paraffin-Embedded Tissue Samples Using Single-Strand-Specific Nuclease P1

DNA from formalin-fixed paraffin-embedded (FFPE) tissues, which are frequently utilized in cancer research, is significantly affected by chemical degradation. It was suggested that approaches that are based on duplex sequencing can significantly improve the accuracy of mutation detection in FFPE-derived DNA. However, the original duplex sequencing method cannot be utilized for the analysis of formalin-fixed paraffin-embedded (FFPE) tissues, as FFPE DNA contains an excessive number of damaged bases, and these lesions are converted to false double-strand nucleotide substitutions during polymerase-driven DNA end repair process. To resolve this drawback, we replaced DNA polymerase by a single strand-specific nuclease P1. Nuclease P1 was shown to efficiently remove RNA from DNA preparations, to fragment the FFPE-derived DNA and to remove 5'/3'-overhangs. To assess the performance of duplex sequencing-based methods in FFPE-derived DNA, we constructed the Bottleneck Sequencing System (BotSeqS) libraries from five colorectal carcinomas (CRCs) using either DNA polymerase or nuclease P1. As expected, the number of identified mutations was approximately an order of magnitude higher in libraries prepared with DNA polymerase vs. nuclease P1 (626 ± 167/Mb vs. 75 ± 37/Mb, paired t-test p-value 0.003). Furthermore, the use of nuclease P1 but not polymerase-driven DNA end repair allowed a reliable discrimination between CRC tumors with and without hypermutator phenotypes. The utility of newly developed modification was validated in the collection of 17 CRCs and 5 adjacent normal tissues. Nuclease P1 can be recommended for the use in duplex sequencing library preparation from FFPE-derived DNA.

Colorectal Cancer Is Associated with the Presence of Cancer Driver Mutations in Normal Colon

Although somatic mutations in colorectal cancer are well characterized, little is known about the accumulation of cancer mutations in the normal colon before cancer. Here, we have developed and applied an ultrasensitive, single-molecule mutational test based on CRISPR-DS technology, which enables mutation detection at extremely low frequency (<0.001) in normal colon from patients with and without colorectal cancer. This testing platform revealed that normal colon from patients with and without colorectal cancer carries mutations in common colorectal cancer genes, but these mutations are more abundant in patients with cancer. Oncogenic KRAS mutations were observed in the normal colon of about one third of patients with colorectal cancer but in none of the patients without colorectal cancer. Patients with colorectal cancer also carried more TP53 mutations than patients without cancer and these mutations were more pathogenic and formed larger clones, especially in patients with early-onset colorectal cancer. Most mutations in the normal colon were different from the driver mutations in tumors, suggesting that the occurrence of independent clones with pathogenic KRAS and TP53 mutations is a common event in the colon of individuals who develop colorectal cancer. These results indicate that somatic evolution contributes to clonal expansions in the normal colon and that this process is enhanced in individuals with cancer, particularly in those with early-onset colorectal cancer.

SIGNIFICANCE

This work suggests prevalent somatic evolution in the normal colon of patients with colorectal cancer, highlighting the potential of using ultrasensitive gene sequencing to predict disease risk.

A Collaborative Initiative to Establish Genomic Biomarkers for Assessing Tumorigenic Potential to Reduce Reliance on Conventional Rodent Carcinogenicity Studies

There is growing recognition across broad sectors of the scientific community that use of genomic biomarkers has the potential to reduce the need for conventional rodent carcinogenicity studies of industrial chemicals, agrochemicals, and pharmaceuticals through a weight-of-evidence approach. These biomarkers fall into 2 major categories: (1) sets of gene transcripts that can identify distinct tumorigenic mechanisms of action; and (2) cancer driver gene mutations indicative of rapidly expanding growth-advantaged clonal cell populations. This call-to-action article describes a collaborative approach launched to develop and qualify biomarker gene expression panels that measure widely accepted molecular pathways linked to tumorigenesis and their activation levels to predict tumorigenic doses of chemicals from short-term exposures. Growing evidence suggests that application of such biomarker panels in short-term exposure rodent studies can identify both tumorigenic hazard and tumorigenic activation levels for chemical-induced carcinogenicity. In the future, this approach will be expanded to include methodologies examining mutations in key cancer driver gene mutation hotspots as biomarkers of both genotoxic and nongenotoxic chemical tumor risk. Analytical, technical, and biological validation studies of these complementary genomic tools are being undertaken by multisector and multidisciplinary collaborative teams within the Health and Environmental Sciences Institute. Success from these efforts will facilitate the transition from current heavy reliance on conventional 2-year rodent carcinogenicity studies to more rapid animal- and resource-sparing approaches for mechanism-based carcinogenicity evaluation supporting internal and regulatory decision-making.

Accurate Detection of Subclonal Variants in Paired Diagnosis-Relapse Acute Myeloid Leukemia Samples by Next Generation Duplex Sequencing

Mutations characterize diverse human cancers; there is a positive correlation between elevated mutation frequency and tumor progression. One exception is acute myeloid leukemia (AML), which has few clonal single nucleotide mutations. We used highly sensitive and accurate Duplex Sequencing (DS) to show now that AML, in addition, has an extensive repertoire of variants with low allele frequencies, < 1%, which is below the accurate detection limit of most other sequencing methodologies. The subclonal variants are unique to each individual and change in composition, frequency, and sequence context from diagnosis to relapse. Their functional significance is apparent by the observation that many are known variants and cluster within functionally important protein domains. Subclones provide a reservoir of variants that could expand and contribute to the development of drug resistance and relapse. In accord, we accurately identified subclonal variants in AML driver genes NRAS and RUNX1 at allele frequencies between 0.1% and 0.3% at diagnosis, which expanded to comprise a major fraction (14-53%) of the blast population at relapse. Early and accurate detection of subclonal variants with low allele frequency thus offers the opportunity for early intervention, prior to detection of clinical relapse, to improve disease outcome and enhance patient survival.

Genome-wide detection of ultralow-frequency substitution mutations in cultures of mouse lymphoma L5178Y cells and Caenorhabditis elegans worms by PacBio sequencing

Many conventional genetic toxicology assays require specialized cell cultures or animals and can only detect mutations that inactivate the function of a reporter gene. These limitations make such assays incompatible with many toxicological models but could be overcome by the development of techniques capable of directly detecting genome-wide somatic mutations through DNA sequencing. PacBio sequencing can generate almost error-free consensus reads by repeatedly inspecting both DNA strands from circularized molecules (a method known as PacBio HiFi). In this study, we show that PacBio HiFi can detect genome-wide ultralow-frequency substitution mutations in cultures of mouse lymphoma L5178Y cells and Caenorhabditis elegans worms. The mutation frequencies (MFs) of unexposed samples in both models were ~1 × 10^-7 mutations per base pair. Compared to these controls, PacBio HiFi detected MF increases of 23-fold in cultures of L5178Y cells exposed to 5 mM ethyl methanosulfonate (EMS) for 4 h, and 5-, 12-, and 29-fold in cultures of C. elegans worms exposed to 12.5, 25, and 50 mM EMS for 4 h, respectively. In both models, the mutation spectra of controls were diverse, while those derived from EMS-exposed samples were dominated by C:G → T:A transitions. To validate these results, clone sequencing analyses were performed on the same cultures of L5178Y cells. The results obtained by clone sequencing and PacBio HiFi were almost identical. Our results suggest that PacBio sequencing could be used for the detection, quantitation, and characterization of mutations in any DNA-containing sample, including those that are not compatible with conventional mutation detection approaches.

Human genetic risk of treatment with antiviral nucleoside analog drugs that induce lethal mutagenesis: The special case of molnupiravir

This review considers antiviral nucleoside analog drugs, including ribavirin, favipiravir, and molnupiravir, which induce genome error catastrophe in SARS-CoV or SARS-CoV-2 via lethal mutagenesis as a mode of action. In vitro data indicate that molnupiravir may be 100 times more potent as an antiviral agent than ribavirin or favipiravir. Molnupiravir has recently demonstrated efficacy in a phase 3 clinical trial. Because of its anticipated global use, its relative potency, and the reported in vitro "host" cell mutagenicity of its active principle, β-d-N4-hydroxycytidine, we have reviewed the development of molnupiravir and its genotoxicity safety evaluation, as well as the genotoxicity profiles of three congeners, that is, ribavirin, favipiravir, and 5-(2-chloroethyl)-2'-deoxyuridine. We consider the potential genetic risks of molnupiravir on the basis of all available information and focus on the need for additional human genotoxicity data and follow-up in patients treated with molnupiravir and similar drugs. Such human data are especially relevant for antiviral NAs that have the potential of permanently modifying the genomes of treated patients and/or causing human teratogenicity or embryotoxicity. We conclude that the results of preclinical genotoxicity studies and phase 1 human clinical safety, tolerability, and pharmacokinetics are critical components of drug safety assessments and sentinels of unanticipated adverse health effects. We provide our rationale for performing more thorough genotoxicity testing prior to and within phase 1 clinical trials, including human PIG-A and error corrected next generation sequencing (duplex sequencing) studies in DNA and mitochondrial DNA of patients treated with antiviral NAs that induce genome error catastrophe via lethal mutagenesis.

Signatures Beyond Oncogenic Mutations in Cell-Free DNA Sequencing for Non-Invasive, Early Detection of Cancer

Early detection of cancer saves lives, but an effective detection strategy in public health settings requires a delicate balance - periodic screening should neither miss rapidly progressing disease nor fail to detect rare tumors at unusual locations; on the other hand, even a modest false positive rate carries risks of over-diagnosis and over-treatment of relatively indolent non-malignant disease. Genomic profiling of cell-free DNA from liquid biopsy using massively parallel sequencing is emerging as an attractive, non-invasive screening platform for sensitive detection of multiple types of cancer in a single assay. Genomic data from cell-free DNA can not only identify oncogenic mutation status, but also additional molecular signatures related to potential tissue of origin, the extent of clonal growth, and malignant disease states. Utilization of the full potential of the molecular signatures from cfDNA sequencing data can guide clinical management strategies for targeted follow-ups using imaging or molecular marker-based diagnostic platforms and treatment options.

PacBio sequencing detects genome-wide ultra-low-frequency substitution mutations resulting from exposure to chemical mutagens

Genetic toxicology uses several assays to identity mutagens and protects the public. Most of these assays, however, rely on reporter genes, can only measure mutation indirectly based on phenotype, and often require specific cell lines or animal models-features that impede their integration with existing and emerging toxicological models, such as organoids. In this study, we show that PacBio Single-Molecule, Real-Time (PB SMRT) sequencing identified substitution mutations caused by chemical mutagens in Escherichia coli by generating nearly error-free consensus reads after repeatedly inspecting both strands of circular DNA molecules. Using DNA from E. coli exposed to ethyl methanosulfonate (EMS) or N-ethyl-N-nitrosourea (ENU), PB SMRT sequencing detected mutation frequencies (MFs) and spectra comparable to those obtained by clone-sequencing from the same exposures. The optimized background MF of PB SMRT sequencing was ≤ 1 × 10^-7 mutations per base pair (mut/bp).

Mutational signatures in T-lymphocytes of rats treated with N-propyl-N-nitrosourea and procarbazine

We have used whole genome sequencing (WGS) to determine mutational signatures induced in the T-cells of rats treated in vivo with N-propyl-N-nitrosourea (PNU) or procarbazine (PCZ). The signatures from the treated rats were different from the signature of background mutations. The main component of the spontaneous T-cell mutational signature was C➔T transition with all other single base substitutions evenly distributed. The PNU-induced mutational signature showed relatively equal contributions from C➔T and T➔C transitions, and T➔A transversions. The PCZ-induced signature was characterized by T➔C transitions, T➔A and, to a smaller extent, T➔G transversions. C➔G transversions were infrequent in either the PNU or PCZ signatures. WGS not only allowed mutational signature detection, but also measured quantitative responses to mutagen treatment: 10-40× increases in the number of mutations per clone were detected in T-cell clones from treated rats. The overall strand specificity of induced mutations for annotated rat genes was comparable to the strand specificity of mutations determined previously for the endogenous X-linked Pig-a gene. Our results provide valuable reference data for future applications of WGS in safety research and risk assessment.

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.