Fidelity of two retroviral reverse transcriptases during DNA-dependent DNA synthesis in vitro

DNA polymerase activity on synthetic N3'→P5' phosphoramidate DNA templates

Genetic polymers that could plausibly govern life in the universe might inhabit a broad swath of chemical space. A subset of these genetic systems can exchange information with RNA and DNA and could therefore form the basis for model protocells in the laboratory. N3'→P5' phosphoramidate (NP) DNA is defined by a conservative linkage substitution and has shown promise as a protocellular genetic material, but much remains unknown about its functionality and fidelity due to limited enzymatic tools. Conveniently, we find widespread NP-DNA-dependent DNA polymerase activity among reverse transcriptases, an observation consistent with structural studies of the RNA-like conformation of NP-DNA duplexes. Here, we analyze the consequences of this unnatural template linkage on the kinetics and fidelity of DNA polymerization activity catalyzed by wild-type and variant reverse transcriptases. Template-associated deficits in kinetics and fidelity suggest that even highly conservative template modifications give rise to error-prone DNA polymerase activity. Enzymatic copying of NP-DNA sequences is nevertheless an important step toward the future study and engineering of this synthetic genetic polymer.

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.

Deoxyuridine in DNA has an inhibitory and promutagenic effect on RNA transcription by diverse RNA polymerases

dUTP is a close structural congener of dTTP and can be readily incorporated into DNA opposite to adenine during DNA replication leading to non-mutagenic dU/A base pairs ('uracilation'). We find that dU/A pairs located within DNA transcriptional templates optimized for either T7 RNA polymerase (T7 RNAP) or human RNA polymerase II (pol II) have inhibitory and mutagenic effects on transcription. The data for T7 RNAP establishes that even a single dU/A pair can inhibit promoter binding and transcription initiation up to 30-fold, and that inhibitory effects on transcription elongation are also possible. Sequencing of the mRNA transcribed from uniformly uracilated DNA templates by T7 RNAP indicated an increased frequency of transversion and insertion mutations compared to all T/A templates. Strong effects of dU/A pairs on cellular transcription activity and fidelity were also observed with RNA pol II using uracil base excision repair (UBER)-deficient human cells. At the highest levels of template uracilation, transcription by RNA pol II was completely blocked. We propose that these effects arise from the decreased thermodynamic stability and increased dynamics of dU/A pairs in DNA. The potential implications of these findings on gene regulation and disease are discussed.

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.

HIV Reverse Transcriptase and Protease Genes Variability Can Be a Biomarker Associated with HIV and Hepatitis B or C Coinfection

Variability of the HIV reverse transcriptase (RT) and protease (PR) genes has been used as indicators of drug resistance and as a mean to evaluate phylogenetic relationships among circulating virus. However, these studies have been carried in HIV mono-infected populations. The goal of this study was to evaluate, for the first time, the HIV PR and RT sequences from HIV/HBV and HIV/HCV co-infected patients. HIV PR and RT genes were amplificated and sequenced to resistance analysis. The bioinformatics analysis was performed to infer about sequences clustering and molecular evolution. The results showed that the most frequent amino acid substitutions in RT were L214F (67.6%), I135T (55.9%), and in PR was V15I (41.2%). The molecular clock analysis showed that the HIV circulating in co-infected patients were separated in two clusters in the years 1999-2000. Some patients included as HIV mono-infected according patients' medical records and inside the co-infected cluster were, in fact, co-infected by PCR analysis. Analysis of the decision trees showed susceptibility to lamivudine and emtricitabine were important attribute to characterize co-infected patients. In conclusion, the results obtained in this study suggest, for the first time, that HIV RT and PR genes variability could be a genetic biomarker to coinfection.

Lethal mutagenesis of an RNA plant virus via lethal defection

Lethal mutagenesis is an antiviral therapy that relies on increasing the viral mutation rate with mutagenic nucleoside or base analogues. Currently, the molecular mechanisms that lead to virus extinction through enhanced mutagenesis are not fully understood. Increasing experimental evidence supports the lethal defection model of lethal mutagenesis of RNA viruses, where replication-competent-defectors drive infective virus towards extinction. Here, we address lethal mutagenesis in vivo using 5-fluorouracil (5-FU) during the establishment of tobacco mosaic virus (TMV) systemic infections in N. tabacum. The results show that 5-FU decreased the infectivity of TMV without affecting its viral load. Analysis of molecular clones spanning two genomic regions showed an increase of the FU-related base transitions A → G and U → C. Although the mutation frequency or the number of mutations per molecule did not increase, the complexity of the mutant spectra and the distribution of the mutations were altered. Overall, our results suggest that 5-FU antiviral effect on TMV is associated with the perturbation of the mutation-selection balance in the genomic region of the RNA-dependent RNA polymerase (RdRp). Our work supports the lethal defection model for lethal mutagenesis in vivo in a plant RNA virus and opens the way to study lethal mutagens in plant-virus systems.

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.

A multiplex RNA quantification method to determine the absolute amounts of mRNA without reverse transcription

We have developed a highly sensitive microarray-based method that determines the absolute amounts of mRNA in a total RNA sample in a multiplex manner without reverse transcription. This direct mRNA measurement promotes high-throughput testing and reduces bias in transcriptome analyses. Furthermore, quantification of the absolute amount of mRNA allows transcriptome analysis without common controls or additional, complicated normalization. The method, called Photo-DEAN, was validated using chemically synthesized RNAs of known quantities and mouse liver total RNA samples. We found that the absolute amounts of mRNA were successfully measured without the cDNA synthesis step, with a sensitivity of 15 zmol achieved in 7 h.

Pausing kinetics dominates strand-displacement polymerization by reverse transcriptase

Reverse transcriptase (RT) catalyzes the conversion of the viral RNA into an integration-competent double-stranded DNA, with a variety of enzymatic activities that include the ability to displace a non-template strand concomitantly with polymerization. Here, using high-resolution optical tweezers to follow the activity of the murine leukemia Virus RT, we show that strand-displacement polymerization is frequently interrupted. Abundant pauses are modulated by the strength of the DNA duplex ∼8 bp ahead, indicating the existence of uncharacterized RT/DNA interactions, and correspond to backtracking of the enzyme, whose recovery is also modulated by the duplex strength. Dissociation and reinitiation events, which induce long periods of inactivity and are likely the rate-limiting step in the synthesis of the genome in vivo, are modulated by the template structure and the viral nucleocapsid protein. Our results emphasize the potential regulatory role of conserved structural motifs, and may provide useful information for the development of potent and specific inhibitors.

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.

RNA editing generates cellular subsets with diverse sequence within populations

RNA editing is a mutational mechanism that specifically alters the nucleotide content in transcribed RNA. However, editing rates vary widely, and could result from equivalent editing amongst individual cells, or represent an average of variable editing within a population. Here we present a hierarchical Bayesian model that quantifies the variance of editing rates at specific sites using RNA-seq data from both single cells, and a cognate bulk sample to distinguish between these two possibilities. The model predicts high variance for specific edited sites in murine macrophages and dendritic cells, findings that we validated experimentally by using targeted amplification of specific editable transcripts from single cells. The model also predicts changes in variance in editing rates for specific sites in dendritic cells during the course of LPS stimulation. Our data demonstrate substantial variance in editing signatures amongst single cells, supporting the notion that RNA editing generates diversity within cellular populations.

RNA editing rate detected from bulk RNA-seq data can vary widely. Here, by constructing a hierarchical Bayesian model, the authors report substantial variance in editing signatures detected by RNA-seq data from both single cells and a cognate bulk sample.

Anti-HIV cytotoxicity enzyme inhibition and molecular docking studies of quinoline based chalcones as potential non-nucleoside reverse transcriptase inhibitors (NNRT)

A series of fourteen (A1 - A14) qunioline based chalcones were screened for reverse transcriptase inhibitors (RT) and found potentially active against RT. Bioassay, theoretical and dockings studies with RT (the enzyme required for reverse transcription of viral RNA) results showed that the type and positions of the substituents seemed to be critical for their inhibition against RT. The bromo and chloro substituted chalcone displayed high degree of inhibition against RT. The A4 andA6 showed high interaction with RT, contributing high free binding energy (ΔG -9.30 and -9.13kcal) and RT inhibition value (IC50 0.10μg/ml and 0.11μg/ml).

Triticum mosaic virus exhibits limited population variation yet shows evidence of parallel evolution after replicated serial passage in wheat

An infectious cDNA clone of Triticum mosaic virus (TriMV) (genus Poacevirus; family Potyviridae) was used to establish three independent lineages in wheat to examine intra-host population diversity levels within protein 1 (P1) and coat protein (CP) cistrons over time. Genetic variation was assessed at passages 9, 18 and 24 by single-strand conformation polymorphism, followed by nucleotide sequencing. The founding P1 region genotype was retained at high frequencies in most lineage/passage populations, while the founding CP genotype disappeared after passage 18 in two lineages. We found that rare TriMV genotypes were present only transiently and lineages followed independent evolutionary trajectories, suggesting that genetic drift dominates TriMV evolution. These results further suggest that experimental populations of TriMV exhibit lower mutant frequencies than that of Wheat streak mosaic virus (genus Tritimovirus; family Potyviridae) in wheat. Nevertheless, there was evidence for parallel evolution at a synonymous site in the TriMV CP cistron.

Current perspectives on HIV-1 antiretroviral drug resistance

Current advancements in antiretroviral therapy (ART) have turned HIV-1 infection into a chronic and manageable disease. However, treatment is only effective until HIV-1 develops resistance against the administered drugs. The most recent antiretroviral drugs have become superior at delaying the evolution of acquired drug resistance. In this review, the viral fitness and its correlation to HIV-1 mutation rates and drug resistance are discussed while emphasizing the concept of lethal mutagenesis as an alternative therapy. The development of resistance to the different classes of approved drugs and the importance of monitoring antiretroviral drug resistance are also summarized briefly.

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.

Experimental transmission of enzootic nasal adenocarcinoma in sheep

Enzootic nasal adenocarcinoma (ENA) is a contagious neoplasm of the secretory epithelial cells of the nasal mucosa of sheep and goats. It is associated with the betaretrovirus, enzootic nasal tumor virus (ENTV), but a causative relationship has yet to be demonstrated. In this study, 14-day-old lambs were experimentally infected via nebulization with cell-free tumor filtrates derived from naturally occurring cases of ENA. At 12 weeks post-infection (wpi), one of the five infected lambs developed clinical signs, including continuous nasal discharge and open mouth breathing, and was euthanized. Necropsy revealed the presence of a large bilateral tumor occupying the nasal cavity. At 45 wpi, when the study was terminated, none of the remaining infected sheep showed evidence of tumors either by computed tomography or post-mortem examination. ENTV-1 proviral DNA was detected in the nose, lung, spleen, liver and kidney of the animal with experimentally induced ENA, however there was no evidence of viral protein expression in tissues other than the nose. Density gradient analysis of virus particles purified from the experimentally induced nasal tumor revealed a peak reverse transcriptase (RT) activity at a buoyant density of 1.22 g/mL which was higher than the 1.18 g/mL density of peak RT activity of virus purified from naturally induced ENA. While the 1.22 g/mL fraction contained primarily immature unprocessed virus particles, mature virus particles with a similar morphology to naturally occurring ENA could be identified by electron microscopy. Full-length sequence analysis of the ENTV-1 genome from the experimentally induced tumor revealed very few nucleotide changes relative to the original inoculum with only one conservative amino acid change. Taken together, these results demonstrate that ENTV-1 is associated with transmissible ENA in sheep and that under experimental conditions, lethal tumors are capable of developing in as little as 12 wpi demonstrating the acutely oncogenic nature of this ovine betaretrovirus.

Altered error specificity of RNase H-deficient HIV-1 reverse transcriptases during DNA-dependent DNA synthesis

Asp(443) and Glu(478) are essential active site residues in the RNase H domain of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT). We have investigated the effects of substituting Asn for Asp(443) or Gln for Glu(478) on the fidelity of DNA-dependent DNA synthesis of phylogenetically diverse HIV-1 RTs. In M13mp2 lacZα-based forward mutation assays, HIV-1 group M (BH10) and group O RTs bearing substitutions D443N, E478Q, V75I/D443N or V75I/E478Q showed 2.0- to 6.6-fold increased accuracy in comparison with the corresponding wild-type enzymes. This was a consequence of their lower base substitution error rates. One-nucleotide deletions and insertions represented between 30 and 68% of all errors identified in the mutational spectra of RNase H-deficient HIV-1 group O RTs. In comparison with the wild-type RT, these enzymes showed higher frameshift error rates and higher dissociation rate constants (koff) for DNA/DNA template-primers. The effects on frameshift fidelity were similar to those reported for mutation E89G and suggest that in HIV-1 group O RT, RNase H inactivation could affect template/primer slippage. Our results support a role for the RNase H domain during plus-strand DNA polymerization and suggest that mutations affecting RNase H function could also contribute to retrovirus variability during the later steps of reverse transcription.

Conserved hydrogen bonds and water molecules in MDR HIV-1 protease substrate complexes

The success of highly active antiretroviral therapy (HAART) in anti-HIV therapy is severely compromised by the rapidly developing drug resistance. HIV-1 protease inhibitors, part of HAART, are losing their potency and efficacy in inhibiting the target. Multi-drug resistant (MDR) 769 HIV-1 protease (resistant mutations at residues 10, 36, 46, 54, 62, 63, 71, 82, 84, 90) was selected for the present study to understand the binding to its natural substrates. The nine crystal structures of MDR769 HIV-1 protease substrate hepta-peptide complexes were analyzed in order to reveal the conserved structural elements for the purpose of drug design against MDR HIV-1 protease. Our structural studies demonstrated that highly conserved hydrogen bonds between the protease and substrate peptides, together with the conserved crystallographic water molecules, played a crucial role in the substrate recognition, substrate stabilization and protease stabilization. In addition, the absence of the key flap-ligand bridging water molecule might imply a different catalytic mechanism of MDR769 HIV-1 protease compared to that of wild type (WT) HIV-1 protease.

Application of metatranscriptomics to soil environments

The activities of soil microbial communities are of critical importance to terrestrial ecosystem functioning. The mechanisms that determine the interactions between soil microorganisms, their environment and neighbouring organisms, however, are poorly understood. Due to advances in sequencing technologies, an increasing number of metagenomics studies are being conducted on samples from diverse environments including soils. This information has not only increased our awareness of the functional potential of soil microbial communities, but also constitutes powerful reference material for soil metatranscriptomics studies. Metatranscriptomics provides a snapshot of transcriptional profiles that correspond to discrete populations within a microbial community at the time of sampling. This information can indicate the potential activities of complex microbial communities and the mechanisms that regulate them. Here we summarise the technical challenges for metatranscriptomics applied to soil environments and discuss approaches for gaining biologically meaningful insight into these datasets.

Thermostable DNA polymerase from a viral metagenome is a potent RT-PCR enzyme

Viral metagenomic libraries are a promising but previously untapped source of new reagent enzymes. Deep sequencing and functional screening of viral metagenomic DNA from a near-boiling thermal pool identified clones expressing thermostable DNA polymerase (Pol) activity. Among these, 3173 Pol demonstrated both high thermostability and innate reverse transcriptase (RT) activity. We describe the biochemistry of 3173 Pol and report its use in single-enzyme reverse transcription PCR (RT-PCR). Wild-type 3173 Pol contains a proofreading 3′-5′ exonuclease domain that confers high fidelity in PCR. An easier-to-use exonuclease-deficient derivative was incorporated into a PyroScript RT-PCR master mix and compared to one-enzyme (Tth) and two-enzyme (MMLV RT/Taq) RT-PCR systems for quantitative detection of MS2 RNA, influenza A RNA, and mRNA targets. Specificity and sensitivity of 3173 Pol-based RT-PCR were higher than Tth Pol and comparable to three common two-enzyme systems. The performance and simplified set-up make this enzyme a potential alternative for research and molecular diagnostics.

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

Although recent reports have provided strong evidence to suggest that xenotropic murine leukemia virus-related virus (XMRV) is unlikely to be the causative agent of prostate cancer and chronic fatigue syndrome, this recombinant retrovirus can nonetheless infect human cells in vitro and induce a chronic infection in macaques. In the present study, we determined the accuracy of DNA synthesis of the reverse transcriptases (RTs) of XMRV and Moloney murine leukemia virus (MoMLV) using a combination of pre-steady-state kinetics of nucleotide incorporation and an M13mp2-based forward mutation assay. The results obtained were compared with those previously reported for the HIV type 1 BH10 strain (HIV-1(BH10)) RT. MoMLV and XMRV RTs were 13.9 and 110 times less efficient [as determined by the catalytic rate constant of the nucleotide incorporation reaction ((pol))/equilibrium constant (K(d))] than the HIV-1(BH10) RT in incorporating correct nucleotides. Misinsertion and mispair extension kinetic studies demonstrated that MoMLV RT was more accurate than the HIV-1(BH10) RT. In comparison with the MoMLV RT, the XMRV RT showed decreased mispair extension fidelity and was less faithful when misincorporating C or A opposite A. However, the XMRV RT showed stronger selectivity against G in misinsertion fidelity assays. Forward mutation assays revealed that XMRV and MoMLV RTs had similar accuracy of DNA-dependent DNA synthesis, but were > 13 times more faithful than the HIV-1(BH10) enzyme. The mutational spectra of XMRV and MoMLV RTs were similar in having a relatively higher proportion of frameshifts and transversions compared with the HIV-1(BH10) RT. However, the XMRV polymerase was less prone to introduce large deletions and one-nucleotide insertions.

SNP Discovery through Next-Generation Sequencing and Its Applications

The decreasing cost along with rapid progress in next-generation sequencing and related bioinformatics computing resources has facilitated large-scale discovery of SNPs in various model and nonmodel plant species. Large numbers and genome-wide availability of SNPs make them the marker of choice in partially or completely sequenced genomes. Although excellent reviews have been published on next-generation sequencing, its associated bioinformatics challenges, and the applications of SNPs in genetic studies, a comprehensive review connecting these three intertwined research areas is needed. This paper touches upon various aspects of SNP discovery, highlighting key points in availability and selection of appropriate sequencing platforms, bioinformatics pipelines, SNP filtering criteria, and applications of SNPs in genetic analyses. The use of next-generation sequencing methodologies in many non-model crops leading to discovery and implementation of SNPs in various genetic studies is discussed. Development and improvement of bioinformatics software that are open source and freely available have accelerated the SNP discovery while reducing the associated cost. Key considerations for SNP filtering and associated pipelines are discussed in specific topics. A list of commonly used software and their sources is compiled for easy access and reference.

Nine crystal structures determine the substrate envelope of the MDR HIV-1 protease

Under drug selection pressure, emerging mutations render HIV-1 protease drug resistant, leading to the therapy failure in anti-HIV treatment. It is known that nine substrate cleavage site peptides bind to wild type (WT) HIV-1 protease in a conserved pattern. However, how the multidrug-resistant (MDR) HIV-1 protease binds to the substrate cleavage site peptides is yet to be determined. MDR769 HIV-1 protease (resistant mutations at residues 10, 36, 46, 54, 62, 63, 71, 82, 84, and 90) was selected for present study to understand the binding to its natural substrates. MDR769 HIV-1 protease was co-crystallized with nine substrate cleavage site hepta-peptides. Crystallographic studies show that MDR769 HIV-1 protease has an expanded substrate envelope with wide open flaps. Furthermore, ligand binding energy calculations indicate weaker binding in MDR769 HIV-1 protease-substrate complexes. These results help in designing the next generation of HIV-1 protease inhibitors by targeting the MDR HIV-1 protease.

Thermostable HIV-1 group O reverse transcriptase variants with the same fidelity as murine leukaemia virus reverse transcriptase

Wild-type HIV-1 group O RT (reverse transcriptase) shows increased thermostability in comparison with HIV-1 group M subtype B RT and MLV (murine leukaemia virus) RT. However, its utility in the amplification of RNA targets is limited by the reduced accuracy of lentiviral RTs compared with oncoretroviral RTs (i.e. MLV RT). The effects of the mutations K65R, R78A and K65R/V75I on the fidelity of HIV-1 group O RTs were studied using gel-based and M13mp2 lacZ forward-mutation fidelity assays. Forward-mutation assays demonstrated that mutant RTs K65R, R78A and K65R/V75I showed >9-fold increased accuracy in comparison with the wild-type enzyme and were approximately two times more faithful than the MLV RT. Compared with MLV RT, all of the tested HIV-1 group O RT variants showed decreased frameshift fidelity. However, K65R RT showed a higher tendency to introduce one-nucleotide deletions in comparison with other HIV-1 group O RT variants. R78A had a destabilizing effect on the RT, either in the presence or absence of V75I. At temperatures above 52 °C, K65R and K65R/V75I retained similar levels of DNA polymerase activity to the wild-type HIV-1 group O RT, but were more efficient than HIV-1 group M subtype B and MLV RTs. K65R, K65R/V75I and R78A RTs showed decreased misinsertion and mispair extension fidelity in comparison with the wild-type enzyme for most base pairs studied. These assays revealed that nucleotide selection is mainly governed by kpol (pol is polymerization) in the case of K65R, whereas both kpol and Kd affect nucleotide discrimination in the case of K65R/V75I.

Single-molecule direct RNA sequencing without cDNA synthesis

Methods for in-depth genome-wide characterization of transcriptomes and quantification of transcript levels using various microarray and next-generation sequencing technologies have emerged as valuable tools for understanding cellular physiology and human disease biology and have begun to be utilized in various clinical diagnostic applications. Current methods, however, typically require RNA to be converted to complementary DNA prior to measurements. This step has been shown to introduce many biases and artifacts. In order to best characterize the 'true' transcriptome, the single-molecule direct RNA sequencing (DRS) technology was developed. This review focuses on the underlying principles behind the DRS, sample preparation steps, and the current and novel avenues of research and applications DRS offers.

The reverse transcriptase encoded by the non-LTR retrotransposon R2 is as error-prone as that encoded by HIV-1

Reverse transcriptases (RTs) encoded by a wide range of mobile retroelements have had a major impact on the structure and function of genomes. Among the most abundant elements in eukaryotes are the non long terminal repeat (LTR) retrotransposons. Here we compare the dNTP concentration requirements and error rates of the RT encoded by the non-LTR retrotransposon R2 of Bombyx mori with the well-characterized RTs of retroviruses. Surprisingly, R2 was found to have properties more similar to those of lentiviral RTs, such as human immunodeficiency virus type 1 (HIV-1), than to those of oncoretroviral RTs, such as murine leukemia virus. Like HIV-1 RT, R2 RT was able to synthesize DNA at low dNTP concentrations, suggesting that R2 is able to retrotranspose in nondividing cells. R2 RT also showed levels of misincorporation in biased dNTP pools and replication error rates in M13 lacZα forward mutation assays, similar to HIV-1 RT. Most of the R2 base substitutions in the forward mutation assay were caused by the misincorporation of dTMP. Analogous to HIV-1, the high error rate of R2 RT appears to be a result of its ability to extend mismatches once generated. We suggest that the low fidelity of R2 RT is a by-product of the flexibility of its active site/dNTP binding pocket required for the target-primed reverse transcription reaction used by R2 for retrotransposition. Finally, we discuss that in spite of the high R2 RT error rate, the long-term nucleotide substitution rate for R2 is not significantly above that associated with cellular DNA replication, based on the frequency of R2 retrotranspositions determined in natural populations.

RNA sequencing: advances, challenges and opportunities

In the few years since its initial application, massively parallel cDNA sequencing, or RNA-seq, has allowed many advances in the characterization and quantification of transcriptomes. Recently, several developments in RNA-seq methods have provided an even more complete characterization of RNA transcripts. These developments include improvements in transcription start site mapping, strand-specific measurements, gene fusion detection, small RNA characterization and detection of alternative splicing events. Ongoing developments promise further advances in the application of RNA-seq, particularly direct RNA sequencing and approaches that allow RNA quantification from very small amounts of cellular materials.

RNA sequencing: advances, challenges and opportunities

In the few years since its initial application, massively parallel cDNA sequencing, or RNA-seq, has allowed many advances in the characterization and quantification of transcriptomes. Recently, several developments in RNA-seq methods have provided an even more complete characterization of RNA transcripts. These developments include improvements in transcription start site mapping, strand-specific measurements, gene fusion detection, small RNA characterization and detection of alternative splicing events. Ongoing developments promise further advances in the application of RNA-seq, particularly direct RNA sequencing and approaches that allow RNA quantification from very small amounts of cellular materials.

Strand-specific transcriptome profiling with directly labeled RNA on genomic tiling microarrays

Background

With lower manufacturing cost, high spot density, and flexible probe design, genomic tiling microarrays are ideal for comprehensive transcriptome studies. Typically, transcriptome profiling using microarrays involves reverse transcription, which converts RNA to cDNA. The cDNA is then labeled and hybridized to the probes on the arrays, thus the RNA signals are detected indirectly. Reverse transcription is known to generate artifactual cDNA, in particular the synthesis of second-strand cDNA, leading to false discovery of antisense RNA. To address this issue, we have developed an effective method using RNA that is directly labeled, thus by-passing the cDNA generation. This paper describes this method and its application to the mapping of transcriptome profiles.

Results

RNA extracted from laboratory cultures of Porphyromonas gingivalis was fluorescently labeled with an alkylation reagent and hybridized directly to probes on genomic tiling microarrays specifically designed for this periodontal pathogen. The generated transcriptome profile was strand-specific and produced signals close to background level in most antisense regions of the genome. In contrast, high levels of signal were detected in the antisense regions when the hybridization was done with cDNA. Five antisense areas were tested with independent strand-specific RT-PCR and none to negligible amplification was detected, indicating that the strong antisense cDNA signals were experimental artifacts.

Conclusions

An efficient method was developed for mapping transcriptome profiles specific to both coding strands of a bacterial genome. This method chemically labels and uses extracted RNA directly in microarray hybridization. The generated transcriptome profile was free of cDNA artifactual signals. In addition, this method requires fewer processing steps and is potentially more sensitive in detecting small amount of RNA compared to conventional end-labeling methods due to the incorporation of more fluorescent molecules per RNA fragment.

Transcriptome profiling using single-molecule direct RNA sequencing

Methods for in-depth characterization of transcriptomes and quantification of transcript levels have emerged as valuable tools for understanding cellular physiology and human disease biology, and have begun to be utilized in various clinical diagnostic applications. Today, current methods utilized by the scientific community typically require RNA to be converted to cDNA prior to comprehensive measurements. However, this cDNA conversion process has been shown to introduce many biases and artifacts that interfere with the proper characterization and quantitation of transcripts. We have developed a direct RNA sequencing (DRS) approach, in which, unlike other technologies, RNA is sequenced directly without prior conversion to cDNA. The benefits of DRS include the ability to use minute quantities (e.g. on the order of several femtomoles) of RNA with minimal sample preparation, the ability to analyze short RNAs which pose unique challenges for analysis using cDNA-based approaches, and the ability to perform these analyses in a low-cost and high-throughput manner. Here, we describe the strategies and procedures we employ to prepare various RNA species for analysis with DRS.

Intra-host diversities of the receptor-binding domain of stork faeces-derived avian H5N1 viruses and its significance as predicted by molecular dynamic simulation

Virus evolution facilitates the emergence of viruses with unpredictable impacts on human health. This study investigated intra-host variations of the receptor-binding domain (RBD) of the haemagglutinin (HA) gene of the avian H5N1 viruses obtained from the 2004 and 2005 epidemics. The results showed that the mutation frequency of the RBD ranged from 0.3 to 0.6 %. The mutations generated one consensus and several minor populations. The consensus population of the 2004 epidemic was transmitted to the 2005 outbreak with increased frequency (39 and 45 %, respectively). Molecular dynamics simulation was applied to predict the significance of the variants. The results revealed that the consensus sequence (E218K/V248I) interacted unstably with sialic acid (SA) with an α2,6 linkage (SAα2,6Gal). Although the mutated K140R/E218K/V248I and Y191C/E218K/V248I sequences decreased the HA binding capacity to α2,3-linked SA, they were shown to bind α2,6-linked SA with increased affinity. Moreover, the substitutions at aa 140 and 191 were positive-selection sites. These data suggest that the K140R and Y191C mutations may represent a step towards human adaptation of the avian H5N1 virus.

Mutation rates and intrinsic fidelity of retroviral reverse transcriptases

Retroviruses are RNA viruses that replicate through a DNA intermediate, in a process catalyzed by the viral reverse transcriptase (RT). Although cellular polymerases and host factors contribute to retroviral mutagenesis, the RT errors play a major role in retroviral mutation. RT mutations that affect the accuracy of the viral polymerase have been identified by in vitro analysis of the fidelity of DNA synthesis, by using enzymological (gel-based) and genetic assays (e.g., M13mp2 lacZ forward mutation assays). For several amino acid substitutions, these observations have been confirmed in cell culture using viral vectors. This review provides an update on studies leading to the identification of the major components of the fidelity center in retroviral RTs.

Unequal distribution of RT-PCR artifacts along the E1-E2 region of Hepatitis C virus

Although viral variability studies have focused traditionally on consensus sequences, the relevance of molecular clone sequences for studying viral evolution at the intra-host level is being increasingly recognized. However, for this approach to be reliable, RT-PCR artifacts do not have to contribute excessively to the observed variability. Molecular clone sequences were obtained from an in vitro transcript to estimate the maximum error rate associated to RT-PCR for the Hepatitis C virus (HCV) E1-E2 region. On average, the frequency of RT-PCR errors was one order of magnitude lower than the level of intra-host genetic variability observed in samples from an HCV outbreak. However, RT-PCR errors were not distributed evenly along the E1-E2 region and were concentrated heavily in the hypervariable region 2 (HVR 2). Although it is concluded that RT-PCR molecular clone sequences are reliable, these results warn against extrapolation of RT-PCR error rates to different genome regions. The data suggest that the RNA sequence context or secondary structure can determine the fidelity of in vitro transcription or reverse transcription. Potentially, these factors might also modify the fidelity of the viral polymerase.

Direct RNA sequencing

Our understanding of human biology and disease is ultimately dependent on a complete understanding of the genome and its functions. The recent application of microarray and sequencing technologies to transcriptomics has changed the simplistic view of transcriptomes to a more complicated view of genome-wide transcription where a large fraction of transcripts emanates from unannotated parts of genomes, and underlined our limited knowledge of the dynamic state of transcription. Most of this broad body of knowledge was obtained indirectly because current transcriptome analysis methods typically require RNA to be converted to complementary DNA (cDNA) before measurements, even though the cDNA synthesis step introduces multiple biases and artefacts that interfere with both the proper characterization and quantification of transcripts. Furthermore, cDNA synthesis is not particularly suitable for the analysis of short, degraded and/or small quantity RNA samples. Here we report direct single molecule RNA sequencing without prior conversion of RNA to cDNA. We applied this technology to sequence femtomole quantities of poly(A)(+) Saccharomyces cerevisiae RNA using a surface coated with poly(dT) oligonucleotides to capture the RNAs at their natural poly(A) tails and initiate sequencing by synthesis. We observed transcript 3' end heterogeneity and polyadenylated small nucleolar RNAs. This study provides a path to high-throughput and low-cost direct RNA sequencing and achieving the ultimate goal of a comprehensive and bias-free understanding of transcriptomes.

RNA-dependent DNA polymerases

Reverse transcriptases (RTs) are multifunctional enzymes, but are mainly used as RNA-directed DNA polymerases in first-strand cDNA synthesis. Specifically, oligodeoxynucleotides are used as primers for extension on RNA templates. The DNA synthesized from an RNA template is referred to as complementary DNA (cDNA) and is often used as a template for PCR or converted to dsDNA for cloning. This unit describes appropriate reaction conditions for RTs from Moloney murine leukemia virus (MMLV) and avian myeloblastosis virus (AMV), along with applications such as synthesizing cDNA, 3' fill-in reactions, and labeling the 3' terminus of DNA fragments with 5' protruding ends, and DNA sequencing.

Quantitative analysis of clinically relevant mutations occurring in lymphoid cells harboring gamma-retrovirus-encoded hsvtk suicide genes

The in vivo regulation of T lymphocyte activity by the activation of a suicide mechanism is an essential paradigm for the safety of adoptive cell therapies. In light of reports showing that gamma-retroviral vector-encoded herpes simplex virus thymidine kinase (hsvtk) undergoes recombination, we undertook a thorough investigation of the genomic stability of SFG-based vectors using two variants of the wild-type hsvtk gene. In a large panel of independent clones, we demonstrate that both hsvtk genes undergo recombination with molecular signatures indicative of template switching in GC-rich regions displaying homology at the deletion junctions or RNA splicing. In the absence of ganciclovir selection, the frequency of recombination is 3% per retroviral replication cycle. Our results underscore the importance of the five nucleotide difference between the two hsvtk genes that account for the presence of recombinogenic hot spots in one variant and not the other, indicating that the probability of RNA splicing is influenced by minute nucleotide changes in sequences adjacent to the splice donor and acceptor sites. Furthermore, our mutational analysis in an unbiased panel of human lymphoid cells (that is, without immune or ganciclovir-mediated selective pressure) provides a robust in vitro assay to predict and quantify clinically relevant mutations in hsvtk suicide genes, which can be applied to studying and improving the stability of any transgene expressed in gamma-retroviral or lentiviral vectors.

New approaches to HIV protease inhibitor drug design II: testing the substrate envelope hypothesis to avoid drug resistance and discover robust inhibitors

PURPOSE OF REVIEW

Drug resistance results when the balance between the binding of inhibitors and the turnover of substrates is perturbed in favor of the substrates. Resistance is quite widespread to the HIV-1 protease inhibitors permitting the protease to process its 10 different substrates. This processing of the substrates permits the virus HIV-1 to mature and become infectious. The design of HIV-1 protease inhibitors that closely fit within the substrate-binding region is proposed to be a strategy to avoid drug resistance.

RECENT FINDINGS

Cocrystal structures of HIV-1 protease with its substrates define an overlapping substrate-binding region or substrate envelope. Novel HIV-1 protease inhibitors that were designed to fit within this substrate envelope were found to retain high binding affinity and have a flat binding profile against a panel of drug-resistant HIV-1 proteases.

SUMMARY

The avoidance of drug resistance needs to be considered in the initial design of inhibitors to quickly evolving targets such as HIV-1 protease. Using a detailed knowledge of substrate binding appears to be a promising strategy for achieving this goal to obtain robust HIV-1 protease inhibitors.

Retroviral reverse transcriptases (other than those of HIV-1 and murine leukemia virus): a comparison of their molecular and biochemical properties

This chapter reviews most of the biochemical data on reverse transcriptases (RTs) of retroviruses, other than those of HIV-1 and murine leukemia virus (MLV) that are covered in detail in other reviews of this special edition devoted to reverse transcriptases. The various RTs mentioned are grouped according to their retroviral origins and include the RTs of the alpharetroviruses, lentiviruses (both primate, other than HIV-1, and non-primate lentiviruses), betaretroviruses, deltaretroviruses and spumaretroviruses. For each RT group, the processing, molecular organization as well as the enzymatic activities and biochemical properties are described. Several RTs function as dimers, primarily as heterodimers, while the others are active as monomeric proteins. The comparisons between the diverse properties of the various RTs show the common traits that characterize the RTs from all retroviral subfamilies. In addition, the unique features of the specific RTs groups are also discussed.

Mapping of the ICOS binding surface of murine B7h using an unbiased, cellular library of B7h mutants created by cyclical packaging rescue

Functional studies of immunologically relevant molecules often involve time-consuming generation and cloning of gene mutations prior to introduction into mammalian cells. We describe here an alternative mutagenesis approach that relies solely on transfer of helper-free retroviral supernatants to rapidly create in virtually any cell line of interest a large cellular library that retrovirally expresses a defined number of independent point mutations in a gene of interest. Using this rapid non-cloning approach, we generated a 3T3 cellular library retrovirally expressing 2 x 10(5) mutants of the murine costimulatory B7h gene. Screening of this unbiased cellular library identified six residues of murine B7h that are critical for binding to the ICOS receptor. These residues are located on the same strands of human B7h that were identified by targeted mutagenesis [Chattopadhyay, K., Bhatia, S., Fiser, A., Almo, S.C., Nathenson, S.G. (2006). Structural basis of inducible costimulator ligand costimulatory function: determination of the cell surface oligomeric state and functional mapping of the receptor-binding site of the protein. J. Immunol. 177, 3920], indicating that the ICOS receptor-binding interface is similar in mouse and human B7h. Based on this proof-of-principle study, CPR-based mutagenesis is applicable to studies of gene function in a variety of mammalian cells.

Trans-lesion synthesis and RNaseH activity by reverse transcriptases on a true abasic RNA template

While much is known about abasic DNA, the biological impact of abasic RNA is largely unexplored. To test the mutagenic potential of this RNA lesion in the context of retroviruses, we synthesized a 31-mer oligoribonucleotide containing an abasic (rAS) site and used it as a template for studying DNA primer extension by HIV-1, avian myeloblastosis virus (AMV) and moloney murine leukemia virus (MMLV) reversed transcriptases (RT). We found that trans-lesion synthesis readily takes place with HIV-1 RT and to a lesser extent with AMV RT while MMLV RT aborts DNA synthesis. The preference of dNTP incorporation follows the order A∼G > C∼T and thus obeys to the ‘A-rule’. In the case of HIV-1 RT, we measured the kinetic data of dNTP incorporation and compared it to abasic DNA. We found that A-incorporation is only 2-fold slower relative to a matched (undamaged) RNA template while it is 7-fold slower in the case of DNA. Furthermore, there is less discrimination in incorporation between the four dNTPs in the case of abasic RNA compared to abasic DNA. These experiments clearly point to a higher promiscuity of lesion bypass on abasic RNA. Given their known higher chemical stability, such rAS sites can clearly contribute to (retro)viral evolution.

Improved repeat identification and masking in Dipterans

Repetitive sequences are a major constituent of many eukaryote genomes and play roles in gene regulation, chromosome inheritance, nuclear architecture, and genome stability. The identification of repetitive elements has traditionally relied on in-depth, manual curation and computational determination of close relatives based on DNA identity. However, the rapid divergence of repetitive sequence has made identification of repeats by DNA identity difficult even in closely related species. Hence, the presence of unidentified repeats in genome sequences affects the quality of gene annotations and annotation-dependent analyses (e.g. microarray analyses). We have developed an enhanced repeat identification pipeline using two approaches. First, the de novo repeat finding program PILER-DF was used to identify interspersed repetitive elements in several recently finished Dipteran genomes. Repeats were classified, when possible, according to their similarity to known elements described in Repbase and GenBank, and also screened against annotated genes as one means of eliminating false positives. Second, we used a new program called RepeatRunner, which integrates results from both RepeatMasker nucleotide searches and protein searches using BLASTX. Using RepeatRunner with PILER-DF predictions, we masked repeats in thirteen Dipteran genomes and conclude that combining PILER-DF and RepeatRunner greatly enhances repeat identification in both well-characterized and un-annotated genomes.

Cis-preferential LINE-1 reverse transcriptase activity in ribonucleoprotein particles

LINE-1 retrotransposons (L1s) constitute approximately 17% of human DNA, and their activity continues to affect genome evolution. Retrotransposition-competent human L1s encode two proteins required for their mobility (ORF1p and ORF2p); however, biochemical activities associated with ORF2p have been difficult to detect in cells. Here, we show for the first time the colocalization of L1 RNA, ORF1p and ORF2p to a putative ribonucleoprotein retrotransposition intermediate. We further demonstrate that ORF2p preferentially uses its encoding RNA as a template for reverse transcription. Thus, our data provide the first biochemical evidence supporting the cis-preferential action of the L1 reverse transcriptase.

[RNA viruses and mutations]

Actively replicating RNA viruses in nature are continually changing their genetic information by spontaneous mutations. These changes often result in alterations in immune-sensitivity, drug-sensitivity, cell-tropism, and host-range, causing uncontrollability of the pathogen and emerging/re-emerging infections. To better understand the virus changes and develop effective methods to control the moving targets, it is essential to obtain information on changes in viral genomes and proteins. Although information on genetic changes is being accumulated very rapidly, assessment of changes in protein structure and function still requires time-consuming works. In this review, we will overview mutation studies of human immunodeficiency virus and other RNA viruses. In addition, we will introduce recent advances in the computational science and its application on mutation studies and drug development.

Comparison of DNA polymerase activities between recombinant feline immunodeficiency and leukemia virus reverse transcriptases

In this study, we present enzymatic differences found between recombinant RTs derived from feline leukemia virus and feline immunodeficiency virus. Firstly, FIV RT showed low steady state K(m) values for dNTPs compared to FeLV RT. Consistent with this, FIV RT synthesized DNA more efficiently than FeLV RT at low dNTP concentrations. We observed similar concentration-dependent activity differences between other lentiviral (HIV-1 and SIV) and non-lentiviral (MuLV and AMV) RTs. Second, FeLV RT showed less efficient misincorporation with biased dNTP pools and mismatch primer extension capabilities, compared to FIV RT. In M13mp2 lacZalpha forward mutation assays, FeLV RT displayed approximately 11-fold higher fidelity than FIV RT. Finally, FeLV RT was less sensitive to 3TCTP and ddATP than FIV RT. This study represents the comprehensive enzymatic characterization of RTs from a lentivirus and a non-lentivirus retrovirus from the same host species. The data presented here support enzymatic divergences seen among retroviral RTs.