Discrete measurements of RNA polymerase and reverse transcriptase fidelity reveal evolutionary tuning.

Direct methods for determining the fidelity of DNA polymerases are robust, with relatively little sample manipulation before sequencing. In contrast, methods for measuring RNA polymerase and reverse transcriptase fidelities are complicated by additional preparation steps that introduce ambiguity and error. Here, we describe a sequencing method, termed Roll-Seq, for simultaneously determining the individual fidelities of RNA polymerases and reverse transcriptases (RT) using Pacific Biosciences Single Molecule Real-Time sequencing. By employing reverse transcriptases with high rolling-circle activity, Roll-Seq generates long concatemeric cDNA from a circular RNA template. To discern the origin of a mutation, errors are recorded and determined to occur within a single concatemer (reverse transcriptase error) or all concatemers (RNA polymerase error) over the cDNA strand. We used Roll-Seq to measure the fidelities of T7 RNA polymerases, a Group II intron-encoded RT (Induro), and two LINE RTs (Fasciolopsis buski R2-RT and human LINE-1). Substitution rates for Induro and R2-RT are the same for cDNA and second strand synthesis while LINE-1 has 2.5-fold lower fidelity when performing second strand synthesis. Deletion and insertion rates increase for all RTs during second strand synthesis. In addition, we find that a structured RNA template impacts fidelity for both RNA polymerase and RT. The accuracy and precision of Roll-Seq enable this method to be applied as a complementary analysis to structural and mechanistic characterization of RNA polymerases and reverse transcriptases or as a screening method for RNAP and RT fidelity.

Base modifications affecting RNA polymerase and reverse transcriptase fidelity.

Ribonucleic acid (RNA) is capable of hosting a variety of chemically diverse modifications, in both naturally-occurring post-transcriptional modifications and artificial chemical modifications used to expand the functionality of RNA. However, few studies have addressed how base modifications affect RNA polymerase and reverse transcriptase activity and fidelity. Here, we describe the fidelity of RNA synthesis and reverse transcription of modified ribonucleotides using an assay based on Pacific Biosciences Single Molecule Real-Time sequencing. Several modified bases, including methylated (m6A, m5C and m5U), hydroxymethylated (hm5U) and isomeric bases (pseudouridine), were examined. By comparing each modified base to the equivalent unmodified RNA base, we can determine how the modification affected cumulative RNA polymerase and reverse transcriptase fidelity. 5-hydroxymethyluridine and N6-methyladenosine both increased the combined error rate of T7 RNA polymerase and reverse transcriptases, while pseudouridine specifically increased the error rate of RNA synthesis by T7 RNA polymerase. In addition, we examined the frequency, mutational spectrum and sequence context of reverse transcription errors on DNA templates from an analysis of second strand DNA synthesis.

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.

Detection and Quantitation of DNA Damage on a Genome-wide Scale Using RADAR-seq.

The formation and persistence of DNA damage can impact biological processes such as DNA replication and transcription. To maintain genome stability and integrity, organisms rely on robust DNA damage repair pathways. Techniques to detect and locate DNA damage sites across a genome enable an understanding of the consequences of DNA damage as well as how damage is repaired, which can have key diagnostic and therapeutic implications. Importantly, advancements in technology have enabled the development of high-throughput sequencing-based DNA damage detection methods. These methods require DNA enrichment or amplification steps that limit the ability to quantitate the DNA damage sites. Further, each of these methods is typically tailored to detect only a specific type of damage. RAre DAmage and Repair (RADAR) sequencing is a DNA sequencing workflow that overcomes these limitations and enables detection and quantitation of DNA damage sites in any organism on a genome-wide scale. RADAR-seq works by replacing DNA damage sites with a patch of modified bases that can be directly detected by Pacific Biosciences Single-Molecule Real Time sequencing. Here, we present three protocols that enable detection of thymine dimers and ribonucleotides in bacterial and archaeal genomes. Basic Protocol 1 enables construction of a reference genome required for RADAR-seq analyses. Basic Protocol 2 describes how to locate, quantitate, and compare thymine dimer levels in Escherichia coli exposed to varying amounts of UV light. Basic Protocol 3 describes how to locate, quantitate, and compare ribonucleotide levels in wild-type and ΔRNaseH2 Thermococcus kodakarensis. Importantly, all three protocols provide in-depth steps for data analysis. Together they serve as proof-of-principle experiments that will allow users to adapt the protocols to locate and quantitate a wide variety of DNA damage sites in any organism. © 2022 New England Biolabs. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Constructing a reference genome utilizing SMRT sequencing Basic Protocol 2: Mapping and quantitating genomic thymine dimer formation in untreated versus UV-irradiated E. coli using RADAR-seq Basic Protocol 3: Mapping and quantitating genomic ribonucleotide incorporation in wildtype versus ΔRNaseH2 T. kodakarensis using RADAR-seq.

N1-methyl-pseudouridine is incorporated with higher fidelity than pseudouridine in synthetic RNAs.

In vitro transcribed synthetic messenger RNAs (mRNAs) represent a novel therapeutic modality. To overcome the inherent immunogenicity, as well as to increase the therapeutic efficacy of the molecules, uridine analogs-such as pseudouridine (Ψ) and N1-methyl-pseudouridine (m1Ψ), are incorporated in the synthetic mRNA. To decipher the fidelity with which these modifications are incorporated during the in vitro transcription (IVT) process, we compared the incorporation fidelity of uridine analogs with different RNA polymerases. We demonstrate that m1Ψ is incorporated with higher fidelity than Ψ. The fidelity of nucleotide incorporation differs between RNA polymerases; however, the spectrum of mutations observed between the RNAPs is similar. We also show that the array of nucleotide misincorporation is not dependent on the template DNA sequence context and that the distribution of these misincorporated nucleotides is not localized to any specific region along the length of the RNA. Based on our findings, we introduce a novel method to improve uridine analog incorporation fidelity during IVT. Our proof-of-concept experiments for higher-fidelity incorporation of uridine analogs during IVT provide guidelines when choosing RNAPs for the generation of modified uridine-containing mRNAs in vitro.

A bi-specific lectin from the mushroom Boletopsis grisea and its application in glycoanalytical workflows.

The BLL lectin from the edible Japanese "Kurokawa" mushroom (Boletopsis leucomelaena) was previously reported to bind to N-glycans harboring terminal N-acetylglucosamine (GlcNAc) and to induce apoptosis in a leukemia cell line. However, its gene has not been reported. In this study, we used a transcriptomics-based workflow to identify a full-length transcript of a BLL functional ortholog (termed BGL) from Boletopsis grisea, a close North American relative of B. leucomelaena. The deduced amino acid sequence of BGL was an obvious member of fungal fruit body lectin family (Pfam PF07367), a highly conserved group of mushroom lectins with a preference for binding O-glycans harboring the Thomsen-Friedenreich antigen (TF-antigen; Galß1,3GalNAc-α-) and having two ligand binding sites. Functional characterization of recombinant BGL using glycan microarray analysis and surface plasmon resonance confirmed its ability to bind both the TF-antigen and ß-GlcNAc-terminated N-glycans. Structure-guided mutagenesis of BGL's two ligand binding clefts showed that one site is responsible for binding TF-antigen structures associated with O-glycans, whereas the second site specifically recognizes N-glycans with terminal ß-GlcNAc. Additionally, the two sites show no evidence of allosteric communication. Finally, mutant BGL proteins having single functional bindings site were used to enrich GlcNAc-capped N-glycans or mucin type O-glycopeptides from complex samples in glycomics and glycoproteomics analytical workflows.

Therminator DNA Polymerase: Modified Nucleotides and Unnatural Substrates.

A variant of 9°N DNA polymerase [Genbank ID (AAA88769.1)] with three mutations (D141A, E143A, A485L) and commercialized under the name "Therminator DNA polymerase" has the ability to incorporate a variety of modified nucleotide classes. This Review focuses on how Therminator DNA Polymerase has enabled new technologies in synthetic biology and DNA sequencing. In addition, we discuss mechanisms for increased modified nucleotide incorporation.

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.

Pol V-Mediated Translesion Synthesis Elicits Localized Untargeted Mutagenesis during Post-replicative Gap Repair.

In vivo, replication forks proceed beyond replication-blocking lesions by way of downstream repriming, generating daughter strand gaps that are subsequently processed by post-replicative repair pathways such as homologous recombination and translesion synthesis (TLS). The way these gaps are filled during TLS is presently unknown. The structure of gap repair synthesis was assessed by sequencing large collections of single DNA molecules that underwent specific TLS events in vivo. The higher error frequency of specialized relative to replicative polymerases allowed us to visualize gap-filling events at high resolution. Unexpectedly, the data reveal that a specialized polymerase, Pol V, synthesizes stretches of DNA both upstream and downstream of a site-specific DNA lesion. Pol V-mediated untargeted mutations are thus spread over several hundred nucleotides, strongly eliciting genetic instability on either side of a given lesion. Consequently, post-replicative gap repair may be a source of untargeted mutations critical for gene diversification in adaptation and evolution.

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.

Directed evolution of DNA polymerase, RNA polymerase and reverse transcriptase activity in a single polypeptide.

DNA polymerases enable key technologies in modern biology but for many applications, native polymerases are limited by their stringent substrate recognition. Here we describe short-patch compartmentalized self-replication (spCSR), a novel strategy to expand the substrate spectrum of polymerases in a targeted way. spCSR is based on the previously described CSR, but unlike CSR only a short region (a "patch") of the gene under investigation is diversified and replicated. This allows the selection of polymerases under conditions where catalytic activity and processivity are compromised to the extent that full self-replication is inefficient. We targeted two specific motifs involved in substrate recognition in the active site of DNA polymerase I from Thermus aquaticus (Taq) and selected for incorporation of both ribonucleotide- (NTP) and deoxyribonucleotide-triphosphates (dNTPs) using spCSR. This allowed the isolation of multiple variants of Taq with apparent dual substrate specificity. They were able to synthesize RNA, while still retaining essentially wild-type (wt) DNA polymerase activity as judged by PCR. One such mutant (AA40: E602V, A608V, I614M, E615G) was able to incorporate both NTPs and dNTPs with the same catalytic efficiency as the wt enzyme incorporates dNTPs. AA40 allowed the generation of mixed RNA-DNA amplification products in PCR demonstrating DNA polymerase, RNA polymerase as well as reverse transcriptase activity within the same polypeptide. Furthermore, AA40 displayed an expanded substrate spectrum towards other 2'-substituted nucleotides and was able to synthesize nucleic acid polymers in which each base bore a different 2'-substituent. Our results suggest that spCSR will be a powerful strategy for the generation of polymerases with altered substrate specificity for applications in nano- and biotechnology and in the enzymatic synthesis of antisense and RNAi probes.

Examining Sources of Error in PCR by Single-Molecule Sequencing.

Next-generation sequencing technology has enabled the detection of rare genetic or somatic mutations and contributed to our understanding of disease progression and evolution. However, many next-generation sequencing technologies first rely on DNA amplification, via the Polymerase Chain Reaction (PCR), as part of sample preparation workflows. Mistakes made during PCR appear in sequencing data and contribute to false mutations that can ultimately confound genetic analysis. In this report, a single-molecule sequencing assay was used to comprehensively catalog the different types of errors introduced during PCR, including polymerase misincorporation, structure-induced template-switching, PCR-mediated recombination and DNA damage. In addition to well-characterized polymerase base substitution errors, other sources of error were found to be equally prevalent. PCR-mediated recombination by Taq polymerase was observed at the single-molecule level, and surprisingly found to occur as frequently as polymerase base substitution errors, suggesting it may be an underappreciated source of error for multiplex amplification reactions. Inverted repeat structural elements in lacZ caused polymerase template-switching between the top and bottom strands during replication and the frequency of these events were measured for different polymerases. For very accurate polymerases, DNA damage introduced during temperature cycling, and not polymerase base substitution errors, appeared to be the major contributor toward mutations occurring in amplification products. In total, we analyzed PCR products at the single-molecule level and present here a more complete picture of the types of mistakes that occur during DNA amplification.

Correction: Examining Sources of Error in PCR by Single-Molecule Sequencing.

[This corrects the article DOI: 10.1371/journal.pone.0169774.].

Intravesicle isothermal DNA replication.

BACKGROUND: Bacterial and viral DNA replication was previously reconstituted in vitro from component parts 1234. Significant advances in building minimal cell-like structures also have been made recently 567. Combining the two approaches would further attempts to build a minimal cell-like structure capable of undergoing evolution by combining membrane encapsulation and genome replication. Towards this end, we attempted to use purified genomic replication protein components from thermophilic bacterial sources to copy strands of DNA isothermally within lipid vesicles. FINDINGS: Bacterial replication components (such as helicases and DNA polymerases) are compatible with methods for the generation of lipid vesicles. Encapsulation inside phospholipid vesicles does not inhibit the activity of bacterial DNA genome replication machinery. Further the described system is efficient at isothermally amplifying short segments of DNA within phospholipid vesicles. CONCLUSIONS: Herein we show that bacterial isothermal DNA replication machinery is functional inside of phospholipid vesicles, suggesting that replicating cellular mimics can be built from purified bacterial components.

Isolation of receptor-ligand pairs by capture of long-lived multivalent interaction complexes.

We have combined phage display and array screening for the rapid isolation of pairs of interacting polypeptides. Our strategy, named SAC (selection by avidity capture), is based on the avidity effect, the formation of highly stable complexes formed by multivalent interactions; in our case, between a receptor (multivalently displayed on phage) and a ligand (coexpressed as a multimeric fusion protein). Capture of the long-lived interaction complex allows the isolation of phage bearing cognate interaction pairs, as we demonstrate for a range of interactions, including Ab-antigen pairs and the rapamycin-dependent interaction of FKBP-12 and FRAP. Cognate phage are enriched by SAC up to 1000-fold and interacting pairs can be identified by array screening. Application of SAC to Ab-antigen interactions as a model system yielded over 140 specific Abs to a single antigen and 92 Abs to three different fetal human brain antigens in a single round of SAC each. Our results suggest that SAC should prove useful for the identification and study of receptor-ligand interactions in particular among extracellular proteins, as well as for the rapid generation of specific Abs to multiple antigens.