Genotypic tropism testing by massively parallel sequencing: qualitative and quantitative analysis

Development and Validation of a Genotypic Assay to Quantify CXCR4- and CCR5-Tropic Human Immunodeficiency Virus Type-1 (HIV-1) Populations and a Comparison to Trofile®

HIV-1 typically infects cells via the CD4 receptor and CCR5 or CXCR4 co-receptors. Maraviroc is a CCR5-specific viral entry inhibitor; knowledge of viral co-receptor specificity is important prior to usage. We developed and validated an economical V3-env Illumina-based assay to detect and quantify the frequency of viruses utilizing each co-receptor. Plasma from 54 HIV+ participants (subtype B) was tested. The viral template cDNA was generated from plasma RNA with unique molecular identifiers (UMIs). The sequences were aligned and collapsed by the UMIs with a custom bioinformatics pipeline. Co-receptor usage, determined by codon analysis and online phenotype predictors PSSM and Geno2pheno, were compared to existing Trofile® data. The cost of V3-UMI was tallied. The sequences interpreted by Geno2pheno using the most conservative cut-off, a 2% false-positive-rate (FPR), predicted CXCR4 usage with the greatest sensitivity (76%) and specificity (100%); PSSM and codon analysis had similar sensitivity and lower specificity. Discordant Trofile® and genotypic results were more common when participants had specimens from different dates analyzed by either assay. V3-UMI reagents cost USD$62/specimen. A batch of ≤20 specimens required 5 h of technical time across 1.5 days. V3-UMI predicts HIV tropism at a sensitivity and specificity similar to those of Trofile®, is relatively inexpensive, and could be performed by most central laboratories. The adoption of V3-UMI could expand HIV drug therapeutic options in lower-resource settings that currently do not have access to phenotypic HIV tropism testing.

In-depth virological and immunological characterization of HIV-1 cure after CCR5Δ32/Δ32 allogeneic hematopoietic stem cell transplantation

Despite scientific evidence originating from two patients published to date that CCR5Δ32/Δ32 hematopoietic stem cell transplantation (HSCT) can cure human immunodeficiency virus type 1 (HIV-1), the knowledge of immunological and virological correlates of cure is limited. Here we characterize a case of long-term HIV-1 remission of a 53-year-old male who was carefully monitored for more than 9 years after allogeneic CCR5Δ32/Δ32 HSCT performed for acute myeloid leukemia. Despite sporadic traces of HIV-1 DNA detected by droplet digital PCR and in situ hybridization assays in peripheral T cell subsets and tissue-derived samples, repeated ex vivo quantitative and in vivo outgrowth assays in humanized mice did not reveal replication-competent virus. Low levels of immune activation and waning HIV-1-specific humoral and cellular immune responses indicated a lack of ongoing antigen production. Four years after analytical treatment interruption, the absence of a viral rebound and the lack of immunological correlates of HIV-1 antigen persistence are strong evidence for HIV-1 cure after CCR5Δ32/Δ32 HSCT.

A Review of Next Generation Sequencing Methods and its Applications in Laboratory Diagnosis

Next-generation sequencing (NGS) is a new technology used to detect the sequence of DNA and RNA and to detect mutations or variations of significance. NGS generates large quantities of sequence data within a short time duration. The various types of sequencing includes Sanger Sequencing, Pyrosequencing, Sequencing by Synthesis (Illumina), Ligation (SoLID), Single molecule Fluorescent Sequencing (Helicos), Single molecule Real time Sequencing (Pacbio), Semiconductor sequencing (Ion torrent technology), Nanopore sequencing and fourth generation sequencing. These methods of sequencing have been modified and improved over the years such that it has become cost effective and accessible to diagnostic laboratories. Management of Outbreaks, rapid identification of bacteria, molecular case finding, taxonomy, detection of the zoonotic agents and guiding prevention strategies in HIV outbreaks are just a few of the many applications of Next Generation sequencing in clinical microbiology.

Monophylogenetic HIV-1C epidemic in Ethiopia is dominated by CCR5-tropic viruses-an analysis of a prospective country-wide cohort

Background

CCR5 coreceptor using HIV-1 subtype C (HIV-1C) has been reported to dominate the Ethiopian epidemic. However, almost all data have been obtained from two large cities in the central and north-west regions and recent data is lacking.

Methods

Plasma were obtained from 420 treatment-naïve patients recruited 2009–2011 to a large country-wide Ethiopian cohort. The V3 region was sequenced and the co-receptor tropism was predicted by the clinical and clonal models of the geno2pheno tool at different false positive rates (fpr) and for subtype. In an intention to treat analysis the impact of baseline tropism on outcome of antiretroviral therapy was evaluated.

Results

V3 loop sequencing was successful in 352 (84%) patients. HIV-1C was found in 350 (99.4%) and HIV-1A in two (0.6%) patients. When comparing the geno2pheno fpr10% clonal and clinical models, 24.4% predictions were discordant. X4-virus was predicted in 17.0 and 19.0%, respectively, but the predictions were concordant in only 6%. At fpr5%, concordant X4-virus predictions were obtained in 3.1%. The proportion of X4-tropic virus (clonal fpr10%) increased from 5.6 to 17.3% (p < 0.001) when 387 Ethiopian V3 loop sequences dated from 1984 to 2003 were compared with ours. In an intention to treat analysis, 67.9% reached treatment success at month 6 and only 50% at month 12. Only age and not tropism predicted therapy outcome and no difference was found in CD4+ cell gain between R5-tropic and X4-tropic infected patients. At viral failure, R5 to X4 switch was rare while X4 to R5 switch occurred more frequently (month 6: p = 0.006; month 12: p = 0.078).

Conclusion

The HIV-1C epidemic is monophylogenetic in all regions of Ethiopia and R5-tropic virus dominates, even in patients with advanced immunodeficiency, although the proportion of X4-tropic virus seems to have increased over the last two decades. Geno2pheno clinical and clonal prediction models show a large discrepancy at fpr10%, but not at fpr5%. Hence further studies are needed to assess the utility of genotypic tropism testing in HIV-1C. In ITT analysis only age and not tropism influenced the outcome.

Statistical correlation of nonconservative substitutions of HIV gp41 variable amino acid residues with the R5X4 HIV-1 phenotype

BACKGROUND

The interaction of the envelope glycoprotein of HIV-1 (gp120/gp41) with coreceptor molecules has important implications for specific cellular targeting and pathogenesis. Experimental and theoretical evidences have shown a role for gp41 in coreceptor tropism, although there is no consensus about the positions involved. Here we analyze the association of physicochemical properties of gp41 amino acid residues with viral tropism (X4, R5, and R5X4) using a large set of HIV-1 sequences. Under the assumption that conserved regions define the complex structural features essential for protein function, we focused our search only on amino acids in the gp41 variable regions.

METHODS

Gp41 amino acid sequences of 2823 HIV-1 strains from all clades with known coreceptor tropism were retrieved from Los Alamos HIV Database. Consensus sequences were constructed for homologous sequences (those obtained from the same patient and having the same tropism) in order to avoid bias due to sequence overrepresentation, and the variability (entropy) per site was determined. Comparisons of hydropathy index (HI) and charge (Q) of amino acid residues at highly variable positions between coreceptor groups were performed using two non-parametrical tests and Benjamini-Hochberg correction. Pearson's correlation analysis was performed to determine covariance of HI and Q values.

RESULTS

Calculation of variability per site rendered 58 highly variable amino acid positions. Of these, statistical analysis rendered significantly different HI or Q only for the R5 vs. R5X4 comparison at twelve positions: 535, 602, 619, 636, 640, 641, 658, 662, 667, 723, 756 and 841. The largest differences in particular amino acid frequencies between coreceptor groups were found at 619, 636, 640, 641, 662, 723 and 756. A hydrophobic tendency of residues 619, 640, 641, 723 and 756, along with a hydrophilic/charged tendency at residues 636 and 662 was observed in R5X4 with respect to R5 sequences. HI of position 640 covariated with that of 602, 619, 636, 662, and 756.

CONCLUSIONS

Variability and significant correlations of physicochemical properties with viral phenotype suggest that substitutions at residues in the loop (602 and 619), the HR2 (636, 640, 641, 662), and the C-terminal tail (723, 756) of gp41 may contribute to phenotype of R5X4 strains.

Comparison of 454 Ultra-Deep Sequencing and Allele-Specific Real-Time PCR with Regard to the Detection of Emerging Drug-Resistant Minor HIV-1 Variants after Antiretroviral Prophylaxis for Vertical Transmission

Background

Pregnant HIV-infected women were screened for the development of HIV-1 drug resistance after implementation of a triple-antiretroviral transmission prophylaxis as recommended by the WHO in 2006. The study offered the opportunity to compare amplicon-based 454 ultra-deep sequencing (UDS) and allele-specific real-time PCR (ASPCR) for the detection of drug-resistant minor variants in the HIV-1 reverse transcriptase (RT).

Methods

Plasma samples from 34 Tanzanian women were previously analysed by ASPCR for key resistance mutations in the viral RT selected by AZT, 3TC, and NVP (K70R, K103N, Y181C, M184V, T215Y/F). In this study, the RT region of the same samples was investigated by amplicon-based UDS for resistance mutations using the 454 GS FLX System.

Results

Drug-resistant HIV-variants were identified in 69% (20/29) of women by UDS and in 45% (13/29) by ASPCR. The absolute number of resistance mutations identified by UDS was twice that identified by ASPCR (45 vs 24). By UDS 14 of 24 ASPCR-detected resistance mutations were identified at the same position. The overall concordance between UDS and ASPCR was 61.0% (25/41). The proportions of variants quantified by UDS were approximately 2–3 times lower than by ASPCR. Amplicon generation from samples with viral loads below 20,000 copies/ml failed more frequently by UDS compared to ASPCR (limit of detection = 650 copies/ml), resulting in missing or insufficient sequence coverage.

Conclusions

Both methods can provide useful information about drug-resistant minor HIV-1 variants. ASPCR has a higher sensitivity than UDS, but is restricted to single resistance mutations.

In contrast, UDS is limited by its requirement for high viral loads to achieve sufficient sequence coverage, but the sequence information reveals the complete resistance patterns within the genomic region analysed. Improvements to the UDS limit of detection are in progress, and UDS could then facilitate monitoring of drug-resistant minor variants in the HIV-1 quasispecies.

HIV-1 Coreceptor Usage Assessment by Ultra-Deep Pyrosequencing and Response to Maraviroc

Background

Maraviroc is an HIV entry inhibitor that alters the conformation of CCR5 and is poorly efficient in patients infected by viruses that use CXCR4 as an entry coreceptor. The goal of this study was to assess the capacity of ultra-deep pyrosequencing (UDPS) and different data analysis approaches to characterize HIV tropism at baseline and predict the therapeutic outcome on maraviroc treatment.

Methods

113 patients with detectable HIV-1 RNA on HAART were treated with maraviroc. The virological response was assessed at months 1, 3 and 6. The sequence of the HIV V3 loop was determined at baseline and prediction of maraviroc response by different software and interpretation algorithms was analyzed.

Results

UDPS followed by analysis with the Pyrotrop software or geno2pheno algorithm provided better prediction of the response to maraviroc than Sanger sequencing. We also found that the H34Y/S substitution in the V3 loop was the strongest individual predictor of maraviroc response, stronger than substitutions at positions 11 or 25 classically used in interpretation algorithms.

Conclusions

UDPS is a powerful tool that can be used with confidence to predict maraviroc response in HIV-1-infected patients. Improvement of the predictive value of interpretation algorithms is possible and our results suggest that adding the H34S/Y substitution would substantially improve the performance of the 11/25/charge rule.

HIV-1 tropism dynamics and phylogenetic analysis from longitudinal ultra-deep sequencing data of CCR5- and CXCR4-using variants

OBJECTIVE

Coreceptor switch from CCR5 to CXCR4 is associated with HIV disease progression. The molecular and evolutionary mechanisms underlying the CCR5 to CXCR4 switch are the focus of intense recent research. We studied the HIV-1 tropism dynamics in relation to coreceptor usage, the nature of quasispecies from ultra deep sequencing (UDPS) data and their phylogenetic relationships.

METHODS

Here, we characterized C2-V3-C3 sequences of HIV obtained from 19 patients followed up for 54 to 114 months using UDPS, with further genotyping and phylogenetic analysis for coreceptor usage. HIV quasispecies diversity and variability as well as HIV plasma viral load were measured longitudinally and their relationship with the HIV coreceptor usage was analyzed. The longitudinal UDPS data were submitted to phylogenetic analysis and sampling times and coreceptor usage were mapped onto the trees obtained.

RESULTS

Although a temporal viral genetic structuring was evident, the persistence of several viral lineages evolving independently along the infection was statistically supported, indicating a complex scenario for the evolution of viral quasispecies. HIV X4-using variants were present in most of our patients, exhibiting a dissimilar inter- and intra-patient predominance as the component of quasispecies even on antiretroviral therapy. The viral populations from some of the patients studied displayed evidences of the evolution of X4 variants through fitness valleys, whereas for other patients the data favored a gradual mode of emergence.

CONCLUSIONS

CXCR4 usage can emerge independently, in multiple lineages, along the course of HIV infection. The mode of emergence, i.e. gradual or through fitness valleys seems to depend on both virus and patient factors. Furthermore, our analyses suggest that, besides becoming dominant after population-level switches, minor proportions of X4 viruses might exist along the infection, perhaps even at early stages of it. The fate of these minor variants might depend on both viral and host factors.

Evolution of HIV-1 quasispecies and coreceptor use in cell reservoirs of patients on suppressive antiretroviral therapy

OBJECTIVES

To track changes in the V3 env region of HIV-1 quasispecies and determine virus coreceptor use in cell reservoirs of patients on long-term suppressive antiretroviral therapy (ART).

PATIENTS AND METHODS

Ten patients whose plasma viraemia had been suppressed for a median of 5.5 years were followed for 5 years. The V3 env regions of viruses in peripheral blood mononuclear cells were analysed by ultra-deep sequencing (UDS). HIV-1 tropism was predicted using the geno2pheno 5.75 algorithm and a phenotypic assay.

RESULTS

The UDS and phenotypic assay data were concordant for predicting HIV-1 tropism. CXCR4-using viruses detected by UDS accounted for 14.7%-76.5% of the virus populations in samples from five patients at enrolment. Five patients harboured pure R5 virus populations and no X4 viruses emerged during the 5 years. The selection pressures estimated by the dN/dS ratio were acting on the V3 region to produce diversification of the quasispecies in CXCR4-infected patients and purification of the quasispecies in R5-infected patients on effective ART.

CONCLUSIONS

UDS showed that the virus quasispecies in cell reservoirs of patients on long-term suppressive ART continued to evolve. CXCR4-using variants became more diversified. Analysis of the selection pressures on the virus quasispecies could provide a clearer picture of virus persistence in patients on effective ART.

An international multicenter study on HIV-1 drug resistance testing by 454 ultra-deep pyrosequencing

The detection of mutant spectra within the viral quasispecies is critical for therapeutic management of HIV-1 infections. Routine clinical application of ultrasensitive genotyping requires reproducibility and concordance within and between laboratories. The goal of the study was to evaluate a new protocol on HIV-1 drug resistance testing by 454 ultra-deep pyrosequencing (454-UDS) in an international multicenter study. Sixteen blinded HIV-1 subtype B samples were provided for 454-UDS as both RNA and cDNA with viral titers of 88,600-573,000 HIV-1 RNA copies/ml. Eight overlapping amplicons spanning protease (PR) codons 10-99 and reverse transcriptase (RT) codons 1-251 were generated using molecular barcoded primers. 454-UDS was performed using the 454 Life Sciences/Roche GS FLX platform. PR and RT sequences were analyzed using 454 Life Sciences Amplicon Variant Analyzer (AVA) software. Quantified variation data were analyzed for intra-laboratory reproducibility and inter-laboratory concordance. Routine population sequencing was performed using the ViroSeq HIV-1 genotyping system. Eleven laboratories and the reference laboratory 454 Life Sciences sequenced the HIV-1 sample set. Data presented are derived from seven laboratories and the reference laboratory since severe study protocol execution errors occurred in four laboratories leading to exclusion. The median sequencing depth across all sites was 1364 reads per position (IQR=809-2065). 100% of the ViroSeq-reported mutations were also detected by 454-UDS. Minority HIV-1 drug resistance mutations, defined as HIV-1 drug resistance mutations identified at frequencies of 1-25%, were only detected by 454-UDS. Analysis of 10 preselected majority and minority mutations were consistently found across sites. The analysis of drug-resistance mutations detected between 1 and 10% demonstrated high intra- and inter-laboratory consistency in frequency estimates for both RNA and prepared cDNA samples, indicating robustness of the method. HIV-1 drug resistance testing using 454 ultra-deep pyrosequencing results in an accurate and highly reproducible, albeit complex, approach to the analysis of HIV-1 mutant spectra, even at frequencies well below those detected by routine population sequencing.

Personalized HIV therapy to control drug resistance

The therapy of HIV patients is characterized by both the high genomic diversity of the virus population harbored by the patient and a substantial volume of therapy options. The virus population is unique for each patient and time point. The large number of therapy options makes it difficult to select an optimal or near optimal therapy, especially with therapy-experienced patients. In the past decade, computer-based support for therapy selection, which assesses the level of viral resistance against drugs has become a mainstay for HIV patients. We discuss the properties of available systems and the perspectives of the field.

Comparison of genotypic and phenotypic HIV type 1 tropism assay: results from the screening samples of Cenicriviroc Study 202, a randomized phase II trial in treatment-naive subjects

Cenicriviroc is a once-daily oral CCR5/CCR2 antagonist in development for treatment of HIV infection. CVC Study 202 (652-2-202; NCT01338883) excluded treatment-naive subjects demonstrated to harbor non-R5 (CXCR4-tropic or dual-mixed) tropic HIV-1 by either genotypic or phenotypic tropism testing. Here we compare the results of genotypic and phenotypic tropism testing in Study 202. A total of 304 subjects screened had paired genotypic and phenotypic results. Genotypic tropism testing (GTT) incorporated triplicate population sequencing using the geno2pheno algorithm and the PSSM algorithm, followed by ultradeep sequencing (UDS) for samples with R5 results. All samples were further evaluated with a phenotypic test, the enhanced-sensitivity Trofile assay (ESTA). Concordance between GTT and ESTA was 80% and increased to 84% when only geno2pheno was used for triplicate population sequencing. GTT (geno2pheno) classified 18% of the samples as non-R5 compared to 16% by ESTA. Only one-third of samples with non-R5 results by either test were classified as non-R5 by both tests. Median CD4((+)) cell counts were lower in patients with concordant non-R5 results by UDS and ESTA than in subjects with an R5 result by either assay (p=0.0004). UDS detected non-R5 virus in an additional 27/304 subjects (median 15% non-R5, interquartile range: 3.7-62%) with R5 results by ESTA. In conclusion, the geno2pheno algorithm improves concordance of GTT with a clinically validated phenotypic tropism assay as does the use of UDS. These findings provide support for recent guidelines indicating that genotypic tropism testing may be considered as an alternative to phenotypic testing.

Bioinformatic analysis of HIV-1 entry and pathogenesis

The evolution of human immunodeficiency virus type 1 (HIV-1) with respect to co-receptor utilization has been shown to be relevant to HIV-1 pathogenesis and disease. The CCR5-utilizing (R5) virus has been shown to be important in the very early stages of transmission and highly prevalent during asymptomatic infection and chronic disease. In addition, the R5 virus has been proposed to be involved in neuroinvasion and central nervous system (CNS) disease. In contrast, the CXCR4-utilizing (X4) virus is more prevalent during the course of disease progression and concurrent with the loss of CD4(+) T cells. The dual-tropic virus is able to utilize both co-receptors (CXCR4 and CCR5) and has been thought to represent an intermediate transitional virus that possesses properties of both X4 and R5 viruses that can be encountered at many stages of disease. The use of computational tools and bioinformatic approaches in the prediction of HIV-1 co-receptor usage has been growing in importance with respect to understanding HIV-1 pathogenesis and disease, developing diagnostic tools, and improving the efficacy of therapeutic strategies focused on blocking viral entry. Current strategies have enhanced the sensitivity, specificity, and reproducibility relative to the prediction of co-receptor use; however, these technologies need to be improved with respect to their efficient and accurate use across the HIV-1 subtypes. The most effective approach may center on the combined use of different algorithms involving sequences within and outside of the env-V3 loop. This review focuses on the HIV-1 entry process and on co-receptor utilization, including bioinformatic tools utilized in the prediction of co-receptor usage. It also provides novel preliminary analyses for enabling identification of linkages between amino acids in V3 with other components of the HIV-1 genome and demonstrates that these linkages are different between X4 and R5 viruses.

Using ultradeep pyrosequencing to study HIV-1 coreceptor usage in primary and dual infection

HIV-1 dual infection (DI) and CXCR4 (X4) coreceptor usage are associated with accelerated disease progression but frequency and dynamics of coreceptor usage during DI is unknown. Ultradeep sequencing was used to interrogate for DI and infer coreceptor usage in longitudinal blood samples of 102 subjects. At baseline, X4 usage was high (23 subjects harbored X4 variants) and was not associated with infection duration or DI. Coreceptor usage changed over time in 12 of 47 participants, and X4 usage emerged in 4 of 41 monoinfections vs 2 of 5 superinfections (P = .12), suggesting a weak statistical trend toward occurrence of superinfection and acquiring X4 usage.

Next-generation sequencing to assess HIV tropism

PURPOSE OF REVIEW

Deep sequencing of the V3 region of the HIV envelope gene can detect minority non-R5 variants in patients with high sensitivity and specificity. As next-generation sequencing approaches have matured, the clinical utility of deep sequencing for HIV tropism has entered the clinic. Accurate and sensitive tropism testing is essential for successful treatment with the CCR5 antagonist class of antiretrovirals.

RECENT FINDINGS

This review will focus on five aspects of next-generation sequencing for assessing HIV tropism: some background on the necessity of deep sequencing versus other tropism methods; the methodological process of 454 sequencing and analysis; other next-generation sequencing technologies; the diagnostic performance of deep sequencing relative to other tropism assays; and the use of deep sequencing in clinical practice.

SUMMARY

This method has emerged quickly as both a research and clinical tool because of its high concordance with commonly used phenotypic tropism assays and its ability to predict virological response to CCR5 antagonist-containing regimens.

Sexual transmission of predicted CXCR4-tropic HIV-1 likely originating from the source partner's seminal cells

We present a case of sexual transmission of HIV-1 predicted to have CXCR4-tropism during male-to-male sexual exposure. Phylogenetic analyses exclude cell-free virus in the seminal plasma of the source partner and possibly point to the seminal cells as the origin of the transmission event.

Higher prevalence of predicted X4-tropic strains in perinatally infected older children with HIV-1 subtype C in India

BACKGROUND

Coreceptor switch from CCR5 to CXCR4 is considered to be less common in HIV-1 subtype C even in advanced stages of infection. In this study, we have examined viral genotypic coreceptor tropism and its clinical, virological, and host genetic determinants among perinatally infected children in India.

METHODS

Genotypic coreceptor tropism analysis was conducted on env V3 sequences using Geno2pheno with a threshold of 10% false-positive rate. A total of 473 sequences were obtained from 72 isolates amplified from children aged 2-17 years. Factors associated with viral tropism in subtype C infections were studied using logistic regression.

RESULTS

Among the samples, 98.6% (71 of 72) were HIV-1 subtype C. Coreceptor tropism analysis determined 81.7% (58 of 71) as R5 tropic, 9.9% (7 of 71) as X4 tropic, and 8.5% (6 of 71) as R5/X4 tropic or dual-tropic HIV-1 strains. Children with X4 or R5/X4 strains were more likely to be older than those with R5-tropic strains (P < 0.05), have lower CD4 counts (P < 0.05), and have viral populations with greater intrapopulation viral divergence (P < 0.01). Older age was a significant independent predictor for X4 or R5/X4 tropism in these children (P < 0.05). None were identified as being heterozygous or homozygous for the CCR5[INCREMENT]32 deletion.

CONCLUSIONS

The high prevalence of X4 and R5/X4 tropic strains among older perinatally infected children with HIV-1 subtype C in India indicate that this phenomenon is not uncommon as previously thought and suggest that coreceptor transition can occur with longer duration of infection and greater disease progression in this population of perinatally infected children living with HIV-1 subtype C.

Geno2pheno[454]: a Web server for the prediction of HIV-1 coreceptor usage from next-generation sequencing data

Inferring HIV-1 coreceptor usage from a genotype is becoming more and more important for the appropriate treatment of long-term patients. While results are already encouraging where standard bulk-nucleic acid sequencing methods are used, they are limited with respect to the detection of minor variants. In contrast, next-generation sequencing methods (ultradeep sequencing, pyrosequencing) are capable of sequencing virus quasispecies at very low quantities. However, as well as being very expensive, these methods generate vast amounts of data such that sequence analysis has to be automated by computer assistance. Here, we describe the geno2pheno[454] system which handles all processing and prediction steps involved in the prediction of coreceptor usage from massively parallel sequencing data. The system is split into a JAVA preprocessor which is run locally on the client side and a Web server which generates the prediction results. Predictions are based on the same prediction method as used in the geno2pheno[coreceptor] tool.

Ultra-deep sequencing for the analysis of viral populations

Next-generation sequencing allows for cost-effective probing of virus populations at an unprecedented level of detail. The massively parallel sequencing approach can detect low-frequency mutations and it provides a snapshot of the entire virus population. However, analyzing ultra-deep sequencing data obtained from diverse virus populations is challenging because of PCR and sequencing errors and short read lengths, such that the experiment provides only indirect evidence of the underlying viral population structure. Recent computational and statistical advances allow for accommodating some of the confounding factors, including methods for read error correction, haplotype reconstruction, and haplotype frequency estimation. With these methods ultra-deep sequencing can be more reliably used to analyze, in a quantitative manner, the genetic diversity of virus populations.