High-resolution mapping of DNA polymerase fidelity using nucleotide imbalances and next-generation sequencing

Correction of non-random mutational biases along a linear bacterial chromosome by the mismatch repair endonuclease NucS

The linear chromosome of Streptomyces exhibits a highly compartmentalized structure with a conserved central region flanked by variable arms. As double strand break (DSB) repair mechanisms play a crucial role in shaping the genome plasticity of Streptomyces, we investigated the role of EndoMS/NucS, a recently characterized endonuclease involved in a non-canonical mismatch repair (MMR) mechanism in archaea and actinobacteria, that singularly corrects mismatches by creating a DSB. We showed that Streptomyces mutants lacking NucS display a marked colonial phenotype and a drastic increase in spontaneous mutation rate. In vitro biochemical assays revealed that NucS cooperates with the replication clamp to efficiently cleave G/T, G/G and T/T mismatched DNA by producing DSBs. These findings are consistent with the transition-shifted mutational spectrum observed in the mutant strains and reveal that NucS-dependent MMR specific task is to eliminate G/T mismatches generated by the DNA polymerase during replication. Interestingly, our data unveil a crescent-shaped distribution of the transition frequency from the replication origin towards the chromosomal ends, shedding light on a possible link between NucS-mediated DSBs and Streptomyces genome evolution.

Expedient production of site specifically nucleobase-labelled or hypermodified RNA with engineered thermophilic DNA polymerases

Innovative approaches to controlled nucleobase-modified RNA synthesis are urgently needed to support RNA biology exploration and to synthesize potential RNA therapeutics. Here we present a strategy for enzymatic construction of nucleobase-modified RNA based on primer-dependent engineered thermophilic DNA polymerases - SFM4-3 and TGK. We demonstrate introduction of one or several different base-modified nucleotides in one strand including hypermodified RNA containing all four modified nucleotides bearing four different substituents, as well as strategy for primer segment removal. We also show facile site-specific or segmented introduction of fluorophores or other functional groups at defined positions in variety of RNA molecules, including structured or long mRNA. Intriguing translation efficacy of single-site modified mRNAs underscores the necessity to study isolated modifications placed at designer positions to disentangle their biological effects and enable development of improved mRNA therapeutics. Our toolbox paves the way for more precise dissecting RNA structures and functions, as well as for construction of diverse types of base-functionalized RNA for therapeutic applications and diagnostics.

Neq2X7: a multi-purpose and open-source fusion DNA polymerase for advanced DNA engineering and diagnostics PCR

Thermostable DNA polymerases, such as Taq isolated from the thermophilic bacterium Thermus aquaticus, enable one-pot exponential DNA amplification known as polymerase chain reaction (PCR). However, properties other than thermostability - such as fidelity, processivity, and compatibility with modified nucleotides - are important in contemporary molecular biology applications. Here, we describe the engineering and characterization of a fusion between a DNA polymerase identified in the marine archaea Nanoarchaeum equitans and a DNA binding domain from the thermophile Sulfolobus solfataricus. The fusion creates a highly active enzyme, Neq2X7, capable of amplifying long and GC-rich DNA, unaffected by replacing dTTP with dUTP in PCR, and tolerant to various known PCR inhibitors. This makes it an attractive DNA polymerase for use, e.g., with uracil excision (USER) DNA assembly and for contamination-free diagnostics. Using a magnification via nucleotide imbalance fidelity assay, Neq2X7 was estimated to have an error rate lower than 2 ∙ 10-5 bp-1 and an approximately 100x lower fidelity than the parental variant Neq2X, indicating a trade-off between fidelity and processivity - an observation that may be of importance for similarly engineered DNA polymerases. Neq2X7 is easy to produce for routine application in any molecular biology laboratory, and the expression plasmid is made freely available.

DNA Polymerases for Whole Genome Amplification: Considerations and Future Directions

In the same way that specialized DNA polymerases (DNAPs) replicate cellular and viral genomes, only a handful of dedicated proteins from various natural origins as well as engineered versions are appropriate for competent exponential amplification of whole genomes and metagenomes (WGA). Different applications have led to the development of diverse protocols, based on various DNAPs. Isothermal WGA is currently widely used due to the high performance of Φ29 DNA polymerase, but PCR-based methods are also available and can provide competent amplification of certain samples. Replication fidelity and processivity must be considered when selecting a suitable enzyme for WGA. However, other properties, such as thermostability, capacity to couple replication, and double helix unwinding, or the ability to maintain DNA replication opposite to damaged bases, are also very relevant for some applications. In this review, we provide an overview of the different properties of DNAPs widely used in WGA and discuss their limitations and future research directions.

CReSIL: accurate identification of extrachromosomal circular DNA from long-read sequences

Extrachromosomal circular DNA (eccDNA) of chromosomal origin is found in many eukaryotic species and cell types, including cancer, where eccDNAs with oncogenes drive tumorigenesis. Most studies of eccDNA employ short-read sequencing for their identification. However, short-read sequencing cannot resolve the complexity of genomic repeats, which can lead to missing eccDNA products. Long-read sequencing technologies provide an alternative to constructing complete eccDNA maps. We present a software suite, Construction-based Rolling-circle-amplification for eccDNA Sequence Identification and Location (CReSIL), to identify and characterize eccDNA from long-read sequences. CReSIL's performance in identifying eccDNA, with a minimum F1 score of 0.98, is superior to the other bioinformatic tools based on simulated data. CReSIL provides many useful features for genomic annotation, which can be used to infer eccDNA function and Circos visualization for eccDNA architecture investigation. We demonstrated CReSIL's capability in several long-read sequencing datasets, including datasets enriched for eccDNA and whole genome datasets from cells containing large eccDNA products. In conclusion, the CReSIL suite software is a versatile tool for investigating complex and simple eccDNA in eukaryotic cells.

Prediction and validation of hematopoietic stem and progenitor cell off-target editing in transplanted rhesus macaques

The programmable nuclease technology CRISPR-Cas9 has revolutionized gene editing in the last decade. Due to the risk of off-target editing, accurate and sensitive methods for off-target characterization are crucial prior to applying CRISPR-Cas9 therapeutically. Here, we utilized a rhesus macaque model to compare the predictive values of CIRCLE-seq, an in vitro off-target prediction method, with in silico prediction (ISP) based solely on genomic sequence comparisons. We use AmpliSeq HD error-corrected sequencing to validate off-target sites predicted by CIRCLE-seq and ISP for a CD33 guide RNA (gRNA) with thousands of off-target sites predicted by ISP and CIRCLE-seq. We found poor correlation between the sites predicted by the two methods. When almost 500 sites predicted by each method were analyzed by error-corrected sequencing of hematopoietic cells following transplantation, 19 off-target sites revealed insertion or deletion mutations. Of these sites, 8 were predicted by both methods, 8 by CIRCLE-seq only, and 3 by ISP only. The levels of cells with these off-target edits exhibited no expansion or abnormal behavior in vivo in animals followed for up to 2 years. In addition, we utilized an unbiased method termed CAST-seq to search for translocations between the on-target site and off-target sites present in animals following transplantation, detecting one specific translocation that persisted in blood cells for at least 1 year following transplantation. In conclusion, neither CIRCLE-seq or ISP predicted all sites, and a combination of careful gRNA design, followed by screening for predicted off-target sites in target cells by multiple methods, may be required for optimizing safety of clinical development.

Accurate and scalable variant calling from single cell DNA sequencing data with ProSolo

Accurate single cell mutational profiles can reveal genomic cell-to-cell heterogeneity. However, sequencing libraries suitable for genotyping require whole genome amplification, which introduces allelic bias and copy errors. The resulting data violates assumptions of variant callers developed for bulk sequencing. Thus, only dedicated models accounting for amplification bias and errors can provide accurate calls. We present ProSolo for calling single nucleotide variants from multiple displacement amplified (MDA) single cell DNA sequencing data. ProSolo probabilistically models a single cell jointly with a bulk sequencing sample and integrates all relevant MDA biases in a site-specific and scalable-because computationally efficient-manner. This achieves a higher accuracy in calling and genotyping single nucleotide variants in single cells in comparison to state-of-the-art tools and supports imputation of insufficiently covered genotypes, when downstream tools cannot handle missing data. Moreover, ProSolo implements the first approach to control the false discovery rate reliably and flexibly. ProSolo is implemented in an extendable framework, with code and usage at: https://github.com/prosolo/prosolo.

Recording Temporal Signals with Minutes Resolution Using Enzymatic DNA Synthesis

Employing DNA as a high-density data storage medium has paved the way for next-generation digital storage and biosensing technologies. However, the multipart architecture of current DNA-based recording techniques renders them inherently slow and incapable of recording fluctuating signals with subhour frequencies. To address this limitation, we developed a simplified system employing a single enzyme, terminal deoxynucleotidyl transferase (TdT), to transduce environmental signals into DNA. TdT adds nucleotides to the 3'-ends of single-stranded DNA (ssDNA) in a template-independent manner, selecting bases according to inherent preferences and environmental conditions. By characterizing TdT nucleotide selectivity under different conditions, we show that TdT can encode various physiologically relevant signals such as Co2+, Ca2+, and Zn2+ concentrations and temperature changes in vitro. Further, by considering the average rate of nucleotide incorporation, we show that the resulting ssDNA functions as a molecular ticker tape. With this method we accurately encode a temporal record of fluctuations in Co2+ concentration to within 1 min over a 60 min period. Finally, we engineer TdT to allosterically turn off in the presence of a physiologically relevant concentration of calcium. We use this engineered TdT in concert with a reference TdT to develop a two-polymerase system capable of recording a single-step change in the Ca2+ signal to within 1 min over a 60 min period. This work expands the repertoire of DNA-based recording techniques by developing a novel DNA synthesis-based system that can record temporal environmental signals into DNA with a resolution of minutes.

Kinetic assays of DNA polymerase fidelity: A theoretical perspective beyond Michaelis-Menten kinetics

The high fidelity of DNA polymerase (DNAP) is critical for the faithful replication of DNA. There are several quantitative approaches to measure DNAP fidelity. Directly counting the error frequency in the replication products gives the true fidelity but it turns out very hard to implement in practice. Two biochemical kinetic approaches, the steady-state assay and the transient-state assay, were then suggested and widely adopted. In these assays, the error frequency is indirectly estimated by using kinetic theories combined with the measured apparent kinetic rates. However, whether it is equivalent to the true fidelity has never been clarified theoretically, and in particular there are different strategies using these assays to quantify the proofreading efficiency of DNAP but often lead to inconsistent results. In this paper, we make a comprehensive examination on the theoretical foundation of the two kinetic assays, based on the theory of DNAP fidelity recently proposed by us. Our studies show that while the conventional kinetic assays are generally valid to quantify the discrimination efficiency of DNAP, they are valid to quantify the proofreading efficiency of DNAP only when the kinetic parameters satisfy some constraints which will be given explicitly in this paper. These results may inspire more carefully-designed experiments to quantify DNAP fidelity.

Dual Deep Sequencing Improves the Accuracy of Low-Frequency Somatic Mutation Detection in Cancer Gene Panel Testing

Cancer gene panel testing requires accurate detection of somatic mosaic mutations, as the test sample consists of a mixture of cancer cells and normal cells; each minor clone in the tumor also has different somatic mutations. Several studies have shown that the different types of software used for variant calling for next generation sequencing (NGS) can detect low-frequency somatic mutations. However, the accuracy of these somatic variant callers is unknown. We performed cancer gene panel testing in duplicate experiments using three different high-fidelity DNA polymerases in pre-capture amplification steps and analyzed by three different variant callers, Strelka2, Mutect2, and LoFreq. We selected six somatic variants that were detected in both experiments with more than two polymerases and by at least one variant caller. Among them, five single nucleotide variants were verified by CEL nuclease-mediated heteroduplex incision with polyacrylamide gel electrophoresis and silver staining (CHIPS) and Sanger sequencing. In silico analysis indicated that the FBXW7 and MAP3K1 missense mutations cause damage at the protein level. Comparing three somatic variant callers, we found that Strelka2 detected more variants than Mutect2 and LoFreq. We conclude that dual sequencing with Strelka2 analysis is useful for detection of accurate somatic mutations in cancer gene panel testing.

Development of a next generation DNA sequencing-based multi detection assay for detecting and identifying Leishmania parasites, blood sources, plant meals and intestinal microbiome in phlebotomine sand flies

Leishmaniasis is a disease caused by Leishmania parasites transmitted by phlebotomine sand flies (Diptera: Psychodidae). Human infections with different Leishmania species cause characteristic clinical manifestations; cutaneous or visceral leishmaniasis. Here we describe the development and application of a Miseq Next GenerationSequencing (NGS)-based Multi Detection Assay (MDA) designed to characterize metagenomics parameters pertinent to the sand fly vectors which may affect their vectorial capacity for Leishmania. For this purpose, we developed a MDA by which, DNA fragments were amplified through polymerase chain reactions (PCR) and then sequenced by MiSeq/NGS. PCR amplification was achieved using some published and some new primers designed specifically for identifying Leishmania spp. (ITS1), sand fly spp. (cytochrome oxidase I), vertebrate blood (Cytochrome b), plant DNA ribulose-1,5-bisphosphate carboxylase large subunit gene (rbcL), and prokaryotic micobiome (16 s rRNA). This MDA/NGS analysis was performed on two species of wild-caught sand flies that transmit different Leishmania spp. in two ecologically distinct, but geographically neighboring locations. The results were analyzed to identify, quantitate and correlate the measured parameters in order to assess their putative importance in the transmission dynamics of leishmaniasis.

Calling Variants in the Clinic: Informed Variant Calling Decisions Based on Biological, Clinical, and Laboratory Variables

Deep sequencing genomic analysis is becoming increasingly common in clinical research and practice, enabling accurate identification of diagnostic, prognostic, and predictive determinants. Variant calling, distinguishing between true mutations and experimental errors, is a central task of genomic analysis and often requires sophisticated statistical, computational, and/or heuristic techniques. Although variant callers seek to overcome noise inherent in biological experiments, variant calling can be significantly affected by outside factors including those used to prepare, store, and analyze samples. The goal of this review is to discuss known experimental features, such as sample preparation, library preparation, and sequencing, alongside diverse biological and clinical variables, and evaluate their effect on variant caller selection and optimization.

The live cell DNA stain SiR-Hoechst induces DNA damage responses and impairs cell cycle progression

SiR-Hoechst (SiR-DNA) is a far-red fluorescent DNA probe being used widely for time-lapse imaging of living cells that is reported to be minimally toxic at concentrations as high as 10-25 µM. However, measuring nuclear import of Cyclin B1, inhibition of mitotic entry, and the induction of γH2AX foci in cultured human cells reveals that SiR-Hoechst induces DNA damage responses and G2 arrest at concentrations well below 1 µM. SiR-Hoechst is useful for live cell imaging, but it should be used with caution and at the lowest practicable concentration.