Mutational spectra for polycyclic aromatic hydrocarbons in the supF target gene

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.

Genome-wide somatic mutation analysis via Hawk-Seq™ reveals mutation profiles associated with chemical mutagens

It is difficult to identify mutagen-induced genome-wide somatic mutations using next generation sequencing; hence, mutagenic features of each mutagen and their roles in cancer development require further elucidation. We described Hawk-Seq™, a highly accurate genome sequencing method and the optimal conditions, for using it to construct libraries that would enable the accurate (c.a. 1 error/10⁷-10⁸ bp) and efficient survey of genome-wide mutations. Genomic mutations in gpt delta mice or Salmonella typhimurium TA100 exposed to methylnitrosourea (MNU), ethylnitrosourea (ENU), diethylnitrosamine (DEN), benzo[a]pyrene (BP), and aristolochic acid (AA) were profiled using Hawk-Seq™ to analyse positions, substitution patterns, or frequencies. The resultant vast mutation data provided high-resolution mutational signatures, including for minor mutational fractions (e.g. G:C>A:T by AA), which enabled the clarification of the mutagenic features of all mutagens. The 96-type mutational signatures of MNU, AA, and BP indicate their partial similarity to signature 11, 22, and 4 or 29, respectively. Meanwhile, signatures attributable to ENU and DEN were highly similar to each other, but not to signature 11, suggesting that the mechanisms of these agents differed from those of typical alkylating agents. Thus, Hawk-Seq™ can clarify genome-wide chemical mutagenicity profiles at extraordinary resolutions, thereby providing insight into mutagen mechanisms and their roles in cancer development.

A genome-wide mutation analysis method enabling high-throughput identification of chemical mutagen signatures

Trinucleotide mutational signatures extracted from cancer genomes provide clues useful in understanding the roles of mutagens and mutagenic mechanisms in cancer development. The lack of a simple method for genome-wide analysis of alterations induced by mutagens hampers the identification of trinucleotide signatures of mutagen exposure and evaluation of their relationships with human cancers. Here, we describe a novel approach to facilitate analysis of chemically induced mutations in bacterial cells by detection of increased frequencies of base substitutions after mutagen exposure, using paired-end overlapping next-generation sequencing. DNA samples from Salmonella typhimurium strain TA100, exposed to three alkylating agents, ethylnitrosourea (ENU), methylnitrosourea (MNU), and ethyl methansulphonate (EMS), were analysed. The G:C > A:T mutation frequency was increased in all samples, whereas A:T base pair substitution frequencies were increased specifically in samples exposed to ENU, consistent with previous reports. Mutation patterns in the context of 96 possible trinucleotide formats in these samples exhibited a sharp peak corresponding to an NpCpY consensus sequence, which is similar to the mutational signature of alkylating agents in human cancer. These results indicate that our approach can be useful in facilitating the understanding of mechanisms underlying chemical mutagenicity and for identification of unknown causal mutagens in human cancer.

Analysis and Identification of 2'-Deoxyadenosine-Derived Adducts in Lung and Liver DNA of F-344 Rats Treated with the Tobacco-Specific Carcinogen 4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone and Enantiomers of its Metabolite 4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanol

4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) and its metabolite 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanol (NNAL) are carcinogenic in animal models and are believed to play an important role in human lung carcinogenesis for cigarette smokers. Cytochrome P450-mediated metabolism of these tobacco-specific nitrosamines produces reactive species that alkylate DNA in the form of pyridyloxobutyl (POB)- or pyridylhydroxybutyl (PHB)-DNA adducts. Understanding the formation mechanism and overall levels of these adducts can potentially enhance cancer prevention methods through the identification of particularly susceptible smokers. Previous studies have identified and measured a panel of POB- and PHB-DNA base adducts of dGuo, dCyd, and Thd; however, dAdo adducts have yet to be determined. In this study, we complete this DNA adduct panel by identifying and quantifying levels of NNK- and NNAL-derived dAdo adducts in vitro and in vivo. To accomplish this, we synthesized standards for expected dAdo-derived DNA adducts and used isotope-dilution LC-ESI⁺-MS/MS to identify POB adducts formed in vitro from the reaction of 4-(acetoxymethylnitrosamino)-1-(3-pyridyl)-1-butanone (NNKOAc) with calf thymus DNA. Adduct levels were then quantified in lung and liver DNA of rats chronically treated with NNK or NNAL for 50 weeks using similar LC-MS detection methods. The in vitro studies identified N⁶-POB-dAdo and N¹-POB-dIno as products of the reaction of NNKOAc with DNA, which supports our proposed mechanism of formation. Though both N⁶-dAdo and N¹-dIno adducts were found in vitro, only N⁶-dAdo adducts were found in vivo, implying possible intervention by DNA repair mechanisms. Analogous to previous studies, levels of N⁶-POB-dAdo and N⁶-PHB-dAdo varied both with tissue and treatment type. Despite the adduct levels being relatively modest compared to most other POB- and PHB-DNA adducts, they may play a biological role and could be used in future studies as NNK- and NNAL-specific DNA damage biomarkers.

Limited mutagenicity of electronic cigarettes in mouse or human cells in vitro

OBJECTIVES

Electronic cigarettes (e-cig), which are promoted as safe alternatives to tobacco cigarettes or as aides to smoking cessation, are becoming increasingly popular among adult chronic smokers and adolescents experimenting with tobacco products. Despite the known presence of toxicants and carcinogens in e-cig liquid and vapor, the possible carcinogenic effects of e-cig use in humans are unknown.

MATERIALS AND METHODS

We have utilized two validated in vitro model systems to investigate whether e-cig vapor induces mutation in mouse or human cells. We have exposed transgenic mouse fibroblasts in vitro to e-cig vapor extracts prepared from three popular brands, and determined the induction of mutagenesis in a reporter gene, the cII transgene. Furthermore, we have treated the pSP189 plasmid with e-cig vapor extract, transfected human fibroblast cells with the e-cig-treated plasmid, and screened for the induced mutations in the supF gene.

RESULTS AND CONCLUSION

We observed no statistically significant increases in relative mutant frequency in the cII transgene or supF gene in the e-cig treated mouse or human cells, respectively. Our data indicate that e-cig vapor extracts from the selected brands and at concentrations tested in this study have limited mutagenicity in both mouse and human cells in vitro.

Assessment of heritable genetic effects using new genetic tools and sentinels in an era of personalized medicine

The challenge of estimating human health effects from damage to the germ line may be met in the genomic era. Understanding the genetic, as opposed to postconception developmental basis of birth defects is critical to their use in monitoring heritable genetic damage. The causes of common birth defects are analyzed here: mendelian genetic, multigenic, developmental, inherited, or combinational. Only a small fraction of these (noninherited, mendelian genetic) are likely to be informative relative to germ cell mutagenesis, and these won't be discernible against the general background of birth defects. Targeted genetic testing as part of personalized medicine could be integrated into a strategy for assessing germ cell alterations in populations. Thus, "sentinel mutations," as originally proposed by Mulvihill and Ceizel, need not be restricted to X-linked or dominant mutations or conditions visible at birth. Several new sentinels related to personalized medicine are proposed, based on health impact (likelihood of monitoring), frequency, and genetic target suitability (responsiveness to diverse mutational mechanisms). Candidates could include CYP genes (related to metabolism of xenobiotics) important in optimizing drug doses and avoiding adverse reactions. High frequency LDLR mutations (related to familial high cholesterol) predict myocardial infarction in approximately50% of individuals. The more common recessive genetic diseases (cystic fibrosis, phenylketonuria, and others) monitored in newborn screening programs could be informative given parental analysis. New opportunities for genetic analyses need to be coupled with epidemiological studies on environmental exposures. These could focus on adverse outcomes related to tobacco, the mostubiquitous and potent environmental mutagen.

Mutation spectra induced by adozelesin in the supF gene of human XP-A fibroblasts

Adozelesin is a synthetic analog of the antitumor antibiotic, CC-1065, which alkylates N3 of adenine in the minor DNA groove in a sequence-specific manner. Here we tested the mutation spectra induced by adozelesin in the supF gene of human XP-A fibroblasts using a shuttle vector assay. Adozelesin primarily induces mutations via an A --> T transversion and a single base insertion. The A --> T transversion (43/59) was observed at the adenine alkylation site in the 5'-ATTTA* sequence (A* is the site of alkylation). The single base insertion (5/59) was observed at the 3'-side of the covalently modified adenine in the 5'-CTAAA* sequence. The results of this study suggest that the DNA alkylating sequence of adozelesin influences the type of DNA mutation.

New approaches to assessing the effects of mutagenic agents on the integrity of the human genome

Heritable genetic alterations, although individually rare, have a substantial collective health impact. Approximately 20% of these are new mutations of unknown cause. Assessment of the effect of exposures to DNA damaging agents, i.e. mutagenic chemicals and radiations, on the integrity of the human genome and on the occurrence of genetic disease remains a daunting challenge. Recent insights may explain why previous examination of human exposures to ionizing radiation, as in Hiroshima and Nagasaki, failed to reveal heritable genetic effects. New opportunities to assess the heritable genetic damaging effects of environmental mutagens are afforded by: (1) integration of knowledge on the molecular nature of genetic disorders and the molecular effects of mutagens; (2) the development of more practical assays for germline mutagenesis; (3) the likely use of population-based genetic screening in personalized medicine.

Characterization of naphtho[1,2-a]pyrene and naphtho[1,2-e]pyrene DNA adducts in C3H10T1/2 fibroblasts

Polycyclic aromatic hydrocarbons (PAHs) are a class of carcinogenic chemicals that are ubiquitous in the environment. Fjord-region naphthopyrene isomers are structurally similar to the potent fjord-region PAH carcinogen dibenzo[a,l]pyrene and thus have the potential to be potent carcinogens. Naphtho[1,2-a]pyrene (N[1,2-a]P) exhibited similar bacterial mutagenicity and morphological cell transforming activity when compared to benzo[a]pyrene (B[a]P), whereas the structural isomer, naphtho[1,2-e]pyrene (N[1,2-e]P) was inactive is these bioassays. In this study, we examined the formation of DNA adducts in C3H10T1/2Cl8 (C3H10T1/2) mouse embryo fibroblasts exposed to N[1,2-a]P or N[1,2-e]P and their respective dihydrodiols. The DNA adducts were characterized by co-chromatography with reaction products from anti-N[1,2-a]P diol epoxide (DE) or anti-N[1,2-e]PDE and polydeoxyadenosine (dAdo) or oligodeoxyguanosine (dGuo). C3H10T1/2 fibroblasts exposed to N[1,2-a]P or N[1,2-a]P-9,10-diol produced both anti-N[1,2-a]P-DE-dAdo and -dGuo adducts with total DNA adduction levels of 22.2 to 33.3 pmol DNA adducts/mug DNA. C3H10T1/2 fibroblasts exposed to N[1,2-e]P produced 2 major and 1 minor adducts. C3H10T1/2 fibroblasts exposed to N[1,2-e]P-11,12-diol produced 2 major adducts. All of the identified adducts were anti-N[1,2-e]PDE-dGuo and -dAdo adducts. While the total DNA adduct level in N[1,2-e]P-11,12-diol-treated fibroblasts was extremely high, 105.9 pmol DNA adducts/mug DNA, the level in N[1,2-e]P-treated fibroblasts was 1.47 pmol DNA adducts/microg DNA. We conclude that lack of biological activity of N[1,2-e]P may be related to its inability to form sufficient amounts of N[1,2-e]P-11,12-diol, which would then be metabolized to sufficient amounts of anti-N[1,2-e]PDE needed to transform these fibroblasts.

The mutagenicity of urban particulate matter in an enzyme free system is associated with the generation of reactive oxygen species

Urban particulate matter (UPM) contributes to lung cancer incidence. UPM has been shown to be genotoxic to mammalian cells and to induce mutations in the Ames assay. Here, we have studied the induction of mutations generated by direct acting mutagenic components of UPM, using the supF forward mutation assay. Plasmid pSP189 was exposed to UPM in aqueous solution in the presence of sucrose buffer, to reduce strand breaks. The mutation frequency induced by 1 microg/microl UPM was 4.99 mutants per 10(4) colonies. This was reduced to 0.84 and 1.48 mutants per 10(4) colonies by addition of mannitol (1 mM) or EDTA (1 mM), respectively. A large percentage of mutant plasmids contained frameshift mutations (57%), and 31% of mutant plasmids contained multiple mutations. Of the base substitution mutations, 88% were at GC pairs, with twice as many transversions as transitions. The types of mutations induced, the reduction of mutagenicity by the inclusion of the free radical scavenger, mannitol, or the metal chelator, EDTA, and the sequence context of the induced mutations all support the conclusion that the majority of mutations were induced by reactive oxygen species generated by metal ions present in the UPM. Most mutation studies with UPM have focused on organic carcinogens present on UPM. Our results highlight the potential contribution of metal ions to the mutagenicity of UPM.

Co-mutagenic activity of arsenic and benzo[a]pyrene in mouse skin

Exposure to inorganic arsenic in drinking water is linked to skin, lung and bladder cancer in humans. The mechanism of arsenic-induced cancer is not clear, but exposure to arsenic and polycyclic arylhydrocarbons (PAH) is more carcinogenic than exposure to either type of carcinogen alone. Arsenic can also generate reactive oxygen species, suggesting that oxidation of DNA may play a role in carcinogenesis. Oxidization of guanosines in polyG tracts is known to cause frameshift mutations, and such events can be detected in situ using the G11 placental alkaline phosphatase (PLAP) transgenic mouse model, which reports frameshift mutations in a run of 11 G:C basepairs by generating cells containing heat-resistant alkaline phosphatase activity. PAH can also induce frameshift mutations. In the study described here, FVB/N mice carrying the G11 PLAP transgene were crossed to C57Bl/6 mice. Half of the hybrid mice were given drinking water with sodium arsenite (10 mg/L) for 10 weeks. Half of the arsenic treated mice were also exposed to benzo[a]pyrene (BaP) by skin painting (500 nmol/week) for 8 weeks. Another group of mice was exposed to BaP but not arsenic. The effect on frameshift mutation was assessed by staining sections of skin tissue to detect cells with PLAP activity. Arsenic alone had no significant effect. On average, mice given BaP alone had approximately three times more PLAP-positive (PLAP+) cells. By contrast, mice exposed to both arsenic and BaP exhibited 10-fold more PLAP+ cells in the skin, and these cells were often arranged in large clusters, suggesting derivation from stem cells. Whereas combined treatment produced more PLAP+ cells, stable BaP adduct levels and arsenic burdens were not higher in mice exposed to both agents compared to mice exposed to either one agent or the other.

Role of base sequence context in conformational equilibria and nucleotide excision repair of benzo[a]pyrene diol epoxide-adenine adducts

We investigate the influence of base sequence context on the conformations of the 10S (+)- and 10R (-)-trans-anti-[BP]-N(6)-dA adducts through molecular dynamics (MD) simulations with free energy calculations, and relate the structural findings to results of nucleotide excision repair (NER) assays in human cell extracts. In previous studies, these adducts were studied in the CA*A sequence context, and here we report results for the CA*C sequence. Our simulations indicate that the base sequence context affects the syn-anti conformational equilibrium in the 10S (+) adduct by modulating the barrier heights between these states on the energy surface, with a higher barrier in the CA*C case. Our nucleotide excision repair assay finds greater NER susceptibilities in the 10S (+) adduct for the CA*C sequence context. A structural rationale ties together these results. A sequence specific hydrogen bond, accompanied by a significantly increased roll and consequent bending in the 10S (+) adduct, has been found in our simulations for the CA*C sequence, which could account for the enhanced nucleotide excision repair as well as the syn-anti equilibrium difference we observe in this isomer and sequence. Such sequence specific differential repair could contribute to the existence of mutational hotspots and thereby contribute to the complexity of cancer initiation.

Cyclohexene ring and Fjord region twist inversion in stereoisomeric DNA adducts of enantiomeric benzo[c]phenanthrene diol epoxides

The sterically hindered, nonplanar fjord region polycyclic aromatic hydrocarbons (PAHs) have been of great interest because of the exceptionally high mutagenic and tumorigenic activity of certain of their metabolically activated diol epoxides. Benzo[c]phenanthrene (B[c]Ph), a representative fjord region PAH, is metabolically activated to a pair of enantiomers, 1S,2R,3R,4S-3,4-dihydroxy-1,2-epoxy-1,2,3,4-tetrahydrobenzo[c]phenanthrene, (+)-anti-B[c]PhDE, and the corresponding 1R,2S,3S,4R enantiomer, (-)-anti-B[c]PhDE. Both of these can bind covalently to the amino group of purines in DNA via trans addition. In the present work we carry out an extensive computational investigation of the 1R(+) and 1S(-)-trans-anti-B[c]Ph adducts to the base guanine, with the goal of delineating the conformational possibilities for the fjord region and the adjacent cyclohexene-type benzylic ring and their relevance to DNA duplexes. We created 10 369 starting structures for each adduct and minimized the energy using AMBER 5.0. A limited set of conformational families is computed, in which the R isomer structures are near mirror images of the S isomer. The benzylic rings are essentially all half-chair-type. Cyclohexene-type ring inversion as well as fjord region twist inversion are possible for each isomer and are correlated. DNA duplexes modified by fjord region adducts select conformers from the allowed families that optimize stacking interactions, which contributes to the stability of the carcinogen-intercalated DNA duplex structures [Cosman et al. (1993) Biochemistry 32, 12488-12497; Cosman et al. (1995) Biochemistry 34, 1295-1307; Suri et al. (1999) J. Mol. Biol. 292, 289-307; Lin et al. (2001) J. Mol. Biol. 306, 1059-1080]. In turn, this stability could contribute to the resistance to repair by the human nucleotide excision system observed in fjord region adducts [Buterin et al. (2000) Cancer Res. 60, 1849-1856].

Methods for testing compounds for DNA adduct formation

This article, based on a presentation on DNA adduct detection given at a Genetic Toxicology Association workshop, is an overview of methods used for testing compounds for DNA adduct formation. A DNA adduct study may be initiated on a case by case basis when there are conflicting results within the standard battery of genetic toxicology tests or when tumors are detected in the animal bioassay for nongenotoxic compounds. Methods for adduct detection include the 32P-postlabeling assay, the use of radioactive test chemicals, physicochemical methods, and immunoassays. Of these, the 32P-postlabeling assay and the use of radiochemicals are discussed in greater detail, since only these two methods are readily applicable to test a compound for the formation of uncharacterized DNA adducts. The other methods are applicable to those adducts that have been chemically characterized or that contain a fluorophore or electrochemically active groups. Evaluation of mutagenic and carcinogenic risk from DNA adducts would require the understanding of various parameters, including the chemical nature, quantity and stability of adducts, proliferation rates for target cells to fix adducts into mutations, mutagenic and repair efficiencies of adducts, and the extent of modifications in critical genes. Since such data cannot be readily obtainable, the toxicological risk from uncharacterized adducts is difficult to assess.

Mutations induced by reactive nitrogen oxide species in the supF forward mutation assay

Nitric oxide is an important bioregulatory molecule with a range of physiological functions. Nitric oxide can also react with oxygen species to produce a range of reactive nitrogen oxides that can damage DNA and lead to mutations of the DNA base sequence. The mutagenicity of a variety of reactive nitrogen oxide species and related DNA damaging agents in the supF assay are reviewed here, in the context of recent reports that relate to the nature of the DNA lesions responsible for the induced mutations. Mutations induced by nitric oxide in the supF assay are compared to those induced by N(2)O(3), nitrous acid, peroxynitrite and different reactive oxygen species. The effect of replication of the damaged pSP189 plasmid in human cells or Escherichia coli cells is also considered.