Skip to main content

Application of a case–control study design to investigate genotypic signatures of HIV-1 transmission

Abstract

Background

The characterization of HIV-1 transmission strains may inform the design of an effective vaccine. Shorter variable loops with fewer predicted glycosites have been suggested as signatures enriched in envelope sequences derived during acute HIV-1 infection. Specifically, a transmission-linked lack of glycosites within the V1 and V2 loops of gp120 provides greater access to an α4β7 binding motif, which promotes the establishment of infection. Also, a histidine at position 12 in the leader sequence of Env has been described as a transmission signature that is selected against during chronic infection. The purpose of this study is to measure the association of the presence of an α4β7 binding motif, the number of N-linked glycosites, the length of the variable loops, and the prevalence of histidine at position 12 with HIV-1 transmission. A case–control study design was used to measure the prevalence of these variables between subtype B and C transmission sequences and frequency-matched randomly-selected sequences derived from chronically infected controls.

Results

Subtype B transmission strains had shorter V3 regions than chronic strains (p = 0.031); subtype C transmission strains had shorter V1 loops than chronic strains (p = 0.047); subtype B transmission strains had more V3 loop glycosites (p = 0.024) than chronic strains. Further investigation showed that these statistically significant results were unlikely to be biologically meaningful. Also, there was no difference observed in the prevalence of a histidine at position 12 among transmission strains and controls of either subtype.

Conclusions

Although a genetic bottleneck is observed after HIV-1 transmission, our results indicate that summary characteristics of Env hypothesised to be important in transmission are not divergent between transmission and chronic strains of either subtype. The success of a transmission strain to initiate infection may be a random event from the divergent pool of donor viral sequences. The characteristics explored through this study are important, but may not function as genotypic signatures of transmission as previously described.

Background

Multiple research groups have invested substantial effort into defining signatures within the envelope protein (Env) of HIV-1 that are highly predictive of viral transmission [113]. A genetic bottleneck occurs at HIV-1 transmission, and early replicating viral populations remain homogenous [3], so that investigating the transmitted strains across many individuals may reveal a transmission signature and point to a process important in establishing infection [1, 2, 47]. These studies not only provide insight into the biological processes occurring during mucosal transmission, they may also inform the development of an effective vaccine against HIV-1. The discovery of a transmission signature could reveal motifs that may be successfully targeted during vaccination which could lead to a more protective response compared to vaccines developed from chronically replicating viral sequences [1417].

A number of different approaches have been used to identify transmission signatures, investigating whether shorter variable loops with fewer predicted N-linked glycosylation sites (glycosites) are characteristics enriched in Env from viruses derived during acute infection. Comparing sequences from chronically infected donors and early sequences of infected recipients is a more powerful way to study transmission, although sample sizes often remain small. Using eight donor-recipient pairs, Derdeyn et al. showed significant differences in subtype C Env sequences, where early V1-V4 loops were shorter (p = 0.02) and had fewer glycosites (p = 0.037) [4]. Also using donor-recipient pairs, Frost et al. found no difference between envelope length or the number of glycosites in sequences derived from eight men who have sex with men during acute infection and their transmission pairs [18]. However, in the study Frost et al. state that only half of the donors had confirmed chronic infection, and that four were recently infected donor subjects [18], which could bias the results. In ten self-reported transmission pairs of subtype B, Liu et al. showed that gp120 sequences had shorter variable regions and reduced numbers of glycosites following transmission to the new host [19]. However, they state that their results were dependent upon how pairs were chosen for inclusion, and that no other study had yet found these reductions to be measured consistently across various acute infections. Chohan et al. compared 35 early and 51 late sequences of subtype A and 13 early and 82 chronic sequences of subtype B (chronic sequences were selected from the Los Alamos HIV Sequence Database (LANL) [20]). They showed early sequences had shorter loops (p = 0.008) and fewer glycosites (p = 0.017) for subtype A, but not subtype B (p = 0.340; p = 0.640, respectively) [5].

Other studies have suggested transmission-linked glycosites in the V1 and V2 loops surrounding the α4β7 binding site are significant in terms of transmission [11, 15, 21, 22]. The V2 loop of gp120 contains an α4β7 binding site, defined as LDV/I at position 180 (according to the HXB2 reference sequence) [11, 2123]. α4β7 is a gut homing receptor colocalized with CCR5 and CD4 on CD4+ T cells, and projects out further from the cell surface than CCR5 or CD4 [11, 2124]. α4β7 +/CD4+ T cells infected with HIV-1 will migrate to the gut, where HIV-1 preferentially replicates and causes irrevocable damage to the immune system [11, 2124]. Although not necessary for HIV-1 entry into CD4+ T cells, Env lacking specific glycosites within the V1 and V2 loops of gp120 have greater access to bind α4β7, an interaction that has been shown to promote a greater susceptibility to infection [11, 15, 22]. This increased capacity for infection may be an important transmission signature [11, 15]. However, no statistical analyses have been used to measure the association of the absence of glycans in close proximity to the α4β7 binding site and transmission.

Recently, models scanning the entire Env sequence for patterns in amino acids have been used to search for transmission signatures in subtype B Env sequences [1, 2]. Gnanakaran et al., through a retrospective hypothesis-raising ‘test’ analysis, found that a histidine (His) at position 12 in the Env leader sequence was associated with transmission because it was selected against in chronic infection [1, 2]. Asmal et al. supported this phenotypically, showing that Env sequences containing a positively charged amino acid at position 12 were associated with higher levels of expression due to more efficient trafficking of the nascent polypeptide to the endoplasmic reticulum [2].

Here, we applied an epidemiological approach to investigating the transmission signature, using a case–control study design to measure associations between HIV-1 transmission strains (‘cases’) and chronic infection strains (‘controls’). Case–control studies are often used in epidemiology and measure exposure variables between cases with the outcome of interest, in this study, a transmission event, and randomly selected population controls [25, 26]. The transmission strains used in this study were consensus sequences (obtained through backwards phylogeny as the most recent common ancestor of acute infection) from patients infected with one viral sequence of either subtype B or C [6, 7].

Exposure variables (biological characteristics hypothesized to be associated with transmission) explored in this study were the amino acid length and number of glycosites (all PNGs are shown to be glycosylated in gp120 [27, 28]) of each variable loop. Where previous studies examined the loops together [4, 5], we chose to measure the hypervariable region [27] of each variable loop (see methods for loop definitions) independently to see if the biological characteristics in one loop had a significant affect where another loop did not. The numbers of glycosites in the conserved regions of Env were also explored. Other variables measured included the prevalence of the α4β7 binding site, as well as the His at position 12 in the leader sequence of Env, between transmission and chronic strains. In this study, significant associations from both a biological and statistical perspective between variables and HIV-1 transmission would be considered evidence of transmission signatures. Rather than looking for new potential transmission signatures, the aim of this study was to measure associations between previously suggested transmission signatures and their relationship with a successful transmission event.

Results

Amino acid lengths of the variable loops

When measuring the amino acid (AA) lengths, only subtype C V1 loops and subtype B V3 loops were significantly different between transmission strains and chronic controls (p = 0.047; p = 0.031, respectively; Table 1). Subtype C V1 loops were shorter in transmission strains, with a median difference of 2 AA; and the median difference between subtype B V3 lengths was zero, as demonstrated in the box plots (Figure 1). In additional analyses based on these significant findings, and because the sample size was large enough, the mean differences in amino acid length were explored. Subtype C V1 loops from transmitted strains were, on average, 3.2 AA shorter than controls; and the mean difference between V3 loops of subtype B transmission strains and controls was 0.14 AA (Additional file 1: Table S1). There were no other significant differences in AA length for the other loops or in either subtype (Figure 1).

Table 1 Comparison of median amino acid lengths between transmission strains and controls
Figure 1
figure 1

Box plots displaying amino acid length of each loop between chronic sequences and transmission strains.

Glycosites within the variable loops and conserved regions

In a logistic regression analysis exploring predictors of whether or not a strain was a transmission case or chronic control, no significant results were observed when measuring the association of the number of glycosites between cases and controls in the univariate analyses for the variable loops (Table 2) or conserved regions (Table 3). In models allowing for differences by subtype (using a statistical interaction term), odds ratios for the number of glycosites changed by over 10% in the V1 and V4 loops, but not in the V2 loop (Table 2). In the V4 loop, interaction was further confirmed by the test for homogeneity (p = 0.001; Table 2). This was not observed in the V1 or V2 loops. In the V3 loop, logistic regression could not be performed since all of subtype B transmission strains contained the one glycosite, as did all subtype C controls. Overall, 98.5% of transmission strains and 96.2% of controls contained this glycosite in the V3 loop (Table 4). The Wilcoxon rank sum test showed a significant difference (p = 0.024) between cases and controls when considering subtype B sequences with 1 versus 0 glycosites (Table 2), but the median difference of the number of glycosites in the V3 loop between transmission strains and controls was zero. There were no other significant findings in the variable loops when considering the number of glycosites, a predictor of transmission. Furthermore, there were no significant findings in the conserved regions (Table 3).

Table 2 Logistic Regression for predicting transmission based on number of glycosites in the variable loops
Table 3 Logistic Regression for predicting transmission based on number of conserved glycosites
Table 4 Baseline characteristics of cases (transmission strains) and controls (chronic strains)

Exposure of the α4β7 binding site

To predict the capacity of gp120 to bind α4β7, we first determined the prevalence of this binding site in the viral sequences. Of the 133 transmission strains, 69.2% contained the LDV/I α4β7 binding site, and 100% contained the aspartic acid (D). Of the 133 controls, 70.1% contained the LDV/I, and 98.5% contained the aspartic acid at this site (Table 4). As described above, there was no difference in the number of glycosites between transmission strains and controls, and thus no difference in exposure for the binding interaction between gp120 and the α4β7 receptor on CD4+ T cells.

Histidine at position 12

As a His at position 12 has been suggested as a transmission signature [1, 2], we measured the prevalence of this motif among transmission strains and controls. 75.6% of subtype B transmission strains contained a His at position 12, as did 71.8% of the controls Additional file 1:Table S1. Using logistic regression, the odds ratio changed by over 10% from the univariate to interaction model and so subtype-specific odds ratios were reported (Table 5). There were no significant differences when measuring if a His at position 12 was a predictor of transmission. Further, the p-values from the Wilcoxon rank sum test comparing sequences with and without the His at position 12 showed no difference between transmission strains and controls for either subtype (Table 5).

Table 5 Logistic Regression for predicting transmission based on His at position 12.

Discussion

This study found that, with one exception, there were no biologically meaningful differences in amino acid length or the number of glycosites in the variable loops or conserved regions, or the prevalence of His at position 12, between transmission strains and controls in either subtype. The exception was that the V1 loops of subtype C transmission strains were shorter than controls, based on a median difference of 2, and a mean difference of 3.2 AA. Our results, derived from a case–control study design with the largest currently available data set (133 “founder” sequences and matched controls) indicate that if there is indeed a “transmission signature” for HIV-1 then it is not manifest in the previously hypothesized summary measures that we have investigated.

As in all case–control studies, we cannot exclude potential confounding effects due to unrecorded factors such as route of transmission, coinfection with other STIs, or phylogenetic relationships between samples [25, 26]. Another limitation is that during the random control selection, it is possible some transmission strains may have been selected. Such selections, expected to be exceedingly unlikely due to the diversity of circulating viral sequences within an individual, would bias our results towards the null hypotheses of no difference between transmission and chronic strains. Another potential limitation of this study is that the sequences used for cases and controls have not been validated for functionality. This is an issue that is inherent in studies examining large numbers of sequences at a population level. Other information that would greatly benefit investigations of this type include the sequence variation of viruses in the transmitting fluids of the infected donor, and how variation from this quasispecies decreases in the new host at the mucosa and through early selection pressure to persist beyond this genetic bottleneck. Unfortunately this information is currently unobtainable.

Patients that displayed multivariant transmission were excluded from this study because it is currently unknown what factors lead to multiple viruses establishing infection. Rieder et al. investigated envelope sequences from 145 patients during acute infection. They found that concomitant STI, gender, and sexual practice were not associated with transmission of heterogeneous virus populations, suggesting that transmission of multiple HIV-1 variants is not dependent on mucosal factors, but a more complex combination of factors that have yet to be determined [29].

Current research shows that having a His at position 12 in the Env leader sequence is enriched in transmission strains [1, 2]. However, our study showed there was no difference in the prevalence of His between transmission strains and controls of subtype B or C. Asmal et al. also indicated that His is the most common AA at position 12 in subtype B and most other subtypes [2]. However, our study revealed that the overall prevalence of His at position 12 in subtype C is small (4.5%; Additional file 2: Table S2), whilst glutamine, an AA without charge, is most common at position 12 among transmission strains and controls of subtype C (85.5% and 76.4%, respectively).

Gnanakaran et al. showed that consensus sequences of Subtype B transmission strains contained a His at position 12 more often than consensus sequences of chronic strains (74% versus 57%) [1]. Of note, in this study, there was no statistical difference between transmission strains and controls of subtype B (72% versus 76%; Additional file 2: Table S2). The differences between the prevalence of His at position 12 in these chronic populations could be due to control selection where we selected controls from the Los Alamos database (n = 78) and Gnanakaran et al. used intra-patient consensus sequences that had all been sequenced using single genome amplification (n = 43) [1]. However, in their study, Gnanakaran et al. also used another set of controls from chronically infected patients selected from the Los Alamos database to increase their statistical power [1]. After investigation of this control population, we found that 86% of these chronic sequences contained a His at position 12, supporting our findings.

The various methods used in this study allow different statistical parameters to be explored. This is advantageous when searching for biological meaning behind statistical significance. The Wilcoxon rank sum test concludes a significant difference in subtype B V3 AA length by determining the probability that differences in median length occurred by chance between transmission strains and controls from chronic infection. However, further investigation reveals that the median difference is zero. This is also reflected in the t-test (where we probe mean, rather than median) where we detected that on average, the V3 loops of transmission strains are just 0.14 AA shorter. This value does not represent a full AA, and therefore the means of the groups are biologically similar, if not indistinguishable. This deeper investigation enables us to interpret biological meaning, rather than just presenting significant p-values from the Wilcoxon rank sum test. In subtype C, the V1 loops of transmission strains are on average, 3.2 AA acids shorter than chronic controls. While this is a more relevant number consisting of greater than one residue, phenotypic studies must consider if three residues contributes to a functional difference in envelope behavior that would be important during transmission. There were no biologically relevant differences in the number of glycosites between transmission strains and chronic controls for either subtype.

Conclusions

When gp120 has access to bind α4β7, there exists a capacity for a population of CD4+ T cells to be highly susceptible to infection with HIV-1 [11, 15, 2123]. However, no differences in the number of glycosites in the V1 and V2 loops were observed between transmission strains and controls of either subtype. These results suggest that although sequences lacking these glycosites can bind α4β7, it may not be an event strongly associated with transmission. Likewise, the His at position 12 in the leader sequence does impact envelope expression and thus virion infectivity [1, 2]. However, the results from this study suggest the prevalence of His is not different between transmission strains and controls.

Although important, the characteristics explored through this study may not function as genotypic signatures characteristic of HIV-1 during transmission, and may not differ throughout pathogenesis. Although a genetic bottleneck is observed after transmission, the selection of the transmitting strain from the divergent pool of donor viral sequences may be a random event. The “random transmission” of HIV-1 has been hypothesized by Hedskog et al., who provide a compelling argument that in the context of CCR5, CXCR4, and R5X4 coreceptor use at transmission, the viral characteristics that dominate the viral population of the donor will be transmitted [30].

Because of discordant findings across research studies, further genotypic, phenotypic, and structural studies with larger sample sizes are warranted and may require the exploration of the subtypes independently for vaccine design.

Methods

To investigate genotypic signatures of HIV-1 transmission, a case–control study design was used to measure the association of AA lengths and the number of glycosites in each of the variable loops of gp120, as well as the conserved regions. The prevalence of the α4β7 binding site (LDV/I) was also investigated due to the transmission-linked associations of glycosites surrounding this site in the V2 loop of gp120 [11, 2123]. In addition, the prevalence of a His at position 12 was determined, and the association of this residue as a transmission signature was measured. The total sample size of 266 sequences was powered at 80% to reject the null hypothesis (of no difference) at 95% precision.

133 transmission strain cases (78 subtype B, 55 subtype C) were acquired from Keele et al. and Abrahams et al.[6, 7]. These envelope sequences, derived during acute infection, were obtained from individuals infected with only one transmitting virus, as determined using a mathematical model of backwards phylogeny [6, 7]. Keele et al. and Abrahams et al. were able to define the most recent common ancestor of infection (the consensus founder sequence) for each case and it is these sequences that are used as the transmission strains in this study. To date, only Keele et al. and Abrahams et al. have defined transmission strains in this way. Therefore, we have used the entire founder “case population” currently available for this type of study.

133 control sequences were derived from the plasma of individuals with chronic infection, acquired from LANL [20]. A random numbers table was generated using the Stat Trek Random Number Generator http://stattrek.com to guarantee the random selection of controls. Each selected control sequence was investigated and submitted to rigorous exclusion criteria, contributing to the strong internal validity of this study (Figure 2). For example, if a randomly selected sequence happened to come from a patient whom we had already derived a control sequence from, that newly selected sequence was discarded and another random selection from the LANL database made. This process ensured the statistical independence of our control samples and that all sequences from a particular patient in the LANL database had the same chance of being selected as a control sequence. To adjust for confounding variables, controls were frequency-matched on HIV-1 subtype and geographical location. However, only the subtype designation was considered in the analyses due to the collinearity between location and subtype. 78 subtype B chronic infection controls were selected from the USA/Trinidad and Tobago and 55 subtype C sequences were selected from South Africa/Malawi. Demographic details including behaviour status were not available for most control sequences and so were not controlled for in the analyses.

Figure 2
figure 2

Case and Control Selection.

Case and control DNA sequences, as well as the HXB2 reference envelope sequence [27], were converted to AA sequences and then aligned using a progressive multiple alignment method (multialign, MATLAB 2010b, The MathWorks Inc., Natick MA, USA). AA positions in the aligned sequences were numbered relative to HXB2 according to the convention of Korber et al.[31]. N-glycosylation sites were then determined on the original sequence for each variable loop identified through the alignment. LDV/I residues and position 12 Histidine residues were counted in each of the 266 sequences. AA lengths and glycosites were counted separately for each of the variable loops (V1-V4) and conserved regions. In this study, the hypervariable regions of the gp120 variable loops were used and defined as Env137-151 for the V1 loop, Env161-195 for the V2 loop, Env300-328 for the V3 loop, and Env393-414 for the V4 loop according to the HXB2 sequence [27].

The Wilcoxon rank sum test was used to determine differences in the AA lengths among the variable loops between cases and controls, analysed separately for subtype B and subtype C. The null hypothesis of the rank sum test was that when values are ranked in ascending order, there is no difference in median values between both groups. Box plots were created to present the observations. Where results displayed p-values less than or equal to 0.05, t-tests were also performed to explore the mean difference.

Logistic regression was used to assess the association of the number of glycosites on transmission. The number of glycosites was considered a categorical variable, based on the distribution and median values of the number of glycosites in each variable loop and each conserved region among the sequences (Figure 3; characteristics for conserved regions not shown). For example, the V1 and V4 loops have the largest number of glycosites present (1 to 8 and 0 to 4, respectively); the median value is an independent category, with categories above and below the median value. The range in the V2 loop (0 to 3) is smaller. The median value is included in the lower category and values above the median are in a higher category. Both cases and controls have either 0 or 1 glycosite in the V3 loop.

Figure 3
figure 3

Histograms of baseline frequencies of the number of glycosites in each variable loop (for subtype-specific characteristics, see Additional file 2 : Table S2).

Because subtypes are genetically distinct, subtype specific odds ratios were measured in a model that tested for interaction between the subtypes. When compared to the univariate models, a change in odds ratio greater than 10% for either subtype was used as a condition to include the interaction and report both univariate and subtype-specific odds ratios (Table 2). The null hypotheses were that the number of glycosites is not associated with HIV-1 transmission in any loop or conserved region, and that subtype does not modify this effect (interaction).

Logistic regression was used to assess the relationship between the His at position 12 and transmission. The null hypothesis with logistic regression was that a sequence containing a His at position 12 is not associated with HIV-1 transmission, and that subtype does not affect this association.

For both logistic models, Wilcoxon rank sum tests were used to test for differences in the number of glycosites and the presence of His at position 12 for comparison with previously published work. All analyses were performed using Stata version 10.1 (Stata Corp, TX, USA).

Author’s contributions

TM designed the study, collected the data, performed the statistical analyses, interpreted the data and prepared the manuscript. JMM aligned the sequences and participated in data collection and provided helpful discussion regarding data analysis. RJC provided meaningful discussion and helped to interpret the data. DFJP conceived of the virological aspects of the study and provided helpful discussion. JMcC oversaw and contributed to the design of the study and helped to draft the manuscript. All authors read and approved the final manuscript. Both corresponding authors contributed equally to this paper.

Accession Numbers

Control sequences: all from LANL [AY352275, AY357342, EF593188, DQ853447, DQ410638, DQ410615, DQ410601, DQ410579, DQ410566, DQ410557, DQ410531, DQ410522, DQ410509, DQ410474, DQ410422, DQ410400, DQ410361, DQ410321, DQ410275, DQ410237, DQ410224, DQ410214, DQ410170, DQ410119, GU728030, DQ410096, DQ410083, DQ410058, AY835781, AY835775, AY835771, AY835763, AY835759, AY842843, AY842821, AY842807, U69584, AF025763, AF025760, AF025750, U26546, FJ469771, FJ469768, FJ469755, FJ469749, FJ469746, FJ469735, FJ469732, FJ469712, FJ469700, FJ653616, FJ653586, FJ653567, FJ653545, FJ653486, FJ653473, FJ653427, FJ653397, FJ653378, FJ653309, FJ653240, FJ653208, FJ653135, FJ653113, FJ798580, FJ798423, FJ798397, GU938295, GU727972, GU728371, EU604604, GU728337, EU604554, GU728368, GU728375, GU728244, GU728195, GU728369, AB485644, AF110967, AF110977, AF290027, AF391244, AF391246, AF443087, AF443097, AF443099, AY463234, AY529660, AY772691, AY878056, AY901970, AY901972, DQ093591, DQ164122, DQ275658, DQ275660, DQ351224, DQ351233, DQ369982, DQ369984, DQ382365, DQ382379, DQ396367, DQ396377, DQ396399, DQ396399, DQ978981, FJ846650, GU080179, GU216703, GU216785, GU216805, GU216824, GU216838, GU329068, GU329085, GU329098, GU329118, GU329147, GU329160, GU329169, GU329191, GU329208, GU329221, GU329241, GU329261, GU329273, GU329294, GU329314, GU329377, GU329389, GU329409]

Transmission strain consensus sequences: Obtained through personal correspondence with Brandon F. Keele and Melissa-Rose Abrahams.

Abbreviations

PNGs:

predicted N-linked glycosylation sites; glycosites

LANL:

Los Alamos National Laboratory

HIV:

Sequence Database

STI:

sexually transmitted infection.

References

  1. Gnanakaran S, Bharracharya T, Daniels M, Keele BF, Hraber PT, Lapedes AS, Shen T, Gaschen B, Krishnamoorthy M, Li H, Decker JM, Salazar-Gonzalez JF, Wang S, Jiang C, Gao F, Swanstrom F, Anderson J, Ping L-H, Cohen MS, Markowitz M, Goepfert PA, Saag MS, Eron JJ, Hicks CB, Blattner WA, Tomaras GD, Asmal M, Letvin NL, Gilbert PB, DeCamp AC, Magaret CA, Schief WR, Ban Y-E A, Zhang M, Soderberg KA, Sodroski JG, Haynew BF, Shaw GM, Hahn BH, Korber B: Recurrent signature patterns in HIV-1 B clade envelope glycoproteins associated with either early or chronic infections. PLoS Pathog. 2011, 7 (9): e1002209-10.1371/journal.ppat.1002209.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Asmal M, Hellmann I, Lie W, Keele BF, Perelson AS, Bhattacharya T, Gnanakaran S, Daniels M, Haynes BF, Korber BT, Haha BH, Shaw GM, Letvin NL: A Signature in HIV-1 envelope leader peptide associated with transition from acute to chronic infection impacts envelope processing and infectivity. PLoS ONE. 2011, 6: e23673-10.1371/journal.pone.0023673.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Haaland RE, Hawkins PA, Salazar-Gonzalez J, Johnson A, Tichacek A, Karita E, Manigart O, Mulenga J, Keele BF, Shaw GM, Hahn BH, Allen S, Derdeyn CA, Hunter E: Inflammatory genital infections mitigate a severe genetic bottleneck in heterosexual transmission of subtype A and C HIV-1. PLoS Pathog. 2009, 5: e1000274-10.1371/journal.ppat.1000274.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Derdeyn CA, Kecker JM, Bibollet-Ruche F, Mokili JL, Muldoon M, Denham SA, Heil ML, Kasolo F, Musonda R, Hahn B, Shaw G, Korber BT, Allen S, Hunter E: Envelope-constrained neutralization-sensitive HIV-1 after heterosexual transmission. Science. 2004, 303: 2019-2022. 10.1126/science.1093137.

    Article  CAS  PubMed  Google Scholar 

  5. Chohan B, Lang D, Sagar M, Korber B, Lavreys L, Richardson B, Overbaugh J: Selection for human immunodeficiency virus type 1 evnelope glycosylation variants with shorter V1-V2 loop sequences occurs during transmission of certain genetic subtypes and may impact viral RNA levels. Virology. 2005, 79 (10): 6528-6531. 10.1128/JVI.79.10.6528-6531.2005. Slides with audio can be accessed: http://pag.aids2010.org/flash/?pid=100868.

    Article  CAS  Google Scholar 

  6. Keele BF, Giorgi EE, Salazar-Gonzalez JF, Decker JM, Pham KT, Salazar MG, Sun C, Grayson T, Wang S, Li H, Wei X, Jiang C, Kirchherr JL, Gao F, Anderson JA, Ping L-H, Swanstrom R, Tomaras GD, Blattner WA, Goepfert PA, Kilby JM, Saag MS, Delwart EL, Busch MP, Cohen MS, Montefiori DC, Haynes BF, Gaschen B, Athreya GS, Lee HY, Woo N, Seoighe C, Perelson AS, Bhattacharya T, Korber BT, Hahn BH, Shaw GM: Identification and characterization of transmitted and early founder virus envelopes in primary HIV-1 infection. Proc Natl Acad Sci. 2008, 105: 7552-7557. 10.1073/pnas.0802203105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Abrahams MR, Anderson JA, Giorgi EE, Seoighe C, Mlisana K, Ping L-H, Athreya GS, Treurnicht FK, Keele BF, Wood N, Salazar-Gonzalez JF, Bhattacharya T, Chu H, Hoffman I, Galvin S, Mapanje C, Kazembe P, Thebus R, Fiscus S, Hide W, Cohen MS, Abdool Kaim S, Haynes BF, Shaw GM, Hahn BH, Korber BT, Swanstrom R, Williamson C: Quantitating the multiplicity of infection with human immunodeficiency virus type 1 subtype C reveals a non-poisson distriution of transmitted variants. J Virol. 2009, 83 (8): 3556-10.1128/JVI.02132-08.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Salazar-Gonzalez JF, Salazar MG, Keele BF, Learn GH, Giorgi EE, Li H, Decker JM, Wang S, Baalwa J, Kraus MH, Parrish NF, Ochsenbauer-Jambor C, Kappes JC, Saag MS, Cohen MS, Mulenga J, Derdeyn CA, Allen S, Hunter E, Markowitz M, Hraber R, Perelson AS, Bhattacharya T, Haynes BF, Korber BT, Hahn BH, Shaw GM: Genetic identity, biological phenotype, and evolutionary pathways of transmitted/founder viruses in acute and early HIV-1 infection. J Exp Med. 2009, 206 (6): 1273-1289. 10.1084/jem.20090378.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Sagar M, Laeyendecker O, Lee S, Gamiel J, Wawr MJ, Gray RH, Serwadda D, Sewankambo NK, Shepherd JC, Toma J, Huang W, Quinn T: Selection of HIV variants with signature genotypic characteristics during hetersexual transmission. JID. 2009, 199: 580-589. 10.1086/596557.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Sagar M, Kirkegaard E, Lavreys L, Overbaugh J: Diversity in HIV-1 Envelope V1-V3 sequences early in infection reflects sequences diversity throughout the HIV-1 genome but down not predict the extent of sequence diversity during chronic infection. AIDS Research and Human Retroviruses. 2006, 22 (5): 430-437. 10.1089/aid.2006.22.430.

    Article  CAS  PubMed  Google Scholar 

  11. Nawaz F, Cicala C, Van Ryk D, Block KE, Jelicic K, McNally JP, Ogundare O, Pascuccio M, Patel N, Wie D, Fauci AS, Arthos J: The genotype of early-transmitting HIV gp120s promotes α4β7-reactivity, revealing α4β7+/CD4+ T cells as key targets in mucosal transmission. PLoS Pathog. 2011, 7 (2): e1001301-10.1371/journal.ppat.1001301.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Ritola K, Pilcher C, Fiscus S, Hoffman NG, Nelson JAE, Kitrinos KM, Hicks CB, Eron JJ, Swanstorm R: Multiple V1/V2 env variants re frequently present during primary infection with human immunodeficiency virus type 1. Virology. 2004, 78 (20): 11208-11218. 10.1128/JVI.78.20.11208-11218.2004.

    Article  CAS  Google Scholar 

  13. Wood N, Bhattacharya T, Keele BF, Giorgi E, Lie M, Gaschen B, Daniels M, Ferrari G, Haynes BF, McMichael A, Shaw G, Hahn BA, Korber B, Seoighe C: HIV evolution in early infection: selection pressures, patterns of insertion and deletion, and the impact of APOBEC. PLoS Pathogens. 2009, 5 (5): e1000414-10.1371/journal.ppat.1000414.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Keele BF, Derdeyn C: Genetic and antigenic features of the transmitted virus. Curr Opin HIV AIDS. 2009, 4: 352-357. 10.1097/COH.0b013e32832d9fef.

    Article  PubMed  Google Scholar 

  15. Anthony Fauci: New Concepts in HIV/AIDS Pathogenesis: Implications for Interventions. XVIII International AIDS Conference. 2010, Vienna, Slides with audio can be accessed:http://pag.aids2010.org/flash/?pid=100868

    Google Scholar 

  16. Derdeyn CA, Hunter E: Viral characteristics of transmitted HIV. Curr Opin HIV AIDS. 2009, 3: 16-21.

    Article  Google Scholar 

  17. Keele BF: Identifying and characterizing recently transmitted viruses. Curr Opin HIV AIDS. 2010, 5: 327-334. 10.1097/COH.0b013e32833a0b9b.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Frost SDW, Liu Y, Kosakovsky Pond SL, Chappey C, Wrin T, Petropoulos CJ, Little SJ, Richman DD: J Virol. 2005, 79 (10): 6523-10.1128/JVI.79.10.6523-6527.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Liu Y, Curlin ME, Diem K, Zhao H, Ghosh AK, Zhu H, Woodward AS, Maenza J, Stevens CE, Stekler J, Collier AC, Genowati I, Deng W, Zioni R, Corey L, Zhu T, Mullins JI: Env length and N-linked glycosylation following transmission of human immunodeficiency virus Type 1 subtype B viruses. Virology. 2008, 347: 229-233.

    Article  Google Scholar 

  20. Los Alamos National Laboratories. Los Alamos HIV Sequence Database. http://www.hiv.lanl.gov/ Operated by Los Alamos National Security, LLC, for the U.S. Department of Energy's National Nuclear Security Administration

  21. Cicala C, Martinelli E, McNally J, Goode DJ, Gopaul R, Hiatt J, Jelicic K, Kottilil S, Macleod K, O’Shea A, Patel N, Van Ryk D, Wie D, Pascuccio M, Ling Y, McKinnon L, Izulla P, Kimani J, Kaul R, Fauci AS, Arthos J: The integrin α4β7 forms a complex with cell-surface CD4 and defines a T-cell subset that is highly susceptible to infection by HIV-1. December. 2009, 106 (49): 20877-20882.

    CAS  Google Scholar 

  22. Cicala C, Arthos J, Fauci A: HIV-1 envelope, integrins and co-receptor use in mucosal transmission of HIV. J Transl Med. 2011, 9 (Suppl 1): S2-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Arthos J, Cicala C, Martinelli E, Macleod K, Van Ryk D, Wie D, Xiao Z, Veenstra TD, Conrad TP, Lempicki RA, McLaughlin S, Pascuccio M, Gopaul R, McNally J, Cruz CC, Censoplano N, Chung E, Reitano KN, Kottilil S, Goode DJ, Fauci AS: HIV-1 envelope protein binds to and signals through integrin α4β7, the gut mucosal homing receptor for peripheral T cells. Nat Immunol. 2008, 9 (3): 301-309. 10.1038/ni1566.

    Article  CAS  PubMed  Google Scholar 

  24. von Andrian UH: CR Mackay. T-cell function and migration. Two sides of the same coin. N Engl J Med. 2000, 343: 1020-1034. 10.1056/NEJM200010053431407.

    Article  CAS  PubMed  Google Scholar 

  25. Webb P, Bain C, Pirozzo S: Essential Epidemiology: and introduction for Students and Health Professionals. 2005, Cambridge University Press, New York

    Google Scholar 

  26. Kirkwood BR, Sterne JAC: Medical Statistics. 2003, Blackwell Science, Inc, Oxford, 2

    Google Scholar 

  27. Leonard CK, Spellman MW, Riddle L, Harris RJ, Thomas JN, Gregory TJ, et al: Assignment of intrachain disulfide bonds and characterization of potential glycosylation sites of the type 1 recombinant human immunodeficiency virus envelope glycoprotein (gp120) expressed in chinese hamster ovary cells. J Biological Chemistry. 1990, 265 (18): 10373-10382.

    CAS  Google Scholar 

  28. Binley JM, Ban Y-E A, Crooks ET, Eggink D, Osawa K, Schief WR, Sanders RW: Role of complex carbohydrates in human immunodeficiency virus type 1 infection and resistance to antibody neutralization. J of Virology. 2010, 84 (11): 5637-5655. 10.1128/JVI.00105-10.

    Article  CAS  Google Scholar 

  29. Rieder P, Joos B, Scherrer AU, Kuster H, Braun D, Grube C, Niederost B, Leemann C, Gianella S, Metzner KJ, Boni J, Weber R, Gunthard HF: Characterization of human immunodficiency virus type 1 (HIV-1) diversity and tropism in 145 patients with primary HIV-1 infection. Clin Infect Dis. 2011, 53 (12): 1271-9. 10.1093/cid/cir725.

    Article  CAS  PubMed  Google Scholar 

  30. Hedskog C, Mild M, Albert J: Transmission of the X4 Phenotype of HIV-1: Is there evidence against the “Random Transmission” hypothesis?. J Infect Dis. 2012, 205: 163-5. 10.1093/infdis/jir719.

    Article  PubMed  Google Scholar 

  31. Korber BT, Foley BT, Kuiken CL, Pillai SK, Sodroski JG: Numbering Positions in HIV Relative to HXB2CG. Human Retroviruses and AIDS. Edited by: Korber B, Kuiken CL, Foley B, Hahn B, McCutchan F, Mellors JW, Sodroski J. 1998, Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, pp. III-102-111

    Google Scholar 

Download references

Acknowledgements

We thank Melissa Abrahams for her help in obtaining the transmission sequences from her study. We thank James Arthos for his discussions regarding α4β7. This work was supported by NH&MRC Program Grant 510448. JMcC is supported by an ARC Future Fellowship (110100250).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Damian F J Purcell or James M McCaw.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material

12977_2012_2887_MOESM1_ESM.doc

Additional file 1: Table S1. T-tests performed on groups resulting in significant findings from Table 1. (DOC 28 KB)

Additional file 2: Table S2. Subtype-specific baseline characteristics. (DOC 162 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Mota, T.M., Murray, J.M., Center, R.J. et al. Application of a case–control study design to investigate genotypic signatures of HIV-1 transmission. Retrovirology 9, 54 (2012). https://doi.org/10.1186/1742-4690-9-54

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1742-4690-9-54

Keywords