Appraising the performance of genotyping tools in the prediction of coreceptor tropism in HIV-1 subtype C viruses

Performance of Geno2Pheno[coreceptor] to infer coreceptor use in human immunodeficiency virus type 1 (HIV-1) subtype A

BACKGROUND

Assessment of human immunodeficiency virus type 1 (HIV-1) coreceptor usage is required prior to treatment with the CCR5 antagonist maraviroc to exclude the presence of CXCR4-using (X4) strains. Genotype-based interpretation systems are mostly designed on subtype B and have been reported to be less accurate for subtype A/CRF02_AG.

OBJECTIVES

To evaluate the performance of the widely used Geno2Pheno[coreceptor] (G2P[c]) algorithm for prediction of coreceptor usage with subtype A/CRF02_AG vs. subtype B.

STUDY DESIGN

Co-receptor tropism of 24 subtype A/CRF02_AG and 24 subtype B viruses was measured phenotypically by a homebrew single-cycle assay and genotypically by using G2P[c]. Samples with discrepant genotype-phenotype results were analyzed by next generation sequencing (NGS) and interpreted by the NGS Geno2Pheno algorithm (G2P[454]).

RESULTS

At 10% false positive rate (FPR), the G2P[c]/phenotype discordance rate was 12.5% (n = 3) for subtype A/CRF02_AG and 8.3% (n = 2) for subtype B. Minority X4 species escaping detection by bulk sequencing but documented by NGS explained the two subtype B and possibly one subtype A/CRF02_AG discordant case. The other two subtype A/CRF02_AG miscalled by G2P[c] could be explained by X4 overcalling at borderline FPR and/or by algorithm failure.

DISCUSSION

Our study did not demonstrate relevantly higher G2P[c] inaccuracy with subtype A/CRF02_AG with respect to subtype B. Genotype/phenotype discordances can be due to different reasons, including but not limited to, algorithm inaccuracy. Very large genotype/phenotype correlation panels are required to detect and explain the reason for any consistent difference in genotypic tropism prediction for subtype A/CRF02_AG vs. subtype B.

Prediction of coreceptor usage by five bioinformatics tools in a large Ethiopian HIV-1 subtype C cohort

Background

Genotypic tropism testing (GTT) has been developed largely on HIV-1 subtype B. Although a few reports have analysed the utility of GTT in other subtypes, more studies using HIV-1 subtype C (HIV-1C) are needed, considering the huge contribution of HIV-1C to the global epidemic.

Methods

Plasma was obtained from 420 treatment-naïve HIV-1C infected Ethiopians recruited 2009–2011. The V3 region was sequenced and the coreceptor usage was predicted by five tools: Geno2Pheno clinical–and clonal–models, PhenoSeq-C, C-PSSM and Raymond’s algorithm. The impact of baseline tropism on antiretroviral treatment (ART) outcome was evaluated.

Results

Of 352 patients with successful baseline V3 sequences, the proportion of predicted R5 virus varied between the methods by 12.5% (78.1%-90.6%). However, only 58.2% of the predictions were concordant and only 1.7% were predicted to be X4-tropic across the five methods. Compared pairwise, the highest concordance was between C-PSSM and Geno2Pheno clonal (86.4%). In bivariate intention to treat (ITT) analysis, R5 infected patients achieved treatment success more frequently than X4 infected at month six as predicted by Geno2Pheno clinical (77.8% vs 58.7%, P = 0.004) and at month 12 by C-PSSM (61.9% vs 46.6%, P = 0.038). However, in the multivariable analysis adjusted for age, gender, baseline CD4 and viral load, only tropism as predicted by C-PSSM showed an impact on month 12 (P = 0.04, OR 2.47, 95% CI 1.06–5.79).

Conclusion

Each of the bioinformatics models predicted R5 tropism with comparable frequency but there was a large discordance between the methods. Baseline tropism had an impact on outcome of first line ART at month 12 in multivariable ITT analysis but only based on prediction by C-PSSM which thus possibly could be used for predicting outcome of ART in HIV-1C infected Ethiopians.

Partial HIV C2V3 envelope sequence analysis reveals association of coreceptor tropism, envelope glycosylation and viral genotypic variability among Kenyan patients on HAART

Background

HIV-1 is highly variable genetically and at protein level, a property it uses to subvert antiviral immunity and treatment. The aim of this study was to assess if HIV subtype differences were associated with variations in glycosylation patterns and co-receptor tropism among HAART patients experiencing different virologic treatment outcomes.

Methods

A total of 118 HIV env C2V3 sequence isolates generated previously from 59 Kenyan patients receiving highly active antiretroviral therapy (HAART) were examined for tropism and glycosylation patterns. For analysis of Potential N-linked glycosylation sites (PNGs), amino acid sequences generated by the NCBI’s Translate tool were applied to the HIVAlign and the N-glycosite tool within the Los Alamos Database. Viral tropism was assessed using Geno2Pheno (G2P), WebPSSM and Phenoseq platforms as well as using Raymond’s and Esbjörnsson’s rules. Chi square test was used to determine independent variables association and ANOVA applied on scale variables.

Results

At respective False Positive Rate (FPR) cut-offs of 5% (p = 0.045), 10% (p = 0.016) and 20% (p = 0.005) for CXCR4 usage within the Geno2Pheno platform, HIV-1 subtype and viral tropism were significantly associated in a chi square test. Raymond’s rule (p = 0.024) and WebPSSM (p = 0.05), but not Phenoseq or Esbjörnsson showed significant associations between subtype and tropism. Relative to other platforms used, Raymond’s and Esbjörnsson’s rules showed higher proportions of X4 variants, while WebPSSM resulted in lower proportions of X4 variants across subtypes. The mean glycosylation density differed significantly between subtypes at positions, N277 (p = 0.034), N296 (p = 0.036), N302 (p = 0.034) and N366 (p = 0.004), with HIV-1D most heavily glycosylated of the subtypes. R5 isolates had fewer PNGs than X4 isolates, but these differences were not significant except at position N262 (p = 0.040). Cell-associated isolates from virologic treatment success subjects were more glycosylated than cell-free isolates from virologic treatment failures both for the NXT (p = 0.016), and for all the patterns (p = 0.011).

Conclusion

These data reveal significant associations of HIV-1 subtype diversity, viral co-receptor tropism, viral suppression and envelope glycosylation. These associations have important implications for designing therapy and vaccines against HIV. Heavy glycosylation and preference for CXCR4 usage of HIV-1D may explain rapid disease progression in patients infected with these strains.

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.

Features of Recently Transmitted HIV-1 Clade C Viruses that Impact Antibody Recognition: Implications for Active and Passive Immunization

The development of biomedical interventions to reduce acquisition of HIV-1 infection remains a global priority, however their potential effectiveness is challenged by very high HIV-1 envelope diversity. Two large prophylactic trials in high incidence, clade C epidemic regions in southern Africa are imminent; passive administration of the monoclonal antibody VRC01, and active immunization with a clade C modified RV144-like vaccines. We have created a large representative panel of C clade viruses to enable assessment of antibody responses to vaccines and natural infection in Southern Africa, and we investigated the genotypic and neutralization properties of recently transmitted clade C viruses to determine how viral diversity impacted antibody recognition. We further explore the implications of these findings for the potential effectiveness of these trials. A panel of 200 HIV-1 Envelope pseudoviruses was constructed from clade C viruses collected within the first 100 days following infection. Viruses collected pre-seroconversion were significantly more resistant to serum neutralization compared to post-seroconversion viruses (p = 0.001). Over 13 years of the study as the epidemic matured, HIV-1 diversified (p = 0.0009) and became more neutralization resistant to monoclonal antibodies VRC01, PG9 and 4E10. When tested at therapeutic levels (10ug/ml), VRC01 only neutralized 80% of viruses in the panel, although it did exhibit potent neutralization activity against sensitive viruses (IC50 titres of 0.42 μg/ml). The Gp120 amino acid similarity between the clade C panel and candidate C-clade vaccine protein boosts (Ce1086 and TV1) was 77%, which is 8% more distant than between CRF01_AE viruses and the RV144 CRF01_AE immunogen. Furthermore, two vaccine signature sites, K169 in V2 and I307 in V3, associated with reduced infection risk in RV144, occurred less frequently in clade C panel viruses than in CRF01_AE viruses from Thailand. Increased resistance of pre-seroconversion viruses and evidence of antigenic drift highlights the value of using panels of very recently transmitted viruses and suggests that interventions may need to be modified over time to track the changing epidemic. Furthermore, high divergence such as that observed in the older clade C epidemic in southern Africa may impact vaccine efficacy, although the correlates of infection risk are yet to be defined in the clade C setting. Findings from this study of acute/early clade C viruses will aid vaccine development, and enable identification of new broad and potent antibodies to combat the HIV-1 C-clade epidemic in southern Africa.

Risk of breast cancer with CXCR4-using HIV defined by V3 loop sequencing

OBJECTIVE

Evaluate the risk of female breast cancer associated with HIV-CXCR4 (X4) tropism as determined by various genotypic measures.

METHODS

A breast cancer case-control study, with pairwise comparisons of tropism determination methods, was conducted. From the Women's Interagency HIV Study repository, one stored plasma specimen was selected from 25 HIV-infected cases near the breast cancer diagnosis date and 75 HIV-infected control women matched for age and calendar date. HIV-gp120 V3 sequences were derived by Sanger population sequencing (PS) and 454-pyro deep sequencing (DS). Sequencing-based HIV-X4 tropism was defined using the geno2pheno algorithm, with both high-stringency DS [false-positive rate (3.5) and 2% X4 cutoff], and lower stringency DS (false-positive rate, 5.75 and 15% X4 cutoff). Concordance of tropism results by PS, DS, and previously performed phenotyping was assessed with kappa (κ) statistics. Case-control comparisons used exact P values and conditional logistic regression.

RESULTS

In 74 women (19 cases, 55 controls) with complete results, prevalence of HIV-X4 by PS was 5% in cases vs 29% in controls (P = 0.06; odds ratio, 0.14; confidence interval: 0.003 to 1.03). Smaller case-control prevalence differences were found with high-stringency DS (21% vs 36%, P = 0.32), lower stringency DS (16% vs 35%, P = 0.18), and phenotyping (11% vs 31%, P = 0.10). HIV-X4 tropism concordance was best between PS and lower stringency DS (93%, κ = 0.83). Other pairwise concordances were 82%-92% (κ = 0.56-0.81). Concordance was similar among cases and controls.

CONCLUSIONS

HIV-X4 defined by population sequencing (PS) had good agreement with lower stringency DS and was significantly associated with lower odds of breast cancer.

Performance of Genotypic Tools for Prediction of Tropism in HIV-1 Subtype C V3 Loop Sequences

Currently, there is no consensus on the genotypic tools to be used for tropism analysis in HIV-1 subtype C strains. Thus, the aim of the study was to evaluate the performance of the different V3 loop-based genotypic algorithms available. We compiled a dataset of 645 HIV-1 subtype C V3 loop sequences of known coreceptor phenotypes (531 R5-tropic/non-syncytium-inducing and 114 X4-tropic/R5X4-tropic/syncytium-inducing sequences) from the Los Alamos database (http://www.hiv.lanl.gov/) and previously published literature. Coreceptor usage was predicted based on this dataset using different software-based machine-learning algorithms as well as simple classical rules. All the sophisticated machine-learning methods showed a good concordance of above 85%. Geno2Pheno (false-positive rate cutoff of 5-15%) and CoRSeqV3-C were found to have a high predicting capability in determining both HIV-1 subtype C X4-tropic and R5-tropic strains. The current sophisticated genotypic tropism tools based on V3 loop perform well for tropism prediction in HIV-1 subtype C strains and can be used in clinical settings.

Profile of HIV type 1 coreceptor tropism among Kenyan patients from 2009 to 2010

A switch of HIV coreceptor usage from CCR5 to CXCR4 occurs in AIDS pathogenesis and may play a critical role in the use of entry inhibitors. To determine the potential usefulness of maraviroc and other CCR5 antagonists among drug-naive and experienced patients in Kenya, the env-C2-V3 gene was successfully sequenced in samples from 176 (98 men and 78 female) consenting subjects between January 2009 and December 2012. In silico CPSSM, webPSSM/, and (ds) Kernel tools were used in predicting coreceptor usage. On the basis of the env V3 loop sequences, 84.1% (148) were reported with R-5 tropism, 4.5% (5) were dual tropic, while 13.4% (23) were of X4 tropism. However, similar to previous studies conducted in Kenya on genetic diversity, HIV-1 subtype A1 (73.9%; 130/176) still remains the most dominant subtype. The high levels of R5 tropism among the studied Kenyan infected populations suggested the potential use of CCR5 antagonists as new therapeutic options in Kenya.