Next-generation sequencing technologies can be used to analyse genetically heterogeneous samples at unprecedented detail. The high coverage achievable with these methods enables the detection of many low-frequency variants. However, sequencing errors complicate the analysis of mixed populations and result in inflated estimates of genetic diversity. We developed a probabilistic Bayesian approach to minimize the effect of errors on the detection of minority variants. We applied it to pyrosequencing data obtained from a 1.5-kb-fragment of the HIV-1 gag/pol gene in two control and two clinical samples. The effect of PCR amplification was analysed. Error correction resulted in a two- and five-fold decrease of the pyrosequencing base substitution rate, from 0.05% to 0.03% and from 0.25% to 0.05% in the non-PCR and PCR-amplified samples, respectively. We were able to detect viral clones as rare as 0.1% with perfect sequence reconstruction. Probabilistic haplotype inference outperforms the counting-based calling method in both precision and recall. Genetic diversity observed within and between two clinical samples resulted in various patterns of phenotypic drug resistance and suggests a close epidemiological link. We conclude that pyrosequencing can be used to investigate genetically diverse samples with high accuracy if technical errors are properly treated.

Therapeutic success of anti-HIV therapies is limited by the development of drug resistant viruses. These genetic variants display complex mutational patterns in their pol gene, which codes for protease and reverse transcriptase, the molecular targets of current antiretroviral therapy. Genotypic resistance testing depends on the ability to interpret such sequence data, whereas phenotypic resistance testing directly measures relative in vitro susceptibility to a drug. From a set of 650 matched genotype-phenotype pairs we construct regression models for the prediction of phenotypic drug resistance from genotypes. Since the range of resistance factors varies considerably between different drugs, two scoring functions are derived from different sets of predicted phenotypes. Firstly, we compare predicted values to those of samples derived from 178 treatment-naive patients and report the relative deviance. Secondly, estimation of the probability density of 2000 predicted phenotypes gives rise to an intrinsic definition of a susceptible and a resistant subpopulation. Thus, for a predicted phenotype, we calculate the probability of membership in the resistant subpopulation. Both scores provide standardized measures of resistance that can be calculated from the genotype and are comparable between drugs. The geno2pheno system makes these genotype interpretations available via the Internet (http://www.genafor.org/).

Next-generation sequencing (NGS) technologies enable new insights into the diversity of virus populations within their hosts. Diversity estimation is currently restricted to single-nucleotide variants or to local fragments of no more than a few hundred nucleotides defined by the length of sequence reads. To study complex heterogeneous virus populations comprehensively, novel methods are required that allow for complete reconstruction of the individual viral haplotypes. Here, we show that assembly of whole viral genomes of ∼8600 nucleotides length is feasible from mixtures of heterogeneous HIV-1 strains derived from defined combinations of cloned virus strains and from clinical samples of an HIV-1 superinfected individual. Haplotype reconstruction was achieved using optimized experimental protocols and computational methods for amplification, sequencing and assembly. We comparatively assessed the performance of the three NGS platforms 454 Life Sciences/Roche, Illumina and Pacific Biosciences for this task. Our results prove and delineate the feasibility of NGS-based full-length viral haplotype reconstruction and provide new tools for studying evolution and pathogenesis of viruses.

We introduce a mixture model of trees to describe evolutionary processes that are characterized by the ordered accumulation of permanent genetic changes. The basic building block of the model is a directed weighted tree that generates a probability distribution on the set of all patterns of genetic events. We present an EM-like algorithm for learning a mixture model of K trees and show how to determine K with a maximum likelihood approach. As a case study, we consider the accumulation of mutations in the HIV-1 reverse transcriptase that are associated with drug resistance. The fitted model is statistically validated as a density estimator, and the stability of the model topology is analyzed. We obtain a generative probabilistic model for the development of drug resistance in HIV that agrees with biological knowledge. Further applications and extensions of the model are discussed.

Herpes simplex virus (HSV) glycoprotein B (gB) is an integral part of the multicomponent fusion system required for virus entry and cell-cell fusion. Here we investigated the mechanism of viral neutralization by the monoclonal antibody (MAb) 2c, which specifically recognizes the gB of HSV type 1 (HSV-1) and HSV-2. Binding of MAb 2c to a type-common discontinuous epitope of gB resulted in highly efficient neutralization of HSV at the postbinding/prefusion stage and completely abrogated the viral cell-to-cell spread in vitro. Mapping of the antigenic site recognized by MAb 2c to the recently solved crystal structure of the HSV-1 gB ectodomain revealed that its discontinuous epitope is only partially accessible within the observed multidomain trimer conformation of gB, likely representing its postfusion conformation. To investigate how MAb 2c may interact with gB during membrane fusion, we characterized the properties of monovalent (Fab and scFv) and bivalent [IgG and F(ab) 2 ] derivatives of MAb 2c. Our data show that the neutralization capacity of MAb 2c is dependent on cross-linkage of gB trimers. As a result, only bivalent derivatives of MAb 2c exhibited high neutralizing activity in vitro. Notably, bivalent MAb 2c not only was capable of preventing mucocutaneous disease in severely immunodeficient NOD/SCID mice upon vaginal HSV-1 challenge but also protected animals even with neuronal HSV infection. We also report for the first time that an anti-gB specific monoclonal antibody prevents HSV-1-induced encephalitis entirely independently from complement activation, antibodydependent cellular cytotoxicity, and cellular immunity. This indicates the potential for further development of MAb 2c as an anti-HSV drug.Herpes simplex virus (HSV) is a neuroinvasive human pathogen that critically depends on efficient infection of distinct target cells within a host. At the time of primary lytic infection, HSV replicates in peripheral mucocutaneous tissues and releases virions. A decisive characteristic of HSV infections in animals and humans is the establishment of a lifelong latency. HSV spreads from infected epithelial cells to axons of sensory neurons innervating the site of the primary infection, followed by retrograde transport to the respective dorsal root ganglia (12). Recurrent infections result from reactivation in neuronal cells, followed by virus replication and anterograde transport to cells at peripheral sites innervated by the respective neurons. Transmission between cells without diffusion through the extracellular environment represents a major route for HSV to spread between tissues and is thus a very efficient way for circumventing immunological barriers of the humoral immune response. Regardless of the dissemination pathway, however, fusion of the viral envelope with host membranes for delivery of the viral genome across the cellular lipid bilayer is essential for viral replication. In contrast to that of most other enveloped viruses, entry of herpesviruses into mammalian cells requires a multicomponent sy...

ObjectivesIt is still debated if pre-existing minority drug-resistant HIV-1 variants (MVs) affect the virological outcomes of first-line NNRTI-containing ART.MethodsThis Europe-wide case–control study included ART-naive subjects infected with drug-susceptible HIV-1 as revealed by population sequencing, who achieved virological suppression on first-line ART including one NNRTI. Cases experienced virological failure and controls were subjects from the same cohort whose viraemia remained suppressed at a matched time since initiation of ART. Blinded, centralized 454 pyrosequencing with parallel bioinformatic analysis in two laboratories was used to identify MVs in the 1%–25% frequency range. ORs of virological failure according to MV detection were estimated by logistic regression.ResultsTwo hundred and sixty samples (76 cases and 184 controls), mostly subtype B (73.5%), were used for the analysis. Identical MVs were detected in the two laboratories. 31.6% of cases and 16.8% of controls harboured pre-existing MVs. Detection of at least one MV versus no MVs was associated with an increased risk of virological failure (OR = 2.75, 95% CI = 1.35–5.60, P = 0.005); similar associations were observed for at least one MV versus no NRTI MVs (OR = 2.27, 95% CI = 0.76–6.77, P = 0.140) and at least one MV versus no NNRTI MVs (OR = 2.41, 95% CI = 1.12–5.18, P = 0.024). A dose–effect relationship between virological failure and mutational load was found.ConclusionsPre-existing MVs more than double the risk of virological failure to first-line NNRTI-based ART.

The preeminent mode of action of the broad-spectrum antiviral nucleoside ribavirin in the therapy of chronic hepatitis C is currently unresolved. Particularly under contest are possible mutagenic effects of ribavirin that may lead to viral extinction by lethal mutagenesis of the hepatitis C virus (HCV) genome. We applied ultradeep sequencing to determine ribavirin-induced sequence changes in the HCV coding region (nucleotides [nt] 330 to 9351) of patients treated with 6-week ribavirin monotherapy (n ‫؍‬ 6) in comparison to placebo (n ‫؍‬ 6). Baseline HCV RNA levels maximally declined on average by ؊0.8 or ؊0.1 log 10 IU/ml in ribavirin-versus placebo-treated patients. No general increase in rates of nucleotide substitutions in ribavirin-treated patients was observed. However, more HCV genome positions with high G-to-A and C-to-U transition rates were detected between baseline and treatment week 6 in ribavirin-treated patients in comparison to placebo-treated patients (rate of 0.0041 transitions per base pair versus rate of 0.0022 transitions per base pair; P ‫؍‬ 0.049). Similarly, the sensitive detection of low-frequency minority variants by statistical filtering indicated significantly more positions with G-to-A and C-to-U transitions in ribavirin-treated patients than in placebo-treated patients (rate of 0.0331 transitions versus rate of 0.0186 transitions per G/C-containing position at baseline; P ‫؍‬ 0.018). In contrast, non-ribavirin-associated A-to-G and U-to-C transitions were not enriched in the ribavirin group (P ‫؍‬ 0.152). We conclude that ribavirin exerts a mutagenic effect on the virus in patients with chronic hepatitis C by facilitating G-to-A and C-to-U nucleotide transitions.

No abstract