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

The Distinct Roles of NEIL1 and XPA in Limiting Aflatoxin B1-Induced Mutagenesis in Mice

Dietary exposure to aflatoxin B1 (AFB1) is a risk factor for the development of hepatocellular carcinomas. Following metabolic activation, AFB1 reacts with guanines to form covalent DNA adducts, which induce high-frequency G > T transversions. The molecular signature associated with these mutational events aligns with the single-base substitution signature 24 (SBS24) in the Catalog of Somatic Mutations in Cancer database. Deficiencies in either base excision repair due to the absence of Nei-like DNA glycosylase 1 (NEIL1) or nucleotide excision repair due to the absence of xeroderma complementation group A protein (XPA) contribute to hepatocellular carcinomas in murine models. In the current study, ultra-low error duplex sequencing was used to characterize mutational profiles in liver DNAs of NEIL1-deficient, XPA-deficient, and DNA repair-proficient mice following neonatal injection of 1 mg/kg AFB1. Analyses of AFB1-induced mutations showed high cosine similarity to SBS24 regardless of repair proficiency status. The absence of NEIL1 resulted in an approximately 30% increase in the frequency of mutations, with the distribution suggesting preferential NEIL1-dependent repair of AFB1 lesions in open chromatin regions. A trend of increased mutagenesis was also observed in the absence of XPA. Consistent with the role of XPA in transcription-coupled repair, mutational profiles in XPA-deficient mice showed disruption of the transcriptional bias in mutations associated with SBS24. Implications: Our findings define the roles of DNA repair pathways in AFB1-induced mutagenesis and carcinogenesis in murine models, with these findings having implications in human health for those with base excision repair and nucleotide excision repair deficiencies.

Dose-Related Mutagenic and Clastogenic Effects of Benzo[b]fluoranthene in Mouse Somatic Tissues Detected by Duplex Sequencing and the Micronucleus Assay

Polycyclic aromatic hydrocarbons (PAHs) are common environmental pollutants that originate from the incomplete combustion of organic materials. We investigated the clastogenicity and mutagenicity of benzo[b]fluoranthene (BbF), one of 16 priority PAHs, in MutaMouse males after a 28 day oral exposure. BbF causes robust dose-dependent increases in micronucleus frequency in peripheral blood, indicative of chromosome damage. Duplex sequencing (DS), an error-corrected sequencing technology, reveals that BbF induces dose-dependent increases in mutation frequencies in bone marrow (BM) and liver. Mutagenicity is increased in intergenic relative to genic regions, suggesting a role for transcription-coupled repair of BbF-induced DNA damage. At higher doses, the maximum mutagenic response to BbF is higher in liver, which has a lower mitotic index but higher metabolic capacity than BM; however, mutagenic potency is comparable between the two tissues. BbF induces primarily C:G > A:T mutations, followed by C:G > T:A and C:G > G:C, indicating that BbF metabolites mainly target guanines and cytosines. The mutation spectrum of BbF correlates with cancer mutational signatures associated with tobacco exposure, supporting its contribution to the carcinogenicity of combustion-derived PAHs in humans. Overall, BbF's mutagenic effects are similar to benzo[a]pyrene, a well-studied mutagenic PAH. Our work showcases the utility of DS for effective mutagenicity assessment of environmental pollutants.

Effect of sequencing platforms on the sensitivity of chemical mutation detection using Hawk-Seq™

BACKGROUND

Error-corrected next-generation sequencing (ecNGS) technologies have enabled the direct evaluation of genome-wide mutations after exposure to mutagens. Previously, we reported an ecNGS methodology, Hawk-Seq™, and demonstrated its utility in evaluating mutagenicity. The evaluation of technical transferability is essential to further evaluate the reliability of ecNGS-based assays. However, cutting-edge sequencing platforms are continually evolving, which can affect the sensitivity of ecNGS. Therefore, the effect of differences in sequencing instruments on mutation data quality should be evaluated.

RESULTS

We assessed the performance of four sequencing platforms (HiSeq2500, NovaSeq6000, NextSeq2000, and DNBSEQ-G400) with the Hawk-Seq™ protocol for mutagenicity evaluation using DNA samples from mouse bone marrow exposed to benzo[a]pyrene (BP). The overall mutation (OM) frequencies per 106 bp in vehicle-treated samples were 0.22, 0.36, 0.46, and 0.26 for HiSeq2500, NovaSeq6000, NextSeq2000, and DNBSEQ-G400, respectively. The OM frequency of NextSeq2000 was significantly higher than that of HiSeq2500, suggesting the difference to be based on the platform. The relatively higher value in NextSeq2000 was a consequence of the G:C to C:G mutations in NextSeq2000 data (0.67 per 106 G:C bp), which was higher than the mean of the four platforms by a ca. of 0.25 per 106 G:C bp. A clear dose-dependent increase in G:C to T:A mutation frequencies was observed in all four sequencing platforms after BP exposure. The cosine similarity values of the 96-dimensional trinucleotide mutation patterns between HiSeq and the three other platforms were 0.93, 0.95, and 0.92 for NovaSeq, NextSeq, and DNBSeq, respectively. These results suggest that all platforms can provide equivalent data that reflect the characteristics of the mutagens.

CONCLUSIONS

All platforms sensitively detected mutagen-induced mutations using the Hawk-Seq™ analysis. The substitution types and frequencies of the background errors differed depending on the platform. The effects of sequencing platforms on mutagenicity evaluation should be assessed before experimentation.

Frequency and spectrum of mutations in human sperm measured using duplex sequencing correlate with trio-based de novo mutation analyses

De novo mutations (DNMs) are drivers of genetic disorders. However, the study of DNMs is hampered by technological limitations preventing accurate quantification of ultra-rare mutations. Duplex Sequencing (DS) theoretically has < 1 error/billion base-pairs (bp). To determine the DS utility to quantify and characterize DNMs, we analyzed DNA from blood and spermatozoa from six healthy, 18-year-old Swedish men using the TwinStrand DS mutagenesis panel (48 kb spanning 20 genic and intergenic loci). The mean single nucleotide variant mutation frequency (MF) was 1.2 × 10- 7 per bp in blood and 2.5 × 10- 8 per bp in sperm, with the most common base substitution being C > T. Blood MF and substitution spectrum were similar to those reported in blood cells with an orthogonal method. The sperm MF was in the same order of magnitude and had a strikingly similar spectrum to DNMs from publicly available whole genome sequencing data from human pedigrees (1.2 × 10- 8 per bp). DS revealed much larger numbers of insertions and deletions in sperm over blood, driven by an abundance of putative extra-chromosomal circular DNAs. The study indicates the strong potential of DS to characterize human DNMs to inform factors that contribute to disease susceptibility and heritable genetic risks.

Power analyses to inform Duplex Sequencing study designs for MutaMouse liver and bone marrow

Regulatory genetic toxicology testing is essential for identifying potentially mutagenic hazards. Duplex Sequencing (DS) is an error-corrected next-generation sequencing technology that provides substantial advantages for mutation analysis over conventional mutagenicity assays including: improved accuracy of mutation detection, ability to measure changes in mutation spectrum, and applicability across diverse biological models. To apply DS for regulatory toxicology testing, power analyses are required to determine suitable sample sizes and study designs. In this study, we explored study designs to achieve sufficient power for various effect sizes in chemical mutagenicity assessment. We collected data from MutaMouse bone marrow and liver samples that were analyzed by DS using TwinStrand's Mouse Mutagenesis Panel. Average duplex reads achieved in two separates studies on liver and bone marrow were 8.4 × 108 (± 7.4 × 107) and 9.5 × 108 (± 1.0 × 108), respectively. Baseline mean mutation frequencies (MF) were 4.6 × 10-8 (± 6.7 × 10-9) and 4.6 × 10-8 (± 1.1 × 10-8), with estimated standard deviations for the animal-to-animal random effect of 0.15 and 0.20, for liver and bone marrow, respectively. We conducted simulation analyses based on these empirically derived parameters. We found that a sample size of four animals per group is sufficient to obtain over 80% power to detect a two-fold change in MF relative to baseline. In addition, we estimated the minimal total number of informative duplex bases sequenced with different sample sizes required to retain power for various effect sizes. Our work provides foundational data for establishing suitable study designs for mutagenicity testing using DS.

CarcSeq detection of lorcaserin-induced clonal expansion of Pik3ca H1047R mutants in rat mammary tissue

Lorcaserin is a 5-hydroxytryptamine 2C (serotonin) receptor agonist and a nongenotoxic rat carcinogen, which induced mammary tumors in male and female rats in a 2-yr bioassay. Female Sprague Dawley rats were treated by gavage daily with 0, 30, or 100 mg/kg lorcaserin, replicating bioassay dosing but for shorter duration, 12 or 24 wk. To characterize exposure and eliminate possible confounding by a potentially genotoxic degradation product, lorcaserin and N-nitroso-lorcaserin were quantified in dosing solutions, terminal plasma, mammary, and liver samples using ultra-high-performance liquid chromatography-electrospray tandem mass spectrometry. N-nitroso-lorcaserin was not detected, supporting lorcaserin classification as nongenotoxic carcinogen. Mammary DNA samples (n = 6/dose/timepoint) were used to synthesize PCR products from gene segments encompassing hotspot cancer driver mutations, namely regions of Apc, Braf, Egfr, Hras, Kras, Nfe2l2, Pik3ca, Setbp1, Stk11, and Tp53. Mutant fractions (MFs) in the amplicons were quantified by CarcSeq, an error-corrected next-generation sequencing approach. Considering all recovered mutants, no significant differences between lorcaserin dose groups were observed. However, significant dose-responsive increases in Pik3ca H1047R mutation were observed at both timepoints (ANOVA, P < 0.05), with greater numbers of mutants and mutants with greater MFs observed at 24 wk as compared with 12 wk. These observations suggest lorcaserin promotes outgrowth of spontaneously occurring Pik3ca H1047R mutant clones leading to mammary carcinogenesis. Importantly, this work reports approaches to analyze clonal expansion and demonstrates CarcSeq detection of the carcinogenic impact (selective Pik3ca H0147R mutant expansion) of a nongenotoxic carcinogen using a treatment duration as short as 3 months.

Dose-related Mutagenic and Clastogenic Effects of Benzo[b]fluoranthene in Mouse Somatic Tissues Detected by Duplex Sequencing and the Micronucleus Assay

Polycyclic aromatic hydrocarbons (PAHs) are common environmental pollutants that originate from the incomplete combustion of organic materials. We investigated the clastogenicity and mutagenicity of benzo[ b ]fluoranthene (BbF), one of 16 priority PAHs, in MutaMouse males after a 28-day oral exposure. BbF causes robust dose-dependent increases in micronucleus frequency in peripheral blood, indicative of chromosome damage. Duplex Sequencing (DS), an error-corrected sequencing technology, reveals that BbF induces dose-dependent increases in mutation frequencies in bone marrow (BM) and liver. Mutagenicity is increased in intergenic relative to genic regions, suggesting a role for transcription-coupled repair of BbF-induced DNA damage. At higher doses, the maximum mutagenic response to BbF is higher in liver, which has a lower mitotic index but higher metabolic capacity than BM; however, mutagenic potency is comparable between the two tissues. BbF induces primarily C:G>A:T mutations, followed by C:G>T:A and C:G>G:C, indicating that BbF metabolites mainly target guanines and cytosines. The mutation spectrum of BbF correlates with cancer mutational signatures associated with tobacco exposure, supporting its contribution to the carcinogenicity of combustion-derived PAHs in humans. Overall, BbF's mutagenic effects are similar to benzo[ a ]pyrene, a well-studied mutagenic PAH. Our work showcases the utility of DS for effective mutagenicity assessment of environmental pollutants.

Synopsis

We used Duplex Sequencing to study the mutagenicity of benzo[ b ]fluoranthene across the mouse genome. Dose-dependent changes in mutation frequency and spectrum quantify its role in PAH-induced carcinogenicity.

Assessing the genotoxicity of N-nitrosodiethylamine with three in vivo endpoints in male Big Blue® transgenic and wild-type C57BL/6N mice

The detection of N-nitrosamines in drug products has raised global regulatory interest in recent years due to the carcinogenic potential of some nitrosamines in animals and a need to identify a testing strategy has emerged. Ideally, methods used would allow for the use of quantitative analysis of dose-response data from in vivo genotoxicity assays to determine a compound-specific acceptable intake for novel nitrosamines without sufficient carcinogenicity data. In a previous study we compared the dose-response relationships of N-nitrosodiethylamine (NDEA) in three in vivo genotoxicity endpoints in rats. Here we report a comparison of NDEA's genotoxicity profile in mice. Big Blue® mice were administered NDEA at doses of 0.001, 0.01, 0.1, 1 and 3 mg/kg/day by oral gavage for 28 days followed by 3 days of expression. Statistically significant increases in the NDEA induced mutations were detected by both the transgenic rodent mutation assay (TGR) using the cII endpoint and by duplex sequencing in the liver but not bone marrow of mice. In addition, administration of NDEA for two consecutive days in male C57BL/6N mice caused elevated DNA damage levels in the liver as measured by % tail DNA in comet assay. The benchmark dose (BMD) analysis shows a BMDL50 of 0.03, 0.04 and 0.72 mg/kg/day for TGR, duplex sequencing and comet endpoints, respectively. Overall, this study demonstrated a similar genotoxicity profile of NDEA between mice and rats and provides a reference that can be used to compare the potential potency of other novel nitrosamines for the induction of gene mutations.

Genome-scale mutational signature analysis in archived fixed tissues

Mutation spectra and mutational signatures in cancerous and non-cancerous tissues can be identified by various established techniques of massively parallel sequencing (or next-generation sequencing) including whole-exome or whole-genome sequencing, and more recently by error-corrected/duplex sequencing. One rather underexplored area has been the genome-scale analysis of mutational signatures as markers of mutagenic exposures, and their impact on cancer driver events applied to formalin-fixed or alcohol-fixed paraffin embedded archived biospecimens. This review showcases successful applications of the next-generation sequencing methodologies in archived fixed tissues, including the delineation of the specific tissue fixation-related DNA damage manifesting as artifactual signatures, distinguishable from the true signatures that arise from biological mutagenic processes. Overall, we discuss and demonstrate how next-generation sequencing techniques applied to archived fixed biospecimens can enhance our understanding of cancer causes including mutagenic effects of extrinsic cancer risk agents, and the implications for prevention efforts aimed at reducing avoidable cancer-causing exposures.

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.

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.

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.

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.

Strategic Assessment of Boron-Enriched Carbon Dots/Naproxen: Diagnostic, Toxicity, and In Vivo Therapeutic Evaluation

Cancer is a significant global public health concern, ranking as the leading cause of mortality worldwide. This study thoroughly explores boron-doped carbon dots (B-CDs) through a simple/rapid microwave-assisted approach and their versatile applications in cancer therapy. The result was highly uniform particles with an average diameter of approximately 4 nm. B-CDs exhibited notable properties, including strong fluorescence with a quantum yield of 33%. Colloid stability tests revealed their robustness within a pH range of 6-12, NaCl concentrations up to 0.5 M, and temperatures ranging from 30 to 60 °C. The study also delved into the kinetics of naproxen release from B-CDs as a drug delivery system. The loading efficacy of naproxen exceeded 55.56%. Under varying pH conditions, the release of naproxen from B-CDs conformed to the Peppas-Sahlin model, demonstrating the potential of Naproxen-loaded CDs for cancer drug delivery. In vitro cytotoxicity assessments, conducted using the CCK-8 Assay and flow cytometry, consistently indicated low toxicity with average cell viability exceeding 80%. An in vivo toxicity test on female mice administered 20 mg/kg of B-CDs for 31 days revealed reversible histological changes in the liver and kidneys, while the pancreas remained unaffected. Importantly, B-CDs did not impact the mice's physical behavior, body weight, or survival. In vivo experiments targeting benzo(a)pyrene-induced fibrosarcoma demonstrated the efficacy of B-CDs as naproxen carriers in the treatment of cancer. This in vivo study provides a thorough comprehension of B-CDs synthesis and toxicity and their potential applications in cancer therapy and drug delivery systems.

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.

Detection of in vivo mutagenicity in rat liver samples using error-corrected sequencing techniques

BACKGROUND

Mutagenicity, the ability of chemical agents to cause mutations and potentially lead to cancer, is a critical aspect of substance safety assessment for protecting human health and the environment. Metabolic enzymes activate multiple mutagens in living organisms, thus in vivo animal models provide highly important information for evaluating mutagenicity in human. Rats are considered suitable models as they share a similar metabolic pathway with humans for processing toxic chemical and exhibit higher responsiveness to chemical carcinogens than mice. To assess mutagenicity in rats, transgenic rodents (TGRs) are widely used for in vivo gene mutation assays. However, such assays are labor-intensive and could only detect transgene mutations inserted into the genome. Therefore, introducing a technology to directly detect in vivo mutagenicity in rats would be necessary. The next-generation sequencing (NGS) based error-corrected sequencing technique is a promising approach for such purposes.

RESULTS

We investigated the applicability of paired-end and complementary consensus sequencing (PECC-Seq), an error-corrected sequencing technique, for detecting in vivo mutagenicity in the rat liver samples. PECC-Seq allows for the direct detection of ultra-rare somatic mutations in the genomic DNA without being constrained by the genomic locus, tissue, or organism. We tested PECC-Seq feasibility in rats treated with diethylnitrosamine (DEN), a mutagenic compound. Interestingly, the mutation and mutant frequencies between PECC-Seq and the TGR assay displayed a promising correlation. Our results also demonstrated that PECC-Seq could successfully detect the A:T > T:A mutation in rat liver samples, consistent with the TGR assay. Furthermore, we calculated the trinucleotide mutation frequency and proved that PECC-Seq accurately identified the DEN treatment-induced mutational signatures.

CONCLUSIONS

Our study provides the first evidence of using PECC-Seq for in vivo mutagenicity detection in rat liver samples. This approach could provide a valuable alternative to conventional TGR assays as it is labor- and time-efficient and eliminates the need for transgenic rodents. Error-corrected sequencing techniques, such as PECC-Seq, represent promising approaches for enhancing mutagenicity assessment and advancing regulatory science.

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.

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.

Duplex sequencing provides detailed characterization of mutation frequencies and spectra in the bone marrow of MutaMouse males exposed to procarbazine hydrochloride

Mutagenicity testing is an essential component of health safety assessment. Duplex Sequencing (DS), an emerging high-accuracy DNA sequencing technology, may provide substantial advantages over conventional mutagenicity assays. DS could be used to eliminate reliance on standalone reporter assays and provide mechanistic information alongside mutation frequency (MF) data. However, the performance of DS must be thoroughly assessed before it can be routinely implemented for standard testing. We used DS to study spontaneous and procarbazine (PRC)-induced mutations in the bone marrow (BM) of MutaMouse males across a panel of 20 diverse genomic targets. Mice were exposed to 0, 6.25, 12.5, or 25 mg/kg-bw/day for 28 days by oral gavage and BM sampled 42 days post-exposure. Results were compared with those obtained using the conventional lacZ viral plaque assay on the same samples. DS detected significant increases in mutation frequencies and changes to mutation spectra at all PRC doses. Low intra-group variability within DS samples allowed for detection of increases at lower doses than the lacZ assay. While the lacZ assay initially yielded a higher fold-change in mutant frequency than DS, inclusion of clonal mutations in DS mutation frequencies reduced this discrepancy. Power analyses suggested that three animals per dose group and 500 million duplex base pairs per sample is sufficient to detect a 1.5-fold increase in mutations with > 80% power. Overall, we demonstrate several advantages of DS over classical mutagenicity assays and provide data to support efforts to identify optimal study designs for the application of DS as a regulatory test.

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.

Possible Genetic Risks from Heat-Damaged DNA in Food

The consumption of foods prepared at high temperatures has been associated with numerous health risks. To date, the chief identified source of risk has been small molecules produced in trace levels by cooking and reacting with healthy DNA upon consumption. Here, we considered whether the DNA in food itself also presents a hazard. We hypothesize that high-temperature cooking may cause significant damage to the DNA in food, and this damage might find its way into cellular DNA by metabolic salvage. We tested cooked and raw foods and found high levels of hydrolytic and oxidative damage to all four DNA bases upon cooking. Exposing cultured cells to damaged 2'-deoxynucleosides (particularly pyrimidines) resulted in elevated DNA damage and repair responses in the cells. Feeding a deaminated 2'-deoxynucleoside (2'-deoxyuridine), and DNA containing it, to mice resulted in substantial uptake into intestinal genomic DNA and promoted double-strand chromosomal breaks there. The results suggest the possibility of a previously unrecognized pathway whereby high-temperature cooking may contribute to genetic risks.

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.