Single-strand specific nuclease enhances accuracy of error-corrected sequencing and improves rare mutation-detection sensitivity

Genome-wide direct quantification of in vivo mutagenesis using high-accuracy paired-end and complementary consensus sequencing

Error-corrected next-generation sequencing (ecNGS) is an emerging technology for accurately measuring somatic mutations. Here, we report paired-end and complementary consensus sequencing (PECC-Seq), a high-accuracy ecNGS approach for genome-wide somatic mutation detection. We characterize a novel 2-aminoimidazolone lesion besides 7,8-dihydro-8-oxoguanine and the resulting end-repair artifacts originating from NGS library preparation that obscure the sequencing accuracy of NGS. We modify library preparation protocol for the enzymatic removal of end-repair artifacts and improve the accuracy of our previously developed duplex consensus sequencing method. Optimized PECC-Seq shows an error rate of <5 × 10-8 with consensus bases compressed from approximately 25 Gb of raw sequencing data, enabling the accurate detection of low-abundance somatic mutations. We apply PECC-Seq to the quantification of in vivo mutagenesis. Compared with the classic gpt gene mutation assay using gpt delta transgenic mice, PECC-Seq exhibits high sensitivity in quantitatively measuring dose-dependent mutagenesis induced by Aristolochic acid I (AAI). Moreover, PECC-Seq specifically characterizes the distinct genome-wide mutational signatures of AAI, Benzo[a]pyrene, N-Nitroso-N-ethylurea and N-nitrosodiethylamine and reveals the mutational signature of Quinoline in common mouse models. Overall, our findings demonstrate that high-accuracy PECC-Seq is a promising tool for genome-wide somatic mutagenesis quantification and for in vivo mutagenicity testing.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Variations in the enzymatic activity of S1-type nucleases results from differences in their active site structures

BACKGROUND

S1-like nucleases are widespread enzymes commonly used in biotechnology and molecular biology. Although it is commonly believed that they are mainly Zn2+-dependent acidic enzymes, we have found that numerous members of this family deviate from this rule. Therefore, in this work, we decided to check how broad is the range of non‑zinc-dependent S1-like nucleases and what is the molecular basis of their activities.

METHODS

S1-like nucleases chosen for analysis were achieved through heterologous expression in appropriate eukaryotic hosts. To characterize nucleases' active-site properties, point mutations were introduced in selected positions. The enzymatic activities of wild-type and mutant nucleases were tested by in-gel nuclease activity assay.

RESULTS

We discovered that S1-like nucleases encoded by non-vascular plants and single-celled protozoa, like their higher plant homologues, exhibit a large variety of catalytic properties. We have shown that these individual properties are determined by specific non-conserved active site residues.

CONCLUSIONS

Our findings demonstrate that mutations that occur during evolution can significantly alter the catalytic properties of S1-like nucleases. As a result, different ions can compete for particular S1-type nucleases' active sites. This phenomenon undermines the existing classification of S1-like nucleases.

GENERAL SIGNIFICANCE

Our findings have numerous implications for applications and understanding the S1-like nucleases' biological functions. For example, new biotechnological applications should take into account their unexpected catalytic properties. Moreover, these results demonstrate that the trinuclear zinc-based model commonly used to characterize the catalytic activities of S1-like nucleases is insufficient to explain the actions of non‑zinc-dependent members of this family.

Preparation of Duplex Sequencing Libraries for Archival Paraffin-Embedded Tissue Samples Using Single-Strand-Specific Nuclease P1

DNA from formalin-fixed paraffin-embedded (FFPE) tissues, which are frequently utilized in cancer research, is significantly affected by chemical degradation. It was suggested that approaches that are based on duplex sequencing can significantly improve the accuracy of mutation detection in FFPE-derived DNA. However, the original duplex sequencing method cannot be utilized for the analysis of formalin-fixed paraffin-embedded (FFPE) tissues, as FFPE DNA contains an excessive number of damaged bases, and these lesions are converted to false double-strand nucleotide substitutions during polymerase-driven DNA end repair process. To resolve this drawback, we replaced DNA polymerase by a single strand-specific nuclease P1. Nuclease P1 was shown to efficiently remove RNA from DNA preparations, to fragment the FFPE-derived DNA and to remove 5'/3'-overhangs. To assess the performance of duplex sequencing-based methods in FFPE-derived DNA, we constructed the Bottleneck Sequencing System (BotSeqS) libraries from five colorectal carcinomas (CRCs) using either DNA polymerase or nuclease P1. As expected, the number of identified mutations was approximately an order of magnitude higher in libraries prepared with DNA polymerase vs. nuclease P1 (626 ± 167/Mb vs. 75 ± 37/Mb, paired t-test p-value 0.003). Furthermore, the use of nuclease P1 but not polymerase-driven DNA end repair allowed a reliable discrimination between CRC tumors with and without hypermutator phenotypes. The utility of newly developed modification was validated in the collection of 17 CRCs and 5 adjacent normal tissues. Nuclease P1 can be recommended for the use in duplex sequencing library preparation from FFPE-derived DNA.