Background: Despite the existence of over 50 subtypes and circulating recombinant forms of HIV-1, subtype C dominates the heterosexual pandemic causing approximately 56% of all infections.
Objective: To evaluate whether viral genetic factors may contribute to the observed subtype-C predominance.
Methods: Chimeric viruses were generated using V1–V3 envelope fragments from a subtype-A/C dually infected woman with preferential genital replication of subtype C. Viral adaptation, spread and cell fusion ability were evaluated in vitro using peripheral blood mononuclear cells and HeLa-CD4-CCR5 cell lines, sequencing and cloning. Structural modeling was performed using a crystal structure of gp120-CD4-X5. Phylogenetic analysis was done using subtype-A, subtype-B and subtype-C sequences from blood and cervix of 37 infected women and database sequences.
Results: We identified two envelope motifs, compact V1–V2 loops and V3-316T, which are found at high frequency throughout subtype-C evolution and affect gp120 interactions with CD4 and CCR5, respectively. When a V1-Δ5 deletion or V3-A316T was incorporated into subtype A, each increased viral fusion and spread several fold in peripheral blood mononuclear cell and cell lines with low CCR5 expression. Structural modeling suggested the formation of an additional hydrogen bond between V3 and CCR5. Moreover, we found preferential selection of HIV with 316T and/or extremely short V1–V2 loops in cervices of three women infected with subtypes A/C, B or C.
Conclusion: As CD4+-CCR5+-T cells are key targets for genital HIV infection and cervical selection can favor compact V1–V2 loops and 316T, which increase viral infectivity, we propose that these conserved subtype-C motifs may contribute to transmission and spread of this subtype.
aLaboratory of Persistent Viral Diseases, Rocky Mountain Laboratories, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Hamilton, Montana, USA
bMRC Human Immunology Unit, Weatherall Institute of Molecular Medicine, University of Oxford, John Radcliffe Hospital, UK
cDivision of Structural Biology, Wellcome Trust Centre for Human Genetics, Oxford, UK
dHRC Pathogen Bioinformatics Unit, South African National Bioinformatics Institute, University of the Western Cape, Cape Town, South Africa
eDepartment of Infectious Diseases, Rigshospitalet, Copenhagen, Denmark.
Received 7 October, 2008
Revised 13 January, 2009
Accepted 2 February, 2009
Correspondence to Astrid K.N. Iversen, MRC HIU, WIMM, University of Oxford, John Radcliffe Hospital, Oxford OX3 9DS, UK. Tel: +44 0 1865 222498; fax: +44 0 1865 222502; e-mail: firstname.lastname@example.org
The current HIV pandemic is dominated by subtype C, which accounts for approximately 56% of the estimated 33 million infections worldwide, despite the concurrent existence of over 50 other subtypes and circulating recombinant forms (CRFs) . First identified in Ethiopia, subtype C has spread through southern and eastern Africa during the last two decades, invading areas formerly dominated by subtypes A and D and is currently causing devastating epidemics in parts of sub-Saharan Africa, Brazil, India and China [2–5].
These transmissions are largely through heterosexual contact or from mother to child. Although founder effects probably can explain the more recent subtype-C epidemics in India, Brazil and China, subtyping of early HIV sequences from southern and eastern Africa strongly suggest that subtype C was introduced into these areas subsequent to subtypes A and D (www.hiv.lanl.gov). As no epidemiological studies have shown any subtype-related differences between heterosexual transmission networks, it has been suggested that subtype-C virus may have a transmission advantage over subtypes A and D [6–8]. Moreover, epidemiological studies [9,10] have demonstrated differences between these subtypes with respect to vaginal shedding and in-utero transmission, subtype C associating with increased risks of both. However, no functional mechanism explaining these differences has so far been identified.
Although many studies have successfully focused on identifying human genetic factors influencing HIV-1 disease progression, these seem to account for just around 10% of the variability in disease progression rates and only some have been demonstrated to affect transmission . Viral genetic variation is also likely to influence pathogenesis and transmission, and short V1–V4 loops, and more neutralization-sensitive viruses, have been suggested to be selected during, or soon after, many sexual transmissions caused by subtypes A and C, but not B [12–16].
Here we asked whether conserved, subtype-specific sequence motifs in subtype C could affect viral replication and transmission. Our results highlight the functional importance of two common subtype-C motifs, short V1–V2 loops and V3-316T, and demonstrate that these are found at high frequency throughout subtype C evolution and may be specifically selected in the genital tract of infected women.
HIV-1-infected women were recruited at Rigshospitalet, Denmark, between 1995 and 1996. Most were untreated; however, some had received monotherapy [azidothymidine (AZT) (8/37), didanosine (ddI) (1/37)] or dual therapy [AZT and ddI (2/37)] for over 6 months [17,18]. The study was approved by the Danish Board of Medical Ethics and patients gave informed consent .
Limiting dilution, nested PCR and sequencing
DNA was extracted from peripheral blood mononuclear cell (PBMC) and cervical cell pellets using the QIAamp Blood kit (Qiagen, Valencia, California, USA). V1–V3, C2–C4 and gp160 envelope fragments were amplified by limiting-dilution nested PCR using Advantage 2 Polymerase mix (Clontech, Mountain View, California, USA) and sequenced (sequences are available upon request) .
Construction of recombinant HIV-1 clones
Infectious HIV-A and HIV-C chimeric clones with V1–V3 env fragments obtained from blood or cervix of a dually infected patient were generated as described previously for HIV clones 81A-4 and 49-5 [20–23].
PCR amplification, sequencing and cloning of PBMC-adapted virus
Virus-containing supernatants were samples throughout the PBMC infection experiment. Viral RNA was purified using the QIAamp Viral RNA Mini Kit (Qiagen, Valencia, California, USA), reverse transcribed and amplified using the SuperScript One-Step RT-PCR Kit for Long Templates (Invitrogen, Paisley, UK), sequenced and cloned. The reconstructed PBMC-adapted chimeric subtype-A clones were named A20-V1Δ5, A2-T146I, A1-A316T and A1-N302K, according to envelope region and type of mutation.
Virus stocks and infection
Proviral HIV constructs were transfected into RC49 HeLa cells and viral stocks were generated for PBMC, HeLa cell and fusion experiments; titers were expressed as focus-forming units (FFUs) per ml [24,25]. Different virus dilutions [multiplicity of infection (MOI) 0.004–0.1] were used to infect PBMC and cultures were continued for 16 days; p24 levels were measured using an antigen capture assay (HIV-1 p24CA) as recommended (AIDS Vaccine Program, National Cancer Institute (NCI)-Frederick Cancer Research and Development Center). Infected HeLa cell monolayers were stained for p24 and foci counted as described previously [24,25].
Fusion experiments and immunohistochemical detection
293T cells, which do not express CD4 or CCR5, were transfected with subtype-A, subtype-C or adapted-A chimeric HIV constructs and relative envelope expression levels compared using pooled antisera from HIV-infected patients and fluorescently labeled antihuman IgG (Invitrogen). Data were acquired on a FACSCalibur flow cytometer (Becton Dickinson, Franklin Lakes, USA). Transfected 293T cells were subsequently overlaid onto JC37, or RC49 monolayers for 24 h and stained for p24 antigen.
Phylogenetic trees were constructed by maximum likelihood using PAUP* 4.0b10 (Sinauer Ltd., Sunderland, Massachusetts, USA) and phyML software (www.atgc.lirmm.fr/phyml/) [26,27]. Nonparametric bootstrap support estimates were generated using 1000 replicates under neighbor-joining tree algorithm using the maximum likelihood substitution model. Tree figures were produced using FigTree1.1 (tree.bio.ed.ac.uk/software/figtree/). Distance analysis was done using MEGA and Kimura-2 parameters; standard errors were estimated using 100 bootstrap runs . The variation analyses were computed using a counting algorithm (http://bioafrica.mrc.ac.za/JoeTool/counting.html). Frequencies of 316A and T were calculated using 1671 subtype-A and 4165 subtype-C sequences. To decrease the risk of sampling bias in these database-derived sequences, we used only one sequence from each patient or transmission cluster from as many countries as possible (i.e. multiple sequences from known transmission clusters such as wife/husband or mother/baby were excluded). Sequences with translational problems (i.e. frameshift or stop codons) were furthermore excluded as these may either represent nonfunctional virus or could indicate poor sequence quality.
Molecular modeling and disorder analysis
The A316T and N302K single amino acid point mutations were modeled on the crystal structure of gp120-CD4-X5  (PDB ID 2B4C) by manually substituting the individual side chains using Coot , ensuring by visual inspection that the resulting model maintained good local stereochemistry (Supplementary Fig. S1). The increase in local positive charge at the region of the N302K mutation was confirmed by mapping the molecular electrostatic potential onto the solvent-accessible surface of gp120 (PDB ID 2B4C) and the gp120 N302K model (not shown) using the Adaptive Poisson-Boltzmann Solver (APBS)  plug-in for PyMOL (www.pymol.sourceforge.net/; DeLano Scientific) with default parameters. The propensity of gp120 residues to be disordered was predicted from the sequences of primary isolates and adapted strains using the neural-network pattern recognition program regional order neural network (RONN)  with default parameters. Molecular graphics were prepared using PyMOL (DeLano Scientific).
Generation of infectious HIV molecular clones
We previously reported preferential genital replication of HIV subtype C in an asymptomatic, untreated woman dually infected by subtypes A and C in blood; her blood-derived and cervix-derived C sequences were indistinguishable and subtype A appeared to have infected her prior to subtype C .
We hypothesized that preferential cervical subtype-C infection was associated with viral genetic traits favoring genital infection. To compare the biological characteristics of HIV sequences found in blood and cervix, we amplified and sequenced single-copy V1–V3 envelope fragments from these tissues, as this region contains key determinants of cell tropism, infectivity and chemokine sensitivity [34–36] (Fig. 1a). Sequences were inserted into the pNL4-3ΔV1–V3 vector to generate infectious molecular clones (subtype C: C57, C64, C66, C68, C69 and C85; subtype A: A1, A2 and A20). Infection studies using HeLa cells expressing CD4 and CXCR4 plus or minus CCR5 showed that all virus clones required the CCR5 coreceptor for infection. Furthermore, all viruses were able to infect and spread slowly in monocyte-derived macrophage cultures (data not shown).
Bioinformatic analysis demonstrated that these sequences were common subtype representatives as they clustered within their respective subtypes on phylogenetic trees, and the distance and amino acid variation between them and other subtype-A and subtype-C sequences were comparable to the medium intrasubtype distance and variation (supplementary Table S1, data not shown). Thus, the viral sequences from our patient may provide information on how subtypes A and C compare functionally in vivo.
In-vitro PBMC replication and cell fusion
In PBMC cultures, we noted marked subtype-specific differences in live virus production. Whereas subtype-C titers peaked at 4–6 days postinfection, subtype-A titers rose more slowly over the course of 16 days (Fig. 1b). Throughout the PBMC infection experiment, viral supernatants were removed and used to infect JC37 cells, which expressed low levels of CCR5 and high levels of CD4. Regardless of sampling time, supernatants from the subtype-C infections produced large foci (10–20 nuclei/focus; Fig. 1d). In contrast, supernatants from the subtype-A infections at days 4–8 generated mostly small foci (1–2 nuclei/focus). However, subtype-A foci increased markedly in size at days 12–16, coincident with increased replication in PBMC, suggesting that adaptation occurred during in-vitro culture of these viruses.
We examined the ability of the viral envelope protein to fuse with cell lines expressing CD4 and CCR5 by transfecting 293T (CD4−, CCR5−) cells with plasmid DNA of subtype-A and subtype-C infectious clones. The transfected cells expressed similar cell surface viral envelope levels by flow cytometry analysis at 48 h posttransfection. When these cells were overlaid for 24 h on RC49 cells (CD4low, CCR5high), both subtypes A and C induced a high level of fusion (Fig. 1e). By contrast, when they were overlaid onto JC37 cells (CD4high, CCR5low), subtype-A viruses generated small HIV-positive fusion foci ranging in size from 1 to 9 nuclei (mean ∼2 nuclei per focus), whereas subtype-C foci contained from 1 to more than 100 nuclei per focus (mean ∼20 nuclei). Thus, HeLa cells with low CCR5 levels were fused more efficiently by subtype-C viruses than by subtype A viruses, but both subtypes resulted in extensive fusion on cells expressing high CCR5 levels.
Selection of adapted subtype-A mutants in PBMC
During replication in PBMC (Fig. 1b), subtype-A viruses attained the high titers seen early in subtype C infection only slowly, and the early and late supernatants from subtype A infections differed in their ability to fuse JC37 cells, suggesting that mutant subtype A viruses had been selectively amplified during in-vitro culture. To search for evidence of adaptation, we reverse transcribed, amplified and sequenced viral subtype-A RNA from day 16 and found four distinct envelope sequence changes targeting the V1 and V3 loops. Using these sequences, we generated new chimeric viruses, A1-A316T, A1-N302K, A2-T146I and A20-V1Δ5, which, apart from the indicated changes, were otherwise identical to their parental virus (Fig. 1a). These adapted viruses displayed a distinct functional phenotype, which differed from their parental subtype-A clones but was similar to the chimeric subtype-C clones, that is, they spread rapidly in PBMC and produced large foci in JC37 cells (Fig. 1c,f), cell types that both express low levels of CCR5.
The influence of CCR5 expression levels on virus titers
To determine whether CCR5 expression levels influenced viral replication in CD4+, CCR5+ HeLa cells, subtype-A, subtype-C and PBMC-adapted subtype-A clones were used to infect JC37 and RC49 cells. In RC49 (CD4low, CCR5high) cells, all viruses gave high titers (2.5 × 105 − 3.7 × 106 FFU/ml) at 5 days postinfection (Table 1). In contrast, in JC37 (CD4high, CCR5low) cells, all three subtype-A clones produced low-virus titers (1.9 × 103 − 1.7 × 105 FFU/ml), whereas subtype-C and all four PBMC-adapted subtype-A clones produced high titers similar to those seen on RC49 cells.
Structural modeling of envelope variation
Structural modeling was performed to help understand the molecular consequences of the four subtype-A in-vitro adaptations (Supplementary Fig. S1). The V1–V2 envelope region lies close to the CD4-binding site, whereas the V3 region interacts with two parts of CCR5 (Supplementary Fig. S1a). The V3-A316T substitution may strengthen HIV-CCR5 interactions as the additional polar group is in an appropriate position to form a new hydrogen bond to CCR5, most likely independently of the surrounding genomic context (Supplementary Fig. S1b). Similarly, the N302K substitution adds another positive charge to the V3 stem, presumably strengthening CCR5 binding, as a crystal structure of V3 and the N-terminal negatively charged CCR5 tail has shown it to interact directly with position 302 .
The interactions between the CD4 receptor and the V1–V2 envelope region are less precisely understood as no high-resolution structures of V1–V2-containing gp120 envelopes exist. The loss of a V1-glycosylation site in A2-T146I and the A20-V1Δ5 deletion might enhance CD4 receptor contact, perhaps by decreasing steric hindrance. We used a neural-network technique to analyze the propensity for disorder of the envelope proteins, as disordered regions can combine high-binding specificity with modest affinity [32,38] (Supplementary Fig. S1c). We found significant differences in disorder probability between subtypes A, A20-V1Δ5 and C in the V1–V2 region, suggesting variation in the way these sequences might bind to CD4 and/or conformationally mask the CD4-binding site; however, both wild-type-A and wild-type-C chimeric viruses utilized the CD4 receptor efficiently in our functional assay.
Analysis of in-vivo occurrence of mutations selected in vitro
We examined the frequency of the four subtype-A in-vitro adaptations in natural HIV infections. Two substitutions were found occasionally in natural infections in some subtypes, however, rarely (146I, 2%) or never (302K) in subtype A (Table 2). By contrast, V3-316T was seen in 68% of all subtype-C sequences, but at much lower frequencies, if at all, in the other subtypes (25–0%) (Table 2, Fig. 2). Furthermore, the V1-Δ5 subtype-A deletion, which shortens the loop, resembled the short V1–V2 loops typically found in subtype C but not in subtypes A, B and D (Table 2). The high frequency of V3-316T and short V1–V2 loops in subtype C, which both increase HIV spread in vitro, suggest these motifs may affect HIV replication in infected patients.
Phylogenetic analysis of 316T selection
While deletions cannot be incorporated in phylogenetic analysis at present, we could examine V3-316T selection. We generated a phylogenetic tree based on subtype-C sequences from 1985 to 2005 and superimposed the nature of the amino acid at position 316 (Fig. 2a, Supplementary Fig. S2). The earliest sequences are from the Democratic Republic of Congo (DRC) where HIV is thought to have evolved and diversified. The phylogenetic tree demonstrates that both 316A and 316T could be found in the DRC sequences at the base of the tree. However, the deep lineages and branches of the tree mostly carry 316T, suggesting that this residue has been generally advantageous to subtype-C spread. Whereas 316V can be found sporadically, alanine seems to spread locally, especially in more recent sequences from South Africa and Botswana. We studied the distribution of 316A and 316T in larger alignments of shorter envelope sequences dating from 1990 to 2005 to see if this observation could be substantiated (Fig. 2b). We found that the frequency of 316T indeed had decreased from around 80% in 1990 to approximately 68% in 2005. This suggests that there may be alternate selective advantages for 316A in some situations in vivo and/or that other changes sometimes are sufficient to favor transmission of subtype C.
Analysis of HIV-316T distribution in blood and cervix
To analyze the prevalence of 316T in blood and cervix, we obtained samples from 36 additional women infected by HIV subtypes A (n = 3), B (n = 28) or C (n = 5). Whereas single-copy envelope sequences could be amplified from PBMC from all women, no HIV was found in the cervices of four .
Most subtype-B infected women carried virus with V3-316A, but four had either 316T, 316T/M, 316S or 316V virus in both samples (Fig. 3a). Moreover, we found V3-316 compartmentalization between blood and cervix in one subtype-B-infected woman (1/28, 4%); all genital viruses carried 316T (5/5), whereas only 316A could be found in PBMC (7/7) (Fisher's exact test, P < 0.001; Fig. 3b).
Of the five subtype-C-infected women, only one carried HIV-316A, and no virus could be amplified from her genital sample; the rest had either 316T or 316T/M in both PBMC and cervix. In one of these patients with HIV-316T, we were able to analyze gp160 from vagina as well as cervix and PBMC (Fig. 3c). Although 28% of the cervical sequences had very short V1–V2 loops (55 amino acids), such short loops were not found in blood; most vaginal sequences had intermediate V1–V2 loops. Overall, V1–V2 loop lengths differed significantly between genitalia and blood when grouped into either short (≤69 amino acids) or long (>69 amino acids) loops (Fisher's exact test, P = 0.02).
While our results demonstrate that subtype-A and subtype-B virus with 316A can replicate in the cervix, they also highlight that selective forces in the genital tract may favor replication of virus with 316T even in subtype-B-infected women. This is reminiscent of the preferential cervical replication of subtype-C virus in our dually infected woman. Furthermore, we demonstrate that in a woman carrying only HIV-316T, selective forces in the genital tract may favor shorter V1–V2 loops than selection in blood.
Subtype C is responsible for an ever-increasing fraction of new HIV infections worldwide, especially in areas where heterosexual transmission predominates, and it now prevails even in regions where subtype A or D or both used to be common [1,2,39,40]. This suggests that there are biological factors favoring transmission of subtype C, but no plausible mechanism has been suggested.
Our study of infectious, chimeric HIV clones carrying V1–V3 sequences from a woman dually infected with subtypes A and C, but with preferential cervical subtype-C replication, combined with the analysis of HIV sequences from blood and cervix from 36 additional women, lead to four findings suggesting that subtype C generally has a relative replication advantage compared with most subtype-A and subtype-B viruses in the female genital tract.
First, we identify a V3-envelope signature sequence (316T) at high frequency in subtype C, which might strengthen the interaction between HIV and CCR5, the main coreceptor in the female genital tract [41–45]. We find that the effect of 316T may be independent of the surrounding genomic context and show the effect of an A316T substitution in subtype A. Second, we find that selection for 316T can occur even in the cervix of a minority of untreated subtype-B-infected women. Third, we demonstrate that a subtype-A V1-deletion augments virus production in cell lines expressing low levels of CD4 or CCR5 and propose that the shortening of the V1 loop increases viral interactions with the CD4 receptor. Fourth, we find that genital selection in a subtype-C-infected woman favors outgrowth of very short V1–V2 loops, despite subtype C already having the most compact V1–V2 loops of all HIV subtypes.
We found that subtype C had a greater ability to replicate in PBMC cultures than subtype A and mediated more extensive virus-induced fusion of cells expressing low CCR5 levels. However, we observed marked changes in subtype-A replication patterns over time and identified adaptations in either the V1 (T114I, V1-Δ5) or V3 (N302K, A316T) region, which are known to interact with the CD4 and CCR5 receptors, respectively. These adaptations enhanced virus replication and fusion comparable to what we found using the wild-type subtype-C clones. Whereas two were found rarely, if at all, in HIV sequences obtained from infected persons, two others (A316T and V1-Δ5) resembled sequence motifs in the subtype C sequences from our patient and in most subtype C sequences from natural infections. In other subtypes, and in most CRFs, alanine is preferred at position 316, with the notable exception of CRF_02_AG, which carries threonine. Analogous to subtype-C spread in southern and eastern Africa, CRF_02_AG has spread throughout western Africa, primarily via heterosexual contact [1,46].
Both subtype-A and subtype-C chimeric viruses interacted efficiently with the CD4 receptor; however, the V1-Δ5 adaptation augmented this interaction and increased infectivity several-fold in both cell lines and PBMC. The CD4-binding region is predicted to be shielded by the V1–V2 loop and shorter loops have been associated with increased susceptibility to antibody neutralization [13–15]. As short V1–V4 loops and more neutralization-sensitive HIV variants have been found in early subtype-A and subtype-C infections, it was suggested that these were either preferentially transmitted or rapidly selected in the recipient [12–15]. We demonstrate that genital selection pressures in subtype-C-infected women may sometimes favor virus with more compact V1–V2 loops than what is commonly found in blood. Consequently, whether transmission of such variants is preferential or at times simply occurs by chance, remains to be elucidated.
Structural modeling suggested that the A316T effect might be independent of the surrounding genomic background, as a new hydrogen bond may form between the polar threonine and CCR5. The importance of position 316 in gp120-CCR5 interactions is indicated by the induction of drug-resistance mutations at this site by two different CCR5 inhibitors (A316T, maraviroc; A316V, AD101) [47,48]. CCR5-binding studies [47,48] have demonstrated that these polar and nonpolar substitutions crucially affect HIV-CCR5 interactions through different mechanisms.
Our phylogenetic analysis demonstrated that 316T must have been favorable for subtype-C transmission. However, it also revealed that 316A seems to spread locally, especially in recent sequences from South Africa and Botswana. This implies that viruses with 316T incur a viral fitness cost even in subtype C, as has been directly demonstrated for maraviroc-resistant subtype-B virus with 316T in in-vitro replication assays . Additionally, in some situations in vivo, there could be selective advantages associated with 316A and/or other changes resulting in effective transmission. The A316T substitution could also affect virus susceptibility to host antibodies or cytotoxic T cell (CTL) responses, as the V3 crown is a prime target for neutralizing antibodies and encompasses CTL epitopes restricted by at least five different human leukocyte antigen (HLA) alleles.
Whereas HIV carrying 316A is fully replication competent and can be found in cervix, virus with 316T may sometimes have a relative replication advantage as CCR5 expression has been found on cells sampled from vagina, ectocervix and endocervix [41–45,49], the proportion of cervicovaginal T cells expressing CCR5 is markedly expanded as compared with peripheral blood , HIV primarily propagates in vagina by T-cell infection  and cell-to-cell transfer is the predominant mode of HIV spread . This relative replication advantage may be greatest in women with low genital CCR5-expression levels such as premenopausal women with normal genital mucosa, and less in women with high expression levels due to, for instance, sexually transmitted diseases (STD), cellular infiltrates, menopause or oral contraceptive use [42,52]. Indeed, our patients with preferential cervical replication of HIV subtype C or B with 316T had normal genital epithelium and no STD .
We propose that the efficient HIV-316T-CCR5 interaction may increase viral transmissibility. That a polymorphism at a single amino acid site can have such profound consequences agrees with the first report of species-specific selection in primate lentiviruses, which demonstrates that efficient spread of SIV cpz in human cells depends on a single amino acid change in Gag . Our proposal is indirectly supported by transmission studies in discordant couples with known CCR5 genotype; if the uninfected partner was heterozygous for the Δ32-CCR5 allele and thus had lower levels of CCR5 on T cells and monocytes , HIV transmission decreased during both heterosexual  and homosexual [43,44] contact, although one study found no effect .
Our results raise the possibility of subtype-specific differences in the efficacy of topical genital CCR5 inhibitors and highlight the need for clinical trials on patients infected with different subtypes. Furthermore, they suggest that the maraviroc-induced drug mutation A316T could result in viruses that are more easily transmitted through sexual contact.
In conclusion, our study suggests that subtype C has evolved in ways compatible with increased sexual transmissibility. Although other envelope regions, additional viral genes and long terminal repeat promotor differences  in combination with epidemiological dynamics, such as founder effects, also must be important factors, these subtype-C motifs may contribute to its current dominance in southern and eastern Africa, which facilitated subsequent worldwide spread.
We thank the patients, Andrew Rambaut, Kim Hasenkrug, Rachel Lacasse, Angela Vincent and Nick Willcox for helpful comments; Anita Mora for graphics, Kevin Braughton and Ron Messer for technical assistance.
A.K.N.I. initiated and designed the overall study; A.K.N.I. and B.C. planned and supervised different parts of the project; A.K.N.I., B.C. and B.L.W. analyzed data and wrote the article; B.L.W., A.E.A. and A.K.N.I. performed experiments; S.C.G., E.Y.J. and D.I.S. performed the structural modeling and RONN analyses; T.d.O. and A.K.N.I. performed the phylogenetic analysis; A.K.N.I. and P.S. provided patients, clinical data and reagents; and A.J.M. contributed reagents. All authors discussed the results and approved the final version of the article.
Funding was provided by the Danish AIDS Foundation, the Nuffield Dominions Trust, Cancer Research UK, the EU (LSHG-CT-2006-031220) and MRC, UK. The authors have no competing financial interests. Supplementary documentation can be obtained on request.
1. Osmanov S, Pattou C, Walker N, Schwardlander B, Esparza J. Estimated global distribution and regional spread of HIV-1 genetic subtypes in the year 2000. J Acquir Immune Defic Syndr 2002; 29:184–190.
2. Morison L, Weiss HA, Buve A, Caraël M, Abega SC, Kaona F, et al
. Commercial sex and the spread of HIV in four cities in sub-Saharan Africa. AIDS 2001; 15(Suppl 4):S61–S69.
3. Shankarappa R, Chatterjee R, Learn GH, Neogi D, Ding M, Roy P, et al
. Human immunodeficiency virus type 1 env sequences from Calcutta in eastern India: identification of features that distinguish subtype C sequences in India from other subtype C sequences. J Virol 2001; 75:10479–10487.
4. Yu XF, Wang X, Mao P, Wang S, Li Z, Zhang J, et al
. Characterization of HIV type 1 heterosexual transmission in Yunnan, China. AIDS Res Hum Retroviruses 2003; 19:1051–1055.
5. Soares EA, Martinez AM, Souza TM, Santos AF, Da Hora V, Silveira J, et al
. HIV-1 subtype C dissemination in southern Brazil. AIDS 2005; 19(Suppl 4):S81–S86.
6. Carswell W. HIV in South Africa. Lancet 1993; 342:132.
7. Hunter DJ. AIDS in sub-Saharan Africa: the epidemiology of heterosexual transmission and the prospects for prevention. Epidemiology 1993; 4:63–72.
8. Morison L, Buve A, Zekeng L, Heyndrickx L, Anagonou S, Musonda R, et al
. HIV-1 subtypes and the HIV epidemics in four cities in sub-Saharan Africa. AIDS 2001; 15(Suppl 4):S109–S116.
9. Renjifo B, Gilbert P, Chaplin B, Msamanga G, Mwakagile D, Fawzi W, et al
. Preferential in-utero transmission of HIV-1 subtype C as compared to HIV-1 subtype A or D. AIDS 2004; 18:1629–1636.
10. John-Stewart GC, Nduati RW, Rousseau CM, Mbori-Ngacha DA, Richardson BA, Rainwater S, et al
. Subtype C is associated with increased vaginal shedding of HIV-1. J Infect Dis 2005; 192:492–496.
11. O'Brien SJ, Nelson GW. Human genes that limit AIDS. Nat Genet 2004; 36:565–574.
12. Chohan B, Lang D, Sagar M, Korber B, Lavreys L, Richardson B, et al
. Selection for human immunodeficiency virus type 1 envelope glycosylation variants with shorter V1-V2 loop sequences occurs during transmission of certain genetic subtypes and may impact viral RNA levels. J Virol 2005; 79:6528–6531.
13. Li B, Decker JM, Johnson RW, Bibollet-Ruche F, Wei X, Mulenga J, et al
. Evidence for potent autologous neutralizing antibody titers and compact envelopes in early infection with subtype C human immunodeficiency virus type 1. J Virol 2006; 80:5211–5218.
14. Sagar M, Wu X, Lee S, Overbaugh J. Human immunodeficiency virus type 1 V1-V2 envelope loop sequences expand and add glycosylation sites over the course of infection, and these modifications affect antibody neutralization sensitivity. J Virol 2006; 80:9586–9598.
15. Derdeyn CA, Decker JM, Bibollet-Ruche F, Mokili JL, Muldoon M, Denham SA, et al
. Envelope-constrained neutralization-sensitive HIV-1 after heterosexual transmission. Science 2004; 303:2019–2022.
16. Frost SD, Liu Y, Pond SL, Chappey C, Wrin T, Petropoulos CJ, et al
. Characterization of human immunodeficiency virus type 1 (HIV-1) envelope variation and neutralizing antibody responses during transmission of HIV-1 subtype B. J Virol 2005; 79:6523–6527.
17. Iversen AK, Larsen AR, Jensen T, Fugger L, Balslev U, Wahl S, et al
. Distinct determinants of human immunodeficiency virus type 1 RNA and DNA loads in vaginal and cervical secretions. J Infect Dis 1998; 177:1214–1220.
18. Iversen AK, Fugger L, Eugen-Olsen J, Balslev U, Jensen T, Wahl S, et al
. Cervical human immunodeficiency virus type 1 shedding is associated with genital beta-chemokine secretion. J Infect Dis 1998; 178:1334–1342.
19. Iversen AK, Learn GH, Fugger L, Gerstoft J, Mullins JI, Skinhoj P. Presence of multiple HIV subtypes and a high frequency of subtype chimeric viruses in heterosexually infected women. J Acquir Immune Defic Syndr 1999; 22:325–332.
20. Adachi A, Gendelman HE, Koenig S, Folks T, Willey R, Rabson A, et al
. Production of acquired immunodeficiency syndrome-associated retrovirus in human and nonhuman cells transfected with an infectious molecular clone. J Virol 1986; 59:284–291.
21. Chesebro B, Wehrly K, Nishio J, Perryman S. Macrophage-tropic human immunodeficiency virus isolates from different patients exhibit unusual V3 envelope sequence homogeneity in comparison with T-cell-tropic isolates: definition of critical amino acids involved in cell tropism. J Virol 1992; 66:6547–6554.
22. Chesebro B, Nishio J, Perryman S, Cann A, O'Brien W, Chen IS, et al
. Identification of human immunodeficiency virus envelope gene sequences influencing viral entry into CD4-positive HeLa cells, T-leukemia cells, and macrophages. J Virol 1991; 65:5782–5789.
23. Toohey K, Wehrly K, Nishio J, Perryman S, Chesebro B. Human immunodeficiency virus envelope V1 and V2 regions influence replication efficiency in macrophages by affecting virus spread. Virology 1995; 213:70–79.
24. Platt EJ, Wehrly K, Kuhmann SE, Chesebro B, Kabat D. Effects of CCR5 and CD4 cell surface concentrations on infections by macrophagetropic isolates of human immunodeficiency virus type 1. J Virol 1998; 72:2855–2864.
25. Walter BL, Wehrly K, Swanstrom R, Platt E, Kabat D, Chesebro B. Role of low CD4 levels in the influence of human immunodeficiency virus type 1 envelope V1 and V2 regions on entry and spread in macrophages. J Virol 2005; 79:4828–4837.
26. Swofford DL. PAUP* 4.0: phylogenetic analysis using parsimony (*and other methods)
. Sunderland, MA; 1999.
27. Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 2003; 52:696–704.
28. Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007; 24:1596–1599.
29. Huang CC, Tang M, Zhang MY, Majeed S, Montabana E, Stanfield RL, et al
. Structure of a V3-containing HIV-1 gp120 core. Science 2005; 310:1025–1028.
30. Emsley P, Cowtan K. Coot: model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr 2004; 60:2126–2132.
31. Baker NA, Sept D, Joseph S, Holst MJ, McCammon JA. Electrostatics of nanosystems: application to microtubules and the ribosome. Proc Natl Acad Sci U S A 2001; 98:10037–10041.
32. Yang ZR, Thomson R, McNeil P, Esnouf RM. RONN: the bio-basis function neural network technique applied to the detection of natively disordered regions in proteins. Bioinformatics 2005; 21:3369–3376.
33. Iversen AK, Learn GH, Skinhoj P, Mullins JI, McMichael AJ, Rambaut A. Preferential detection of HIV subtype C' over subtype A in cervical cells from a dually infected woman. AIDS 2005; 19:990–993.
34. Berger EA, Murphy PM, Farber JM. Chemokine receptors as HIV-1 coreceptors: roles in viral entry, tropism, and disease. Annu Rev Immunol 1999; 17:657–700.
35. Cormier EG, Dragic T. The crown and stem of the V3 loop play distinct roles in human immunodeficiency virus type 1 envelope glycoprotein interactions with the CCR5 coreceptor. J Virol 2002; 76:8953–8957.
36. Nabatov AA, Pollakis G, Linnemann T, Kliphius A, Chalaby MI, Paxton WA. Intrapatient alterations in the human immunodeficiency virus type 1 gp120 V1V2 and V3 regions differentially modulate coreceptor usage, virus inhibition by CC/CXC chemokines, soluble CD4, and the b12 and 2G12 monoclonal antibodies. J Virol 2004; 78:524–530.
37. Huang CC, Lam SN, Acharya P, Tang M, Xiang SH, Hussan SS, et al
. Structures of the CCR5 N terminus and of a tyrosine-sulfated antibody with HIV-1 gp120 and CD4. Science 2007; 317:1930–1934.
38. Schulz GE. Molecular mechanism of biological recognition. New York: Elsevier/North-Holland Biomedical Press; 1979.
39. Janssens W, Buve A, Nkengasong JN. The puzzle of HIV-1 subtypes in Africa. AIDS 1997; 11:705–712.
40. Walker PR, Worobey M, Rambaut A, Holmes EC, Pybus OG. Epidemiology: sexual transmission of HIV in Africa. Nature 2003; 422:679.
41. Hladik F, Lentz G, Delpit E, McElroy A, McElrath MJ. Coexpression of CCR5 and IL-2 in human genital but not blood T cells: implications for the ontogeny of the CCR5+ Th1 phenotype. J Immunol 1999; 163:2306–2313.
42. Patterson BK, Landay A, Andersson J, Brown C, Behbahani H, Jiyamapa D, et al
. Repertoire of chemokine receptor expression in the female genital tract: implications for human immunodeficiency virus transmission. Am J Pathol 1998; 153:481–490.
43. Hladik F, Liu H, Speelmon E, Livingston-Rosanoff D, Wilson S, Sakchalathorn P, et al
. Combined effect of CCR5-delta32 heterozygosity and the CCR5 promoter polymorphism-2459 A/G on CCR5 expression and resistance to human immunodeficiency virus type 1 transmission. J Virol 2005; 79:11677–11684.
44. Thomas SM, Tse DB, Ketner DS, Rochford G, Meyer DA, Zade DD, et al
. CCR5 expression and duration of high risk sexual activity among HIV-seronegative men who have sex with men. AIDS 2006; 20:1879–1883.
45. Hoffman TL, MacGregor RR, Burger H, Mick R, Doms RW, Collman RG. CCR5 genotypes in sexually active couples discordant for human immunodeficiency virus type 1 infection status. J Infect Dis 1997; 176:1093–1096.
46. Laurent C, Bourgeois A, Faye MA, Mougnutou R, Seydi M, Gueye M, et al
. No difference in clinical progression between patients infected with the predominant human immunodeficiency virus type 1 circulating recombinant form (CRF) 02_AG strain and patients not infected with CRF02_AG, in western and west-central Africa: a four-year prospective multicenter study. J Infect Dis 2002; 186:486–492.
47. Westby M, Smith-Burchnell C, Mori J, Lewis M, Mosley M, Stockdale M, et al
. Reduced maximal inhibition in phenotypic susceptibility assays indicates that viral strains resistant to the CCR5 antagonist maraviroc utilize inhibitor-bound receptor for entry. J Virol 2007; 81:2359–2371.
48. Kuhmann SE, Pugach P, Kunstman KJ, Taylor J, Stanfield RL, Snyder A, et al
. Genetic and phenotypic analyses of human immunodeficiency virus type 1 escape from a small-molecule CCR5 inhibitor. J Virol 2004; 78:2790–2807.
49. McClure CP, Tighe PJ, Robins RA, Bansal D, Bowman CA, Kingston M, et al
. HIV coreceptor and chemokine ligand gene expression in the male urethra and female cervix. AIDS 2005; 19:1257–1265.
50. Hladik F, Sakchalathorn P, Ballweber L, Lentz G, Fialkow M, Eschenbach D, et al
. Initial events in establishing vaginal entry and infection by human immunodeficiency virus type-1. Immunity 2007; 26:257–270.
51. Sourisseau M, Sol-Foulon N, Porrot F, Blanchet F, Schwartz O. Inefficient human immunodeficiency virus replication in mobile lymphocytes. J Virol 2007; 81:1000–1012.
52. Prakash M, Patterson S, Gotch F, Kapembwa MS. Ex vivo analysis of HIV-1 co-receptors at the endocervical mucosa of women using oral contraceptives. BJOG 2004; 111:1468–1470.
53. Wain LV, Bailes E, Bibollet-Ruche F, Decker JM, Keele BF, Van Heuverswyn F, et al
. Adaptation of HIV-1 to its human host. Mol Biol Evol 2007; 24:1853–1860.
54. Lockett SF, Alonso A, Wyld R, Martin MP, Robertson JR, Gore SM, et al
. Effect of chemokine receptor mutations on heterosexual human immunodeficiency virus transmission. J Infect Dis 1999; 180:614–621.
55. Centlivre M, Sala M, Wain-Hobson S, Berkhout B. In HIV-1 pathogenesis the dice is cast during primary infection. AIDS 2007; 21:1–11.
Keywords:© 2009 Lippincott Williams & Wilkins, Inc.
CCR5; cervix; evolution; HIV-1; phylogeny; subtype A; subtype C; transmission