Chemically induced mutations in a MutaMouse reporter gene inform mechanisms underlying human cancer mutational signatures

Duplex sequencing identifies unique characteristics of ENU-induced mutations in male mouse germ cells†

Germ cell mutagenicity testing is increasingly required for chemical risk assessment. Duplex sequencing is rapidly gaining acceptance as a method to assess in vivo mutagenesis, and as a valid alternative to transgenic rodent mutation models such as the MutaMouse. We used a duplex sequencing panel of 20 genomic targets and the transgenic rodent assay to measure mutations in the germ cells of MutaMouse males exposed to 0, 1, 2, or 5 mg/kg N-ethyl-N-nitrosourea for 28 days. Germ cells from the seminiferous tubules were collected 28 days post-exposure. The transgenic rodent assay showed a significant increase in mutant frequencies at the high (P < 0.001) and medium (P = 0.01) N-ethyl-N-nitrosourea doses relative to controls, while duplex sequencing revealed a significant increase (P < 0.001) in N-ethyl-N-nitrosourea-induced mutations only at the high dose. Duplex sequencing mutation frequencies were lower in genic than in intergenic targets, suggesting a protective role for transcription-coupled repair. Interestingly, we observed several unique germ cell characteristics with respect to duplex sequencing data from rodent somatic tissues: 1) larger inter-animal variability in clonally expanded mutations that affects the ability to detect significant increases in mutation frequency; 2) a target on chromosome 2 showing much higher susceptibility to spontaneous and chemical-induced mutagenesis than other targets; and 3) a mutation spectrum consistent with that observed in the offspring of N-ethyl-N-nitrosourea-treated males but not with the spectrum in bone marrow of directly-exposed males. These results suggest that duplex sequencing is a promising approach for characterizing germ cell mutagenesis and that mutagenic mechanisms operating in germ cells differ from those in somatic tissues.

Exploring the Interconnections Between Mitochondrial Dysfunction and Polycystic Ovary Syndrome: A Comprehensive Integrated Analysis

Polycystic ovary syndrome (PCOS) is a leading cause of anovulatory infertility and is strongly linked to mitochondrial dysfunction (MD) in reproductive-age women. MD contributes to excessive reactive oxygen species (ROS) accumulation, exacerbating disease progression. This study aimed to identify key MD-related genes (MDRGs) involved in PCOS through bioinformatics analyses and experimental validation. Two PCOS transcriptome datasets (GSE34526 and GSE5850) were analyzed to identify differentially expressed genes (DEGs), which were then intersected with MDRGs to obtain MD-related DEGs (MDDEGs). Functional enrichment (GO, KEGG, GSEA) and protein-protein interaction (PPI) network analyses identified eight hub MDDEGs (MMP9, PPP1 CA, PSMD12, LIFR, PRKAA1, ITGAM, SUCLA2, GPBAR1). A rat PCOS model was established to validate hub gene expression via RT-qPCR, western blotting, and immunohistochemistry. The experimental data confirmed that seven hub genes exhibited consistent expression patterns with GSE34526 (P < 0.05), while only PRKAA1 and LIFR matched GSE5850 findings. Additionally, ROC analysis for the five most significant genes (LIFR, PBK, PRKAA1, RCAN1, MMP9) demonstrated promising diagnostic value (AUC > 0.85). This study highlights the role of MD in shaping the immune microenvironment of PCOS and identifies novel molecular targets for potential therapeutic interventions.

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.

The transgenic MutaMouse hepatocyte mutation assay in vitro: Mutagenicity and mutation spectra of six substances with different mutagenic mechanisms

Mutagenicity testing is a component of the hazard assessment of industrial chemicals, biocides, and pesticides. Mutations induced by test substances can be detected by in vitro and in vivo methods that have been adopted as OECD Test Guidelines. One of these in vivo methods is the Transgenic Rodent Assay (TGRA), OECD test guideline no. 488. An analogous in vitro TGRA has been described, but experience with this test method is limited. In this study, six in vivo TGRA positive mutagens were tested in the in vitro TGRA based on primary MutaMouse hepatocytes. In addition to the functional read-out of the lacZ reporter gene, induced mutations were analysed by next-generation sequencing (NGS). Five of the six in vivo TGRA positive mutagens (N-ethyl-N-nitrosourea (ENU), ethyl methanesulfonate (EMS), mitomycin C (MMC), benzo[a]pyrene (B[a]P), and azathioprine (AZA), but not cyproterone acetate) mutated the lacZ gene in vitro. NGS identified mutations which matched the mutagenic mechanisms described in the literature. The alkylating agent ENU induced a greater proportion of A:T to T:A transversions than did the other alkylating agent, EMS, whereas EMS increased smaller deletions (1-4 bp). G:C to T:A transversions accounted for the majority of mutations identified after treatments with MMC and B[a]P, both of which form monoadducts at the guanine N2 position. AZA induced mainly G:C to A:T transitions, explained by the structural similarity of one of its metabolites to guanine. An increased proportion of mid-size changes (0.3-2.5 kb) was detected only for the crosslinking mutagen MMC. The in vitro TGRA based on primary MutaMouse hepatocytes is a promising in vitro assay for the assessment of mutation induction, reflecting many aspects of the corresponding in vivo TGRA and allowing for mutation spectra analysis to evaluate the induced mutations.

AI and ML-based risk assessment of chemicals: predicting carcinogenic risk from chemical-induced genomic instability

Chemical risk assessment plays a pivotal role in safeguarding public health and environmental safety by evaluating the potential hazards and risks associated with chemical exposures. In recent years, the convergence of artificial intelligence (AI), machine learning (ML), and omics technologies has revolutionized the field of chemical risk assessment, offering new insights into toxicity mechanisms, predictive modeling, and risk management strategies. This perspective review explores the synergistic potential of AI/ML and omics in deciphering clastogen-induced genomic instability for carcinogenic risk prediction. We provide an overview of key findings, challenges, and opportunities in integrating AI/ML and omics technologies for chemical risk assessment, highlighting successful applications and case studies across diverse sectors. From predicting genotoxicity and mutagenicity to elucidating molecular pathways underlying carcinogenesis, integrative approaches offer a comprehensive framework for understanding chemical exposures and mitigating associated health risks. Future perspectives for advancing chemical risk assessment and cancer prevention through data integration, advanced machine learning techniques, translational research, and policy implementation are discussed. By implementing the predictive capabilities of AI/ML and omics technologies, researchers and policymakers can enhance public health protection, inform regulatory decisions, and promote sustainable development for a healthier future.

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.

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.

N-Nitrosodimethylamine investigations in Muta™Mouse define point-of-departure values and demonstrate less-than-additive somatic mutant frequency accumulations

The N-nitrosamine, N-nitrosodimethylamine (NDMA), is an environmental mutagen and rodent carcinogen. Small levels of NDMA have been identified as an impurity in some commonly used drugs, resulting in several product recalls. In this study, NDMA was evaluated in an OECD TG-488 compliant Muta™Mouse gene mutation assay (28-day oral dosing across seven daily doses of 0.02-4 mg/kg/day) using an integrated design that assessed mutation at the transgenic lacZ locus in various tissues and at the endogenous Pig-a gene-locus, along with micronucleus frequencies in peripheral blood. Liver pathology was determined together with NDMA exposure in blood and liver. The additivity of mutation induction was assessed by including two acute single-dose treatment groups (i.e. 5 and 10 mg/kg dose on Day 1), which represented the same total dose as two of the repeat dose treatment groups. NDMA did not induce statistically significant increases in mean lacZ mutant frequency (MF) in bone marrow, spleen, bladder, or stomach, nor in peripheral blood (Pig-a mutation or micronucleus induction) when tested up to 4 mg/kg/day. There were dose-dependent increases in mean lacZ MF in the liver, lung, and kidney following 28-day repeat dosing or in the liver and kidney after a single dose (10 mg/kg). No observed genotoxic effect levels (NOGEL) were determined for the positive repeat dose-response relationships. Mutagenicity did not exhibit simple additivity in the liver since there was a reduction in MF following NDMA repeat dosing compared with acute dosing for the same total dose. Benchmark dose modelling was used to estimate point of departure doses for NDMA mutagenicity in Muta™Mouse and rank order target organ tissue sensitivity (liver > kidney or lung). The BMD50 value for liver was 0.32 mg/kg/day following repeat dosing (confidence interval 0.21-0.46 mg/kg/day). In addition, liver toxicity was observed at doses of ≥ 1.1 mg/kg/day NDMA and correlated with systemic and target organ exposure. The integration of these results and their implications for risk assessment are discussed.

Somatic mutation patterns at Ig and Non-Ig Loci

The reverse transcriptase (RT) model of immunoglobulin (Ig) somatic hypermutation (SHM) has received insufficient scientific attention. This is understandable given that DNA deamination mediated by activation-induced deaminase (AID), the initiating step of Ig SHM, has dominated experiments since 2002. We summarise some key history of the RT Ig SHM model dating to 1987. For example, it is now established that DNA polymerase η, the sole DNA repair polymerase involved in post-replication short-patch repair, is an efficient cellular RT. This implies that it is potentially able to initiate target site reverse transcription by RNA-directed DNA repair at AID-induced lesions. Recently, DNA polymerase θ has also been shown to be an efficient cellular RT. Since DNA polymerase θ plays no significant role in Ig SHM, it could serve a similar RNA-dependent DNA polymerase role as DNA polymerase η at non-Ig loci in the putative RNA-templated nucleotide excision repair of bulky adducts and other mutagenic lesions on the transcribed strand. A major yet still poorly recognised consequence of the proposed RT process in Ig SHM is the generation of significant and characteristic strand-biased mutation signatures at both deoxyadenosine/deoxythymidine and deoxyguanosine/deoxycytidine base pairs. In this historical perspective, we highlight how diagnostic strand-biased mutation signatures are detected in vivo during SHM at both Ig loci in germinal centre B lymphocytes and non-Ig loci in cancer genomes. These strand-biased signatures have been significantly obscured by technical issues created by improper use of the polymerase chain reaction technique. A heightened awareness of this fact should contribute to better data interpretation and somatic mutation pattern recognition both at Ig and non-Ig loci.

The functional mutational landscape of the lacZ gene

The lacZ gene of Escherichia coli encodes β-galactosidase (β-gal), a lactose metabolism enzyme of the lactose operon. Previous chemical modification or site-directed mutagenesis experiments have identified 21 amino acids that are essential for β-gal catalytic activity. We have assembled over 10,000 lacZ mutations from published studies that were collected using a positive selection assay to identify mutations in lacZ that disrupted β-gal function. We analyzed 6,465 independent lacZ mutations that resulted in 2,732 missense mutations that impaired β-gal function. Those mutations affected 492 of the 1,023 lacZ codons, including most of the 21 previously known residues critical for catalytic activity. Most missense mutations occurred near the catalytic site and in regions important for subunit tetramerization. Overall, our work provides a comprehensive and detailed map of the amino acid residues affecting the structure and catalytic activity of the β-gal enzyme.

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 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.

Patterns of Somatic Variants in Colorectal Adenoma and Carcinoma Tissue and Matched Plasma Samples from the Hungarian Oncogenome Program

Analysis of circulating cell-free DNA (cfDNA) of colorectal adenoma (AD) and cancer (CRC) patients provides a minimally invasive approach that is able to explore genetic alterations. It is unknown whether there are specific genetic variants that could explain the high prevalence of CRC in Hungary. Whole-exome sequencing (WES) was performed on colon tissues (27 AD, 51 CRC) and matched cfDNAs (17 AD, 33 CRC); furthermore, targeted panel sequencing was performed on a subset of cfDNA samples. The most frequently mutated genes were APC, KRAS, and FBN3 in AD, while APC, TP53, TTN, and KRAS were the most frequently mutated in CRC tissue. Variants in KRAS codons 12 (AD: 8/27, CRC: 11/51 (0.216)) and 13 (CRC: 3/51 (0.06)) were the most frequent in our sample set, with G12V (5/27) dominance in ADs and G12D (5/51 (0.098)) in CRCs. In terms of the cfDNA WES results, tumor somatic variants were found in 6/33 of CRC cases. Panel sequencing revealed somatic variants in 8 out of the 12 enrolled patients, identifying 12/20 tumor somatic variants falling on its targeted regions, while WES recovered only 20% in the respective regions in cfDNA of the same patients. In liquid biopsy analyses, WES is less efficient compared to the targeted panel sequencing with a higher coverage depth that can hold a relevant clinical potential to be applied in everyday practice in the future.

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.

Cascade synthetic strategies opening access to medicinal-relevant aliphatic 3- and 4-membered N-heterocyclic scaffolds

Cascade reactions are often 'employed' by nature to construct structurally diverse nitrogen-containing heterocycles in a highly stereoselective fashion, i.e., secondary metabolites important for pharmacy. Nitrogen-containing heterocycles of three- and four-membered rings, as standalone and bicyclic compounds, inhibit different enzymes and are pharmacophores of approved drugs or drug candidates considered in many therapies, e.g. anticancer, antibacterial or antiviral. Domino transformations are in most cases in line with modern green chemistry concepts due to atom economy, one-pot procedures often without use the protective groups, time-saving and at markedly lower costs than multistep transformations. The tandem approaches can help to obtain novel N-heterocyclic scaffolds, functionalized according to structural requirements of the target in cells, taking into account the nature of functional group and stereochemistry. On the other hand cascade strategies allow to modify small N-heterocyclic rings in a systematic way, which is beneficial for structure-activity relationship (SAR) analyses. This review is focused on the biological relevance of the N-heterocyclic scaffolds with smaller 3- and 4-membered rings among approved drugs and leading structures of drug candidates. The cascade synthetic strategies offering N-heterocyclic scaffolds, at relatively good yields and high stereoselectivity, are discussed here. The review covers mainly years from 2015 to 2021.

Folate Intake Alters Mutation Frequency and Profiles in a Tissue- and Dose-Specific Manner in MutaMouse Male Mice

BACKGROUND

While cancer is common, its incidence varies widely by tissue. These differences are attributable to variable risk factors, such as environmental exposure, genetic inheritance, and lifetime number of stem cell divisions in a tissue. Folate deficiency is generally associated with increased risk for colorectal cancer (CRC) and acute lymphocytic leukemia (ALL). Conversely, high folic acid (FA) intake has also been associated with higher CRC risk.

OBJECTIVE

Our objective was to compare the effect of folate intake on mutant frequency (MF) and types of mutations in the colon and bone marrow of mice.

METHODS

Five-week-old MutaMouse male mice were fed a deficient (0 mg FA/kg), control (2 mg FA/kg), or supplemented (8 mg FA/kg) diet for 20 wk. Tissue MF was assessed using the lacZ mutant assay and comparisons made by 2-factor ANOVA. LacZ mutant plaques were sequenced using next-generation sequencing, and diet-specific mutation profiles within each tissue were compared by Fisher's exact test.

RESULTS

In the colon, the MF was 1.5-fold and 1.3-fold higher in mice fed the supplemented diet compared with mice fed the control (P = 0.001) and deficient (P = 0.008) diets, respectively. This contrasted with the bone marrow MF in the same mice where the MF was 1.7-fold and 1.6-fold higher in mice fed the deficient diet compared with mice fed the control (P = 0.02) and supplemented (P = 0.03) diets, respectively. Mutation profiles and signatures (mutation context) were tissue-specific.

CONCLUSIONS

Our data indicate that dietary folate intake affects mutagenesis in a tissue- and dose-specific manner in mice. Mutation profiles were generally tissue- but not dose-specific, suggesting that altered cellular folate status appears to interact with endogenous mutagenic mechanisms in each tissue to create a permissive context in which specific mutation types accumulate. These data illuminate potential mechanisms underpinning differences in observed associations between folate intake/status and cancer.

The 28 + 28 day design is an effective sampling time for analyzing mutant frequencies in rapidly proliferating tissues of MutaMouse animals

The Organisation for Economic Co-Operation and Development Test Guideline 488 (TG 488) uses transgenic rodent models to generate in vivo mutagenesis data for regulatory submission. The recommended design in TG 488, 28 consecutive daily exposures with tissue sampling three days later (28 + 3d), is optimized for rapidly proliferating tissues such as bone marrow (BM). A sampling time of 28 days (28 + 28d) is considered more appropriate for slowly proliferating tissues (e.g., liver) and male germ cells. We evaluated the impact of the sampling time on mutant frequencies (MF) in the BM of MutaMouse males exposed for 28 days to benzo[a]pyrene (BaP), procarbazine (PRC), isopropyl methanesulfonate (iPMS), or triethylenemelamine (TEM) in dose-response studies. BM samples were collected + 3d, + 28d, + 42d or + 70d post exposure and MF quantified using the lacZ assay. All chemicals significantly increased MF with maximum fold increases at 28 + 3d of 162.9, 6.6, 4.7 and 2.8 for BaP, PRC, iPMS and TEM, respectively. MF were relatively stable over the time period investigated, although they were significantly increased only at 28 + 3d and 28 + 28d for TEM. Benchmark dose (BMD) modelling generated overlapping BMD confidence intervals among the four sampling times for each chemical. These results demonstrate that the sampling time does not affect the detection of mutations for strong mutagens. However, for mutagens that produce small increases in MF, sampling times greater than 28 days may produce false-negative results. Thus, the 28 + 28d protocol represents a unifying protocol for simultaneously assessing mutations in rapidly and slowly proliferating somatic tissues and male germ cells.

Direct quantification of in vivo mutagenesis and carcinogenesis using duplex sequencing

Error-corrected next-generation sequencing (ecNGS) can be used to rapidly detect and quantify the in vivo mutagenic impact of environmental exposures or endogenous processes in any tissue, from any species, at any genomic location. The greater speed, higher scalability, richer data outputs, and cross-species and cross-locus applicability of ecNGS compared to existing methods make it a powerful new tool for mutational research, regulatory safety testing, and emerging clinical applications.

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