Simian immunodeficiency virus (SIV) envelope quasispecies transmission and evolution in infant rhesus macaques after oral challenge with uncloned SIVmac251: increased diversity is associated with neutralizing antibodies and improved survival in previously immunized animals
Virology Journal volume 2, Article number: 11 (2005)
Oral infection of infant macaques with simian immunodeficiency virus (SIV) is a useful animal model to test interventions to reduce postnatal HIV transmission via breast-feeding. We previously demonstrated that immunization of infant rhesus macaques with either modified vaccinia virus Ankara (MVA) expressing SIV Gag, Pol and Env, or live-attenuated SIVmac1A11 resulted in lower viremia and longer survival compared to unimmunized controls after oral challenge with virulent SIVmac251 (Van Rompay et al., J. Virology 77:179–190, 2003). Here we evaluate the impact of these vaccines on oral transmission and evolution of SIV envelope variants.
Limiting dilution analysis of SIV RNA followed by heteroduplex mobility assays of the V1–V2 envelope (env) region revealed two major env variants in the uncloned SIVmac251 inoculum. Plasma sampled from all infants 1 week after challenge contained heterogeneous SIV env populations including one or both of the most common env variants in the virus inoculum; no consistent differences in patterns of env variants were found between vaccinated and unvaccinated infants. However, SIV env variant populations diverged in most vaccinated monkeys 3 to 5 months after challenge, in association with the development of neutralizing antibodies.
These patterns of viral envelope diversity, immune responses and disease course in SIV-infected infant macaques are similar to observations in HIV-infected children, and underscore the relevance of this pediatric animal model. The results also support the concept that neonatal immunization with HIV vaccines might modulate disease progression in infants infected with HIV by breast-feeding.
The continued need for breast-feeding in developing countries due to nutritional or socio-economic reasons poses a considerable risk for postnatal mother-to-child transmission of HIV, and breastfeeding is estimated to account for 33–50% of infant HIV infections worldwide [1–5]. This dilemma underscores the need for a vaccine that, when administered shortly after birth to the infant, could protect against HIV transmission via breast-feeding. The ultimate goal of a neonatal HIV vaccine is to prevent infection; however, vaccination of newborns of HIV-infected women early in life may elicit HIV-specific immune responses that substantially reduce infant disease progression in the event that breast milk transmission occurs.
Advances in the understanding of the mechanisms of oral transmission of HIV variants may aid the development of an effective infant HIV-1 vaccine. Recent studies have demonstrated that infants of HIV-infected women can be infected with single or multiple HIV variants [6, 7] shortly before or during the birth process. However, little is known regarding the diversity of HIV transmitted by breastfeeding. These questions are difficult to address in human studies because the characteristics of HIV variants in breast-milk at the time of transmission are unknown. In addition, it is often difficult to obtain virus from infants at early times after HIV infection. Finally, the presence in infants of different levels of transplacentally transferred HIV-specific maternal antibodies with differing anti-viral properties complicates assessments of HIV variant transmission.
Longitudinal studies of HIV-infected adults have shown that the rate of disease progression is inversely related to the rate of evolution of HIV envelope quasispecies [8, 9]. Also, without antiviral treatment, virus-specific immune responses are directly related to HIV quasispecies evolution . The reported relationship between HIV envelope variant evolution and disease progression in HIV-infected infants and children is contradictory. Some studies have found greater HIV envelope variant evolution in rapid progressors [11–13] while other investigations have found that slowly progressing HIV-infected children have greater HIV quasispecies divergence or diversity over time [14, 15]. However, all of these retrospective studies necessarily evaluated HIV variant evolution in a limited number of serial blood samples during the first months of life from a small number of HIV-infected children (two to six per cohort). More recently, a longitudinal study of 10 perinatally HIV-infected children found that changes in HIV envelope quasispecies during the first year of life were associated with a better clinical outcome . A few reports have described a correlation between nascent HIV-specific immune responses, the evolution of HIV variants and disease progression in HIV-infected infants [16, 17].
Simian immunodeficiency virus (SIV) infection of infant macaques is a useful and relevant animal model of pediatric HIV infection for rapidly testing the efficacy of pediatric HIV vaccine and drug interventions [18–20]. This SIV/infant macaque model was previously used to assess the efficacy of two vaccines, (i) modified vaccinia virus Ankara (MVA) expressing SIV Gag, Pol and Env (MVA-SIVgpe) and (ii) live-attenuated SIVmac1A11, against oral challenge with virulent uncloned SIVmac251. We reported an improved clinical outcome (i.e., disease-free survival) for vaccinated compared with unvaccinated infants, which was associated with reduced plasma SIV RNA and sustained SIV-specific humoral immune responses . Here in this report, we used a heteroduplex mobility assay (HMA) to evaluate the genetic diversity in the V1–V2 envelope (env) region of SIV variants present in the SIVmac251 virus inoculum and compare the transmission and evolution of the SIV env quasispecies in plasma following oral inoculation of these vaccinated and unvaccinated infant macaques. Three major questions were addressed: (i) Compared to the SIVmac251 virus inoculum, are few SIV envelope variants transmitted orally?, (ii) Is the lower viremia and better clinical outcome of vaccinated infants related to the initial genetic diversity of SIV env quasispecies?, and, (iii) Is the evolution of SIV envelope quasispecies during the course of infection associated with the development of SIV neutralizing antibody? We demonstrate that while the vaccines did not modulate oral transmission of viral variants, an association was found between vaccine-induced enhanced antiviral immune responses, increased env diversity, and a slower disease course. These findings in vaccinated infant macaques are similar to observations in HIV-infected children with slow disease progression and support the relevance of the SIV infant macaque model for developing neonatal vaccine strategies to prevent pediatric HIV infection and AIDS.
Characterization of variants in SIVmac251-5/98 virus stock
HMA analysis revealed that the undiluted SIVmac251-5/98 virus stock was comprised of a diverse population of V1–V2 env variants. To determine the most common variant(s) in the virus stock, six independent serial dilution experiments were conducted. Viral RNA was isolated from 1 ml of virus stock and 10-fold dilution series (undiluted to 10-9) of the RNA were prepared from 6 separate aliquots of virus stock. The resulting RNA was analyzed by RT-PCR and HMA. Figure 1 shows the results of 4 of these 6 separate virus stock dilution/HMA experiments. The observation that multiple heteroduplex bands were observed through the 10-5 or 10-6 dilutions of viral RNA indicates that the undiluted SIVmac251-5/98 stock contains multiple env variants at high frequency. An RT-PCR endpoint (i.e. dilution to a single variant) was not reached in 3 of the 6 dilution experiments. An example of this is shown in dilution series A (Figure 1). In the other 3 dilution series (Fig. 1, series B, C, and F), the last dilution that yielded an RT-PCR product consisted of a homogeneous population of envelope variants represented by one main variant (homoduplex band). This endpoint variant was designated the virus stock endpoint variant (VSEV). The fact that endpoint variants were reached at different dilutions for each dilution series is probably due to the variability at each step of these independently performed experiments.
The VSEV in dilution series B and F had different mobilities on the HMA gel (Fig. 1). Dilution series C resulted in two endpoint variants, one at 10-6 and the other at 10-7; the positions of these two VSEV corresponded to one of each of the two VSEV in dilution series B and F. Thus, the dilution of the virus stock to an RT-PCR endpoint resulted in 4 independent variants (represented by homoduplex bands) that migrated to two different positions on the HMA gels. Based on these positions, the homoduplex bands that migrated furthest were referred to as VSEV-1 and the variants that migrated a shorter distance were designated VSEV-2 (Fig. 1). To confirm that the four endpoint homoduplexes represented only two variants, an HMA mixture experiment was performed, in which all pairwise combinations of the virus stock endpoint variants were mixed prior to HMA . These experiments demonstrated that the two variants designated VSEV-1 are indeed similar (i.e., = 1–2% difference in nucleotides with no insertion/deletion), as their mixtures resulted in the formation of a single homoduplex band on an HMA gel; similarly, the two variants referred to as VSEV-2 are similar (Fig. 2). In contrast, the formation of heteroduplexes and two main homoduplexes in the mixtures of VSEV-1 and VSEV-2 demonstrate that these 2 variants are significantly different from each other (Fig. 2). Thus, VSEV-1 and VSEV-2 are 2 distinct variants that exist at similar frequencies and represent the most common variants in the undiluted SIVmac251-5/98 virus stock. These results are consistent with observations of the virus stock from which SIVmac251-5/98 was made .
Experimental design of animal experiments and summary of outcome
Nineteen newborn rhesus macaques were divided into 5 experimental vaccine groups (table 1). Group 1 (n = 5) consisted of unimmunized control animals. Group 2 (n = 2), group 3 (n = 4) and group 4 (n = 4) were vaccinated with MVA-SIVgpe at 0 and 3 weeks of age; group 4 had maternally-derived SIV antibodies (due to immunization of their mothers with inactivated SIV). Group 5 (n = 4) was immunized with live-attenuated SIVmac1A11 at 0 and 3 weeks of age. As described elsewhere , except for group 2, all other groups were inoculated orally with SIVmac251-5/98 at 4 weeks of age; all these animals became persistently viremic, but the immunized animals had lower virus levels, enhanced antiviral immune responses and a delayed disease course in comparison to the unimmunized animals. Four of the 5 unimmunized infected animals developed AIDS within 14 weeks of age, while the fifth animal needed euthanasia at 28 weeks. Four MVA-SIVgpe-vaccinated SIVmac251-5/98-infected animals developed AIDS by 19 to 27 weeks of age (2 animals of groups 3 and 4 each; table 1). The remaining eight vaccinated SIVmac251-5/98-infected infants, including all four SIVmac1A11-vaccinated animals, were clinically stable at the end of the observation period (28 weeks of age).
Detection of SIV envelope variants in plasma of neonates early after oral inoculation
The genetic diversity of SIV env variant populations in the plasma of the infant monkeys one week after oral inoculation with SIVmac251 was analyzed by HMA (Fig. 3). Each plasma sample was analyzed in replicates (≥ 2) to assure reproducibility of the gel banding patterns. As indicated by the presence of heteroduplex bands, all infants were infected with multiple SIV env variants, indicating that the SIVmac251-5/98 virus stock contained several variants capable of establishing infection by the oral route. However, there were differences in HMA banding patterns. In each group, some animals had several strong heteroduplex bands; this pattern of variant transmission was referred to as infection pattern A (e.g. Fig. 3, animal 31319). In contrast, one or two infants in each group were infected with a genetically more homogenous variant population, consisting of one major variant (homoduplex band), while heteroduplex bands were less pronounced. These monkeys infected with genetically more homogeneous viral populations harbored one of two main env variants, distinguished by different electrophoretic mobilities of the homoduplexes representing these variants. These more homogeneous variant populations were referred to as infection patterns B and C (e.g. Fig. 3, animals 31325 and 31608, respectively). Infection pattern C contained a homoduplex band that migrated slightly slower than the homoduplex band characterizing infection pattern B. One newborn in each vaccine group was infected with a SIV variant of transmission pattern B. Infection pattern C was detected in one newborn of each group except the SIVmac1A11 vaccinates (table 1, group 5). Therefore, no substantial difference was observed among the different vaccine groups in viral genetic diversity in plasma collected 1 week after virus inoculation. However, all but one infant (31540) infected with more homogenous populations of env variants (infection patterns B or C) had 10- to 100-fold lower virus levels one week after SIVmac251 challenge than all but one infant (31378) infected with more heterogeneous populations of SIV variants (transmission pattern A, Table 1). This association of homogeneous viral variants with reduced SIV RNA in plasma at 1 week after infection was statistically significant (P ≤ 0.05; one-sided Fisher's Exact test) but did not persist. From week 2 after challenge throughout the duration of the study, plasma SIV RNA levels showed no correlation with the initial SIV variant pattern detected in plasma. The rate of disease progression in these animals was also not associated with the initial envelope variant transmission patterns (table 1). Further, there was no correlation between the presence of the MHC type I alleles Mamu-A*01 or Mamu-B*01 and the viral variant infection patterns, levels of SIV RNA in plasma, or disease progression (table 1).
To determine which SIV envelope variant was present in the highest frequency in each infection pattern, serial end-point dilution experiments were performed with RNA isolated from plasma collected one week after SIVmac251 challenge, followed by RT-PCR and HMA. Similar to the methods described above, mixture experiments were then performed, including with VSEV-1 and VSEV-2. These experiments demonstrated that 1 week after infection, the most common variants in animals with the more homogenous transmission patterns B and C were similar (i.e., less than 1–2 % difference based on the absence of heteroduplex bands) to VSEV-1 and VSEV-2, respectively (data not shown). The most common variants by end-point dilution in the 11 monkeys with transmission pattern A and A/C were similar to VSEV-1 (5 animals), or VSEV-2 (5 animals) or both (1 animal).
Greater quasispecies diversity in vaccinated compared to control infants during chronic SIV infection
HMA was used to analyze the evolution of genetic diversity of V1–V2 env populations in plasma of the monkeys during the course of infection (1 week after oral SIVmac251-5/98 challenge until euthanasia) (Fig. 4). Results from two standard measures of the nucleotide sequence heterogeneity of V1–V2 env plasma variants derived from the HMA analyses are shown in Fig. 5: (i) entropy (E), an estimate of the overall viral RNA sequence complexity for each sample and, (ii) median mobility shift (MMS), a measure of the SIV quasispecies sequence divergence reflected by the degree of base-pair mismatch after DNA strand re-annealing of strands of envelope variants .
The diversity of SIV env quasispecies in plasma varied among animals at the first sample (1 week after challenge) as indicated by the gel banding pattern (Fig. 4) and entropy measures (Fig. 5). Entropy for SIV env populations was high (> 0.9) for all unvaccinated animals (Group 1) and for 7 of the 12 vaccinated animals (Fig. 5). Initial entropy < 0.9 for vaccinated animals was associated with lower SIV RNA in plasma at 1 week after challenge (P < 0.05; Fisher's Exact Test). No consistent pattern of entropy over the 24 week course of infection was observed; in two of the five controls (31321, 31608) and three of the 12 vaccinates (31533, 31732, 31780) entropy decreased near the time of euthanasia. Overall, there was no association of SIV envelope diversity as measured by entropy with either viral RNA levels or virus-specific neutralizing antibodies in plasma (see below and Fig. 5).
The sequence divergence of SIV envelope variants in plasma of each animal over time was estimated by the MMS, shown in Fig. 5. Four of the five unvaccinated animals had initial MMS values ≥ 0.5 which decreased at varying rates until the time of euthanasia; the remaining control animal (31325) had an initial MMS < 0.3 which did not change significantly over the course of infection (Fig. 5). Thus, in unvaccinated infants the population of SIV env variants in plasma exhibited either no sequence divergence or increasing sequence similarity over time; this observation is consistent with the absence of sustained SIV-specific immune responses in these animals (; see below). There was no association between MMS values and SIV RNA plasma levels for these unimmunized animals.
For 3 of the 4 vaccinated animals that developed AIDS within the observation period of 28 weeks (animals 31732, 31533, 31542), we also observed little change or a decrease of genetic divergence (i.e., as measured by stable or decreasing MMS values) of plasma env variant quasispecies. In contrast, diversification in plasma SIV env variant populations (i.e., a significant increase in MSS values) was observed by 3 to 5 months of infection in 4 of the 8 vaccinated monkeys (31540, 31833, 31856 and 31778) that were still relatively healthy at 28 weeks (Fig. 5). This diversification corresponded to the detection of different patterns of heteroduplex bands and/or fainter homoduplex bands over time (Fig. 4). Although increased diversification seemed to correlate with improved disease-free survival, this diversification was not associated with any obvious changes in plasma virus levels.
SIV neutralizing antibodies in vaccinates correlate with evolution of SIV quasispecies diversity
The possibility that SIV envelope-specific immune responses were associated with the observed plasma SIV RNA levels or evolution of SIV env quasispecies was evaluated by measuring levels of plasma antibodies that neutralized the homologous challenge virus, SIVmac251-5/98, during the course of infection (Fig. 5). SIV neutralizing antibodies were detected in none of the unvaccinated control animals, but in all except one (31777) of the 12 vaccinated animals within 16 to 20 weeks after infection (Fig. 5). Although the presence of SIV neutralizing antibodies was associated with increased survival of the vaccinated animals, no obvious relationship was detected between the SIV neutralizing antibody levels and either SIV RNA plasma levels or entropy over time in vaccinated animals. However, in the 5 animals with increasing sequence divergence (i.e., increasing MMS values; animals 31540, 31526, 31833, 31856 and 31778), neutralizing antibodies were detected around the time that MMS values increased, and the neutralizing antibody response was sustained (i.e., detectable in ≥ 3 plasma samples) in these 5 animals (Fig. 5). In contrast, animals with stable or declining MMS values had sustained (31533, 31542, 31779), transiently detected (31378, 31732, 31780) or undetectable (31319, 31321, 31322, 31325, 31608, 31777) anti-SIV neutralizing antibodies. Thus, a sustained SIV neutralizing antibody response was associated with increased divergence of SIV envelope variants in plasma (P = 0.009; Fisher's Exact test).
The present study is among the most comprehensive longitudinal studies describing SIV envelope variation in vivo following mucosal SIV infection of infant macaques. In this study, we examined the extent of genetic diversity of the SIV envelope variant pool in the plasma of infant macaques that were inoculated orally at 4 weeks of age with an uncloned, genetically diverse virus stock SIVmac251-5/98. In addition, this is the first study to evaluate whether the transmission and evolution of viral variants was modulated by two different SIV vaccines, MVA-SIVgpe and SIVmac1A11, or the presence of maternally-derived anti-SIV antibodies.
HMA analysis revealed that the animals became infected with multiple SIV envelope variant populations, but which predominantly consisted of one of two single envelope variants that were very similar to the two most common variants in the SIVmac251-5/98 stock. These results are consistent with reports of mother-to-infant HIV transmission of multiple variants [23–25], single variants [14, 26, 27] or both [6, 28–31], but inconsistent with studies reporting vertical transmission of single, minor variants [10, 26, 32–34] from the mothers' virus population. This discrepancy could be explained by differences in the HIV inoculum regarding dose, virulence and genetic diversity compared to SIV. In the present study, macaques were inoculated with a relatively high dose of SIVmac251-5/98, while infection of human infants is likely to occur due to exposure to lower amounts of virus. An inherent limitation of studies of vertical transmission of HIV is that the exact timing of infection is usually unknown, and therefore the mothers' population of viral variants at the time of transmission and the source (e.g., breast-milk) and dose of virus is unknown. Our observation that oral exposure of 17 infant macaques to the same dose of the same virus stock resulted in different transmission patterns further underscores the complexity of studying variant transmission in humans, and suggests that the different outcomes observed for vertical transmission of HIV may not necessarily reflect "selection" of HIV variants but may be more a stochastic event. In this context, studies looking at the effect of heterogeneity of viral variants in the HIV-1 infected mother and the rate of vertical transmission have also shown conflicting results [6, 26, 35]. Also, the HIV studies mentioned focused on prenatal or intra-partum transmission, whereas our study modeled postnatal HIV transmission via breastfeeding by oral inoculation of 1-month old infant macaques with SIVmac251-5/98. The route(s) of infection in utero or during birth for individual infants and source of virus (cell-free or cell-associated) is usually unknown, and therefore different mechanisms may be responsible for viral transmission via these routes . Consistent with this view, others have reported that more SIV variants were detected in orally infected newborn macaques than in infants born to SIV-infected female macaques for which transmission occurred in utero  or during the late breast-feeding period .
Neither of the SIV vaccines used in this experiment (MVA-SIVgpe and SIVmac1A11), nor the presence of maternal antibodies in one of the MVA-SIVgpe immunized groups altered which envelope variants were transmitted because in each group, some monkeys became infected with more heterogeneous and others with more homogeneous virus populations. It is possible that neither MVA-SIVgpe nor SIVmac1A11 elicited immune responses that effectively targeted the predominant SIV env variants in the SIVmac251-5/98 stock, or that anti-envelope immune responses were elicited against regions of the envelope other than V1–V2. It is also possible that vaccine-induced immune mechanisms at the time and/or site(s) of initial infection were not potent enough to modulate the variant transmission patterns.
Viral levels in plasma of monkeys with more homogeneous populations of SIV env variants tended to be lower one week after oral inoculation with SIVmac251-5/98. The higher initial virus levels in infants infected with multiple variants may reflect higher replication capacities of diverse variant populations compared to those comprised of one main variant, especially in the initial target cells during the first days of infection. We have observed this previously for adult macaques inoculated intravaginally . From the second week after SIVmac251-5/98 inoculation onwards, however, there was no correlation between viral genetic complexity (measured by entropy) or divergence (measured by MMS) and plasma SIV RNA levels. Thus, once systemic infection was established, virus replication attained similar levels regardless of the initial diversity, and there was no difference in AIDS-free survival times.
Based on the measurement of MMS values, we observed little change or a decrease in genetic divergence of plasma SIV env variant quasispecies in all unvaccinated and most vaccinated animals that developed AIDS within the observation period of 28 weeks. In contrast, diversification in plasma SIV envelope variant populations was observed in 4 of the 8 vaccinated monkeys that were still relatively healthy at 28 weeks. This increased divergence of plasma viral variants at ~3 to 5 months after infection was generally associated with more sustained levels of SIV-specific neutralizing antibodies, and also of SIV Gag and Env-specific antibodies (measured by ELISA, as shown previously ). Similar associations between viral genetic divergence, immune parameters and/or disease progression have been described in HIV-infected adults and children [6, 8, 12, 13, 15, 16, 29, 38–44], and recently also in juvenile macaques following intravenous or intra-rectal SIVsm inoculation . Our studies extend these observations by demonstrating that this correlation of more sustained immune responses, enhanced viral divergence and slower disease progression is also observed in infant macaques following oral SIV infection. Together, these results suggest that the rate of virus evolution is determined by a combination of the extent of virus replication (which induces random mutations due to the error-prone reverse transcriptase) and selection pressures such as antiviral immune responses that promote the outgrowth of new variants. The generation of increasingly divergent viral variants ("immune escape mutants") reflects attempts of the immune system, albeit only partially effective, to control virus replication. In contrast, high viremia and little evolution of viral envelope variants is associated with severe immunodeficiency (and thus little immune selection pressure) and rapid disease progression.
The patterns of SIV env variant transmission and evolution in infant macaques that were inoculated orally with the same SIVmac251-5/98 stock reflect the range of results that is observed in mother-to-infant transmission of HIV, where the dose and genetic diversity of the virus at the time of transmission are unknown. While the vaccines tested here did not modulate oral transmission of viral variants, an association was found between vaccination and enhanced antiviral immune responses, increased env diversity, and a slower disease course. These findings are similar to observations in HIV-infected children with slow disease progression and underscore the relevance of the infant macaque model for developing neonatal vaccine strategies to prevent pediatric HIV infection and AIDS . These results also support the concept that neonatal immunization could prevent rapid disease progression in infants who become HIV-infected by breast-feeding.
Materials and Methods
Infant immunizations, virus inoculations, and sample collection
All newborn rhesus macaques (Macaca mulatta) were from the HIV-2, SIV, type D retrovirus, and simian T-cell lymphotropic virus type 1-free colony at the California National Primate Research Center. Newborn monkeys were hand-reared in a primate nursery, and all animals were housed in accordance with American Association for Accreditation of Laboratory Animal Care standards. We adhered to the "Guide for Care and Use of Laboratory Animals" . When necessary, animals were immobilized with 10 mg/kg ketamine hydrochloride (Parke-Davis, Morris Plains, NJ) injected intramuscularly (IM). EDTA-anticoagulated blood samples were collected regularly for monitoring virologic and immunologic parameters as described previously .
Four newborn macaques had maternally derived SIV antibodies, because their mothers had been immunized and boosted during three or four consecutive pregnancies with whole-inactivated SIVmac251 plus Montanide ISA 51 adjuvant (Seppic, Fairfield, NJ), administered intramuscularly as previously described . One of two SIV vaccines was administered to newborn monkeys at birth and 3 weeks of age: Modified Vaccinia virus Ankara expressing SIVmac239 gag, pol, and env (MVA-SIVgpe) was given to 8 newborn monkeys, including the 4 with maternal antibodies. SIVmac1A11 was given to 4 newborn monkeys. Details about these vaccines are described elsewhere .
At 4 weeks of age, these 17 monkeys were inoculated orally with 2 doses (24 hours apart) of uncloned virulent SIVmac251. Ketamine anesthesia was used for each inoculation. Each dose consisted of 1 ml of undiluted SIVmac251 of a stock designated by lot number -5/98, and was administered atraumatically by dispensing virus slowly into the mouth with a syringe. The SIVmac251-5/98 virus stock used in this study was derived from a previous SIVmac251 stock (lot 8/95) that was serially passaged intravenously in rhesus macaques as described . This SIVmac251-5/98 stock contained 1 × 105 50% tissue culture infective doses (TCID50) and 1.4 × 109 copies of RNA per ml (determined by bDNA assay).
Quantitation of plasma viral RNA
Viral RNA in plasma was quantified using a branched DNA (bDNA) signal amplification assay specific for SIV, with conditions as described previously .
RNA isolation and RT-PCR
RNA was extracted from plasma samples (100–140 μl) using a viral RNA isolation kit (Qiagen, Inc., Valencia, CA) following the manufacturer's protocol. A 590 bp fragment encompassing the V1–V2 region of SIV env was then amplified in a nested RT-PCR assay as previously described .
Analysis of SIV variants by heteroduplex mobility assay (HMA)
Genetic diversity in viral variant populations was analyzed using a modification of the HMA methods described elsewhere [22, 48]. In brief, V1–V2 env fragments were generated by RT-PCR as described above, and the presence of sufficient product was confirmed on a 1.5% agarose gel. The RT-PCR products were then mixed with 1.5 μl of 10× annealing buffer (1 M NaCl, 100 mM Tris, 20 mM EDTA), denatured at 94°C for 2 minutes and placed immediately on wet ice to promote heteroduplex formation. Samples were then run on non-denaturing 5% polyacrylamide gels and stained with ethidium bromide (0.5 μg/ml). Reverse-images of the stained gels were photographed with a digital imaging system (Alpha Innotech Corporation, San Leandro, CA). All gel images were color reversed to enhance visualization of banding patterns (e.g. black bands on white background). The number of heteroduplex bands observed is a measure of SIV envelope diversity in each monkey's virus population (i.e. a large number of bands on a gel corresponds to a large number of V1–V2 variants in the sample). The RT-PCR and HMA analysis on plasma samples was performed in replicates (of at least 2) to assure reproducibility of the gel banding patterns.
To further characterize SIV envelope variants, we used an additional form of HMA analysis that assesses the relative genetic similarity of specific viral variants by comparing the HMA patterns that result from combinations of these variants. Mixtures of SIV RNA from plasma samples were analyzed by HMA to allow estimation of sequence similarity between two different homogeneous virus populations (i.e. the most common inoculum variants and/or plasma variants from infected monkeys). This "mixture analysis" is based on the HMA HIV subtyping protocol developed by Delwart et al.  and was performed as described previously . Briefly, mixtures of equal volumes of SIV V1–V2 env PCR product amplified from two different samples were mixed and subjected to HMA analysis as described above. Mixtures that result in a single homoduplex band are 98–100% identical  in the nucleotide sequence of the PCR fragment analyzed (V1–V2 env region). Mixtures that result in the formation of heteroduplexes are comprised of variant populations with nucleotide sequences that differ by more than 1–2% or have an insertion/deletion (i.e., a single codon length difference will also cause a gel shift) . We have validated this HMA method for SIVmac251 in a previous study .
Calculation of entropy and median mobility shift
All measures of entropy (E) and median mobility shift (MMS) were estimated according to methods described by Delwart et al. . Images from HMA gels were captured with a CCD camera as binary TIFF files and color reversed to enhance visualization of banding patterns (e.g. black bands on white background). Each TIFF gel image file was then opened using Adobe Photoshop Version 6.0 (Adobe Inc., San Jose, CA) and edited to ensure that the lightest inter-lane areas of the gel image had "0" signal intensity (as read by the NIH Image Program and required for the Hdent program described below). Digitized gel lanes were scanned by using the plot profile function of the NIH Image Program (available at http://rsb.info.nih.gov/nih-image). Lane scans within the same gel were of equal length (i.e. same number of pixels) and were recorded from positions immediately below the single-stranded DNA position to immediately below the homoduplex. The signal intensity at each pixel along the scan was transferred to a Microsoft Excel (Richmond, Wash.) file. Because different numbers of pixels per lane were acquired from different gels, each gel was standardized by partitioning into 191 divisions, the smallest number of pixels in the scans under study. This allowed the maximum distinction of fine banding patterns, while permitting unbiased comparison between gels.
The quasispecies diversity for each sample was estimated by calculating a normalized Shannon entropy, a measure of the breadth or spread of the signal distribution in each HMA gel lane, using the HDent program (available at http://www.ma.umist.ac.uk/mrm) as described by Delwart et al. . The Shannon entropy (S) is defined as: S = - Σ (from i = 1 to N)P(i)ln [P(i)], where N is the number of partitions in a lane, and P(i) is the fraction of the total signal in partition i. The maximum possible entropy is ln(N), and we defined the normalized entropy as S/ln(N). The normalized entropy has a range of 0 to 1, where 0 reflects no diversity (all of the signal is in a single partition), and 1 reflects maximum entropy, in which the signal is evenly distributed throughout all partitions in the lane. Thus, entropy is large for lanes with many, closely spaced or overlapping bands and small for lanes with only one band or a few, narrow bands.
Shannon entropy estimates quasispecies genetic diversity by measuring the pattern of SIV V1–V2 env heteroduplex distribution in an HMA gel lane rather than the specific electrophoretic mobility of heteroduplexes. However, the electrophoretic mobility of heteroduplexes through a polyacrylamide gel is proportional to the sequence differences in reannealed DNA strands [38, 48, 49]. The degree of SIV quasispecies envelope sequence divergence among the V1–V2 env variants present in plasma samples was estimated by calculating a median mobility shift (MMS) for each HMA gel lane using the HDent program. The MMS is a measure of the midpoint of the total signal in an HMA gel lane that has values between 0 and 1, where 0 corresponds to the bottom of an HMA gel lane (i.e. nearest homoduplex bands) and 1 corresponds to the top of the lane. Thus, a MMS score of 1 reflects maximum sequence diversity (i.e. all heteroduplexes bands have maximum mobility reduction and no visible homoduplexes); a MMS value of 0 reflects maximum sequence similarity (> 98%) where all signal for a lane is in homoduplex bands and there are no visible heteroduplexes.
Assessment of MHC class I alleles
DNA extracted from lymphoid cells (with QIAamp® DNA mini kit, QIAgen, Valencia, CA) was used to screen for the presence of the rhesus macaque major histocompatibility complex (MHC) class I alleles Mamu A*01 and Mamu B*01, using a PCR-based technique [50, 51].
Neutralizing antibody titers in EDTA-anticoagulated plasma were measured according to methods described previously , except that CEM-R5 cells (i.e. CEM×174 cells expressing CCR5 by transfection; generously provided by James Robinson) were used. Neutralizing antibody titers were expressed as the reciprocal of the plasma dilution at which 50% of cells were protected from virus-induced killing as measured by neutral red uptake. The virus consisted of SIVmac251-5/98 briefly propagated in human PBMC.
Fisher's exact test, performed with Instat v. 2.03 (GraphPad Software, Inc., San Diego, CA), was used to evaluate possible association of SIV env V1–V2 variants detected by HMA in plasma of monkeys one week after oral SIVmac251 inoculation with the levels of plasma viral RNA at this same time point. To determine potential linear associations of Entropy values or MMS values over time and SIV neutralizing antibody levels over time, linear regression was performed using Prism v.3.0 (GraphPad Software Inc., San Diego CA). For all statistical comparisons, a P value less than 0.05 was considered significant.
Kreiss J: Breastfeeding and vertical transmission of HIV-1. Acta Paediatr Suppl 1997, 421: 113-117.
Leroy V, Newell ML, Dabis F, Peckham C, Van de Perre P, Bulterys M, Kind C, Simonds RJ, Wiktor S, Msellati P: International mulicentre pooled analysis of late postnatal mother-to-child transmission of HIV-1 infection. Ghent International Working Group on Mother-to-Child Transmission of HIV. Lancet 1998, 352: 597-600. 10.1016/S0140-6736(98)01419-6
Nduati R, John G, Mbori-Ngacha D, Richardson B, Overbaugh J, Mwatha A, Ndinya-Achola J, Bwayo J, Onyango FE, Hughes J, Kreiss J: Effect of breastfeeding and formula feeding on transmission of HIV-1. JAMA 2000, 283: 1167-1174. 10.1001/jama.283.9.1167
Bertolli J, St Louis ME, Simonds RJ, Nieburg P, Kamenga M, Brown C, Tarande M, Quinn T, Ou CY: Estimating the timing of mother-to-child transmission of human immunodeficiency virus in a breast-feeding population in Kinshasa, Zaire. J Infect Dis 1996, 174: 722-726.
De Cock K, Fowler MG, Mercier E, de Vincenzi I, Saba J, Hoff E, Alnwick DJ, Rogers M, Shaffer N: Prevention of mother-to-child HIV transmission in resource-poor countries. Translating research into policy and practice. JAMA 2000, 283: 1175-1182. 10.1001/jama.283.9.1175
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: 2194-2203. 10.1128/JVI.75.5.2194-2203.2001
Essajee SM, Pollack H, Rochford G, Oransky I, Krasinski K, Borkowsky W: Early changes in quasispecies repertoire in HIV-infected infants: correlation with disease progression. AIDS Res Hum Retroviruses 2000, 16: 1949-57. 10.1089/088922200750054675
Delwart EL, Pan H, Sheppard HW, Wolpert D, Neumann AU, Korber B, Mullins JI: Slower evolution of human immunodeficiency virus type 1 quasispecies during progression to AIDS. J Virol 1997, 71: 7498-7508.
Shankarappa R, Margolick JB, Gange SJ, Rodrigo AG, Upchurch D, Farzadegan H, Gupta P, Rinaldo CR, Learn GH, He X, Huang XL, Mullins JI: Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection. J Virol 1999, 73: 10489-502.
Wolinsky SM, Carla MW, Korber BTM, Hutto C, Parks WP, Rosenblum LL, Kunstman KJ, Furtado MR, Munoz JL: Selective transmission of human immunodeficiency virus type 1 variant from mothers to infants. Science 1992, 255: 1134-1137.
Strunnikova N, Ray SC, Lancioni C, Nguyen M, Viscidi RP: Evolution of human immunodeficiency virus type 1 in relation to disease progression in children. J Hum Virol 1998, 1: 224-39.
Strunnikova N, Ray S, Livingston R, Rubalcaba E, Viscidi R: Convergent evolution within the V3 loop domain of human immunodeficiency virus type 1 in association with disease progression. J Virol 1995, 69: 7548-7558.
Hutto C, Zhou Y, He J, Geffin R, Hill M, Scott W, Wood C: Longitudinal studies of viral sequence, viral phenotype, and immunologic parameters of human immunodeficiency virus type 1 infection in perinatally infected twins with discordant diseases courses. J Virol 1996, 70: 3589-3598.
Salvatori F, Masiero S, Giaquinto C, Wade CM, Leigh-Brown AJ, Chieco-Bianchi L, De Rossi A: Evolution of human immunodeficiency virus type 1 in perinatally infected infants with rapid and slow progression to disease. J Virol 1997, 71: 4694-4706.
Ganeshan S, Dickover RE, Korber BTM, Bryson YJ, Wolinsky S: Human immunodeficiency virus type 1 genetic evolution in children with different rates of development of disease. J Virol 1997, 71: 663-677.
Halapi E, Leitner T, Jansson M, Scarlatti G, Orlandi P, Plebani A, Romiti L, Albert J, Wigzell H, Rossi P: Correlation between HIV sequence evolution, specific immune response and clinical outcome in vertically infected infants. AIDS 1997, 11: 1709-1717. 10.1097/00002030-199714000-00007
Geffin R, Hutto C, Andrew C, Scott GB: A longitudinal assessment of autologous neutralizing antibodies in children perinatally infected with human immunodeficiency virus type 1. Virology 2003, 310: 207-215.
Marthas ML, Van Rompay KKA, Otsyula M, Miller CJ, Canfield D, Pedersen NC, McChesney MB: Viral factors determine progression to AIDS in simian immunodeficiency virus-infected newborn rhesus macaques. J Virol 1995, 69: 4198-4205.
Van Rompay KKA, Berardi CJ, Dillard-Telm S, Tarara RP, Canfield DR, Valverde CR, Montefiori DC, Stefano Cole K, Montelaro RC, Miller CJ, Marthas ML: Passive immunization of newborn rhesus macaques prevents oral simian immunodeficiency virus infection. J Infect Dis 1998, 177: 1247-1259.
Van Rompay KKA, McChesney MB, Aguirre NL, Schmidt KA, Bischofberger N, Marthas ML: Two low doses of tenofovir protect newborn macaques against oral simian immunodeficiency virus infection. J Infect Dis 2001, 184: 429-438. 10.1086/322781
Van Rompay KKA, Greenier JL, Cole KS, Earl P, Moss B, Steckbeck JD, Pahar B, Rourke T, Montelaro RC, Canfield DR, Tarara RP, Miller CJ, McChesney MB, Marthas ML: Immunization of newborn rhesus macaques with simian immunodeficiency virus (SIV) vaccines prolongs survival after oral challenge with virulent SIVmac251. J Virol 2003, 77: 179-190. 10.1128/JVI.77.1.179-190.2003
Greenier JL, Miller CJ, Lu D, Dailey PJ, Lü FX, Kunstman KJ, Wolinsky SM, Marthas ML: Route of simian immunodeficiency virus inoculation determines the complexity but not the identity of viral variant populations that infect rhesus macaques. J Virol 2001, 75: 3753-3765. 10.1128/JVI.75.8.3753-3765.2001
Lamers SL, Sleasman JW, She JX, Barrie KA, Pomeroy SM, Barrett DJ, Goodenow MM: Persistence of multiple maternal genotypes of human immunodeficiency virus type I in infants infected by vertical transmission. J Clin Invest 1994, 93: 380-390.
Briant L, Wade CM, Puel J, Leigh Brown AJ, Guyader M: Analysis of envelope sequence variants suggests multiple mechanisms of mother-to-child transmission of human immunodeficiency virus type 1. J Virol 1995, 69: 3778-3788.
Wike CM, Korber BT, Daniels MR, Hutto C, Munoz J, Furtado M, Parks W, Saah A, Bulterys M, Kurawige JB, Wolinsky SM: HIV-1 sequence variation between isolates from mother-infant transmission pairs. AIDS Res Hum Retroviruses 1992, 8: 1297-1300.
Kliks S, Contag CH, Corliss H, Learn G, Rodrigo A, Wara D, Mullins JI, Levy JA: Genetic analysis of viral variants selected in transmission of human immunodeficiency viruses in newborns. AIDS Res Hum Retroviruses 2000, 16: 1223-1233. 10.1089/08892220050116998
Arroyo MA, Tien H, Pagán M, Swanstrom R, Hilyer GV, Cadilla CL, Meléndez-Guerrero LM: Virologic risk factors for vertical transmission of HIV type 1 in Puerto Rico. AIDS Res Hum Retroviruses 2002, 18: 447-460. 10.1089/088922202753614218
Pasquier C, Cayrou C, Blancher A, Tourne-Petheil C, Berrebi A, Tricoire J, Puel J, Izopet J: 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: 8493-8501.
Sutthent R, Foongladda S, Chearskul S, Wanprapa N, Likanonskul S, Kositanont U, Riengrojpitak S, Sahaphong S, Wasi C: V3 sequence diversity of HIV-1 subtype E in infected mothers and their infants. J Acquir Immune Defic Syndr Hum Retrovirol 1998, 18: 323-331.
Sutthent R, Foongladda S, Chearskul S, Wanaprapa N, Likanonskul S, Kositanont U, Riengrojpitak S, Sahapong S, Wasi C: Maternal and viral factors in vertical transmission of HIV-1 subtype E. Southeast Asian J Trop Med Public Health 1997, 28: 689-98.
Renjifo B, Chung M, Gilbert P, Mwakagile D, Msamanga G, Fawzi W, Essex M: In-utero transmission of quasispecies among human immunodeficiency virus type 1 genotypes. Virology 2003, 307: 278-282. 10.1016/S0042-6822(02)00066-1
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: 1001-1012.
Contag CH, Ehrnst A, Duda J, Bohlin AB, Lindgren S, Learn GH, Mullins JI: Mother-to-infant transmission of human immunodeficiency virus type 1 involving five envelope sequence subtypes. J Virol 1997, 71: 1292-1300.
Mulder-Kampinga G, Kuiken C, Dekker J, Scherpbier H, Boer K, Goudsmit J: Genomic human immunodeficiency virus type 1 RNA variation in mother and child following intra-uterine virus transmission. J Gen Virol 1993, 74: 1747-1756.
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: 1481-1489. 10.1089/088922200750006001
Amedee AM, Lacour N, Martin LN, Clements JE, Bohm RB Jr, Davison B, Harrison R, Murphey-Corb M: Genotypic analysis of infant macaques infected transplacentally and orally. J Med Primatol 1996, 25: 225-235.
Amedee AM, Lacour N, Ratterree M: Mother-to-infant transmission of SIV via breast-feeding in rhesus macaques. J Med Primatol 2003, 32: 187-193. 10.1034/j.1600-0684.2003.00024.x
Delwart E, Sheppard H, Walker B, Goudsmit J, Mullins J: Human immunodeficiency virus type 1 evolution in vivo tracked by DNA heteroduplex mobility assays. J Virol 1994, 68: 6672-6683.
Lukashov VV, Kuiken CL, Goudsmit J: Intrahost human immunodeficiency virus type 1 evolution is related to length of the immunocompetent period. J Virol 1995, 69: 6911-6916.
Mani I, Gilbert P, Sankalé J-L, Eisen G, Mboup S, Kanki PJ: Intrapatient diversity and its correlation with viral setpoint in human immunodeficiency virus type 1 CRF02_A/G-IbNG infection. J Virol 2002, 76: 10745-10755. 10.1128/JVI.76.21.10745-10755.2002
Markham RB, Wang W-C, Weisstein AE, Wang Z, Munoz A, Templeton A, Margolick J, Vlahov D, Quinn T, Farzadegan H, Yu X-F: Patterns of HIV-1 evolution in individuals with differing rates of CD4 T cell decline. Proc Natl Acad Sci USA 1998, 95: 12568-12573. 10.1073/pnas.95.21.12568
Wolinsky SM, Korber BTM, Neumann AU, Daniels M, Kunstman KJ, Whetsell AJ, Furtado MR, Cao Y, Ho DD, Safrit JT: Adaptive evolution of human immunodeficiency virus-type 1 during the natural course of infection. Science 1996, 272: 537-542.
Wolfs TF, de Jong JJ, Van den Berg H, Tijnagel JM, Krone WJ, Goudsmit J: Evolution of sequences encoding the principal neutralization epitope of human immunodeficiency virus 1 is host dependent, rapid, and continuous. Proc Natl Acad Sci USA 1990, 87: 9938-9942.
Liu SL, Schacker T, Musey L, Shriner D, McElrath MJ, Corey L, Mullins JI: Divergent patterns of progression to AIDS after infection from the same source: human immunodeficiency virus type 1 evolution and antiviral responses. J Virol 1997, 71: 4284-4295.
Rybarczyk BJ, Montefiori D, Johnson PR, West A, Johnston RE, Swanstrom R: Correlation between env V1/V2 region diversification and neutralizing antibodies during primary infection by simian immunodeficiency virus sm in rhesus macaques. J Virol 2004, 78: 3561-3571. 10.1128/JVI.78.7.3561-3571.2004
Van Rompay KKA, Abel K, Lawson JR, Singh RP, Schmidt KA, Evans T, Earl P, Harvey D, Franchini G, Tartaglia J, Montefiori D, Hattangadi S, Moss B, Marthas ML: Attenuated poxvirus-based SIV vaccines given in infancy partially protect infant and juvenile macaques against repeated oral challenge with virulent SIV. J Acquir Immune Defic Syndr 2005, 38: 124-134. 10.1097/00126334-200502010-00002
National Research Council: Guide for the care and use of laboratory animals. Washington, D. C.: National Academy Press; 1996.
Delwart E, Shaper E, McCutchan F, Louwagie J, Grez M, Rubsamen-Waigmann H, Mullins J: Genetic relationships determined by a heteroduplex mobility assay: analysis of HIV env genes. Science 1993, 262: 1257-1261.
Delwart EL, Herring B, Rodrigo AG, Mullins JI: Genetic subtyping of human immunodeficiency virus using a heteroduplex mobility assay. PCR Methods Appl 1995, 4: S202-S216.
Evans DT, Knapp LA, Jing P, Mitchen JL, Dykhuizen M, Montefiori DC, Pauza CD, Watkins DI: Rapid and slow progressors differ by a single MHC class I haplotype in a family of MHC-defined rhesus macques infected with SIV. Immunol Lett 1999, 66: 53-59. 10.1016/S0165-2478(98)00151-5
Knapp LA, Lehmann E, Piekarczyk MS, Urvater JA, Watkins DI: A high frequency of Mamu-1*01 in the rhesus macaque detected by polymerase chain reaction with sequence-specific primers and direct sequencing. Tissue Antigens 1997, 50: 657-661.
Montefiori DC, Baba TW, Li A, Bilska M, Ruprecht RM: Neutralizing and infection-enhancing antibody responses do not correlate with the differential pathogenicity of SIVmac239Δ3 in adult and infant rhesus monkeys. J Immunol 1996, 157: 5528-5535.
We thank D. Bennett, D. Brandt, I. Bolton, L. Brignolo, K. Christe, L. Hirst, A. Spinner, W. von Morgenland and the California National Primate Research Center Colony Services for expert technical assistance; M. Ma for assistance with image analysis of HMA gels; E. Delwart (Univ. of California, San Francisco) for useful discussions and suggestions. We thank Shilpa Hattangadi and Lynn Frampton for construction of recombinant MVAs. This work was supported by Public Health Science grant RR00169 from the National Center for Research Resources, NIH/NIAID grants AI39109 and AI46320 (MLM), and Elizabeth Glaser Scientist award #8-97 (MLM) from the Elizabeth Glaser Pediatric AIDS Foundation.
The author(s) declare they have no competing interests.
JG carried out the HMA studies and drafted the manuscript; KVR participated in the design and coordination of the study, acquisition and analysis of data, and helped draft the manuscript; DM analyzed and interpreted all neutralizing antibody assays; PE and BM designed and provided the MVA constructs, participated in the experimental design and manuscript writing; MM designed and coordinated the study, assisted in the data analyses and helped draft the manuscript.
About this article
Cite this article
Greenier, J.L., Van Rompay, K.K., Montefiori, D. et al. Simian immunodeficiency virus (SIV) envelope quasispecies transmission and evolution in infant rhesus macaques after oral challenge with uncloned SIVmac251: increased diversity is associated with neutralizing antibodies and improved survival in previously immunized animals. Virol J 2, 11 (2005). https://doi.org/10.1186/1743-422X-2-11