- Short report
- Open Access
Co-receptor tropism prediction among 1045 Indian HIV-1 subtype C sequences: Therapeutic implications for India
AIDS Research and Therapy volume 7, Article number: 24 (2010)
Understanding co-receptor tropism of HIV-1 strains circulating in India will provide key analytical leverage for assessing the potential usefulness of newer antiretroviral drugs such as chemokine co-receptor antagonists among Indian HIV-infected populations. The objective of this study was to determine using in silico methods, HIV-1 tropism among a large number of Indian isolates both from primary clinical isolates as well as from database-derived sequences.
R5-tropism was seen in 96.8% of a total of 1045 HIV-1 subtype C Indian sequences. Co-receptor prediction of 15 primary clinical isolates detected two X4-tropic strains using the C-PSSM matrix. R5-tropic HIV-1 subtype C V3 sequences were conserved to a greater extent than X4-tropic strains. X4-tropic strains were obtained from subjects who had a significantly longer time since HIV diagnosis (96.5 months) compared to R5-tropic strains (20.5 months).
High prevalence of R5 tropism and greater homogeneity of the V3 sequence among HIV-1 subtype C strains in India suggests the potential benefit of CCR5 antagonists as a therapeutic option in India.
After the discovery of the CD4 molecule as the major cellular receptor for HIV entry [1, 2], multiple studies suggested the presence of a secondary cellular receptor for HIV entry into the human CD4 cell [3, 4]. These co-receptors, particularly the chemokine receptors CCR5 and CXCR4, have been the subject of much research attempting to elucidate viral entry mechanisms, disease progression, antiretroviral therapy and vaccine development. Based on co-receptor usage, viral strains are classified into R5-tropic (those that use CCR5 for cellular entry), X4-tropic (those that use CXCR4) and dual tropic strains (those that use both co-receptors) . Co-receptor tropism of individual viral strains can be delineated using reporter cells expressing different coreceptors; however such cell-based assays are labor-intensive, expensive and not appropriate for high throughput screening . As an alternative, in silico strategies using computer simulation and bioinformatics have been developed to predict viral co-receptor usage from env gene sequence information [5–7]. Of late, the in silico approaches have been gaining popularity given the simplicity of this strategy and the fact that env sequences are increasingly becoming available globally.
The simplest method used for delineating HIV-1 tropism is known as the 'charge rule' , which relies on the charge of amino acids at positions 11 or 25 within the V3 loop when aligned against a consensus. Presence of positively charged amino acids (i.e. arginine, lysine, or histidine) in these positions typically is indicative of X4-tropism, while presence of other amino acid residues is associated with R5 phenotype . Currently a number of tools are available online to predict the co-receptor usage on the basis of the V3 sequence. Such tools include among others, (i) Geno2Pheno which predicts whether the corresponding virus is capable of using CXCR4 as a co-receptor (R5/X4 or X4 variants) or not (R5 variants) , (ii) the distant segments (ds)Kernel which include relative positional information of segments in a string of symbols which detects R5-, X4- and R5X4-tropic strains, and (iii) WebPSSM using CPSSM, a genotypic predictor based on position-specific scoring matrices (PSSM) which detects R5- or X4- tropic strains specially designed and validated for HIV-1 subtype C . Dual-tropic strains in C-PSSM are grouped with the X4- data set . Till date, molecular epidemiological information from India has indicated that > 96% HIV-1 circulating strains are HIV-1 subtype C (Geographic search interface, Los Alamos database, accessed on February 2010) . While X4-tropic HIV-1 subtype C strains have been widely reported from Africa [14–16], the presence of CXCR4 as a co-receptor to facilitate entry into the host cell is uncommon among Indian subtype C strains [17, 18]. R5-tropic viruses constitute by far the predominant strains in India although recent reports indicate the occasional presence of HIV-1 subtype C X4-tropic strains [19–21].
We aimed to characterize co-receptor tropism of HIV-1 subtype C strains isolated from a clinical cohort in southern India, using three different online bioinformatics tools. Furthermore, we aimed to validate this strategy and expand our understanding of co-receptor tropism preference among Indian strains by extending this analysis to a total of 1030 V3 sequences of Indian origin available at Los Alamos databank.
Study population and sample collection
A single peripheral blood sample was collected from 15 ART-naïve patients (10 males, 5 females) attending the Infectious Disease Clinic at St. John's Medical Hospital, Bangalore, between 1 and 30 November 2009. Patient characteristics are described in Table 1. Routine CD4 count was performed using a dual-platform flow cytometer (FACSCalibur, BD, USA). Genomic DNA from whole blood was extracted using a commercial kit (QIAamp Blood DNA kit, Qiagen, Germany).
Polymerase chain reaction and sequencing
The env gene portion encoding the V3-V5 region was amplified by nested polymerase chain reaction (PCR) from whole blood genomic DNA using iNtRON Taq Polymerase (Intron Biotech, South Korea). Primers were designed based on a consensus Indian sequence and modified from previously published reports [22, 23]. The first round of PCR was carried out with a forward primer, FP1: 5-CACCGGCTTAGGCATCTCCTATGGCAGGAAGAA-3 and reverse primer RP1: 5-TAACCCTTCCAGGTACCCCCTTTTCTTTTA-3. The nested PCR was carried out with the forward primer, FP2: 5'-tgtaaaacgacggccagtCTGTTAAATGGCAGTCTAGC and reverse primer, RP2: 5'-caggaaacagctatgaccCACTTCTCCAATTGTCCCTCA. Primers FP2 and RP2 contain the M13 universal primer sequence (lower case), which was used for population based sequencing.
HIV subtyping was carried out using three different tools i.e. REGA subtyping tools v2.0 http://www.bioafrica.net/rega-genotype/html/subtypinghiv.html, NCBI Viral Genotyping tools http://www.ncbi.nlm.nih.gov/projects/genotyping/formpage.cgi and RIP 3.0 http://www.hiv.lanl.gov/content/sequence/RIP/RIP.html.
Indian V3 sequence and co-receptor prediction
Sequences from our primary clinical isolates (n = 15) were pooled with 1030 Indian V3 sequences (additional file 1; accession numbers and sequence information) available from the Los Alamos database accessed on 7 Feb 2010. All available HIV-1 subtype C V3 sequences (n = 1114) were downloaded. Sequences containing a premature stop codon (n = 84) were excluded from the study. In all, 1045 sequences were analyzed in silico for co-receptor tropism using three different tools; (i) C-PSSM http://indra.mullins.microbiol.washington.edu/webpssm/, (ii) Geno2pheno [co-receptor] http://coreceptor.bioinf.mpi-inf.mpg.de/ and (iii) (ds)Kernel http://genome.ulaval.ca/hiv-dskernel/. To compare V3 characteristics of the R5 and X4 tropic strains, the consensus sequences of R5 and X4 tropic strains detected in C-PSSM were obtained using WebLogo http://weblogo.berkeley.edu/logo.cgi and consensus maker tool present in Los Alamos database http://www.hiv.lanl.gov/content/sequence/CONSENSUS/consensus.html. The mutation patterns of the cohort V3 loop sequences were compared with the consensus HIV-1 subtype C sequence and presented as additional file 2.
This study was approved by the Institutional Ethical Review Board of St. John's Medical College and Hospital. Written informed consent was obtained from each participant prior to sample collection.
All the 15 clinical isolates were detected as HIV-1 Subtype C; using three different subtyping tools, REGA subtyping tools v 2, NCBI Viral Genotyping tools and RIP 3.
Table 1 summarizes the predicted co-receptor tropism of each of the individual viral strains derived from our clinical cohort. While all the three tools predicted 13 out of the 15 viruses to be R5-tropic, there was disagreement among them with respect to two viral strains. The viral strain SJNAHS04 was predicted to be an X4 virus by C-PSSM while the other two tools characterized it to be a R5 virus. Similarly, the strain SJNAHS13 was predicted to be an R5 virus by (ds)Kernel while the other two tools found it to be an X4 virus. To see how reliable the three bioinformatics tools were in predicting co-receptor tropism of subtype C strains, we applied each of these tools to 1030 Indian env V3 sequences available at the Los Alamos database. This analysis confirmed a high magnitude of R5-tropism in the Indian env sequences predicted by all the three tools, 97%, 99% and 99.6% of R5-tropic sequences by CPSSM, Geno2Pheno and (ds) Kernel, respectively (Figure 1). As multiple sequences may be derived from the same patient, a second analysis was done after eliminating sequences that were derived from the same individual, and the result was similar to that obtained when all included sequences were analyzed.
The score in C-PSSM used for coreceptor tropism prediction was previously validated using sequences of the known syncytium inducing (SI) phenotypes on the MT2 cell line . C-PSSM score for SJNAHS04 (known history of HIV positivity > 5 years) was -16.09 and SJNAHS13 (known history of HIV positivity > 10 years) was -19.35; both were well above the prediction cutoff of -21.64 for X4-tropic viruses . The remaining isolates characterized to be R5-tropic were obtained from subjects with a recent history of HIV diagnosis within the past 5 years. The mean duration since detection of HIV-1 infection was longer for X4-tropic strains compared to R5-tropic strains (96.5 months and 20.5 months respectively), although only 2 X4-tropic strains were available for this analysis.
Sequence characteristics of Indian V3 sequence
Given that the emergence of X4 viruses is correlated with disease progression in subtype B infection, and that, although rare, several dual-tropic and X4-tropic viruses have been reported in subtype C from within and outside India, we sought to understand the relative magnitude of genetic variation between the V3 amino acid sequences of the predicted X4- and R5-tropic strains of the Indian origin. Consensus sequence logos for V3 amino acid sequences of 1012 R5-and 33 X4-tropic strains were determined using WebLogo v.3 http://weblogo.berkeley.edu/logo.cgi. This analysis identified a high degree of conservation within the key amino acid residues of V3 loop of the R5-tropic (Figure 2A), but not X4-tropic strains (Figure 2B). The V3 loop amino acid residues of the X4-tropic strains were highly variable (Figure 2B). An analysis of the V3 sequences of Indian subtype C, subtype B and subtype A/A1 irrespective of co-receptor tropism detected greater magnitude of sequence diversity in subtype B and A/A1 compared to V3 sequences of Indian HIV-1 subtype C (see additional file 3).
The overwhelming majority of Indian HIV-1 subtype C strains, 96.8% (1012 out of 1045) were predicted to be R5-tropic according to C-PSSM analysis. There were discrepancies in predicting HIV-1 subtype C X4-tropism as most of the tools were developed for HIV-1 subtype B strains. The C-PSSM represents an improvement over currently available methods for predicting X4-viruses in subtype C populations; this method has an estimated sensitivity of 81.8% and specificity of 93.3%, values considerably superior compared to other tools . Investigators studying co-receptor tropism prediction among African subtype C strains reported 100% concordance of C-PSSM with phenotypic assay in detecting CXCR4 . Our data supports the predominance of R5 phenotype in subtype C infected patients in India. Our data also revealed a correlation between X4-tropism and the duration since first diagnosis although this is to be interpreted with caution given that only 2 X4-tropic strains were studied. Both patients harboring probable X4-tropic virus were diagnosed with HIV-1 infection for longer than 5 years. Two previous reports of X4-tropic HIV-1 subtype C viruses from India have not commented on the length of infection [19, 20]. The presence of X4-tropic strains is well known to be significantly associated with longer duration of HIV infection . Variations at position 16 and 18 of the V3 loop in R5 viruses have been reported to lead to X4 tropism . In all the Indian R5 tropic HIV-1 subtype C sequences including cohort sequences (additional file 2), the GPGQ crown motif was significantly conserved, but in Indian X4 tropic strains, the 18th position was variable (Figure 2B). The overall conserved nature of HIV-1 subtype C V3 sequence may reduce the possibility of co-receptor switch in subtype C viruses and may partially explain the low prevalence of X4-tropic strains. An additional contributory factor for the high prevalence of R5-tropic strains may be the presence of a large pool of CCR5 positive CD4 cells in the Indian population which allows for the R5-strains to have improved replication fitness . CCR5 expression on CD4+ cells of HIV-1 infected individuals is higher among the Indian population (26.8%)  compared to the population in the USA (13.2%) . Two European studies showed varying levels of expression of CCR5 on CD4+ cells from HIV-infected individuals; 28% in Netherlands  and 17% in Italy .
The implementation of antiretroviral therapy (ART) in resource-limited settings requires use of standard first- and second-line therapies. CCR5 receptor antagonists such as maraviroc is a potential future option for second-line therapy in populations where R5-tropic strains predominate [33, 34]. The high proportion of R5-tropic strains and decreased evidence of co-receptor switch in HIV-1 subtype C viruses in India support the proposition that CCR5-antagonists may be promising drugs for future HIV treatment although concerns about potential overgrowth of X4-tropic strains need to be adequately addressed.
The present study, the first of its kind from India where a large number of env sequences were subjected to in silico co-receptor prediction analysis, revealed high prevalence of R5-tropism and greater homogeneity within the V3 loop sequences of HIV-1 subtype C Indian strains. Although prediction tools may not entirely substitute experimental evaluation, the simplicity of in silico strategies highlighted in this study can be a major advantage for coreceptor tropism prediction in resource-constrained settings. Furthermore, our findings also allude to the possibility of including CCR5 antagonists to the anti-retroviral repertoire with additional necessary precautions. The therapeutic implications of our findings are of global relevance and will facilitate further research on HIV-1 co-receptor usage and viral diversity.
Conflicts of interests
The authors declare that they have no competing interests.
Dalgleish A, Beverley P, Clapham P, Crawford D, Greaves M, Weiss R: The CD4 (T4) antigen is an essential component of the receptor for the AIDS retrovirus. Nature. 1984, 312: 763-767. 10.1038/312763a0
Klatzmann D, Champagne E, Chamaret S, Gruest J, Guetard D, Hercend T, Gluckman JC, Montagnier L: T-lymphocyte T4 molecule behaves as the receptor for human retrovirus LAV. Nature. 1984, 312: 767-768. 10.1038/312767a0
Feng Y, Broder CC, Kennedy PE, Berger EA: HIV-1 entry cofactor: Functional cDNA cloning of a seven-transmembrane, G protein-coupled receptor. Science. 1996, 272: 872-877. 10.1126/science.272.5263.872
Deng H, Liu R, Ellmeier W, Choe S, Unutmaz D, Burkhart M, Di Marzio P, Marmon S, Sutton RE, Hill CM, Davis CB, Peiper SC, Schall TJ, Littman DR, Landau NR: Identification of a major co-receptor for primary isolates of HIV-1. Nature. 1996, 381: 661-666. 10.1038/381661a0
Jensen M, 't Wout A: Predicting HIV-1 coreceptor usage with sequence analysis. AIDS Rev. 2003, 5: 104-112.
Skrabal K, Low AJ, Dong W, Sing T, Cheung PK, Mammano F, Harrigan PR: Determining human immunodeficiency virus coreceptor use in a clinical setting: degree of correlation between two phenotypic assays and a bioinformatic model. J Clin Microbiol. 2007, 45: 279-284. 10.1128/JCM.01118-06
Poveda E, Briz V, Quinones-Mateu M, Soriano V: HIV tropism: diagnostic tools and implications for disease progression and treatment with entry inhibitors. AIDS. 2006, 20: 1359-1367. 10.1097/01.aids.0000233569.74769.69
Raymond S, Delobel P, Mavigner M, Ferradini L, Cazabat M, Souyris C, Sandres-Sauné K, Pasquier C, Marchou B, Massip P, Izopet J: Prediction of HIV type 1 subtype C tropism by genotypic algorithms built from subtype B viruses. J Acquir Immune Defic Syndr. 2010, 53: 167-75. 10.1097/QAI.0b013e3181c8413b
Pillai S, Good B, Richman D, Corbeil J: A new perspective on V3 phenotype prediction. AIDS Res Hum Retroviruses. 2003, 19: 145-149. 10.1089/088922203762688658
Sing T, Low AJ, Beerenwinkel N, Sander O, Cheung PK, Domingues FS, Büch J, Däumer M, Kaiser R, Lengauer T, Harrigan PR: Predicting HIV co-receptor usage based on genetic and clinical covariates. Antivir Ther. 2007, 12: 1097-1106.
Boisvert S, Marchand M, Laviolette F, Corbeil J: HIV-1 coreceptor usage prediction without multiple alignments: an application of string kernels. Retrovirology. 2008, 5: 110- 10.1186/1742-4690-5-110
Jensen MA, Coetzer M, van 't Wout AB, Morris L, Mullins JI: A reliable phenotype predictor for human immunodeficiency virus type 1 subtype C based on Envelope V3 sequences. J Virol. 2006, 80: 4698-4704. 10.1128/JVI.80.10.4698-4704.2006
Geography Search Interface, Los Alamos database. http://www.hiv.lanl.gov/components/sequence/HIV/geo/geo.comp
Papathanasopoulos MA, Cilliers T, Morris L, Mokili JL, Dowling W, Birx DL, McCutchan FE: Full-length genome analysis of HIV-1 subtype C utilizing CXCR4 and intersubtype recombinants isolated in South Africa. AIDS Res Hum Retroviruses. 2002, 18: 879-886. 10.1089/08892220260190362
Abebe A, Demissie D, Goudsmit J, Brouwer M, Kuiken CL, Pollakis G, Schuitemaker H, Fontanet AL, Rinke de Wit TF: HIV-1 subtype C syncytium- and non-syncytiuminducing phenotypes and coreceptor usage among Ethiopian patients with AIDS. AIDS. 1999, 13: 1305-1311. 10.1097/00002030-199907300-00006
van Rensburg EJ, Smith TL, Zeier M, Robson B, Sampson C, Treurnicht F, Engelbrecht S: Change in co-receptor usage of current South African HIV-1 subtype C primary isolates. AIDS. 2002, 16: 2479-2480. 10.1097/00002030-200212060-00015
Cecilia D, Kulkarni SS, Tripathy SP, Gangakhedkar RR, Paranjape RS, Gadkari DA: Absence of coreceptor switch with disease progression in human immunodeficiency virus infections in India. Virology. 2000, 271: 253-258. 10.1006/viro.2000.0297
Sundaravaradan V, Das SR, Ramakrishnan R, Sehgal S, Gopalan S, Ahmad N, Jameel S: Role of HIV-1 subtype C envelope V3 to V5 regions in viral entry, co-receptor utilization and replication efficiency in primary T-lymphocytes and monocyte derived macrophages. Virol J. 2007, 4: 126- 10.1186/1743-422X-4-126
Kandathil AJ, Kannangai R, Abraham OC, Pulimood SA, Jensen MA, Sridharan G: HIV-1 with Predicted CXCR4 Genotype Identified in Clade C from India. Mol Diag Ther. 2009, 13: 19-24.
Neogi U, Sood V, Chowdhury A, Das S, Ramachandran VG, Sreedhar VK, Wanchu A, Ghosh N, Banerjea AC: Genetic analysis of HIV-1 Circulating Recombinant Form 02_AG, B and C subtype-specific envelope sequences from Northern India and their predicted co-receptor usage. AIDS Res Ther. 2009, 6: 28- 10.1186/1742-6405-6-28
Dash PK, Siddappa NB, Mangaiarkarasi A, Mahendarkar AV, Roshan P, Anand KK, Mahadevan A, Satishchandra P, Shankar SK, Prasad VR, Ranga U: Exceptional molecular and coreceptor-requirement properties of molecular clones isolated from an Human Immunodeficiency Virus Type-1 subtype C infection. Retrovirology. 2008, 5: 25- 10.1186/1742-4690-5-25
Bachmann MH, Delwart EL, Shpaer EG, Lingenfelter P, Singal R, Mullins JI: WHO Network on HIV Isolation and Characterization. Rapid genetic characterization of HIV-1 from four WHO-sponsored vaccine evaluation sites using a heteroduplex mobility assay. AIDS Res Hum Retroviruses. 1994, 10: 1343-1351. 10.1089/aid.1994.10.1345.
Gharu L, Ringe R, Pandey S, Paranjape R, Bhattacharya J: HIV-1 clade C env clones obtained from an Indian patient exhibiting expanded co-receptor tropism are presented with naturally occurring unusual amino acid substitutions in V3 loop. Virus Res. 2009, 144: 306-314. 10.1016/j.virusres.2009.04.020
REGA subtyping tools v2.0. http://www.bioafrica.net/rega-genotype/html/subtypinghiv.html
NCBI Viral Genotyping tools. http://www.ncbi.nlm.nih.gov/projects/genotyping/formpage.cgi
Recombination Identification Program (RIP 3.0)http://www.hiv.lanl.gov/content/sequence/RIP/RIP.html
Regoes RR, Bonhoeffer S: The HIV co-receptor switch: a population dynamical perspective. Trends Microbiol. 2005, 13: 269-77. 10.1016/j.tim.2005.04.005
Coetzer M, Nedellec R, Salkowitz J, McLaughlin S, Liu Y, Heath L, Mullins JI, Mosier DE: Evolution of CCR5 use before and during co-receptor switching. J Virol. 2008, 82: 11758-66. 10.1128/JVI.01141-08
Ramalingam S, Kannangai R, Vijayakumar TS, Subramanian S, Abraham OC, Rupali P, Jesudason MV, Sridharan G: Increased number of CCR5+ CD4 T cells among south Indian adults probably associated with the low frequency of X4 phenotype of HIV-1 in India. Indian J Med Res. 2002, 116: 90-95.
Ostrowski MA, Justement SJ, Catanzaro A, Hallahan CA, Ehler LA, Mizell SB, Kumar PN, Mican JA, Chun TW, Fauci AS: Expression of chemokine receptors CXCR4 and CCR5 in HIV-1-infected and uninfected individuals. J Immunol. 1998, 161: 3195-3201.
de Roda Husman AM, Blaak H, Brouwer M, Schuitemaker H: CC Chemokine Receptor 5 Cell-Surface expression in relation to CC chemokine receptor 5 genotype and the clinical course of HIV-1 infection. J Immunol. 1999, 163: 4597-4603.
Pierdominici M, Giovannetti A, Ensoli F, Mazzetta F, Marziali M, De Cristofaro MR, Santini-Muratori D, Leti W, Aiuti F: Changes in CCR5 and CXCR4 expression and [beta]-Chemokine production in HIV-1-infected patients treated with highly active antiretroviral therapy. JAIDS. 2002, 29: 122-131.
Castonguay L, Wengg Y, Adolfsen W, Di Salvo J, Kilburn R, Caldwell C, Daugherty BL, Finke PE, Hale JJ, Lynch CL, Mills SG, MacCoss M, Springer MS, Julie A: Binding of 2-aryl-4(piperidin-1yl)butanmines and 1, 3, 4-trisubstitued pyrrolidines to human CCR5: a molecular modelling-guide mutagenesis study of the binding pocket. Biochemistry. 2003, 42: 1544-1550. 10.1021/bi026639s
Soriano V, Clo-Federico P, Kaiser R, Calvez V, Gatell JM, di Perri G, Pillay D, Rockstroh J, Geretti AM: When and how to use maraviroc in HIV-infected patients. AIDS. 2009, 23: 2377-2385. 10.1097/QAD.0b013e328332d32d
We thank the staff of the Infectious Disease Clinic, St. John's Hospital, Bangalore, for assistance with patient recruitment. We are thankful to Mr. Soham Gupta and Ms. Pravat Nalini for their critical review of this manuscript. This study was partially funded by European Union Framework Program 7. Part of this work was presented at AIDS 2010, XVIII International AIDS Conference, Austria, Vienna, 18-23 July 2010.
UN designed the study, performed all laboratory tests and bioinformatics analysis. UN, ADC and AS drafted the manuscript. UN, AS, GDS, PBS obtained ethical approval and helped with subject recruitment. AS, VSK, GDS and PBS provided clinical expertise and UR provided scientific guidance. All authors have read and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Neogi, U., Prarthana, S.B., D'Souza, G. et al. Co-receptor tropism prediction among 1045 Indian HIV-1 subtype C sequences: Therapeutic implications for India. AIDS Res Ther 7, 24 (2010). https://doi.org/10.1186/1742-6405-7-24
- Subtyping Tool
- Individual Viral Strain
- Primary Clinical Isolate
- REGA Subtyping Tool
- Loop Amino Acid Residue