Intrinsic DNA synthesis fidelity of xenotropic murine leukemia virus-related virus reverse transcriptase

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.

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.

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.

Evolution of HIV-1 within untreated individuals and at the population scale in Uganda

HIV-1 undergoes multiple rounds of error-prone replication between transmission events, resulting in diverse viral populations within and among individuals. In addition, the virus experiences different selective pressures at multiple levels: during the course of infection, at transmission, and among individuals. Disentangling how these evolutionary forces shape the evolution of the virus at the population scale is important for understanding pathogenesis, how drug- and immune-escape variants are likely to spread in populations, and the development of preventive vaccines. To address this, we deep-sequenced two regions of the HIV-1 genome (p24 and gp41) from 34 longitudinally-sampled untreated individuals from Rakai District in Uganda, infected with subtypes A, D, and inter-subtype recombinants. This dataset substantially increases the availability of HIV-1 sequence data that spans multiple years of untreated infection, in particular for different geographical regions and viral subtypes. In line with previous studies, we estimated an approximately five-fold faster rate of evolution at the within-host compared to the population scale for both synonymous and nonsynonymous substitutions, and for all subtypes. We determined the extent to which this mismatch in evolutionary rates can be explained by the evolution of the virus towards population-level consensus, or the transmission of viruses similar to those that establish infection within individuals. Our findings indicate that both processes are likely to be important.

Transcriptional inaccuracy threshold attenuates differences in RNA-dependent DNA synthesis fidelity between retroviral reverse transcriptases

In M13mp2 lacZα forward mutation assays measuring intrinsic fidelity of DNA-dependent DNA synthesis, wild-type human immunodeficiency virus type 1 (HIV-1) RTs of group M/subtype B previously showed >10-fold higher error rates than murine leukaemia virus (MLV) and avian myeloblastosis virus (AMV) RTs. An adapted version of the assay was used to obtain error rates of RNA-dependent DNA synthesis for several RTs, including wild-type HIV-1BH10, HIV-1ESP49, AMV and MLV RTs, and the high-fidelity mutants of HIV-1ESP49 RT K65R and K65R/V75I. Our results showed that there were less than two-fold differences in fidelity between the studied RTs with error rates ranging within 2.5 × 10-5 and 3.5 × 10-5. These results were consistent with the existence of a transcriptional inaccuracy threshold, generated by the RNA polymerase while synthesizing the RNA template used in the assay. A modest but consistent reduction of the inaccuracy threshold was achieved by lowering the pH and Mg2+ concentration of the transcription reaction. Despite assay limitations, we conclude that HIV-1BH10 and HIV-1ESP49 RTs are less accurate when copying DNA templates than RNA templates. Analysis of the RNA-dependent mutational spectra revealed a higher tendency to introduce large deletions at the initiation of reverse transcription by all HIV-1 RTs except the double-mutant K65R/V75I.

Viral reverse transcriptases

Reverse transcriptases (RTs) play a major role in the replication of Retroviridae, Metaviridae, Pseudoviridae, Hepadnaviridae and Caulimoviridae. RTs are enzymes that are able to synthesize DNA using RNA or DNA as templates (DNA polymerase activity), and degrade RNA when forming RNA/DNA hybrids (ribonuclease H activity). In retroviruses and LTR retrotransposons (Metaviridae and Pseudoviridae), the coordinated action of both enzymatic activities converts single-stranded RNA into a double-stranded DNA that is flanked by identical sequences known as long terminal repeats (LTRs). RTs of retroviruses and LTR retrotransposons are active as monomers (e.g. murine leukemia virus RT), homodimers (e.g. Ty3 RT) or heterodimers (e.g. human immunodeficiency virus type 1 (HIV-1) RT). RTs lack proofreading activity and display high intrinsic error rates. Besides, high recombination rates observed in retroviruses are promoted by poor processivity that causes template switching, a hallmark of reverse transcription. HIV-1 RT inhibitors acting on its polymerase activity constitute the backbone of current antiretroviral therapies, although novel drugs, including ribonuclease H inhibitors, are still necessary to fight HIV infections. In Hepadnaviridae and Caulimoviridae, reverse transcription leads to the formation of nicked circular DNAs that will be converted into episomal DNA in the host cell nucleus. Structural and biochemical information on their polymerases is limited, although several drugs inhibiting HIV-1 RT are known to be effective against the human hepatitis B virus polymerase. In this review, we summarize current knowledge on reverse transcription in the five virus families and discuss available biochemical and structural information on RTs, including their biosynthesis, enzymatic activities, and potential inhibition.

Fidelity of classwide-resistant HIV-2 reverse transcriptase and differential contribution of K65R to the accuracy of HIV-1 and HIV-2 reverse transcriptases

Nucleoside reverse transcriptase (RT) inhibitors constitute the backbone of current therapies against human immunodeficiency virus type 1 and type 2 (HIV-1 and HIV-2, respectively). However, mutational pathways leading to the development of nucleoside analogue resistance are different in both types of HIV. In HIV-2, resistance to all approved nucleoside analogues is conferred by the combination of RT substitutions K65R, Q151M and M184V. Nucleotide incorporation kinetic analyses of mutant and wild-type (WT) HIV-2 RTs show that the triple-mutant has decreased catalytic efficiency due to the presence of M184V. Although similar effects were previously reported for equivalent mutations in HIV-1 RT, the HIV-2 enzymes were catalytically less efficient. Interestingly, in highly divergent HIV-1 RTs, K65R confers several-fold increased accuracy of DNA synthesis. We have determined the intrinsic fidelity of DNA synthesis of WT HIV-2 RT and mutants K65R and K65R/Q151M/M184V. Our results show that those changes in HIV-2 RT have a relatively small impact on nucleotide selectivity. Furthermore, we found that there were less than two-fold differences in error rates obtained with forward mutation assays using mutant and WT HIV-2 RTs. A different conformation of the β3-β4 hairpin loop in HIV-1 and HIV-2 RTs could probably explain the differential effects of K65R.

Preparation and characterization of the RNase H domain of Moloney murine leukemia virus reverse transcriptase

Moloney murine leukemia virus reverse transcriptase (MMLV RT) contains fingers, palm, thumb, and connection subdomains as well as an RNase H domain. The DNA polymerase active site resides in the palm subdomain, and the RNase H active site is located in the RNase H domain. The RNase H domain contains a positively charged α-helix called the C helix (H(594)GEIYRRR(601)), that is thought to be involved in substrate recognition. In this study, we expressed three versions of the RNase H domain in Escherichia coli, the wild-type domain (WT) (residues Ile498-Leu671) and two variants that lack the regions containing the C helix (Ile593-Leu603 and Gly595-Thr605, which we called ΔC1 and ΔC2, respectively) with a strep-tag at the N-terminus and a deca-histidine tag at the C-terminus. These peptides were purified from the cells by anion-exchange, Ni(2+) affinity, and Strep-Tactin affinity column chromatography, and then the tags were removed by proteolysis. In an RNase H assay using a 25-bp RNA-DNA heteroduplex, WT, ΔC1, and ΔC2 produced RNA fragments ranging from 7 to 16 nucleotides (nt) whereas the full-length MMLV RT (Thr24-Leu671) produced 14-20-nt RNA fragments, suggesting that elimination of the fingers, palm, thumb, and connection subdomains affects the binding of the RNase H domain to the RNA-DNA heteroduplex. The activity levels of WT, ΔC1, and ΔC2 were estimated to be 1%, 0.01%, and 0.01% of full-length MMLV RT activity, indicating that the C helix is important, but not critical, for the activity of the isolated RNase H domain.

Major groove binding track residues of the connection subdomain of human immunodeficiency virus type 1 reverse transcriptase enhance cDNA synthesis at high temperatures

At high temperatures, RNA denaturation can improve the efficiency and specificity of reverse transcription. Refined structures and molecular models of HIV-1 reverse transcriptases (RTs) from phylogenetically distant clades (i.e., group M subtype B and group O) revealed a major interaction between the template-primer and the Arg³⁵⁸-Gly³⁵⁹-Ala³⁶⁰ triad in the large subunit of HIV-1M/B RT. However, fewer contacts were predicted for the equivalent Lys³⁵⁸-Ala³⁵⁹-Ser³⁶⁰ triad of HIV-1O RT and the nucleic acid. An engineered HIV-1O K358R/A359G/S360A RT showed increased cDNA synthesis efficiency above 68 °C, as determined by qualitative and quantitative reverse transcription polymerase chain reactions. In comparison with wild-type HIV-1O RT, the mutant enzyme showed higher thermal stability but retained wild-type RNase H activity. Mutations that increased the accuracy of HIV-1M/B RTs were tested in combination with the K358R/A359G/S360A triple mutation. Some of them (e.g., F61A, K65R, K65R/V75I, and V148I) had a negative effect on reverse transcription efficiency above 65 °C. RTs with improved DNA binding affinities also showed higher cDNA synthesis efficiencies at elevated temperatures. Two of the most thermostable RTs (i.e., mutants T69SSG/K358R/A359G/S360A and K358R/A359G/S360A/E478Q) showed moderately increased fidelity in forward mutation assays. Our results demonstrate that the triad of Arg³⁵⁸, Gly³⁵⁹, and Ala³⁶⁰ in the major groove binding track of HIV-1 RT is a major target for RT stabilization, and most relevant for improving reverse transcription efficiency at high temperatures.

Temperature effects on the fidelity of a thermostable HIV-1 reverse transcriptase

Transcriptomics and gene expression analysis are largely dependent of the availability of efficient thermostable reverse transcriptases (RTs). However, the intrinsic fidelity of DNA synthesis catalyzed by retroviral RTs is low. Reported error rates are in the range 1.2 × 10(-5)-6.7 × 10(-4), with oncoretroviral RTs being the most faithful enzymes. Wild-type HIV-1 group O (HIV-1O) RT is a thermostable polymerase that is able to synthesize cDNA at temperatures as high as 70 °C. At 37 °C, its error rate has been estimated at 5.8 × 10(-5) in M13mp2 lacZ-based forward mutation assays. However, at higher temperatures (e.g. 50 and 55 °C), the accuracy of HIV-1O RT is increased by approximately two- to five-fold. At 55 °C, the HIV-1O RT error rate (1.3 × 10(-5)) was similar to that shown by the AffinityScript (Agilent Technologies Inc., La Jolla, CA, USA) RT, a commercially available thermostable murine leukaemia virus RT. At higher temperatures, the increased accuracy of the HIV-1 enzyme results from a lower base substitution error rate, although it shows a higher tendency to introduce frameshifts. Kinetic studies carried out with model template-primers suggest minor differences in nucleotide discrimination, although, at higher temperatures, HIV-1O RT showed a reduced ability to extend mispaired template-primers.