Next-generation sequencing-based analysis of reverse transcriptase fidelity

Enhancing thermostability of Moloney murine leukemia virus reverse transcriptase through greedy combination of multiple mutant residues

Reverse transcription is crucial in bioengineering and biomedical fields, particularly for genome sequencing and virus diagnosis. Enhancing the thermostability of reverse transcriptase can significantly improve its efficiency and accuracy by enabling it to function at higher temperatures, thereby reducing RNA secondary structures and minimizing interference from contaminating enzymes, particularly in clinical samples. Here, using a combinatorial strategy, a variant of Moloney Murine Leukemia Virus reverse transcriptase (MMLV RT) with improved activity across a wide temperature range (30-50 °C) was identified and maintained 100% activity after incubation at 50 °C for 10 min. Eleven hot-spot residues were mutated in various combinations, and the mutant proteins were rapidly expressed in a cell-free system for reverse transcription activity testing. Variant M5, which carries five mutated residues (E47K/E280R/T284R/L413G/D631V), exhibited enhanced thermostability and activity compared to any variant with a single residue mutation. Using purified recombinant protein for precise characterization, the melting temperature (Tm) of M5 increased by 4.7 °C when assembled with a nucleotide template-primer (T/P). Consequently, the half-life of M5 at 50 °C extended to approximately 60 min, in contrast to less than 4 min for the wild type. These findings demonstrate that the epistasis of combining multiple mutant residues holds excellent potential for significantly enhancing enzyme activity, even with existing knowledge. This heat-stable MMLV RT variant M5 may potentially improve efficiency and accuracy in molecular biology research and clinical diagnostics.

Next-Generation Sequencing Methods to Determine the Accuracy of Retroviral Reverse Transcriptases: Advantages and Limitations

Retroviruses, like other RNA viruses, mutate at very high rates and exist as genetically heterogeneous populations. The error-prone activity of viral reverse transcriptase (RT) is largely responsible for the observed variability, most notably in HIV-1. In addition, RTs are widely used in biotechnology to detect RNAs and to clone expressed genes, among many other applications. The fidelity of retroviral RTs has been traditionally analyzed using enzymatic (gel-based) or reporter-based assays. However, these methods are laborious and have important limitations. The development of next-generation sequencing (NGS) technologies opened the possibility of obtaining reverse transcription error rates from a large number of sequences, although appropriate protocols had to be developed. In this review, we summarize the developments in this field that allowed the determination of RNA-dependent DNA synthesis error rates for different RTs (viral and non-viral), including methods such as PRIMER IDs, REP-SEQ, ARC-SEQ, CIR-SEQ, SMRT-SEQ and ROLL-SEQ. Their advantages and limitations are discussed. Complementary DNA (cDNA) synthesis error rates obtained in different studies, using RTs and RNAs of diverse origins, are presented and compared. Future improvements in methodological pipelines will be needed for the precise identification of mutations in the RNA template, including modified bases.

Environmental RNA metabarcoding for ballast water microbial diversity: Minimizing false positives

While maritime transport boosts global trade by shipping bulk goods, it raises concerns about the spread of harmful bacteria via ballast water. Moreover, the dark and cold environments of ballast tanks often harbor extracellular DNA from dead organisms, leading to false positives in traditional environmental DNA (eDNA) metabarcoding analyses. Here, we alternatively employed environmental RNA (eRNA) metabarcoding to assess its potential for reducing false positive in ballast water monitoring. We collected eDNA and eRNA in parallel from ballast water before and after disinfection in three vessels. High-throughput sequencing of the 16S rRNA V4-V5 regions and cDNA counterparts was conducted to compare bacterial community composition. Our findings showed that over 80 % of the top 150 abundant amplicon sequence variants (ASVs) were detected by both eRNA and eDNA metabarcoding. Samples sequenced separately using DNA and RNA consistently clustered together, indicating similar community recovery efficacy. However, 42 % of ASVs were detected exclusively in DNA, resulting in significantly higher bacterial diversity compared to RNA, which suggests false positives in the DNA method. In treated samples with higher dead cell counts, the RNA method showed significantly lower bacterial diversity, indicating its effectiveness in detecting live bacteria. In summary, eRNA metabarcoding offers comparable recovery efficiency while maintaining a lower false-positive rate than eDNA metabarcoding.

HIV-1 Reverse Transcriptase Error Rates and Transcriptional Thresholds Based on Single-strand Consensus Sequencing of Target RNA Derived From In Vitro-transcription and HIV-infected Cells

Nucleotide incorporation and lacZ-based forward mutation assays have been widely used to determine the accuracy of reverse transcriptases (RTs) in RNA-dependent DNA polymerization reactions. However, they involve quite complex and laborious procedures, and cannot provide accurate error rates. Recently, NGS-based methods using barcodes opened the possibility of detecting all errors introduced by the RT, although their widespread use is limited by cost, due to the large size of libraries to be sequenced. In this study, we describe a novel and relatively simple NGS assay based on single-strand consensus sequencing that provides robust results with a relatively small number of raw sequences (around 60 Mb). The method has been validated by determining the error rate of HIV-1 (BH10 strain) RT using the HIV-1 protease-coding sequence as target. HIV-1 reverse transcription error rates in standard conditions (37 °C/3 mM Mg2+) using an in vitro-transcribed RNA were around 7.3 × 10-5. In agreement with previous reports, an 8-fold increase in RT's accuracy was observed after reducing Mg2+ concentration to 0.5 mM. The fidelity of HIV-1 RT was also higher at 50 °C than at 37 °C (error rate 1.5 × 10-5). Interestingly, error rates obtained with HIV-1 RNA from infected cells as template of the reverse transcription at 3 mM Mg2+ (7.4 × 10-5) were similar to those determined with the in vitro-transcribed RNA, and were reduced to 1.8 × 10-5 in the presence of 0.5 mM Mg2+. Values obtained at low magnesium concentrations were modestly higher than the transcription error rates calculated for human cells, thereby suggesting a realistic transcriptional threshold for our NGS-based error rate determinations.

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 using 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.

Digital RNA sequencing using unique molecular identifiers enables ultrasensitive RNA mutation analysis

Mutation analysis is typically performed at the DNA level since most technical approaches are developed for DNA analysis. However, some applications, like transcriptional mutagenesis, RNA editing and gene expression analysis, require RNA analysis. Here, we combine reverse transcription and digital DNA sequencing to enable low error digital RNA sequencing. We evaluate yield, reproducibility, dynamic range and error correction rate for seven different reverse transcription conditions using multiplexed assays. The yield, reproducibility and error rate vary substantially between the specific conditions, where the yield differs 9.9-fold between the best and worst performing condition. Next, we show that error rates similar to DNA sequencing can be achieved for RNA using appropriate reverse transcription conditions, enabling detection of mutant allele frequencies <0.1% at RNA level. We also detect mutations at both DNA and RNA levels in tumor tissue using a breast cancer panel. Finally, we demonstrate that digital RNA sequencing can be applied to liquid biopsies, analyzing cell-free gene transcripts. In conclusion, we demonstrate that digital RNA sequencing is suitable for ultrasensitive RNA mutation analysis, enabling several basic research and clinical applications.

Inosine-Induced Base Pairing Diversity during Reverse Transcription

A-to-I editing catalyzed by adenosine deaminase acting on RNAs impacts numerous physiological and biochemical processes that are essential for cellular functions and is a big contributor to the infectivity of certain RNA viruses. The outcome of this deamination leads to changes in the eukaryotic transcriptome functionally resembling A-G transitions since inosine preferentially pairs with cytosine. Moreover, hyper-editing or multiple A to G transitions in clusters were detected in measles virus. Inosine modifications either directly on viral RNA or on cellular RNA can have antiviral or pro-viral repercussions. While many of the significant roles of inosine in cellular RNAs are well understood, the effects of hyper-editing of A to I on viral polymerase activity during RNA replication remain elusive. Moreover, biological strategies such as molecular cloning and RNA-seq for transcriptomic interrogation rely on RT-polymerase chain reaction with little to no emphasis placed on the first step, reverse transcription, which may reshape the sequencing results when hypermodification is present. In this study, we systematically explore the influence of inosine modification, varying the number and position of inosines, on decoding outcomes using three different reverse transcriptases (RTs) followed by standard Sanger sequencing. We find that inosine alone or in clusters can differentially affect the RT activity. To gain structural insights into the accommodation of inosine in the polymerase site of HIV-1 reverse transcriptase (HIV-1-RT) and how this structural context affects the base pairing rules for inosine, we performed molecular dynamics simulations of the HIV-1-RT. The simulations highlight the importance of the protein-nucleotide interaction as a critical factor in deciphering the base pairing behavior of inosine clusters. This effort sets the groundwork for decrypting the physiological significance of inosine and linking the fidelity of reverse transcriptase and the possible diverse transcription outcomes of cellular RNAs and/or viral RNAs where hyper-edited inosines are present in the transcripts.

M-MuLV reverse transcriptase: Selected properties and improved mutants

Reverse transcriptases (RTs) are enzymes synthesizing DNA using RNA as the template and serving as the standard tools in modern biotechnology and molecular diagnostics. To date, the most commonly used reverse transcriptase is the enzyme from Moloney murine leukemia virus, M-MuLV RT. Since its discovery, M-MuLV RT has become indispensable for modern RNA studies; the range of M-MuLV RT applications is vast, from scientific tasks to clinical testing of human pathogens. This review will give a brief description of the structure, thermal stability, processivity, and fidelity, focusing on improving M-MuLV RT for practical usage.

Physiological magnesium concentrations increase fidelity of diverse reverse transcriptases from HIV-1, HIV-2, and foamy virus, but not MuLV or AMV

Reverse transcriptases (RTs) are typically assayed using optimized Mg2+ concentrations (~5-10 mM) several-fold higher than physiological cellular free Mg2+ (~0.5 mM). Recent analyses demonstrated that HIV-1, but not Moloney murine leukaemia (MuLV) or avain myeloblastosis (AMV) virus RTs has higher fidelity in low Mg2+. In the current report, lacZα-based α-complementation assays were used to measure the fidelity of several RTs including HIV-1 (subtype B and A/E), several drug-resistant HIV-1 derivatives, HIV-2, and prototype foamy virus (PFV), all which showed higher fidelity using physiological Mg2+, while MuLV and AMV RTs demonstrated equivalent fidelity in low and high Mg2+. In 0.5 mM Mg2+, all RTs demonstrated approximately equal fidelity, except for PFV which showed higher fidelity. A Next Generation Sequencing (NGS) approach that used barcoding to determine mutation profiles was used to examine the types of mutations made by HIV-1 RT (type B) in low (0.5 mM) and high (6 mM) Mg2+ on a lacZα template. Unlike α-complementation assays which are dependent on LacZα activity, the NGS assay scores mutations at all positions and of every type. Consistent with α-complementation assays, a ~four-fold increase in mutations was observed in high Mg2+. These findings help explain why HIV-1 RT displays lower fidelity in vitro (with high Mg2+ concentrations) than other RTs (e.g. MuLV and AMV), yet cellular fidelity for these viruses is comparable. Establishing in vitro conditions that accurately represent RT's activity in cells is pivotal to determining the contribution of RT and other factors to the mutation profile observed with HIV-1.

High Throughput Sequencing for the Detection and Characterization of RNA Viruses

This review aims to assess and recommend approaches for targeted and agnostic High Throughput Sequencing of RNA viruses in a variety of sample matrices. HTS also referred to as deep sequencing, next generation sequencing and third generation sequencing; has much to offer to the field of environmental virology as its increased sequencing depth circumvents issues with cloning environmental isolates for Sanger sequencing. That said however, it is important to consider the challenges and biases that method choice can impart to sequencing results. Here, methodology choices from RNA extraction, reverse transcription to library preparation are compared based on their impact on the detection or characterization of RNA viruses.

Improvement of Moloney murine leukemia virus reverse transcriptase thermostability by introducing a disulfide bridge in the ribonuclease H region

Moloney murine leukemia virus (MMLV) reverse transcriptase (RT) is widely used in research and clinical diagnosis. Improvement of MMLV RT thermostability has been an important topic of research for increasing the efficiency of cDNA synthesis. In this study, we attempted to increase MMLV RT thermostability by introducing a disulfide bridge in its RNase H region using site-directed mutagenesis. Five variants were designed, focusing on the distance between the two residues to be mutated into cysteine. The variants were expressed in Escherichia coli and purified. A551C/T662C was determined to be the most thermostable variant.

Navigating the Pandemic Response Life Cycle: Molecular Diagnostics and Immunoassays in the Context of COVID-19 Management

Coronavirus disease 2019 (COVID-19) is an infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). To counter COVID-19 spreading, an infrastructure to provide rapid and thorough molecular diagnostics and serology testing is the cornerstone of outbreak and pandemic management. We hereby review the clinical insights with regard to using molecular tests and immunoassays in the context of COVID-19 management life cycle: the preventive phase, the preparedness phase, the response phase and the recovery phase. The spatial and temporal distribution of viral RNA, antigens and antibodies during human infection is summarized to provide a biological foundation for accurate detection of the disease. We shared the lessons learned and the obstacles encountered during real world high-volume screening programs. Clinical needs are discussed to identify existing technology gaps in these tests. Leverage technologies, such as engineered polymerases, isothermal amplification, and direct amplification from complex matrices may improve the productivity of current infrastructure, while emerging technologies like CRISPR diagnostics, visual end point detection, and PCR free methods for nucleic acid sensing may lead to at-home tests. The lessons learned, and innovations spurred from the COVID-19 pandemic could upgrade our global public health infrastructure to better combat potential outbreaks in the future.

Alteration of enzymes and their application to nucleic acid amplification (Review)

Since the discovery of polymerase chain reaction (PCR) in 1985, several methods have been developed to achieve nucleic acid amplification, and are currently used in various fields including clinical diagnosis and life science research. Thus, a wealth of information has accumulated regarding nucleic acid-related enzymes. In this review, some nucleic acid-related enzymes were selected and the recent advances in their modification along with their application to nucleic acid amplification were described. The discussion also focused on optimization of the corresponding reaction conditions. Using newly developed enzymes under well-optimized reaction conditions, the sensitivity, specificity, and fidelity of nucleic acid tests can be improved successfully.

Discovery and evolution of RNA and XNA reverse transcriptase function and fidelity

The ability of reverse transcriptases (RTs) to synthesize a complementary DNA from natural RNA and a range of unnatural xeno nucleic acid (XNA) template chemistries, underpins key methods in molecular and synthetic genetics. However, RTs have proven challenging to discover and engineer, in particular for the more divergent XNA chemistries. Here we describe a general strategy for the directed evolution of RT function for any template chemistry called compartmentalized bead labelling and demonstrate it by the directed evolution of efficient RTs for 2'-O-methyl RNA and hexitol nucleic acids and the discovery of RTs for the orphan XNA chemistries D-altritol nucleic acid and 2'-methoxyethyl RNA, for which previously no RTs existed. Finally, we describe the engineering of XNA RTs with active exonucleolytic proofreading as well as the directed evolution of RNA RTs with very high complementary DNA synthesis fidelities, even in the absence of proofreading.

Reverse Transcriptase: From Transcriptomics to Genome Editing

Reverse transcriptases (RTs) are enzymes that can generate a complementary strand of DNA (cDNA) from RNA. Coupled with PCR, RTs have been widely used to detect RNAs and to clone expressed genes. Classical retroviral RTs have been improved by protein engineering. These enzymes and newly characterized RTs are key elements in the development of next-generation sequencing techniques that are now being applied to the study of transcriptomics. In addition, engineered RTs fused to a CRISPR/Cas9 nickase have recently shown great potential as tools to manipulate eukaryotic genomes. In this review, we discuss the properties and uses of wild type and engineered RTs in biotechnological applications, from conventional RT-PCR to recently introduced prime editing.

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 (m⁶A, m⁵C and m⁵U), hydroxymethylated (hm⁵U) 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 N⁶-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.

Decreased Km to dNTPs is an essential M-MuLV reverse transcriptase adoption required to perform efficient cDNA synthesis in One-Step RT-PCR assay

Personalized medicine and advanced diagnostic tools based on RNA analysis are focusing on fast and direct One-Step RT-PCR assays. First strand complementary DNA (cDNA) synthesized by the reverse transcriptase (RT) is exponentially amplified in the end-point or real-time PCR. Even a minor discrepancy in PCR conditions would result in big deviations during the data analysis. Thus, One-Step RT-PCR composition is typically based on the PCR buffer. In this study, we have used compartmentalized ribosome display technique for in vitro evolution of the Moloney Murine Leukemia Virus reverse transcriptase (M-MuLV RT) that would be able to perform efficient full-length cDNA synthesis in PCR buffer optimized for Thermus aquaticus DNA polymerase. The most frequent mutations found in a selected library were analyzed. Aside from the mutations, which switch off RNase H activity of RT and are beneficial for the full-length cDNA synthesis, we have identified several mutations in the active center of the enzyme (Q221R and V223A/M), which result in 4-5-fold decrease of Km for dNTPs (<0.2 mM). The selected mutations are in surprising agreement with the natural evolution process because they transformed the active center from the oncoretroviral M-MuLV RT-type to the lenitiviral enzyme-type. We believe that this was the major and essential phenotypic adjustment required to perform fast and efficient cDNA synthesis in PCR buffer at 0.2-mM concentration of each dNTP.

High sensitive RNA detection by one-step RT-PCR using the genetically engineered variant of DNA polymerase with reverse transcriptase activity from hyperthermophilies

One-step RT-PCR has not been widely used even though some thermostable DNA polymerases with reverse transcriptase (RT) activity were developed from bacterial and archaeal polymerases, which is owing to low cDNA synthesis activity from RNA. In the present study, we developed highly-sensitive one-step RT-PCR using the single variant of family A DNA polymerase with RT activity, K4pol_L329A (L329A), from the hyperthermophilic bacterium Thermotoga petrophila K4 or the 16-tuple variant of family B DNA polymerase with RT activity, RTX, from the hyperthermophilic archaeon Thermococcus kodakarensis. Optimization of reaction condition revealed that the activities for cDNA synthesis and PCR of K4pol_L329A and RTX were highly affected by the concentrations of MgCl₂ and Mn(OCOCH₃)₂ as well as those of K4pol_L329A or RTX. Under the optimized condition, 300 copies/μl of target RNA in 10 μl reaction volumes were successfully detected by the one-step RT-PCR with K4pol_L329A or RTX, which was almost equally sensitive enough compared with the current RT-PCR condition using retroviral RT and thermostable DNA polymerase. Considering that K4pol_L329A and RTX are stable even at 90-100°C, our results suggest that the one-step RT-PCR with K4pol_L329A or RTX is more advantageous than the current one.

RNA/DNA structures recognized by RNase H2

Ribonuclease H (RNase H) [EC 3.1.26.4] is an enzyme that specifically degrades RNA from RNA/DNA hybrids. Since its discovery in 1969, the enzyme has been extensively studied for its catalytic mechanism and physiological role. RNase H has been classified into two major families, Type 1 and Type 2. Type 1 enzymes are designated RNase HI in prokaryotes and RNase H1 in eukaryotes, while Type 2 enzymes are designated RNase HII in prokaryotes and RNase H2 in eukaryotes. Type 2 enzymes are able to cleave the 5′-phosphodiester bond of one ribonucleotide embedded in a DNA double strand. Recent studies have shown that RNase H2 is involved in excision of a single ribonucleotide embedded in genomic DNA and removal of an R-loop formed in cells. It is also involved in double-strand break of DNA and its repair. In this review, we aim to outline the structures recognized by RNase H2.

Application of long single-stranded DNA donors in genome editing: generation and validation of mouse mutants

BACKGROUND

Recent advances in clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) genome editing have led to the use of long single-stranded DNA (lssDNA) molecules for generating conditional mutations. However, there is still limited available data on the efficiency and reliability of this method.

RESULTS

We generated conditional mouse alleles using lssDNA donor templates and performed extensive characterization of the resulting mutations. We observed that the use of lssDNA molecules as donors efficiently yielded founders bearing the conditional allele, with seven out of nine projects giving rise to modified alleles. However, rearranged alleles including nucleotide changes, indels, local rearrangements and additional integrations were also frequently generated by this method. Specifically, we found that alleles containing unexpected point mutations were found in three of the nine projects analyzed. Alleles originating from illegitimate repairs or partial integration of the donor were detected in eight projects. Furthermore, additional integrations of donor molecules were identified in four out of the seven projects analyzed by copy counting. This highlighted the requirement for a thorough allele validation by polymerase chain reaction, sequencing and copy counting of the mice generated through this method. We also demonstrated the feasibility of using lssDNA donors to generate thus far problematic point mutations distant from active CRISPR cutting sites by targeting two distinct genes (Gckr and Rims1). We propose a strategy to perform extensive quality control and validation of both types of mouse models generated using lssDNA donors.

CONCLUSION

lssDNA donors reproducibly generate conditional alleles and can be used to introduce point mutations away from CRISPR/Cas9 cutting sites in mice. However, our work demonstrates that thorough quality control of new models is essential prior to reliably experimenting with mice generated by this method. These advances in genome editing techniques shift the challenge of mutagenesis from generation to the validation of new mutant models.