Evolution of subtype C HIV-1 Env in a slowly progressing Zambian infant
© Zhang et al; licensee BioMed Central Ltd. 2005
Received: 30 June 2005
Accepted: 07 November 2005
Published: 07 November 2005
Given the high prevalence of mother to child infection, the development of a better understanding of African subtype C HIV-1 transmission and natural evolution is of significant importance. In this study, we genotypically and phenotypically characterized subtype C viruses isolated over a 67-month follow-up period from an in utero-infected Zambian infant. Changes in genotype and phenotype were correlated to alterations of the host humoral immune response.
A comparison of baseline maternal and infant samples indicated that the infant sequences are monophyletic and contain a fraction of the diversity observed in the mother. This finding suggests that selective transmission occurred from mother to child. Peaks in infant HIV-1 Env genetic diversity and divergence were noted at 48 months, but were not correlated with changes in co-receptor usage or syncytia phenotype. Phylogenetic analyses revealed an accumulation of mutations over time, as well as the reappearance of ancestral lineages. In the infant C2-V4 region of Env, neither the median number of putative N-glycosylation sites or median sequence length showed consistent increases over time. The infant possessed neutralizing antibodies at birth, but these decreased in effectiveness or quantity with time. De novo humoral responses were detected in the child after 12 months, and corresponded with an increase in Env diversity.
Our study demonstrates a correlation between HIV-1 Env evolution and the humoral immune response. There was an increase in genetic diversification in the infant viral sequences after 12 months, which coincided with increases in neutralizing antibody titers. In addition, episodes of viral growth and successive immune reactions in the first 5–6 years were observed in this slow progressor infant with delayed onset of AIDS. Whether this pattern is typical of slow progressing subtype C HIV-1 infected infant needs to be further substantiated.
Subtype C human immunodeficiency virus type 1 (HIV-1) accounts for over 56% of HIV-1 infections [1–3]. Globally, HIV-1 infection is one of the leading causes of childhood morbidity and mortality. HIV-1 infected children account for 20% of all HIV-1 related deaths; 7% of individuals living with HIV-1 infection, and 16% of new HIV-1 infections annually . In sub-Saharan Africa, HIV-1 subtype C is responsible for approximately 50% of infections and a significant number of infections are in infants and children. Transmission of HIV-1 from infected mothers to their infants is the primary mode of HIV-1 infection in children and can occur in utero, intrapartum, or postnatally through breast milk. The use of antiretroviral regimens has successfully reduced the rate of HIV-1 infection in infants in the developed world to approximately 1%; nevertheless, such regimens have only recently become available in many of the developing nations where HIV-1 mother to child transmission (MTCT) is most significant .
HIV-1 MTCT is complex, and its determinants are not completely understood. Several factors, including high maternal viral load, maternal env gene homogeneity, and rapid viral replication kinetics, have been correlated with perinatal HIV-1 transmission [6–8]. In addition, advanced maternal disease status, lack of drug therapy, and lack of breast-feeding alternatives contribute to increased MTCT . Moreover, several studies have demonstrated the transmission of minor [9–12], major [9, 11], and multiple [9, 13, 14] HIV-1 genotypes from mother to infant. Our understanding of perinatal transmission and disease progression in infants is mainly derived from studies of subtype B infected individuals. The applicability of such findings to other subtypes remains to be substantiated.
The natural history of subtype C HIV-1 infection has not been extensively studied in children. It is known that infant disease survival times are considerably shorter than those of HIV-infected adults, and that without treatment, most HIV-1 infected African children die before their third birthday . Given the expanding distribution of subtype C infections, a complete understanding of virus transmission and natural evolution is increasingly important.
HIV-1 transmission is, in part, a function of the receptor binding by the envelope glycoprotein (Env) that mediates virus-cell fusion. Alteration of Env has been linked to expanded host range, alternative co-receptor usage and in vitro syncytium induction and associated with viral pathogenesis and disease progression [16–25]. Accumulating evidence suggests that subtype C Env displays biological properties, such as near-exclusive CCR-5 utilization, that distinguish it from other subtypes. In addition, the subtype C Env glycoprotein, third variable region (V3) is more conserved than the previously defined "constant" regions [26, 27]. Whether differences in cellular tropism, transmission and pathogenetic outcome observed between subtype C and other subtypes correlate with the Env glycoprotein biological or genetic properties need to be examined. In addition, whether there exist differences in Env evolution in infected children based on viral subtype, remains to be determined. Recently it has been suggested that particular changes in env in Zambian adults correlated with heterosexual transmission. Viruses with shorter Env length, and fewer putative N-linked glycosylation sites (PNGS) were suggested to be more susceptible to neutralizing antibodies, yet more efficient at transmission . Similar correlates have not been reported for transmission to children.
In the present study, we investigated the longitudinal variation of the viruses in a subtype C HIV-1 infected Zambian mother/infant pair (MIP 1157). This pair was antiretroviral therapy naïve over a six-year follow-up period. The extended follow-up enabled us to examine the interplay between humoral immune selection and virus evolution. We describe changes in the infant Env C2-V4 region over the follow-up period, and correlate these changes with alterations in viral phenotype and host humoral immune response. Our findings indicate that genetic diversification in the infant Env gene increased after 12 months, and is correlated with increases in neutralizing antibody titers.
HIV-1 infected mother-infant pair
We characterized HIV-1 transmission and longitudinal evolution of the HIV-1 envelope glycoprotein in a Zambian mother and infant pair (MIP 1157) for more than 6 years. The mother and child are anti-retroviral naïve and remain clinically asymptomatic. Infant 1157 was infected in utero since HIV-1 sequences were detected by DNA PCR of infant blood samples collected at birth. The baby was delivered naturally, healthy and with normal birth weight, and was breast-fed until 20 months of age. The child remains clinically asymptomatic throughout the follow-up study period and his CD4 counts was 658 cells/μl at 6 years old. The child has been evaluated at the study clinic where blood specimens were collected every 6 months for the first 24 months and at 12-month intervals thereafter. The prolonged survival of this infected child is unusual since most untreated HIV-1 infected African children do not survive beyond the first three years of life. The extended follow-up of infant1157 provided us with an opportunity to investigate correlates of virus transmission in the Env glycoprotein and to track genetic variation and evolution of this gene over time.
MIP1157 viruses use CCR5 as co-receptor and belong to subtype C
All viral isolates recovered from MIP 1157 replicated efficiently in PBMC and monocyte-derived macrophages (MDM), but failed to grow in MT-2 and C8166 T-cell lines. Viral isolates did not induce syncytia in infected PBMC and MDM. We evaluated viral co-receptor usage in cell lines that co-express CD4 with a single co-receptor. All isolates failed to infect CXCR4-expressing CEMx174-GFP cells, and similarly, none of the viruses grew in cells expressing only CCR3 (HOS-CD4-CCR3) (data not shown). In addition, 1157 viruses failed to replicate in PBMC homozygous for the Δ32 deletion variant of CCR5. In contrast, cells expressing normal CCR5 and CD4 (GHOST-CD4-CCR5) were readily infected, suggesting that 1157 HIV-1 isolates primarily use CCR5 as co-receptor (data not shown). This is in agreement with infectivity assays demonstrating that only primary PBMC and MDM support viral growth. Phylogenetic analyses clustered all 1157 env sequences with subtype C.
Viral variations in the different mother and infant populations
Infant at birth
2 (0 – 7)
1 (0 – 5)
15 (14 – 15)
183 (183 – 183)
Infant 6 months
4 (0 – 11)
3 (0 – 9)
14 (13 – 15)
183 (175 – 183)
Infant 12 months
6 (0 – 15)
4 (0 – 11)
13 (13 – 15)
174 (174 – 183)
Infant 18 months
15 (0 – 29)
11 (0 – 23)
14 (13 – 15)
179 (175 – 183)
Infant 24 months
14 (1 – 21)
11 (0 – 18)
12 (11 – 16)
177 (174 – 183)
Infant 29 months
16 (1 – 26)
12 (0 – 19)
12.5 (11 – 15)
182 (173 – 183)
Infant 36 months
13 (0 – 27)
8 (0 – 19)
12 (10 – 14)
176 (173 – 183)
Infant 48 months
25 (2 – 37)
16 (2 – 26)
13 (12 – 14)
182 (176 – 185)
Infant 67 months
35 (0 – 57)
24 (0 – 37)
13 (12 – 15)
183 (180 – 185)
Mother at delivery
5 (1 – 10)
3 (0 – 7)
15 (14 – 15)
183 (183 – 183)
Mother 12 months
8 (1 – 23)
5 (0 – 13)
15 (13 – 15)
183 (183 – 183)
Mother 18 months
5 (0 – 13)
2 (0 – 9)
15 (13 – 15)
183 (183 – 183)
Mother 24 months
2 (0 – 7)
1 (0 – 5)
14 (13 – 14)
183 (183 – 183)
Longitudinal variation in envsequences
A recent report suggested that subtype C viruses transmitted between members of Zambian discordant couples possess envelope glycoproteins that are under-glycosylated, neutralization sensitive and contain short loop structures . To explore the potential role of specific sequence characteristics in virus transmission between mother and child, we compared the sequence length polymorphism and variation in the number of PNGS for baseline maternal and infant Env C2-V4 sequences. There are 15 PNGS in this region of 1157 Env. Maternal and infant baseline sequences are all of the same length, and showed little variation in the PNGS (Table 1 and Figure 4). In the mother, there were 4 sequences, out of 26, that lost a PNGS, and the position at which this site was ablated was not conserved among any of the four. In the infant, 6 of 48 sequences lost a single PNGS, but in parallel with the mother, there was no conservation in the position of that loss. Moreover, only one variable PNGS was shared between the mother and the infant.
A similar evaluation of the C2-V4 length polymorphism and PNGS alteration was carried out on subsequent infant samples to assess the longitudinal variation in these two parameters (Table 1 and Figure 4). Length polymorphism was only observed in infant sequences where putative insertions and/or deletions occur in a subset of sequences at amino-acid positions 106–109 and 166–180. Maternal sequences remained of constant length, 183 amino acids, throughout the follow-up. All transmitted sequences in the infant were initially of the same length (also 183 amino acids). Length polymorphism in the region spanning amino acids 166–180 appeared 6 months postpartum, whereas polymorphism in the region spanning 106–109 was first observed at 12 months. The longest sequences, isolated at 48 and 67 months, were 185 amino acids in length, whereas the shortest sequences, 173 amino acids, were isolated at 29 and 36 months. All infant PNGS present at baseline remain present in a fraction of sequences from subsequent timepoints; however, only 3 sites remained fixed over the entire course of infection. The largest PNGS variation was observed at positions 7, 104, and 177, which oscillate between high and low prevalence (Figure 4, months 24, 36 and 48; position 7). In addition, there are 2 sites gained, one at position 72 at 6 months, and the other at position 109 at 12 months. Both of these polymorphisms are low in frequency and both are adjacent to another PNGS.
Neutralization of infant HIV-1 isolates
Longitudinal evolution of HIV-1 subtype C has rarely been evaluated in infected children. The survival of infant 1157 for more than 6 years post-infection provided us with an opportunity to track genetic variation and phenotypic evolution in the viral envelope glycoprotein over that period. In addition, we were able to examine correlations between these viral properties and the humoral immune response of the child. Detection of HIV-1 sequences in PBMC collected from the child at birth indicated in utero infection. The pattern of genetic variation shown by phylogenetic analysis at baseline is compatible with an episode of selective transmission, as reported in previous studies [9–12, 27, 29]. In utero infection of infants has been reported to result in more rapid disease progression [30–32]; however, the extended survival of infant 1157 suggests the route of infection alone is not predictive of disease progression in subtype C infected children.
HIV-1 has replication and mutation rates that generate high numbers of progeny and significant genetic variation. The env gene has been calculated to diverge at a rate of about 1% per year . The patterns of HIV-1 evolution in infected individuals, even for subtype B viruses, are ambiguous. Delwart et al. reported several-fold higher diversity at the early stage versus the late stage of infection . In contrast, other studies have shown that viral sequences in env are more homogenous early in infection and diversify with disease progression and decline in CD4+ T cell counts [33, 35–40]. Here we show that birth env sequences in the recipient child were highly homogenous, as indicated by env diversity, and were closely related to, but encompassed only a subset of the contemporary maternal variation. Genetic analysis at multiple timepoints showed that diversity in env as well as divergence from the initial infecting species increased with time up to 48 months. This increase in diversity and divergence correlated with parallel increases in non-synonymous changes. Whether such an increase is unique to this case needs to be further substantiated. Our findings contrast with those from studies of subtype B infected adults where, in patients infected with viruses that undergo co-receptor switching, the peak of diversity correlated with the development of CXCR4 utilization and the peak of divergence correlated with the maximal prevalence of CXCR4 utilizing species . These phenomena are not relevant to 1157 since no alternative co-receptor usage was detected in either the mother or the child. Although X4-utilizing subtype C strains have been described [41–44], they are unusual, thus pointing to distinct evolutionary pressures on the various subtypes. Since subtype C infected individuals possess X4-expressing cells, it is likely that immunological and viral replicative selection in these individuals do not force or allow subtype C to efficiently utilize these targets or other constraints make such utilization significantly unfavorable.
Interestingly, we have observed the apparent reappearance of earlier lineages at the 67-month time point, and this is probably correlated to the decrease in viral genetic divergence at the same time point. Our observation would indicate that viral sequences, presumably emerging from latently infected cells, can reintroduce ancestral lineages and thus could lead to the decrease in divergence. It is tempting to speculate that such reintroduction might coincide with the waning of the immune response to these 'earlier' viruses in much the same way as antiretroviral therapy interruption often results in repopulation of the patient with drug-sensitive ancestral strains. Alternatively, the host environment may have altered in such a fashion that an ancestral variant becomes more viable due to higher replication fitness and decay of immune selection.
Our sequence analysis also revealed a substantial amount of variation (mutation, deletion and insertion) in envC3 and V4 regions in infant samples, implying that C3 or V4 domain is a likely target of immunological or replicative selective pressure during subtype C virus evolution and disease progression in children. The significance of C3 and V4 variation is currently under investigation. It is important to recognize that definitions of the constant and variable domains in Env are derived primarily from studies of subtype B viruses, and the patterns of sequence diversity in those isolates may not be reflected in other subtypes such as subtype C.
Our neutralization assays support the concept that the humoral immune response developed in parallel with the evolving HIV-1 envelope sequences and constitutes part of the selective pressure on the gene [45, 46]. The persistence of high level neutralizing antibodies against early infant viral isolates indicated that the infant immune system is capable of developing and maintaining strong responses to eliminate the initially transmitted and replicating virus (Figure 6A). It has been shown that neutralization escape mutants with reduced sensitivity to autologous sera emerge rapidly in HIV-1 infected adults [46–48], but patients subsequently developed additional neutralizing antibodies to the 'escape' viruses after a delay . The initial effectiveness of the infant sera is likely due to a significant contribution by maternal antibodies to neutralization titer. Nevertheless, the child does not receive the full repertoire of maternal neutralizing antibody since a disparity was observed between the effectiveness of maternal and infant baseline neutralization titers. This idea is reinforced by the fact that the maternal baseline serum continues to be effective against the infant viruses for the duration of infection; whereas the ability of the infant serum to neutralize contemporary viruses is reduced after the early timepoints. Moreover, differences in the susceptibility of viral isolates to be neutralized by antibodies was independent of the replication rates, since the 6, 12 and 48-month viral isolates replicated with nearly identical kinetics.
The observed viral diversity increase at 12-months might coincide with the diminution of maternal antibody effectiveness. However, the increasing titer of antibodies beyond 12 month implied the development of de novo infant humoral immune responses against the diversifying population. This response, as might be anticipated, is always in reaction to the viral alterations, not in anticipation of it. This conclusion is supported by the finding that despite an apparent failure of the humoral immunity to control HIV-1 replication through neutralizing antibodies at 48 months, infant 1157 mounted an effective neutralizing response to that virus at subsequent timepoints (67 month) (Figure 6C) and this coincided with a decrease in viral diversity (Figure 2). However, the role of cell-mediated immunity in controlling viral replication cannot be determined for this infant since viable cells were not available.
It has been suggested that a more antigenically diverse virus population would correlate to a broader immune reactivity, a slower rate of disease progression [50, 51], and selection of neutralization escape mutants in HIV-1 infected individuals, including long-term non-progressors [47, 52–54]. Our study, even though with only one mother infant pair, appears to support this hypothesis but further analysis involving a larger number of patients, including rapid and slow progressors, followed longitudinally will be needed to substantiate this observation. A more complete understanding of the mechanisms of humoral immune escape with a more precise definition of the regions in Env where such mutations cluster is likely to impact vaccine design.
It has recently been observed that viruses with shorter V1-V4 Env length, and fewer glycans are more susceptible to neutralizing antibodies, but mediate more efficient transmission in discordant couples . Assuming this concept, one would expect to see a relative lengthening of Env, and an increase in the number of glycans with time. Our analysis of MIP 1157 longitudinally, which was based on C2-V4 sequences, cannot be used for direct comparison for transmission, we did, however, observe increases in variation at PNGS and in sequence length over time. The variation in Env over the follow-up period frequently resulted in the deletion, addition, or relocation of potential N-glycans, suggesting a role of N-glycans for immune selection in the HIV-1 evolution. The hot spots of N-glycan variation were particularly evident in the C2 and C3 regions. Similar changes in potential glycosylation sites have been hypothesized to modify a "glycan shield" for evading neutralizing antibodies .
We have demonstrated that genetic diversification in the infant sequences increased after 12 months, and this coincided with increases in neutralizing antibody titers. In addition, episodes of viral growth and successive immune reactions in the first 5–6 years were observed in this slow progressor infant with delayed onset of AIDS. Longitudinal studies such as the one described here underscore the dynamic and complex interactions of viral populations and immune responses. Whether this pattern of viral host interaction is typical of slow progressing infected infant needs to be further substantiated.
Patient population and sample collection
The mother-infant pair (MIP) 1157 characterized in this study was recruited to investigate the routes of transmission of HIV-1. Venous blood was obtained from the mother before delivery and from the infant within 24 hours of delivery. Follow-up blood specimens were obtained when the pair returned for visits at 6, 12, 18, 24, 29, 36, 48 and 67-months after delivery. The HIV-1 serological status of the mother was determined by two rapid assays, Capillus (Cambridge Biotech, Ireland) and Determine (Abbott laboratories, USA), on the initial blood samples. The positive serological result was confirmed by immunofluorescence assay (IFA), as previously described . The status of HIV-1 infection in the infant was determined by performing viral isolation from the infant's peripheral blood mononuclear cells (PBMC) and by PCR on DNA isolated on the day of birth.
HIV-1 was isolated sequentially over a 67-month post-delivery period by standard co-culture procedures. Donor HIV-1-negative PBMC were purified using Lymphoprep (Life Technology). The purified lymphocytes were then propagated in RPMI 1640 medium containing 10% heat-inactivated fetal bovine serum (FBS) and 5 μg/ml of phytohemagglutinin (Sigma) for 40 h before co-culturing with MIP 1157 PBMC or whole blood at a combined final concentration of 2 × 106cells /ml. Equal numbers of fresh uninfected PHA-stimulated PBMC were added to the culture weekly. Virus production was monitored using a commercial ELISA to measure HIV-1 p24 antigen levels (Coulter immunology, FL). Virus stocks were prepared when p24 antigen concentration exceeded 10 ng /ml (about 7–10 days). Viral isolates were recovered from 6-month maternal and 6, 12, 18, 24, 29, and 48-month infant samples
Phenotype, syncytium-inducing (SI) or non-syncytium-inducing (NSI), was determined by infecting MT-2 cells in a 12-well tissue culture plate (5 × 105 cells / well) with 5 ng of p24 virus stock per well. Cell cultures were observed daily for syncytia formation, over a course of 10 days. Levels of p24 antigen were determined in supernatants collected on day 2, 4, 7, and 10 post-infection. Virus was scored as SI if syncytia formation and increasing level of p24 antigen were observed within the 10-day period, and as NSI if syncytia failed to form within that time.
To define the viral tropism, primary monocyte-derived macrophages (MDM), and MT-2 or C8166 T-cell lines were infected with the virus stocks using standard methods. Primary monocytes were obtained from gradient-purified PBMC by adherence to plastic culture dishes . Adherent cells were cultured for 7 to 10 days in RPMI 1640 medium containing 10% FBS and 10 ng/ml of granulocyte-macrophage, colony-stimulating factor (GIBCO) to promote differentiation of monocytes to macrophages. Differentiated macrophages, or T-cell lines, were infected with 5 to 10 ng of HIV-1 p24 antigen per 5 × 105 cells and incubated for 4 to 5 h at 37°C. Subsequently, the infected cells were washed twice with phosphate buffered saline (PBS) and resuspended in fresh culture medium. Culture supernatants were removed at 3, 7, and 14 days post-infection and assayed for HIV-1 p24 antigen. A culture well was considered virus-positive if increasing level of p24 antigen was observed.
Chemokine co-receptor usage
Determination of co-receptor usage was carried out using cell lines obtained through the NIH AIDS Research and Reference Reagent Program, Division of AIDS, NIAID, NIH from Dr. Nathaniel Landau that express specific co-receptors (CEMx174-GFP cells [CXCR4], Ghost-CCR5 cells [CCR5] and HOS-CD4-CCR3 cells [CCR3]). PBMC from an individual homozygous for CCR5 mutation Δ32 were obtained from Dr. James Hoxie (University of Pennsylvania). To test for co-receptor usage, the CCR5-Δ32 PBMC and the three co-receptor-specific cell lines were seeded at a density of 1 × 106 cells/ml into 24-well culture plates. The cells were infected with 5 ng /ml of HIV-1 p24. The infected CEMx174-GFP and Ghost-CCR-5 cells were observed microscopically on day 2–3 post-infection for green fluorescent protein (GFP) expression. Wells exhibiting a count of GFP-expressing cells greater than or equal to 3-fold the negative control wells were scored as positive. Uninfected control wells produced only one to two GFP expressing cells per well. HIV-1 strains SF2 and NL4-3 were used as positive controls for viruses that use CXCR4, and HIV-1 strain SF128A was used as positive control for CCR5 utilization. Positive control viruses consistently gave 7-fold, or greater, GFP- expressing cells than the background control.
Infection of CCR5-Δ32 PBMC and HOS-CD4-CCR3 cells was monitored by measuring HIV-1 p24 antigen production in culture supernatants. HIV-1 p24 was measured at 3 days post-infection for HOS-CD4-CCR3 cells, and at 3, 7 and 10 days post-infection for the CCR5-Δ32 PBMC. Cultures were considered positive for viral growth if more than 100 pg/ml of p24 was detected.
Viral isolates replication kinetics
Equivalent infectious units, 100 TCID50, of the infant viral isolates obtained at 6, 12 and 48-month follow up were added to triplicate wells in a 12-well plate containing 6 × 106 PHA-stimulated PBMC from a HIV-1 seronegative blood donor. The laboratory isolate 128A was used as a replication kinetics control. After incubation at 37°C for 6 hours, cells were washed 3 times with PBS and refilled with fresh medium. All infected cultures were sampled and supplemented with a 50% volume of fresh culture medium at day 1, 3, 5, 7, 9, 11, 13 and 15. Viral replication kinetics in PBMC was determined by measuring RT units in culture supernatants at day 0, 1, 3, 5, 7, 9, 11, 13 and 15-postinfection. Viruses were lysed using 10% Triton X-100 (1% final concentration) in RPMI medium supplemented with 10% FBS, and RT units was measured using EnzChek Reverse Transcriptase Assay Kit (Invitrogen, Eugene, Oregon). The assay was performed in triplicate.
Virus neutralization assay
Plasma neutralization activity was determined through infections of TZM-bl cells (NIH AIDS Research and Reference Reagent Program catalogy no. 8129, TZM-bl) as described in Wei et al (2003) with modifications. TZM-bl cells stably express high levels of CD4, CCR5 and CXCR4. The cells contain HIV-1 LTR promoter cassettes that express luciferase and β-galactosidase in response to stimulation with HIV-1 Tat. TZM-bl cells were plated at a density of 6 × 103/well in 96-well tissue culture plates (Falcon) and cultured overnight in DMEM supplemented with 10% FBS. Test plasma was heat-inactivated at 56°C for 30 min, spun at 3,000 × g for 5 min and diluted 1:20, 1:100 and 1:500 in DMEM plus 6% FBS. Viral aliquots of 100 TCID50/ml were prepared in DMEM supplemented with 6% FBS and 80 μg/ml DEAE dextran, to a combined total volume of 100 μl. The virus aliquots (100 μl) were combined with 100 μl of the different test plasma dilutions and the mixture was incubated for 1h at 37°C. Following incubation, the virus-plasma mixture was added to TZM- bl cells and incubated at 37°C for two days. Following two washes with PBS, the level of virus infection was measured by luciferase activity. Cells were lysed using Luciferase Assay Reagent (Promega, Madison, WI) and the luciferase activity was measured using a LUMIstar luminometer (BMG Lab Technologies, Offenburg, Germany). The assay was performed in triplicate. Controls included cells infected by virus inoculated with medium or normal human plasma instead of the test plasma. Effective neutralization by the plasma would reduce the level of luciferase versus controls lacking plasma. The luciferase activity in the control wells without plasma was defined as 100 %, and the neutralization titer of the test plasma was calculated relative to this value.
Polymerase chain reaction, gene cloning, sequencing and subtype identification
Genomic DNA was purified from patient PBMC using the ISOQUICK kit (ORCA Research, Inc.). Primers used for amplification of the subtype C env gene were designed based on a reference alignment of all HIV-1 subtypes obtained from the Los Alamos HIV-1 Sequence Database http://hiv-web.lanl.gov. Nested PCR, as previously described , was used to amplify a 617 bp fragment containing the C2-V4 region of the maternal and infant env gene from samples collected at birth through 67-month. The amplified fragments were cloned into the pGEM-T Easy vector (Promega) and sequenced in both directions with dideoxy terminators (ABI BigDye Kit).
Sequence alignment and analyses
Nucleotide sequences were translated and the amino acid sequences aligned using Clustal W, as implemented in BioEdit (Hall 1999), and further refined manually, preserving the reading frame.
Bayesian (BA), maximum likelihood (ML), and neighbor-joining (NJ) phylogenetic analyses were implemented to: 1) determine subtype affiliation (aligning infant clones to a reference panel from the HIV database at Los Alamos National Laboratory http://www.hiv.lanl.gov; 2) assess the transmission pattern between 1157 mother and infant (samples collected at birth from the mother and the infant were evaluated), 3) visualize how viral populations change with time (all samples collected from mother and infant were included). Bayesian searches were run in MrBayes version 3.04  for 5 × 106 generations. Trees were sampled every 1000 generations. Convergence was reached after 1.5 × 106, and a majority rule consensus was built based on the last 3000 trees. ML searches were implemented in Treefinder  using a general time-reversible model of nucleotide substitution for each codon position. NJ analyses were performed in PAUP* ver4.10b , estimating distances by maximum likelihood, after selecting the best fitting model of nucleotide substitution in Modeltest ver 3.5 . Support for the nodes in ML and NJ was evaluated by running 1000 pseudoreplicates in bootstrap analyses.
Genetic diversity, temporal structure and estimates of selective strength
Samples were grouped into contemporaneous sets according to the time of collection. Sites with alignment gaps, corresponding to putative insertions or deletions, were excluded from comparisons. Variation in the pattern of genetic diversity for each time-point was explored using MEGA3,  and DNAsp ver 4.06 , Datamonkey  and PAML ver 3.13. , and the number and location of putative glycosylation sites (PNGS) were estimated using N-GlycoSite from the Los Alamos National Laboratory database http://hiv-web.lanl.gov. Viral diversity and viral divergence from the initial population were calculated for each timepoint. Viral diversity estimates were based on π, the average nucleotide differences between sequences per site. Changes in viral diversity were used to estimate relative changes in the effective viral population size, assuming a similar mutation rate at each timepoint. Viral divergence was calculated as the average genetic distance to the viral population in the infant at birth.
The relative strength of positive selection was evaluated in a ML framework using PAML ver 3.13. . To estimate the overall dN/dS for each set of contemporaneous sequences we implemented the one-rate codon model described by Yang and Nielsen (1998) . The number of synonymous and non-synonymous substitutions observed per site in the infant Env domains (C2-V4) was calculated in Datamonkey so as to visualize the distribution of synonymous and non-synonymous substitutions along the sequence and to identify regions that accumulate most variation ('hot spots').
This study was supported by PHS grants HD39620, CA76958, PO1 AI 48240, Fogarty Training grant TW01429, and NCRR COBRE grant RR15635 to C.W, PO1 AI34266 to R.M.R. We thank Dr. James Hoxie and for providing us with CCR5 Δ 32 cells; Dr. Nathaniel Landau for HOS-CD4-CCR3 cells through the NIH AIDS research and reference reagent Program, Division of AIDS.
- Essex M: State of the HIV pandemic. J Hum Virol. 1998, 1 (7): 427-429.PubMedGoogle Scholar
- Janssens W, Buve A, Nkengasong JN: The puzzle of HIV-1 subtypes in Africa. AIDS. 1997, 11 (6): 705-712. 10.1097/00002030-199706000-00002.View ArticlePubMedGoogle Scholar
- Rodenburg CM, Li Y, Trask SA, Chen Y, Decker J, Robertson DL, Kalish ML, Shaw GM, Allen S, Hahn BH, Gao F: Near full-length clones and reference sequences for subtype C isolates of HIV type 1 from three different continents. AIDS Res Hum Retroviruses. 2001, 17 (2): 161-168. 10.1089/08892220150217247.View ArticlePubMedGoogle Scholar
- Luzuriaga K, Sullivan JL: Pediatric HIV-1 infection: advances and remaining challenges. AIDS Rev. 2002, 4 (1): 21-26.PubMedGoogle Scholar
- Edgeworth RL, Ugen KE: Immunopathological factors for vertical transmission of HIV-1. Pathobiology. 2000, 68 (2): 53-67. 10.1159/000028115.View ArticlePubMedGoogle Scholar
- Arroyo MA, Tien H, Pagan M, Swanstrom R, Hillyer GV, Cadilla CL, Melendez-Guerrero LM: Virologic risk factors for vertical transmission of HIV type 1 in Puerto Rico. AIDS Res Hum Retroviruses. 2002, 18 (6): 447-460. 10.1089/088922202753614218.View ArticlePubMedGoogle Scholar
- Garcia PM, Kalish LA, Pitt J, Minkoff H, Quinn TC, Burchett SK, Kornegay J, Jackson B, Moye J, Hanson C, Zorrilla C, Lew JF: Maternal levels of plasma human immunodeficiency virus type 1 RNA and the risk of perinatal transmission. Women and Infants Transmission Study Group. N Engl J Med. 1999, 341 (6): 394-402. 10.1056/NEJM199908053410602.View ArticlePubMedGoogle Scholar
- John GC, Nduati RW, Mbori-Ngacha DA, Richardson BA, Panteleeff D, Mwatha A, Overbaugh J, Bwayo J, Ndinya-Achola JO, Kreiss JK: Correlates of mother-to-child human immunodeficiency virus type 1 (HIV-1) transmission: association with maternal plasma HIV-1 RNA load, genital HIV-1 DNA shedding, and breast infections. J Infect Dis. 2001, 183 (2): 206-212. 10.1086/317918.View ArticlePubMedGoogle Scholar
- 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-2203. 10.1128/JVI.75.5.2194-2203.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Wolinsky SM, Wike CM, Korber BT, Hutto C, Parks WP, Rosenblum LL, Kunstman KJ, Furtado MR, Munoz JL: Selective transmission of human immunodeficiency virus type-1 variants from mothers to infants. Science. 1992, 255 (5048): 1134-1137.View ArticlePubMedGoogle Scholar
- Scarlatti G, Leitner T, Halapi E, Wahlberg J, Marchisio P, Clerici-Schoeller MA, Wigzell H, Fenyo EM, Albert J, Uhlen M, Rossi P: 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-1725.PubMed CentralView ArticlePubMedGoogle Scholar
- 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-1012.PubMed CentralPubMedGoogle Scholar
- 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 (11): 8493-8501.PubMed CentralPubMedGoogle Scholar
- 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 (1): 380-390.PubMed CentralView ArticlePubMedGoogle Scholar
- Chakraborty R: HIV-1 infection in children: a clinical and immunologic overview. Curr HIV Res. 2005, 3 (1): 31-41. 10.2174/1570162052773022.View ArticlePubMedGoogle Scholar
- Hwang SS, Boyle TJ, Lyerly HK, Cullen BR: Identification of the envelope V3 loop as the primary determinant of cell tropism in HIV-1. Science. 1991, 253 (5015): 71-74.View ArticlePubMedGoogle Scholar
- Fouchier RA, Groenink M, Kootstra NA, Tersmette M, Huisman HG, Miedema F, Schuitemaker H: Phenotype-associated sequence variation in the third variable domain of the human immunodeficiency virus type 1 gp120 molecule. J Virol. 1992, 66 (5): 3183-3187.PubMed CentralPubMedGoogle Scholar
- de Jong JJ, Goudsmit J, Keulen W, Klaver B, Krone W, Tersmette M, de Ronde A: Human immunodeficiency virus type 1 clones chimeric for the envelope V3 domain differ in syncytium formation and replication capacity. J Virol. 1992, 66 (2): 757-765.PubMed CentralPubMedGoogle Scholar
- Milich L, Margolin B, Swanstrom R: V3 loop of the human immunodeficiency virus type 1 Env protein: interpreting sequence variability. J Virol. 1993, 67 (9): 5623-5634.PubMed CentralPubMedGoogle Scholar
- Stamatatos L, Cheng-Mayer C: Evidence that the structural conformation of envelope gp120 affects human immunodeficiency virus type 1 infectivity, host range, and syncytium-forming ability. J Virol. 1993, 67 (9): 5635-5639.PubMed CentralPubMedGoogle Scholar
- Richman DD, Bozzette SA: The impact of the syncytium-inducing phenotype of human immunodeficiency virus on disease progression. J Infect Dis. 1994, 169 (5): 968-974.View ArticlePubMedGoogle Scholar
- Cocchi F, DeVico AL, Garzino-Demo A, Cara A, Gallo RC, Lusso P: The V3 domain of the HIV-1 gp120 envelope glycoprotein is critical for chemokine-mediated blockade of infection. Nat Med. 1996, 2 (11): 1244-1247. 10.1038/nm1196-1244.View ArticlePubMedGoogle Scholar
- Connor RI, Sheridan KE, Ceradini D, Choe S, Landau NR: Change in coreceptor use coreceptor use correlates with disease progression in HIV-1 – infected individuals. J Exp Med. 1997, 185 (4): 621-628. 10.1084/jem.185.4.621.PubMed CentralView ArticlePubMedGoogle Scholar
- Hayman A, Moss T, Arnold C, Naylor-Adamson L, Balfe P: Disease progression in heterosexual patients infected with closely related subtype B strains of HIV type 1 with differing coreceptor usage properties. AIDS Res Hum Retroviruses. 2004, 20 (4): 365-371. 10.1089/088922204323048113.View ArticlePubMedGoogle Scholar
- Nabatov AA, Pollakis G, Linnemann T, Kliphius A, Chalaby MI, Paxton WA: Intrapatient alterations in the human immunodeficiency virus type 1 gp120 V1V2 and V3 regions differentially modulate coreceptor usage, virus inhibition by CC/CXC chemokines, soluble CD4, and the b12 and 2G12 monoclonal antibodies. J Virol. 2004, 78 (1): 524-530. 10.1128/JVI.78.1.524-530.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Oelrichs RB, Shrestha IL, Anderson DA, Deacon NJ: The explosive human immunodeficiency virus type 1 epidemic among injecting drug users of Kathmandu, Nepal, is caused by a subtype C virus of restricted genetic diversity. J Virol. 2000, 74 (3): 1149-1157. 10.1128/JVI.74.3.1149-1157.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang H, Orti G, Du Q, He J, Kankasa C, Bhat G, Wood C: Phylogenetic and phenotypic analysis of HIV type 1 env gp120 in cases of subtype C mother-to-child transmission. AIDS Res Hum Retroviruses. 2002, 18 (18): 1415-1423. 10.1089/088922202320935492.View ArticlePubMedGoogle Scholar
- Derdeyn CA, Decker JM, Bibollet-Ruche F, Mokili JL, Muldoon M, Denham SA, Heil ML, Kasolo F, Musonda R, Hahn BH, Shaw GM, Korber BT, Allen S, Hunter E: Envelope-constrained neutralization-sensitive HIV-1 after heterosexual transmission. Science. 2004, 303 (5666): 2019-2022. 10.1126/science.1093137.View ArticlePubMedGoogle Scholar
- Verhofstede C, Demecheleer E, De Cabooter N, Gaillard P, Mwanyumba F, Claeys P, Chohan V, Mandaliya K, Temmerman M, Plum J: 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-3057. 10.1128/JVI.77.5.3050-3057.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Blanche S, Rouzioux C, Moscato ML, Veber F, Mayaux MJ, Jacomet C, Tricoire J, Deville A, Vial M, Firtion G, Crepy AD, Douard D, Robin M, Courptin C, Ciraru-Vigneron N, Diest F, Griscelli C: A prospective study of infants born to women seropositive for human immunodeficiency virus type 1. HIV Infection in Newborns French Collaborative Study Group. N Engl J Med. 1989, 320 (25): 1643-1648.View ArticlePubMedGoogle Scholar
- Blanche S, Tardieu M, Duliege A, Rouzioux C, Le Deist F, Fukunaga K, Caniglia M, Jacomet C, Messiah A, Griscelli C: Longitudinal study of 94 symptomatic infants with perinatally acquired human immunodeficiency virus infection. Evidence for a bimodal expression of clinical and biological symptoms. Am J Dis Child. 1990, 144 (11): 1210-1215.View ArticlePubMedGoogle Scholar
- Mayaux MJ, Burgard M, Teglas JP, Cottalorda J, Krivine A, Simon F, Puel J, Tamalet C, Dormont D, Masquelier B, Doussin A, Rouzioux C, Blanche S: Neonatal characteristics in rapidly progressive perinatally acquired HIV-1 disease. The French Pediatric HIV Infection Study Group. Jama. 1996, 275 (8): 606-610. 10.1001/jama.275.8.606.View ArticlePubMedGoogle Scholar
- 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 (12): 10489-10502.PubMed CentralPubMedGoogle Scholar
- 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 (10): 7498-7508.PubMed CentralPubMedGoogle Scholar
- 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 disease courses. J Virol. 1996, 70 (6): 3589-3598.PubMed CentralPubMedGoogle Scholar
- Markham RB, Wang WC, Weisstein AE, Wang Z, Munoz A, Templeton A, Margolick J, Vlahov D, Quinn T, Farzadegan H, Yu XF: Patterns of HIV-1 evolution in individuals with differing rates of CD4 T cell decline. Proc Natl Acad Sci USA. 1998, 95 (21): 12568-12573. 10.1073/pnas.95.21.12568.PubMed CentralView ArticlePubMedGoogle Scholar
- Nowak MA, Bangham CR: Population dynamics of immune responses to persistent viruses. Science. 1996, 272 (5258): 74-79.View ArticlePubMedGoogle Scholar
- Salvatori F, Masiero S, Giaquinto C, Wade CM, 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 (6): 4694-4706.PubMed CentralPubMedGoogle Scholar
- Joos B, Trkola A, Fischer M, Kuster H, Rusert P, Leemann C, Boni J, Oxenius A, Price DA, Phillips RE, Wong JK, Hirschel B, Weber R, Gunthard HF: Low human immunodeficiency virus envelope diversity correlates with low in vitro replication capacity and predicts spontaneous control of plasma viremia after treatment interruptions. J Virol. 2005, 79 (14): 9026-9037. 10.1128/JVI.79.14.9026-9037.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Troyer RM, Collins KR, Abraha A, Fraundorf E, Moore DM, Krizan RW, Toossi Z, Colebunders RL, Jensen MA, Mullins JI, Vanham G, Arts EJ: Changes in human immunodeficiency virus type 1 fitness and genetic diversity during disease progression. J Virol. 2005, 79 (14): 9006-9018. 10.1128/JVI.79.14.9006-9018.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Abebe A, Demissie D, Goudsmit J, Brouwer M, Kuiken CL, Pollakis G, Schuitemaker H, Fontanet AL, Rinke de Wit TF: HIV-1 subtype C syncytium- and non-syncytium-inducing phenotypes and coreceptor usage among Ethiopian patients with AIDS. AIDS. 1999, 13 (11): 1305-1311. 10.1097/00002030-199907300-00006.View ArticlePubMedGoogle Scholar
- Cilliers T, Nhlapo J, Coetzer M, Orlovic D, Ketas T, Olson WC, Moore JP, Trkola A, Morris L: The CCR5 and CXCR4 Coreceptors Are Both Used by Human Immunodeficiency Virus Type 1 Primary Isolates from Subtype C. J Virol. 2003, 77 (7): 4449-4456. 10.1128/JVI.77.7.4449-4456.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Johnston ER, Zijenah LS, Mutetwa S, Kantor R, Kittinunvorakoon C, Katzenstein DA: High Frequency of Syncytium-Inducing and CXCR4-Tropic Viruses among Human Immunodeficiency Virus Type 1 Subtype C-Infected Patients Receiving Antiretroviral Treatment. J Virol. 2003, 77 (13): 7682-7688. 10.1128/JVI.77.13.7682-7688.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Pollakis G, Abebe A, Kliphuis A, Chalaby MIM, Bakker M, Mengistu Y, Brouwer M, Goudsmit J, Schuitemaker H, Paxton WA: Phenotypic and Genotypic Comparisons of CCR5- and CXCR4-Tropic Human Immunodeficiency Virus Type 1 Biological Clones Isolated from Subtype C-Infected Individuals. J Virol. 2004, 78 (6): 2841-2852. 10.1128/JVI.78.6.2841-2852.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Korber B, Gaschen B, Yusim K, Thakallapally R, Kesmir C, Detours V: Evolutionary and immunological implications of contemporary HIV-1 variation. Br Med Bull. 2001, 58: 19-42. 10.1093/bmb/58.1.19.View ArticlePubMedGoogle Scholar
- Richman DD, Wrin T, Little SJ, Petropoulos CJ: Rapid evolution of the neutralizing antibody response to HIV type 1 infection. Proc Natl Acad Sci USA. 2003, 100 (7): 4144-4149. 10.1073/pnas.0630530100.PubMed CentralView ArticlePubMedGoogle Scholar
- Arendrup M, Nielsen C, Hansen JE, Pedersen C, Mathiesen L, Nielsen JO: Autologous HIV-1 neutralizing antibodies: emergence of neutralization-resistant escape virus and subsequent development of escape virus neutralizing antibodies. J Acquir Immune Defic Syndr. 1992, 5 (3): 303-307.View ArticlePubMedGoogle Scholar
- Wei X, Decker JM, Wang S, Hui H, Kappes JC, Wu X, Salazar-Gonzalez JF, Salazar MG, Kilby JM, Saag MS, Komarova NL, Nowak MA, Hahn BH, Kwong PD, Shaw GM: Antibody neutralization and escape by HIV-1. Nature. 2003, 422 (6929): 307-312. 10.1038/nature01470.View ArticlePubMedGoogle Scholar
- Khouri YF, McIntosh K, Cavacini L, Posner M, Pagano M, Tuomala R, Marasco WA: Vertical Transmission of HIV-1. Correlation with maternal viral load and plasma levels of CD4 binding site anti-gp120 antibodies. J Clin Invest. 1995, 95 (2): 732-737.PubMed CentralView ArticlePubMedGoogle Scholar
- Ganeshan S, Dickover RE, Korber BT, Bryson YJ, Wolinsky SM: Human immunodeficiency virus type 1 genetic evolution in children with different rates of development of disease. J Virol. 1997, 71 (1): 663-677.PubMed CentralPubMedGoogle Scholar
- 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 (14): 1709-1717. 10.1097/00002030-199714000-00007.View ArticlePubMedGoogle Scholar
- Bradney AP, Scheer S, Crawford JM, Buchbinder SP, Montefiori DC: Neutralization escape in human immunodeficiency virus type 1-infected long-term nonprogressors. J Infect Dis. 1999, 179 (5): 1264-1267. 10.1086/314711.View ArticlePubMedGoogle Scholar
- Lathey JL, Pratt RD, Spector SA: Appearance of autologous neutralizing antibody correlates with reduction in virus load and phenotype switch during primary infection with human immunodeficiency virus type 1. J Infect Dis. 1997, 175 (1): 231-232.View ArticlePubMedGoogle Scholar
- McGrath KM, Hoffman NG, Resch W, Nelson JA, Swanstrom R: Using HIV-1 sequence variability to explore virus biology. Virus Res. 2001, 76 (2): 137-160. 10.1016/S0168-1702(01)00271-4.View ArticlePubMedGoogle Scholar
- Mantina H, Kankasa C, Klaskala W, Brayfield B, Campbell J, Du Q, Bhat G, Kasolo F, Mitchell C, Wood C: Vertical transmission of Kaposi's sarcoma-associated herpesvirus. Int J Cancer. 2001, 94 (5): 749-752. 10.1002/ijc.1529.View ArticlePubMedGoogle Scholar
- Cheng-Mayer C, Shioda T, Levy JA: Host range, replicative, and cytopathic properties of human immunodeficiency virus type 1 are determined by very few amino acid changes in tat and gp120. J Virol. 1991, 65 (12): 6931-6941.PubMed CentralPubMedGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25 (24): 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Jobb G, von Haeseler A, Strimmer K: TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics. BMC Evol Biol. 2004, 4 (1): 18-. 10.1186/1471-2148-4-18.PubMed CentralView ArticlePubMedGoogle Scholar
- Swofford DL: PAUP* Phylogenetic Analysis using parsimony and other methods. Version 4.0 beta edn. 1998, Sunderland, MA: Sinauer AssociatesGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14 (9): 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA: Molecular Evolutionary Genetics Analysis. 1.01 edn. 1993, University Park, PA 16802: The Pennsylvania State UniversityGoogle Scholar
- Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19 (18): 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMedGoogle Scholar
- Kosakovsky Pond SL, Frost SDW: Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics. 2005, bti320-Google Scholar
- Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13 (5): 555-556.PubMedGoogle Scholar
- Nielsen R, Yang Z: Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics. 1998, 148 (3): 929-936.PubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.