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

Insights to HIV-1 coreceptor usage by estimating HLA adaptation with Bayesian generalized linear mixed models

The mechanisms triggering the human immunodeficiency virus type I (HIV-1) to switch the coreceptor usage from CCR5 to CXCR4 during the course of infection are not entirely understood. While low CD4+ T cell counts are associated with CXCR4 usage, a predominance of CXCR4 usage with still high CD4+ T cell counts remains puzzling. Here, we explore the hypothesis that viral adaptation to the human leukocyte antigen (HLA) complex, especially to the HLA class II alleles, contributes to the coreceptor switch. To this end, we sequence the viral gag and env protein with corresponding HLA class I and II alleles of a new cohort of 312 treatment-naive, subtype C, chronically-infected HIV-1 patients from South Africa. To estimate HLA adaptation, we develop a novel computational approach using Bayesian generalized linear mixed models (GLMMs). Our model allows to consider the entire HLA repertoire without restricting the model to pre-learned HLA-polymorphisms. In addition, we correct for phylogenetic relatedness of the viruses within the model itself to account for founder effects. Using our model, we observe that CXCR4-using variants are more adapted than CCR5-using variants (p-value = 1.34e-2). Additionally, adapted CCR5-using variants have a significantly lower predicted false positive rate (FPR) by the geno2pheno[coreceptor] tool compared to the non-adapted CCR5-using variants (p-value = 2.21e-2), where a low FPR is associated with CXCR4 usage. Consequently, estimating HLA adaptation can be an asset in predicting not only coreceptor usage, but also an approaching coreceptor switch in CCR5-using variants. We propose the usage of Bayesian GLMMs for modeling virus-host adaptation in general.

Phenotypic and Genotypic Co-receptor Tropism Testing in HIV-1 Epidemic Region of Tanzania Where Multiple Non-B Subtypes Co-circulate

HIV human immunodeficiency virus type I (HIV-1) entry inhibitor potency is dependent on viral co-receptor tropisms and thereby tropism determination is clinically important. However, phenotypic tropisms of HIV-1 non-B subtypes have been poorly investigated and the genotypic prediction algorithms remain insufficiently validated. To clarify this issue, we recruited 52 treatment-naïve, HIV-1-infected patients in Tanzania, where multiple HIV-1 non-B subtypes co-circulate. Sequence analysis of 93 infectious envelope clones isolated from their plasma viral RNA revealed the co-circulation of subtypes A1, C, D, and inter-subtype recombinant forms (isRFs). Phenotypic tropism assays revealed that lentivirus reporters pseudotyped with 75 (80.6%) and 5 (5.4%) envelope clones could establish infection toward U87.CD4 cells expressing CCR5 (R5) and CXCR4 (X4), respectively; whereas the remaining 13 (14%) clones could infect both cells. Genotypic analyses by widely used algorithms including V3 net charge, Geno2pheno, WebPSSM, and PhenoSeq showed that almost all phenotypic X4-tropic clones and only 15 of 75 phenotypic R5-tropic clones were concordantly predicted. However, the remaining 60 phenotypic R5-tropic clones were discordantly predicted by at least one algorithm. In particular, 2 phenotypic R5-tropic clones were discordantly predicted by all algorithms tested. Taken together, the results demonstrate the limitation of currently available genotypic algorithms for predicting co-receptor inference among co-circulating multiple non-B subtypes and emerging isRFs. Also, the phenotypic tropism dataset presented here could be valuable for retraining of the widely used genotypic prediction algorithms to enhance their performance.

Longitudinal analysis of subtype C envelope tropism for memory CD4+ T cell subsets over the first 3 years of untreated HIV-1 infection

Background

HIV-1 infects a wide range of CD4⁺ T cells with different phenotypic properties and differing expression levels of entry coreceptors. We sought to determine the viral tropism of subtype C (C-HIV) Envelope (Env) clones for different CD4⁺ T cell subsets and whether tropism changes during acute to chronic disease progression. HIV-1 envs were amplified from the plasma of five C-HIV infected women from three untreated time points; less than 2 months, 1-year and 3-years post-infection. Pseudoviruses were generated from Env clones, phenotyped for coreceptor usage and CD4⁺ T cell subset tropism was measured by flow cytometry.

Results

A total of 50 C-HIV envs were cloned and screened for functionality in pseudovirus infection assays. Phylogenetic and variable region characteristic analysis demonstrated evolution in envs between time points. We found 45 pseudoviruses were functional and all used CCR5 to mediate entry into NP2/CD4/CCR5 cells. In vitro infection assays showed transitional memory (TM) and effector memory (EM) CD4⁺ T cells were more frequently infected (median: 46% and 25% of total infected CD4⁺ T cells respectively) than naïve, stem cell memory, central memory and terminally differentiated cells. This was not due to these subsets contributing a higher proportion of the CD4⁺ T cell pool, rather these subsets were more susceptible to infection (median: 5.38% EM and 2.15% TM cells infected), consistent with heightened CCR5 expression on EM and TM cells. No inter- or intra-participant changes in CD4⁺ T cell subset tropism were observed across the three-time points.

Conclusions

CD4⁺ T cell subsets that express more CCR5 were more susceptible to infection with C-HIV Envs, suggesting that these may be the major cellular targets during the first 3 years of infection. Moreover, we found that viral tropism for different CD4⁺ T cell subsets in vitro did not change between Envs cloned from acute to chronic disease stages. Finally, central memory, naïve and stem cell memory CD4⁺ T cell subsets were susceptible to infection, albeit inefficiently by Envs from all time-points, suggesting that direct infection of these cells may help establish the latent reservoir early in infection.

HIVCoR: A sequence-based tool for predicting HIV-1 CRF01_AE coreceptor usage

Determination of HIV-1 coreceptor usage is strongly recommended before starting the coreceptor-specific inhibitors for HIV treatment. Currently, the genotypic assays are the most interesting tools due to they are more feasible than phenotypic assays. However, most of prediction models were developed and validated by data set of HIV-1 subtype B and C. The present study aims to develop a powerful and reliable model to accurately predict HIV-1 coreceptor usage for CRF01_AE subtype called HIVCoR. HIVCoR utilized random forest and support vector machine as the prediction model, together with amino acid compositions, pseudo amino acid compositions and relative synonymous codon usage frequencies as the input feature. The overall success rate of 93.79% was achieved from the external validation test on the objective benchmark dataset. Comparison results indicated that HIVCoR was superior to other bioinformatics tools and genotypic predictors. For the convenience of experimental scientists, a user-friendly webserver has been established at http://codes.bio/hivcor/.

The Molecular Basis of pH-Modulated HIV gp120 Binding Revealed

Decades of research has yet to provide a vaccine for HIV, the virus which causes AIDS. Recent theoretical research has turned attention to mucosa pH levels over systemic pH levels. Previous research in this field developed a computational approach for determining pH sensitivity that indicated higher potential for transmission at mucosa pH levels present during intercourse. The process was extended to incorporate a principal component analysis (PCA)-based machine learning technique for classification of gp120 proteins against a known transmitted variant called Biomolecular Electro-Static Indexing (BESI). The original process has since been extended to the residue level by a process we termed Electrostatic Variance Masking (EVM) and used in conjunction with BESI to determine structural differences present among various subspecies across Clades A1 and C. Results indicate that structures outside of the core selected by EVM may be responsible for binding affinity observed in many other studies and that pH modulation of select substructures indicated by EVM may influence specific regions of the viral envelope protein (Env) involved in protein-protein interactions.

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.

Individuals with HIV-1 Subtype C Infection and Cryptococcal Meningitis Exhibit Viral Genetic Intermixing of HIV-1 Between Plasma and Cerebrospinal Fluid and a High Prevalence of CXCR4-Using Variants

The genotypic properties of human immunodeficiency virus type 1 (HIV-1) subtype C in individuals presenting with cryptococcal meningitis (CM) are not well established. Employing single-genome amplification as well as bulk PCR, cloning and sequencing strategies, we evaluated the genetic properties of HIV-1 subtype C env in 16 antiretroviral therapy-naive study participants with CM. Eleven of the 16 participants had matched blood plasma and cerebrospinal fluid (CSF) evaluated, with the rest having either a plasma or CSF sample evaluated. Before antiretroviral therapy initiation, matched plasma and CSF-derived env sequences of all 11 participants displayed genetic intermixing between the two compartments. Overall, 7 of the 16 (∼43.8%) participants harbored CXCR4-using variants in plasma and/or CSF, according to coreceptor usage prediction algorithms. This study suggests that HIV-1 subtype C genetic intermixing between peripheral blood and the central nervous system is common in individuals presenting with CM, and that CXCR4 usage is present in one or both compartments in approximately 44% of individuals.

Novel Strategy To Adapt Simian-Human Immunodeficiency Virus E1 Carrying env from an RV144 Volunteer to Rhesus Macaques: Coreceptor Switch and Final Recovery of a Pathogenic Virus with Exclusive R5 Tropism

The phase III RV144 human immunodeficiency virus (HIV) vaccine trial conducted in Thailand remains the only study to show efficacy in decreasing the HIV acquisition risk. In Thailand, circulating recombinant forms of HIV clade A/E (CRF01_AE) predominate; in such viruses, env originates from clade E (HIV-E). We constructed a simian-human immunodeficiency virus (SHIV) chimera carrying env isolated from an RV144 placebo recipient in the SHIV-1157ipd3N4 backbone. The latter contains long terminal repeats (LTRs) with duplicated NF-κB sites, thus resembling HIV LTRs. We devised a novel strategy to adapt the parental infectious molecular clone (IMC), R5 SHIV-E1, to rhesus macaques: the simultaneous depletion of B and CD8⁺ cells followed by the intramuscular inoculation of proviral DNA and repeated administrations of cell-free virus. High-level viremia and CD4⁺ T-cell depletion ensued. Passage 3 virus unexpectedly caused acute, irreversible CD4⁺ T-cell loss; the partially adapted SHIV had become dual tropic. Virus and IMCs with exclusive R5 tropism were reisolated from earlier passages, combined, and used to complete adaptation through additional macaques. The final isolate, SHIV-E1p5, remained solely R5 tropic. It had a tier 2 neutralization phenotype, was mucosally transmissible, and was pathogenic. Deep sequencing revealed 99% Env amino acid sequence conservation; X4-only and dual-tropic strains had evolved independently from an early branch of parental SHIV-E1. To conclude, our primate model data reveal that SHIV-E1p5 recapitulates important aspects of HIV transmission and pathobiology in humans.IMPORTANCE Understanding the protective principles that lead to a safe, effective vaccine against HIV in nonhuman primate (NHP) models requires test viruses that allow the evaluation of anti-HIV envelope responses. Reduced HIV acquisition risk in RV144 has been linked to nonneutralizing IgG antibodies with a range of effector activities. Definitive experiments to decipher the mechanisms of the partial protection observed in RV144 require passive-immunization studies in NHPs with a relevant test virus. We have generated such a virus by inserting env from an RV144 placebo recipient into a SHIV backbone with HIV-like LTRs. The final SHIV-E1p5 isolate, grown in rhesus monkey peripheral blood mononuclear cells, was mucosally transmissible and pathogenic. Earlier SHIV-E passages showed a coreceptor switch, again mimicking HIV biology in humans. Thus, our series of SHIV-E strains mirrors HIV transmission and disease progression in humans. SHIV-E1p5 represents a biologically relevant tool to assess prevention strategies.

Contribution of the gp120 V3 loop to envelope glycoprotein trimer stability in primate immunodeficiency viruses

The V3 loop of the human immunodeficiency virus type 1 (HIV-1) gp120 exterior envelope glycoprotein (Env) becomes exposed after CD4 binding and contacts the coreceptor to mediate viral entry. Prior to CD4 engagement, a hydrophobic patch located at the tip of the V3 loop stabilizes the non-covalent association of gp120 with the Env trimer of HIV-1 subtype B strains. Here, we show that this conserved hydrophobic patch (amino acid residues 307, 309 and 317) contributes to gp120-trimer association in HIV-1 subtype C, HIV-2 and SIV. Changes that reduced the hydrophobicity of these V3 residues resulted in increased gp120 shedding and decreased Env-mediated cell-cell fusion and virus entry in the different primate immunodeficiency viruses tested. Thus, the hydrophobic patch is an evolutionarily conserved element in the tip of the gp120 V3 loop that plays an essential role in maintaining the stability of the pre-triggered Env trimer in diverse primate immunodeficiency viruses.

•

The V3-loop of HIV-1 gp120 contributes to Env trimer stability and viral entry.

•

The hydrophobic patch in the tip of the V3 loop is critical for pre-triggered Env trimer stability.

•

The hydrophobic patch is a conserved motif in primate immunodeficiency viruses.

HIV-1 and SIV Predominantly Use CCR5 Expressed on a Precursor Population to Establish Infection in T Follicular Helper Cells

Background

T follicular helper (Tfh) cells are increasingly recognized as a major reservoir of HIV infection that will likely need to be addressed in approaches to curing HIV. However, Tfh express minimal CCR5, the major coreceptor for HIV-1, and the mechanism by which they are infected is unclear. We have previously shown that macaque Tfh lack CCR5, but are infected in vivo with CCR5-using SIV at levels comparable to other memory CD4⁺ T cells. Similarly, human splenic Tfh cells are highly infected with HIV-1 DNA. Therefore, we set out to examine the mechanism of infection of Tfh cells.

Methodology

Tfh and other CD4⁺ T cell subsets from macaque lymph nodes and spleens, splenic Tfh from HIV⁺ subjects, and tonsillar Tfh from HIV-uninfected subjects were isolated by cell sorting prior to cell surface and molecular characterization. HIV proviral gp120 sequences were submitted to genotypic and phenotypic tropism assays. Entry of CCR5- and CXCR4-using viruses into Tfh from uninfected tonsillar tissue was measured using a fusion assay.

Results

Phylogenetic analysis, genotypic, and phenotypic analysis showed that splenic Tfh cells from chronic HIV⁺ subjects were predominantly infected with CCR5-using viruses. In macaques, purified CCR5⁺PD-1^{intermediate(int)+} memory CD4⁺ T cells were shown to include pre-Tfh cells capable of differentiating in vitro to Tfh by upregulation of PD-1 and Bcl6, confirmed by qRT-PCR and single-cell multiplex PCR. Infected PD-1^int cells survive, carry SIV provirus, and differentiate into PD-1^hi Tfh after T cell receptor stimulation, suggesting a pathway for SIV infection of Tfh. In addition, a small subset of macaque and human PD-1^hi Tfh can express low levels of CCR5, which makes them susceptible to infection. Fusion assays demonstrated CCR5-using HIV-1 entry into CCR5⁺ Tfh and pre-Tfh cells from human tonsils.

Conclusion

The major route of infection of Tfh in macaques and humans appears to be via a CCR5-expressing pre-Tfh population. As the generation of Tfh are important for establishing effective immune responses during primary infections, Tfh are likely to be an early target of HIV-1 following transmission, creating an important component of the reservoir that has the potential to expand over time.

Sequential CCR5-Tropic HIV-1 Reactivation from Distinct Cellular Reservoirs following Perturbation of Elite Control

Background

HIV Elite Controllers may reveal insights into virus persistence given they harbour small reservoir sizes, akin to HIV non-controllers treated early with combination antiretroviral therapy. Both groups of patients represent the most promising candidates for interventions aimed at sustained remission or ‘cure’. Analytic treatment interruption (ATI) in the latter group leads to stochastic rebound of virus, though it is unclear whether loss of elite control is also associated with similar rebound characteristics.

Methods

We studied three discrete periods of virus rebound during myeloma related immune disruption over 2.5 years in an elite controller who previously underwent autologous stem cell transplantation (ASCT) in the absence of any antiretroviral therapy. Single genome sequencing of the V1-V4 region of env in PBMC and plasma was performed and phylogenies reconstructed. Average pairwise distance (APD) was calculated and non-parametric methods used to assess compartmentalisation. Coreceptor usage was predicted based on genotypic algorithms.

Results

122 single genome sequences were obtained (median 26 sequences per rebound). The initial rebounding plasma env sequences following ASCT represented two distinct lineages, and clustered with proviral DNA sequences isolated prior to ASCT. One of the lineages was monophyletic, possibly indicating reactivation from clonally expanded cells. The second rebound occurred 470 days after spontaneous control of the first rebound and was phylogenetically distinct from the first, confirmed by compartmentalisation analysis, with a different cellular origin rather than ongoing replication. By contrast, third rebound viruses clustered with second rebound viruses, with evidence for ongoing evolution that was associated with lymphopenia and myeloma progression. Following ASCT a shift in tropism from CXCR4-tropic viruses to a CCR5-tropic population was observed to persist through to the third rebound.

Conclusions

Our data highlight similarities in the viral reservoir between elite and non-controllers undergoing ATI following allogeneic transplantation. The lack of propagation of CXCR4 using viruses following transplantation warrants further study.

HIV-1 Group O Genotypes and Phenotypes: Relationship to Fitness and Susceptibility to Antiretroviral Drugs

Despite only 30,000 group O HIV-1 infections, a similar genetic diversity is observed among the O subgroups H (head) and T (tail) (previously described as subtypes A, B) as in the 9 group M subtypes (A-K). Group O isolates bearing a cysteine at reverse transcriptase (RT) position 181, predominantly the H strains are intrinsically resistant to non-nucleoside reverse transcriptase inhibitors (NNRTIs). However, their susceptibility to newer antiretroviral drugs such as etravirine, maraviroc, raltegravir (RAL), and elvitegravir (EVG) remains relatively unknown. We tested a large collection of HIV-1 group O strains for their susceptibility to four classes of antiretroviral drugs namely nucleoside RT, non-nucleoside RT, integrase, and entry inhibitors knowing in advance the intrinsic resistance to NNRTIs. Drug target regions were sequenced to determine various polymorphisms and were phylogenetically analyzed. Replication kinetics and fitness assays were performed in U87-CD4(+)CCR5 and CXCR4 cells and peripheral blood mononuclear cells. With all antiretroviral drugs, group O HIV-1 showed higher variability in IC50 values than group M HIV-1. The mean IC50 values for entry and nucleoside reverse transcriptase inhibitor (NRTI) were similar for group O and M HIV-1 isolates. Despite similar susceptibility to maraviroc, the various phenotypic algorithms failed to predict CXCR4 usage based on the V3 Env sequences of group O HIV-1 isolates. Decreased sensitivity of group O HIV-1 to integrase or NNRTIs had no relation to replicative fitness. Group O HIV-1 isolates were 10-fold less sensitive to EVG inhibition than group M HIV-1. These findings suggest that in regions where HIV-1 group O is endemic, first line treatment regimens combining two NRTIs with RAL may provide more sustained virologic responses than the standard regimens involving an NNRTI or protease inhibitors.

Highly Accurate Structure-Based Prediction of HIV-1 Coreceptor Usage Suggests Intermolecular Interactions Driving Tropism

HIV-1 entry into host cells is mediated by interactions between the V3-loop of viral glycoprotein gp120 and chemokine receptor CCR5 or CXCR4, collectively known as HIV-1 coreceptors. Accurate genotypic prediction of coreceptor usage is of significant clinical interest and determination of the factors driving tropism has been the focus of extensive study. We have developed a method based on nonlinear support vector machines to elucidate the interacting residue pairs driving coreceptor usage and provide highly accurate coreceptor usage predictions. Our models utilize centroid-centroid interaction energies from computationally derived structures of the V3-loop:coreceptor complexes as primary features, while additional features based on established rules regarding V3-loop sequences are also investigated. We tested our method on 2455 V3-loop sequences of various lengths and subtypes, and produce a median area under the receiver operator curve of 0.977 based on 500 runs of 10-fold cross validation. Our study is the first to elucidate a small set of specific interacting residue pairs between the V3-loop and coreceptors capable of predicting coreceptor usage with high accuracy across major HIV-1 subtypes. The developed method has been implemented as a web tool named CRUSH, CoReceptor USage prediction for HIV-1, which is available at http://ares.tamu.edu/CRUSH/.

Analysis of Immunological, Viral, Genetic, and Environmental Factors That Might Be Associated with Decreased Susceptibility to HIV Infection in Serodiscordant Couples in Florianópolis, Southern Brazil

Individuals who have been exposed to human immunodeficiency virus (HIV) and have not been infected might possess natural resistance mechanisms. An understanding of the sociodemographic and immunological conditions that influence resistance to HIV is a challenge, and very little is known about the role of intrinsic antiviral factors that restrict HIV infection. The aim of this study was to analyze potential factors responsible for resistance to HIV infection in serodiscordant couples by comparing HIV-exposed seronegative individuals (HESN) to HIV-seropositive individuals treated with antiretroviral therapy (HIV-ART) along with healthy controls (HC). The results revealed one HLA-B*27 and two HLA-B*57 individuals among the HESN; a CCR5Δ32 heterozygous deletion was observed in one serodiscordant couple, while the homozygous genotype for this variant was not observed. There were no differences in the basal mRNA expression of APOBEC3G, CFLAR, TRIM5α, LEDGF/p75, BST-2, or SAMHD1 in CD4(+) T lymphocyte- and monocyte-enriched populations among the three groups, and lower HBD-3 concentrations were observed in saliva from HIV-ART compared to HESN and HC. The most prevalent HIV-1 subtype was C or C-containing recombinant forms. Six HIV-ART individuals and one HIV-ART individual were infected with the R5 HIV and X4 HIV strains, respectively. The ability to control infection or delay disease progression is probably defined by a balance between viral and host factors, and further evaluation should be performed in larger cohorts. Our data suggest that susceptibility to HIV infection varies among individuals and strengthens the multifactorial characteristics underlying the resistance mechanisms in HIV.

Defining the fitness of HIV-1 isolates with dual/mixed co-receptor usage

Background

CCR5-using (r5) HIV-1 predominates during asymptomatic disease followed by occasional emergence of CXCR4-using (x4) or dual tropic (r5x4) virus. We examined the contribution of the x4 and r5 components to replicative fitness of HIV-1 isolates.

Methods

Dual tropic r5x4 viruses were predicted from average HIV-1 env sequences of two primary subtype C HIV-1 isolates (C19 and C27) and from two patient plasma samples (B12 and B19). Chimeric Env viruses with an NL4-3 backbone were constructed from the B12 and B19 env sequences. To determine replicative fitness, these primary and chimeric dual tropic HIV-1 were then competed against HIV-1 reference isolates in U87.CD4 cells expressing CXCR4 or CCR5 or in PBMCs ± entry inhibitors. Contribution of the x4 and r5 clones within the quasispecies of these chimeric or primary HIV-1 isolates were then compared to the frequency of x4, r5, and dual tropic clones within the quasispecies as predicted by phenotypic assays, clonal sequencing, and 454 deep sequencing.

Results

In the primary HIV-1 isolates (C19 and C27), subtype C dual tropic clones dominated over x4 clones while pure r5 clones were absent. In two subtype B chimeric viruses (B12 and B19), r5 clones were >100-fold more abundant than x4 or r5/x4 clones. The dual tropic C19 and C27 HIV-1 isolates outcompeted r5 primary HIV-1 isolates, B2 and C3 in PBMCs. When AMD3100 was added or when only U87.CD4.CCR5 cells were used, the B2 and C3 reference viruses now out-competed the r5 component of the dual tropic C19 and C27. In contrast, the same replicative fitness was observed with dualtropic B12 and B19 HIV-1 isolates relative to x4 HIV-1 A8 and E6 or the r5 B2 and C3 viruses, even when the r5 or x4 component was inhibited by maraviroc (or AMD3100) or in U87.CD4.CXCR4 (or CCR5) cells.

Conclusions

In the dual tropic HIV-1 isolates, the x4 replicative fitness is higher than r5 clones but the x4 or x4/r5 clones are typically at low frequency in the intrapatient virus population. Ex vivo HIV propagation promotes outgrowth of the x4 clones and provides an over-estimate of x4 dominance in replicative fitness within dual tropic viruses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12981-015-0066-7) contains supplementary material, which is available to authorized users.

Effect of lysine to arginine mutagenesis in the V3 loop of HIV-1 gp120 on viral entry efficiency and neutralization

HIV-1 infection is characterized by an ongoing replication leading to T-lymphocyte decline which is paralleled by the switch from CCR5 to CXCR4 coreceptor usage. To predict coreceptor usage, several computer algorithms using gp120 V3 loop sequence data have been developed. In these algorithms an occupation of the V3 positions 11 and 25, by one of the amino acids lysine (K) or arginine (R), is an indicator for CXCR4 usage. Amino acids R and K dominate at these two positions, but can also be identified at positions 9 and 10. Generally, CXCR4-viruses possess V3 sequences, with an overall positive charge higher than the V3 sequences of R5-viruses. The net charge is calculated by subtracting the number of negatively charged amino acids (D, aspartic acid and E, glutamic acid) from the number of positively charged ones (K and R). In contrast to D and E, which are very similar in their polar and acidic properties, the characteristics of the R guanidinium group differ significantly from the K ammonium group. However, in coreceptor predictive computer algorithms R and K are both equally rated. The study was conducted to analyze differences in infectivity and coreceptor usage because of R-to-K mutations at the V3 positions 9, 10 and 11. V3 loop mutants with all possible RRR-to-KKK triplets were constructed and analyzed for coreceptor usage, infectivity and neutralization by SDF-1α and RANTES. Virus mutants R₉R₁₀R₁₁ showed the highest infectivity rates, and were inhibited more efficiently in contrast to the K₉K₁₀K₁₁ viruses. They also showed higher efficiency in a virus-gp120 paired infection assay. Especially V3 loop position 9 was relevant for a switch to higher infectivity when occupied by R. Thus, K-to-R exchanges play a role for enhanced viral entry efficiency and should therefore be considered when the viral phenotype is predicted based on V3 sequence data.

Reliable genotypic tropism tests for the major HIV-1 subtypes

Over the past decade antiretroviral drugs have dramatically improved the prognosis for HIV-1 infected individuals, yet achieving better access to vulnerable populations remains a challenge. The principal obstacle to the CCR5-antagonist, maraviroc, from being more widely used in anti-HIV-1 therapy regimens is that the pre-treatment genotypic "tropism tests" to determine virus susceptibility to maraviroc have been developed primarily for HIV-1 subtype B strains, which account for only 10% of infections worldwide. We therefore developed PhenoSeq, a suite of HIV-1 genotypic tropism assays that are highly sensitive and specific for establishing the tropism of HIV-1 subtypes A, B, C, D and circulating recombinant forms of subtypes AE and AG, which together account for 95% of HIV-1 infections worldwide. The PhenoSeq platform will inform the appropriate use of maraviroc and future CCR5 blocking drugs in regions of the world where non-B HIV-1 predominates, which are burdened the most by the HIV-1 pandemic.

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.

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.

Coreceptor usage, diversity, and divergence in drug-naive and drug-exposed individuals from Malawi, infected with HIV-1 subtype C for more than 20 years

There are few cohorts of individuals who have survived infection with HIV-1 for more than 20 years, reported and followed in the literature, and even fewer from Africa. Here we present data on a cohort of subtype C-infected individuals from rural northern Malawi. By sequencing multiple clones from long-term survivors at different time points, and using multiple genotyping approaches, we show that 5 of the 11 individuals are predicted as CXCR4 using (by ≥3/5 predictors) but only one individual is predicted as CXCR4 using by all five algorithms. Using any one genotyping approach overestimates the number of predicted CXCR4 sequences. Patterns of diversity and divergence were variable between the HIV-1 long-term survivors with some individuals showing very small amounts of variation and change, and others showing a greater amount; both patterns are consistent with what has been described in the literature.

Alterations in HIV-1 gp120 V3 region are necessary but not sufficient for coreceptor switching in CRF07_BC in China

The most predominant HIV-1 strains in China's current epidemic is the Circulating Recombinant Form 07_BC (CRF07_BC). CRF07_BC is mainly considered as a CCR5-tropic (R5) virus, since CXCR4-tropic (X4) viruses have thus far not been found in this subtype, and the molecular determinants of coreceptor switching remain unknown. To investigate the mechanisms underlying coreceptor requirement in CRF07_BC viruses, we characterized a panel of pNL4-3-based chimeric viruses with mutated V3 loop regions derived from an HIV-1 CRF07_BC infectious clone pXJDC13. Among 17 chimeric viruses, seven were dual-tropic and induced syncytium formation in MT-2 cells. Two amino acid insertions between positions 13 and 14, as well as arginine substitution at position 11 or 16 (IG insertion and P16R mutation or MG insertion and S11R mutation), conferred the chimeric viruses CXCR4-tropic features, which were same as subtype C X4 viruses. Next, to construct CRF07_BC X4 variants, mutated V3 loops were cloned into the CRF07_BC infectious clone pXJDC13. These V3 loops, which in the pNL4-3 backbone conferred chimeric viruses with CXCR4-using ability, abrogated infectivity completely in the CRF07_BC pXJDC13 genetic background. Similarly, IG insertion or MG insertion and S11R mutation dramatically diminished or completely abolished viral infectivity in other envelopes of subtype C or CRF07_BC. These results suggest that the effects of IG insertion and P16R mutation or MG insertion and S11R mutation on CXCR4 usage are context dependent, and additional mutations elsewhere in the envelope are needed to compensate for these fitness-reducing alterations.

Viral tropism and antiretroviral drug resistance in HIV-1 subtype C-infected patients failing highly active antiretroviral therapy in Johannesburg, South Africa

Reports show that up to 30% of antiretroviral drug-naive patients in Johannesburg have CXCR4-utilizing HIV-1 subtype C. We assessed whether HIV-1 subtype C-infected individuals failing highly active antiretroviral therapy (HAART) have a higher proportion of CXCR4-utilizing viruses compared to antiretroviral drug-naive patients. The V3 loop was sequenced from plasma from 100 randomly selected HAART-failing patients, and tropism was established using predictive algorithms. All patients harbored HIV-1 subtype C with at least one antiretroviral drug resistance mutation. Viral tropism prediction in individuals failing HAART revealed similar proportions (29%) of X4-utilizing viruses compared to antiretroviral drug-naive patients (30%). Findings are in contrast to reports from Durban in which 60% of HAART-failing subjects harbored X4/dual/mixed-tropic viruses. Despite differences in proportions of X4-tropism within South Africa, the high proportion of thymidine analogue mutations (TAMs) and CXCR4-utilizing HIV-1 highlights the need for intensified monitoring of HAART patients and the predicament of diminishing drug options, including CCR5 antagonists, for patients failing therapy.

Quantifying susceptibility of CD4+ stem memory T-cells to infection by laboratory adapted and clinical HIV-1 strains

CD4⁺ T cells are principal targets for human immunodeficiency virus type 1 (HIV-1) infection. CD4⁺ T cell subsets are heterogeneous cell populations, divided by functional and phenotypic differences into naïve and memory T cells. The memory CD4⁺ T cells are further segregated into central, effector and transitional memory cell subsets by functional, phenotypic and homeostatic characteristics. Defining the distribution of HIV-1 infection in different T cell subsets is important, as this can play a role in determining the size and composition of the viral reservoir. Both central memory and transitional memory CD4⁺ T cells have been described as long-lived viral reservoirs for HIV. Recently, the newly described stem memory T cell subset has also been implicated as a long-lived HIV reservoir. Using green fluorescent protein (GFP) reporter strains of HIV-1 and multi parameter flow cytometry, we developed an assay to simultaneously quantify the susceptibility of stem memory (TSCM), central memory, effector memory, transitional memory and naïve CD4⁺ T cell subsets, to HIV-1 infection in vitro. We show that TSCM are susceptible to infection with laboratory adapted and clinical HIV-1 strains. Our system facilitates the quantitation of HIV-1 infection in alternative T cell subsets by CCR5- and CXCR4-using viruses across different HIV-1 subtypes, and will be useful for studies of HIV-1 pathogenesis and viral reservoirs.

Longitudinal Analysis of CCR5 and CXCR4 Usage in a Cohort of Antiretroviral Therapy-Naïve Subjects with Progressive HIV-1 Subtype C Infection

HIV-1 subtype C (C-HIV) is responsible for most HIV-1 cases worldwide. Although the pathogenesis of C-HIV is thought to predominantly involve CCR5-restricted (R5) strains, we do not have a firm understanding of how frequently CXCR4-using (X4 and R5X4) variants emerge in subjects with progressive C-HIV infection. Nor do we completely understand the molecular determinants of coreceptor switching by C-HIV variants. Here, we characterized a panel of HIV-1 envelope glycoproteins (Envs) (n = 300) cloned sequentially from plasma of 21 antiretroviral therapy (ART)-naïve subjects who experienced progression from chronic to advanced stages of C-HIV infection, and show that CXCR4-using C-HIV variants emerged in only one individual. Mutagenesis studies and structural models suggest that the evolution of R5 to X4 variants in this subject principally involved acquisition of an “Ile-Gly” insertion in the gp120 V3 loop and replacement of the V3 “Gly-Pro-Gly” crown with a “Gly-Arg-Gly” motif, but that the accumulation of additional gp120 “scaffold” mutations was required for these V3 loop changes to confer functional effects. In this context, either of the V3 loop changes could confer possible transitional R5X4 phenotypes, but when present together they completely abolished CCR5 usage and conferred the X4 phenotype. Our results show that the emergence of CXCR4-using strains is rare in this cohort of untreated individuals with advanced C-HIV infection. In the subject where X4 variants did emerge, alterations in the gp120 V3 loop were necessary but not sufficient to confer CXCR4 usage.