- Open Access
Breast milk and in utero transmission of HIV-1 select for envelope variants with unique molecular signatures
Retrovirology volume 14, Article number: 6 (2017)
Mother-to-child transmission of human immunodeficiency virus-type 1 (HIV-1) poses a serious health threat in developing countries, and adequate interventions are as yet unrealized. HIV-1 infection is frequently initiated by a single founder viral variant, but the factors that influence particular variant selection are poorly understood.
Our analysis of 647 full-length HIV-1 subtype C and G viral envelope sequences from 22 mother–infant pairs reveals unique genotypic and phenotypic signatures that depend upon transmission route. Relative to maternal strains, intrauterine HIV transmission selects infant variants that have shorter, less-glycosylated V1 loops that are more resistant to soluble CD4 (sCD4) neutralization. Transmission through breastfeeding selects for variants with fewer potential glycosylation sites in gp41, are more sensitive to the broadly neutralizing antibodies PG9 and PG16, and that bind sCD4 with reduced cooperativity. Furthermore, experiments with Affinofile cells indicate that infant viruses, regardless of transmission route, require increased levels of surface CD4 receptor for productive infection.
These data provide the first evidence for transmission route-specific selection of HIV-1 variants, potentially informing therapeutic strategies and vaccine designs that can be tailored to specific modes of vertical HIV transmission.
Despite focused efforts upon preventative measures, pediatric human immunodeficiency virus-type 1 (HIV-1) infections through mother-to-child transmission (MTCT) continue to challenge clinicians and strain healthcare systems, with 1.5 million HIV-positive women giving birth and 240,000 children acquiring the virus in 2013 . HIV-1 MTCT can occur through three distinct routes: in utero (transplacental passage), intrapartum, and postpartum through breastfeeding. Studies of HIV-infected infants indicate that disease progression can vary dramatically among these transmission modes, with infants infected in utero showing the shortest times of survival [2, 3]. An understanding of the genotypic and phenotypic factors that uniquely distinguish these routes of HIV transmission is crucial for developing targeted preventative strategies including vaccines.
A large majority of HIV-1 infections that occur through sexual or vertical transmission are established by a single “founder” variant, with primary infection being characterized by a virus population with significantly less sequence diversity compared to donor viruses [4–17]. The founding variant(s) typically utilize the chemokine receptor CCR5 for entry, and may represent a minority population from the donor [4, 8–10, 14, 18–22]. Studies of heterosexual transmission (subtypes A, C and D) [5, 18, 23, 24] and of MTCT (subtype CRF01_AE)  have identified certain “molecular signatures” in founder variants, such as shorter envelope variable loops, loss of potential N-glycosylation sites (PNGs), or selection of particular PNGs in transmitted strains, but no such signatures have been observed following transmission of subtype B variants [23, 25, 26]. Although previous MTCT studies had limitations of small sample size, a limited focus to specific regions within the envelope, and/or the inability to discriminate between different routes of MTCT [4, 6, 7, 9, 12, 13, 27–33], all of these studies reported detection of a single or few founder variants. Some (but not all) MTCT studies have reported genetic differences between viral strains transmitted in utero and intrapartum [6, 12, 34]. Investigations of the role of maternal neutralizing antibodies in shaping the founder virus population during MTCT have yielded conflicting results: while some data show that selected founder viruses represent maternal antibody escape variants [30, 34], other studies find no such relationship .
As the major viral protein assembly on the virion surface, the HIV envelope protein spike holds the dual distinction of being both the ‘key’ that the virus uses to unlock susceptible target cells in the exposed infant as well as being the primary target of neutralizing antibodies. The HIV-1 envelope spike assembly is composed of a trimer of glycosylated gp120/gp41 heterodimers, and is responsible for mediating entry into host cells via interactions with the CD4 receptor and a co-receptor (typically CCR5). The gp120 subunit lies entirely on the extracellular surface, while the standard envelope topology model describes gp41 as having an extracellular domain, a single membrane-spanning domain, and a cytoplasmic domain (CD) . Analysis of gp120 sequences from diverse viral strains identifies five regions of high sequence variability (variable domains V1–V5) interspersed within five relatively constant domains (C1–C5) [36, 37] (Additional file 1: Figure S1). The tremendous sequence diversity within the gp120 variable regions enables escape from neutralizing antibodies . In addition to CD4 and CCR5, the gp120 subunit can engage other cell-surface proteins that may facilitate mucosal transmission, such as dendritic cell-specific intercellular adhesion molecule-3-grabbing non-integrin (DC-SIGN)  and the integrin α4β7 [40, 41]. Envelope recognition of the CD4 receptor, co-receptor, and other key cell surface proteins, all while simultaneously escaping antibody surveillance, suggests functional sequence constraints that govern the selection of the founding variant(s) from the swarm of donor isolates during transmission.
MTCT differs significantly from other HIV transmission modes: first, MTCT occurs in the presence of maternal antibodies to HIV, with infant levels of maternal antibodies being greater than or equivalent to those in the mother at birth, and then decreasing over the first few months of life . Second, the infant immune system differs in the quantity, distribution, and activation state of CD4+ T cells compared to adults . Third, infants share half their genetic identity with their mothers, leading to greater immunologic overlap between donor and recipient than occurs during sexual or parenteral transmission. Moreover, MTCT in utero, intrapartum and through breastfeeding each involves unique physical exposures, biological barriers, and developmental states of the infant that could influence transmission dynamics and founding strain selection. Given such distinctive transmission features, as well as the critical role of HIV envelope protein in establishing infection and escaping the host immune system, we performed a large, systematic analysis of the genetic and phenotypic properties of HIV envelope variants transmitted from mother to child by different routes, with the aim of identifying viral characteristics that might confer selective advantages for mode-specific transmission.
Transmission timing and phylogenetic linkage of HIV strains within mother-to-child transmission pairs
We examined 22 MTCT pairs participating in the Zambia Exclusive Breastfeeding Study (ZEBS) . Six mothers from this cohort transmitted in utero (IUT), as their infants were HIV-positive at birth by PCR tests, while 13 mothers transmitted virus through breastfeeding (BMT), since their infants were HIV-negative at birth and at 1 month of age, but HIV-positive after 42 days. Three mothers were classified as “indeterminate” transmitters (IND), as their infants were HIV-negative at birth but HIV-positive after 1 month: as such, the precise timing of when these infants acquired virus (either postpartum or during labor and delivery) is unclear. To confirm transmission linkage, we cloned, sequenced and compared 647 envelope genes following polymerase chain reaction (PCR) amplification. In all, we obtained 245 envelope sequences from infants from the first available HIV-1 PCR-positive infant blood sample, and 402 sequences from mothers (Fig. 1). Envelope sequences of maternal origin were typically taken from maternal blood and breast milk samples collected either at study entry (for IUT and IND pairs) or from maternal samples collected within ~35 days prior to the infant’s first positive PCR test (for BMT pairs) (Table 1). In all, we obtained 228 envelope sequences from the IUT group, 365 sequences from the BMT group, and 54 sequences from the IND group (Fig. 1). From within the IUT group, 162 envelope sequences were amplified from blood samples and 66 from milk, while from within the BMT group 309 and 56 envelope sequences were obtained from blood and milk samples, respectively (Table 2). When gp160 sequences from the entire cohort were placed into a phylogenetic tree, each transmission pair formed a monophyletic cluster, with infant envelope sequences forming a sub-cluster within each larger maternal cluster (Fig. 2a). From this analysis, we determined that 21 transmission pairs were infected with subtype C virus, while one pair (Pair 8 from the BMT group) was infected with HIV-1 subtype G (Fig. 1).
Founder populations and reduction in viral diversity upon mother-to-child transmission
Envelope amino acid sequences obtained from infant samples were significantly less diverse than maternal sequences (Fig. 2b), but did not differ significantly in gp160 length (mean of 858.2 residues for all infant strains in the cohort, and 857.4 residues for the maternal strains; p = 0.760) or in number of PNGs (mean of 29.6 and 28.7 PNGs for infant and maternal strains, respectively; p = 0.192), as determined by comparison of means derived from generalized estimating equations (GEE) modeling (Fig. 3a, b). The number of founder variants establishing infection was judged by analyzing phylogenies, viral diversity, insertion and deletion patterns (indels) and the number and linkage of phylogenetically informative mutations in each infant’s virus population (see Additional file 10: Table S1; Additional file 2: Figures S2, Additional file 3: Figure S3). At least 12/22 infections (9/13 BMT, 2/6 IUT, and 1/3 IND) were consistent with a single founder variant. The incidence of putative multi-founder infections was lower in the BMT group than in the IUT group (31% of breast milk transmissions vs. 67% of in utero transmissions), but this difference was not statistically significant (p = 0.319, Fisher’s exact test).
Distinct envelope protein characteristics of maternal and infant virus variants
All infant variants used CCR5 for entry, as determined from infection assays using pseudotyped HIV virions bearing the infant envelope sequences and TZM-bl cells in the presence of high TAK-779 concentrations, and/or from infection of GHOST cells co-expressing CD4 and either CCR5 or CXCR4 (data not shown).
To determine if a critical set of epitopes within HIV-1 envelopes are being selected during MTCT in this cohort, we tested a broad array of envelope-specific inhibitors in infection assays against pseudotyped HIV-1 virions bearing viral envelopes of either infant or maternal origin. These inhibitors included broadly neutralizing monoclonal antibodies (bnAbs) such as the V1/V2 loop-specific PG9 and PG16 [45–47]; 2F5 and 4E10, which bind to the membrane-proximal external region (MPER) of gp41 [48–51]; 2G12, which recognizes carbohydrate moieties on the outer domain of gp120 [48, 52, 53]; and b12, which targets a CD4 binding epitope on gp120 [54–58]. We also tested the inhibitory activities of soluble sCD4 (a 26 kDa protein) , the small molecule CCR5 inhibitor TAK-779 [60, 61], and the 36-amino acid fusion inhibitor enfuvirtide (T-20) , which we used to probe HIV-1 envelope interactions with the CD4 receptor, the CCR5 co-receptor, and the envelope capacity to trigger membrane fusion, respectively. These experiments allowed (a) probing of neutralization sensitivity using standardized reagents, (b) examination of the accessibility of specific neutralizing epitopes (e.g., the gp41 MPER), and (c) assessment of the potential usefulness of available bnAbs in the setting of MTCT (either as a prophylactic intervention or a vaccine template).
In all, we performed inhibition assays with 322 different envelopes from this cohort: 223 obtained from maternal samples and 99 from infant samples. Susceptibility to the bnAbs 4E10, 2F5, b12, 2G12, and T-20 was similar between infant and maternal variants (data not shown). Infant variants also showed similar susceptibility as the maternal variants to PG16 (p = 0.150) and TAK-779 (p = 0.429), but were more sensitive to PG9 (mean IC50 of 0.17 µg/mL for infant strains and 0.38 µg/mL for maternal strains; p = 0.035) and more resistant than the maternal variants to sCD4 (mean IC50 of 13.7 µg/mL for infant strains and 8.8 µg/mL for maternal strains; p = 0.0007) (Fig. 3c–f). When the data were stratified by route of transmission, a unique transmission signature was identified for each route, as described below. To better examine genotypic and phenotypic differences in strains having precisely defined modes of transmission, we focused our analysis on strains solely from the IUT and BMT groups, and did not explore further the envelope variants from the IND group.
‘Transmission signatures’ of in utero and breast milk infection
In contrast to the overall comparisons between maternal and infants viruses described above, when stratified by transmission route, in utero transmission selected for gp160 variants that were shorter (mean length of 852.0 residues in infant variants and 857.6 in maternal variants; p = 0.008) and encoded fewer PNGs (mean of 27.8 sites in infant variants and 29.7 in maternal variants; p = 0.001) (Fig. 4a, b). When such analyses were confined to the V1–V4 region for the IUT group, the infant envelope variants had shorter V1–V4 lengths than maternal variants (mean of 277.0 residues in infant strains and 282.7 residues in maternal strains; p = 0.001) and fewer PNGs (mean of 18.8 sites within V1–V4 in infant strains, and 20.7 sites in V1–V4 of maternal strains; p = 0.001). This observed decrease in V1–V4 length and number of PNGs is similar to that encountered in studies of heterosexual transmission of HIV-1 subtype C , but neither was detected in our breast milk transmission group.
Examination of each individual variable and constant domains (V1–V5 and C1–C5) indicated that infant variants within the IUT group had shorter V1 loops than those found in their maternal counterparts (mean V1 length of 20.1 residues in infant isolates and 24.1 residues in maternal isolates; p = 0.001). Strains from infected infants also had fewer V1 PNGs (mean of 2.4 sites in infant isolates and 3.0 sites in maternal isolates; p = 0.017) and fewer V5 PNGs (mean of 1.2 sites in infant isolates and 1.5 sites in maternal isolates; p = 0.017) (Fig. 4c). No differences in length or number of PNGs for gp160, V1–V4, or V1 regions were observed in the BMT group (Fig. 4d). Infant isolates from the BMT group had fewer PNG sites only within gp41 (mean of 3.9 PNGs in gp41 of infant strains and 4.3 PNGs in gp41 from maternal strains; p = 0.017), while no such difference in gp41 PNGs was observed in the IUT group.
Inhibitor IC50 comparison
In utero transmission failed to select infant-derived variants that differed significantly in bnAb sensitivity from maternal-derived variants, or in sensitivity to the inhibitors TAK-779 (Fig. 5) and T-20 (data not shown). However, GEE analysis of the mean sCD4 IC50’s showed that infant-derived variants from the IUT group were more resistant to sCD4 than corresponding maternal variants (mean sCD4 IC50 of 12.3 µg/mL for virus with infant-derived envelopes and 6.2 µg/mL for virus with maternal-derived envelopes; p = 0.0002) (Fig. 5). We further noted that when confining comparisons to within mother–infant pairs, the infant isolates showed higher median sCD4 IC50 values than those seen for corresponding maternal isolates in 5/6 IUT pairs (Additional file 4: Figure S4a). However, no statistically-significant difference in sCD4 IC50 for maternal and infant isolates within any given IUT mother–infant pair was confirmed (as determined by Wilcoxon rank-sum tests, followed by p value correction for multiple comparisons), perhaps owing to the relatively small number of envelopes being examined within each mother–infant pair. No difference in mean sCD4 sensitivity between infant- and maternal-derived isolates was observed in the BMT group, either in the aggregate or on a pair-by-pair basis (Additional file 4: Figure S4b).
Within the BMT group, infant variants were more susceptible than maternal variants to neutralization by bnAbs PG9 (mean IC50 of 0.16 µg/mL for virus with infant-derived envelopes and 0.68 µg/mL for virus with maternal-derived envelopes; p = 0.0004) and PG16 (mean IC50 of 0.04 µg/mL for virus with infant-derived envelopes and 0.21 µg/mL for virus with maternal-derived envelopes; p = 0.013), a difference that is not observed in the IUT group (Fig. 5). When comparing within BMT mother–infant pairs, the infant variants showed lower median PG9 and PG16 IC50 values than the corresponding maternal variants in 8/13 pairs (Additional file 4: Figure S4b). Using Wilcoxon rank-sum tests, no statistically significant differences in PG9 or PG16 IC50 values for the maternal and infant isolates within any of the BMT pairs were discerned, though this analysis may have been limited by the relatively low number of isolates analyzed within each BMT pair.
In utero MTCT may depend upon initial infection of placental macrophages (e.g., Hofbauer cells), which have low surface densities of CD4 molecules [31, 63, 64]. To further explore the CD4 usage phenotype in our cohort, we investigated whether infant- and maternal-derived envelopes productively engage CD4 with differing efficiencies, and whether any difference in CD4 usage depends upon route of transmission. Towards this end, we used Affinofile cells  to test viral entry properties of 151 molecular envelope clones obtained from the 6 IUT and 11 of the 13 BMT pairs in our cohort (10 BMT pairs infected with HIV-1 subtype C and 1 with subtype G).
Affinofile cells can be induced to express varying levels of surface CD4 in a controllable fashion, allowing for standardized and reproducible infection of cells across a wide, physiologically-relevant range of CD4 surface densities. Affinofile cells were treated with ponasterone to maximally induce CCR5 expression, followed by addition of varying doses of doxycycline to yield controlled levels of CD4 expression. Infectivity curves were constructed by varying CD4 expression/cell over a 100-fold range, infecting cells with 2000 infectious units of virus pseudotyped with either infant or maternal envelope variants, and subsequently determining for each variant the number of CD4 molecules/cell (expressed as the number of α-CD4 antibody binding sites (ABS)/cell) required to achieve 20% of maximal entry (EC20). Analysis of EC20’s showed that infant-derived variants from the IUT group required more CD4 molecules than their maternal counterparts (mean EC20 of 4596 ABS/cell for virus with infant-derived envelopes and 3817 ABS/cell for virus with maternal-derived envelopes; p < 0.0001) (Fig. 6, Additional file 5: Figure S5). Pseudotyped virions from the BMT group showed similar behavior: infant variants in the BMT group required more CD4 molecules to achieve 20% of maximal entry than corresponding maternal strains (mean EC20 of 4409 ABS/cell for infant strains and 3387 ABS/cell for maternal variants; p = 0.001) (Fig. 6, Additional file 5: Figure S5). These results indicate that envelopes from the infant-derived strains in both IUT and BMT groups use CD4 less efficiently than corresponding maternal strains, as the infant strains require more CD4 molecules than their maternal counterparts to achieve comparable levels of target cell infection. Additionally, we evaluated CCR5 usage in a small subset of mother–infant pairs, 2 IUT and 4 BMT. Despite the small numbers, infant variants in the IUT group did require less CCR5 molecules for entry than corresponding maternal strains (mean EC20 of 411 α-CCR5 antibody binding sites (ABS)/cell) required to achieve 20% of maximal entry (EC20ABS/cell for infant strains and 493 ABS/cell for maternal variants; p < 0.001) (Additional file 6: Figures S6, Additional file 7: Figure S7). This difference was not seen in the BMT group (mean EC20 of 370 ABS/cell for infant strains and 553 ABS/cell for maternal variants; p = 0.153) nor when all 6 subjects were evaluated together (Additional file 6: Figures S6, Additional file 7: Figure S7).
Dose–response curve analysis
To further understand inhibition mechanisms for HIV-1 neutralization, and to ascertain functional differences in virions that are transmitted either in utero or through breast milk feeding, we fit inhibition curves of PG9 and sCD4 to the median effect model [66, 67]. In the median effect model, the slope parameter (m) provides a measure of allosteric cooperativity in inhibitor binding for multivalent targets, such as the HIV Env trimer. Slopes greater than 1 (m > 1) indicate positive cooperativity, where the binding of one inhibitor molecule to an Env spike allosterically increases the affinity of other binding sites on the trimeric Env. This favors the binding of additional inhibitor molecules and results in a steep dose-dependent rise in inhibitor activity and a high slope. Slopes less than 1 (m < 1) indicate negative cooperativity, where the inhibitor allosterically decreases the affinity of other binding sites, resulting in a shallow rise in inhibition and a low slope. Slopes equal to 1 (m = 1) indicate non-cooperativity, where an inhibitor has no allosteric effects on the Env structure.
Mean slope parameters obtained for PG9 against infant-derived and maternal-derived virions were both <1 and not significantly different from one another, either within the IUT or within the BMT groups (Additional file 8: Figure S8). Mean slope parameters for sCD4 inhibition of infant and maternal variants were also <1, and not significantly different within IUT group (mean slope of 0.59 for infant isolates and 0.68 for maternal isolates within the IUT group; p = 0.345) (Fig. 7). However, the sCD4 slopes for infant isolates in the BMT group were significantly lower than their maternal counterparts (mean slope of 0.52 for infant isolates and 0.76 for maternal isolates within the BMT group; p = 0.0003) (Additional file 9: Figure S9). Slopes that are <1, as seen in these sCD4 inhibition experiments, suggest that binding of an initial sCD4 molecule discourages binding of additional sCD4 molecules to the trimeric HIV-1 spike complex. Such behavior is characteristic of sCD4 inhibition for infant and maternal variants in both IUT and BMT groups in this cohort, and this negative cooperativity becomes more pronounced in the infant isolates from the BMT group relative to their maternal counterparts.
Our detailed investigation of viral isolates obtained from maternal–infant transmission pairs revealed that genotypic and phenotypic features of HIV-1 envelopes acquired in utero (transplacental transmission) can be clearly distinguished from those acquired from breast milk (mucosal transmission). In utero transmission selects for HIV-1 subtype C infant isolates with shorter, less-glycosylated envelopes: specifically, the V1 domain of in utero-acquired envelope sequences is shorter, and the V1 and V5 domains contain fewer PNGs than corresponding maternal isolates. These observations suggest that envelopes with more compact V1 domains and less-glycosylated V1 and V5 domains impart a selective advantage to virions in establishing in utero infections. A recent study shows that HIV-1 subtype C isolates with compact V1/V2 regions have enhanced fusion efficiency and envelope incorporation into virions , and therefore the shorter V1 region within in utero-acquired strains in our present study may confer improved viral fusion kinetics and envelope uptake, potentially a basis for selection during intrauterine transmission. In contrast, breast milk transmission selects for HIV-1 variants that do not differ in overall length or glycosylation from maternal variants, but have fewer PNG sites in gp41, and are more sensitive to neutralization by the bnAbs PG9 and PG16. These particular genotypic and phenotypic signatures are unique to either intrauterine or breast milk transmission, but not both (Fig. 8).
Placental macrophages represent a population of HIV-susceptible fetal cells in close proximity to the maternal circulatory system, and express relatively low amounts of CD4 . Previous work suggests that a higher affinity for CD4 correlates with increased susceptibility to sCD4 . We therefore hypothesized that infection of placental macrophages would require variants with high CD4 affinity to scavenge the relatively scarce amounts of available CD4, and anticipated that envelope isolates from in utero-infected infants would thus be more sensitive to sCD4 inhibition. Contrary to our expectation, our data indicate that infant strains that emerged following in utero transmission display weaker interactions with the CD4 receptor and a lower binding affinity than maternal strain counterparts. Infection assays with TZM-bl cells show that infant envelope variants from the IUT group have an approximately twofold higher mean sCD4 IC50 value than that of corresponding maternal variants. These findings suggest that these infant-derived envelopes bind to cognate CD4 receptors on the cell surface with lower affinity than those from paired maternal isolates. In addition, our experiments with Affinofile cells show that all infant isolates require more surface CD4 than their maternal counterparts to achieve comparable levels of target cell infection. These data indicate that the transmitted infant variants exhibit a decreased CD4 usage efficiency reflected by both a decreased affinity and a higher stoichiometric requirement for entry. Conceivably, despite having lower CD4 affinity these infant envelope variants may have high CD4 avidity to achieve tight CD4 binding in this setting, exploiting interactions of multiple trimeric spike protein subunits with clusters of CD4 receptors in a contact zone (a so-called “entry claw”) [70, 71] to successfully attach to and infect target cells. Here, an elevated CD4 surface density promotes high HIV-1 envelope-CD4 avidity and ensuing infection. One additional possibility that is consistent with our sCD4 data is that the target fetal cells initially infected during in utero transmission actually express a relatively higher amount of CD4, which may result if a transient inflammatory process occurs during gestation, resulting in elevated CD4 expression and facilitating MTCT despite transmitted isolates having a relatively low CD4 affinity.
Alternatively, in utero-transmitted isolates with weakened CD4 affinity might actually possess a reduced reliance upon CD4, perhaps by assuming envelope conformations that are predisposed to immediate CCR5 engagement. In an infection model that requires little CD4, a reduced CD4 affinity would pose no great fitness penalty: conceivably, any additional surface CD4 density beyond a certain threshold level could help promote infection, but this excess CD4 would not be strictly required. Our dose–response curve analysis of the sCD4 inhibition data using the median effect framework  is consistent with this later interpretation. The mean sCD4 slope parameter for in utero-transmitted infant variants, though not significantly different from that of paired maternal variants, is <1 (mean slope of ~0.6). This low slope parameter suggests that initial envelope trimer engagement of a CD4 molecule discourages subsequent binding of additional CD4 molecules (i.e., negative cooperativity), perhaps by the envelope adopting a conformation that occludes or disturbs the remaining CD4 binding sites on the trimer. Such a conformation might confer a selective advantage by better exposing the CCR5 co-receptor binding site, providing easier access for CCR5 co-receptor and relying less upon the usual CD4 binding event(s) to engage CCR5 and subsequently trigger viral fusion. The low slope parameter observed in our study might also imply slower envelope activation kinetics, which would provide additional time to assemble an entry-competent complex with co-receptor following an initial CD4 binding event.
While most HIV-1 isolates are CD4-dependent, some variants with low or no requirement for CD4 have been described (reviewed in Ref. ). Brain-derived, macrophage-tropic HIV-1 envelope variants with low CD4 requirements and increased fusogenicity have been isolated , and macrophage tropism of simian immunodeficiency virus (SIV) is associated with CD4 independence [74, 75] and reduced CD4 affinity . In addition, data from other Affinofile cell experiments suggest that HIV-1 macrophage tropism requires little CD4 (but is not completely CD4-independent) . Given our data and observations from studies described above, we surmise that selective pressures during in utero transmission yield infant isolates with low CD4 requirements for productive infection, and that are well-poised to infect placental macrophages (e.g., Hofbauer cells) that have lower cell surface CD4 densities, perhaps through exploiting higher CD4 avidity, direct co-receptor binding, and/or having enhanced fusogenic properties.
Our studies provide other insights into factors that select the HIV variants that establish infection. Replication of viral lineages that are genetically distinct from those circulating in the blood can be found within the central nervous system [78, 79], in the genital tract , and within the placenta . Conversely, we and others have previously shown that virus found in maternal blood and in breast milk and are not genetically distinct from one another [81–84], suggesting that the observed genetic restriction during BMT must occur within recipient infant tissues, perhaps involving specific cellular factors.
HIV variants from infants infected in utero in our study were similar to those described in heterosexual transmission of HIV-1 subtypes A and C [5, 23], i.e., having shorter variable loops with fewer PNGs than those found in the transmitting partner. The finding that variants that traverse the placenta have similar genotypic characteristics to those that breach the male and female genital tracts suggests that these different anatomical compartments impose common selective pressures that favor (in each setting) transmission of virions with more compact envelope variable regions having fewer PNGs (Fig. 8). Although the origin of such selective pressures is currently unknown, candidates under consideration include cellular lectins (reviewed in Ref. ), whose presence could potentially favor transmission of under-glycosylated envelope isolates that can escape the innate immune system, or the cellular integrin α4β7, which may facilitate mucosal transmission through strong interactions with less-glycosylated V1/V2 regions . Finally, a strong requirement for virions with elevated envelope incorporation and efficient viral fusion properties during heterosexual and in utero transmission might lead to the selection of isolates having compact V1/V2 regions, which have previously been shown to display this phenotype .
Unlike other forms of HIV transmission, MTCT occurs in the presence of HIV specific antibody. At birth, IgG antibody levels in the infant are equal or greater than those in the mother, but rapidly decline during the first few months of life [42, 87]. Although some studies indicate that founder variants in infants have escaped maternal neutralizing antibody , other studies have shown that neither breadth nor potency of passively acquired maternal antibody protect against breast milk transmission  and that maternal antibody does not strongly influence founder selection . Nevertheless, the search for broadly neutralizing antibodies (bnAbs) that can serve either as a basis for vaccine design or as passively-administered therapeutics remains an active area of research . The increased susceptibility to PG9 and PG16 of subtype C and G isolates from infants infected through breast milk in our study is consistent with observations from a previous report, which also showed that variants from infants infected through breast milk are sensitive to PG9 and PG16 . In addition, subtype CRF01_AE isolates from infants infected perinatally also display increased PG9 and PG16 sensitivity over corresponding maternal strains . Collectively, these findings support the exploration of bnAbs such as PG9 and PG16 in MTCT clinical trials.
In summary, our data indicate that different MTCT routes select genotypically and phenotypically distinct viral variants. These findings may have implications for both HIV pathophysiology and the development of successful MTCT intervention strategies. Past studies suggest that infants infected perinatally experience more rapid disease progression than those infected later by breastfeeding, with the accelerated pathogenesis being attributed to differences in immunologic maturity at the time of infection, and/or to detrimental effects of HIV infection upon infant growth and development [91, 92]. Alternatively, our data suggest that characteristic viral isolates for each transmission setting might infiltrate target tissues through unique mechanisms that may contribute to the observed differences in disease progression. Further comparison analyses of the HIV strains that are selected by different MTCT routes could elucidate these mechanisms, as well as lead to new therapeutic strategies or vaccine designs aimed at blocking vertical transmission in a mode-specific fashion.
We selected 22 women who were enrolled in the Zambia Exclusive Breastfeeding Study (ZEBS)  and who were HIV-positive and pregnant at the time of enrollment. The mothers and infants in our study population received a single-dose of nevirapine peripartum, as per Zambian government guidelines at the time of study enrollment. Infants were also given cotrimoxazole prophylaxis from 6 weeks to 12 months of age . The mode of viral transmission for each mother–infant pair was determined by the timing of the first HIV-1 PCR-positive test for the infant: infants who were showed a PCR-positive test at birth were deemed to have acquired virus in utero (n = 6), while those infants who were PCR-negative at birth and at 1 month of age but PCR-positive after 42 days were considered to be infected through breastfeeding (n = 13). Infants whose HIV-1 PCR tests were negative at birth but positive after 1 month were considered to have acquired virus either intrapartum or early postpartum (n = 3).
Peripheral blood and whole breast milk was collected at regular intervals starting at study enrollment (for maternal blood), and extending to 24 months post-delivery. Peripheral blood mononuclear cells (PBMCs) were separated by Ficoll–Hypaque centrifugation (MP Biomedicals) and preserved with RNAlater (Applied Biosystems) or DNAzol (Invitrogen), followed by freezing of both cells and cell-free plasma. Cell and supernatant/lipid fractions from breast milk were separated by centrifugation and then frozen.
Viral RNA copy numbers were determined from maternal plasma samples using the Roche Amplicor v1.5 standard or ultrasensitive methods (Roche Diagnostics). Full-length gp160 sequences were amplified and cloned from either reverse-transcribed plasma RNA or from DNA extracted from PBMCs or breast milk cells using nested PCR with Vif1/OFM19 (outer) and EnvA/EnvN (inner) primer pairs, as previously described [5, 93]. Multiple independent PCR reactions were performed at or near limiting dilution (i.e. 0.1–10 template copies per reaction) to favor amplification of a single viral envelope sequence. The env sequences from the 647 clones isolated in this cohort were obtained through bidirectional Sanger sequencing, and initially analyzed using Sequencher (Genecodes), MacVector (MacVector Inc), and the LANL website tools (http://www.hiv.lanl.gov/content/sequence/LOCATE/locate.html). Examination of env sequences proceeded as has been described previously . In brief, all chromatograms from Sanger sequencing runs were visually inspected during assembly, and any chromatograms with dual peaks were excluded from further analysis. When multiple clones that resulted from a single PCR reaction were not identical, suggesting that more than one amplifiable template was present during PCR amplification, only one representative clone was selected for further analysis.
Examination of env sequences proceeded much as has been described previously . In brief, all chromatograms from Sanger sequencing runs were visually inspected during assembly, and any chromatograms with dual peaks were excluded from further analysis. All sequences were compared against the HIV-1 database using ViroBLAST  to rule out cross-contamination. Sequences were aligned in MUSCLE v3.7  and refined manually using Geneious (Biomatters, Auckland). A maximum likelihood tree was calculated in PhyML v3.0  using the online tool DIVEIN (http://indra.mullins.microbiol.washington.edu/DIVEIN/) (doi: 10.2144/000113370), which implements the evolutionary model GTR + I + G. The tree was rooted with the subtype B reference sequence HXB2 and represented using FigTree version 1.3.1 (http://tree.bio.ed.ac.uk/software/figtree). Potential N-linked glycosylation sites (PNGs) were identified using the N-Glycosite tool on the LANL website (http://www.hiv.lanl.gov/content/sequence/GLYCOSITE/glycosite.html) .
The number of founder variants establishing infection was inferred by detailed assessment of the infant virus populations, including the number of phylogenetically informative mutation sites (i.e., specific mutations found in more than one sequence, or “InSites”) as determined by DIVEIN; the number of observed sublineages (multiple sequences sharing the same InSites) and the number of mutations defining the sublineage; whether insertions or deletions (InDels) were detected and shared across different sequences; whether probable recombination was observed; and the number of maternal sequences found to represent the most recent common ancestor (MRCA) associated with transmission. The amount of time between PCR-based identification of infant infection and virus population sampling was also considered in making this determination (Additional file 10: Table S1; Additional file 2: Figures S2, Additional file 3: Figure S3).
Cells and reagents
293T retroviral packaging cells were obtained from either the American Type Culture Collection (ATCC) or from Clontech (Lenti-X 293T cells), and TZM-bl cells were obtained through the National Institutes of Health (NIH) AIDS Reagent Program (catalog #8129), courtesy of Dr. John C. Kappes, Dr. Xiaoyun Wu and Tranzyme Inc. TZM-bl is a HeLa clone that constitutively expresses high surface levels of CD4, CCR5, and CXCR4 and contains both luciferase and β-galactosidase reporter genes under control of the HIV-1 LTR promoter [98–101]. TZM-bl and 293T packaging cells were maintained in Dulbecco’s modified Eagle’s medium (DMEM) (Fisher Scientific) supplemented with 10% fetal bovine serum (Gemini Bio-products), 100 U/mL penicillin-100 µg/µL streptomycin (Gibco), and 2 mM l-glutamine (Gibco). GHOST-R5 (expressing CCR5) and GHOST-X4 (expressing CXCR4) cells were also obtained from the NIH AIDS Reagent Program (catalog #3944 and #3685, respectively), courtesy of Dr. Vineet N. KewalRamani and Dr. Dan R. Littman, and were maintained in DMEM with 4.5 g/L glucose (Fisher Scientific), 10% FBS, 500 µg/mL G418 (Invitrogen), 100 µg/mL hygromycin (Invitrogen), and 1 µg/mL puromycin (Sigma-Aldrich). Affinofile cells have been previously described ; briefly, these cells contain independently inducible CD4 and CCR5 expression constructs under the regulatory control of tetracycline and ponasterone, respectively, as well as an HIV-promoter-driven, secreted Gaussia luciferase reporter. Affinofile cells were maintained in Affinofile media (AfM) containing Advanced DMEM (Gibco, Carlsbad, CA) supplemented with 2% dialyzed FBS (Hyclone, Waltham, MA), 2 mM l-glutamine (HiMedia), 100 U/mL penicillin G/100 µg/µL streptomycin (Mediatech), and 50 µg/mL blasticidin (Mediatech). QuantiBRITE™ beads for flow cytometry experiments with Affinofile cells were purchased from BD Biosciences (catalog #340495; Franklin Lakes, NJ). Phycoerythrin (PE)-conjugated mouse anti-human CD4 (clone RPA-T4, catalog #555347), anti-human CCR5 (clone 2D7, catalog #555993), and isotype control (clone MOPC-21, catalog #555749) monoclonal antibodies (mAbs) with a 1:1 antibody:PE ratio were also purchased from BD Biosciences.
The following drugs and antibodies were obtained through the NIH AIDS Reagent Program: HIV-1 gp41 mAbs 4E10 (catalog #10091) and 2F5 (catalog #1475), courtesy of Dr. Hermann Katinger [48–51]; HIV-1 gp120 mAb IgG1 b12 (catalog #2640) from Drs. Dennis Burton and Carlos Barbas [54–57]; TAK-779 from Takeda Chemical Industries, Ltd. (catalog #4983) [60, 61]; T-20 from Roche (catalog #9845); and sCD4-183 (catalog #7356) from Pharmacia, Inc.  The mAbs 4E10 and IgG1 b12 were also purchased directly from Polymun Scientific. The bnAbs PG9 and PG16 were kindly provided through the International AIDS Vaccine Initiative (IAVI), courtesy of Dr. Dennis Burton and the Protocol G team.
Pseudovirus was produced by co-transfection of 293T cells with a plasmid containing a cloned maternal or infant env sequence and the pSG3∆Env backbone (obtained through the NIH AIDS Reagent Program; catalog #11051) using Fugene 6 (Roche Applied Science). The culture medium was clarified by centrifugation at 48–72 h post-transfection, and then aliquoted and stored at −80 °C. One aliquot was thawed and titered on TZM-bl cells. The percentage of functional clones relative to the total number of clones amplified varied from 6 to 90% among the transmission pairs.
Infectivity assays were performed in a 96-well plate format using TZM-bl cells as previously described . Fifty percent inhibitory concentrations (IC50’s) for each tested inhibitor were interpolated from data points immediately above and below the 50% infectivity level in log-linear plots of dose–response curves using Excel (Microsoft). Each infectivity data point in the dose response curve represented the mean of at least two separate tests of a given inhibitor concentration and showed less then 30% variation among these tests.
On Day 0, Affinofile cells were trypsinized (TrypLE; Lifetech, Carlsbad, CA), washed once with AfM, counted, and diluted to 200,000 cells/mL prior to induction. An aliquot of uninduced cells treated with carrier controls (ethanol and deionized water) was replated, and the remaining cells were treated with 2 µM ponasterone (Invitrogen/Life Technologies) to induce high levels of CCR5 expression. CCR5-induced cells were then split into six pools for treatment with doses of doxycycline (Calbiochem, Billerica, MA) ranging between 0.1 and 1.6 ng/mL to induce CD4 expression. Doxycycline-induced cells were plated at 20,000 cells/well in 96-well plates, with a small remainder reserved for culture in a T-25 flask for subsequent flow cytometry analysis.
On Day 1, media was aspirated from all wells and 2000 IU virus was added to duplicate wells of a 96-well plate after being diluted in AfM without diethylaminoethyl (DEAE)-dextran to a total volume of 200 µL, which yielded 10,000–20,000 relative light units (RLUs) per well after a 48 h incubation. AfM without virus was added to uninfected cell control wells. The induced and uninduced Affinofile cells were trypsinized in T25 flasks and washed. Receptor expression levels were subsequently measured by quantitative fluorescence-activated flow cytometry (qFACS) using PE-conjugated α-CD4 and α-CCR5 mAbs, PE-conjugated QuantiBRITE™ beads, and a LSR II flow cytometer (Becton–Dickinson). A standard curve generated from measurements with the QuantiBRITE™ beads having known PE levels allowed calculation of cell surface concentrations of CD4 (expressed as antibody binding sites (ABS)/cell), taking into account that the α-CD4 mAbs used in Affinofile experiments have a 1:1 PE:antibody molar ratio.
On Day 3, 30 µL of luciferase-containing supernatant was transferred from wells containing infected or mock-infected cells to a 96-well white opaque plate and read on a Lumistar Optima instrument (BMG Labtech, Cary, NC) using a Gaussia-luciferase kit (New England Biolabs, Ipswich, MA) per manufacturer’s instructions. RLU values were compared between duplicate wells and any pairs with >25% variation were discarded as invalid. Data from an entire sample were discarded if the 1.6 ng/mL doxycycline wells (high control) were either invalid due to unacceptable variance, or less than threefold higher RLU readout than the cell control well values. RLU values for each pair of wells were averaged and the mean cell-control value subtracted before normalization (as a percentage) to the 1.6 ng/mL doxycycline RLU value for that sample. Percent infectivity values were then pooled across a minimum of three independent experiments for each sample and plotted in Excel (Microsoft, Redmond, WA) against the absolute number of CD4 molecules/cell (as determined through comparison to the QuantiBRITE™-derived standard curve) to construct a single infectivity curve for each isolate across a ~100-fold range of CD4 expression levels. Each infectivity curve was used to calculate the number of CD4 molecules required to achieve 20% of maximum entry (EC20) for a given variant. The 20% value was chosen based on an early analysis of control viruses, before the bulk of the data were analyzed. These assay parameters consistently resulted in detectable, low-level viral entry and good assay-to-assay reproducibility.
CCR5 infectivity curves were generated using the same protocol, except CD4 was induced to saturating levels with 3 ng/mL doxycycline and CCR5 levels were adjusted using a ponasterone concentration range of 3–0.05 µM, giving a 100-fold range of CCR5 surface density (approximately 100–10,000 CCR5 ABS/cell).
Median effect analysis
The median effect framework provides characterization of viral inhibitors through a logarithmic transformation of dose response curves from infection assays [66, 67], as described in Eq. 1:
where f a is the fraction of affected virus, D is the inhibitor concentration, D m represents the inhibitor concentration that causes 50% of the maximum inhibitory effect (equivalently, the inhibitor IC50), and m is a slope parameter that is analogous to the Hill coefficient  and provides a measure of the cooperativity of the system. The median effect model was fit to viral inhibition data using linear regression with in-house perl scripts, providing slope and IC50 values of a given inhibitor for each tested maternal and infant isolate. In all cases, median effect fits were determined from the average of two experimental replicates for each neutralization curve, and isolates for which the squared correlation coefficient (R2) was 0.9 or better were selected for further analysis. Those isolates that could not be neutralized by 50% within the range of inhibitor concentrations used in viral infection assays were not analyzed further.
We used GEE modeling to compare mean values of genotypic and phenotypic features between mother and infant sequences. GEE is a parametric method that takes into account differences in the number of isolates tested per patient, as well as the relatedness of linked maternal and infant isolates. All inhibitor IC50 datasets were log10-transformed prior to GEE analysis to yield a normal distribution, and an exchangeable correlation matrix was used for GEE calculations. Data analysis was conducted using SAS 9.1 (Cary, NC, USA) and R . Diversity comparisons were done using the method of Gilbert et al. , and proportion of putative multi-variant transmissions were compared using Fisher’s exact test. All reported p values are based on GEE testing unless otherwise specified, and were adjusted for multiple comparisons using the method of Benjamini and Hochberg to control for the false discovery rate .
Nucleotide sequences associated with this manuscript have been submitted to GenBank with accession numbers: GU939049–GU939104, GU939106, GU939107, GU939124–GU939142, HM036739–HM036745, HM036747–HM036749, HM036751–HM036789, HM036791, HM036793, HM036794, HM036797–HM036800, HM036802–HM036831, HM036983–HM037012, HM037015–HM037021, HM037023–HM037037 and KY229251–KY229682. Alignments are available at: https://mullinslab.microbiol.washington.edu/publications/nakamura_2016/.
UNAIDS: The Gap Report 2014: children and pregnant women living with HIV. Joint United Nations Programme on HIV/AIDS (UNAIDS); 2014.
Marinda E, Humphrey JH, Iliff PJ, Mutasa K, Nathoo KJ, Piwoz EG, et al. Child mortality according to maternal and infant HIV status in Zimbabwe. Pediatr Infect Dis J. 2007;26(6):519–26. doi:10.1097/01.inf.0000264527.69954.4c.
Little K, Thorne C, Luo C, Bunders M, Ngongo N, McDermott P, et al. Disease progression in children with vertically-acquired HIV infection in sub-Saharan Africa: reviewing the need for HIV treatment. Curr HIV Res. 2007;5(2):139–53.
Ahmad N, Baroudy BM, Baker RC, Chappey C. Genetic analysis of human immunodeficiency virus type 1 envelope V3 region isolates from mothers and infants after perinatal transmission. J Virol. 1995;69(2):1001–12.
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(5666):2019–22. doi:10.1126/science.1093137.
Dickover RE, Garratty EM, Plaeger S, Bryson YJ. Perinatal transmission of major, minor, and multiple maternal human immunodeficiency virus type 1 variants in utero and intrapartum. J Virol. 2001;75(5):2194–203. doi:10.1128/JVI.75.5.2194-2203.2001.
Kampinga GA, Simonon A, Van de Perre P, Karita E, Msellati P, Goudsmit J. Primary infections with HIV-1 of women and their offspring in Rwanda: findings of heterogeneity at seroconversion, coinfection, and recombinants of HIV-1 subtypes A and C. Virology. 1997;227(1):63–76. doi:10.1006/viro.1996.8318.
Keele BF, Giorgi EE, Salazar-Gonzalez JF, Decker JM, Pham KT, Salazar MG, et al. Identification and characterization of transmitted and early founder virus envelopes in primary HIV-1 infection. Proc Natl Acad Sci USA. 2008;105(21):7552–7. doi:10.1073/pnas.0802203105.
Kwiek JJ, Russell ES, Dang KK, Burch CL, Mwapasa V, Meshnick SR, et al. The molecular epidemiology of HIV-1 envelope diversity during HIV-1 subtype C vertical transmission in Malawian mother–infant pairs. AIDS. 2008;22(7):863–71. doi:10.1097/QAD.0b013e3282f51ea0.
Pasquier C, Cayrou C, Blancher A, Tourne-Petheil C, Berrebi A, Tricoire J, et al. Molecular evidence for mother-to-child transmission of multiple variants by analysis of RNA and DNA sequences of human immunodeficiency virus type 1. J Virol. 1998;72(11):8493–501.
Salazar-Gonzalez JF, Bailes E, Pham KT, Salazar MG, Guffey MB, Keele BF, et al. Deciphering human immunodeficiency virus type 1 transmission and early envelope diversification by single-genome amplification and sequencing. J Virol. 2008;82(8):3952–70. doi:10.1128/JVI.02660-07.
Samleerat T, Braibant M, Jourdain G, Moreau A, Ngo-Giang-Huong N, Leechanachai P, et al. Characteristics of HIV type 1 (HIV-1) glycoprotein 120 env sequences in mother–infant pairs infected with HIV-1 subtype CRF01_AE. J Infect Dis. 2008;198(6):868–76. doi:10.1086/591251.
Tovanabutra S, de Souza M, Sittisombut N, Sriplienchan S, Ketsararat V, Birx DL, et al. HIV-1 genetic diversity and compartmentalization in mother/infant pairs infected with CRF01_AE. AIDS. 2007;21(8):1050–3. doi:10.1097/QAD.0b013e32810c8cf3.
Wolinsky SM, Wike CM, Korber BT, Hutto C, Parks WP, Rosenblum LL, et al. Selective transmission of human immunodeficiency virus type-1 variants from mothers to infants. Science. 1992;255(5048):1134–7.
Zhu T, Mo H, Wang N, Nam DS, Cao Y, Koup RA, et al. Genotypic and phenotypic characterization of HIV-1 patients with primary infection. Science. 1993;261(5125):1179–81.
Delwart EL, Sheppard HW, Walker BD, Goudsmit J, Mullins JI. Human immunodeficiency virus type 1 evolution in vivo tracked by DNA heteroduplex mobility assays. J Virol. 1994;68(10):6672–83.
Gottlieb GS, Heath L, Nickle DC, Wong KG, Leach SE, Jacobs B, et al. HIV-1 variation before seroconversion in men who have sex with men: analysis of acute/early HIV infection in the multicenter AIDS cohort study. J Infect Dis. 2008;197(7):1011–5. doi:10.1086/529206.
Sagar M, Laeyendecker O, Lee S, Gamiel J, Wawer MJ, Gray RH, et al. Selection of HIV variants with signature genotypic characteristics during heterosexual transmission. J Infect Dis. 2009;199(4):580–9. doi:10.1086/596557.
van’t Wout AB, Kootstra NA, Mulder-Kampinga GA, Albrecht-van Lent N, Scherpbier HJ, Veenstra J, et al. Macrophage-tropic variants initiate human immunodeficiency virus type 1 infection after sexual, parenteral, and vertical transmission. J Clin Invest. 1994;94(5):2060–7. doi:10.1172/JCI117560.
Zhang H, Rola M, West JT, Tully DC, Kubis P, He J, et al. Functional properties of the HIV-1 subtype C envelope glycoprotein associated with mother-to-child transmission. Virology. 2010;. doi:10.1016/j.virol.2009.12.019.
Zhang H, Tully DC, Hoffmann FG, He J, Kankasa C, Wood C. Restricted genetic diversity of HIV-1 subtype C envelope glycoprotein from perinatally infected Zambian infants. PLoS ONE. 2010;5(2):e9294. doi:10.1371/journal.pone.0009294.
Herbeck JT, Rolland M, Liu Y, McLaughlin S, McNevin J, Zhao H, et al. Demographic processes affect HIV-1 evolution in primary infection before the onset of selective processes. J Virol. 2011;85(15):7523–34. doi:10.1128/JVI.02697-10.
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(10):6528–31. doi:10.1128/JVI.79.10.6528-6531.2005.
Ping LH, Joseph SB, Anderson JA, Abrahams MR, Salazar-Gonzalez JF, Kincer LP, et al. Comparison of viral Env proteins from acute and chronic infections with subtype C human immunodeficiency virus type 1 identifies differences in glycosylation and CCR5 utilization and suggests a new strategy for immunogen design. J Virol. 2013;87(13):7218–33. doi:10.1128/JVI.03577-12.
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(10):6523–7. doi:10.1128/JVI.79.10.6523-6527.2005.
Liu Y, Curlin ME, Diem K, Zhao H, Ghosh AK, Zhu H, et al. Env length and N-linked glycosylation following transmission of human immunodeficiency virus Type 1 subtype B viruses. Virology. 2008;374(2):229–33. doi:10.1016/j.virol.2008.01.029.
Matala E, Crandall KA, Baker RC, Ahmad N. Limited heterogeneity of HIV type 1 in infected mothers correlates with lack of vertical transmission. AIDS Res Hum Retroviruses. 2000;16(15):1481–9. doi:10.1089/088922200750006001.
Scarlatti G, Leitner T, Halapi E, Wahlberg J, Marchisio P, Clerici-Schoeller MA, et al. Comparison of variable region 3 sequences of human immunodeficiency virus type 1 from infected children with the RNA and DNA sequences of the virus populations of their mothers. Proc Natl Acad Sci USA. 1993;90(5):1721–5.
Verhofstede C, Demecheleer E, De Cabooter N, Gaillard P, Mwanyumba F, Claeys P, et al. Diversity of the human immunodeficiency virus type 1 (HIV-1) env sequence after vertical transmission in mother-child pairs infected with HIV-1 subtype A. J Virol. 2003;77(5):3050–7.
Wu X, Parast AB, Richardson BA, Nduati R, John-Stewart G, Mbori-Ngacha D, et al. Neutralization escape variants of human immunodeficiency virus type 1 are transmitted from mother to infant. J Virol. 2006;80(2):835–44. doi:10.1128/JVI.80.2.835-844.2006.
Kumar SB, Handelman SK, Voronkin I, Mwapasa V, Janies D, Rogerson SJ, et al. Different regions of HIV-1 subtype C env are associated with placental localization and in utero mother-to-child transmission. J Virol. 2011;85(14):7142–52. doi:10.1128/JVI.01955-10.
Russell ES, Kwiek JJ, Keys J, Barton K, Mwapasa V, Montefiori DC, et al. The genetic bottleneck in vertical transmission of subtype C HIV-1 is not driven by selection of especially neutralization-resistant virus from the maternal viral population. J Virol. 2011;85(16):8253–62. doi:10.1128/JVI.00197-11.
Russell ES, Ojeda S, Fouda GG, Meshnick SR, Montefiori D, Permar SR, et al. Short communication: HIV type 1 subtype C variants transmitted through the bottleneck of breastfeeding are sensitive to new generation broadly neutralizing antibodies directed against quaternary and CD4-binding site epitopes. AIDS Res Hum Retroviruses. 2013;29(3):511–5. doi:10.1089/AID.2012.0197.
Dickover R, Garratty E, Yusim K, Miller C, Korber B, Bryson Y. Role of maternal autologous neutralizing antibody in selective perinatal transmission of human immunodeficiency virus type 1 escape variants. J Virol. 2006;80(13):6525–33. doi:10.1128/JVI.02658-05.
Checkley MA, Luttge BG, Freed EO. HIV-1 envelope glycoprotein biosynthesis, trafficking, and incorporation. J Mol Biol. 2011;410(4):582–608. doi:10.1016/j.jmb.2011.04.042.
Starcich BR, Hahn BH, Shaw GM, McNeely PD, Modrow S, Wolf H, et al. Identification and characterization of conserved and variable regions in the envelope gene of HTLV-III/LAV, the retrovirus of AIDS. Cell. 1986;45(5):637–48.
Willey RL, Rutledge RA, Dias S, Folks T, Theodore T, Buckler CE, et al. Identification of conserved and divergent domains within the envelope gene of the acquired immunodeficiency syndrome retrovirus. Proc Natl Acad Sci USA. 1986;83(14):5038–42.
Pantophlet R, Burton DR. GP120: target for neutralizing HIV-1 antibodies. Annu Rev Immunol. 2006;24:739–69. doi:10.1146/annurev.immunol.24.021605.090557.
Geijtenbeek TB, Kwon DS, Torensma R, van Vliet SJ, van Duijnhoven GC, Middel J, et al. DC-SIGN, a dendritic cell-specific HIV-1-binding protein that enhances trans-infection of T cells. Cell. 2000;100(5):587–97.
Arthos J, Cicala C, Martinelli E, Macleod K, Van Ryk D, Wei D, et al. 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–9. doi:10.1038/ni1566.
Cicala C, Martinelli E, McNally JP, Goode DJ, Gopaul R, Hiatt J, et al. 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. Proc Natl Acad Sci USA. 2009;106(49):20877–82. doi:10.1073/pnas.0911796106.
Stiehm ER. Immunologic disorders in infants and children. 3rd ed. Philadelphia: Saunders; 1989.
Mo H, Monard S, Pollack H, Ip J, Rochford G, Wu L, et al. Expression patterns of the HIV type 1 coreceptors CCR5 and CXCR4 on CD4 + T cells and monocytes from cord and adult blood. AIDS Res Hum Retroviruses. 1998;14(7):607–17.
Kuhn L, Aldrovandi GM, Sinkala M, Kankasa C, Semrau K, Mwiya M, et al. Effects of early, abrupt weaning on HIV-free survival of children in Zambia. N Engl J Med. 2008;359(2):130–41. doi:10.1056/NEJMoa073788.
Walker LM, Phogat SK, Chan-Hui PY, Wagner D, Phung P, Goss JL, et al. Broad and potent neutralizing antibodies from an African donor reveal a new HIV-1 vaccine target. Science. 2009;326(5950):285–9. doi:10.1126/science.1178746.
Pancera M, McLellan JS, Wu X, Zhu J, Changela A, Schmidt SD, et al. Crystal structure of PG16 and chimeric dissection with somatically related PG9: structure–function analysis of two quaternary-specific antibodies that effectively neutralize HIV-1. J Virol. 2010;84(16):8098–110. doi:10.1128/JVI.00966-10.
McLellan JS, Pancera M, Carrico C, Gorman J, Julien JP, Khayat R, et al. Structure of HIV-1 gp120 V1/V2 domain with broadly neutralizing antibody PG9. Nature. 2011;480(7377):336–43. doi:10.1038/nature10696.
Buchacher A, Predl R, Strutzenberger K, Steinfellner W, Trkola A, Purtscher M, et al. Generation of human monoclonal antibodies against HIV-1 proteins; electrofusion and Epstein–Barr virus transformation for peripheral blood lymphocyte immortalization. AIDS Res Hum Retroviruses. 1994;10(4):359–69.
Purtscher M, Trkola A, Grassauer A, Schulz PM, Klima A, Dopper S, et al. Restricted antigenic variability of the epitope recognized by the neutralizing gp41 antibody 2F5. AIDS. 1996;10(6):587–93.
Purtscher M, Trkola A, Gruber G, Buchacher A, Predl R, Steindl F, et al. A broadly neutralizing human monoclonal antibody against gp41 of human immunodeficiency virus type 1. AIDS Res Hum Retroviruses. 1994;10(12):1651–8.
Stiegler G, Kunert R, Purtscher M, Wolbank S, Voglauer R, Steindl F, et al. A potent cross-clade neutralizing human monoclonal antibody against a novel epitope on gp41 of human immunodeficiency virus type 1. AIDS Res Hum Retroviruses. 2001;17(18):1757–65. doi:10.1089/08892220152741450.
Trkola A, Pomales AB, Yuan H, Korber B, Maddon PJ, Allaway GP, et al. Cross-clade neutralization of primary isolates of human immunodeficiency virus type 1 by human monoclonal antibodies and tetrameric CD4-IgG. J Virol. 1995;69(11):6609–17.
Calarese DA, Scanlan CN, Zwick MB, Deechongkit S, Mimura Y, Kunert R, et al. Antibody domain exchange is an immunological solution to carbohydrate cluster recognition. Science. 2003;300(5628):2065–71. doi:10.1126/science.1083182.
Burton DR, Barbas CF 3rd, Persson MA, Koenig S, Chanock RM, Lerner RA. A large array of human monoclonal antibodies to type 1 human immunodeficiency virus from combinatorial libraries of asymptomatic seropositive individuals. Proc Natl Acad Sci USA. 1991;88(22):10134–7.
Barbas CF 3rd, Bjorling E, Chiodi F, Dunlop N, Cababa D, Jones TM, et al. Recombinant human Fab fragments neutralize human type 1 immunodeficiency virus in vitro. Proc Natl Acad Sci USA. 1992;89(19):9339–43.
Roben P, Moore JP, Thali M, Sodroski J, Barbas CF 3rd, Burton DR. Recognition properties of a panel of human recombinant Fab fragments to the CD4 binding site of gp120 that show differing abilities to neutralize human immunodeficiency virus type 1. J Virol. 1994;68(8):4821–8.
Burton DR, Pyati J, Koduri R, Sharp SJ, Thornton GB, Parren PW, et al. Efficient neutralization of primary isolates of HIV-1 by a recombinant human monoclonal antibody. Science. 1994;266(5187):1024–7.
Zhou T, Xu L, Dey B, Hessell AJ, Van Ryk D, Xiang SH, et al. Structural definition of a conserved neutralization epitope on HIV-1 gp120. Nature. 2007;445(7129):732–7. doi:10.1038/nature05580.
Garlick RL, Kirschner RJ, Eckenrode FM, Tarpley WG, Tomich CS. Escherichia coli expression, purification, and biological activity of a truncated soluble CD4. AIDS Res Hum Retroviruses. 1990;6(4):465–79.
Baba M, Nishimura O, Kanzaki N, Okamoto M, Sawada H, Iizawa Y, et al. A small-molecule, nonpeptide CCR5 antagonist with highly potent and selective anti-HIV-1 activity. Proc Natl Acad Sci USA. 1999;96(10):5698–703.
Dragic T, Trkola A, Thompson DA, Cormier EG, Kajumo FA, Maxwell E, et al. A binding pocket for a small molecule inhibitor of HIV-1 entry within the transmembrane helices of CCR5. Proc Natl Acad Sci USA. 2000;97(10):5639–44. doi:10.1073/pnas.090576697.
Wild C, Greenwell T, Matthews T. A synthetic peptide from HIV-1 gp41 is a potent inhibitor of virus-mediated cell-cell fusion. AIDS Res Hum Retroviruses. 1993;9(11):1051–3.
Torres G, Garcia V, Sanchez E, Segarra A, Patterson BK, Melendez-Guerrero LM. Expression of the HIV-1 co-receptors CCR5 and CXCR4 on placental macrophages and the effect of IL-10 on their expression. Placenta. 2001;22(Suppl A):S29–33. doi:10.1053/plac.2001.0652.
Al-Husaini AM. Role of placenta in the vertical transmission of human immunodeficiency virus. J Perinatol. 2009;29(5):331–6. doi:10.1038/jp.2008.187.
Chikere K, Webb NE, Chou T, Borm K, Sterjovski J, Gorry PR, et al. Distinct HIV-1 entry phenotypes are associated with transmission, subtype specificity, and resistance to broadly neutralizing antibodies. Retrovirology. 2014;11:48. doi:10.1186/1742-4690-11-48.
Chou TC, Talalay P. Quantitative analysis of dose–effect relationships: the combined effects of multiple drugs or enzyme inhibitors. Adv Enzyme Regul. 1984;22:27–55.
Shen L, Peterson S, Sedaghat AR, McMahon MA, Callender M, Zhang H, et al. Dose-response curve slope sets class-specific limits on inhibitory potential of anti-HIV drugs. Nat Med. 2008;14(7):762–6. doi:10.1038/nm1777.
Cavrois M, Neidleman J, Santiago ML, Derdeyn CA, Hunter E, Greene WC. Enhanced fusion and virion incorporation for HIV-1 subtype C envelope glycoproteins with compact V1/V2 domains. J Virol. 2014;88(4):2083–94. doi:10.1128/JVI.02308-13.
Thali M, Olshevsky U, Furman C, Gabuzda D, Li J, Sodroski J. Effects of changes in gp120-CD4 binding affinity on human immunodeficiency virus type 1 envelope glycoprotein function and soluble CD4 sensitivity. J Virol. 1991;65(9):5007–12.
Chojnacki J, Staudt T, Glass B, Bingen P, Engelhardt J, Anders M, et al. Maturation-dependent HIV-1 surface protein redistribution revealed by fluorescence nanoscopy. Science. 2012;338(6106):524–8. doi:10.1126/science.1226359.
Sougrat R, Bartesaghi A, Lifson JD, Bennett AE, Bess JW, Zabransky DJ, et al. Electron tomography of the contact between T cells and SIV/HIV-1: implications for viral entry. PLoS Pathog. 2007;3(5):e63. doi:10.1371/journal.ppat.0030063.
Bhattacharya J, Peters PJ, Clapham PR. CD4-independent infection of HIV and SIV: implications for envelope conformation and cell tropism in vivo. AIDS. 2003;17(Suppl 4):S35–43.
Rossi F, Querido B, Nimmagadda M, Cocklin S, Navas-Martin S, Martin-Garcia J. The V1–V3 region of a brain-derived HIV-1 envelope glycoprotein determines macrophage tropism, low CD4 dependence, increased fusogenicity and altered sensitivity to entry inhibitors. Retrovirology. 2008;5:89. doi:10.1186/1742-4690-5-89.
Puffer BA, Pohlmann S, Edinger AL, Carlin D, Sanchez MD, Reitter J, et al. CD4 independence of simian immunodeficiency virus Envs is associated with macrophage tropism, neutralization sensitivity, and attenuated pathogenicity. J Virol. 2002;76(6):2595–605.
Means RE, Matthews T, Hoxie JA, Malim MH, Kodama T, Desrosiers RC. Ability of the V3 loop of simian immunodeficiency virus to serve as a target for antibody-mediated neutralization: correlation of neutralization sensitivity, growth in macrophages, and decreased dependence on CD4. J Virol. 2001;75(8):3903–15. doi:10.1128/JVI.75.8.3903-3915.2001.
Ryzhova E, Whitbeck JC, Canziani G, Westmoreland SV, Cohen GH, Eisenberg RJ, et al. Rapid progression to simian AIDS can be accompanied by selection of CD4-independent gp120 variants with impaired ability to bind CD4. J Virol. 2002;76(15):7903–9.
Joseph SB, Arrildt KT, Swanstrom AE, Schnell G, Lee B, Hoxie JA, et al. Quantification of entry phenotypes of macrophage-tropic HIV-1 across a wide range of CD4 densities. J Virol. 2014;88(4):1858–69. doi:10.1128/JVI.02477-13.
Sturdevant CB, Dow A, Jabara CB, Joseph SB, Schnell G, Takamune N, et al. Central nervous system compartmentalization of HIV-1 subtype C variants early and late in infection in young children. PLoS Pathog. 2012;8(12):e1003094. doi:10.1371/journal.ppat.1003094.
Schnell G, Joseph S, Spudich S, Price RW, Swanstrom R. HIV-1 replication in the central nervous system occurs in two distinct cell types. PLoS Pathog. 2011;7(10):e1002286. doi:10.1371/journal.ppat.1002286.
Bednar MM, Hauser BM, Ping LH, Dukhovlinova E, Zhou S, Arrildt KT, et al. R5 macrophage-tropic HIV-1 in the male genital tract. J Virol. 2015;89(20):10688–92. doi:10.1128/JVI.01842-15.
Heath L, Conway S, Jones L, Semrau K, Nakamura K, Walter J, et al. Restriction of HIV-1 genotypes in breast milk does not account for the population transmission genetic bottleneck that occurs following transmission. PLoS ONE. 2010;5(4):e10213. doi:10.1371/journal.pone.0010213.
Gray RR, Salemi M, Lowe A, Nakamura KJ, Decker WD, Sinkala M, et al. Multiple independent lineages of HIV-1 persist in breast milk and plasma. AIDS. 2011;25(2):143–52. doi:10.1097/QAD.0b013e328340fdaf.
Salazar-Gonzalez JF, Salazar MG, Learn GH, Fouda GG, Kang HH, Mahlokozera T, et al. Origin and evolution of HIV-1 in breast milk determined by single-genome amplification and sequencing. J Virol. 2011;85(6):2751–63. doi:10.1128/JVI.02316-10.
Gantt S, Carlsson J, Heath L, Bull ME, Shetty AK, Mutsvangwa J, et al. Genetic analyses of HIV-1 env sequences demonstrate limited compartmentalization in breast milk and suggest viral replication within the breast that increases with mastitis. J Virol. 2010;84(20):10812–9. doi:10.1128/JVI.00543-10.
Ballegaard V, Haugaard AK, Garred P, Nielsen SD, Munthe-Fog L. The lectin pathway of complement: advantage or disadvantage in HIV pathogenesis? Clin Immunol. 2014;154(1):13–25. doi:10.1016/j.clim.2014.06.002.
Nawaz F, Cicala C, Van Ryk D, Block KE, Jelicic K, McNally JP, et al. The genotype of early-transmitting HIV gp120 s promotes α4β7-reactivity, revealing α4β +7 /CD4+ T cells as key targets in mucosal transmission. PLoS Pathog. 2011;7(2):e1001301.
Braibant M, Barin F. The role of neutralizing antibodies in prevention of HIV-1 infection: what can we learn from the mother-to-child transmission context? Retrovirology. 2013;10:103. doi:10.1186/1742-4690-10-103.
Lynch JB, Nduati R, Blish CA, Richardson BA, Mabuka JM, Jalalian-Lechak Z, et al. The breadth and potency of passively acquired human immunodeficiency virus type 1-specific neutralizing antibodies does not correlate with risk of infant infection. J Virol. 2011;. doi:10.1128/JVI.02216-10.
Burton DR, Ahmed R, Barouch DH, Butera ST, Crotty S, Godzik A, et al. A blueprint for HIV vaccine discovery. Cell Host Microbe. 2012;12(4):396–407. doi:10.1016/j.chom.2012.09.008.
Thenin S, Samleerat T, Tavernier E, Ngo-Giang-Huong N, Jourdain G, Lallemant M, et al. Envelope glycoproteins of human immunodeficiency virus type 1 variants issued from mother–infant pairs display a wide spectrum of biological properties. Virology. 2012;426(1):12–21. doi:10.1016/j.virol.2012.01.017.
Marston M, Becquet R, Zaba B, Moulton LH, Gray G, Coovadia H, et al. Net survival of perinatally and postnatally HIV-infected children: a pooled analysis of individual data from sub-Saharan Africa. Int J Epidemiol. 2011;40(2):385–96. doi:10.1093/ije/dyq255.
Becquet R, Marston M, Dabis F, Moulton LH, Gray G, Coovadia HM, et al. Children who acquire HIV infection perinatally are at higher risk of early death than those acquiring infection through breastmilk: a meta-analysis. PLoS ONE. 2012;7(2):e28510. doi:10.1371/journal.pone.0028510.
Nakamura KJ, Gach JS, Jones L, Semrau K, Walter J, Bibollet-Ruche F, et al. 4E10-resistant HIV-1 isolated from four subjects with rare membrane-proximal external region polymorphisms. PLoS ONE. 2010;5(3):e9786. doi:10.1371/journal.pone.0009786.
Deng W, Nickle DC, Learn GH, Maust B, Mullins JI. ViroBLAST: a stand-alone BLAST web server for flexible queries of multiple databases and user’s datasets. Bioinformatics. 2007;23(17):2334–6. doi:10.1093/bioinformatics/btm331.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7. doi:10.1093/nar/gkh340.
Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52(5):696–704.
Zhang M, Gaschen B, Blay W, Foley B, Haigwood N, Kuiken C, et al. Tracking global patterns of N-linked glycosylation site variation in highly variable viral glycoproteins: HIV, SIV, and HCV envelopes and influenza hemagglutinin. Glycobiology. 2004;14(12):1229–46. doi:10.1093/glycob/cwh106.
Takeuchi Y, McClure MO, Pizzato M. Identification of gammaretroviruses constitutively released from cell lines used for human immunodeficiency virus research. J Virol. 2008;82(24):12585–8. doi:10.1128/JVI.01726-08.
Wei X, Decker JM, Liu H, Zhang Z, Arani RB, Kilby JM, et al. Emergence of resistant human immunodeficiency virus type 1 in patients receiving fusion inhibitor (T-20) monotherapy. Antimicrob Agents Chemother. 2002;46(6):1896–905.
Derdeyn CA, Decker JM, Sfakianos JN, Wu X, O’Brien WA, Ratner L, et al. Sensitivity of human immunodeficiency virus type 1 to the fusion inhibitor T-20 is modulated by coreceptor specificity defined by the V3 loop of gp120. J Virol. 2000;74(18):8358–67.
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(4):2855–64.
Hill AV. The possible effects of the aggregation of the molecules of haemoglobin on its dissociation curves. J Physiol. 1910;40:4–7.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014.
Gilbert PB, Rossini AJ, Shankarappa R. Two-sample tests for comparing intra-individual genetic sequence diversity between populations. Biometrics. 2005;61(1):106–17. doi:10.1111/j.0006-341X.2005.020719.x.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995;57:289–300.
KJN participated in study design and data analysis, performed the major experiments, and drafted the manuscript. ERS and KS performed statistical analyses. TAW performed data analysis and TAW and NHT interpreted data and participated in drafting the manuscript. LH and JIM performed multiple sequence alignment of the viral envelope sequences, constructed the phylogenetic tree, and characterized the number of founder variants, variable loop length and PNG content within envelope sequences. NEW participated in the design of Affinofile cell experiments and median effect analysis. BL participated in the design of Affinofile cell experiments and contributed Affinofile cells to the study. CK, DMT, and LK participated in study design. GMA conceived of the study, participated in study design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.
First and foremost, we would like to thank the participants in the ZEBS trial, without whom none of this work would have been possible. We would also like to thank the staff of the NIH AIDS Reference Reagent Repository (Drs. D. Finzi, O. Sharma, M. Sharp, and J. Fitzgibbon) for outstanding support. We appreciate the insightful discussions and comments provided by Drs. Steven D. Douglas, J. Michael McCune, Maureen M. Goodenow, and John W. Sleasman. We recognize Colin Sevareid for performing experiments. We also acknowledge S. Conway, L. Jones, and J. Hong for technical support.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
All women participating in the study provided written informed consent. ZEBS was approved by human subjects committees at the investigators’ institutions in the United States (Boston University, Columbia University, and Children’s Hospital Los Angeles), as well as by the University of Zambia Research Ethics Committee.
This work was supported by the National Institutes of Health Cellular, Biochemical, and Molecular Sciences Training Program Grant (T32 067587), the National Institute of Child Health and Human Development (RO1 HD 39611, RO1 HD 40777), and International Maternal Pediatric Adolescent AIDS Trials Group (U01 A1068632) (GMA), as well as by the National Institute of Allergy and Infectious Diseases (R37 AI047734) and the Seattle Center for AIDS Research Molecular Profiling and Computational Biology Core (P30 AI027757) (JIM). GMA is an Elizabeth Glaser Pediatric AIDS Foundation Scientist. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Additional file 1: Figure S1.Domains of the HIV-1 envelope proteins. SP, signal peptide; C1–C5, conserved domains; V1–V5, variable domains; ECD, gp41 extracellular domain; TM, gp41 transmembrane domain; CD, gp41 cytoplasmic domain.
Additional file 2: Figure S2.Pairwise comparisons of genetic diversity within infant founder virus populations. All pairwise comparisons are shown scaled to the same genetic divergence in each panel on the X-axis, with bin intervals of 0.002, but with differentially scaled Y axes that reflect the variable number of sequences obtained from each infant.
Additional file 3: Figure S3.Phylogenetically informative mutation sites and assessment of single versus multiple virus founders in infants. (a) Number of InSites in the infant virus founder populations and the assessment of single versus multiple founders. Note that tally of these sites were only one of several criteria used to assess funder numbers (see Additional file 10: Table S1 for full list of criteria). (b, c) Phylogenetic trees and alignment of phylogenetically informative sites for sequences from the two infants in panel (a) with 6 such sites. These were the only two infants for which the data across all criteria were not unambiguous. While both had evidence of multiple variants, the presence of two founders defined by GGC versus AAA (the latter in the lower infant clade) was clearer in B19, whereas the multiple variants in B16 could have resulted from immune escape when sampled at 88 days of virus positivity.
Additional file 4: Figure S4.sCD4, TAK-779, PG9, and PG16 IC50 data for maternal- and infant-derived envelopes, presented on a pair-by-pair basis. (a) the IUT group and (b) the BMT group. Median IC50 values are indicated by a horizontal bar. Wilcoxon rank-sum tests were used to compare IC50 values obtained from maternal and infant isolates from a given maternal–infant pair, with resulting p values being corrected for multiple comparisons using the method of Benjamini and Hochberg . Any p values showing a significant difference between IC50’s from the maternal- and infant-derived variants (p ≤ 0.05) are indicated in red.
Additional file 5: Figure S5.Representative CD4 infectivity curves using Affinofile cells for IUT (top) and BMT (bottom) maternal–infant pairs. Affinofile cells were induced to generate a 100-fold range of CD4 surface density (ABS/cell) and infected with 2000 IU pseudotyped virus. Percent infection was measured as the percent luciferase relative to infected and maximally induced Affinofile cells. Data shown are representative curves among 3–4 experimental replicates.
Additional file 6: Figure S6.Surface CCR5 molecules needed to produce 20% of maximum infectivity (EC20), by subject and transmission type. The number of envelope sequences tested for each group (n) is shown. The grey box indicates the 95% confidence interval and the whiskers indicate the 99% confidence interval.
Additional file 7: Figure S7.CCR5 infectivity curves using Affinofile cells for IUT (top) and BMT (bottom) maternal–infant pairs. Affinofile cells were induced to generate a 100-fold range of CCR5 surface density (ABS/cell) and infected with 2000 IU pseudotyped virus. Percent infection was measured as the percent luciferase relative to infected and maximally induced Affinofile cells. Data shown are the average of three replicates.
Additional file 8: Figure S8.Representative PG9 neutralization curves for IUT (top) and BMT (bottom) maternal–infant pairs. Data shown are representative curves among 3–4 experimental replicates.
Additional file 9: Figure S9.Representative sCD4 neutralization curves for IUT (left) and BMT (right) maternal–infant pairs. Data shown are representative curves among 3–4 experimental replicates.
Additional file 10: Table S1.Identification of founder variants establishing infection. Infant infections that were inferred to be associated with single (S) or multiple (M) founder variants are shown, along with data used to make these determinations. Given the varying time between the first PCR-positive HIV test and virus population sampling in the infected infant, the wide range in the number of maternal and infant viral sequences sampled, and the unknown impact of early immune responses in selecting for mutations in the virus population, multiple criteria were considered when assessing whether transmission was associated with single or multiple founder variants. These criteria included (i) the number of phylogenetically informative mutation sites (InSites) computed by DIVEIN (doi: 10.2144/000113370), (ii) the minimum number of viral sublineages (multiple sequences sharing the same InSites), (iii) the estimated number of shared mutations defining these sublineages, (iv) whether insertions or deletions (InDels) were detected and shared across different sequences, (v) whether probable recombination was observed between sequences, and (vi) the number of maternal sequences found to represent the most recent common ancestor (MRCA) of the infant’s infection in phylogenetic trees. Multiple founders were inferred when the infant virus population emerged from more than a single MRCA in the mother, and/or when sublineages were detected that differed by 4 or more mutations (including InDels). Supporting these inferences were the detection of recombinants between sublineage sequences. The number of sublineages and mutations defining these sublineages were in some cases not precisely defined due to recombination. While in utero transmission showed a higher incidence of multiple founder virus outgrowth (67% in the IUT group versus 31% in the BMT group), our sampling was insufficient to determine a statistically-significant difference (p = 0.319, Fisher’s exact test).
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Nakamura, K.J., Heath, L., Sobrera, E.R. et al. Breast milk and in utero transmission of HIV-1 select for envelope variants with unique molecular signatures. Retrovirology 14, 6 (2017). https://doi.org/10.1186/s12977-017-0331-z
- Mother-to-child transmission
- Broadly neutralizing antibodies