Search

Understanding the extent and limitation of viral genome evolution can provide insight about potential drug and vaccine targets. Influenza B Viruses (IBVs) infect humans in a seasonal manner and causes significant morbidity and mortality. IBVs are negative-sense single-stranded RNA viruses with a segmented genome and can be divided into two antigenically distinct lineages. The two lineages have been circulating and further evolving for almost four decades. The immune response to IBV infection can lead to antibodies that target the strain causing the infection. Some antibodies are cross-reactive and are able to bind strains from both lineages but, because of antigenic drift and immunodominance, both lineages continue to evolve and challenge human health. Here we investigate changes in the genomes of an IBVs from each lineage after passage in tissue culture in the presence of human sera containing polyclonal antibodies directed toward antigenically and temporally distinct viruses. Our previous analysis of the fourth segment, which encodes the major surface protein HA, revealed a pattern of change in which signature sequences from one lineage mutated to the signature sequences of the other lineage. Here we analyze genes from the other genomic segments and observe that most of the quasispecies’ heterogeneity occurs at the same loci in each lineage. The nature of the variants at these loci are investigated and possible reasons for this pattern are discussed. This work expands our understanding of the extent and limitations of genomic change in IBV.

Rubella viruses (RV) have been found in an association with granulomas in children with primary immune deficiencies (PID). Here, we report the recovery and characterization of infectious immunodeficiency-related vaccine-derived rubella viruses (iVDRV) from diagnostic skin biopsies of four patients. Sequence evolution within PID hosts was studied by comparison of the complete genomic sequences of the iVDRVs with the genome of the vaccine virus RA27/3. The degree of divergence of each iVDRV correlated with the duration of persistence indicating continuous intrahost evolution. The evolution rates for synonymous and nonsynonymous substitutions were estimated to be 5.7 x 10−3 subs/site/year and 8.9 x 10−4 subs/site/year, respectively. Mutational spectra and signatures indicated a major role for APOBEC cytidine deaminases and a secondary role for ADAR adenosine deaminases in generating diversity of iVDRVs. The distributions of mutations across the genes and 3D hotspots for amino acid substitutions in the E1 glycoprotein identified regions that may be under positive selective pressure. Quasispecies diversity was higher in granulomas than in recovered infectious iVDRVs. Growth properties of iVDRVs were assessed in WI-38 fibroblast cultures. None of the iVDRV isolates showed complete reversion to wild type phenotype but the replicative and persistence characteristics of iVDRVs were different from those of the RA27/3 vaccine strain, making predictions of iVDRV transmissibility and teratogenicity difficult. However, detection of iVDRV RNA in nasopharyngeal specimen and poor neutralization of some iVDRV strains by sera from vaccinated persons suggests possible public health risks associated with iVDRV carriers. Detection of IgM antibody to RV in sera of two out of three patients may be a marker of virus persistence, potentially useful for identifying patients with iVDRV before development of lesions. Studies of the evolutionary dynamics of iVDRV during persistence will contribute to development of infection control strategies and antiviral therapies., Author summary Primary immunodeficiency diseases (PID) are caused by genetic defects and lead to serious problems including chronic granulomas (abnormal collections (nodules) of inflammatory cells), sometimes lasting for decades and sometimes leading to severe ulcers. Initial reports (2014–2016), including our report of a blinded study using ultrasensitive virus detection in biopsies, proved the association between granuloma of the skin in PID patients and rubella virus. The viruses in these reports and the current report were derived from a widely used vaccine strain of the rubella virus. Work reported here shows that these vaccine-derived viruses are biologically different from the vaccine virus and that their genomes have changed. Genomic changes could be analyzed largely because the exact sequence of starting vaccine virus genome was known. These genomic differences are likely generated via mechanisms similar to those occurring during normal circulation of wild type rubella. We present data that newly recognized mechanisms for generation of sequence diversity in viruses (because of cellular deaminases) likely occurs in the generation of these vaccine-derived rubella viruses. Thousands of PID patients in the United States are likely shedding these vaccine-derived rubella viruses. Our work presented here characterizing viruses in diagnostic specimens highlights at least two areas where insufficient work has been done: 1) research on the properties of rubella virus (limited understanding of the antibody binding sites on the virus); 2) controlled research studies to assess the public health impact of viruses in populations with high immunity.

The impacts of invertebrate RNA virus population dynamics on virulence and infection outcomes are poorly understood. Deformed wing virus (DWV), the main viral pathogen of honey bees, negatively impacts bee health, which can lead to colony death. Despite previous reports on the reduction of DWV diversity following the arrival of the parasitic mite Varroa destructor, the key DWV vector, we found high genetic diversity of DWV in infested United States honey bee colonies. Phylogenetic analysis showed that divergent US DWV genotypes are of monophyletic origin and were likely generated as a result of diversification after a genetic bottleneck. To investigate the population dynamics of this divergent DWV, we designed a series of novel infectious cDNA clones corresponding to coexisting DWV genotypes, thereby devising a reverse-genetics system for an invertebrate RNA virus quasispecies. Equal replication rates were observed for all clone-derived DWV variants in single infections. Surprisingly, individual clones replicated to the same high levels as their mixtures and even the parental highly diverse natural DWV population, suggesting that complementation between genotypes was not required to replicate to high levels. Mixed clone–derived infections showed a lack of strong competitive exclusion, suggesting that the DWV genotypes were adapted to coexist. Mutational and recombination events were observed across clone progeny, providing new insights into the forces that drive and constrain virus diversification. Accordingly, our results suggest that Varroa influences DWV dynamics by causing an initial selective sweep, which is followed by virus diversification fueled by negative frequency-dependent selection for new genotypes. We suggest that this selection might reflect the ability of rare lineages to evade host defenses, specifically antiviral RNA interference (RNAi). In support of this hypothesis, we show that RNAi induced against one DWV strain is less effective against an alternate strain from the same population., Deformed wing virus, a key pathogen of honey bees, shows rapid diversification after genetic bottlenecks; a novel reverse-genetic system provides insights into the forces that shape virus diversity, suggesting that virus quasi-species diversification may be driven by selection of genotypes capable of evading host RNAi defences.

HIV-1 transmission is associated with a severe bottleneck in which a limited number of variants from a pool of genetically diverse quasispecies establishes infection. The IAVI protocol C cohort of discordant couples, female sex workers, other heterosexuals and men who have sex with men (MSM) present varying risks of HIV infection, diverse HIV-1 subtypes and represent a unique opportunity to characterize transmitted/founder viruses (TF) where disease outcome is known. To identify the TF, the HIV-1 repertoire of 38 MSM participants’ samples was sequenced close to transmission (median 21 days post infection, IQR 18–41) and assessment of multivariant infection done. Patient derived gag genes were cloned into an NL4.3 provirus to generate chimeric viruses which were characterized for replicative capacity (RC). Finally, an evaluation of how the TF virus predicted disease progression and modified the immune response at both acute and chronic HIV-1 infection was done. There was higher prevalence of multivariant infection compared with previously described heterosexual cohorts. A link was identified between multivariant infection and replicative capacity conferred by gag, whereby TF gag tended to be of lower replicative capacity in multivariant infection (p = 0.02) suggesting an overall lowering of fitness requirements during infection with multiple variants. Notwithstanding, multivariant infection was associated with rapid CD4+ T cell decline and perturbances in the CD4+ T cell and B cell compartments compared to single variant infection, which were reversible upon control of viremia. Strategies aimed at identifying and mitigating multivariant infection could contribute toward improving HIV-1 prognosis and this may involve strategies that tighten the stringency of the transmission bottleneck such as treatment of STI. Furthermore, the sequences and chimeric viruses help with TF based experimental vaccine immunogen design and can be used in functional assays to probe effective immune responses against TF., Author summary The development of a safe and effective HIV-1 vaccine alongside cure strategies is a major public health concern and requires in-depth knowledge of the HIV-1 populations that establish infection. Here we sequenced the HIV genetic populations present during the acute stage of HIV-1 infection in 38 Men who have sex with men (MSM) and identified the transmitted/founder variant’s sequence. We then cloned the gag gene from each patient’s transmitted/founder gag sequences in both single and multiple variant infection into a common, lab-adapted virus and measured the replicative capacity it conferred on this virus. Finally, cellular immune responses were compared between single variant and multiple variant infection at 0–3, 6–9 and 24–30 months after infection. We observed a higher frequency of multivariant infection in MSM than has been previously described in heterosexually infected participants, and this was associated with faster decline of CD4+ T cells and perturbances in the CD4+ T cell and the B cell compartment. Moreover, the replicative capacity conferred by gag was lower in multivariant infection, suggesting a less stringent transmission bottleneck that allowed for less fit variants to establish infection. These results suggest that strategies to tighten the stringency of the transmission bottleneck may be of benefit to patients by reducing the likelihood of multivariant transmission and potentially slowing down disease progression.

HIV-1 can downregulate HLA-C on infected cells, using the viral protein Vpu, and the magnitude of this downregulation varies widely between primary HIV-1 variants. The selection pressures that result in viral downregulation of HLA-C in some individuals, but preservation of surface HLA-C in others are not clear. To better understand viral immune evasion targeting HLA-C, we have characterized HLA-C downregulation by a range of primary HIV-1 viruses. 128 replication competent viral isolates from 19 individuals with effective anti-retroviral therapy, show that a substantial minority of individuals harbor latent reservoir virus which strongly downregulates HLA-C. Untreated infections display no change in HLA-C downregulation during the first 6 months of infection, but variation between viral quasispecies can be detected in chronic infection. Vpu molecules cloned from plasma of 195 treatment naive individuals in chronic infection demonstrate that downregulation of HLA-C adapts to host HLA genotype. HLA-C alleles differ in the pressure they exert for downregulation, and individuals with higher levels of HLA-C expression favor greater viral downregulation of HLA-C. Studies of primary and mutant molecules identify 5 residues in the transmembrane region of Vpu, and 4 residues in the transmembrane domain of HLA-C, which determine interactions between Vpu and HLA. The observed adaptation of Vpu-mediated downregulation to host genotype indicates that HLA-C alleles differ in likelihood of mediating a CTL response that is subverted by viral downregulation, and that preservation of HLA-C expression is favored in the absence of these responses. Finding that latent reservoir viruses can downregulate HLA-C could have implications for HIV-1 cure therapy approaches in some individuals.

Dengue virus affects millions of people worldwide each year. To date, there is no drug for the treatment of dengue-associated disease. Nucleosides are effective antivirals and work by inhibiting the accurate replication of the viral genome. Nucleobases offer a cheaper alternative to nucleosides for broad antiviral applications. Metabolic activation of nucleobases involves condensation with 5-phosphoribosyl-1-pyrophosphate to give the corresponding nucleoside-5’-monophosphate. This could provide an alternative to phosphorylation of a nucleoside, a step that is often rate limiting and inefficient in activation of nucleosides. We evaluated more than 30 nucleobases and corresponding nucleosides for their antiviral activity against dengue virus. Five nucleobases and two nucleosides were found to induce potent antiviral effects not previously described. Our studies further revealed that nucleobases were usually more active with a better tissue culture therapeutic index than their corresponding nucleosides. The development of viral lethal mutagenesis, an antiviral approach that takes into account the quasispecies behavior of RNA viruses, represents an exciting prospect not yet studied in the context of dengue replication. Passage of the virus in the presence of the nucleobase 3a (T-1105) and corresponding nucleoside 3b (T-1106), favipiravir derivatives, induced an increase in apparent mutations, indicating lethal mutagenesis as a possible antiviral mechanism. A more concerted and widespread screening of nucleobase libraries is a very promising approach to identify dengue virus inhibitors including those that may act as viral mutagens., Author summary Dengue virus is a world-wide public health menace estimated to infect hundreds of millions of people per year. Vaccines to prevent dengue virus infection have had limited success due in part to the requirement to elicit effective immune responses against the four dengue serotypes. There is an urgent unmet need for anti-dengue virus therapies. Nucleosides are effective antiviral small molecules which usually work by inhibiting the accurate replication of the viral genome. Typically, nucleosides must be converted within the cell to their triphosphate form to inhibit virus replication, thus inefficient phosphorylation often leads to suboptimal activity. We screened a small library of nucleobases that require an activation pathway different from nucleosides to achieve the same active form. We identified some known and previously undescribed dengue virus nucleobase inhibitors and their corresponding nucleosides. Our investigation of the mechanism of action of one nucleobase and its corresponding nucleoside found evidence for enhanced mutagenesis of the dengue virus genome in the presence of the compounds in cell culture. A wide screening of nucleobases libraries is a promising strategy to discover dengue virus inhibitors including potential viral mutagens.

Broadly neutralizing antibodies (bNAbs) are rarely elicited by current human immunodeficiency virus type 1 (HIV-1) vaccine designs, but the presence of bNAbs in naturally infected individuals may be associated with high plasma viral loads, suggesting that the magnitude, duration, and diversity of viral exposure may contribute to the development of bNAbs. Here, we report the isolation and characterization of a panel of human monoclonal antibodies (mAbs) from two subjects who developed broadly neutralizing autologous antibody responses during HIV-1 infection. In both subjects, we identified collections of mAbs that exhibited specificity only to a few autologous envelopes (Envs), with some mAbs exhibiting specificity only to a subset of Envs within the quasispecies of a particular sample at one time point. Neutralizing antibodies (NAbs) isolated from these subjects mapped mostly to epitopes in the Env V3 loop region and the CD4 binding site. None of the individual neutralizing mAbs recovered exhibited the cumulative breadth of neutralization present in the serum of the subjects. Surprisingly, however, the activity of polyclonal mixtures comprising individual mAbs that each possessed limited neutralizing activity, could achieve increased breadth of neutralizing activity against autologous isolates. While a single broadly neutralizing antibody targeting one epitope can mediate neutralization breadth, the findings presented here suggest that a cooperative polyclonal process mediated by diverse antibodies with more limited breadth targeting multiple epitopes also can achieve neutralization breadth against HIV-1.

Occult HBV infection (OBI), defined by the presence of HBV DNA in absence of hepatitis B surface antigen (HBsAg), is a significant concern in the HIV-infected population. Of 441 HIV+/HBsAg- patients analyzed, the overall prevalence of OBI was 6.3% (28/441). OBI was identified in 21 anti-HBc positives (17.8%), as well as among those who lacked any HBV-specific serological markers (2.2%). Comparison with HIV/HBV co-infection revealed that the levels of CD4, ALT, and HBV DNA were significantly lower during occult infection. Discrete differences were also observed with respect to quasispecies divergence. Additionally, subgenotype D1 was most frequent in occult infection, while D2 was widespread during chronic infection. The majority (~90%) of occult D1 sequences had the sQ129R mutation in the surface gene. This study highlights several distinct features of OBI in India and underscores the need for additional HBV DNA screening in HIV-positive individuals.

Purpose Hepatitis C virus (HCV) predominantly infects hepatocytes, although it is known that receptors for viral entry are distributed on a wide array of target cells. Chronic HCV infection is indeed characterized by multiple non-liver manifestations, suggesting a more complex HCV tropism extended to extrahepatic tissues and remains to be fully elucidated. In this study, we investigated the gastrointestinal mucosa (GIM) as a potential extrahepatic viral replication site and its contribution to HCV recurrence. Methods We analyzed GIM biopsies from a cohort of 76 patients, 11 of which were HCV-negative and 65 HCV-positive. Of these, 54 biopsies were from liver-transplanted patients. In 29 cases, we were able to investigate gastrointestinal biopsies from the same patient before and after transplant. To evaluate the presence of HCV, we looked for viral antigens and genome RNA, whilst to assess viral replicative activity, we searched for the replicative intermediate minus-strand RNA. We studied the genetic diversity and the phylogenetic relationship of HCV quasispecies from plasma, liver and gastrointestinal mucosa of HCV-liver-transplanted patients in order to assess HCV compartmentalization and possible contribution of gastrointestinal variants to liver re-infection after transplantation. Results Here we show that HCV infects and replicates in the cells of the GIM and that the favorite hosts were mostly enteroendocrine cells. Interestingly, we observed compartmentalization of the HCV quasispecies present in the gastrointestinal mucosa compared to other tissues of the same patient. Moreover, the phylogenetic analysis revealed a high similarity between HCV variants detected in gastrointestinal mucosa and those present in the re-infected graft. Conclusions Our results demonstrated that the gastrointestinal mucosa might be considered as an extrahepatic reservoir of HCV and that could contribute to viral recurrence. Moreover, the finding that HCV infects and replicates in neuroendocrine cells opens new perspectives on the role of these cells in the natural history of HCV infection.

Altres ajuts: CONACYT/GCPS/44519 Although Structured Treatment Interruptions (STI) are currently not considered an alternative strategy for antiretroviral treatment, their true benefits and limitations have not been fully established. Some studies suggest the possibility of improving the quality of life of patients with this strategy; however, the information that has been obtained corresponds mostly to studies conducted in adults, with a lack of knowledge about its impact on children. Furthermore, mutations associated with antiretroviral resistance could be selected due to sub-therapeutic levels of HAART at each interruption period. Genotyping methods to determine the resistance profiles of the infecting viruses have become increasingly important for the management of patients under STI, thus low-abundance antiretroviral drug-resistant mutations (DRM's) at levels under limit of detection of conventional genotyping (

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.

This study assesses the presence and outcome of genotype mixtures in the polymerase/surface and X/preCore regions of the HBV genome in patients with chronic hepatitis B virus (HBV) infection. Thirty samples from ten chronic hepatitis B patients were included. The polymerase/surface and X/preCore regions were analyzed by deep sequencing (UDPS) in the first available sample at diagnosis, a pre-treatment sample, and a sample while under treatment. HBV genotype was determined by phylogenesis. Quasispecies complexity was evaluated by mutation frequency and nucleotide diversity. The polymerase/surface and X/preCore regions were validated for genotyping from 113 GenBank reference sequences. UDPS yielded a median of 10,960 sequences per sample (IQR 16,645) in the polymerase/surface region and 11,595 sequences per sample (IQR 14,682) in X/preCore. Genotype mixtures were more common in X/preCore (90%) than in polymerase/surface (30%) (p

Made available in DSpace on 2018-12-11T17:26:54Z (GMT). No. of bitstreams: 0 Previous issue date: 2015-12-01 National Institutes of Health The hepatitis C virus (HCV) is a species of diverse genotypes that infect over 170 million people worldwide, causing chronic inflammation, cirrhosis and hepatocellular carcinoma. HCV genotype 3a is common in Brazil, and it is associated with a relatively poor response to current direct-acting antiviral therapies. The HCV NS3 protein cleaves part of the HCV polyprotein, and cellular antiviral proteins. It is therefore the target of several HCV drugs. In addition to its protease activity, NS3 is also an RNA helicase. Previously, HCV present in a relapse patient was found to harbor a mutation known to be lethal to HCV genotype 1b. The point mutation encodes the amino acid substitution W501R in the helicase RNA binding site. To examine how the W501R substitution affects NS3 helicase activity in a genotype 3a background, wild type and W501R genotype 3a NS3 alleles were sub-cloned, expressed in E. coli, and the recombinant proteins were purified and characterized. The impact of the W501R allele on genotype 2a and 3a subgenomic replicons was also analyzed. Assays monitoring helicase-catalyzed DNA and RNA unwinding revealed that the catalytic efficiency of wild type genotype 3a NS3 helicase was more than 600 times greater than the W501R protein. Other assays revealed that the W501R protein bound DNA less than 2 times weaker than wild type, and both proteins hydrolyzed ATP at similar rates. In Huh7.5 cells, both genotype 2a and 3a subgenomic HCV replicons harboring the W501R allele showed a severe defect in replication. Since the W501R allele is carried as a minor variant, its replication would therefore need to be attributed to the trans-complementation by other wild type quasispecies. São Paulo State University - UNESP Department of Biology University of Wisconsin- Milwaukee Department of Chemistry and Biochemistry São José do Rio Preto Medical School Laboratory of Virology São Paulo State University - UNESP Department of Biology National Institutes of Health: R01 AI088001

Background Aim was to determine the dynamics of peripheral blood mononuclear cells (PBMC)- associated total HIV-1 DNA in successfully ART-treated HIV/HCV co-infected patients receiving DAA treatment and to explore possible virological hypotheses underlying the phenomenon. Methods Longitudinal, single-centre study measuring total HIV-1 DNA before the start of DAA, at the end of treatment (EOT), and 3 months after treatment. Univariable and multivariable analyses were used to assess factors associated with HIV-1 DNA increase ≥0.5 Log copies/million PBMC. Episomal 2-LTR forms, residual HIV-1 viremia and proviral DNA quasispecies evolution were also investigated. Results 119 successfully ART-treated HIV/HCV co-infected patients were included. Median baseline HIV-1 DNA was 3.84 Log copies/million PBMC (95%CI 3.49–4.05), and no significant variation with respect to baseline was found at EOT and after 3 months of DAA termination. In 17% of cases an increase ≥0.5 Log copies/million PBMC was observed at EOT compared to baseline. HIV-1 DNA increase was independently associated with lower baseline HIV-1 DNA, longer HIV suppression, raltegravir-based ART and previous exposure to interferon/ribavirin for HCV treatment. In none of the patients with HIV-1 DNA increase, 2-LTR forms were detected at baseline, while in 2 cases 2-LTR forms were found at EOT, without association with residual HIV-1 RNA viremia. No evidence of viral evolution was observed. Conclusions In successfully ART-treated HIV/HCV co-infected patients receiving DAA, PBMC-associated total HIV-1 DNA was quite stable over time, but some patients showed a considerable increase at EOT when compared to baseline. A significantly higher risk of HIV DNA increase was found, in presence of lower cellular HIV reservoir at baseline. Activation of replicative-competent virus generating new rounds of viral replication seems unlikely, while mobilization of cell-associated HIV from tissue reservoirs could be hypothesized.

Background Hepatitis C virus (HCV) infection is a global health problem estimated to affect almost 200 million people worldwide. The aim of this study is to analyze the subtypes and existence of variants resistant to protease inhibitors and their association with potential HCV risk factors among blood donors in Brazil. Methods Repeat anti-HCV reactive blood donors are systematically asked to return for retest, notification, and counseling in which they are interviewed for risk factors for transfusion-transmitted diseases. We analyzed 202 donors who returned for counseling from 2007 to 2010 and presented enzyme immunoassay- and immunoblot-reactive results. The HCV genotypes and resistance mutation analyses were determined by the direct sequencing of the NS5b and NS3 regions, respectively. The HCV viral load was determined using an in-house real-time PCR assay targeting the 5′-NCR. Results HCV subtypes 1b, 1a, and 3a were found in 45.5%, 32.0%, and 18.0% of the donors, respectively. The mean viral load of genotype 1 was significantly higher than that of the genotype 3 isolates. Subtype 1a was more frequent among young donors and 3a was more frequent among older donors. Protease inhibitor-resistant variants were detected in 12.8% of the sequenced samples belonging to genotype 1, and a higher frequency was observed among subtype 1a (20%) in comparison to 1b (8%). There was no difference in the prevalence of HCV risk factors among the genotypes or drug-resistant variants. Conclusions We found a predominance of subtype 1b, with an increase in the frequency of subtype 1a, in young subjects. Mutations conferring resistance to NS3 inhibitors were frequent in treatment-naive blood donors, particularly those infected with subtype 1a. These variants were detected in the major viral population of HCV quasispecies, have replicative capacities comparable to nonresistant strains, and could be important for predicting the response to antiviral triple therapy.

Objective: To assess the presence of two major non-nucleoside reverse transcriptase inhibitors (NNRTI) drug resistance mutations (DRMs), Y181C and K103N, in minor viral quasispecies of treatment naive HIV-1 infected East-African and Swedish patients by allele-specific polymerase chain reaction (AS-PCR).Methods: Treatment naive adults (n = 191) with three epidemiological backgrounds were included: 92 Ethiopians living in Ethiopia; 55 East-Africans who had migrated to Sweden; and 44 Caucasians living in Sweden. The pol gene was analysed by standard population sequencing and by AS-PCR for the detection of Y181C and K103N.Results: The Y181C was detected in the minority quasispecies of six Ethiopians (6.5%), in two Caucasians (4.5%), and in one East-African (1.8%). The K103N was detected in one East-African (1.8%), by both methods. The proportion of mutants ranged from 0.25% to 17.5%. Additional DRMs were found in all three treatment naive patient groups by population sequencing.Conclusions: Major NNRTI mutations can be found by AS-PCR in minor quasispecies of treatment naive HIV-1 infected Ethiopians living in Ethiopia, in East-African and Caucasian patients living in Sweden in whom population sequencing reveal wild-type virus only. Surveys with standard sequencing are likely to underestimate transmitted drug resistance and the presence of resistant minor quasispecies in treatment naive patients should be topic for future large scale studies.

It was brought to the attention of PLOS ONE staff editors that there are similarities in study design and some findings between this article and the following previously published work: Zhang M, Foley B, Schultz AK, Macke JP, Bulla I, Stanke M, Morgenstern B, Korber B, Leitner T: The role of recombination in the emergence of a complex and dynamic HIV epidemic. Retrovirology 2010 Mar, 7:25. doi:10.1186/1742-4690-7-25 [2] This previously published study used the jpHMM program to analyze sequences, including BC recombinants CRF07 and CRF08 and BF recombinants, and thus reported work related to the findings described in the current publication. Jia et al. have used the same program, jpHMM, as well as RDP3, to analyze HIV sequences. The authors wish to clarify that the current article was focused on describing a more accurate recombination pattern of the Reference strains of 6 selected circulating recombinant forms (CRFs) compared to those previously reported, while Zhang et al. reported whether full-length and sequence fragments (but not the Reference strains) of CRFs required subtype assignment. Jia et al. also discuss a comparison between jpHMM-derived recombination structures of Reference strains of the CRFs with the original ones listed on the HIV Circulating Recombinant Forms page of the Los Alamos database (http://www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html). In addition, the following paragraph from the Methods section was duplicated almost verbatim from the above article [2]: “In jpHMM, each HIV-1 subtype is represented by a profile hidden Markov model. All profile models are connected by empirical probabilities, allowing the detection of possible recombinants and related breakpoints by jumping from one profile to another. jpHMM performs best in predicting recombinants that involve subtypes that have had adequate sampling to build well-informed profiles. It is less effective in cases related to subtypes H, J, and K because so few full-length genome sequences are available (N = 4, 3, and 2, respectively). In the present study, jpHMM was used to detect the recombination patterns in recombinants composed exclusively of subtypes B, C, F, and G. Each of these subtypes has enough data to form a good model of sequence variation.” The authors wished to ensure an accurate description of the jpHMM tool by using the description previously published by the developers of the tool; however, the authors recognize the sentences should have been reproduced in quotation marks. Figure 2B of the published article shows some similar findings to those presented in the following thesis, which was not cited: Doctoral Thesis Ming Zhang: Tracking HIV-1 Genetic Variation https://www.uni-goettingen.de/en/zhang-ming-tracking-hiv-1-genetic-variation/73586.html [3] Similarly to the work on the thesis, Figure 2B in this article displays sub-regions trees; however, the authors would like to clarify that the goal is to confirm the subtype adscription of the fragment in nef region of CRF08, while the thesis reports trees constructed for sub-regions delimited by all CRF07 and CRF08 sequences in order to investigate the relationship between CRF07 and CRF08. Additional instances of duplicated text were also identified within the following excerpts from the Introduction and Methods sections: “RDP3 treats every sequence within the analyzed alignment as a potential recombinant and systematically screens sequence triplets or quartets to identify viruses that contain a recombinant and two sequences that could serve as parents while performing a statistical evaluation of recombination signals … Such an approach eliminates the need for reference sequences, which makes analysis of viral quasispecies from epidemiologically unlinked patients more practical.” The above text is similar to that found in “Frequent Intra-Subtype Recombination among HIV-1 Circulating in Tanzania” PLOS ONE DOI: 10.1371/journal.pone.0071131. [4] "as the prototype of subtype I based on gp120 sequences … Full genome sequencing revealed this virus to be a complex mosaic with multiple breakpoints between regions of several distinct subtypes, including A, G, and I … Subsequent analysis with previously unavailable complete genome sequences has revealed that the virus is in fact mosaic with regions associated with subtypes A, G, H, and K and unclassified regions." The above text is similar to that found in “CRF06-cpx: A New Circulating Recombinant Form of HIV-1 in West Africa Involving Subtypes A, G, K, and J.” JAIDS Journal of Acquired Immune Deficiency Syndromes 29:522–530. [5] "This prediction method is based on a precalculated multiple alignment of the major HIV-1 subtypes including CRF01_AE references, and it is more accurate than the competing methods used for phylogenetic breakpoint detection." The above text is similar to that found in “Identification and Genetic Characterization of a Novel CRF22_01A1 Recombinant Form of HIV Type 1 in Cameroon” DOI:10.1089/aid.2009.0197. [6] The authors would like to apologize for all instances of overlapping text.

One of the hurdles to understanding the role of viral quasispecies in RNA virus cross-species transmission (CST) events is the need to analyze a densely sampled outbreak using deep sequencing in order to measure the amount of mutation occurring on a small time scale. In 2009, the California Department of Public Health reported a dramatic increase (350) in the number of gray foxes infected with a rabies virus variant for which striped skunks serve as a reservoir host in Humboldt County. To better understand the evolution of rabies, deep-sequencing was applied to 40 unpassaged rabies virus samples from the Humboldt outbreak. For each sample, approximately 11 kb of the 12 kb genome was amplified and sequenced using the Illumina platform. Average coverage was 17,448 and this allowed characterization of the rabies virus population present in each sample at unprecedented depths. Phylogenetic analysis of the consensus sequence data demonstrated that samples clustered according to date (1995 vs. 2009) and geographic location (northern vs. southern). A single amino acid change in the G protein distinguished a subset of northern foxes from a haplotype present in both foxes and skunks, suggesting this mutation may have played a role in the observed increased transmission among foxes in this region. Deep-sequencing data indicated that many genetic changes associated with the CST event occurred prior to 2009 since several nonsynonymous mutations that were present in the consensus sequences of skunk and fox rabies samples obtained from 20032010 were present at the sub-consensus level (as rare variants in the viral population) in skunk and fox samples from 1995. These results suggest that analysis of rare variants within a viral population may yield clues to ancestral genomes and identify rare variants that have the potential to be selected for if environment conditions change., Author Summary Understanding the role of genetic variants within a viral population is a necessary step toward predicting and treating emerging infectious diseases. The high mutation rate of RNA viruses increases the ability of these viruses to adapt to diverse hosts and cause new human and zoonotic diseases. The genetic diversity of a viral population within a host may allow the virus to adapt to a diverse array of selective pressures and enable cross-species transmission events. In 2009 a large outbreak of rabies in Northern California involved a skunk rabies virus variant that efficiently transmitted within a population of gray foxes, suggesting possible adaptation to a novel host species. To better understand the evolution of rabies virus that enabled this host jump, we applied deep-sequencing analysis to rabies virus samples from the outbreak. Deep-sequencing data indicated that many of the genetic changes associated with host jump occurred prior to 2009, and these mutations were present at very low frequencies in viral populations from samples dating back to 1995. These results suggest deep sequencing is useful for characterization of viral populations, and may provide insight to ancestral genomes and role of rare variants in viral emergence.

Uncultivable HPR0 strains of infectious salmon anaemia viruses (ISAVs) infecting gills are non-virulent putative precursors of virulent ISAVs (vISAVs) causing systemic disease in farmed Atlantic salmon (Salmo salar). The transition to virulence involves two molecular events, a deletion in the highly polymorphic region (HPR) of the hemagglutinin-esterase (HE) gene and a Q266→L266 substitution or insertion next to the putative cleavage site (R267) in the fusion protein (F). We have performed ultra-deep pyrosequencing (UDPS) of these gene regions from healthy fish positive for HPR0 virus carrying full-length HPR sampled in a screening program, and a vISAV strain from an ISA outbreak at the same farming site three weeks later, and compared the mutant spectra. As the UDPS data shows the presence of both HE genotypes at both sampling times, and the outbreak strain was unlikely to be directly related to the HPR0 strain, this is the first report of a double infection with HPR0s and vISAVs. For F amplicon reads, mutation frequencies generating L266 codons in screening samples and Q266 codons in outbreak samples were not higher than at any random site. We suggest quasispecies heterogeneity as well as RNA structural properties are linked to transition to virulence. More specifically, a mechanism where selected single point mutations in the full-length HPR alter the RNA structure facilitating single- or sequential deletions in this region is proposed. The data provides stronger support for the deletion hypothesis, as opposed to recombination, as the responsible mechanism for generating the sequence deletions in HE.

Following long-term infection with virus derived from the pathogenic GL8 molecular clone of feline immunodeficiency virus (FIV), a range of viral variants emerged with distinct modes of interaction with the viral receptors CD134 and CXCR4, and sensitivities to neutralizing antibodies. In order to assess whether this viral diversity would be maintained following subsequent transmission, a synthetic quasispecies was reconstituted comprising molecular clones bearing envs from six viral variants and its replicative capacity compared in vivo with a clonal preparation of the parent virus. Infection with either clonal (Group 1) or diverse (Group 2) challenge viruses, resulted in a reduction in CD4+ lymphocytes and an increase in CD8+ lymphocytes. Proviral loads were similar in both study groups, peaking by 10 weeks post-infection, a higher plateau (set-point) being achieved and maintained in study Group 1. Marked differences in the ability of individual viral variants to replicate were noted in Group 2; those most similar to GL8 achieved higher viral loads while variants such as the chimaeras bearing the B14 and B28 Envs grew less well. The defective replication of these variants was not due to suppression by the humoral immune response as virus neutralising antibodies were not elicited within the study period. Similarly, although potent cellular immune responses were detected against determinants in Env, no qualitative differences were revealed between animals infected with either the clonal or the diverse inocula. However, in vitro studies indicated that the reduced replicative capacity of variants B14 and B28 in vivo was associated with altered interactions between the viruses and the viral receptor and co-receptor. The data suggest that viral variants with GL8-like characteristics have an early, replicative advantage and should provide the focus for future vaccine development.

Hepatitis A virus (HAV) is the most common cause of infectious hepatitis throughout the world, spread largely by the fecal-oral route. To characterize the genetic diversity of the virus circulating in China where HAV in endemic, we selected the outbreak cases with identical sequences in VP1-2A junction region and compiled a panel of 42 isolates. The VP3-VP1-2A regions of the HAV capsid-coding genes were further sequenced and analyzed. The quasispecies distribution was evaluated by cloning the VP3 and VP1-2A genes in three clinical samples. Phylogenetic analysis demonstrated that the same genotyping results could be obtained whether using the complete VP3, VP1, or partial VP1-2A genes for analysis in this study, although some differences did exist. Most isolates clustered in sub-genotype IA, and fewer in sub-genotype IB. No amino acid mutations were found at the published neutralizing epitope sites, however, several unique amino acid substitutions in the VP3 or VP1 region were identified, with two amino acid variants closely located to the immunodominant site. Quasispecies analysis showed the mutation frequencies were in the range of 7.22 x 10(-4) -2.33 x 10(-3) substitutions per nucleotide for VP3, VP1, or VP1-2A. When compared with the consensus sequences, mutated nucleotide sites represented the minority of all the analyzed sequences sites. HAV replicated as a complex distribution of closely genetically related variants referred to as quasispecies, and were under negative selection. The results indicate that diverse HAV strains and quasispecies inside the viral populations are presented in China, with unique amino acid substitutions detected close to the immunodominant site, and that the possibility of antigenic escaping mutants cannot be ruled out and needs to be further analyzed.

Triple therapy of chronic hepatitis C virus (HCV) infection with boceprevir (BOC) or telaprevir (TVR) leads to virologic failure in many patients which is often associated with the selection of resistance-associated variants (RAVs). These resistance profiles are of importance for the selection of potential rescue treatment options. In this study, we sequenced baseline NS3 RAVs population-based and investigated the sensitivity of NS3 phenotypes in an HCV replicon assay together with clinical factors for a prediction of treatment response in a cohort of 165 German and Swiss patients treated with a BOC or TVR-based triple therapy. Overall, the prevalence of baseline RAVs was low, although the frequency of RAVs was higher in patients with virologic failure compared to those who achieved a sustained virologic response (SVR) (7% versus 1%, P = 0.06). The occurrence of RAVs was associated with a resistant NS3 quasispecies phenotype (P

The use of mutagenic drugs to drive HIV-1 past its error threshold presents a novel intervention strategy, as suggested by the quasispecies theory, that may be less susceptible to failure via viral mutation-induced emergence of drug resistance than current strategies. The error threshold of HIV-1, , however, is not known. Application of the quasispecies theory to determine poses significant challenges: Whereas the quasispecies theory considers the asexual reproduction of an infinitely large population of haploid individuals, HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We performed population genetics-based stochastic simulations of the within-host evolution of HIV-1 and estimated the structure of the HIV-1 quasispecies and . We found that with small mutation rates, the quasispecies was dominated by genomes with few mutations. Upon increasing the mutation rate, a sharp error catastrophe occurred where the quasispecies became delocalized in sequence space. Using parameter values that quantitatively captured data of viral diversification in HIV-1 patients, we estimated to be substitutions/site/replication, ∼2–6 fold higher than the natural mutation rate of HIV-1, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs. The latter estimate was weakly dependent on the within-host effective population size of HIV-1. With large population sizes and in the absence of recombination, our simulations converged to the quasispecies theory, bridging the gap between quasispecies theory and population genetics-based approaches to describing HIV-1 evolution. Further, increased with the recombination rate, rendering HIV-1 less susceptible to error catastrophe, thus elucidating an added benefit of recombination to HIV-1. Our estimate of may serve as a quantitative guideline for the use of mutagenic drugs against HIV-1., Author Summary Currently available antiretroviral drugs curtail HIV infection but fail to eradicate the virus. A strategy of intervention radically different from that employed by current drugs has been proposed by the molecular quasispecies theory. The theory predicts that increasing the viral mutation rate beyond a critical value, called the error threshold, would cause a severe loss of genetic information, potentially leading to viral clearance. Several chemical mutagens are now being developed that can increase the mutation rate of HIV-1. Their success depends on reliable estimates of the error threshold of HIV-1, which are currently lacking. The quasispecies theory cannot be applied directly to HIV-1: the theory considers an infinitely large population of asexually reproducing haploid individuals, whereas HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We employed detailed stochastic simulations that overcome the limitations of the quasispecies theory and accurately mimic HIV-1 evolution in vivo. With these simulations, we estimated the error threshold of HIV-1 to be ∼2–6-fold higher than its natural mutation rate, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs.

To address whether sequences of viral gag and env quasispecies collected during the early post-acute period can be utilized to determine multiplicity of transmitted HIV's, recently developed approaches for analysis of viral evolution in acute HIV-1 infection [1], [2] were applied. Specifically, phylogenetic reconstruction, inter- and intra-patient distribution of maximum and mean genetic distances, analysis of Poisson fitness, shape of highlighter plots, recombination analysis, and estimation of time to the most recent common ancestor (tMRCA) were utilized for resolving multiplicity of HIV-1 transmission in a set of viral quasispecies collected within 50 days post-seroconversion (p/s) in 25 HIV-infected individuals with estimated time of seroconversion. The decision on multiplicity of HIV infection was made based on the model's fit with, or failure to explain, the observed extent of viral sequence heterogeneity. The initial analysis was based on phylogeny, inter-patient distribution of maximum and mean distances, and Poisson fitness, and was able to resolve multiplicity of HIV transmission in 20 of 25 (80%) cases. Additional analysis involved distribution of individual viral distances, highlighter plots, recombination analysis, and estimation of tMRCA, and resolved 4 of the 5 remaining cases. Overall, transmission of a single viral variant was identified in 16 of 25 (64%) cases, and transmission of multiple variants was evident in 8 of 25 (32%) cases. In one case multiplicity of HIV-1 transmission could not be determined. In primary HIV-1 subtype C infection, samples collected within 50 days p/s and analyzed by a single-genome amplification/sequencing technique can provide reliable identification of transmission multiplicity in 24 of 25 (96%) cases. Observed transmission frequency of a single viral variant and multiple viral variants were within the ranges of 64% to 68%, and 32% to 36%, respectively.

A large number of medically important viruses, including HIV, hepatitis C virus, and influenza, have RNA genomes. These viruses replicate with extremely high mutation rates and exhibit significant genetic diversity. This diversity allows a viral population to rapidly adapt to dynamic environments and evolve resistance to vaccines and antiviral drugs. For the last 30 years, quasispecies theory has provided a population-based framework for understanding RNA viral evolution. A quasispecies is a cloud of diverse variants that are genetically linked through mutation, interact cooperatively on a functional level, and collectively contribute to the characteristics of the population. Many predictions of quasispecies theory run counter to traditional views of microbial behavior and evolution and have profound implications for our understanding of viral disease. Here, we discuss basic principles of quasispecies theory and describe its relevance for our understanding of viral fitness, virulence, and antiviral therapeutic strategy.

Background Mother-to-child transmission of HIV-1 remains a significant problem in the resource-constrained settings where anti-retroviral therapy is still not widely available. Understanding the earliest events during HIV-1 transmission and characterizing the newly transmitted or founder virus is central to intervention efforts. In this study, we analyzed the viral env quasispecies of six mother-infant transmission pairs (MIPs) and characterized the genetic features of envelope glycoprotein that could influence HIV-1 subtype C perinatal transmission. Methodology and Findings The V1-V5 region of env was amplified from 6 MIPs baseline samples and 334 DNA sequences in total were analyzed. A comparison of the viral population derived from the mother and infant revealed a severe genetic bottleneck occurring during perinatal transmission, which was characterized by low sequence diversity in the infant. Phylogenetic analysis indicates that most likely in all our infant subjects a single founder virus was responsible for establishing infection. Furthermore, the newly transmitted viruses from the infant had significantly fewer potential N-linked glycosylation sites in Env V1-V5 region and showed a propensity to encode shorter variable loops compared to the nontransmitted viruses. In addition, a similar intensity of selection was seen between mothers and infants with a higher rate of synonymous (dS) compared to nonsynonymous (dN) substitutions evident (dN/dS

Background The rapidly expanding availability of de novo sequencing technologies can greatly facilitate efforts to monitor the relatively high mutation rates of influenza A viruses and the detection of quasispecies. Both the mutation rates and the lineages of influenza A viruses are likely to play an important role in the natural history of these viruses and the emergence of phenotypically and antigenically distinct strains. Methodology and Principal Findings We evaluated quasispecies and mixed infections by de novo sequencing the whole genomes of 10 virus isolates, including eight avian influenza viruses grown in embryonated chicken eggs (six waterfowl isolates - five H3N2 and one H4N6; an H7N3 turkey isolate; and a bald eagle isolate with H1N1/H2N1 mixed infection), and two tissue cultured H3N2 swine influenza viruses. Two waterfowl cloacal swabs were included in the analysis. Full-length sequences of all segments were obtained with 20 to 787-X coverage for the ten viruses and one cloacal swab. The second cloacal swab yielded 15 influenza reads of ∼230 bases, sufficient for bioinformatic inference of mixed infections or quasispecies. Genomic subpopulations or quasispecies of viruses were identified in four egg grown avian influenza isolates and one cell cultured swine virus. A bald eagle isolate and the second cloacal swab showed evidence of mixed infections with two (H1 and H2) and three (H1, H3, and H4) HA subtypes, respectively. Multiple sequence differences were identified between cloacal swab and the virus recovered using embryonated chicken eggs. Conclusions We describe a new approach to comprehensively identify mixed infections and quasispecies in low passage influenza A isolates and cloacal swabs and add to the understanding of the ecology of influenza A virus populations.

The HIV-1 epidemic in sub-Saharan Africa is driven largely by heterosexual transmission of non-subtype B viruses, of which subtypes C and A are predominant. Previous studies of subtype B and subtype C transmission pairs have suggested that a single variant from the chronically infected partner can establish infection in their newly infected partner. However, in subtype A infected individuals from a sex worker cohort and subtype B individuals from STD clinics, infection was frequently established by multiple variants. This study examined over 1750 single-genome amplified viral sequences derived from epidemiologically linked subtype C and subtype A transmission pairs very early after infection. In 90% (18/20) of the pairs, HIV-1 infection is initiated by a single viral variant that is derived from the quasispecies of the transmitting partner. In addition, the virus initiating infection in individuals who were infected by someone other than their spouse was characterized to determine if genital infections mitigated the severe genetic bottleneck observed in a majority of epidemiologically linked heterosexual HIV-1 transmission events. In nearly 50% (3/7) of individuals infected by someone other than their spouse, multiple genetic variants from a single individual established infection. A statistically significant association was observed between infection by multiple genetic variants and an inflammatory genital infection in the newly infected individual. Thus, in the vast majority of HIV-1 transmission events in cohabiting heterosexual couples, a single genetic variant establishes infection. Nevertheless, this severe genetic bottleneck can be mitigated by the presence of inflammatory genital infections in the at risk partner, suggesting that this restriction on genetic diversity is imposed in large part by the mucosal barrier., Author Summary Previous studies of HIV transmission have yielded conflicting results regarding the genetic heterogeneity of the virus establishing infection in the newly infected individual. In this study of populations from Zambia and Rwanda that are infected by two distinct viral genetic subtypes, we compared viral sequences that encode the entry-mediating envelope glycoproteins from newly infected individuals (recipients) and their spouses (donors) very early after infection, as well as newly infected individuals infected by someone other than their spouse. In spite of the genetically diverse virus population in the donor, approximately 90% of newly infected individuals were infected by a single viral variant, while the rest were infected by multiple viral variants. The homogeneity of the virus population in the newly infected recipient, as well as the presence, in some cases, of identical virus variants in the donor, allowed us to precisely identify the transmitted variant. We were able to examine the clinical history of each newly infected individual and observed that all individuals infected by multiple variants also showed evidence of inflammatory genital infections. Our results suggest that the genital mucosa provides a natural barrier to infection by multiple genetic variants of HIV-1, but that this barrier can be lowered by inflammatory genital infections.

Mechanisms by which hepatitis C virus (HCV) evades cellular immunity to establish persistence in chronically infected individuals are not clear. Mutations in human leukocyte antigen (HLA) class I-restricted epitopes targeted by CD8+ T cells are associated with persistence, but the extent to which these mutations affect viral fitness is not fully understood. Previous work showed that the HCV quasispecies in a persistently infected chimpanzee accumulated multiple mutations in numerous class I epitopes over a period of 7 years. During the acute phase of infection, one representative epitope in the C-terminal region of the NS3/4A helicase, NS31629-1637, displayed multiple serial amino acid substitutions in major histocompatibility complex (MHC) anchor and T cell receptor (TCR) contact residues. Only one of these amino acid substitutions at position 9 (P9) of the epitope was stable in the quasispecies. We therefore assessed the effect of each mutation observed during in vivo infection on viral fitness and T cell responses using an HCV subgenomic replicon system and a recently developed in vitro infectious virus cell culture model. Mutation of a position 7 (P7) TCR-contact residue, I1635T, expectedly ablated the T cell response without affecting viral RNA replication or virion production. In contrast, two mutations at the P9 MHC-anchor residue abrogated antigen-specific T cell responses, but additionally decreased viral RNA replication and virion production. The first escape mutation, L1637P, detected in vivo only transiently at 3 mo after infection, decreased viral production, and reverted to the parental sequence in vitro. The second P9 variant, L1637S, which was stable in vivo through 7 years of follow-up, evaded the antigen-specific T cell response and did not revert in vitro despite being less optimal in virion production compared to the parental virus. These studies suggest that HCV escape mutants emerging early in infection are not necessarily stable, but are eventually replaced with variants that achieve a balance between immune evasion and fitness for replication., Author Summary Hepatitis C virus (HCV)-associated liver disease is a leading indication for liver transplantation in the United States. With more than 170 million people infected with HCV worldwide and more than 70% of those infected unable to clear the virus, it is of paramount importance to elucidate the factors leading to viral persistence. A well-characterized experimental chimpanzee chronically infected with HCV was found to develop multiple viral sequence variations over the course of 7 years. Serial mutations in an HCV epitope that was originally immunogenic in the host were observed during the course of infection, and we sought to better understand how each mutation affected viral persistence. We recreated the viral variants detected in the chimpanzee and assessed the ability of each variant to replicate, produce progeny virus, and evade the host immune system. We found that certain HCV variants present long after initial infection are able to avoid recognition by host immune cells, but have reduced replication fitness. Importantly, in the absence of immune pressure, these viruses mutate back to variants better able to replicate and produce new viral progeny. Our experiments suggest that rapidly mutating viruses like HCV seek a balance between replication efficiency and the ability to evade the host immune system.

CD8+ cytotoxic T lymphocytes (CTL) can be effective at controlling HIV-1 in humans and SIV in macaques, but their utility is partly offset by mutational escape. The kinetics of CTL escape and reversion of escape mutant viruses upon transmission to MHC-mismatched hosts can help us understand CTL-mediated viral control and the fitness cost extracted by immune escape mutation. Traditional methods for following CTL escape and reversion are, however, insensitive to minor viral quasispecies. We developed sensitive quantitative real-time PCR assays to track the viral load of SIV Gag164–172 KP9 wild-type (WT) and escape mutant (EM) variants in pigtail macaques. Rapid outgrowth of EM virus occurs during the first few weeks of infection. However, the rate of escape plateaued soon after, revealing a prolonged persistence of WT viremia not detectable by standard cloning and sequencing methods. The rate of escape of KP9 correlated with levels of vaccine-primed KP9-specific CD8+ T cells present at that time. Similarly, when non-KP9 responder (lacking the restricting Mane-A*10 allele) macaques were infected with SHIVmn229 stock containing a mixture of EM and WT virus, rapid reversion to WT was observed over the first 2 weeks following infection. However, the rate of reversion to WT slowed dramatically over the first month of infection. The serial quantitation of escape mutant viruses evolving during SIV infection shows that rapid dynamics of immune escape and reversion can be observed in early infection, particularly when CD8 T cells are primed by vaccination. However, these early rapid rates of escape and reversion are transient and followed by a significant slowing in these rates later during infection, highlighting that the rate of escape is significantly influenced by the timing of its occurrence., Author Summary Immune escape from AIDS virus–specific cellular immunity is common. The driving forces behind how quickly cellular immunity forces escape are poorly understood. We developed a novel assay for a common immune escape variant of SIV in macaques. This allowed us to sensitively track the rates of immune escape even when levels of escape mutant or wild-type virus were low. We found that prior immunization of macaques resulted in very rapid immune escape during acute infection. However, when escape starts to occur later, during chronic infection, the rate of immune escape is much more gradual. Thus, both prior vaccination and timing influence the rates of immune escape and provide a fuller picture of the effectiveness of T cell immunity to HIV.

Association between hepatitis C virus (HCV) quasispecies and treatment outcome among patients with chronic hepatitis C has been the subject of many studies. However, these studies focused mainly on viral variable regions (E1 and E2) and usually did not include human immunodeficiency virus (HIV)-positive patients. The aim of the present study was to analyze heterogeneity of the 5' untranslated region (5'UTR) in HCV/HIV coinfected patients treated with interferon and ribavirin. The HCV 5'UTR was amplified from serum and peripheral blood mononuclear cells (PBMC) samples in 37 HCV/HIV coinfected patients treated for chronic hepatitis C. Samples were collected right before treatment, and at 2, 4, 6, 8, 12, 20, 24, 36, 44, 48, 60, and 72 weeks. Heterogeneity of the 5'UTR was analyzed by single strand conformational polymorphism (SSCP), cloning and sequencing. Sustained virological response (SVR) was achieved in 46% of analyzed HCV/HIV co-infected patients. Stable SSCP band pattern was observed in 22 patients (62.9%) and SVR rate among these patients was 23%. Decline in the number of bands and/or shift in band positions were found in 6 patients (17.1%), 5 (83%) of whom achieved SVR (p=0.009). A novel viral genotype was identified in all but one of these patients. In 5 of these 6 patients a new genotype was dominant. 5'UTR heterogeneity may correlate with interferon and ribavirin treatment outcome. In the analyzed group of HCV/HIV coinfected patients, viral quasispecies stability during treatment favored viral persistence, whereas decrease in the number of variants and/or emergence of new variants was associated with SVR. Among injection drug users (IDU) patients, a new genotype may become dominant during treatment, probably due to the presence of mixed infections with various strains, which have different susceptibility to treatment.

RNA viruses have high error rates, and the resulting quasispecies may aid survival of the virus population in the presence of selective pressure. Therefore, it has been theorized that RNA viruses require high error rates for survival, and that a virus with high fidelity would be less able to cope in complex environments. We previously isolated and characterized poliovirus with a mutation in the viral polymerase, 3D-G64S, which confers resistance to mutagenic nucleotide analogs via increased fidelity. The 3D-G64S virus was less pathogenic than wild-type virus in poliovirus-receptor transgenic mice, even though only slight growth defects were observed in tissue culture. To determine whether the high-fidelity phenotype of the 3D-G64S virus could decrease its fitness under a defined selective pressure, we compared growth of the 3D-G64S virus and 3D wild-type virus in the context of a revertible attenuating point mutation, 2C-F28S. Even with a 10-fold input advantage, the 3D-G64S virus was unable to compete with 3D wild-type virus in the context of the revertible attenuating mutation; however, in the context of a non-revertible version of the 2C-F28S attenuating mutation, 3D-G64S virus matched the replication of 3D wild-type virus. Therefore, the 3D-G64S high-fidelity phenotype reduced viral fitness under a defined selective pressure, making it likely that the reduced spread in murine tissue could be caused by the increased fidelity of the viral polymerase., Synopsis RNA viruses have the highest error rates in nature, resulting in the likelihood that each virus differs from other viruses in the population by one or more mutations. The consequence of this “infidelity” is that the viral population as a whole, under selective pressure from the immune system or antiviral drugs, may benefit from adaptive changes in a subset of its members. Therefore, it has been theorized that RNA viruses need high error rates to survive in complex environments. We tested this hypothesis using a drug-resistant poliovirus that contains a mutation in its polymerase that reduces errors during replication. We found that this high-fidelity mutant virus has reduced growth in mice, a complex environment where mutations may be required for growth and spread within the infected animal. At least part of this attenuation is likely due to the high fidelity of this mutant virus, since it was unable to compete with the low-fidelity version of the virus in the context of a defined selective pressure. Therefore, it is likely that mutations do benefit viral populations, especially in complex environments such as an infected animal or human.

Background Epidemiological investigations of infectious disease are mainly dependent on indirect contact information and only occasionally assisted by characterization of pathogen sequence variation from clinical isolates. Direct sequence analysis of the pathogen, particularly at a population level, is generally thought to be too cumbersome, technically difficult, and expensive. We present here a novel application of mass spectrometry (MS)–based technology in characterizing viral sequence variations that overcomes these problems, and we apply it retrospectively to the severe acute respiratory syndrome (SARS) outbreak in Singapore. Methods and Findings The success rate of the MS-based analysis for detecting SARS coronavirus (SARS-CoV) sequence variations was determined to be 95% with 75 copies of viral RNA per reaction, which is sufficient to directly analyze both clinical and cultured samples. Analysis of 13 SARS-CoV isolates from the different stages of the Singapore outbreak identified nine sequence variations that could define the molecular relationship between them and pointed to a new, previously unidentified, primary route of introduction of SARS-CoV into the Singapore population. Our direct determination of viral sequence variation from a clinical sample also clarified an unresolved epidemiological link regarding the acquisition of SARS in a German patient. We were also able to detect heterogeneous viral sequences in primary lung tissues, suggesting a possible coevolution of quasispecies of virus within a single host. Conclusion This study has further demonstrated the importance of improving clinical and epidemiological studies of pathogen transmission through the use of genetic analysis and has revealed the MS-based analysis to be a sensitive and accurate method for characterizing SARS-CoV genetic variations in clinical samples. We suggest that this approach should be used routinely during outbreaks of a wide variety of agents, in order to allow the most effective control., Mass spectrometry-based sequence analysis provides a sensitive and accurate method to characterize genetic variation of the SARS coronavirus in clinical samples

Artificial recombinants can be generated during PCR when more than two genetically distinct templates coexist in a single PCR reaction. These recombinant amplicons can lead to the false interpretation of genetic diversity and incorrect identification of biological phenotypes that do not exist in vivo. We investigated how recombination between 2 or 35 genetically distinct HIV-1 genomes was affected by different PCR conditions using the parallel allele-specific sequencing (PASS) assay and the next generation sequencing method. In a standard PCR condition, about 40% of amplicons in a PCR reaction were recombinants. The high recombination frequency could be significantly reduced if the number of amplicons in a PCR reaction was below a threshold of 10(13)-10(14) using low thermal cycles, fewer input templates, and longer extension time. Heteroduplexes (each DNA strand from a distinct template) were present at a large proportion in the PCR products when more thermal cycles, more templates, and shorter extension time were used. Importantly, the majority of recombinants were identified in heteroduplexes, indicating that the recombinants were mainly generated through heteroduplexes. Since prematurely terminated extension fragments can form heteroduplexes by annealing to different templates during PCR amplification, recombination has a better chance to occur with samples containing different genomes when the number of amplicons accumulate over the threshold. New technologies are warranted to accurately characterize complex quasispecies gene populations.

Trade-offs between throughput, read length, and error rates in high-throughput sequencing limit certain applications such as monitoring viral quasispecies. Here, we describe a molecular-based tag linkage method that allows assemblage of short sequence reads into long DNA fragments. It enables haplotype phasing with high accuracy and sensitivity to interrogate individual viral sequences in a quasispecies. This approach is demonstrated to deduce ∼ 2000 unique 1.3 kb viral sequences from HIV-1 quasispecies in vivo and after passaging ex vivo with a detection limit of ∼ 0.005% to ∼ 0.001%. Reproducibility of the method is validated quantitatively and qualitatively by a technical replicate. This approach can improve monitoring of the genetic architecture and evolution dynamics in any quasispecies population.

Quasispecies variants and recombination were studied longitudinally in an emergent outbreak of beak and feather disease virus (BFDV) infection in the orange-bellied parrot (Neophema chrysogaster). Detailed health monitoring and the small population size (

Detection of recent HIV infections is a prerequisite for reliable estimations of transmitted HIV drug resistance (t-HIVDR) and incidence. However, accurately identifying recent HIV infection is challenging due partially to the limitations of current serological tests. Ambiguous nucleotides are newly emerged mutations in quasispecies, and accumulate by time of viral infection. We utilized ambiguous mutations to establish a measurement for detecting recent HIV infection and monitoring early HIVDR development. Ambiguous nucleotides were extracted from HIV-1 pol-gene sequences in the datasets of recent (HIVDR threshold surveys [HIVDR-TS] in 7 countries; n=416) and established infections (1 HIVDR monitoring survey at baseline; n=271). An ambiguous mutation index of 2.04×10(-3) nts/site was detected in HIV-1 recent infections which is equivalent to the HIV-1 substitution rate (2×10(-3) nts/site/year) reported before. However, significantly higher index (14.41×10(-3) nts/site) was revealed with established infections. Using this substitution rate, 75.2% subjects in HIVDR-TS with the exception of the Vietnam dataset and 3.3% those in HIVDR-baseline were classified as recent infection within one year. We also calculated mutation scores at amino acid level at HIVDR sites based on ambiguous or fitted mutations. The overall mutation scores caused by ambiguous mutations increased (0.54×10(-2)3.48×10(-2)/DR-site) whereas those caused by fitted mutations remained stable (7.50-7.89×10(-2)/DR-site) in both recent and established infections, indicating that t-HIVDR exists in drug-naïve populations regardless of infection status in which new HIVDR continues to emerge. Our findings suggest that characterization of ambiguous mutations in HIV may serve as an additional tool to differentiate recent from established infections and to monitor HIVDR emergence.

Viral variants with decreased susceptibility to HCV protease inhibitors (PIs) occur naturally and preexist at low levels within HCV populations. In patients failing PI monotherapy, single and double mutants conferring intermediate to high-level resistance to PIs have been selected in vivo. The abundance, temporal dynamics and linkage of naturally occurring resistance-associated variants (RAVs), however, have not been characterized in detail. Here, using high-density pyrosequencing, we analyzed HCV NS3 gene segments from 20 subjects with chronic HCV infection, including 12 subjects before and after liver transplantation. Bioinformatics analysis revealed that Q80 substitution was a dominant variant in 40% of the subjects, whereas other RAVs circulate at low levels within quasispecies populations. Low frequency mutation linkage was detectable by Illumina paired-end sequencing in as low as 0.5% of the mock populations constructed from in vitro RNA transcripts but were uncommon in vivo. We show that naturally occurring RAVs are common and can persist long term following liver transplant at low levels not readily detectable by conventional sequencing. Our results indicate that mutation linkage at low levels could be identified using the Illumina paired-end approach. The methods described here should facilitate the analysis of low frequency HCV drug resistance, mutation linkage and evolution, which may inform future therapeutic strategies in patients undergoing direct acting antiviral therapies.

Influenza A(H1N1)pdm09 viruses cause sporadically very severe disease including fatal clinical outcomes associated with pneumonia, viremia and myocarditis. A mutation characterized by the substitution of aspartic acid (wild-type) to glycine at position 222 within the haemagglutinin gene (HA-D222G) was recorded during the 2009 H1N1 pandemic in Germany and other countries with significant frequency in fatal and severe cases. Additionally, A(H1N1)pdm09 viruses exhibiting the polymorphism HA-222D/G/N were detected both in the respiratory tract and in blood. Specimens from mild, fatal and severe cases were collected to study the heterogeneity of HA-222 in A(H1N1)pdm09 viruses circulating in Germany between 2009 and 2011. In order to enable rapid and large scale analysis we designed a pyrosequencing (PSQ) assay. In 2009/2010, the 222D wild-type of A(H1N1)pdm09 viruses predominated in fatal and severe outcomes. Moreover, co-circulating virus mutants exhibiting a D222G or D222E substitution (8/6%) as well as HA-222 quasispecies were identified (10%). Both the 222D/G and the 222D/G/N/V/Y polymorphisms were confirmed by TA cloning. PSQ analyses of viruses associated with mild outcomes revealed mainly the wild-type 222D and no D222G change in both seasons. However, an increase of variants with 222D/G polymorphism (60%) was characteristic for A(H1N1)pdm09 viruses causing fatal and severe cases in the season 2010/2011. Pure 222G viruses were not observed. Our results support the hypothesis that the D222G change may result from adaptation of viral receptor specificity to the lower respiratory tract. This could explain why transmission of the 222G variant is less frequent among humans. Thus, amino acid changes at HA position 222 may be the result of viral intra-host evolution leading to the generation of variants with an altered viral tropism.

The high mutation rate of RNA viruses enables a diverse genetic population of viral genotypes to exist within a single infected host. In-host genetic diversity could better position the virus population to respond and adapt to a diverse array of selective pressures such as host-switching events. Multiple new coronaviruses, including SARS, have been identified in human samples just within the last ten years, demonstrating the potential of coronaviruses as emergent human pathogens. Deep sequencing was used to characterize genomic changes in coronavirus quasispecies during simulated host-switching. Three bovine nasal samples infected with bovine coronavirus were used to infect human and bovine macrophage and lung cell lines. The virus reproduced relatively well in macrophages, but the lung cell lines were not infected efficiently enough to allow passage of non lab-adapted samples. Approximately 12 kb of the genome was amplified before and after passage and sequenced at average coverages of nearly 950×(454 sequencing) and 38,000×(Illumina). The consensus sequence of many of the passaged samples had a 12 nucleotide insert in the consensus sequence of the spike gene, and multiple point mutations were associated with the presence of the insert. Deep sequencing revealed that the insert was present but very rare in the unpassaged samples and could quickly shift to dominate the population when placed in a different environment. The insert coded for three arginine residues, occurred in a region associated with fusion entry into host cells, and may allow infection of new cell types via heparin sulfate binding. Analysis of the deep sequencing data indicated that two distinct genotypes circulated at different frequency levels in each sample, and support the hypothesis that the mutations present in passaged strains were "selected" from a pre-existing pool rather than through de novo mutation and subsequent population fixation.

Background Hepatitis C is a treatment-resistant disease affecting millions of people worldwide. The hepatitis C virus (HCV) genome is a single-stranded RNA molecule. After infection of the host cell, viral RNA is translated into a polyprotein that is cleaved by host and viral proteinases into functional, structural and non-structural, viral proteins. Cleavage of the polyprotein involves the viral NS3/4A proteinase, a proven drug target. HCV mutates as it replicates and, as a result, multiple emerging quasispecies become rapidly resistant to anti-virals, including NS3/4A inhibitors. Methodology/Principal Findings To circumvent drug resistance and complement the existing anti-virals, NS3/4A inhibitors, which are additional and distinct from the FDA-approved telaprevir and boceprevir α-ketoamide inhibitors, are required. To test potential new avenues for inhibitor development, we have probed several distinct exosites of NS3/4A which are either outside of or partially overlapping with the active site groove of the proteinase. For this purpose, we employed virtual ligand screening using the 275,000 compound library of the Developmental Therapeutics Program (NCI/NIH) and the X-ray crystal structure of NS3/4A as a ligand source and a target, respectively. As a result, we identified several novel, previously uncharacterized, nanomolar range inhibitory scaffolds, which suppressed of the NS3/4A activity in vitro and replication of a sub-genomic HCV RNA replicon with a luciferase reporter in human hepatocarcinoma cells. The binding sites of these novel inhibitors do not significantly overlap with those of α-ketoamides. As a result, the most common resistant mutations, including V36M, R155K, A156T, D168A and V170A, did not considerably diminish the inhibitory potency of certain novel inhibitor scaffolds we identified. Conclusions/Significance Overall, the further optimization of both the in silico strategy and software platform we developed and lead compounds we identified may lead to advances in novel anti-virals.

To assess the level of intra-patient diversity and evolution of HIV-1C non-structural genes in primary infection, viral quasispecies obtained by single genome amplification (SGA) at multiple sampling timepoints up to 500 days post-seroconversion (p/s) were analyzed. The mean intra-patient diversity was 0.11% (95% CI; 0.02 to 0.20) for vif, 0.23% (95% CI; 0.08 to 0.38) for vpr, 0.35% (95% CI; -0.05 to 0.75) for vpu, 0.18% (95% CI; 0.01 to 0.35) for tat exon 1 and 0.30% (95% CI; 0.02 to 0.58) for rev exon 1 during the time period 0 to 90 days p/s. The intra-patient diversity increased gradually in all non-structural genes over the first year of HIV-1 infection, which was evident from the vif mean intra-patient diversity of 0.46% (95% CI; 0.28 to 0.64), vpr 0.44% (95% CI; 0.24 to 0.64), vpu 0.84% (95% CI; 0.55 to 1.13), tat exon 1 0.35% (95% CI; 0.14 to 0.56 ) and rev exon 1 0.42% (95% CI; 0.18 to 0.66) during the time period of 181 to 500 days p/s. There was a statistically significant increase in viral diversity for vif (p = 0.013) and vpu (p = 0.002). No associations between levels of viral diversity within the non-structural genes and HIV-1 RNA load during primary infection were found. The study details the dynamics of the non-structural viral genes during the early stages of HIV-1C infection.

Monoclonal antibodies directed against hepatitis C virus (HCV) E2 protein can neutralize cell-cultured HCV and pseudoparticles expressing envelopes derived from multiple HCV subtypes. For example, based on antibody blocking experiments and alanine scanning mutagenesis, it was proposed that the AR3B monoclonal antibody recognized a discontinuous conformational epitope comprised of amino acid residues 396-424, 436-447, and 523-540 of HCV E2 envelope protein. Intriguingly, one of these segments (436-447) overlapped with hypervariable region 3 (HVR3), a domain that exhibited significant intrahost and interhost genetic diversity. To reconcile these observations, amino-acid sequence variability was examined and homology-based structural modelling of E2 based on tick-borne encephalitis virus (TBEV) E protein was performed based on 413 HCV sequences derived from 18 subjects with chronic hepatitis C. Here we report that despite a high degree of amino-acid sequence variability, the three-dimensional structure of E2 is remarkably conserved, suggesting broad recognition of structural determinants rather than specific residues. Regions 396-424 and 523-540 were largely exposed and in close spatial proximity at the surface of E2. In contrast, region 436-447, which overlaps with HVR3, was >35 Å away, and estimates of buried surface were inconsistent with HVR3 being part of the AR3B binding interface. High-throughput structural analysis of HCV quasispecies could facilitate the development of novel vaccines that target conserved structural features of HCV envelope and elicit neutralizing antibody responses that are less vulnerable to viral escape.

Background Pyrosequencing technology has the potential to rapidly sequence HIV-1 viral quasispecies without requiring the traditional approach of cloning. In this study, we investigated the utility of ultra-deep pyrosequencing to characterize genetic diversity of the HIV-1 gag quasispecies and assessed the possible contribution of pyrosequencing technology in studying HIV-1 biology and evolution. Methodology/Principal Findings HIV-1 gag gene was amplified from 96 patients using nested PCR. The PCR products were cloned and sequenced using capillary based Sanger fluorescent dideoxy termination sequencing. The same PCR products were also directly sequenced using the 454 pyrosequencing technology. The two sequencing methods were evaluated for their ability to characterize quasispecies variation, and to reveal sites under host immune pressure for their putative functional significance. A total of 14,034 variations were identified by 454 pyrosequencing versus 3,632 variations by Sanger clone-based (SCB) sequencing. 11,050 of these variations were detected only by pyrosequencing. These undetected variations were located in the HIV-1 Gag region which is known to contain putative cytotoxic T lymphocyte (CTL) and neutralizing antibody epitopes, and sites related to virus assembly and packaging. Analysis of the positively selected sites derived by the two sequencing methods identified several differences. All of them were located within the CTL epitope regions. Conclusions/Significance Ultra-deep pyrosequencing has proven to be a powerful tool for characterization of HIV-1 genetic diversity with enhanced sensitivity, efficiency, and accuracy. It also improved reliability of downstream evolutionary and functional analysis of HIV-1 quasispecies.

Background and Aims The hepatitis C virus (HCV) invariably shows wide heterogeneity in infected patients, referred to as a quasispecies population. Massive amounts of genetic information due to the abundance of HCV variants could be an obstacle to evaluate the viral genetic heterogeneity in detail. Methods Using a newly developed massive-parallel ultra-deep sequencing technique, we investigated the viral genetic heterogeneity in 27 chronic hepatitis C patients receiving peg-interferon (IFN) α2b plus ribavirin therapy. Results Ultra-deep sequencing determined a total of more than 10 million nucleotides of the HCV genome, corresponding to a mean of more than 1000 clones in each specimen, and unveiled extremely high genetic heterogeneity in the genotype 1b HCV population. There was no significant difference in the level of viral complexity between immediate virologic responders and non-responders at baseline (p = 0.39). Immediate virologic responders (n = 8) showed a significant reduction in the genetic complexity spanning all the viral genetic regions at the early phase of IFN administration (p = 0.037). In contrast, non-virologic responders (n = 8) showed no significant changes in the level of viral quasispecies (p = 0.12), indicating that very few viral clones are sensitive to IFN treatment. We also demonstrated that clones resistant to direct-acting antivirals for HCV, such as viral protease and polymerase inhibitors, preexist with various abundances in all 27 treatment-naive patients, suggesting the risk of the development of drug resistance against these agents. Conclusion Use of the ultra-deep sequencing technology revealed massive genetic heterogeneity of HCV, which has important implications regarding the treatment response and outcome of antiviral therapy.

Genetic analysis of hepatitis B virus (HBV) frequently involves study of intra-host variants, identification of which is commonly achieved using short regions of the HBV genome. However, the use of short sequences significantly limits evaluation of genetic relatedness among HBV strains. Although analysis of HBV complete genomes using genetic cloning has been developed, its application is highly labor intensive and practiced only infrequently. We describe here a novel approach to whole genome (WG) HBV quasispecies analysis based on end-point, limiting-dilution real-time PCR (EPLD-PCR) for amplification of single HBV genome variants, and their subsequent sequencing. EPLD-PCR was used to analyze WG quasispecies from serum samples of patients (n = 38) infected with HBV genotypes A, B, C, D, E and G. Phylogenetic analysis of the EPLD-isolated HBV-WG quasispecies showed the presence of mixed genotypes, recombinant variants and sub-populations of the virus. A critical observation was that HBV-WG consensus sequences obtained by direct sequencing of PCR fragments without EPLD are genetically close, but not always identical to the major HBV variants in the intra-host population, thus indicating that consensus sequences should be judiciously used in genetic analysis. Sequence-based studies of HBV WG quasispecies should afford a more accurate assessment of HBV evolution in various clinical and epidemiological settings.

Mucosal transmission of HIV is inefficient. The virus must breach physical barriers before it infects mucosal CD4+ T cells. Low-level viral replication occurs initially in mucosal CD4+ T cells, but within days high-level replication occurs in Peyer's patches, the gut lamina propria and mesenteric lymph nodes. Understanding the early events in HIV transmission may provide valuable information relevant to the development of an HIV vaccine. The viral quasispecies in a donor contracts through a genetic bottleneck in the recipient, such that, in low-risk settings, infection is frequently established by a single founder virus. Early-transmitting viruses in subtypes A and C mucosal transmission tend to encode gp120s with reduced numbers of N-linked glycosylation sites at specific positions throughout the V1-V4 domains, relative to typical chronically replicating isolates in the donor quasispecies. The transmission advantage gained by the absence of these N-linked glycosylation sites is unknown. Using primary α4β7 +/CD4+ T cells and a flow-cytometry based steady-state binding assay we show that the removal of transmission-associated N-linked glycosylation sites results in large increases in the specific reactivity of gp120 for integrin- α4β7. High-affinity for integrin α4β7, although not found in many gp120s, was observed in early-transmitting gp120s that we analyzed. Increased α4β7 affinity is mediated by sequences encoded in gp120 V1/V2. α4β7-reactivity was also influenced by N-linked glycosylation sites located in C3/V4. These results suggest that the genetic bottleneck that occurs after transmission may frequently involve a relative requirement for the productive infection of α4β7 +/CD4+ T cells. Early-transmitting gp120s were further distinguished by their dependence on avidity-effects to interact with CD4, suggesting that these gp120s bear unusual structural features not present in many well-characterized gp120s derived from chronically replicating viruses. Understanding the structural features that characterize early-transmitting gp120s may aid in the design of an effective gp120-based subunit vaccine., Author Summary In the first days following sexual transmission, HIV replication occurs initially at relatively low levels in mucosal tissues because of a paucity of CD4+ T cell targets for the virus to infect. After a period of days, virus accesses specific gut tissues that are enriched in activated CD4+ T cells, where near-exponential replication ensues. The period of time before HIV accesses gut tissues represents a window of opportunity where a microbicide, pre-exposure and/or post-exposure antiretroviral prophylaxis or a vaccine-induced immune response could block infection. We previously reported that the HIV envelope protein gp120 binds to integrin α4β7 on the surface of CD4+ T cells. α4β7 mediates the homing of CD4+ T cells into the gut tissues where HIV can replicate exponentially. Here we report that the genotypic features that distinguish viruses isolated within the first month after infection, termed early-transmitting isolates, promote increased steady-state reactivity with α4β7. This property likely provides these viruses with enhanced transmission-fitness. These results suggest that the infection of α4β7 +/CD4+ T cells can play an important role early in HIV transmission. These findings have potentially important implications in the design of interventions to block the mucosal transmission of HIV.

HIV/HCV coinfected individuals under highly active antiretroviral therapy (HAART) represent an interesting model for the investigation of the role played by the immune system in driving the evolution of the HCV quasispecies. We prospectively studied the intra-host evolution of the HCV heterogeneity in 8 coinfected subjects, selected from a cohort of 32 patients initiating HAART: 5 immunological responders (group A) and 3 immunological non-responders (group B), and in two HCV singly infected controls not assuming drugs (group C). For all these subjects at least two serial samples obtained at the first observation (before HAART) and more than 1 year later, underwent clonal sequence analysis of partial E1/E2 sequences, encompassing the whole HVR1. Evolutionary rates, dated phylogenies and population dynamics were co-estimated by using a Bayesian Markov Chain Monte Carlo approach, and site specific selection pressures were estimated by maximum likelihood-based methods. The intra-host evolutionary rates of HCV quasispecies was 10 times higher in subjects treated with HAART than in controls without immunodeficiency (1.9 and 2.3 × 10(-3) sub/site/month in group A and B and 0.29 × 10(-3) sub/site/month in group C individuals). The within-host Bayesian Skyline plot analysis showed an exponential growth of the quasispecies populations in immunological responders, coinciding with a peak in CD4 cell counts. On the contrary, quasispecies population remained constant in group B and in group C controls. A significant positive selection pressure was detected in a half of the patients under HAART and in none of the group C controls. Several sites under significant positive selection were described, mainly included in the HVR1. Our data indicate that different forces, in addition to the selection pressure, drive an exceptionally fast evolution of HCV during HAART immune restoration. We hypothesize that an important role is played by the enlargement of the viral replicative space.

The goal of this paper is to develop a mathematical model that analyzes the selective advantage of the SOS response in unicellular organisms. To this end, this paper develops a quasispecies model that incorporates the SOS response. We consider a unicellular, asexually replicating population of organisms, whose genomes consist of a single, double-stranded DNA molecule, i.e. one chromosome. We assume that repair of post-replication mismatched base-pairs occurs with probability , and that the SOS response is triggered when the total number of mismatched base-pairs is at least . We further assume that the per-mismatch SOS elimination rate is characterized by a first-order rate constant . For a single fitness peak landscape where the master genome can sustain up to mismatches and remain viable, this model is analytically solvable in the limit of infinite sequence length. The results, which are confirmed by stochastic simulations, indicate that the SOS response does indeed confer a fitness advantage to a population, provided that it is only activated when DNA damage is so extensive that a cell will die if it does not attempt to repair its DNA.

Most RNA viruses lack the mechanisms to recognize and correct mutations that arise during genome replication, resulting in quasispecies diversity that is required for pathogenesis and adaptation. However, it is not known how viruses encoding large viral RNA genomes such as the Coronaviridae (26 to 32 kb) balance the requirements for genome stability and quasispecies diversity. Further, the limits of replication infidelity during replication of large RNA genomes and how decreased fidelity impacts virus fitness over time are not known. Our previous work demonstrated that genetic inactivation of the coronavirus exoribonuclease (ExoN) in nonstructural protein 14 (nsp14) of murine hepatitis virus results in a 15-fold decrease in replication fidelity. However, it is not known whether nsp14-ExoN is required for replication fidelity of all coronaviruses, nor the impact of decreased fidelity on genome diversity and fitness during replication and passage. We report here the engineering and recovery of nsp14-ExoN mutant viruses of severe acute respiratory syndrome coronavirus (SARS-CoV) that have stable growth defects and demonstrate a 21-fold increase in mutation frequency during replication in culture. Analysis of complete genome sequences from SARS-ExoN mutant viral clones revealed unique mutation sets in every genome examined from the same round of replication and a total of 100 unique mutations across the genome. Using novel bioinformatic tools and deep sequencing across the full-length genome following 10 population passages in vitro, we demonstrate retention of ExoN mutations and continued increased diversity and mutational load compared to wild-type SARS-CoV. The results define a novel genetic and bioinformatics model for introduction and identification of multi-allelic mutations in replication competent viruses that will be powerful tools for testing the effects of decreased fidelity and increased quasispecies diversity on viral replication, pathogenesis, and evolution., Author Summary Quasispecies diversity is critical to virus fitness, adaptation, and pathogenesis. However, the relationship of fidelity to population diversity is less studied because viral systems with engineered differences in fidelity and bioinformatic methods that robustly measure and compare fidelity and diversity during replication and passage have not been available. Coronaviruses contain the largest and most complex RNA genomes, and encode multiple novel replicase nonstructural proteins (nsps). We previously demonstrated that murine hepatitis virus nsp14-exonuclease (ExoN) activity is required for replication fidelity. In the present report we have generated nsp14-ExoN inactivation mutants of SARS-coronavirus (S-ExoN) that have stable growth defects and dramatically decreased replication fidelity during replication in culture. We used the S-ExoN mutant viruses to define the diversity and stability of the genome during replication and passage, and to test the capacity of deep sequencing to track virus population diversity over time. The experiments demonstrate that viable S-ExoN mutants accumulate large numbers of predominantly unique mutations across the genome, and that increased diversity is continuous over passage. The results establish methods for direct comparison of consensus genome sequences with total population diversity and the impact on viral growth and adaptation.