Identification and Characterization of Thermostable Y-Family DNA Polymerases η, ι, κ and Rev1 From a Lower Eukaryote, Thermomyces lanuginosus.

Y-family DNA polymerases (pols) consist of six phylogenetically separate subfamilies; two UmuC (polV) branches, DinB (pol IV, Dpo4, polκ), Rad30A/POLH (polη), and Rad30B/POLI (polι) and Rev1. Of these subfamilies, DinB orthologs are found in all three domains of life; eubacteria, archaea, and eukarya. UmuC orthologs are identified only in bacteria, whilst Rev1 and Rad30A/B orthologs are only detected in eukaryotes. Within eukaryotes, a wide array of evolutionary diversity exists. Humans possess all four Y-family pols (pols η, ι, κ, and Rev1), Schizosaccharomyces pombe has three Y-family pols (pols η, κ, and Rev1), and Saccharomyces cerevisiae only has polη and Rev1. Here, we report the cloning, expression, and biochemical characterization of the four Y-family pols from the lower eukaryotic thermophilic fungi, Thermomyces lanuginosus. Apart from the expected increased thermostability of the T. lanuginosus Y-family pols, their major biochemical properties are very similar to properties of their human counterparts. In particular, both Rad30B homologs (T. lanuginosus and human polÉ©) exhibit remarkably low fidelity during DNA synthesis that is template sequence dependent. It was previously hypothesized that higher organisms had acquired this property during eukaryotic evolution, but these observations imply that polι originated earlier than previously known, suggesting a critical cellular function in both lower and higher eukaryotes.

Rapid identification of methylase specificity (RIMS-seq) jointly identifies methylated motifs and generates shotgun sequencing of bacterial genomes.

DNA methylation is widespread amongst eukaryotes and prokaryotes to modulate gene expression and confer viral resistance. 5-Methylcytosine (m5C) methylation has been described in genomes of a large fraction of bacterial species as part of restriction-modification systems, each composed of a methyltransferase and cognate restriction enzyme. Methylases are site-specific and target sequences vary across organisms. High-throughput methods, such as bisulfite-sequencing can identify m5C at base resolution but require specialized library preparations and single molecule, real-time (SMRT) sequencing usually misses m5C. Here, we present a new method called RIMS-seq (rapid identification of methylase specificity) to simultaneously sequence bacterial genomes and determine m5C methylase specificities using a simple experimental protocol that closely resembles the DNA-seq protocol for Illumina. Importantly, the resulting sequencing quality is identical to DNA-seq, enabling RIMS-seq to substitute standard sequencing of bacterial genomes. Applied to bacteria and synthetic mixed communities, RIMS-seq reveals new methylase specificities, supporting routine study of m5C methylation while sequencing new genomes.

Enzymatic methyl sequencing detects DNA methylation at single-base resolution from picograms of DNA.

Bisulfite sequencing detects 5mC and 5hmC at single-base resolution. However, bisulfite treatment damages DNA, which results in fragmentation, DNA loss, and biased sequencing data. To overcome these problems, enzymatic methyl-seq (EM-seq) was developed. This method detects 5mC and 5hmC using two sets of enzymatic reactions. In the first reaction, TET2 and T4-BGT convert 5mC and 5hmC into products that cannot be deaminated by APOBEC3A. In the second reaction, APOBEC3A deaminates unmodified cytosines by converting them to uracils. Therefore, these three enzymes enable the identification of 5mC and 5hmC. EM-seq libraries were compared with bisulfite-converted DNA, and each library type was ligated to Illumina adaptors before conversion. Libraries were made using NA12878 genomic DNA, cell-free DNA, and FFPE DNA over a range of DNA inputs. The 5mC and 5hmC detected in EM-seq libraries were similar to those of bisulfite libraries. However, libraries made using EM-seq outperformed bisulfite-converted libraries in all specific measures examined (coverage, duplication, sensitivity, etc.). EM-seq libraries displayed even GC distribution, better correlations across DNA inputs, increased numbers of CpGs within genomic features, and accuracy of cytosine methylation calls. EM-seq was effective using as little as 100 pg of DNA, and these libraries maintained the described advantages over bisulfite sequencing. EM-seq library construction, using challenging samples and lower DNA inputs, opens new avenues for research and clinical applications.

Nondestructive enzymatic deamination enables single-molecule long-read amplicon sequencing for the determination of 5-methylcytosine and 5-hydroxymethylcytosine at single-base resolution.

The predominant methodology for DNA methylation analysis relies on the chemical deamination by sodium bisulfite of unmodified cytosine to uracil to permit the differential readout of methylated cytosines. Bisulfite treatment damages the DNA, leading to fragmentation and loss of long-range methylation information. To overcome this limitation of bisulfite-treated DNA, we applied a new enzymatic deamination approach, termed enzymatic methyl-seq (EM-seq), to long-range sequencing technologies. Our methodology, named long-read enzymatic modification sequencing (LR-EM-seq), preserves the integrity of DNA, allowing long-range methylation profiling of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) over multikilobase length of genomic DNA. When applied to known differentially methylated regions (DMRs), LR-EM-seq achieves phasing of >5 kb, resulting in broader and better defined DMRs compared with that previously reported. This result showed the importance of phasing methylation for biologically relevant questions and the applicability of LR-EM-seq for long-range epigenetic analysis at single-molecule and single-nucleotide resolution.

Corrigendum: Identification and Characterization of Thermostable Y-Family DNA Polymerases η, ι, κ and Rev1 From a Lower Eukaryote, Thermomyces lanuginosus.

[This corrects the article DOI: 10.3389/fmolb.2021.778400.].

Novel ribonucleotide discrimination in the RNA polymerase-like two-barrel catalytic core of Family D DNA polymerases.

Family D DNA polymerase (PolD) is the essential replicative DNA polymerase for duplication of most archaeal genomes. PolD contains a unique two-barrel catalytic core absent from all other DNA polymerase families but found in RNA polymerases (RNAPs). While PolD has an ancestral RNA polymerase catalytic core, its active site has evolved the ability to discriminate against ribonucleotides. Until now, the mechanism evolved by PolD to prevent ribonucleotide incorporation was unknown. In all other DNA polymerase families, an active site steric gate residue prevents ribonucleotide incorporation. In this work, we identify two consensus active site acidic (a) and basic (b) motifs shared across the entire two-barrel nucleotide polymerase superfamily, and a nucleotide selectivity (s) motif specific to PolD versus RNAPs. A novel steric gate histidine residue (H931 in Thermococcus sp. 9°N PolD) in the PolD s-motif both prevents ribonucleotide incorporation and promotes efficient dNTP incorporation. Further, a PolD H931A steric gate mutant abolishes ribonucleotide discrimination and readily incorporates a variety of 2' modified nucleotides. Taken together, we construct the first putative nucleotide bound PolD active site model and provide structural and functional evidence for the emergence of DNA replication through the evolution of an ancestral RNAP two-barrel catalytic core.

RADAR-seq: A RAre DAmage and Repair sequencing method for detecting DNA damage on a genome-wide scale.

RAre DAmage and Repair sequencing (RADAR-seq) is a highly adaptable sequencing method that enables the identification and detection of rare DNA damage events for a wide variety of DNA lesions at single-molecule resolution on a genome-wide scale. In RADAR-seq, DNA lesions are replaced with a patch of modified bases that can be directly detected by Pacific Biosciences Single Molecule Real-Time (SMRT) sequencing. RADAR-seq enables dynamic detection over a wide range of DNA damage frequencies, including low physiological levels. Furthermore, without the need for DNA amplification and enrichment steps, RADAR-seq provides sequencing coverage of damaged and undamaged DNA across an entire genome. Here, we use RADAR-seq to measure the frequency and map the location of ribonucleotides in wild-type and RNaseH2-deficient E. coli and Thermococcus kodakarensis strains. Additionally, by tracking ribonucleotides incorporated during in vivo lagging strand DNA synthesis, we determined the replication initiation point in E. coli, and its relation to the origin of replication (oriC). RADAR-seq was also used to map cyclobutane pyrimidine dimers (CPDs) in Escherichia coli (E. coli) genomic DNA exposed to UV-radiation. On a broader scale, RADAR-seq can be applied to understand formation and repair of DNA damage, the correlation between DNA damage and disease initiation and progression, and complex biological pathways, including DNA replication.

Comprehensive Profiling of Four Base Overhang Ligation Fidelity by T4 DNA Ligase and Application to DNA Assembly.

Synthetic biology relies on the manufacture of large and complex DNA constructs from libraries of genetic parts. Golden Gate and other Type IIS restriction enzyme-dependent DNA assembly methods enable rapid construction of genes and operons through one-pot, multifragment assembly, with the ordering of parts determined by the ligation of Watson-Crick base-paired overhangs. However, ligation of mismatched overhangs leads to erroneous assembly, and low-efficiency Watson Crick pairings can lead to truncated assemblies. Using sets of empirically vetted, high-accuracy junction pairs avoids this issue but limits the number of parts that can be joined in a single reaction. Here, we report the use of comprehensive end-joining ligation fidelity and bias data to predict high accuracy junction sets for Golden Gate assembly. The ligation profile accurately predicted junction fidelity in ten-fragment Golden Gate assembly reactions and enabled accurate and efficient assembly of a lac cassette from up to 24-fragments in a single reaction.

Response to Comment on "DNA damage is a pervasive cause of sequencing errors, directly confounding variant identification".

Following the Comment of Stewart et al, we repeated our analysis on sequencing runs from The Cancer Genome Atlas (TCGA) using their suggested parameters. We found signs of oxidative damage in all sequence contexts and irrespective of the sequencing date, reaffirming that DNA damage affects mutation-calling pipelines in their ability to accurately identify somatic variations.

Single-stranded binding proteins and helicase enhance the activity of prokaryotic argonautes in vitro.

Prokaryotic argonautes are a unique class of nucleic acid-guided endonucleases putatively involved in cellular defense against foreign genetic elements. While their eukaryotic homologs and Cas protein counterparts require single-stranded RNAs as guides, some prokaryotic argonautes are able to utilize short single-stranded DNAs as guides for sequence-specific endonuclease activity. Many complications currently prevent the use of prokaryotic argonautes for in vivo gene-editing applications; however, they do exhibit potential as a new class of in vitro molecular tools if certain challenges can be overcome, specifically the limitations on substrate accessibility which leads to unequal levels of activity across a broad palate of substrates and the inability to act on double-stranded DNA substrates. Here we demonstrate the use of accessory factors, including thermostable single-stranded DNA binding proteins and UvrD-like helicase, in conjunction with prokaryotic argonautes to significantly improve enzymatic activity and enable functionality with a broader range of substrates, including linear double-stranded DNA substrates. We also demonstrate the use of Thermus thermophilus argonaute with accessory factors as a programmable restriction enzyme to generate long, unique single-stranded overhangs from linear double-stranded substrates compatible with downstream ligation.

A single-molecule sequencing assay for the comprehensive profiling of T4 DNA ligase fidelity and bias during DNA end-joining.

DNA ligases are key enzymes in molecular and synthetic biology that catalyze the joining of breaks in duplex DNA and the end-joining of DNA fragments. Ligation fidelity (discrimination against the ligation of substrates containing mismatched base pairs) and bias (preferential ligation of particular sequences over others) have been well-studied in the context of nick ligation. However, almost no data exist for fidelity and bias in end-joining ligation contexts. In this study, we applied Pacific Biosciences Single-Molecule Real-Time sequencing technology to directly sequence the products of a highly multiplexed ligation reaction. This method has been used to profile the ligation of all three-base 5'-overhangs by T4 DNA ligase under typical ligation conditions in a single experiment. We report the relative frequency of all ligation products with or without mismatches, the position-dependent frequency of each mismatch, and the surprising observation that 5'-TNA overhangs ligate extremely inefficiently compared to all other Watson-Crick pairings. The method can easily be extended to profile other ligases, end-types (e.g. blunt ends and overhangs of different lengths), and the effect of adjacent sequence on the ligation results. Further, the method has the potential to provide new insights into the thermodynamics of annealing and the kinetics of end-joining reactions.

Comparative analysis of the end-joining activity of several DNA ligases.

DNA ligases catalyze the repair of phosphate backbone breaks in DNA, acting with highest activity on breaks in one strand of duplex DNA. Some DNA ligases have also been observed to ligate two DNA fragments with short complementary overhangs or blunt-ended termini. In this study, several wild-type DNA ligases (phage T3, T4, and T7 DNA ligases, Paramecium bursaria chlorella virus 1 (PBCV1) DNA ligase, human DNA ligase 3, and Escherichia coli DNA ligase) were tested for their ability to ligate DNA fragments with several difficult to ligate end structures (blunt-ended termini, 3'- and 5'- single base overhangs, and 5'-two base overhangs). This analysis revealed that T4 DNA ligase, the most common enzyme utilized for in vitro ligation, had its greatest activity on blunt- and 2-base overhangs, and poorest on 5'-single base overhangs. Other ligases had different substrate specificity: T3 DNA ligase ligated only blunt ends well; PBCV1 DNA ligase joined 3'-single base overhangs and 2-base overhangs effectively with little blunt or 5'- single base overhang activity; and human ligase 3 had highest activity on blunt ends and 5'-single base overhangs. There is no correlation of activity among ligases on blunt DNA ends with their activity on single base overhangs. In addition, DNA binding domains (Sso7d, hLig3 zinc finger, and T4 DNA ligase N-terminal domain) were fused to PBCV1 DNA ligase to explore whether modified binding to DNA would lead to greater activity on these difficult to ligate substrates. These engineered ligases showed both an increased binding affinity for DNA and increased activity, but did not alter the relative substrate preferences of PBCV1 DNA ligase, indicating active site structure plays a role in determining substrate preference.

Colorimetric tests for diagnosis of filarial infection and vector surveillance using non-instrumented nucleic acid loop-mediated isothermal amplification (NINA-LAMP).

Accurate detection of filarial parasites in humans is essential for the implementation and evaluation of mass drug administration programs to control onchocerciasis and lymphatic filariasis. Determining the infection levels in vector populations is also important for assessing transmission, deciding when drug treatments may be terminated and for monitoring recrudescence. Immunological methods to detect infection in humans are available, however, cross-reactivity issues have been reported. Nucleic acid-based molecular assays offer high levels of specificity and sensitivity, and can be used to detect infection in both humans and vectors. In this study we developed loop-mediated isothermal amplification (LAMP) tests to detect three different filarial DNAs in human and insect samples using pH sensitive dyes for enhanced visual detection of amplification. Furthermore, reactions were performed in a portable, non-instrumented nucleic acid amplification (NINA) device that provides a stable heat source for LAMP. The efficacy of several strand displacing DNA polymerases were evaluated in combination with neutral red or phenol red dyes. Colorimetric NINA-LAMP assays targeting Brugia Hha I repeat, Onchocerca volvulus GST1a and Wuchereria bancrofti LDR each exhibit species-specificity and are also highly sensitive, detecting DNA equivalent to 1/10-1/5000th of one microfilaria. Reaction times varied depending on whether a single copy gene (70 minutes, O. volvulus) or repetitive DNA (40 min, B. malayi and W. bancrofti) was employed as a biomarker. The NINA heater can be used to detect multiple infections simultaneously. The accuracy, simplicity and versatility of the technology suggests that colorimetric NINA-LAMP assays are ideally suited for monitoring the success of filariasis control programs.

DNA damage is a pervasive cause of sequencing errors, directly confounding variant identification.

Mutations in somatic cells generate a heterogeneous genomic population and may result in serious medical conditions. Although cancer is typically associated with somatic variations, advances in DNA sequencing indicate that cell-specific variants affect a number of phenotypes and pathologies. Here, we show that mutagenic damage accounts for the majority of the erroneous identification of variants with low to moderate (1 to 5%) frequency. More important, we found signatures of damage in most sequencing data sets in widely used resources, including the 1000 Genomes Project and The Cancer Genome Atlas, establishing damage as a pervasive cause of sequencing errors. The extent of this damage directly confounds the determination of somatic variants in these data sets.

Rapid Time Scale Analysis of T4 DNA Ligase-DNA Binding.

DNA ligases, essential to both in vivo genome integrity and in vitro molecular biology, catalyze phosphodiester bond formation between adjacent 3'-OH and 5'-phosphorylated termini in dsDNA. This reaction requires enzyme self-adenylylation, using ATP or NAD+ as a cofactor, and AMP release concomitant with phosphodiester bond formation. In this study, we present the first fast time scale binding kinetics of T4 DNA ligase to both nicked substrate DNA (nDNA) and product-equivalent non-nicked dsDNA, as well as binding and release kinetics of AMP. The described assays utilized a fluorescein-dT labeled DNA substrate as a reporter for ligase·DNA interactions via stopped-flow fluorescence spectroscopy. The analysis revealed that binding to nDNA by the active adenylylated ligase occurs in two steps, an initial rapid association equilibrium followed by a transition to a second bound state prior to catalysis. Furthermore, the ligase binds and dissociates from nicked and nonsubstrate dsDNA rapidly with initial association affinities on the order of 100 nM regardless of enzyme adenylylation state. DNA binding occurs through a two-step mechanism in all cases, confirming prior proposals of transient binding followed by a transition to a productive ligase·nDNA (Lig·nDNA) conformation but suggesting that weaker nonproductive "closed" complexes are formed as well. These observations demonstrate the mechanistic underpinnings of competitive inhibition by rapid binding of nonsubstrate DNA, and of substrate inhibition by blocking of the self-adenylylation reaction through nick binding by deadenylylated ligase. Our analysis further reveals that product release is not the rate-determining step in turnover.

The Inhibitory Effect of Non-Substrate and Substrate DNA on the Ligation and Self-Adenylylation Reactions Catalyzed by T4 DNA Ligase.

DNA ligases are essential both to in vivo replication, repair and recombination processes, and in vitro molecular biology protocols. Prior characterization of DNA ligases through gel shift assays has shown the presence of a nick site to be essential for tight binding between the enzyme and its dsDNA substrate, with no interaction evident on dsDNA lacking a nick. In the current study, we observed a significant substrate inhibition effect, as well as the inhibition of both the self-adenylylation and nick-sealing steps of T4 DNA ligase by non-nicked, non-substrate dsDNA. Inhibition by non-substrate DNA was dependent only on the total DNA concentration rather than the structure; with 1 µg/mL of 40-mers, 75-mers, or circular plasmid DNA all inhibiting ligation equally. A >15-fold reduction in T4 DNA ligase self-adenylylation rate when in the presence of high non-nicked dsDNA concentrations was observed. Finally, EMSAs were utilized to demonstrate that non-substrate dsDNA can compete with nicked dsDNA substrates for enzyme binding. Based upon these data, we hypothesize the inhibition of T4 DNA ligase by non-nicked dsDNA is direct evidence for a two-step nick-binding mechanism, with an initial, nick-independent, transient dsDNA-binding event preceding a transition to a stable binding complex in the presence of a nick site.

A high-throughput assay for the comprehensive profiling of DNA ligase fidelity.

DNA ligases have broad application in molecular biology, from traditional cloning methods to modern synthetic biology and molecular diagnostics protocols. Ligation-based detection of polynucleotide sequences can be achieved by the ligation of probe oligonucleotides when annealed to a complementary target sequence. In order to achieve a high sensitivity and low background, the ligase must efficiently join correctly base-paired substrates, while discriminating against the ligation of substrates containing even one mismatched base pair. In the current study, we report the use of capillary electrophoresis to rapidly generate mismatch fidelity profiles that interrogate all 256 possible base-pair combinations at a ligation junction in a single experiment. Rapid screening of ligase fidelity in a 96-well plate format has allowed the study of ligase fidelity in unprecedented depth. As an example of this new method, herein we report the ligation fidelity of Thermus thermophilus DNA ligase at a range of temperatures, buffer pH and monovalent cation strength. This screen allows the selection of reaction conditions that maximize fidelity without sacrificing activity, while generating a profile of specific mismatches that ligate detectably under each set of conditions.

Modifiers of solid RNP granules control normal RNP dynamics and mRNA activity in early development.

Ribonucleoproteins (RNPs) often coassemble into supramolecular bodies with regulated dynamics. The factors controlling RNP bodies and connections to RNA regulation are unclear. During Caenorhabditis elegans oogenesis, cytoplasmic RNPs can transition among diffuse, liquid, and solid states linked to mRNA regulation. Loss of CGH-1/Ddx6 RNA helicase generates solid granules that are sensitive to mRNA regulators. Here, we identified 66 modifiers of RNP solids induced by cgh-1 mutation. A majority of genes promote or suppress normal RNP body assembly, dynamics, or metabolism. Surprisingly, polyadenylation factors promote RNP coassembly in vivo, suggesting new functions of poly(A) tail regulation in RNP dynamics. Many genes carry polyglutatmine (polyQ) motifs or modulate polyQ aggregation, indicating possible connections with neurodegenerative disorders induced by CAG/polyQ expansion. Several RNP body regulators repress translation of mRNA subsets, suggesting that mRNAs are repressed by multiple mechanisms. Collectively, these findings suggest new pathways of RNP modification that control large-scale coassembly and mRNA activity during development.

Genome Filtering for New DNA Biomarkers of Loa loa Infection Suitable for Loop-Mediated Isothermal Amplification.

Loa loa infections have emerged as a serious public health problem in patients co-infected with Onchocerca volvulus or Wuchereria bancrofti because of severe adverse neurological reactions after treatment with ivermectin. Accurate diagnostic tests are needed for careful mapping in regions where mass drug administration is underway. Loop-mediated isothermal amplification (LAMP) has become a widely adopted screening method because of its operational simplicity, rapidity and versatility of visual detection readout options. Here, we present a multi-step bioinformatic pipeline to generate diagnostic candidates suitable for LAMP and experimentally validate this approach using one of the identified candidates to develop a species-specific LAMP assay for L. loa. The pipeline identified ~140 new L. loa specific DNA repeat families as putative biomarkers of infection. The consensus sequence of one family, repeat family 4 (RF4), was compiled from ~ 350 sequences dispersed throughout the L. loa genome and maps to a L. loa-specific region of the long terminal repeats found at the boundaries of Bel/Pao retrotransposons. PCR and LAMP primer sets targeting RF4 specifically amplified L. loa but not W. bancrofti, O. volvulus, Brugia malayi, human or mosquito DNA. RF4 LAMP detects the DNA equivalent of one microfilaria (100 pg) in 25-30 minutes and as little as 0.060 pg of L. loa DNA (~1/1600th of a microfilaria) purified from spiked blood samples in approximately 50 minutes. In summary, we have successfully employed a bioinformatic approach to mine the L. loa genome for species-specific repeat families that can serve as new DNA biomarkers for LAMP. The RF4 LAMP assay shows promise as a field tool for the implementation and management of mass drug administration programs and warrants further testing on clinical samples as the next stage in development towards this goal.

Visual detection of isothermal nucleic acid amplification using pH-sensitive dyes.

Nucleic acid amplification is the basis for many molecular diagnostic assays. In these cases, the amplification product must be detected and analyzed, typically requiring extended workflow time, sophisticated equipment, or both. Here we present a novel method of amplification detection that harnesses the pH change resulting from amplification reactions performed with minimal buffering capacity. In loop-mediated isothermal amplification (LAMP) reactions, we achieved rapid (<30 min) and sensitive (<10 copies) visual detection using pH-sensitive dyes. Additionally, the detection can be performed in real time, enabling high-throughput or quantitative applications. We also demonstrate this visual detection for another isothermal amplification method (strand-displacement amplification), PCR, and reverse transcription LAMP (RT-LAMP) detection of RNA. The colorimetric detection of amplification presented here represents a generally applicable approach for visual detection of nucleic acid amplification, enabling molecular diagnostic tests to be analyzed immediately without the need for specialized and expensive instrumentation.