From RNA to quasispecies: a DNA polymerase with proofreading activity is highly recommended for accurate assessment of viral diversity

Quick and Easy Assembly of a One-Step qRT-PCR Kit for COVID-19 Diagnostics Using In-House Enzymes

One-step reverse-transcription quantitative polymerase chain reaction (qRT-PCR) is the most widely applied method for COVID-19 diagnostics. Notwithstanding the facts that one-step qRT-PCR is well suited for the diagnosis of COVID-19 and that there are many commercially available one-step qRT-PCR kits in the market, their high cost and unavailability due to airport closures and shipment restriction became a major bottleneck that had driven the desire to produce the key components of such kits locally. Here, we provide a simple, economical, and powerful one-step qRT-PCR kit based on patent-free, specifically tailored versions of Moloney murine leukemia virus reverse transcriptase and Thermus aquaticus DNA polymerase and termed R3T (Rapid Research Response Team) one-step qRT-PCR. We also demonstrate the robustness of our enzyme production strategies and provide the optimal reaction conditions for their efficient augmentation in a one-step approach. Our kit was routinely able to reliably detect as low as 10 copies of the synthetic RNAs of SARS-CoV-2. More importantly, our kit successfully detected COVID-19 in clinical samples of broad viral titers with similar reliability and selectivity to that of the Invitrogen SuperScript III Platinum One-step qRT-PCR and TaqPath one-step RT-qPCR kits. Overall, our kit has shown robust performance in both laboratory settings and the Saudi Ministry of Health-approved testing facility.

Tick-Borne Encephalitis Virus Adaptation in Different Host Environments and Existence of Quasispecies

A highly virulent strain (Hypr) of tick-borne encephalitis virus (TBEV) was serially subcultured in the mammalian porcine kidney stable (PS) and Ixodes ricinus tick (IRE/CTVM19) cell lines, producing three viral variants. These variants exhibited distinct plaque sizes and virulence in a mouse model. Comparing the full-genome sequences of all variants, several nucleotide changes were identified in different genomic regions. Furthermore, different sequential variants were revealed to co-exist within one sample as quasispecies. Interestingly, the above-mentioned nucleotide changes found within the whole genome sequences of the new variants were present alongside the nucleotide sequence of the parental strain, which was represented as a minority quasispecies. These observations further imply that TBEV exists as a heterogeneous population that contains virus variants pre-adapted to reproduction in different environments, probably enabling virus survival in ticks and mammals.

PCR Amplification Strategies Towards Full-length HIV-1 Genome Sequencing

The advent of next-generation sequencing has enabled greater resolution of viral diversity and improved feasibility of full viral genome sequencing allowing routine HIV-1 full genome sequencing in both research and diagnostic settings. Regardless of the sequencing platform selected, successful PCR amplification of the HIV-1 genome is essential for sequencing template preparation. As such, full HIV-1 genome amplification is a crucial step in dictating the successful and reliable sequencing downstream. Here we reviewed existing PCR protocols leading to HIV-1 full genome sequencing. In addition to the discussion on basic considerations on relevant PCR design, the advantages as well as the pitfalls of the published protocols were reviewed.

Polymerase Mechanism-Based Method of Viral Attenuation

Vaccines remain the most effective way of preventing infection and spread of infectious diseases. These prophylactics have been used for centuries but still to this day only three main design strategies exist: (1) live attenuated virus (LAV) vaccines, (2) killed or inactivated virus vaccines, (3) and subunit vaccines of the three, the most efficacious vaccines remain LAVs. LAVs replicate in relevant tissues, elicit strong cellular and humoral responses, and often confer lifelong immunity. While this vaccine strategy has produced the majority of successful vaccines in use today, there are also important safety concerns to consider with this approach. In the past, the development of LAVs has been empirical. Blind passage of viruses in various cell types results in the accumulation of multiple attenuating mutations leaving the molecular mechanisms of attenuation unknown. Also, due to the high error rate of RNA viruses and selective pressures of the host environment, these LAVs, derived from such viruses, can potentially revert back to wild-type virulence. This not only puts the vaccinee at risk, but if shed can put those that are unvaccinated at risk as well. While these vaccines have been successful there still remains a need for a rational design strategy by which to create additional LAVs.One approach for rational vaccine design involves increasing the fidelity of the viral RdRp. Increased fidelity decreases the viral mutational frequency thereby reducing the genetic variation the virus needs in order to evade the host imposed bottlenecks to infection. While polymerase mutants exist which decrease viral mutation frequency the mutations are not in conserved regions of the polymerase, which doesn't lend itself toward using a common mutant approach toward developing a universal vaccine strategy for all RNA viruses. We have identified a conserved lysine residue in the active site of the PV RdRp that acts as a general acid during nucleotide incorporation. Mutation from a lysine to an arginine results in a high fidelity polymerase that replicates slowly thus creating an attenuated virus that is genetically stable and less likely to revert to a wild-type phenotype. This chapter provides detailed methods in which to identify the conserved lysine residue and evaluating fidelity and attenuation in cell culture (in vitro) and in the PV transgenic murine model (in vivo).

Performance assessment of the Illumina massively parallel sequencing platform for deep sequencing analysis of viral minority variants

Massively parallel sequencing (MPS) technology has opened new avenues to study viral dynamics and treatment-induced resistance mechanisms of infections such as human immunodeficiency virus (HIV) and hepatitis C virus (HCV). Whereas the Roche/454 platform has been used widely for the detection of low-frequent drug resistant variants, more recently developed short-read MPS technologies have the advantage of delivering a higher sequencing depth at a lower cost per sequenced base. This study assesses the performance characteristics of Illumina MPS technology for the characterization of genetic variability in viral populations by deep sequencing. The reported results from MPS experiments comprising HIV and HCV plasmids demonstrate that a 0.5-1% lower limit of detection can be achieved readily with Illumina MPS while retaining good accuracy also at low frequencies. Deep sequencing of a set of clinical samples (12 HIV and 9 HCV patients), designed at a similar budget for both MPS platforms, reveals a comparable lower limit of detection for Illumina and Roche/454. Finally, this study shows the possibility to apply Illumina's paired-end sequencing as a strategy to assess linkage between different mutations identified in individual viral subspecies. These results support the use of Illumina as another MPS platform of choice for deep sequencing of viral minority species.

Deconvoluting the composition of low-frequency hepatitis C viral quasispecies: comparison of genotypes and NS3 resistance-associated variants between HCV/HIV coinfected hemophiliacs and HCV monoinfected patients in Japan

Pre-existing low-frequency resistance-associated variants (RAVs) may jeopardize successful sustained virological responses (SVR) to HCV treatment with direct-acting antivirals (DAAs). However, the potential impact of low-frequency (∼0.1%) mutations, concatenated mutations (haplotypes), and their association with genotypes (Gts) on the treatment outcome has not yet been elucidated, most probably owing to the difficulty in detecting pre-existing minor haplotypes with sufficient length and accuracy. Herein, we characterize a methodological framework based on Illumina MiSeq next-generation sequencing (NGS) coupled with bioinformatics of quasispecies reconstruction (QSR) to realize highly accurate variant calling and genotype-haplotype detection. The core-to-NS3 protease coding sequences in 10 HCV monoinfected patients, 5 of whom had a history of blood transfusion, and 11 HCV/HIV coinfected patients with hemophilia, were studied. Simulation experiments showed that, for minor variants constituting more than 1%, our framework achieved a positive predictive value (PPV) of 100% and sensitivities of 91.7–100% for genotyping and 80.6% for RAV screening. Genotyping analysis indicated the prevalence of dominant Gt1a infection in coinfected patients (6/11 vs 0/10, p = 0.01). For clinical samples, minor genotype overlapping infection was prevalent in HCV/HIV coinfected hemophiliacs (10/11) and patients who experienced whole-blood transfusion (4/5) but none in patients without exposure to blood (0/5). As for RAV screening, the Q80K/R and S122K/R variants were particularly prevalent among minor RAVs observed, detected in 12/21 and 6/21 cases, respectively. Q80K was detected only in coinfected patients, whereas Q80R was predominantly detected in monoinfected patients (1/11 vs 7/10, p < 0.01). Multivariate interdependence analysis revealed the previously unrecognized prevalence of Gt1b-Q80K, in HCV/HIV coinfected hemophiliacs [Odds ratio = 13.4 (3.48–51.9), p < 0.01]. Our study revealed the distinct characteristics of viral quasispecies between the subgroups specified above and the feasibility of NGS and QSR-based genetic deconvolution of pre-existing minor Gts, RAVs, and their interrelationships.

The Quantification of Representative Sequences pipeline for amplicon sequencing: case study on within-population ITS1 sequence variation in a microparasite infecting Daphnia

Next generation sequencing (NGS) platforms are replacing traditional molecular biology protocols like cloning and Sanger sequencing. However, accuracy of NGS platforms has rarely been measured when quantifying relative frequencies of genotypes or taxa within populations. Here we developed a new bioinformatic pipeline (QRS) that pools similar sequence variants and estimates their frequencies in NGS data sets from populations or communities. We tested whether the estimated frequency of representative sequences, generated by 454 amplicon sequencing, differs significantly from that obtained by Sanger sequencing of cloned PCR products. This was performed by analysing sequence variation of the highly variable first internal transcribed spacer (ITS1) of the ichthyosporean Caullerya mesnili, a microparasite of cladocerans of the genus Daphnia. This analysis also serves as a case example of the usage of this pipeline to study within-population variation. Additionally, a public Illumina data set was used to validate the pipeline on community-level data. Overall, there was a good correspondence in absolute frequencies of C. mesnili ITS1 sequences obtained from Sanger and 454 platforms. Furthermore, analyses of molecular variance (amova) revealed that population structure of C. mesnili differs across lakes and years independently of the sequencing platform. Our results support not only the usefulness of amplicon sequencing data for studies of within-population structure but also the successful application of the QRS pipeline on Illumina-generated data. The QRS pipeline is freely available together with its documentation under GNU Public Licence version 3 at http://code.google.com/p/quantification-representative-sequences.

Evidence for weak genetic recombination at the PTP2 locus of Nosema ceranae

The microsporidian Nosema ceranae is an emergent pathogen that threatens the health of honeybees and other pollinators all over the world. Its recent rapid spread across a wide variety of host species and environments demonstrated an enhanced ability of adaptation, which seems to contradict the lack of evidence for genetic recombination and the absence of a sexual stage in its life cycle. Here we retrieved fresh data of the patterns of genetic variation at the PTP2 locus in naturally infected Apis mellifera colonies, by means of single genome amplification. This technique, designed to prevent the formation of chimeric haplotypes during polymerase chain reaction (PCR), provides more reliable estimates of the diversity levels and haplotype structure than standard PCR-cloning methods. Our results are consistent with low but significant rates of recombination in the history of the haplotypes detected: estimates of the population recombination rate are of the order of 30 and support recent evidence for unexpectedly high levels of variation of the parasites within honeybee colonies. These observations suggest the existence of a diploid stage at some point in the life cycle of this parasite and are relevant for our understanding of the dynamics of its expanding population.

Whole genome sequencing and de novo assembly identifies Sydney-like variant noroviruses and recombinants during the winter 2012/2013 outbreak in England

BACKGROUND

Norovirus is the commonest cause of epidemic gastroenteritis among people of all ages. Outbreaks frequently occur in hospitals and the community, costing the UK an estimated £110 m per annum. An evolutionary explanation for periodic increases in norovirus cases, despite some host-specific post immunity is currently limited to the identification of obvious recombinants. Our understanding could be significantly enhanced by full length genome sequences for large numbers of intensively sampled viruses, which would also assist control and vaccine design. Our objective is to develop rapid, high-throughput, end-to-end methods yielding complete norovirus genome sequences. We apply these methods to recent English outbreaks, placing them in the wider context of the international norovirus epidemic of winter 2012.

METHOD

Norovirus sequences were generated from 28 unique clinical samples by Illumina RNA sequencing (RNA-Seq) of total faecal RNA. A range of de novo sequence assemblers were attempted. The best assembler was identified by validation against three replicate samples and two norovirus qPCR negative samples, together with an additional 20 sequences determined by PCR and fractional capillary sequencing. Phylogenetic methods were used to reconstruct evolutionary relationships from the whole genome sequences.

RESULTS

Full length norovirus genomes were generated from 23/28 samples. 5/28 partial norovirus genomes were associated with low viral copy numbers. The de novo assembled sequences differed from sequences determined by capillary sequencing by <0.003%. Intra-host nucleotide sequence diversity was rare, but detectable by mapping short sequence reads onto its de novo assembled consensus. Genomes similar to the Sydney 2012 strain caused 78% (18/23) of cases, consistent with its previously documented association with the winter 2012 global outbreak. Interestingly, phylogenetic analysis and recombination detection analysis of the consensus sequences identified two related viruses as recombinants, containing sequences in prior circulation to Sydney 2012 in open reading frame (ORF) 2.

CONCLUSION

Our approach facilitates the rapid determination of complete norovirus genomes. This method provides high resolution of full norovirus genomes which, when coupled with detailed epidemiology, may improve the understanding of evolution and control of this important healthcare-associated pathogen.

Ultradeep pyrosequencing of hepatitis C virus hypervariable region 1 in quasispecies analysis

Genetic variability of hepatitis C virus (HCV) determines pathogenesis of infection, including viral persistence and resistance to treatment. The aim of the present study was to characterize HCV genetic heterogeneity within a hypervariable region 1 (HVR1) of a chronically infected patient by ultradeep 454 sequencing strategy. Three independent sequencing error correction methods were applied. First correction method (Method I) implemented cut-off for genetic variants present in less than 1%. In the second method (Method II), a condition to call a variant was bidirectional coverage of sequencing reads. Third method (Method III) used Short Read Assembly into Haplotypes (ShoRAH) program. After the application of these three different algorithms, HVR1 population consisted of 8, 40, and 186 genetic haplotypes. The most sensitive method was ShoRAH, allowing to reconstruct haplotypes constituting as little as 0.013% of the population. The most abundant genetic variant constituted only 10.5%. Seventeen haplotypes were present in a frequency above 1%, and there was wide dispersion of the population into very sparse haplotypes. Our results indicate that HCV HVR1 heterogeneity and quasispecies population structure may be reconstructed by ultradeep sequencing. However, credible analysis requires proper reconstruction methods, which would distinguish sequencing error from real variability in vivo.

Intra-host variation structure of classical swine fever virus NS5B in relation to antiviral therapy

Classical swine fever (CSF) is one of most important diseases of the Suidea with severe social economic consequences in case of outbreaks. Antivirals have been demonstrated, in recent publications, to be an interesting alternative method of fighting the disease. However, classical swine fever virus is an RNA virus which presents a challenge as intra-host variation and the error prone RNA dependent RNA polymerase (RdRp) could lead to the emergence/selection of resistant variants hampering further treatment. Therefore, it was the purpose of this study to investigate the intra-host variation of the RdRp gene, targeted by antivirals, in respect to antiviral treatment. Using the non-unique nucleotide changes, a limited intra-host variation was found in the wild type virus with 2 silent and 2 non-synonymous sites. This number shifted significantly when an antiviral resistant variant was analyzed. In total 22nt changes were found resulting in 14 amino acid changes whereby each genome copy contained at least 2 amino-acid changes in the RdRp. Interestingly, the frequency of the mutations situated in close proximity to a region involved in antiviral resistance in CSFV and bovine viral diarrhea virus (BVDV) was elevated compared to the other mutations. None of the identified mutations in the resistant variant and which could potentially result in antiviral resistance was present in the wild type virus as a non-unique mutation. In view of the spectrum of mutations identified in the resistance associated region and that none of the resistance associated mutations reported for another strain of classical swine fever for the same antiviral were observed in the study, it can be suggested that multiple mutations confer resistance to some degree. Although the followed classical approach allowed the analysis the RdRp as a whole, the contribution of unique mutations to the intra-host variation could not be completely resolved. There was a significant difference in de number of unique mutations found between: 1/wild type virus and the antiviral resistant variant and 2/between both and the number to be expected from the error rate of the RT-PCR process. This indicates that the some of the unique mutations contributed to the intra-host variation and that the antiviral pressure also shifted this pattern. This is important as one of the non-synonymous mutations found in the resistant variant and which was located in the antiviral resistance associated region, was present in the wild type virus as a unique mutation. The findings presented in this study not only show the importance of intra-host variation analysis but also warrants further research certainly in view of the potential inclusion of antivirals in a control/eradication strategy.

Variability and resistance mutations in the hepatitis C virus NS3 protease in patients not treated with protease inhibitors

The goal of treatment of chronic hepatitis C is to achieve a sustained virological response, which is defined as exhibiting undetectable hepatitis C virus (HCV) RNA levels in serum following therapy for at least six months. However, the current treatment is only effective in 50% of patients infected with HCV genotype 1, the most prevalent genotype in Brazil. Inhibitors of the serine protease non-structural protein 3 (NS3) have therefore been developed to improve the responses of HCV-infected patients. However, the emergence of drug-resistant variants has been the major obstacle to therapeutic success. The goal of this study was to evaluate the presence of resistance mutations and genetic polymorphisms in the NS3 genomic region of HCV from 37 patients infected with HCV genotype 1 had not been treated with protease inhibitors. Plasma viral RNA was used to amplify and sequence the HCV NS3 gene. The results indicate that the catalytic triad is conserved. A large number of substitutions were observed in codons 153, 40 and 91; the resistant variants T54A, T54S, V55A, R155K and A156T were also detected. This study shows that resistance mutations and genetic polymorphisms are present in the NS3 region of HCV in patients who have not been treated with protease inhibitors, data that are important in determining the efficiency of this new class of drugs in Brazil.

Quantifying the diversification of hepatitis C virus (HCV) during primary infection: estimates of the in vivo mutation rate

Hepatitis C virus (HCV) is present in the host with multiple variants generated by its error prone RNA-dependent RNA polymerase. Little is known about the initial viral diversification and the viral life cycle processes that influence diversity. We studied the diversification of HCV during acute infection in 17 plasma donors, with frequent sampling early in infection. To analyze these data, we developed a new stochastic model of the HCV life cycle. We found that the accumulation of mutations is surprisingly slow: at 30 days, the viral population on average is still 46% identical to its transmitted viral genome. Fitting the model to the sequence data, we estimate the median in vivo viral mutation rate is 2.5×10⁻⁵ mutations per nucleotide per genome replication (range 1.6–6.2×10⁻⁵), about 5-fold lower than previous estimates. To confirm these results we analyzed the frequency of stop codons (N = 10) among all possible non-sense mutation targets (M = 898,335), and found a mutation rate of 2.8–3.2×10⁻⁵, consistent with the estimate from the dynamical model. The slow accumulation of mutations is consistent with slow turnover of infected cells and replication complexes within infected cells. This slow turnover is also inferred from the viral load kinetics. Our estimated mutation rate, which is similar to that of other RNA viruses (e.g., HIV and influenza), is also compatible with the accumulation of substitutions seen in HCV at the population level. Our model identifies the relevant processes (long-lived cells and slow turnover of replication complexes) and parameters involved in determining the rate of HCV diversification.

Hepatitis C virus (HCV) is a RNA virus that infects over 170 million people across the world. It leads to a chronic infection in the majority of people who are infected (>70%). Most people only discover that they are infected long after initial infection. Thus, it is difficult to study the very early events in infection. Here we study 17 individuals during the earliest possible stages of infection, from before the virus is detectable in the plasma to around 35 days post-infection. We focus on understanding the viral kinetics and the diversification of HCV during this acute phase of infection. During chronic infection HCV is present in the host as a swarm of multiple variants generated by its error prone copying. We studied the early diversification of HCV during acute infection using a new mathematical model of HCV replication. We found that after a phase of fast increase in viral load, accompanied by viral diversification, there is a stabilization of viral load and diversity levels. Using our model, we were able to estimate for the first time the HCV mutation rate during acute infection. We estimated the median in vivo viral mutation rate is 2.5×10⁻⁵ mutations per nucleotide per genome replication (range 1.6–6.2×10⁻⁵), about 5-fold lower than previous estimates. We also used a different approach, based on results of classical genetics, to calculate HCV's mutation rate and obtained consistent results (2.8–3.2×10⁻⁵).

Elucidation of hepatitis C virus transmission and early diversification by single genome sequencing

A precise molecular identification of transmitted hepatitis C virus (HCV) genomes could illuminate key aspects of transmission biology, immunopathogenesis and natural history. We used single genome sequencing of 2,922 half or quarter genomes from plasma viral RNA to identify transmitted/founder (T/F) viruses in 17 subjects with acute community-acquired HCV infection. Sequences from 13 of 17 acute subjects, but none of 14 chronic controls, exhibited one or more discrete low diversity viral lineages. Sequences within each lineage generally revealed a star-like phylogeny of mutations that coalesced to unambiguous T/F viral genomes. Numbers of transmitted viruses leading to productive clinical infection were estimated to range from 1 to 37 or more (median = 4). Four acutely infected subjects showed a distinctly different pattern of virus diversity that deviated from a star-like phylogeny. In these cases, empirical analysis and mathematical modeling suggested high multiplicity virus transmission from individuals who themselves were acutely infected or had experienced a virus population bottleneck due to antiviral drug therapy. These results provide new quantitative and qualitative insights into HCV transmission, revealing for the first time virus-host interactions that successful vaccines or treatment interventions will need to overcome. Our findings further suggest a novel experimental strategy for identifying full-length T/F genomes for proteome-wide analyses of HCV biology and adaptation to antiviral drug or immune pressures.

Hepatitis C virus infects as many as 170 million people worldwide. Globally, there are seven major genotypes of HCV that differ by approximately 30% in nucleotide sequence. Importantly, the natural history of HCV infection is variable, ranging from spontaneous resolution to persistent viremia and chronic disease. Factors responsible for this variability in clinical outcome are unknown but likely involve a combination of viral and host determinants. To this end, a precise molecular identification of transmitted HCV genomes could illuminate key aspects of transmission biology, immunopathogenesis and natural history. We used single genome sequencing of plasma viral RNA to identify transmitted viral genomes and their progeny in 17 subjects with acute infection. Numbers of transmitted viruses leading to productive clinical infection ranged from 1 to 37 or more (median = 4). Surprisingly, we found evidence of high multiplicity acute-to-acute HCV transmission in 3 of 17 subjects, which suggests that clinical transmission of HCV, like that of HIV-1, may be enhanced in early infection when virus titers are highest and neutralizing antibodies are absent. These results provide novel insight into HCV transmission and early virus diversification key to our understanding of virus natural history and response to drug selection and immune pressure.

High-resolution analysis of intrahost genetic diversity in dengue virus serotype 1 infection identifies mixed infections

Little is known about the rate at which genetic variation is generated within intrahost populations of dengue virus (DENV) and what implications this diversity has for dengue pathogenesis, disease severity, and host immunity. Previous studies of intrahost DENV variation have used a low frequency of sampling and/or experimental methods that do not fully account for errors generated through amplification and sequencing of viral RNAs. We investigated the extent and pattern of genetic diversity in sequence data in domain III (DIII) of the envelope (E) gene in serial plasma samples (n = 49) taken from 17 patients infected with DENV type 1 (DENV-1), totaling some 8,458 clones. Statistically rigorous approaches were employed to account for artifactual variants resulting from amplification and sequencing, which we suggest have played a major role in previous studies of intrahost genetic variation. Accordingly, nucleotide sequence diversities of viral populations were very low, with conservative estimates of the average levels of genetic diversity ranging from 0 to 0.0013. Despite such sequence conservation, we observed clear evidence for mixed infection, with the presence of multiple phylogenetically distinct lineages present within the same host, while the presence of stop codon mutations in some samples suggests the action of complementation. In contrast to some previous studies we observed no relationship between the extent and pattern of DENV-1 genetic diversity and disease severity, immune status, or level of viremia.

No evidence of XMRV in prostate cancer cohorts in the Midwestern United States

Background

Xenotropic murine leukemia virus (MLV)-related virus (XMRV) was initially identified in prostate cancer (PCa) tissue, particularly in the prostatic stromal fibroblasts, of patients homozygous for the RNASEL R462Q mutation. A subsequent study reported XMRV antigens in malignant prostatic epithelium and association of XMRV infection with PCa, especially higher-grade tumors, independently of the RNASEL polymorphism. Further studies showed high prevalence of XMRV or related MLV sequences in chronic fatigue syndrome patients (CFS), while others found no, or low, prevalence of XMRV in a variety of diseases including PCa or CFS. Thus, the etiological link between XMRV and human disease remains elusive. To address the association between XMRV infection and PCa, we have tested prostate tissues and human sera for the presence of viral DNA, viral antigens and anti-XMRV antibodies.

Results

Real-time PCR analysis of 110 PCa (Gleason scores >4) and 40 benign and normal prostate tissues identified six positive samples (5 PCa and 1 non-PCa). No statistical link was observed between the presence of proviral DNA and PCa, PCa grades, and the RNASEL R462Q mutation. The amplified viral sequences were distantly related to XMRV, but nearly identical to endogenous MLV sequences in mice. The PCR positive samples were also positive for mouse mitochondrial DNA by nested PCR, suggesting contamination of the samples with mouse DNA. Immuno-histochemistry (IHC) with an anti-XMRV antibody, but not an anti-MLV antibody that recognizes XMRV, sporadically identified antigen-positive cells in prostatic epithelium, irrespectively of the status of viral DNA detection. No serum (159 PCa and 201 age-matched controls) showed strong neutralization of XMRV infection at 1:10 dilution.

Conclusion

The lack of XMRV sequences or strong anti-XMRV neutralizing antibodies indicates no or very low prevalence of XMRV in our cohorts. We conclude that real-time PCR- and IHC-positive samples were due to laboratory contamination and non-specific immune reactions, respectively.

Beyond the consensus: dissecting within-host viral population diversity of foot-and-mouth disease virus by using next-generation genome sequencing

The diverse sequences of viral populations within individual hosts are the starting material for selection and subsequent evolution of RNA viruses such as foot-and-mouth disease virus (FMDV). Using next-generation sequencing (NGS) performed on a Genome Analyzer platform (Illumina), this study compared the viral populations within two bovine epithelial samples (foot lesions) from a single animal with the inoculum used to initiate experimental infection. Genomic sequences were determined in duplicate sequencing runs, and the consensus sequence of the inoculum determined by NGS was identical to that previously determined using the Sanger method. However, NGS revealed the fine polymorphic substructure of the viral population, from nucleotide variants present at just below 50% frequency to those present at fractions of 1%. Some of the higher-frequency polymorphisms identified encoded changes within codons associated with heparan sulfate binding and were present in both foot lesions, revealing intermediate stages in the evolution of a tissue culture-adapted virus replicating within a mammalian host. We identified 2,622, 1,434, and 1,703 polymorphisms in the inoculum and in the two foot lesions, respectively: most of the substitutions occurred in only a small fraction of the population and represented the progeny from recent cellular replication prior to onset of any selective pressures. We estimated the upper limit for the genome-wide mutation rate of the virus within a cell to be 7.8 × 10(-4) per nucleotide. The greater depth of detection achieved by NGS demonstrates that this method is a powerful and valuable tool for the dissection of FMDV populations within hosts.

Homogeneity of Powassan virus populations in naturally infected Ixodes scapularis

Powassan virus (POWV, Flaviviridae: Flavivirus) is the sole North American member of the tick-borne encephalitis complex and consists of two distinct lineages that are maintained in ecologically discrete enzootic transmission cycles. The underlying genetic mechanisms that lead to niche partitioning in arboviruses are poorly understood. Therefore, intra- and interhost genetic diversity was analyzed to determine if POWV exists as a quasispecies in nature and quantify selective pressures within and between hosts. In contrast to previous reports for West Nile virus (WNV), significant intrahost genetic diversity was not observed. However, pN (0.238) and d(N)/d(S) ratios (0.092) for interhost diversity were similar to those of WNV. Combined, these data suggest that purifying selection and/or population bottlenecks constrain quasispecies diversity within ticks. These same selective and stochastic mechanisms appear to drive minor sequence changes between ticks. Moreover, Powassan virus populations seem not to be structured as quasispecies in naturally infected adult deer ticks.

Natural epitope variants of the hepatitis C virus impair cytotoxic T lymphocyte activity

AIM: To understand how interactions between hepatitis C virus (HCV) and the host’s immune system might lead to viral persistence or effective elimination of HCV.

METHODS: Nucleotides 3519-3935 of the non-structural 3 (NS3) region were amplified by using reverse transcription polymerase chain reaction (PCR). PCR products of the HCV NS3 regions were integrated into a PCR^® T7TOPO^® TA vector and then sequenced in both directions using an automated DNA sequencer. Relative major histocompatibility complex binding levels of wild-type and variant peptides were performed by fluorescence polarization-based peptide competition assays. Peptides with wild type and variant sequences of NS3 were synthesized locally using F-moc chemistry and purified by high-performance liquid chromatography. Specific cytotoxic T lymphocytes (CTLs) clones toward HCV NS3 wild-type peptides were generated through limiting dilution cloning. The CTL clones specifically recognizing HCV NS3 wild-type peptides were tested by tetramer staining and flow cytometry. Cytolytic activity of CTL clones was measured using target cells labeled with the fluorescence enhancing ligand, DELFIA EuTDA.

RESULTS: The pattern of natural variants within three human leukocyte antigen (HLA)-A2-restricted NS3 epitopes has been examined in one patient with chronic HCV infection at 12, 28 and 63 mo post-infection. Results obtained may provide convincing evidence of immune selection pressure for all epitopes investigated. Statistical analysis of the extensive sequence variation found within these NS3 epitopes favors a Darwinian selection model of variant viruses. Mutations within the epitopes coincided with the decline of CTL responses, and peptide-binding studies suggested a significant impact of the mutation on T cell recognition rather than peptide presentation by HLA molecules. While most variants were either not recognized or elicited low responses, such could antagonize CTL responses to target cells pulsed with wild-type peptides.

CONCLUSION: Cross-recognition of CTL epitopes from wild-type and naturally-occurring HCV variants may lead to impaired immune responses and ultimately contribute to viral persistence.

Low genetic diversities of rabies virus populations within different hosts in Brazil

The low rates of nonsynonymous evolution observed in natural rabies virus (RABV) isolates are suggested to have arisen in association with the structural and functional constraints operating on the virus protein and the infection strategies employed by RABV within infected hosts to avoid strong selection by the immune response. In order to investigate the relationship between the genetic characteristics of RABV populations within hosts and the virus evolution, the present study examined the genetic heterogeneities of RABV populations within naturally infected dogs and foxes in Brazil, as well as those of bat RABV populations that were passaged once in suckling mice. Sequence analyses of complete RABV glycoprotein (G) genes showed that RABV populations within infected hosts were genetically highly homogeneous whether they were infected naturally or experimentally (nucleotide diversities of 0-0.95x10(-3)). In addition, amino acid mutations were randomly distributed over the entire region of the G protein, and the nonsynonymous/synonymous rate ratios (d(N)/d(S)) for the G protein gene were less than 1. These findings suggest that the low genetic diversities of RABV populations within hosts reflect the stabilizing selection operating on the virus, the infection strategies of the virus, and eventually, the evolutionary patterns of the virus.

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.

Analysis of Foot-and-mouth disease virus nucleotide sequence variation within naturally infected epithelium

An important epidemiological tool in the control of epidemics of Foot-and-mouth disease (FMD) is genetic tracing using complete virus genome sequence data. However to interpret these genetic data, it is important to quantify underlying variation present in FMDV populations from individual tissue samples. Cloned complete capsid sequences from two virus populations from epithelium from a cow (n=26) and from a sheep (n=15) infected during the UK 2001 outbreaks were generated. Genetic diversity of the two virus populations differed significantly, with sequences representing virus from the cow having a mutation frequency of 2.79 x 10(-4) mutations per nucleotide sequenced (mpns) and those from the sheep having 3.94 x 10(-4) mpns (chi(2)=8.24, P=0.004). The dN/dS ratio of sequences from the cow was higher (1.228) than that from the sheep (0.187) although not significantly so. The sequences from the cow epithelium exhibited significantly higher than expected number of changes within neutralising antigenic sites (P=0.0007). The performance of two different reverse transcriptase enzymes was found not to differ with respect to the frequency (P=0.559, chi(2)=0.341) or dN/dS ratio (P=0.863, chi(2)=0.03) of the mutations observed. These data provide insight into the population diversity that exists within a single lesion and help toward understanding the mechanisms that underpin sequence evolution of FMDV.

Genetic heterogeneity of L-Zagreb mumps virus vaccine strain

BACKGROUND

The most often used mumps vaccine strains Jeryl Lynn (JL), RIT4385, Urabe-AM9, L-Zagreb and L-3 differ in immunogenicity and reactogenicity. Previous analyses showed that JL, Urabe-AM9 and L-3 are genetically heterogeneous.

RESULTS

We identified the heterogeneity of L-Zagreb throughout the entire genome. Two major variants were defined: variant A being identical to the consensus sequence of viral seeds and vaccine(s) and variant B which differs from variant A in three nucleotide positions. The difference between viral variants in L-Zagreb strain is insufficient for distinct viral strains to be defined. We demonstrated that proportion of variants in L-Zagreb viral population depends on cell substrate used for viral replication in vitro and in vivo.

CONCLUSION

L-Zagreb strain should be considered as a single strain composed of at least two variant viral genomes.

Viral population estimation using pyrosequencing

The diversity of virus populations within single infected hosts presents a major difficulty for the natural immune response as well as for vaccine design and antiviral drug therapy. Recently developed pyrophosphate-based sequencing technologies (pyrosequencing) can be used for quantifying this diversity by ultra-deep sequencing of virus samples. We present computational methods for the analysis of such sequence data and apply these techniques to pyrosequencing data obtained from HIV populations within patients harboring drug-resistant virus strains. Our main result is the estimation of the population structure of the sample from the pyrosequencing reads. This inference is based on a statistical approach to error correction, followed by a combinatorial algorithm for constructing a minimal set of haplotypes that explain the data. Using this set of explaining haplotypes, we apply a statistical model to infer the frequencies of the haplotypes in the population via an expectation–maximization (EM) algorithm. We demonstrate that pyrosequencing reads allow for effective population reconstruction by extensive simulations and by comparison to 165 sequences obtained directly from clonal sequencing of four independent, diverse HIV populations. Thus, pyrosequencing can be used for cost-effective estimation of the structure of virus populations, promising new insights into viral evolutionary dynamics and disease control strategies.

The genetic diversity of viral populations is important for biomedical problems such as disease progression, vaccine design, and drug resistance, yet it is not generally well understood. In this paper, we use pyrosequencing, a novel DNA sequencing technique, to reconstruct viral populations. Pyrosequencing produces DNA sequences, called reads, in numbers much greater than standard DNA sequencing techniques. However, these reads are substantially shorter and more error-prone than those obtained from standard sequencing techniques. Therefore, pyrosequencing data requires new methods of analysis. Here, we develop mathematical and statistical tools for reconstructing viral populations using pyrosequencing. To this end, we show how to correct errors in the reads and assemble them into the different viral strains present in the population. We apply these methods to HIV-1 populations from drug-resistant patients and show that our techniques produce results quite close to accepted techniques at a lower cost and potentially higher resolution.

Protocol for nearly full-length sequencing of HIV-1 RNA from plasma

Nearly full-length genome sequencing of HIV-1 using peripheral blood mononuclear cells (PBMC) DNA as a template for PCR is now a relatively routine laboratory procedure. However, this has not been the case when using virion RNA as the template and this has made full genome analysis of circulating viruses difficult. Therefore, a well-developed procedure for sequencing of full-length HIV-1 RNA directly from plasma was needed. Plasma from U.S. donors representing a range of viral loads (VL) was used to develop the assay. RNA was extracted from plasma and reverse-transcribed. Two or three overlapping regions were PCR amplified to cover the entire viral genome and sequenced for verification. The success of the procedure was sensitive to VL but was routinely successful for VL greater than 10⁵ and the rate declined in proportion to the VL. While the two-amplicon strategy had an advantage of increasing the possibility of amplifying a single species of HIV-1, the three-amplicon strategy was more successful in amplifying samples with low viral loads. This protocol provides a useful tool for molecular analysis to understand the HIV epidemic and pathogenesis, as well as diagnosis, therapy and future vaccine strategies.

Characterization of mutation spectra with ultra-deep pyrosequencing: application to HIV-1 drug resistance

The detection of mutant spectra within a population of microorganisms is critical for the management of drug-resistant infections. We performed ultra-deep pyrosequencing to detect minor sequence variants in HIV-1 protease and reverse transcriptase (RT) genes from clinical plasma samples. We estimated empirical error rates from four HIV-1 plasmid clones and used them to develop a statistical approach to distinguish authentic minor variants from sequencing errors in eight clinical samples. Ultra-deep pyrosequencing detected an average of 58 variants per sample compared with an average of eight variants per sample detected by conventional direct-PCR dideoxynucleotide sequencing. In the clinical sample with the largest number of minor sequence variants, all 60 variants present in > or =3% of genomes and 20 of 35 variants present in <3% of genomes were confirmed by limiting dilution sequencing. With appropriate analysis, ultra-deep pyrosequencing is a promising method for characterizing genetic diversity and detecting minor yet clinically relevant variants in biological samples with complex genetic populations.

Quasispecies evolution in NS5A region of hepatitis C virus genotype 1b during interferon or combined interferon-ribavirin therapy

AIM: To evaluate the implication of substitutions in the hepatitis C virus (HCV) non-structural 5A (NS5A) protein in the resistance of HCV during mono-interferon (IFN) or combined IFN-ribavirin (IFN-R) therapy. Although NS5A has been reported to interact with the HCV RNA-dependent RNA polymerase, NS5B, as well as with many cellular proteins, the function of NS5A in the life cycle of HCV remains unclear.

METHODS: HCV quasispecies were studied by cloning and sequencing of sequential isolates from patients infected by HCV genotype 1b. Patients were treated by IFN-α2b for 3 mo followed by IFN-α2b alone or combined IFN-R therapy for 9 additional months. Patients were categorized intro two groups based on their response to the treatments: 7 with sustained virological response (SVR) (quasispecies = 150) and 3 non-responders (NR) to IFN-R (quasispecies = 106).

RESULTS: Prior to treatment, SVR patients displayed a lower complexity of quasispecies than NR patients. Most patients had a decrease in the complexity of quasispecies during therapy. Analysis of amino acids substitutions showed that the degree of the complexity of the interferon sensitivity-determining region (ISDR) and the V3 domain of NS5A protein was able to discriminate the two groups of patients. Moreover, SVR patients displayed more variability in the NS5A region than NR patients.

CONCLUSION: These results suggest that detailed molecular analysis of the NS5A region may be important for understanding its function in IFN response during HCV 1b infection.

Anti-retroviral drug resistance-associated mutations among non-subtype B HIV-1-infected Kenyan children with treatment failure

Recently increased availability of anti-retroviral therapy (ART) has mitigated HIV-1/AIDS prognoses especially in resource poor settings. The emergence of ART resistance-associated mutations from non-suppressive ART has been implicated as a major cause of ART failure. Reverse transcriptase inhibitor (RTI)-resistance mutations among 12 non-subtype B HIV-1-infected children with treatment failure were evaluated by genotypically analyzing HIV-1 strains isolated from plasma obtained between 2001 and 2004. A region of pol-RT gene was amplified and at least five clones per sample were analyzed. Phylogenetic analysis revealed HIV-1 subtype A1 (n = 7), subtype C (n = 1), subtype D (n = 3), and CRF02_AG (n = 1). Before treatment, 4 of 12 (33.3%) children had primary RTI-resistance mutations, K103N (n = 3, ages 5-7 years) and Y181C (n = 1, age 1 year). In one child, K103N was found as a minor population (1/5 clones) before treatment and became major (7/7 clones) 8 months after RTI treatment. In 7 of 12 children, M184V appeared with one thymidine-analogue-associated mutation (TAM) as the first mutation, while the remaining 5 children had only TAMs appearing either individually (n = 2), or as TAMs 1 (M41L, L210W, and T215Y) and 2 (D67N, K70R, and K219Q/E/R) appearing together (n = 3). These results suggest that "vertically transmitted" primary RTI-resistance mutations, K103N and Y181C, can persist over the years even in the absence of drug pressure and impact RTI treatment negatively, and that appearing patterns of RTI-resistance mutations among non-subtype B HIV-1-infected children could possibly be different from those reported in subtype B-infected children.

Mutations in the MYB intron I regulatory sequence increase transcription in colon cancers

Although MYB overexpression in colorectal cancer (CRC) is known to be a prognostic indicator for poor survival, the basis for this overexpression is unclear. Among multiple levels of MYB regulation, the most dynamic is the control of transcriptional elongation by sequences within intron 1. The authors have proposed that this regulatory sequence is transcribed into an RNA stem-loop and 19-residue polyuridine tract, and is subject to mutation in CRC. When this region was examined in colorectal and breast carcinoma cell lines and tissues, the authors found frequent mutations only in CRC. It was determined that these mutations allowed increased transcription compared with the wild type sequence. These data suggest that this MYB regulatory region within intron 1 is subject to mutations in CRC but not breast cancer, perhaps consistent with the mutagenic insult that occurs within the colon and not mammary tissue. In CRC, these mutations may contribute to MYB overexpression, highlighting the importance of noncoding sequences in the regulation of key cancer genes.

Escherichia coli DNA polymerase III epsilon subunit increases Moloney murine leukemia virus reverse transcriptase fidelity and accuracy of RT-PCR procedures

In an effort to improve reverse transcriptase (RT) fidelity, we measured the error rate of Moloney murine leukemia virus (MMLV) RT in the presence of several autonomous and DNA polymerase-associated 3'-5' exonucleases using a lacZ forward mutation assay. A number of 3'-5' exonucleases were found to lower the error rate of MMLV RT, including p53, Escherichia coli DNA polymerase III epsilon subunit, and the proofreading activities associated with T4, varphi29, and E. coli pol I DNA polymerases. The bacterial epsilon subunit increased RNA-dependent DNA synthesis fidelity by approximately threefold and was the only 3'-5' exonuclease tested that did not deleteriously affect RT-PCR yields. Further testing showed that RT-PCR mutant frequencies were reduced significantly by performing cDNA synthesis in the presence of epsilon subunit, followed by PCR with a high-fidelity proofreading DNA polymerase. DNA sequence analysis was used to show that the combination of MMLV RT/epsilon subunit and PfuUltra DNA polymerase produces approximately eightfold fewer errors compared with the commonly used combination of MMLV RT and a Taq-based high-fidelity blend, consistent with predictions based on experimentally determined polymerase error rates.

Selection and persistence of non-nucleoside reverse transcriptase inhibitor-resistant HIV-1 in patients starting and stopping non-nucleoside therapy

BACKGROUND

Understanding the selection and decay of drug-resistant HIV-1 variants is important for designing optimal antiretroviral therapy.

OBJECTIVE

To develop a high-throughput, real-time reverse transcriptase (RT) polymerase chain reaction (PCR) assay to quantify non-nucleoside reverse transcriptase inhibitor (NNRTI)-resistant variants K103N (AAT or AAC alleles) at frequencies as low as 0.1%, and to apply this to monitor these variants before, during, and after NNRTI therapy.

METHODS

HIV-1 RNA in longitudinal plasma samples obtained from patients starting and stopping NNRTI therapy was converted to cDNA and the target sequence region amplified and quantified by real-time PCR. Approximately 10 copies/reaction provided a template for a second round of PCR using primers that discriminated between the mutant and wild-type alleles. Amplification specificity was confirmed by thermal denaturation analysis.

RESULTS

Frequencies of 103N similar to assay background (0.029%) were observed in longitudinal samples from 9 of 12 treatment-naive patients; three patients had transient increases in 103N frequency to a range of 0.21-0.48%, which was 7-16.5 times assay background. Analysis of longitudinal plasma samples from six NNRTI-experienced patients showed three patterns: persistence of 103N variants after stopping NNRTI therapy, codon switching of 103N between AAC and AAT during NNRTI therapy, and decay of 103N variants to below assay background after cessation of NNRTI therapy.

CONCLUSIONS

Allele-specific RT-PCR quantified the emergence and decay of drug-resistant variants in patients over a broad range of frequencies (0.1-100%). The rate of decay of K103N variants after stopping NNRTI therapy was highly variable.

Hepatitis C virus in the semen of men coinfected with HIV-1: prevalence and origin

OBJECTIVE

To compare the prevalence of hepatitis C (HCV) RNA in semen from men infected with HCV and those coinfected with HIV-1/HCV and to study the origin of HCV shed in semen.

DESIGN

Two prospective studies (HC EP09 and BINECO) included 120 HCV-positive men, 82 coinfected with HIV-1; all had positive HCV RNA detection in blood.

METHODS

Paired blood and semen samples were collected for HCV RNA detection and quantification in seminal plasma and in blood serum; repeated semen samples were obtained for 45 men. HCV RNA was sought in spermatozoa and non-sperm cells. Phylogenetic analysis of the HVR-1 region of HCV compared the quasispecies in blood serum and seminal plasma of two men.

RESULTS

HCV RNA was more frequently found in the semen of men coinfected with HIV-1 (37.8%) than in those with only HCV infection (18.4%) (P = 0.033). HCV RNA detection in semen was intermittent and was positive in at least one semen sample of 42.8% of HIV-1/HCV-coinfected men who provided repeated samples. Men with HCV-positive semen had significantly higher HCV load in blood than men with HCV-negative semen (P = 0.038). Phylogenetic comparison of HCV quasispecies in blood and in semen showed no evidence of HCV replication in genital leukocytes; however, a phenetic structure was observed between compartments (P < 0.001).

CONCLUSIONS

HCV particles in semen originate from passive passage from blood, with preferential transfer of some variants. Nearly half of HIV-1/HCV-coinfected men may intermittently harbour HCV in their semen. Recommendations of protected sex for HIV-infected individuals should be reinforced.

Genetic variation in West Nile virus from naturally infected mosquitoes and birds suggests quasispecies structure and strong purifying selection

Intrahost genetic diversity was analysed in naturally infected mosquitoes and birds to determine whether West Nile virus (WNV) exists in nature as a quasispecies and to quantify selective pressures within and between hosts. WNV was sampled from ten infected birds and ten infected mosquito pools collected on Long Island, NY, USA, during the peak of the 2003 WNV transmission season. A 1938 nt fragment comprising the 3' 1159 nt of the WNV envelope (E) coding region and the 5' 779 nt of the non-structural protein 1 (NS1) coding region was amplified and cloned and 20 clones per specimen were sequenced. Results from this analysis demonstrate that WNV infections are derived from a genetically diverse population of genomes in nature. The mean nucleotide diversity was 0.016 % within individual specimens and the mean percentage of clones that differed from the consensus sequence was 19.5 %. WNV sequences in mosquitoes were significantly more genetically diverse than WNV in birds. No host-dependent bias for particular types of mutations was observed and estimates of genetic diversity did not differ significantly between E and NS1 coding sequences. Non-consensus clones obtained from two avian specimens had highly similar genetic signatures, providing preliminary evidence that WNV genetic diversity may be maintained throughout the enzootic transmission cycle, rather than arising independently during each infection. Evidence of purifying selection was obtained from both intra- and interhost WNV populations. Combined, these data support the observation that WNV populations may be structured as a quasispecies and document strong purifying natural selection in WNV populations.

High levels of discordance between sequencing and serological subtyping in a predominantly non-B subtype HIV-1 infected cohort

Samples from 457 randomly selected HIV-1 infected patients attending King's College Hospital were analysed using a subtype specific enzyme immunoassay. All serotyped non-Bs that provided unambiguous sequence and for which sufficient sample was available (n=100), which included three serotyped subtype B samples were further analysed by env sequencing and subtyping using neighbour joining phylogenetic analysis, the NCBI Retrovirus Genotyping tool and the Los Alamos BLAST search tool. Of the serotyped viruses, 45% (n=204) samples were subtype B. Specifically serotyped non-B strains (n=130) accounted for 28% of the total, of which the largest proportion were subtype C (n=66). Twenty-seven samples (6%) were classified as non-B, 9% (n=40) were multiply-reactive and 12% were non-reactive (n=56). Of the 100 samples subtyped by sequencing the majority were subtype C (n=32), followed by subtype A (n=20). There was little concordance between the two methods. Although a 100% match was found among the serotyped and sequenced non-B viruses (n=13), only 16 of the sequenced subtype C specimens matched the 29 obtained by serotyping. Of the 20 multiply-reactive samples analysed by serotyping, only 1 sample consisted of a subtype mixture by sequencing. Of the 14 serologically non-reactive samples analysed, all were successfully sequenced, with subtype B strains (57%) the most common. Sequencing 15 samples in both env and pol regions revealed differences in subtype assignment for the same sample in some cases. Only 1/6 env subtype A and 4/5 env subtype C samples were concordant in pol sequence subtype. Differences were also found in subtyping by the different methods used. The overall agreement between the three methods was 89%. Four out of 11 samples agreed between the phylogenetic and Los Alamos methods, 1/11 between phylogenetic and BLAST and 2/11 between Los Alamos and BLAST.

Low human immunodeficiency virus envelope diversity correlates with low in vitro replication capacity and predicts spontaneous control of plasma viremia after treatment interruptions

Genetic diversity of viral isolates in human immunodeficiency virus (HIV)-infected individuals varies substantially. However, it remains unclear whether HIV-related disease progresses more rapidly in patients harboring virus swarms with low or high diversity and, in the same context, whether high or low diversity is required to induce potent humoral and cellular immune responses. To explore whether viral diversity predicts virologic control, we studied HIV-infected patients who received antiretroviral therapy (ART) for years before undergoing structured treatment interruptions (STI). Viral diversity before initiation of ART and the ability of the patients to contain viremia after STI and final cessation of treatment was evaluated. Seven out of 21 patients contained plasma viremia at low levels after the final treatment cessation. Clonal sequences encompassing the envelope C2V3C3 domain derived from plasma prior to treatment, exhibited significantly lower diversity in these patients compared to those derived from patients with poor control of viremia. Viral diversity pre-ART correlated with the viral replication capacity of rebounding virus isolates during STI. Neutralizing antibody activity against autologous virus was significantly higher in patients who controlled viremia and was associated with lower pretreatment diversity. No such association was found with binding antibodies directed to gp120. In summary, lower pretreatment viral diversity was associated with spontaneous control of viremia, reduced viral replication capacity and higher neutralizing antibody titers, suggesting a link between viral diversity, replication capacity, and neutralizing antibody activity.

Molecular study of the quasispecies evolution of a typical pigeon paramyxovirus type 1 after serial passages in pigeons by contact

The quasispecies nature of a typical pigeon paramyxovirus type 1 (pPMV-1) was, for the first time, studied under conditions close to the natural infectious environment. The virus was serially passaged in pigeons by successive contacts. Viral heterogeneity was analysed in the kidneys and brain of five pigeons from the last contact, by reverse transcriptase-polymerase chain reactions performed on RNA directly extracted from the organ and targeting the P and HN genes of the virus. The viral diversity following in vivo passage was found to be different from that in the inoculum, but demonstrated the reality of the quasispecies concept for pPMV-1 strains. Moreover, some aberrant genomic RNAs comprising insertions in the P gene editing site or deletions in the HN gene were also detected, with possible consequences for the pathogenicity and infectivity of the virus.

Hepatitis C virus core protein acts as a trans-modulating factor on internal translation initiation of the viral RNA

Translation initiation of hepatitis C virus (HCV) RNA occurs through an internal ribosome entry site (IRES) located at its 5' end. As a positive-stranded virus, HCV uses the genomic RNA template for translation and replication, but the transition between these two processes remains poorly understood. HCV core protein (HCV-C) has been proposed as a good candidate to modulate such a regulation. However, current data are still the subject of controversy in attributing any potential role in HCV translation to the HCV core protein. Here we demonstrate that HCV-C displays binding activities toward both HCV IRES and the 40 S ribosomal subunit by using centrifugation on sucrose gradients. To gain further insight into these interactions, we investigated the effect of exogenous addition of purified HCV-C on HCV IRES activity by using an in vitro reporter assay. We found that HCV IRES-mediated translation was specifically modulated by HCV-C provided in trans, in a dose-dependent manner, with up to a 5-fold stimulation of the IRES efficiency upon addition of low amounts of HCV-C, followed by a decrease at high doses. Interestingly, mutations within some domains of the IRES as well as the presence of an upstream reporter gene both lead to changes in the expected effects, consistent with the high dependence of HCV IRES function on its overall structure. Collectively, these results indicate that the HCV core protein is involved in a tight modulation of HCV translation initiation, depending on its concentration, and they suggest an important biological role of this protein in viral gene expression.

Complete gene scanning by temperature gradient capillary electrophoresis using the cystic fibrosis transmembrane conductance regulator gene as a model

Many inherited diseases involve large genes with many different mutations. Identifying a wide spectrum of mutations requires an efficient gene-scanning method. By differentiating thermodynamic stability and mobility of heteroduplexes from heterozygous samples, temperature gradient capillary electrophoresis (TGCE) was used to scan the entire coding region of the cystic fibrosis transmembrane conductance regulator gene. An initial panel (29 different mutations) showed 100% agreement between TGCE scanning and previously genotyped results for heterozygous samples. Different peak patterns were observed for single base substitutions and base insertions/deletions. Subsequently, 12 deidentified clinical samples genotyped as wild type for 32 mutations were scanned for the entire 27 exons. Results were 100% concordance with the bidirectional sequence analysis. Ten samples had nucleotide variations including a reported base insertion in intron 14b (2789 + 2insA) resulting in a possible mRNA splicing defect, and an unreported missense mutation in exon 20 (3991 G/A) with unknown clinical significance. This methodology does not require labeled primers or probes for detection and separation through a temperature gradient eliminates laborious temperature optimization required for other technologies. TGCE automation and high-throughput capability can be implemented in a clinical environment for mutation scanning with high sensitivity, thus reducing sequencing cost and effort.

Genetic heterogeneity of the NS3 protease gene in hepatitis C virus genotype 1 from untreated infected patients

NS3 protease is essential for hepatitis C Virus (HCV) replication, and is one of the most promising targets for specific anti-HCV therapy. Its natural polymorphism has not been studied at the quasispecies level. In the present work, the genetic heterogeneity of the NS3 protease gene was analyzed in 17 HCV genotype 1 (5 subtypes 1a and 12 subtypes 1b) samples collected from infected patients before anti-viral therapy. A total of 294 clones were sequenced. Although the protease NS3 is considered to be one of the less variable genes in the HCV genome, variability of both nucleotide and amino acid sequences was found. In variants belonging to 1a and 1b subtypes, 224 and 267 of 543 positions showed one or more nucleotide substitutions, respectively. Forty and 74 of the 181 NS3 amino acid positions showed at least one mutation in HCV-1a and HCV-1b isolates, respectively. Most substitutions were conservative. This substantial polymorphism of the NS3 protease produced by HCV-1a and HCV-1b suggests that, despite the numerous functional and structural constraints, the enzyme is sufficiently flexible to tolerate substitutions.

Multiple, linked human immunodeficiency virus type 1 drug resistance mutations in treatment-experienced patients are missed by standard genotype analysis

To investigate the extent to which drug resistance mutations are missed by standard genotyping methods, we analyzed the same plasma samples from 26 patients with suspected multidrug-resistant human immunodeficiency virus type 1 by using a newly developed single-genome sequencing technique and compared it to standard genotype analysis. Plasma samples were obtained from patients with prior exposure to at least two antiretroviral drug classes and who were on a failing antiretroviral regimen. Standard genotypes were obtained by reverse transcriptase (RT)-PCR and sequencing of the bulk PCR product. For single-genome sequencing, cDNA derived from plasma RNA was serially diluted to 1 copy per reaction, and a region encompassing p6, protease, and a portion of RT was amplified and sequenced. Sequences from 15 to 46 single viral genomes were obtained from each plasma sample. Drug resistance mutations identified by single-genome sequencing were not detected by standard genotype analysis in 24 of the 26 patients studied. Mutations present in less than 10% of single genomes were almost never detected in standard genotypes (1 of 86). Similarly, mutations present in 10 to 35% of single genomes were detected only 25% of the time in standard genotypes. For example, in one patient, 10 mutations identified by single-genome sequencing and conferring resistance to protease inhibitors (PIs), nucleoside analog reverse transcriptase inhibitors, and nonnucleoside reverse transcriptase inhibitors (NNRTIs) were not detected by standard genotyping methods. Each of these mutations was present in 5 to 20% of the 20 genomes analyzed; 15% of the genomes in this sample contained linked PI mutations, none of which were present in the standard genotype. In another patient sample, 33% of genomes contained five linked NNRTI resistance mutations, none of which were detected by standard genotype analysis. These findings illustrate the inadequacy of the standard genotype for detecting low-frequency drug resistance mutations. In addition to having greater sensitivity, single-genome sequencing identifies linked mutations that confer high-level drug resistance. Such linkage cannot be detected by standard genotype analysis.

Viral error catastrophe by mutagenic nucleosides

Riboviruses and retroviruses have the highest rates of mutations of any known organism. Increasing the mutation rate of these viruses could exceed the error threshold for viability of a viral population within a host. Recent experiments with mutagenic nucleoside analogs validate this new approach to treating infection of RNA viruses. Lethal mutagenesis with HIV-infected cells in culture has been documented and has been postulated to be the mechanism for treatment of hepatitis C with ribavirin. We consider the viral dynamics involved in the formation of a quasispecies, the choice of mutagenic nucleoside analogs, and the studies that have demonstrated the feasibility of lethal mutagenesis.

Human immunodeficiency virus type 1 (HIV-1) diversity at time of infection is not restricted to certain risk groups or specific HIV-1 subtypes

African women frequently acquire several genetically distinct human immunodeficiency virus type 1 (HIV-1) variants from a heterosexual partner, whereas the acquisition of multiple variants appears to be rare in men. To determine whether newly infected individuals in other risk groups acquire genetically diverse viruses, we examined the viral envelope sequences in plasma samples from 13 women and 4 men from the United States infected with subtype B viruses and 10 men from Kenya infected with non-subtype B viruses. HIV-1 envelope sequences differed by more than 2% in three U.S. women, one U.S. man, and one Kenyan man near the time of seroconversion. These findings suggest that early HIV-1 genetic diversity is not exclusive to women from Africa or to infection with any particular HIV-1 subtype.

Identification of unique hepatitis C virus quasispecies in the central nervous system and comparative analysis of internal translational efficiency of brain, liver, and serum variants

Reports of cerebral dysfunction in chronic hepatitis C virus (HCV) infection have led to the suggestion that HCV may infect the central nervous system (CNS). We used reverse transcription-PCR, cloning, and sequencing to define quasispecies for the HCV internal ribosomal entry site (IRES) and hypervariable region 1 (HVR1) in autopsy-derived brain, liver, lymph node, and serum samples. There was evidence of tissue compartmentalization of sequences in the brain in two patients, with between 24 and 55% of brain-derived IRES sequences absent from the serum, and significant phylogenetic and phenetic clustering of the brain and lymph node HVR1 sequences. The IRES initiates cap-independent translation of the viral polyprotein. Two unique brain-derived IRES mutations (C(204)-->A and G(243)-->A), which have previously been associated with lymphoid replication and altered translational efficiency in cell culture, were found in one patient. We used a dicistronic reporter vector to test whether brain-derived variants showed altered IRES-mediated translational efficiency, which might favor CNS infection. The translational efficiencies of the brain-derived IRES sequences were generally reduced compared to those of the master serum and liver sequences in rabbit reticulocyte cell lysates and two human cell lines, HuH7 (liver) and CHME3 (microglial). The C(204)-->A and G(243)-->A mutations showed preserved translational efficiency in HuH7 cells but reduced efficiency in CHME3 cells. Our data provide evidence that the CNS is a site of HCV replication, consistent with the recent demonstration of negative-strand HCV RNA in brain, and suggest that IRES polymorphisms may be important as a viral strategy of reduced translation to favor latency in the CNS.

Rapid RT-PCR amplification of full-length poliovirus genomes allows rapid discrimination between wild-type and recombinant vaccine-derived polioviruses

Poliomyelitis outbreaks in areas that were free for a long time of wild-type polioviruses have been reported. Characterization at nucleotide level of the causative agents showed that the isolated viruses were recombinant oral polio vaccine (OPV)-derived polioviruses. To allow rapid identification and detailed analysis of such recombinant polioviruses, a robust full-length reverse transcriptase-PCR (RT-PCR) was developed using SuperScript II (RT) and expand (PCR). Without extensive purification, it was possible to amplify and characterize the full-length genomes of all selected vaccine, wild-type, and recombinant vaccine-derived polioviruses within a week. Endonuclease nuclease analysis (SpeI) of the full-length amplicons allowed easy discrimination between recombinant and non-recombinant polioviruses. Furthermore, sequence analysis of cloned full-length amplicons of a recombinant vaccine-derived poliovirus strain showed that the quasi-species nature of a viral stock is preserved during the RT-PCR procedure. This robust and rapid RT-PCR method will allow rapid characterization of (recombinant) poliovirus strains in case of a local poliomyelitis outbreak, and will help to assess the risk of the appearance of such strains after wild-type poliovirus has been eradicated globally.

Do Taq-generated RT-PCR products from RNA viruses accurately reflect viral genetic heterogeneity?

Since the first report of genetically heterogeneous, or quasispecies, populations of RNA viruses, the genetic heterogeneity of the RNA genomes of major viral pathogens has been extensively studied. These studies aim to provide insights into the evolutionary pressures that act upon viruses, in order to define windows where anti-viral therapies will be most effective, to take prognostic values from viral genetic distributions at a given time, and to aid the development of novel therapeutic compounds that may tilt viral replication towards information loss. Many methodologies are employed to analyse genetic distributions of a virus in a given sample, but all involve the generation, and subsequent analysis, of the sequence information contained in a reverse-transcription-polymerase chain reaction (RT-PCR) product. Despite the fact that the aim of these RT-PCRs is to obtain sequence information from viral genomes, their application to this task is approached without adequate consideration of this end-goal. The establishment of an RT-PCR for a specific viral target genome generally proceeds in the same fashion as one would apply to establishing a PCR to determine the presence or absence of a specific target sequence in a given sample. However, it is becoming increasingly apparent that RT-PCR products generated by amplification with the ubiquitous thermostable DNA polymerase Taq, coupled with standard cloning and sequencing methodologies, has the potential to yield inaccurate and misleading data as pertains to the information content of populations of RNA viral genomes. This review discusses varying approaches employed to analyse heterogeneous populations of hepatitis C virus RNA genomes.

Quasispecies variation of porcine reproductive and respiratory syndrome virus during natural infection

Porcine reproductive and respiratory syndrome virus (PRRSV) displays notorious genetic, antigenic, and clinical variability. Little is known, however, about the nature and extent of viral variation present within naturally infected animals. By amplifying and cloning the open reading frame 5 gene from tonsils of naturally infected swine, and by sequencing individual clones, we characterized viral diversity in nine animals from two farms. All animals harbored multiple PRRSV variants at both the nucleic and the amino acid levels. Structural variation and rates of synonymous and nonsynonymous nucleotide substitution were no different within known epitopes than elsewhere. Analysis of molecular variance indicated that differences between farms, among animals within farms, and within individual animals accounted for 92.94, 3.84, and 3.22% of the total viral genetic variability observed, respectively. PRRSV exists during natural infection as a quasispecies distribution of related genotypes. Positive natural selection for immune evasiveness does not appear to maintain this diversity.