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

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.

Assessment of in vivo chemical mutagenesis by long-read sequencing

Evaluating the mutagenic properties of chemicals is crucial for understanding their potential cancer risks. Recent Illumina-based error-corrected sequencing techniques have enabled the direct detection of mutations induced de novo by mutagens. However, as the Illumina platform lacks intrinsic error-correction capabilities, complex library preparations and bioinformatic processes are necessary to identify these rare mutations. In this study, we evaluated whether long-read PacBio-based HiFi sequencing (HiFi seq), which has integrated error-correction, can detect de novo mutations induced by mutagens in C57BL/6 mouse tissues. Using HiFi seq, dose-dependent increases in mutation frequencies were found in tissues from mice exposed to 7,12-dimethylbenz[a]anthracene, procarbazine, and N-propyl-N-nitrosourea. Furthermore, the mutational signatures derived from these exposures were consistent with those previously reported for these mutagens. This study demonstrates that HiFi seq can complement established mutation detection assays to facilitate the identification of hazardous compounds.

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.

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.

Repeat treatment of organotypic airway cultures with ethyl methanesulfonate causes accumulation of somatic cell mutations without expansion of bronchial-carcinoma-specific cancer driver mutations

The human in vitro organotypic air-liquid-interface (ALI) airway tissue model is structurally and functionally similar to the human large airway epithelium and, as a result, is being used increasingly for studying the toxicity of inhaled substances. Our previous research demonstrated that DNA damage and mutagenesis can be detected in human airway tissue models under conditions used to assess general and respiratory toxicity endpoints. Expanding upon our previous proof-of-principle study, human airway epithelial tissue models were treated with 6.25-100 µg/mL ethyl methanesulfonate (EMS) for 28 days, followed by a 28-day recovery period. Mutagenesis was evaluated by Duplex Sequencing (DS), and clonal expansion of bronchial-cancer-specific cancer-driver mutations (CDMs) was investigated by CarcSeq to determine if both mutation-based endpoints can be assessed in the same system. Additionally, DNA damage and tissue-specific responses were analyzed during the treatment and following the recovery period. EMS exposure led to time-dependent increases in mutagenesis over the 28-day treatment period, without expansion of clones containing CDMs; the mutation frequencies remained elevated following the recovery. EMS also produced an increase in DNA damage measured by the CometChip and MultiFlow assays and the elevated levels of DNA damage were reduced (but not eliminated) following the recovery period. Cytotoxicity and most tissue-function changes induced by EMS treatment recovered to control levels, the exception being reduced proliferating cell frequency. Our results indicate that general, respiratory-tissue-specific and genotoxicity endpoints increased with repeat EMS dosing; expansion of CDM clones, however, was not detected using this repeat treatment protocol. DISCLAIMER: This article reflects the views of its authors and does not necessarily reflect those of the U.S. Food and Drug Administration. Any mention of commercial products is for clarification only and is not intended as approval, endorsement, or recommendation.

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.

Evaluating the mutagenicity of N-nitrosodimethylamine in 2D and 3D HepaRG cell cultures using error-corrected next generation sequencing

Human liver-derived metabolically competent HepaRG cells have been successfully employed in both two-dimensional (2D) and 3D spheroid formats for performing the comet assay and micronucleus (MN) assay. In the present study, we have investigated expanding the genotoxicity endpoints evaluated in HepaRG cells by detecting mutagenesis using two error-corrected next generation sequencing (ecNGS) technologies, Duplex Sequencing (DS) and High-Fidelity (HiFi) Sequencing. Both HepaRG 2D cells and 3D spheroids were exposed for 72 h to N-nitrosodimethylamine (NDMA), followed by an additional incubation for the fixation of induced mutations. NDMA-induced DNA damage, chromosomal damage, and mutagenesis were determined using the comet assay, MN assay, and ecNGS, respectively. The 72-h treatment with NDMA resulted in concentration-dependent increases in cytotoxicity, DNA damage, MN formation, and mutation frequency in both 2D and 3D cultures, with greater responses observed in the 3D spheroids compared to 2D cells. The mutational spectrum analysis showed that NDMA induced predominantly A:T → G:C transitions, along with a lower frequency of G:C → A:T transitions, and exhibited a different trinucleotide signature relative to the negative control. These results demonstrate that the HepaRG 2D cells and 3D spheroid models can be used for mutagenesis assessment using both DS and HiFi Sequencing, with the caveat that severe cytotoxic concentrations should be avoided when conducting DS. With further validation, the HepaRG 2D/3D system may become a powerful human-based metabolically competent platform for genotoxicity testing.

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.

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.

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.