Correlation between genotypic predictions based on V3 sequences and phenotypic determination of HIV-1 tropism

Unique profile of predominant CCR5-tropic in CRF07_BC HIV-1 infections and discovery of an unusual CXCR4-tropic strain

CRF07_BC is one of the most prevalent HIV-1 strains in China, which contributes over one-third of the virus transmissions in the country. In general, CRF07_BC is associated with slower disease progression, while the underlying mechanisms remain unclear. Our study focused on envelope proteins (Env) and its V3 loop which determine viral binding to co-receptors during infection of cells. We studied a large dataset of 3,937 env sequences in China and found that CRF07_BC had a unique profile of predominantly single CCR5 tropism compared with CCR5 and CXCR4 dual tropisms in other HIV-1 subtypes. The percentages of the CXCR4-tropic virus in B (3.7%) and CRF01_AE (10.4%) infection are much higher than that of CRF07_BC (0.1%), which is supported by median false-positive rates (FPRs) of 69.8%, 25.5%, and 13.4% for CRF07_BC, B, and CRF01_AE respectively, with a cutoff FPR for CXCR4-tropic at 2%. In this study, we identified the first pure CXCR4-tropic virus from one CRF07_BC-infected patient with an extremely low CD4⁺T cell count (7 cells/mm³). Structural analysis found that the V3 region of this virus has the characteristic 7T and 25R and a substitution of conserved “GPGQ” crown motif for “GPGH”. This study provided compelling evidence that CRF07_BC has the ability to evolve into CXCR4 strains. Our study also lay down the groundwork for studies on tropism switch, which were commonly done for other HIV-1 subtypes, for the long-delayed CRF07_BC.

In vitro susceptibility to fostemsavir is not affected by long-term exposure to antiviral therapy in MDR HIV-1-infected patients

OBJECTIVES

Fostemsavir is the prodrug of the HIV-1 attachment inhibitor temsavir and is currently under clinical assessment in heavily treatment-experienced patients with limited therapeutic options. We evaluated the genotypic and phenotypic susceptibility to temsavir in a panel of samples collected from patients harbouring MDR strains enrolled in the Italian PRESTIGIO Registry.

METHODS

Plasma samples from 24 patients were used for HIV-1 gp120 sequencing, while viral tropism and susceptibility to temsavir were assessed through a homemade phenotypic assay with pseudotyped viruses expressing patient-derived Env protein.

RESULTS

Of the 24 patients enrolled, 18 (75%) were male, median (IQR) age was 55 years (52-61), time since HIV-1 diagnosis was 27 years (24-30), time on ART was 26 years (23-27) and 11 (46%) had a previous AIDS diagnosis. Exposure to entry inhibitors (maraviroc and/or enfuvirtide) had occurred in 19 (79%) patients. Among 23/24 gp120 sequences obtained, temsavir resistance-associated mutations (RAMs) were detected in three cases (two M426L and one S375N). Pseudotyped viruses were obtained from 23/24 samples and viral tropism was CXCR4-tropic, CCR5-tropic and dual/mixed-tropic in six, nine and eight cases, respectively. Phenotypic susceptibility to temsavir was comparable to the reference WT viruses NL4-3 and AD8 in all samples, irrespective of RAMs. Viral tropism and exposure to entry inhibitors did not impact temsavir susceptibility.

CONCLUSIONS

These data support the use of fostemsavir as a valuable therapy option in patients harbouring MDR virus. The role of laboratory testing in optimal screening of patients eligible for fostemsavir treatment remains to be investigated.

THETA: a new genotypic approach for predicting HIV-1 CRF02-AG coreceptor usage

MOTIVATION

The circulating recombinant form of HIV-1 CRF02-AG is the most frequent non-B subtype in Europe. Anti-HIV therapy and pathophysiological studies on the impact of HIV-1 tropism require genotypic determination of HIV-1 tropism for non-B subtypes. But genotypic approaches based on analysis of the V3 envelope region perform poorly when used to determine the tropism of CRF02-AG. We, therefore, designed an algorithm based on information from the gp120 and gp41 ectodomain that better predicts the tropism of HIV-1 subtype CRF02-AG.

RESULTS

We used a bio-statistical method to identify the genotypic determinants of CRF02-AG coreceptor use. Toulouse HIV Extended Tropism Algorithm (THETA), based on a Least Absolute Shrinkage and Selection Operator method, uses HIV envelope sequence from phenotypically characterized clones. Prediction of R5X4/X4 viruses was 86% sensitive and that of R5 viruses was 89% specific with our model. The overall accuracy of THETA was 88%, making it sufficiently reliable for predicting the tropism of subtype CRF02-AG sequences.

AVAILABILITY AND IMPLEMENTATION

Binaries are freely available for download at https://github.com/viro-tls/THETA. It was implemented in Matlab and supported on MS Windows platform. The sequence data used in this work are available from GenBank under the accession numbers MK618182-MK618417.

Clinical use, efficacy, and durability of maraviroc for antiretroviral therapy in routine care: A European survey

OBJECTIVES

The study aimed to survey maraviroc use and assess effectiveness and durability of maraviroc-containing antiretroviral treatment (ART) in routine practice across Europe.

METHODS

Data were retrieved from 26 cohorts in 8 countries comprising adults who started maraviroc in 2005-2016 and had ≥1 follow-up visit. Available V3 sequences were re-analysed centrally for tropism determination by geno2pheno[coreceptor]. Treatment failure (TF) was defined as either virological failure (viral load >50 copies/mL) or maraviroc discontinuation for any reason over 48 weeks. Predictors of TF were explored by logistic regression analysis. Time to maraviroc discontinuation was estimated by Kaplan-Meier survival analysis.

RESULTS

At maraviroc initiation (baseline), among 1,381 patients, 67.1% had experienced ≥3 ART classes and 45.6% had a viral load <50 copies/mL. Maraviroc was occasionally added to the existing regimen as a single agent (7.3%) but it was more commonly introduced alongside other new agents, and was often (70.4%) used with protease inhibitors. Accompanying drugs comprised 1 (40.2%), 2 (48.6%) or ≥3 (11.2%) ART classes. Among 1,273 patients with available tropism data, 17.6% showed non-R5 virus. Non-standard maraviroc use also comprised reported once daily dosing (20.0%) and a total daily dose of 150mg (12.1%). Over 48 weeks, 41.4% of patients met the definition of TF, although the 1-year estimated retention on maraviroc was 82.1% (95% confidence interval 79.9-84.2). Among 1,010 subjects on maraviroc at week 48, the viral load was >50 copies/mL in 19.9% and >200 copies/mL in 10.7%. Independent predictors of TF comprised a low nadir CD4 count, a detectable baseline viral load, previous PI experience, non-R5 tropism, having ≥3 active drugs in the accompanying regimen, and a more recent calendar year of maraviroc initiation.

CONCLUSIONS

This study reports on the largest observation cohort of patients who started maraviroc across 8 European countries. In this overall highly treatment-experienced population, with a small but appreciable subset that received maraviroc outside of standard treatment guidelines, maraviroc was safe and reasonably effective, with relatively low rates of discontinuation over 48 weeks and only 2 cases of serum transaminase elevations reported as reasons for discontinuation.

Long-term evolution of transmitted CXCR4-using HIV-1 under effective antiretroviral therapy

OBJECTIVE

To study the long-term evolution of the transmitted CXCR4-using viruses. CCR5-using viruses (R5 viruses) predominate during primary HIV-1 infections (PHI) while CXCR4-using viruses are isolated in less than 10% of PHI.

DESIGN

Six patients infected with an R5X4 virus, detected by a sensitive phenotypic assay during PHI, were matched with six patients infected with a pure R5 virus for sex, Fiebig stage, time of antiretroviral initiation and duration of follow-up.

METHODS

We used MiSeq ultra-deep sequencing to determine the composition of the virus quasispecies during PHI and at the end of follow-up (median time of follow-up: 12.5 years).

RESULTS

X4 viruses were detected by genetic analysis in three of six samples from the R5X4 group, accounting for 1.3-100% of the virus quasispecies, during PHI, and in four of six samples (accounting for 6.7-100%) at the end of follow-up. No X4 virus was detected in the R5 group during PHI and in only one patient (accounting for 1.2%) at the end of follow-up. The complexity of the virus quasispecies at the stage of PHI was higher in the R5X4 group than in the R5 group. Complexity increased from PHI to the end of follow-up in the R5 group but remained stable in the R5X4 group.

CONCLUSION

CXCR4-using viruses persisted in the peripheral blood mononuclear cells of several patients on suppressive antiretroviral therapy for a median duration of 12.5 years after PHI. The genetic complexity of HIV-1 evolved differently post-PHI in patients infected with R5X4 viruses from those infected with R5 viruses.

Coreceptor Tropism and Maraviroc Sensitivity of Clonally Derived Ethiopian HIV-1C Strains Using an in-house Phenotypic Assay and Commonly Used Genotypic Methods

OBJECTIVES

Genotypic Tropism Testing (GTT) tools are generally developed based on HIV-1 subtype B (HIV-1B) and used for HIV-1C as well but with a large discordance of prediction between different methods. We used an established phenotypic assay for comparison with GTT methods and for the determination of in vitro maraviroc sensitivity of pure R5-tropic and dual-tropic HIV-1C.

METHODS

Plasma was obtained from 58 HIV-1C infected Ethiopians. Envgp120 was cloned into a luciferase tagged NL4-3 plasmid. Phenotypic tropism was determined by in house method and the V3 sequences were analysed by five GTT methods. In vitro maraviroc sensitivity of R5-tropic and dual-tropic isolates were compared in the TZMbl cell-line.

RESULTS

The phenotypes were classified as R5 in 92.4% and dual tropic (R5X4) in 7.6% of 79 clones. The concordance between phenotype and genotype ranged from 64.7% to 84.3% depending on the GTT method. Only 46.9% of the R5 phenotypes were predicted as R5 by all GTT tools while R5X4 phenotypes were predicted as X4 by four methods, but not by Raymond's method. All six tested phenotypic R5 clones, as well as five of six of dual tropic clones, showed a dose response to maraviroc.

CONCLUSION

There is a high discordance between GTT methods, which underestimates the presence of R5 and overestimates X4 strains compared to a phenotypic assay. Currently available GTT algorithms should be further improved for tropism prediction in HIV-1C. Maraviroc has an in vitro activity against most HIV-1C viruses and could be considered as an alternative regimen in individuals infected with CCR5-tropic HIV-1C viruses.

Transmission patterns of HIV-1 non-R5 strains in Poland

HIV-1 env sequencing enables predictions of viral coreceptor tropism and phylogenetic investigations of transmission events. The aim of the study was to estimate the contribution of non-R5 strains to the viral spread in Poland. Partial proviral env sequences were retrieved from baseline blood samples of patients with newly diagnosed HIV-1 infection between 2008-2014, including 46 patients with recent HIV-1 infection (RHI), and 246 individuals with long-term infection (LTHI). These sequences were subjected to the genotypic coreceptor tropism predictions and phylogenetic analyses to identify transmission clusters. Overall, 27 clusters with 57 sequences (19.5%) were detected, including 15 sequences (26.3%) from patients with RHI. The proportion of non-R5 strains among all study participants was 23.3% (68/292), and was comparable between patients with RHI and LTHI (11/46, 23.9% vs 57/246, 23.2%; p = 1.000). All 11 patients with non-R5 strains and RHI were men having sex with men (MSM). Among these patients, 4 had viral sequences grouped within phylogenetic cluster with another sequence of non-R5 strain obtained from patient with LTHI, indicating potential acquisition of non-R5 HIV-1 for at least 4/46 (8.7%) patients with RHI. We were unable to confirm the contribution of patients with RHI to the forward transmission of non-R5 strains, but a relatively high proportion of non-R5 strains among them deserves attention due to the limited susceptibility to CCR5 antagonists.

Accurate quantification of CXCR4-using HIV-1 variants by Illumina deep-sequencing

: Next-generation sequencing is a sensitive method for determining HIV-1 tropism but there is a lack of data on the quantification of X4 variants. We evaluated MiSeq and 454 GS-Junior platforms for determining HIV-1 tropism and for quantifying X4 variants. Both platforms were 93% concordant for determining HIV-1 tropism and correlated well for determining the proportion of X4 variants (Spearman correlation, ρ = 0.748; P < 0.0001). MiSeq Illumina sequencing seems to be well adapted for characterizing X4-containing samples.

CD4 T cell decline following HIV seroconversion in individuals with and without CXCR4-tropic virus

Background

The natural clinical and immunological courses following HIV seroconversion with CXCR4-tropic or dual-mixed (X4/DM) viruses are controversial. We compared spontaneous immunological outcome in patients harbouring an X4/DM virus at the time of seroconversion with those harbouring a CCR5-tropic (R5) virus.

Methods

Data were included from patients participating in CASCADE, a large cohort collaboration of HIV seroconverters, with ≥2 years of follow-up since seroconversion. The HIV envelope gene was sequenced from frozen plasma samples collected at enrolment, and HIV tropism was determined using Geno2Pheno (false-positive rate 10%). The spontaneous CD4 T cell evolution was compared by modelling CD4 kinetics using linear mixed-effects models with random intercept and random slope.

Results

A total of 1387 patients were eligible. Median time between seroconversion and enrolment was 1 month (range 0-3). At enrolment, 202 of 1387 (15%) harboured an X4/DM-tropic virus. CD4 decrease slopes were not significantly different according to HIV-1 tropism during the first 30 months after seroconversion. No marked change in these results was found after adjusting for age, year of seroconversion and baseline HIV viral load. Time to antiretroviral treatment initiation was not statistically different between patients harbouring an R5 (20.76 months) and those harbouring an X4/DM-tropic virus (22.86 months, logrank test P = 0.32). Conclusions: In this large cohort collaboration, 15% of the patients harboured an X4/DM virus close to HIV seroconversion. Patients harbouring X4/DM-tropic viruses close to seroconversion did not have an increased risk of disease progression, estimated by the decline in CD4 T cell count or time to combined ART initiation.

Pace of Coreceptor Tropism Switch in HIV-1-Infected Individuals after Recent Infection

HIV-1 entry into target cells influences several aspects of HIV-1 pathogenesis, including viral tropism, HIV-1 transmission and disease progression, and response to entry inhibitors. The evolution from CCR5- to CXCR4-using strains in a given human host is still unpredictable. Here we analyzed timing and predictors for coreceptor evolution among recently HIV-1-infected individuals. Proviral DNA was longitudinally evaluated in 66 individuals using Geno2pheno_[coreceptor] Demographics, viral load, CD4⁺ and CD8⁺ T cell counts, CCR5Δ32 polymorphisms, GB virus C (GBV-C) coinfection, and HLA profiles were also evaluated. Ultradeep sequencing was performed on initial samples from 11 selected individuals. A tropism switch from CCR5- to CXCR4-using strains was identified in 9/49 (18.4%) individuals. Only a low baseline false-positive rate (FPR) was found to be a significant tropism switch predictor. No minor CXCR4-using variants were identified in initial samples of 4 of 5 R5/non-R5 switchers. Logistic regression analysis showed that patients with an FPR of >40.6% at baseline presented a stable FPR over time whereas lower FPRs tend to progressively decay, leading to emergence of CXCR4-using strains, with a mean evolution time of 27.29 months (range, 8.90 to 64.62). An FPR threshold above 40.6% determined by logistic regression analysis may make it unnecessary to further determine tropism for prediction of disease progression related to emergence of X4 strains or use of CCR5 antagonists. The detection of variants with intermediate FPRs and progressive FPR decay over time not only strengthens the power of Geno2pheno in predicting HIV tropism but also indirectly confirms a continuous evolution from earlier R5 variants toward CXCR4-using strains.IMPORTANCE The introduction of CCR5 antagonists in the antiretroviral arsenal has sparked interest in coreceptors utilized by HIV-1. Despite concentrated efforts, viral and human host features predicting tropism switch are still poorly understood. Limited longitudinal data are available to assess the influence that these factors have on predicting tropism switch and disease progression. The present study describes longitudinal tropism evolution in a group of recently HIV-infected individuals to determine the prevalence and potential correlates of tropism switch. We demonstrated here that a low baseline FPR determined by the Geno2pheno_[coreceptor] algorithm can predict tropism evolution from CCR5 to CXCR4 coreceptor use.

HIV-1 coreceptor usage in perinatally infected Thai children

HIV-1 coreceptor usage in children who were born to HIV-1 infected mothers in Thailand is not well characterized. Here, the prevalence of coreceptor usage and genotype among HIV-1 infected children in Thailand were observed. Proviral DNA from 284 HIV-1 infected children who received HIV-1 early infant diagnosis between 2007 and 2013 under the National AIDS Program were studied. Genotypic tropism testing was performed based on amplification of the V3 region in a triplicate nested-PCR following by DNA sequencing. HIV-1 coreceptor usage was determined using Geno2pheno_[coreceptor] with a false positive rate of 10%. Samples from 267 children were successfully amplified and coreceptor usage could be determined. Two hundred and thirty-seven (89%) children were infected with CRF01_AE, 29 (11%) were subtype B and 1 was subtype C. CCR5-using variants were found in 148 (55%) children and CXCR4-using variants were observed in 119 (45%) children. No significant differences in coreceptor usage and age, gender, signs of HIV infection, children's or maternal ARV receiving were observed. The only significant difference was found in N-linked glycosylation characteristic. This evidence showed that X4 viruses can be highly observed at an early age of children which has important clinical implications and may limit usage of CCR5 antagonist family.

Performance comparison of next-generation sequencing platforms for determining HIV-1 coreceptor use

The coreceptor used by HIV-1 must be determined before a CCR5 antagonist, part of the arsenal of antiretroviral drugs, is prescribed because viruses that enter cells using the CXCR4 coreceptor are responsible for treatment failure. HIV-1 tropism is also correlated with disease progression and so must be determined for virological studies. Tropism can be determined by next-generation sequencing (NGS), but not all of these new technologies have been fully validated for use in clinical practice. The Illumina NGS technology is used in many laboratories but its ability to predict HIV-1 tropism has not been evaluated while the 454 GS-Junior (Roche) is used for routine diagnosis. The genotypic prediction of HIV-1 tropism is based on sequencing the V3 region and interpreting the results with an appropriate algorithm. We compared the performances of the MiSeq (Illumina) and 454 GS-Junior (Roche) systems with a reference phenotypic assay. We used clinical samples for the NGS tropism predictions and assessed their ability to quantify CXCR4-using variants. The data show that the Illumina platform can be used to detect minor CXCR4-using variants in clinical practice but technical optimization are needed to improve quantification.

No selection of CXCR4-using variants in cell reservoirs of dual-mixed HIV-infected patients on suppressive maraviroc therapy

We used ultradeep sequencing to investigate the evolution of the frequency of CXCR4-using viruses in the peripheral blood mononuclear cells of 22 patients infected with both CCR5 and CXCR4-using viruses treated with the CCR5 antagonist maraviroc for 24 weeks and a stable antiviral therapy. The mean CXCR4-using virus frequency in peripheral blood mononuclear cells was 59% before maraviroc intensification and 52% after 24 weeks of effective treatment, indicating no selection by maraviroc.

Position-specific automated processing of V3 env ultra-deep pyrosequencing data for predicting HIV-1 tropism

HIV-1 coreceptor usage must be accurately determined before starting CCR5 antagonist-based treatment as the presence of undetected minor CXCR4-using variants can cause subsequent virological failure. Ultra-deep pyrosequencing of HIV-1 V3 env allows to detect low levels of CXCR4-using variants that current genotypic approaches miss. However, the computation of the mass of sequence data and the need to identify true minor variants while excluding artifactual sequences generated during amplification and ultra-deep pyrosequencing is rate-limiting. Arbitrary fixed cut-offs below which minor variants are discarded are currently used but the errors generated during ultra-deep pyrosequencing are sequence-dependant rather than random. We have developed an automated processing of HIV-1 V3 env ultra-deep pyrosequencing data that uses biological filters to discard artifactual or non-functional V3 sequences followed by statistical filters to determine position-specific sensitivity thresholds, rather than arbitrary fixed cut-offs. It allows to retain authentic sequences with point mutations at V3 positions of interest and discard artifactual ones with accurate sensitivity thresholds.

Envelope gene evolution and HIV-1 neuropathogenesis

In the era of combined antiretroviral therapy (cART), HIV-associated neurocognitive disorders (HAND) account for 40 to 56% of all HIV⁺ cases. During the acute stage of HIV-1 infection (<6 months), the virus invades and replicates within the central nervous system (CNS). Compared to peripheral tissues, the local CNS cell population expresses distinct levels of chemokine receptors, which levels exert selective pressure on the invading virus. HIV-1 envelope (env) sequences recovered from the brains and cerebrospinal fluid (CSF) of neurocognitively impaired HIV⁺ subjects often display higher nucleotide variability as compared to non-impaired HIV⁺ subjects. Specifically, env evolution provides HIV-1 with the strategies to evade host immune response, to reduce chemokine receptor dependence, to increase co-receptor binding efficiency, and to potentiate neurotoxicity. The evolution of env within the CNS leads to changes that may result in the emergence of novel isolates with neurotoxic and neurovirulent features. However, whether specific factors of HIV-1 evolution lead to the emergence of neurovirulent and neurotropic isolates remains ill-defined. HIV-1 env evolution is an ongoing phenomenon that occurs independently of neurological and neurocognitive disease severity; thus HIV env evolution may play a pivotal and reciprocal role in the etiology of HAND. Despite the use of cART, the reactivation of latent viral reservoirs represents a clinical challenge because of the replenishment of the viral pool that may subsequently lead to persistent infection. Therefore, gaining a more complete understanding of how HIV-1 env evolves over the course of the disease should be considered for the development of future therapies aimed at controlling CNS burden, diminishing persistent viremia, and eradicating viral reservoirs. Here we review the current literature on the role of HIV-1 env evolution in the setting of HAND disease progression and on the impact of cART on the dynamics of viral evolution.

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.

Frequency of coreceptor tropism in PBMC samples from HIV-1 recently infected blood donors by massively parallel sequencing: the REDS II study

Background

The interaction of HIV-1 and target cells involves sequential binding of the viral gp120 Env protein to the CD4 receptor and a chemokine co-receptor (either CCR5 or CXCR4). CCR5 antagonists have proved to be an effective salvage therapy in patients with CCR5 using variants (R5) but not with variants capable of using CXCR4 (×4) phenotype. Thus, it is critically important to determine cellular tropism of a country’s circulating HIV strains to guide a management decision to improve treatment outcome. In this study, we report the prevalence of R5 and ×4 HIV strains in 45 proviral DNA massively parallel sequencing “MPS” data from recently infected Brazilian blood donors.

Methods

The MPS data encompassing the tropism-related V3 loop region of the HIV‐1 env gene was extracted from our recently published HIV-1 genomes sequenced by a paired-end protocol (Illumina). HIV‐1 tropism was inferred using Geno2pheno_[coreceptor] algorithm (3.5 % false-positive rate). V3 net charge and 11/25 rules were also used for coreceptor prediction.

Results

Among the 45 samples for which tropism were determined, 39 were exclusively R5 variants, 5 ×4 variants, and one dual-tropic or mixed (D/M) populations of R5 and ×4 viruses, corresponding to 86.7, 11.1 and 2.2 %, respectively. Thus, the proportion of all blood donors that harbor CXCR4-using virus was 13.3 % including individuals with D/M-tropic viruses.

Conclusions

The presence of CCR5-tropic variants in more than 85 % of our cohort of antiretroviral-naïve blood donors with recent HIV-1 infection indicates a potential benefit of CCR5 antagonists as a therapeutic option in Brazil. Therefore, determination of viral co-receptor tropism is an important diagnostic prerequisite.

Performance of commonly used genotypic assays and comparison with phenotypic assays of HIV-1 coreceptor tropism in acutely HIV-1-infected patients

OBJECTIVES

Although founder viruses in primary HIV-1 infections (PHIs) typically use the CCR5 coreceptor (R5-tropic), 3%-19% of subjects also harbour CXCR4-using viruses (X4-tropic), making tropism determination before CCR5 antagonist usage mandatory. Genotypic methods can be used to accurately determine HIV-1 tropism in chronically infected patients.

METHODS

We compared the results of genotypic methods [geno2pheno, PSSMx4r5 including a novel nucleotide-input version (ntPSSM) and distant segments (ds)Kernel] to predict coreceptor usage in a cohort of 67 PHIs. Specimens with discrepant results were phenotypically tested after cloning the V3 gene region into proviral backbones. Recombinant viruses were used to infect U87 indicator cell lines bearing CD4 and either CCR5 or CXCR4.

RESULTS

Geno2pheno10%, PSSMx4r5 and (ds)Kernel gave identical predictions in 85% of cases. Geno2pheno10% predicted the presence of CXCR4 viruses in 18% of patients. Two patients were predicted to carry X4-tropic viruses by all algorithms and X4-tropic viruses were detected in at least one of the recombinant AD8 or NL4-3 backbone-based assays. Ten samples resulted in discordant predictions with at least one algorithm. Full concordance between tropism prediction by using population sequencing and phenotypic assays was observed only with ntPSSM. Geno2pheno prediction and the phenotypic assay gave the same results in a minority of 'discordant' patients.

CONCLUSIONS

Compared with both PSSMx4r5 versions, (ds)Kernel and our phenotypic assay, geno2pheno10% overestimated the frequency of X4-tropic viruses (18% versus 3%). ntPSSM was able to detect one additional X4 virus compared with (ds)Kernel that was confirmed with the phenotypic assay.

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.

Longitudinal Analysis of Cerebrospinal Fluid and Plasma HIV-1 Envelope Sequences Isolated From a Single Donor with HIV Asymptomatic Neurocognitive Impairment

OBJECTIVE

Combined antiretroviral treatment (cART) has changed the clinical presentation of HIV-associated neurocognitive disorders (HAND) to that of the milder forms of the disease. Asymptomatic neurocognitive impairment (ANI) is now more prevalent and is associated with increased morbidity and mortality risk in HIV-1-infected people. HIV-1 envelope (env) genetic heterogeneity has been detected within the central nervous system (CNS) of individuals with ANI. Changes within env determine co-receptor use, cellular tropism, and neuropathogenesis. We hypothesize that compartmental changes are associated with HIV-1 env C2V4 during ANI and sought to analyze paired HIV-1 env sequences from plasma and cerebrospinal fluid (CSF) of a female subject undergoing long-term cART.

METHODS

Paired plasma and CSF samples were collected at 12-month intervals and HIV-1 env C2V4 was cloned and sequenced.

RESULTS

Phylogenetic analysis of paired samples consistently showed genetic variants unique to the CSF. Phenotypic prediction showed CCR5 (R5) variants for all CSF-derived sequences and showed minor X4 variants (or dual-tropic) in the plasma at later time points. Viral compartmentalization was evident throughout the study, suggesting that the occurrence of distinctive env strains may contribute to the neuropathogenesis of HAND.

CONCLUSIONS

Our study provides new insights about the genetic characteristics within the C2V4 of HIV-1 env that persist after long-term cART and during the course of persistent ANI.

HIV gp120 sequence variability associated with HAND in Hispanic Women

OBJECTIVE

HIV-1 variants with different tropisms are associated with various neuropathologies. This study intends to determine if this correlation is determined by unique viral env sequences. We hypothesize that HIV-1 envelope gene sequence changes are associated with cognition status.

METHODS

Viral RNA was extracted from peripheral blood mononuclear cells (PBMCs) co-cultures derived from HIV-1 infected Hispanic women that had been characterized for HIV associated neurocognitive disorders (HAND).

RESULTS

Analyses of the C2V4 region of HIV gp120 demonstrated that increased sequence diversity correlates with cognition status as sequences derived from subjects with normal cognition exhibited less diversity than sequences derived from subjects with cognitive impairment. In addition, differences in V3 and V4 loop charges were also noted as well as differences in the N-glycosylation of the V4 region.

CONCLUSIONS

Our data suggest that the genetic signature within the C2V4 region may contribute to the pathogenesis of HAND. HIV env sequence characteristics for the isolates grouped in milder forms of HAND can provide insightful information of prognostic value to assess neurocognitive status in HIV⁺ subjects, particularly during the era of highly prevalent milder forms of HAND.

Genotypic prediction of tropism of highly diverse HIV-1 strains from Cameroon

BACKGROUND

The use of CCR5 antagonists involves determination of HIV-1 tropism prior to initiation of treatment. HIV-1 tropism can be assessed either by phenotypic or genotypic methods. Genotypic methods are extensively used for tropism prediction. However, their validation in predicting tropism of viral isolates belonging to group M non-B subtypes remains challenging. In Cameroon, the genetic diversity of HIV-1 strains is the broadest reported worldwide. To facilitate the integration of CCR5 antagonists into clinical practice in this region, there is a need to evaluate the performance of genotypic methods for predicting tropism of highly diverse group M HIV-1 strains.

METHODS

Tropism of diverse HIV-1 strains isolated from PBMCs from Cameroon was determined using the GHOST cell assay. Prediction, based on V3 sequences from matched plasma samples, was determined using bioinformatics algorithms and rules based on position 11/25 and net charge applied independently or combined according to Delobel's and Garrido's rules. Performance of genotypic methods was evaluated by comparing prediction generated with tropism assigned by the phenotypic assay.

RESULTS

Specificity for predicting R5-tropic virus was high, ranging from 83.7% to 97.7% depending on the genotypic methods used. Sensitivity for X4-tropic viruses was fairly low, ranging from 33.3% to 50%. In our study, overall, genotypic methods were less able to accurately predict X4-tropic virus belonging to subtype CRF02_AG. In addition, it was found that of the methods we used the Garrido rule has the highest sensitivity rate of over 50% with a specificity of 93%.

CONCLUSION

Our study demonstrated that overall, genotypic methods were less sensitive for accurate prediction of HIV-1 tropism in settings where diverse HIV-1 strains co-circulate. Our data suggest that further optimization of genotypic methods is needed and that larger studies to determine their utility for tropism prediction of diverse HIV-1 strains may be warranted.

Covariance of charged amino acids at positions 322 and 440 of HIV-1 Env contributes to coreceptor specificity of subtype B viruses, and can be used to improve the performance of V3 sequence-based coreceptor usage prediction algorithms

The ability to determine coreceptor usage of patient-derived human immunodeficiency virus type 1 (HIV-1) strains is clinically important, particularly for the administration of the CCR5 antagonist maraviroc. The envelope glycoprotein (Env) determinants of coreceptor specificity lie primarily within the gp120 V3 loop region, although other Env determinants have been shown to influence gp120-coreceptor interactions. Here, we determined whether conserved amino acid alterations outside the V3 loop that contribute to coreceptor usage exist, and whether these alterations improve the performance of V3 sequence-based coreceptor usage prediction algorithms. We demonstrate a significant covariant association between charged amino acids at position 322 in V3 and position 440 in the C4 Env region that contributes to the specificity of HIV-1 subtype B strains for CCR5 or CXCR4. Specifically, positively charged Lys/Arg at position 322 and negatively charged Asp/Glu at position 440 occurred more frequently in CXCR4-using viruses, whereas negatively charged Asp/Glu at position 322 and positively charged Arg at position 440 occurred more frequently in R5 strains. In the context of CD4-bound gp120, structural models suggest that covariation of amino acids at Env positions 322 and 440 has the potential to alter electrostatic interactions that are formed between gp120 and charged amino acids in the CCR5 N-terminus. We further demonstrate that inclusion of a "440 rule" can improve the sensitivity of several V3 sequence-based genotypic algorithms for predicting coreceptor usage of subtype B HIV-1 strains, without compromising specificity, and significantly improves the AUROC of the geno2pheno algorithm when set to its recommended false positive rate of 5.75%. Together, our results provide further mechanistic insights into the intra-molecular interactions within Env that contribute to coreceptor specificity of subtype B HIV-1 strains, and demonstrate that incorporation of Env determinants outside V3 can improve the reliability of coreceptor usage prediction algorithms.

Determination of viral tropism by genotyping and phenotyping assays in Brazilian HIV-1-infected patients

The clinical application of CCR5 antagonists involves first determining the coreceptor usage by the infecting viral strain. Bioinformatics programs that predict coreceptor usage could provide an alternative method to screen candidates for treatment with CCR5 antagonists, particularly in countries with limited financial resources. Thus, the present study aims to identify the best approach using bioinformatics tools for determining HIV-1 coreceptor usage in clinical practice. Proviral DNA sequences and Trofile results from 99 HIV-1-infected subjects under clinical monitoring were analyzed in this study. Based on the Trofile results, the viral variants present were 81.1% R5, 21.4% R5X4 and 1.8% X4. Determination of tropism using a Geno2pheno_[coreceptor] analysis with a false positive rate of 10% gave the most suitable performance in this sampling: the R5 and X4 strains were found at frequencies of 78.5% and 28.4%, respectively, and there was 78.6% concordance between the phenotypic and genotypic results. Further studies are needed to clarify how genetic diversity amongst virus strains affects bioinformatics-driven approaches for determining tropism. Although this strategy could be useful for screening patients in developing countries, some limitations remain that restrict the wider application of coreceptor usage tests in clinical practice.

HIV-1 coreceptor usage in paired plasma RNA and proviral DNA from patients with acute and chronic infection never treated with antiretroviral therapy

Although an independent evolution of viral quasispecies in different body sites might determine a differential compartmentalization of viral variants, the aim of this paper was to establish whether sequences from peripheral blood mononuclear cells (PBMCs) and plasma provide different or complementary information on HIV tropism in patients with acute or chronic infection. Tropism was predicted using genotypic testing combined with geno2pheno (coreceptor) analysis at a 10% false positive rate in paired RNA and DNA samples from 75 antiretroviral-naïve patients (divided on the basis of avidity index into patients with a recent or long-lasting infection). A high prevalence of R5 HIV strains (97%) was observed in both compartments (plasma and PBMCs) in patients infected recently. By contrast, patients with a long-lasting infection showed a quite different situation in the two compartments, revealing more (46%) X4/DM in PBMCs than patients infected recently (3%) (P = 0.008). As- a knowledge of viral strains in different biological compartments might be crucial to establish a therapeutic protocol, it could be extremely useful to detect not only viral strains in plasma, but also viruses hidden or archived in different cell compartments to better understand disease evolution and treatment efficacy in patients infected with HIV.

Genotypic tropism prediction from paired cell and plasma using single and replicate sequences

HIV-1 tropism determination is necessary prior to CCR5 antagonist use as antiretroviral therapy. Genotypic prediction of coreceptor use is a practical alternative to phenotypic tests. Cell DNA and plasma RNA-based prediction has shown discordance in many studies. We evaluate paired cell and plasma either as single or replicate V3 sequences to assess prediction comparability. The HIV-1 partial env region was sequenced and tropism was predicted using geno2pheno and position-specific scoring matrices (PSSM). Nucleotide ambiguities at V3 were quantified and genetic distance (Protdist) was determined using BioEdit. Wilcoxon signed-rank test, t tests, and Spearman correlation were performed with Prism GraphPad5.0. Results are expressed as medians, with a level of significance of p<0.05, two tailed. Single (n=28) or replicate (n=26) paired cell/plasma sequences were obtained from 54 patients. Although the clonalfalse-positive rate (FPR) value from both compartments strongly correlated (r=0.86 p<0.0001), discordance in tropism prediction was observed in both singles and replicates using geno2pheno or PSSM. Applying clonalFPR(10%) 46% (25/54) were X4 tropic, with a plasma/cell discordance of 11% in singles and 23% in replicates. Genetic distance (p<0.0001) and clonalFPR value dispersion (p=0.003) were significantly higher among replicate sequences from cells. Discordance of viral tropism prediction is not uncommon and the use of replicates does not decrease its occurrence, but improves X4 sensitivity. Sequences from provirus had greater genetic distance and dispersion of clonalFPR values. This may suggest that DNA replicate assays may better represent the diversity of HIV-1 variants, but the clinical significance of these findings needs further evaluation.

Sensitive deep-sequencing-based HIV-1 genotyping assay to simultaneously determine susceptibility to protease, reverse transcriptase, integrase, and maturation inhibitors, as well as HIV-1 coreceptor tropism

With 29 individual antiretroviral drugs available from six classes that are approved for the treatment of HIV-1 infection, a combination of different phenotypic and genotypic tests is currently needed to monitor HIV-infected individuals. In this study, we developed a novel HIV-1 genotypic assay based on deep sequencing (DeepGen HIV) to simultaneously assess HIV-1 susceptibilities to all drugs targeting the three viral enzymes and to predict HIV-1 coreceptor tropism. Patient-derived gag-p2/NCp7/p1/p6/pol-PR/RT/IN- and env-C2V3 PCR products were sequenced using the Ion Torrent Personal Genome Machine. Reads spanning the 3' end of the Gag, protease (PR), reverse transcriptase (RT), integrase (IN), and V3 regions were extracted, truncated, translated, and assembled for genotype and HIV-1 coreceptor tropism determination. DeepGen HIV consistently detected both minority drug-resistant viruses and non-R5 HIV-1 variants from clinical specimens with viral loads of ≥1,000 copies/ml and from B and non-B subtypes. Additional mutations associated with resistance to PR, RT, and IN inhibitors, previously undetected by standard (Sanger) population sequencing, were reliably identified at frequencies as low as 1%. DeepGen HIV results correlated with phenotypic (original Trofile, 92%; enhanced-sensitivity Trofile assay [ESTA], 80%; TROCAI, 81%; and VeriTrop, 80%) and genotypic (population sequencing/Geno2Pheno with a 10% false-positive rate [FPR], 84%) HIV-1 tropism test results. DeepGen HIV (83%) and Trofile (85%) showed similar concordances with the clinical response following an 8-day course of maraviroc monotherapy (MCT). In summary, this novel all-inclusive HIV-1 genotypic and coreceptor tropism assay, based on deep sequencing of the PR, RT, IN, and V3 regions, permits simultaneous multiplex detection of low-level drug-resistant and/or non-R5 viruses in up to 96 clinical samples. This comprehensive test, the first of its class, will be instrumental in the development of new antiretroviral drugs and, more importantly, will aid in the treatment and management of HIV-infected individuals.

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.

Genotypic testing on HIV-1 DNA as a tool to assess HIV-1 co-receptor usage in clinical practice: results from the DIVA study group

PURPOSE

We have developed a sequencing assay for determining the usage of the genotypic HIV-1 co-receptor using peripheral blood mononuclear cell (PBMC) DNA in virologically suppressed HIV-1 infected patients. Our specific aims were to (1) evaluate the efficiency of V3 sequences in B versus non-B subtypes, (2) compare the efficiency of V3 sequences and tropism prediction using whole blood and PBMCs for DNA extraction, (3) compare the efficiency of V3 sequences and tropism prediction using a single versus a triplicate round of amplification.

RESULTS

The overall rate of successful V3 sequences ranged from 100 % in samples with >3,000 copies HIV-1 DNA/10(6) PBMCs to 60 % in samples with <100 copies total HIV-1 DNA /10(6) PBMCs. Analysis of 143 paired PBMCs and whole-blood samples showed successful V3 sequences rates of 77.6 % for PBMCs and 83.9 % for whole blood. These rates are in agreement with the tropism prediction obtained using the geno2pheno co-receptor algorithm, namely, 92.1 % with a false-positive rate (FPR) of 10 or 20 % and of 96.5 % with an FPR of 5.75 %. The agreement between tropism prediction values using single versus triplicate amplification was 98.2 % (56/57) of patients using an FPR of 20 % and 92.9 % (53/57) using an FPR of 10 or 5.75 %. For 63.0 % (36/57) of patients, the FPR obtained via the single amplification procedure was superimposable to all three FPRs obtained by triplicate amplification.

CONCLUSIONS

Our results show the feasibility and consistency of genotypic testing on HIV-1 DNA tropism, supporting its possible use for selecting patients with suppressed plasma HIV-1 RNA as candidates for CCR5-antagonist treatment. The high agreement between tropism prediction by single and triple amplification does not support the use of triplicate amplification in clinical practice.

Determination of the high prevalence of Dual/Mixed- or X4-tropism among HIV type 1 CRF01_AE in Hong Kong by genotyping and phenotyping methods

In Hong Kong, the CCR5 antagonist has recently been introduced into salvage therapy for multiclass drug-resistant HIV-1-infected patients. Coreceptor usage must be determined prior to the usage of the CCR5 antagonist, which does not inhibit X4-tropic viruses. This study aimed to determine the tropism prevalence for HIV-1 subtypes B and CRF01_AE in Hong Kong. In addition, a modified promoter-PCR phenotypic assay was used to validate the genotypic tropism prediction on CRF01_AE. One hundred and five subtype B and 98 CRF01_AE antiretroviral-naive patients were recruited for this study. The viral env V3 region isolated from the patients was sequenced and analyzed by Geno2pheno (FPR=5.75% or 10%, Clonal or Clinical), position-specific scoring matrix (WebPSSM, x4r5 subtype B matrix), and the combination of 11/25 and net charge rules. Fifteen concordant and 22 discordant tropism genotyped CRF01_AE samples were further phenotyped by either enhanced sensitivity Trofile assay or an optimized promoter-PCR phenotypic assay. The prevalence of Dual/Mixed- or X4-tropic virus in antiretroviral-naive subtype CRF01_AE was 39.1%, which was significantly higher than subtype B (p<0.05), regardless of the choices of genotypic algorithms. Our phenotypic data proposed that a better genotypic tropism prediction for HIV-1 CRF01_AE would be using both Geno2pheno (FPR=10%, Clonal) and WebPSSM (x4r5 subtype B matrix) algorithms in combination. The sensitivity and specificity for this combination were 88.9% and 89.3%, respectively. The comparatively high prevalence of Dual/Mixed- or X4-tropic virus in CRF01_AE demonstrated the need for special attention to future treatment strategies.

HIV-1 coreceptor switch during 2 years of structured treatment interruptions

The purpose of this investigation was to determine the impact on human immunodeficiency virus (HIV) tropism of uncontrolled virus exposure during 2 years of intermittent highly active antiretroviral therapy (HAART). The Istituto Superiore di Sanità-Pulsed Antiretroviral Therapy (ISS-PART) randomized study compared the outcome of 2 years of structured treatment interruptions (STIs) versus standard continuous treatment in first-line HAART responder subjects. The STI schedule consisted of five STIs of 1, 1, 2, 2, and 3 months, respectively, separated by four periods of 3-month therapy. In the present study, coreceptor tropism was assessed in 12 patients of the STI arm at different time points over a period of 2 years. Tropism was determined on DNA and RNA by V3 loop region sequencing. The Geno2pheno algorithm (false-positive rate, FPR: 20%) was used for data interpretation. At baseline, 9/12 subjects (75.0%) had CCR5-tropic viruses in their HIV. Three had a CXCR4-tropic virus. Ten patients maintained the same coreceptor in DNA after 2 years, whereas in two patients, a shift occurred (one R5-X4, one X4-R5). In a patient with an R5 virus at baseline, a transient change to X4 tropism was seen in the rebounding virus during STI. Changes in tropism were not associated with the amplitude and duration of virus exposure during STIs, residual viremia at baseline, or the development of resistance mutations in the RT region. Our preliminary results suggest that viral replication, observed after short periods of treatment interruption, is not enough to drive the evolution of HIV tropism.

RECall for automated genotypic tropism testing

Standardization of sequence chromatogram analysis is required for consistent genotypic tropism determination across laboratories. A freely available, fast, and automated chromatogram analysis tool (RECall) provided tropism interpretations equivalent to those of manual sequence editing of 521 V3 loop HIV-1 population sequences, suggesting that RECall can be useful in standardizing genotypic tropism testing across laboratories.

HIV-1 tropism determination using a phenotypic Env recombinant viral assay highlights overestimation of CXCR4-usage by genotypic prediction algorithms for CRF01_AE and CRF02_AG [corrected]

BACKGROUND

Human Immunodeficiency virus type-1 (HIV) entry into target cells involves binding of the viral envelope (Env) to CD4 and a coreceptor, mainly CCR5 or CXCR4. The only currently licensed HIV entry inhibitor, maraviroc, targets CCR5, and the presence of CXCX4-using strains must be excluded prior to treatment. Co-receptor usage can be assessed by phenotypic assays or through genotypic prediction. Here we compared the performance of a phenotypic Env-Recombinant Viral Assay (RVA) to the two most widely used genotypic prediction algorithms, Geno2Pheno[coreceptor] and webPSSM.

METHODS

Co-receptor tropism of samples from 73 subtype B and 219 non-B infections was measured phenotypically using a luciferase-tagged, NL4-3-based, RVA targeting Env. In parallel, tropism was inferred genotypically from the corresponding V3-loop sequences using Geno2Pheno[coreceptor] (5-20% FPR) and webPSSM-R5X4. For discordant samples, phenotypic outcome was retested using co-receptor antagonists or the validated Trofile® Enhanced-Sensitivity-Tropism-Assay.

RESULTS

The lower detection limit of the RVA was 2.5% and 5% for X4 and R5 minority variants respectively. A phenotype/genotype result was obtained for 210 samples. Overall, concordance of phenotypic results with Geno2Pheno[coreceptor] was 85.2% and concordance with webPSSM was 79.5%. For subtype B, concordance with Geno2pheno[coreceptor] was 94.4% and concordance with webPSSM was 79.6%. High concordance of genotypic tools with phenotypic outcome was seen for subtype C (90% for both tools). Main discordances involved CRF01_AE and CRF02_AG for both algorithms (CRF01_AE: 35.9% discordances with Geno2Pheno[coreceptor] and 28.2% with webPSSM; CRF02_AG: 20.7% for both algorithms). Genotypic prediction overestimated CXCR4-usage for both CRFs. For webPSSM, 40% discordance was observed for subtype A.

CONCLUSIONS

Phenotypic assays remain the most accurate for most non-B subtypes and new subtype-specific rules should be developed for non-B subtypes, as research studies more and more draw conclusions from genotypically-inferred tropism, and to avoid unnecessarily precluding patients with limited treatment options from receiving maraviroc or other entry inhibitors.

Variability of HIV-1 genomes among children and adolescents from São Paulo, Brazil

BACKGROUND

Genetic variability is a major feature of the human immunodeficiency virus type 1 (HIV-1) and considered the key factor to frustrating efforts to halt the virus epidemic. In this study, we aimed to investigate the genetic variability of HIV-1 strains among children and adolescents born from 1992 to 2009 in the state of Sao Paulo, Brazil.

METHODOLOGY

Plasma and peripheral blood mononuclear cells (PBMC) were collected from 51 HIV-1-positive children and adolescents on ART followed between September 1992 and July 2009. After extraction, the genetic materials were used in a polymerase chain reaction (PCR) to amplify the viral near full length genomes (NFLGs) from 5 overlapped fragments. NFLGs and partial amplicons were directly sequenced and data were phylogenetically inferred.

RESULTS

Of the 51 samples studied, the NFLGs and partial fragments of HIV-1 from 42 PBMCs and 25 plasma were successfully subtyped. Results based on proviral DNA revealed that 22 (52.4%) patients were infected with subtype B, 16 (38.1%) were infected with BF1 mosaic variants and 4 (9.5%) were infected with sub-subtype F1. All the BF1 recombinants were unique and distinct from any previously identified unique or circulating recombinant forms in South America. Evidence of dual infections was detected in 3 patients coinfected with the same or distinct HIV-1 subtypes. Ten of the 31 (32.2%) and 12 of the 21 (57.1%) subjects with recovered proviral and plasma, respectively, protease sequences were infected with major mutants resistant to protease inhibitors. The V3 sequences of 14 patients with available sequences from PBMC/or plasma were predicted to be R5-tropic virus except for two patients who harbored an X4 strain.

CONCLUSIONS

The high proportion of HIV-1 BF1 recombinant, coinfection rate and vertical transmission in Brazil merits urgent attention and effective measures to reduce the transmission of HIV among spouses and sex partners.

Prevalence of R5 strains in multi-treated HIV subjects and impact of new regimens including maraviroc in a selected group of patients with CCR5-tropic HIV-1 infection

OBJECTIVES

Maraviroc currently represents an important antiretroviral drug for multi-experienced and viremic HIV patients. This study focused on two main points: (1) determining the prevalence of R5 and X4 HIV strains in antiretroviral-experienced patients using two main tests currently in use to determine viral tropism, and (2) the follow-up to 3 years of a limited number of patients who started a new antiretroviral protocol including maraviroc.

METHODS

A group of 56 HIV patients, previously multi-treated, were first analyzed by genotyping assay and Trofile™ to establish their eligibility for maraviroc treatment. In addition, 25 subjects selected to follow a new therapeutic protocol including a CCR5 antagonist were monitored by HIV RNA viral load and CD4+ cell count.

RESULTS

The determination of viral tropism showed a large percentage of patients with an R5 profile (72% by genotyping assay and 74% by Trofile). The follow-up of most (21 out 25) patients who started the new antiretroviral protocol showed an undetectable viral load throughout the observation period, accompanied by a major improvement in CD4 cell count (cells/mm(3)) (baseline: median CD4 cell count 365, interquartile range (IQR) 204-511; 12 months: median value 501, IQR 349-677, p=0.042; 24 months: median value 503, IQR 386-678, p=0.026; 36 months: median value 601, IQR 517-717, p=0.001). Among the four non-responder subjects, two showed a lack of drug compliance and two switched from R5 to X4.

CONCLUSION

Although our patient cohort was small, the results showed a high prevalence of R5 viral strains in multi-experienced patients. As well as showing the advantages of genotyping, which can be performed in plasma samples with low viral load replication, the follow-up of HIV patients selected for an alternative drug protocol, including a CCR5 antagonist, showed a persistent undetectable viral replication and a good recovery of CD4 cell count in most treated HIV patients.

Sensitive cell-based assay for determination of human immunodeficiency virus type 1 coreceptor tropism

CCR5 antagonists are a powerful new class of antiretroviral drugs that require a companion assay to evaluate the presence of CXCR4-tropic (non-R5) viruses prior to use in human immunodeficiency virus (HIV)-infected individuals. In this study, we have developed, characterized, verified, and prevalidated a novel phenotypic test to determine HIV-1 coreceptor tropism (VERITROP) based on a sensitive cell-to-cell fusion assay. A proprietary vector was constructed containing a near-full-length HIV-1 genome with the yeast uracil biosynthesis (URA3) gene replacing the HIV-1 env coding sequence. Patient-derived HIV-1 PCR products were introduced by homologous recombination using an innovative yeast-based cloning strategy. The env-expressing vectors were then used in a cell-to-cell fusion assay to determine the presence of R5 and/or non-R5 HIV-1 variants within the viral population. Results were compared with (i) the original version of Trofile (Monogram Biosciences, San Francisco, CA), (ii) population sequencing, and (iii) 454 pyrosequencing, with the genotypic data analyzed using several bioinformatics tools, i.e., the 11/24/25 rule, Geno2Pheno (2% to 5.75%, 3.5%, or 10% false-positive rate [FPR]), and webPSSM. VERITROP consistently detected minority non-R5 variants from clinical specimens, with an analytical sensitivity of 0.3%, with viral loads of ≥1,000 copies/ml, and from B and non-B subtypes. In a pilot study, a 73.7% (56/76) concordance was observed with the original Trofile assay, with 19 of the 20 discordant results corresponding to non-R5 variants detected using VERITROP and not by the original Trofile assay. The degree of concordance of VERITROP and Trofile with population and deep sequencing results depended on the algorithm used to determine HIV-1 coreceptor tropism. Overall, VERITROP showed better concordance with deep sequencing/Geno2Pheno at a 0.3% detection threshold (67%), whereas Trofile matched better with population sequencing (79%). However, 454 sequencing using Geno2Pheno at a 10% FPR and 0.3% threshold and VERITROP more accurately predicted the success of a maraviroc-based regimen. In conclusion, VERITROP may promote the development of new HIV coreceptor antagonists and aid in the treatment and management of HIV-infected individuals prior to and/or during treatment with this class of drugs.

CoRSeqV3-C: a novel HIV-1 subtype C specific V3 sequence based coreceptor usage prediction algorithm

Background

The majority of HIV-1 subjects worldwide are infected with HIV-1 subtype C (C-HIV). Although C-HIV predominates in developing regions of the world such as Southern Africa and Central Asia, C-HIV is also spreading rapidly in countries with more developed economies and health care systems, whose populations are more likely to have access to wider treatment options, including the CCR5 antagonist maraviroc (MVC). The ability to reliably determine C-HIV coreceptor usage is therefore becoming increasingly more important. In silico V3 sequence based coreceptor usage prediction algorithms are a relatively rapid and cost effective method for determining HIV-1 coreceptor specificity. In this study, we elucidated the V3 sequence determinants of C-HIV coreceptor usage, and used this knowledge to develop and validate a novel, user friendly, and highly sensitive C-HIV specific coreceptor usage prediction algorithm.

Results

We characterized every phenotypically-verified C-HIV gp120 V3 sequence available in the Los Alamos HIV Database. Sequence analyses revealed that compared to R5 C-HIV V3 sequences, CXCR4-using C-HIV V3 sequences have significantly greater amino acid variability, increased net charge, increased amino acid length, increased frequency of insertions and substitutions within the GPGQ crown motif, and reduced frequency of glycosylation sites. Based on these findings, we developed a novel C-HIV specific coreceptor usage prediction algorithm (CoRSeq_V3-C), which we show has superior sensitivity for determining CXCR4 usage by C-HIV strains compared to all other available algorithms and prediction rules, including Geno2pheno_[coreceptor] and WebPSSM_SINSI-C, which has been designed specifically for C-HIV.

Conclusions

CoRSeq_V3-C is now openly available for public use at http://www.burnet.edu.au/coreceptor. Our results show that CoRSeq_V3-C is the most sensitive V3 sequence based algorithm presently available for predicting CXCR4 usage of C-HIV strains, without compromising specificity. CoRSeq_V3-C may be potentially useful for assisting clinicians to decide the best treatment options for patients with C-HIV infection, and will be helpful for basic studies of C-HIV pathogenesis.

High prevalence of CXCR4 usage among treatment-naive CRF01_AE and CRF51_01B-infected HIV-1 subjects in Singapore

Background

Recent studies suggest HIV-1 inter-subtype differences in co-receptor usage. We examined the correlation between HIV-1 subtype and co-receptor usage among treatment-naïve HIV-1 subjects in Singapore. Additionally, we investigated whether the subtype co-receptor association was influenced by stage of infection.

Methods

V3 sequences of HIV-1 envelope protein gp120 were obtained from 110 HIV treatment-naïve patients and genotypic co-receptor tropism determination was performed using Geno2pheno. Two false-positive rate (FPR) cut-offs, 10% and 5.75% were selected for tropism testing.

Results

Subtype assignment of viral strains from 110 HIV-infected individuals based on partial sequencing of HIV-1 pol, gp120 and gp41 were as follows: 27 subtype B, 64 CRF01_AE, 10 CRF51_01B, and 9 other subtypes. At FPR=10%, 10 (100%) CRF51_01B-infected subjects and 26 (40.6%) CRF01_AE-infected subjects had CXCR4-using virus, compared to 7 (25.9%) subtype B subjects and 1 (11.1%) CRF33_01B-infected subject (P < 0.001). At FPR=5.75%, 10 (100%) CRF51_01B-infected subjects and 20 (31.3%) CRF01_AE-infected subjects had CXCR4-using virus, compared to 4 (14.8%) subtype B and 1 (11.1%) CRF33_01B-infected subjects (P < 0.001). Among those with evidence of seroconversion within 2 years prior to study enrolment, 100% of CRF51_01B-infected subjects had CXCR4-using virus, independent of Geno2pheno FPR.

Conclusion

CRF51_01B and CRF01_AE-infected individuals have higher prevalence of CXCR4-usage compared to subtype B infected individuals. Further studies examining these differences could help optimise the use of CCR5-antagonist in populations with these subtypes, and increase our understanding of HIV-1 biology.

The genotypic false positive rate determined by V3 population sequencing can predict the burden of HIV-1 CXCR4-using species detected by pyrosequencing

Objective

The false-positive rate (FPR) is a percentage-score provided by Geno2Pheno-algorithm indicating the likelihood that a V3-sequence is falsely predicted as CXCR4-using. We evaluated the correlation between FPR obtained by V3 population-sequencing and the burden of CXCR4-using variants detected by V3 ultra-deep sequencing (UDPS) and Enhanced-Sensitivity Trofile assay (ESTA).

Methods

54 HIV-1 B-subtype infected-patients (all maraviroc-naïve), with viremia >10,000copies/ml, were analyzed. HIV-tropism was assessed by V3 population-sequencing, UDPS (considering variants with >0.5% prevalence), and ESTA.

Results

By UDPS, CCR5-using variants were detected in 53/54 patients, irrespective of FPR values, and their intra-patient prevalence progressively increased by increasing the FPR obtained by V3 population-sequencing (rho = 0.75, p = 5.0e-8). Conversely, the intra-patient prevalence of CXCR4-using variants in the 54 patients analyzed progressively decreased by increasing the FPR (rho = −0.61; p = 9.3e-6). Indeed, no CXCR4-using variants were detected in 13/13 patients with FPR>60. They were present in 7/18 (38.8%) patients with FPR 20–60 (intra-patient prevalence range: 2.1%–18.4%), in 5/7 (71.4%) with FPR 10–20, in 4/6 (66.7%) with FPR 5–10, and in 10/10(100%) with FPR<5 (intra-patient prevalence range: 12.1%–98.1%).

Conclusions

FPR by V3 population-sequencing can predict the burden of CXCR4-using variants. This information can be used to optimize the management of tropism determination in clinical practice. Due to its low cost and short turnaround time, V3 population-sequencing may represent the most feasible test for HIV-1 tropism determination. More sensitive methodologies (as UDPS) might be useful when V3 population-sequencing provides a FPR >20 (particularly in the range 20–60), allowing a more careful identification of patients harboring CXCR4-using variants.

Genotypic prediction of HIV-1 CRF01-AE tropism

HIV-1 subtype CRF01-AE predominates in south Asia and has spread throughout the world. The virus tropism must be determined before using CCR5 antagonists. Genotypic methods could be used, but the prediction algorithms may be inaccurate for non-B subtypes like CRF01-AE and the correlation with the phenotypic approach has not been assessed. We analyzed 61 CRF01-AE V3 clonal sequences of known phenotype from the GenBank database. The sensitivity of the Geno2pheno10 genotypic algorithm was 91%, but its specificity was poor (54%). In contrast, the combined 11/25 and net charge rule was highly specific (98%) but rather insensitive (64%). We thus identified subtype CRF01-AE determinants in the V3 region that are associated with CXCR4 use and developed a new simple rule for optimizing the genotypic prediction of CRF01-AE tropism. The concordance between the predicted CRF01-AE genotype and the phenotype was 95% for the clonal data set. We then validated this algorithm by analyzing the data from 44 patients infected with subtype CRF01-AE, whose tropism was determined using a recombinant phenotypic entry assay and V3-loop bulk sequencing. The CRF01-AE genotypic tool was 70% sensitive and 96% specific for predicting CXCR4 use, and the concordance between genotype and phenotype was 84%, approaching the concordance obtained for predicting the tropism of HIV-1 subtype B. Genotypic predictions that use a subtype CRF01-AE-specific algorithm appear to be preferable for characterizing coreceptor usage both in pathophysiological studies and for ensuring the appropriate use of CCR5 antagonists.

HIV-1 tropism and liver fibrosis in HIV-HCV co-infected patients

Background and Aims

Hepatic stellate cells, the major producers of extracellular matrix in the liver, and hepatocytes bear CXCR4 and CCR5, the two main co-receptors for entry of the human immunodeficiency virus (HIV). In vitro studies suggest that HIV-envelope proteins can modulate the replication of hepatitis C virus (HCV) and fibrogenesis. We investigated the influence of HIV tropism on liver fibrosis and the concentration of HCV RNA in HIV–HCV co-infected patients.

Methods

We used a phenotypic assay to assess HIV tropism in 172 HCV–HIV co-infected patients: one group (75 patients) had mild fibrosis (score ≤F2) and the other (97 patients) had severe fibrosis (score >F2). We also assessed the relationship between HIV tropism and HCV RNA concentration in all these patients. We also followed 34 of these patients for 3 years to determine the evolution of HIV tropism and liver fibrosis, estimated by liver stiffness.

Results

Initially, most patients (91.8%) received a potent antiretroviral therapy. CXCR4-using viruses were found in 29% of patients. The only factor associated with a CXCR4-using virus infection in multivariate analysis was the nadir of CD4 cells: <200/mm³ (OR: 3.94, 95%CI: 1.39–11.14). The median HCV RNA concentrations in patients infected with R5 viruses, those with dual-mixed viruses and those with X4 viruses, were all similar. The prevalence of CXCR4-using viruses in patients with mild fibrosis (≤F2) (31%) and those with severe fibrosis (F3–F4) (28%, p = 0.6) was similar. Longitudinal analyses showed that the presence of CXCR4-using viruses did not increase the likelihood of fibrosis progression, evaluated by measuring liver stiffness.

Conclusions

The presence of CXCR4-using viruses in patients receiving a potent antiretroviral therapy does not influence HCV RNA concentration or liver fibrosis.

Use of four next-generation sequencing platforms to determine HIV-1 coreceptor tropism

HIV-1 coreceptor tropism assays are required to rule out the presence of CXCR4-tropic (non-R5) viruses prior treatment with CCR5 antagonists. Phenotypic (e.g., Trofile™, Monogram Biosciences) and genotypic (e.g., population sequencing linked to bioinformatic algorithms) assays are the most widely used. Although several next-generation sequencing (NGS) platforms are available, to date all published deep sequencing HIV-1 tropism studies have used the 454™ Life Sciences/Roche platform. In this study, HIV-1 co-receptor usage was predicted for twelve patients scheduled to start a maraviroc-based antiretroviral regimen. The V3 region of the HIV-1 env gene was sequenced using four NGS platforms: 454™, PacBio® RS (Pacific Biosciences), Illumina®, and Ion Torrent™ (Life Technologies). Cross-platform variation was evaluated, including number of reads, read length and error rates. HIV-1 tropism was inferred using Geno2Pheno, Web PSSM, and the 11/24/25 rule and compared with Trofile™ and virologic response to antiretroviral therapy. Error rates related to insertions/deletions (indels) and nucleotide substitutions introduced by the four NGS platforms were low compared to the actual HIV-1 sequence variation. Each platform detected all major virus variants within the HIV-1 population with similar frequencies. Identification of non-R5 viruses was comparable among the four platforms, with minor differences attributable to the algorithms used to infer HIV-1 tropism. All NGS platforms showed similar concordance with virologic response to the maraviroc-based regimen (75% to 80% range depending on the algorithm used), compared to Trofile (80%) and population sequencing (70%). In conclusion, all four NGS platforms were able to detect minority non-R5 variants at comparable levels suggesting that any NGS-based method can be used to predict HIV-1 coreceptor usage.

Primary CXCR4 co-receptor use in acute HIV infection leads to rapid disease progression in the AE subtype

This is a comparative study of HIV co-receptor usage in the early stages of HIV infection between two distinct patient groups, one with a low CD4 count (group 1), and the other with a high CD4 count (group 2). Group 1 progressed to a CD4 count below 200 cells/μL within 2 y, while group 2 had a CD4 count above 500 cells/μL within 2 y. Viral RNA was extracted from the plasma of these patients, and the C2-V5 region of the HIV-1 env genes were cloned and sequenced. The co-receptor usage was predicated based on V3 loop amino acid sequences using Geno2pheno and PSSM programs. Our results indicate that in acute HIV infection of rapid progressors (low CD4 count; group 1), the primary co-receptor usage is CXCR4, while in the high CD4 count group (group 2), the co-receptor usage is predominantly CCR5. One-year follow-up data from these patients showed no obvious change in HIV co-receptor usage in either group. Sequence analysis of patients from both study groups showed prevalence of the AE subtype, and therefore we can speculate that the CXCR4 co-receptor may be the primary HIV-1 co-receptor used in the HIV-1 AE subtype, and may be responsible for rapid HIV-1 disease progression in the MSM cohort.

Selected amino acid changes in HIV-1 subtype-C gp41 are associated with specific gp120(V3) signatures in the regulation of co-receptor usage

The majority of studies have characterized the tropism of HIV-1 subtype-B isolates, but little is known about the determinants of tropism in other subtypes. So, the goal of the present study was to genetically characterize the envelope of viral proteins in terms of co-receptor usage by analyzing 356 full-length env sequences derived from HIV-1 subtype-C infected individuals. The co-receptor usage of V3 sequences was inferred by using the Geno2Pheno and PSSM algorithms, and also analyzed to the "11/25 rule". All reported env sequences were also analyzed with regard to N-linked glycosylation sites, net charge and hydrophilicity, as well as the binomial correlation phi coefficient to assess covariation among gp120(V3) and gp41 signatures and the average linkage hierarchical agglomerative clustering were also performed. Among env sequences present in Los Alamos Database, 255 and 101 sequences predicted as CCR5 and CXCR4 were selected, respectively. The classical V3 signatures at positions 11 and 25, and other specific V3 and gp41 amino acid changes were found statistically associated with different co-receptor usage. Furthermore, several statistically significant associations between V3 and gp41 signatures were also observed. The dendrogram topology showed a cluster associated with CCR5-usage composed by five gp41 mutated positions, A22V, R133M, E136G, N140L, and N166Q that clustered with T2V(V3) and G24T(V3) (bootstrap=1). Conversely, a heterogeneous cluster with CXCR4-usage, involving S11GR(V3), 13-14insIG/LG(V3), P16RQ(V3), Q18KR(V3), F20ILV(V3), D25KRQ(V3), Q32KR(V3) along with A30T(gp41), S107N(gp41), D148E(gp41), A189S(gp41) was identified (bootstrap=0.86). Our results show that as observed for HIV-1 subtype-B, also in subtype-C specific and different gp41 and gp120V3 amino acid changes are associated individually or together with CXCR4 and/or CCR5 usage. These findings strengthen previous observations that determinants of tropism may also reside in the gp41 protein.

Discordance in HIV-1 co-receptor use prediction by different genotypic algorithms and phenotype assay: intermediate profile in relation to concordant predictions

Concordant and discordant genotypic predictions of HIV-1 co-receptor tropism were analyzed. V3 region was sequenced from plasma samples of patients screened for R5 tropism by the Trofile® assay, before CCR5 antagonist prescription. Ten tools including geno2pheno, PSSM, an "11/25" and "net charge" rule, and other published algorithms were used. Patients were grouped according to concordance or discordance between tools and Trofile® result. Trofile® tropism reports from 50 patient samples were R5 in 38 and Dual/Mixed (DM) in 12. Prediction with the genotypic tools were concordant for 23 R5 samples, and discordant for the 15 other ones. From Trofile® DM strains were concordant in 6 and discordant in 6. V3 sequences were not clearly distinct between R5 and DM strains, except a greater diversity in the later. Discordances were found with any tool or combination of them, so that no one can be proposed as better than the others. Predictive values of each algorithm were similar and rather good (efficacy ranged from 74% to 84%), but the rate of non-confirmed prediction is greater when compelling the results of all tools with each individual sample. The mean of quantitative values obtained with one tool when another tool give the opposite prediction were different from those obtained when all tools agree with that prediction. The two discordant groups were often not distinguishable from each other. These results suggest that viruses giving discordant prediction with bioinformatic tools could be functionally distinct and/or in a different evolutionary state compared to those with concordant prediction.