Mutation assay using single-molecule real-time (SMRT(TM)) sequencing technology

DNA mismatch and damage patterns revealed by single-molecule sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other diseases1,2. Most mutations begin as nucleotide mismatches or damage in one of the two strands of the DNA before becoming double-strand mutations if unrepaired or misrepaired3,4. However, current DNA-sequencing technologies cannot accurately resolve these initial single-strand events. Here we develop a single-molecule, long-read sequencing method (Hairpin Duplex Enhanced Fidelity sequencing (HiDEF-seq)) that achieves single-molecule fidelity for base substitutions when present in either one or both DNA strands. HiDEF-seq also detects cytosine deamination-a common type of DNA damage-with single-molecule fidelity. We profiled 134 samples from diverse tissues, including from individuals with cancer predisposition syndromes, and derive from them single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumours deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples that are deficient in only polymerase proofreading. We also define a single-strand damage signature for APOBEC3A. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. As double-strand DNA mutations are only the end point of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable studies of how mutations arise in a variety of contexts, especially in cancer and ageing.

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

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

Unbiased whole genome detection of ultrarare off-target mutations in genome-edited cell populations by HiFi sequencing

DNA base editors (BEs) composed of a nuclease-deficient Cas9 fused to a DNA-modifying enzyme can achieve on-target mutagenesis without creating double-strand DNA breaks (DSBs). As a result, BEs generate far less DNA damage than traditional nuclease-proficient Cas9 systems, which do rely on the creation of DSBs to achieve on-target mutagenesis. The inability of BEs to create DSBs makes the detection of their undesired off-target effects very difficult. PacBio HiFi sequencing can efficiently detect ultrarare mutations resulting from chemical mutagenesis in whole genomes with a sensitivity ~1 × 10-8 mutations per base pair. In this proof-of-principle study, we evaluated whether this technique could also detect the on- and off-target mutations generated by a cytosine-to-thymine (C>T) BE targeting the LacZ gene in Escherichia coli (E. coli). HiFi sequencing detected on-target mutant allele fractions ranging from ~7% to ~63%, depending on the single-guide RNA (sgRNA) used, while no on-target mutations were detected in controls lacking the BE. The presence of the BE resulted in a ~3-fold increase in mutation frequencies compared to controls lacking the BE, irrespective of the sgRNA used. These increases were mostly composed of C:G>T:A substitutions distributed throughout the genome. Our results demonstrate that HiFi sequencing can efficiently identify on- and off-target mutations in cell populations that have undergone genome editing.

Single-strand mismatch and damage patterns revealed by single-molecule DNA sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other genetic diseases1-4. Almost all of these mosaic mutations begin as nucleotide mismatches or damage in only one of the two strands of the DNA prior to becoming double-strand mutations if unrepaired or misrepaired5. However, current DNA sequencing technologies cannot resolve these initial single-strand events. Here, we developed a single-molecule, long-read sequencing method that achieves single-molecule fidelity for single-base substitutions when present in either one or both strands of the DNA. It also detects single-strand cytosine deamination events, a common type of DNA damage. We profiled 110 samples from diverse tissues, including from individuals with cancer-predisposition syndromes, and define the first single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumors deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples deficient in only polymerase proofreading. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. Since the double-strand DNA mutations interrogated by prior studies are only the endpoint of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable new studies of how mutations arise in a variety of contexts, especially in cancer and aging.

Pharmacogenomics in drug-induced cardiotoxicity: Current status and the future

Drug-induced cardiotoxicity (DICT) is an important concern of drug safety in both drug development and clinical application. The clinical manifestations of DICT include cardiomyopathy, arrhythmia, myocardial ischemia, heart failure, and a series of cardiac structural and functional changes. The occurrence of DICT has negative impacts on the life quality of the patients, brings additional social and economic burden. It is important to identify the potential factors and explore the mechanisms of DICT. Traditional cardiovascular risk factors can only partially explain the risk of DICT. Pharmacogenomic studies show accumulated evidence of genetics in DICT and suggest the potential to guide precision therapy to reduce risk of cardiotoxicity. The comprehensive application of technologies such as third-generation sequencing, human induced pluripotent stem (iPS) cells and genome editing has promoted the in-depth understanding of the functional role of susceptible genes in DICT. This paper reviewed drugs that cause DICT, the clinical manifestations and laboratory tests, as well as the related content of genetic variations associated with the risk of DICT, and further discussed the implication of new technologies in pharmacogenomics of DICT.

Evaluation of the mutagenic effects of Molnupiravir and N4-hydroxycytidine in bacterial and mammalian cells by HiFi sequencing

Molnupiravir (MOV) is used to treat COVID-19. In cells, MOV is converted to the ribonucleoside analog N4-hydroxycytidine (NHC) and incorporated into the SARS-CoV-2 RNA genome during its replication, resulting in RNA mutations. The widespread accumulation of such mutations inhibits SARS-CoV-2 propagation. Although safety assessments by many regulatory agencies across the world have concluded that the genotoxic risks associated with the clinical use of MOV are low, concerns remain that it could induce DNA mutations in patients, particularly because numerous in vitro studies have shown that NHC is a DNA mutagen. In this study, we used HiFi sequencing, a technique that can detect ultralow-frequency substitution mutations in whole genomes, to evaluate the mutagenic effects of MOV in E. coli and of MOV and NHC in mouse lymphoma L5178Y cells and human lymphoblastoid TK6 cells. In all models, exposure to these compounds increased genome-wide mutation frequencies in a dose-dependent manner, and these increases were mainly composed of A:T → G:C transitions. The NHC exposure concentrations used for mammalian cells were comparable to those observed in the plasma of humans who received clinical doses of MOV.

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

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

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

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

Hawk-Seq™ differentiates between various mutations in Salmonella typhimurium TA100 strain caused by exposure to Ames test-positive mutagens

A precise understanding of differences in genomic mutations according to the mutagenic mechanisms detected in mutagenicity data is required to evaluate the carcinogenicity of environmental mutagens. Recently, we developed a highly accurate genome sequencing method, ‘Hawk-Seq™’, that enables the detection of mutagen-induced genome-wide mutations. However, its applicability to detect various mutagens and identify differences in mutational profiles is not well understood. Thus, we evaluated DNA samples from Salmonella typhimurium TA100 exposed to 11 mutagens, including alkylating agents, aldehydes, an aromatic nitro compound, epoxides, aromatic amines and polycyclic aromatic hydrocarbons (PAHs). We extensively analysed mutagen-induced mutational profiles and studied their association with the mechanisms of mutagens. Hawk-Seq™ sensitively detected mutations induced by all 11 mutagens, including one that increased the number of revertants by approximately 2-fold in the Ames test. Although the sensitivity for less water-soluble mutagens was relatively low, we increased the sensitivity to obtain high-resolution spectra by modifying the exposure protocol. Moreover, two epoxides indicated similar 6- or 96-dimensional mutational patterns; likewise, three S_N1-type alkylating agents indicated similar mutational patterns, suggesting that the mutational patterns are compound category specific. Meanwhile, an S_N2 type alkylating agent exhibited unique mutational patterns compared to those of the S_N1 type alkylating agents. Although the mutational patterns induced by aldehydes, the aromatic nitro compound, aromatic amines and PAHs did not differ substantially from each other, the maximum total base substitution frequencies (MTSFs) were similar among mutagens in the same structural groups. Furthermore, the MTSF was found to be associated with the carcinogenic potency of some direct-acting mutagens. These results indicate that our method can generate high-resolution mutational profiles to identify characteristic features of each mutagen. The detailed mutational data obtained by Hawk-Seq™ can provide useful information regarding mutagenic mechanisms and help identify its association with the carcinogenicity of mutagens without requiring carcinogenicity data.

My career development with Ames test: A personal recollection

I first became acquainted with the Ames test at the very beginning of my career in 1978, when my task at the National Institute of Health Sciences (Tokyo) was to screen for mutagenicity of food additives used in Japan, using the Ames test. I also used this test to research the metabolic activation mechanisms of chemical carcinogens, in particular, the analgesic drug, phenacetin. This chemical was not mutagenic in Salmonella typhimurium TA100 with standard 9000 × g supernatant of liver homogenates (S9) from rat but was mutagenic with hamster S9. It was revealed that hamster S9 had much higher deacetylation activities than rat S9, which accounts for the species difference. Then, my work was focused on molecular biology. We cloned the genes encoding nitroreductase and acetyltransferase in Salmonella typhimurium TA1538. Plasmids carrying these genes made strain TA98 more sensitive to mutagenic nitroarenes and aromatic amines. Because of their high sensitivity, the resulting strains such as YG1021 and YG1024 are widely used to monitor mutagenic nitroarenes and aromatic amines in complex mixtures. Later, we disrupted the genes encoding DNA polymerases in TA1538 and classified chemical mutagens into four classes depending on their use of different DNA polymerases. I was also involved in the generation of gpt delta transgenic rodent gene mutation assays, which examine the results of the Ames test in vivo. I have unintentionally developed my career under the influence of Dr. Ames and I would like to acknowledge his remarkable achievements in the field of environmental mutagenesis and carcinogenesis.

A decision tree-based integrated testing strategy for tailor-made carcinogenicity evaluation of test substances using genotoxicity test results and chemical spaces

Genotoxicity evaluation has been widely used to estimate the carcinogenicity of test substances during safety evaluation. However, the latest strategies using genotoxicity tests give more weight to sensitivity; therefore, their accuracy has been very low. For precise carcinogenicity evaluation, we attempted to establish an integrated testing strategy for the tailor-made carcinogenicity evaluation of test materials, considering the relationships among genotoxicity test results (Ames, in vitro mammalian genotoxicity and in vivo micronucleus), carcinogenicity test results and chemical properties (molecular weight, logKow and 179 organic functional groups). By analyzing the toxicological information and chemical properties of 230 chemicals, including 184 carcinogens in the Carcinogenicity Genotoxicity eXperience database, a decision tree for carcinogenicity evaluation was optimised statistically. A decision forest model was generated using a machine-learning method-random forest-which comprises thousands of decision trees. As a result, balanced accuracies in cross-validation of the optimised decision tree and decision forest model, considering chemical space (71.5% and 75.5%, respectively), were higher than balanced accuracy of an example regulatory decision tree (54.1%). Moreover, the statistical optimisation of tree-based models revealed significant organic functional groups that would cause false prediction in standard genotoxicity tests and non-genotoxic carcinogenicity (e.g., organic amide and thioamide, saturated heterocyclic fragment and aryl halide). In vitro genotoxicity tests were the most important parameters in all models, even when in silico parameters were integrated. Although external validation is required, the findings of the integrated testing strategies established herein will contribute to precise carcinogenicity evaluation and to determine new mechanistic hypotheses of carcinogenicity.

Development of a supF-based mutation-detection system in the extreme thermophile Thermus thermophilus HB27

Thermus thermophilus (T. thermophilus) HB27 is an extreme thermophile that grows optimally at 65-72 °C. Heat-induced DNA lesions are expected to occur at a higher frequency in the genome of T. thermophilus than in those of mesophiles; however, the mechanisms underlying the maintenance of genome integrity at high temperatures remain poorly understood. The study of mutation spectra has become a powerful approach to understanding the molecular mechanisms responsible for DNA repair and mutagenesis in mesophilic species. Therefore, we developed a supF-based system to detect a broad spectrum of mutations in T. thermophilus. This system was validated by measuring spontaneous mutations in the wild type and a udgA, B double mutant deficient in uracil-DNA glycosylase (UDG) activity. We found that the mutation frequency of the udgA, B strain was 4.7-fold higher than that of the wild type and G:C→A:T transitions dominated, which was the most reasonable for the mutator phenotype associated with the loss of UDG function in T. thermophilus. These results show that this system allowed for the rapid analysis of mutations in T. thermophilus, and may be useful for studying the molecular mechanisms responsible for DNA repair and mutagenesis in this extreme thermophile.

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.

Genome resequencing analysis of Salmonella typhimurium LT-2 strains TA98 and TA100 for the establishment of a next-generation sequencing-based mutagenicity assay

Next-generation sequencing (NGS) is a potentially useful technology to achieve a more precise evaluation of chemical mutagenicity. To establish NGS-based mutagenicity assays, which enable the direct detection of chemically induced mutations in a whole genome manner, the selection of appropriate biological resources and their precise genome sequences are essential. Here, we performed genome re-sequencing analyses of Salmonella typhimurium LT-2 strains TA98 and TA100, which have been frequently used in mutagenicity assays. We identified several strain-specific mutations including those that were relevant to their known phenotypes (his, ΔuvrB and rfa). The details of rfa mutations were first clarified in this study, which was a frameshift variant in rfaF and a missense variant in rfaC in TA98 and TA100, respectively. The uvrB deletion in TA98 was larger than that in TA100, which suggested differences in defects of lipopolysaccharide synthesis between these strains. The re-sequenced genome data of TA98 and TA100 will help us establish NGS-based bacterial mutagenicity assays and understand the biological events seen in them. Copyright © 2017 John Wiley & Sons, Ltd.

Mutation assay using single-molecule real-time (SMRT(TM)) sequencing technology

Introduction

We present here a simple, phenotype-independent mutation assay using a PacBio RSII DNA sequencer employing single-molecule real-time (SMRT) sequencing technology. Salmonella typhimurium YG7108 was treated with the alkylating agent N-ethyl-N-nitrosourea (ENU) and grown though several generations to fix the induced mutations, the DNA was extracted and the mutations were analyzed by using the SMRT DNA sequencer.

Results

The ENU-induced base-substitution frequency was 15.4 per Megabase pair, which is highly consistent with our previous results based on colony isolation and next-generation sequencing. The induced mutation spectrum (95% G:C → A:T, 5% A:T → G:C) is also consistent with the known ENU signature. The base-substitution frequency of the control was calculated to be less than 0.12 per Megabase pair. A current limitation of the approach is the high frequency of artifactual insertion and deletion mutations it detects.

Conclusions

Ultra-low frequency base-substitution mutations can be detected directly by using the SMRT DNA sequencer, and this technology provides a phenotype-independent mutation assay.