Analysis of HIV Type 1 BF Recombinant Sequences from South America Dates the Origin of CRF12

Abstract

HIV-1 epidemics in South America are believed to have originated in part from the subtype B epidemic initiated in the Caribbean/North America region. However, circulation of BF recombinants in similar proportions was extensively reported. Information currently shows that many BF recombinants share a recombination structure similar to that found in the CRF12_BF. In the present study, analyzing a set of 405 HIV sequences, we identified the most likely origin of the BF epidemic in an early event of recombination. We found that the subtype B epidemics in South America analyzed in the present study were initiated by a founder event that occurred in the early 1970s, a few years after the introduction of these strains in the Americas. Regarding the F/BF recombinant epidemics, by analyzing a subtype F genomic segment within the viral gene gag present in the majority of the BF recombinants, we found evidence of a geographic divergence very soon after the introduction of subtype F strains in South America. Moreover, through analysis of a subtype B segment present in all the CRF12_BF-like recombination structure, we estimated the circulation of the subtype B strain that gave rise to that recombinant structure around the same time period estimated for the introduction of subtype F strains. The HIV epidemics in South America were initiated in part through a founder event driven by subtype B strains coming from the previously established epidemic in the north of the continent. A second introduction driven by subtype F strains is likely to have encountered the incipient subtype B epidemic that soon after their arrival recombined with them, originating the BF epidemic in the region. These results may explain why in South America the majority of F sequences are found as BF recombinants.

Introduction

T he human immunodeficiency virus type 1 (HIV-1) is the etiological agent of the acquired immunodeficiency syndrome (AIDS) and the cause of death for 2 million individuals per year.¹ Although AIDS was first identified in 1981^2,3 and HIV was recognized as the etiological agent in 1983,⁴ there is an estimated cryptic time of several decades during which the virus was circulating among the human population. The emergence of the M group that contains the HIV strains responsible for the majority of infections was estimated to have occurred around 1930⁵ and the arrival to America of subtype B strains was estimated around the late 1960s,⁶ several years before the first cases of AIDS were detected.

After being introduced into the American continent, the subtype B strains subsequently dispersed to the rest of America and to the rest of the world, initiating the different subtype B regional epidemics. However, several studies on molecular characterization have shown that HIV-1 epidemics in South America are constituted not only by subtype B strains but also by subtype F strains, as well as recombinants between both subtype B and F.^{7

–18} In particular, the prevalence of BF recombinants is significantly higher than that observed for pure subtype F strains, reaching a prevalence similar to subtype B strains in some regions.¹⁴ Considering the genetic distance between subtypes B and F, the presence of the BF recombinants in South America can only be explained by a second introduction of HIV-1 into the American Continent, independent of what occurred in the Caribbean/U.S. region with subtype B strains and, in particular, driven by subtype F strains. Current evidence about the origin of the BF epidemic was obtained from characterization of the recombinant structures. Those studies found that the majority of the BF recombinant sequences share common breakpoints that are identical to those found in circulating recombinant form 12 (CRF12_BF),^{12

–16,19} which was the first recombinant to be identified in the region, suggesting a common BF recombinant ancestor. Later, new HIV circulating recombinant forms (CRFs) between subtype B and F, such as CRF28_BF and CRF29_BF, were identified in the region.²⁰ In addition, analysis of recombination patterns of BF recombinants circulating in South America has shown that ongoing current recombination drives a dynamic epidemic leading to the rise of new unique recombinant forms.²¹

In the present study we analyzed the origin of the BF HIV epidemic and, in particular, that of the CRF12_BF that previous studies has proposed to be one of the most ancestral recombinant structures in the region. To assess this issue, we analyzed through phylogenetic methods a set of sequences from BF recombinants and subtype B strains that together span a significant proportion of the time period since the regional epidemics were established.

Materials and Methods

Sequences analyzed

Sequences analyzed in the present study were either obtained from the Los Alamos National Laboratory HIV sequence database (http://www.hiv.lanl.gov/) or were previously characterized.²² The set of sequence includes the viral gene gag, the viral gene env, and full-length viral genomes as schematized in Fig. 1 and described as follows. Dataset 1: In this dataset a total of 182 sequences of the gag gene were included. Seventeen were subtype reference sequences obtained from the Los Alamos National Laboratory Data Base and 126 sequences were previously obtained in our laboratory from samples collected from newly diagnosed individuals between 1987 and 2006.²² The remaining 39 were other subtype F or BF recombinant sequences obtained from the Los Alamos National Laboratory Data Base (http://hiv-web.lanl.gov). Dataset 2 : In this dataset a total of 105 sequences of the env gene were included. This dataset contains the one previously analyzed by Gilbert et al. ⁶ and 37 additional env sequences from previously characterized complete subtype B genomes from South America, mainly from Argentina and Brazil because of sequence availability. The overall dataset covers the years 1981 to 2004. Dataset 3 : In this dataset subtype B sequences of the genomic region comprised between positions 3087 and 3597 in HXB2 (within pol gene) were included. Thirty-four of them were extracted from full-length recombinant BF genomes available online, together with other 56 sequences extracted from full-length subtype B genomes found to be related to American strains.

FIG. 1.

Viral genome regions under analysis. The recombination patterns of 10 full-length genomes previously characterized are shown as an example of diversity. The first three structures are recombinant forms circulating in Brazil. The fourth is the circulating recombinant form number 12. The other six structures are other unique recombinant forms identified in Argentina. One complete subtype B genome is represented by the sequence ARMA132. The three viral genome regions under analysis are highlighted in the schematic representation of the HIV genome (solid gray square) and in the sequences (open gray square). The gag region was analyzed in all the BF recombinants, the env region was analyzed in all the subtype B sequences, and the pol region was analyzed only in those BF recombinants that were subtype B in the region under analysis (in the example the CRF40_BF is excluded from this analysis).

For each dataset, alignments were performed by using Clustal W (BioEdit 7.0.4.1 sequence alignment editor²³). Codon optimization was then performed with Gene Cutter [B. Gaschen (Los Alamos National Laboratory); http://www.hiv.lanl.gov/content/hivdb/] and after gap-stripping a neighbor-joining (NJ) tree was constructed for each dataset. Duplicated or resampled sequences were identified as those clustering together with a bootstrap support higher than 99%, evidencing a common origin in the sequence information and reduced genetic distance according to branch lengths. After eliminating duplicated or resampled sequences, quality analysis was assessed by calculating the number of mismatches with a data-derived consensus sequence obtained from the corresponding alignment. Confidence intervals were estimated for the mean mismatched proportion and sequences classified as hypermutated were eliminated. Hypermutated sequences were rich in ambiguities.

Viral subtype characterization

Sequences were aligned with subtype references sequences available (http://www.hiv.lanl.gov/content/sequence/NEWALIGN/align.html) using Clustal W in the BioEdit 7.0.4.1 sequence alignment editor,²³ and codon optimization was performed using Gene Cutter [B. Gaschen (Los Alamos National Laboratory); http://www.hiv.lanl.gov/content/hivdb/] and hand aligning. After gap stripping, NJ trees were constructed under the Kimura two-parameter model with MEGA4.²⁴ All sequences were individually analyzed for similarity with consensus references by SimPlot v2.5 [S.C. Ray (Johns Hopkins University); http://sray.med.som.jhmi.edu/SCRoftware/simplot/]. Sequences bearing similarities to two or more different subtypes were further analyzed by bootscanning also using SimPlot v2.5 with a windows size of 100 bases and a step size of 20 and by visual inspection of alignments in order to identify recombination breakpoints.

The parsimony analysis of datasets

Phylogenetic trees based on parsimony method were obtained using TNT.²⁵ This software implements new technologies for dealing with large datasets^26,27 aimed to provide a thorough exploration of the tree space. This ensures that all the possible phylogenetic hypotheses that could be supported by the data are found.²⁷ Thus, any phylogenetic uncertainty present in the data is summarized in the corresponding strict consensus tree. Tree searches were performed as described elsewhere.²² Briefly, we obtained 100 Random Addition Sequence trees, or Wagner (W) trees,²⁸ and the shorter ones were submitted to branch swapping using three rounds of a combination of Tree Fusing, Ratchet, Tree Drift, and Sectorial Search (W + T-R-Df-SS).²⁶ The best trees found in each W + T-R-Df-SS round were saved for ‘‘feeding’’ the subsequent W + T-R-Df-SS rounds (builder). Ten independent runs of the builder were performed and the best trees obtained from each builder were pooled and subjected to further cycles of T-R-Df-SS to verify that no further improvements of length were possible, at least under this experimental setting. Ambiguously supported branches were automatically collapsed during searches.

Estimation of the time to the most recent common ancestor

The BEAUti/BEAST v1.4.8 package [A. J. Drummond (University of Auckland, Auckland, New Zealand) and A. Rambaut (University of Edinburgh); beast.bio.ed.ac.uk] was used to perform the Bayesian MCMC estimation of the TMRCA for different taxon subsets. The year of sample collection was assigned to each sequence. The nucleotide substitution model was obtained with ModelTest²⁹ and PAUP4.0b10 [D.L. Swofford (Smithsonian Institution, Washington, DC); http://paup.csit.fsu.edu/]. The analysis was performed with a Bayesian Skyline coalescent tree under the Uncorrelated Lognormal Relaxed Molecular Clock model. An initial run of 1 million generations was performed for each analysis in order to optimize the different operators. Then new analyses with 10 million generations were performed in order to improve the optimizations. The final run was performed with 8 independent MCMC analyses, each one of 10 million generations with a burn-in period of 2 million generations for dataset 1 and 3 million generations for dataset 2. The MCMC samples were analyzed with Tracer v1.4.1 [A. Rambaut (University of Edinburgh) and A. J. Drummond (University of Auckland, Auckland, New Zealand); http://beast.bio.ed.ac.uk].

Statistical analysis

Chi-square test and Fisher's exact test were used to test the hypothesis of association between sampling country and phylogenetic clade.

Results

The phylogeography of subtype F

We first analyzed a total of 165 sequences of the gag gene from strains circulating in South America, previously characterized as BF recombinants,²² in order to identify the different recombination patterns. We selected this dataset for analysis because it covers the broader range of years in the BF epidemic. After alignment and optimization the genomic viral region under analysis was 1147 bases in length. Results obtained through recombination analysis performed by bootscanning showed that the BF recombinant HIV sequences had one of either two different recombinant structures. Both structures had only one recombination breakpoint, although located in a different position: one of them had the recombination breakpoint around nucleotide 200 (989 in HXB2) in our alignment and the other one had the recombination breakpoint around nucleotide 500 (1301 in HXB2). The recombination pattern was found to be significantly associated with the country of sampling (p < 0.001, Fisher's exact test). All the sequences with a BF200-like recombinant structure came from Argentina, while 14 out of the 26 (53.8%) BF500-like structures came from Brazil. Comparative analysis of recombination breakpoints showed that the BF200 pattern is present in the recombination structure of circulating recombinant form 12 (CRF12_BF).

Next, in order to track the origin of the subtype F strain that gave rise to these BF recombinant forms, we performed a phylogenetic analysis over the larger region in which all the sequences belonged to the viral subtype F in order to avoid confounding effects due to recombination. The region analyzed was located between nucleotide positions 1401 and 1951 in HXB2 and constitute Dataset 1 (Fig. 1). After trimming the alignment, all the sequences were reanalyzed by bootscanning in order to confirm the absence of recombination breakpoints within the region, and the phylogenetic history was inferred by maximum parsimony and Bayesian methods. By using the first method, a total of 38 optimal trees was found and the strict consensus was constructed. As shown in Fig. 2A, this analysis indicated that there is a monophyletic clade that groups together all the sequences obtained in the American continent, whereas all the sequences collected outside the continent are paraphyletic. Within the American clade, two reciprocally monophyletic clades cluster with significantly different proportions (p < 0.001): the sequences from Brazil and the sequences from Argentina. Within the “Brazilian clade” only four out of the 25 (16%) sequences had a BF200-like structure, which, in addition, clustered in a monophyletic clade. In contrast, within the “Argentinean clade” the BF200-like structure was in the 48.5% (n = 67) of the sequences, the BF500-like structure in 12.3% (n = 17, five of them in a monophyletic clade), and the remaining 39.2% (n = 54) were pure subtype F sequences in the gag gene. Although the Bayesian method could not resolve phylogenetic relationships inside subtype F (Supplementary Fig. S1; Supplementary Data are available online at www.liebertonline.com/aid), the parsimony analysis was supported by a geographic correlate for the phylogenetic reconstruction, finding sequences sampled in Argentina more frequently located within the “Argentinean clade” in the tree and sequences sampled from Brazil more frequently located within the “Brazilian clade” in the tree (Fig. 2A).

FIG. 2.

Phylogenetic reconstructions of regional HIV-1 epidemics. The analysis was performed by a maximum parsimony method. Branches are colored in green for samples of Brazilian origin and in blue for samples of Argentinean origin. (A) Strict consensus constructed from 38 equally optimal trees for the F segment of all the F and BF recombinants forms available. The recombination structures of the sequences analyzed are detailed (green circles: BF500; blue circles: BF200). Branches are colored in black for samples with an origin outside of the American Continent. (B) Strict consensus constructed from the four equally optimal trees found for the subtype B env region. Samples from Haiti are labeled as “HT” and samples from Trinidad and Tobago are labeled as “TT”. Also the “Brazilian clade” and the “Argentinean clade” are detailed. Branches are colored in black for samples from the United States. (C) Strict consensus constructed from the 49 equally optimal trees found for the B segment shared by all the CRF12_BF genomes and the majority of Brazilian BF sequences. The CRF12_BF-like sequences are detailed with open red circles and the CRF12_BF clade is detailed. Branches are colored in black for samples from other regions of the American Continent.

The origin of the subtype B epidemic in South America

To assess the origin of the subtype B epidemic in the region we performed a phylogenetic analysis over the env gene (Dataset 2, positions 6221–8795 in HXB2, Fig. 1). We based our analysis on the dataset previously described by Gilbert et al.⁶ because it is the most complete evidence available for reconstruction of the subtype B epidemic in the American continent and included additional sequences previously obtained by our group and others obtained from Los Alamos Sequence Database. By using the parsimony method for phylogenetic reconstruction we were able to obtain four optimal trees. Different search strategies were not able to obtain shorter trees. As shown in Fig. 2B, we found that the majority of Argentinean and Brazilian sequences clustered in two reciprocally monophyletic clades. As for the subtype F sequences in dataset 1, character optimization suggests a Brazilian origin for one of them and an Argentinean origin for the other. The “Brazilian clade” has 14 sequences from Brazil and 10 sequences from Argentina. The “Argentinean clade” has 11 sequences from Argentina, 2 sequences from Brazil, and 3 sequences from the United States. The Bayesian-based phylogenetic reconstruction found the same two clades with a posterior probability of 1.00 and 0.78, respectively (Supplementary Fig. S2).

Phylogenetic reconstruction of subtype B regions in CRF12_BF resembles the monophyletic origin of subtype F

Considering previous reports that suggest that CRF12_BF is an ancestral recombination pattern we analyzed the larger subtype B segment present in the CRF12_BF and also in the majority of BF recombinants over a total of 102 sequences in order to assess the origin of the subtype B strain that originated the BF recombinant forms in the region. After gap striping, the region analyzed was of 500 bases in length, located within the pol gene between positions 3087 and 3597 in HXB2 (Fig. 1, dataset 3). After performing phylogenetic reconstruction through the parsimony method we obtained 49 optimal trees. As shown in Fig. 2C, the strict consensus has all the 13 CRF12_BF sequences clustering in a clade with other BF recombinants and subtype B sequences suggesting a monophyletic origin of subtype B genomic segments of CRF12_BF strains (“CRF12_BF clade”). The phylogenetic relationship for the remaining sequences could not be resolved and are collapsed in the node of subtype B. As shown in Fig. 3, the Bayesian-based phylogenetic reconstruction found the same clade in the majority-rule consensus tree with a posterior probability of 0.56 (see also Supplementary Fig. S3 and Supplementary Table S1).

FIG. 3.

Comparison of the phylogenetic reconstructions for the CRF12_BF clade obtained by the maximum parsimony and Bayesian method. CRF12_BF-like recombination structures are detailed with black branches and other BF recombinant structures are detailed in white. Subtype B sequences are detailed in gray. The viral subtype of the complete sequence is detailed in the first two letters of the name. The country of origin is detailed next with a two-letters code (AR: Argentina; BR: Brazil; US: United States).

Timing the origin of HIV epidemics in the region

Finally, considering the phylogenies reconstructed above, we aimed to estimate the time to the most recent common ancestor (TMRCA) for each of the subtype B and subtype F epidemics. Because the dating analysis is critically dependent on the monophyly of the clades for which the TMRCA is to be estimated, we performed the analysis only in those cases where that condition was satisfied according to phylogenetic support. The estimations of the TMRCAs and substitution rates were performed with the Bayesian Markov chain Monte Carlo (MCMC) phylogenetic analysis that allows the substitution rate to vary among branches in the tree.

In the case of the subtype F sequences we estimated the TMRCAs for the American subtype F clade, that is, for the ancestral node containing both “Brazilian clade” and “Argentinean clade” in Fig. 2A. Our analysis performed over dataset 1 found the maximum probability for the ancestor to that clade in the late 1960s/beginning of the 1970s (TMRCA: 1969; 95%HPD: 1959–1978, Fig. 4A). The majority-rule consensus tree obtained from the posterior sample of this analysis confirmed the monophyletic nature of this clade (Supplementary Fig. S1). In addition, we estimated the maximum probability for the time of divergence between viral subtypes F1 and F2 to be around in the late 1950s/early 1960s (TMRCA: 1956, 95%HPD: 1935–1972).

FIG. 4.

Estimation of the time to the most recent common ancestor (TMRCA). The density plots for the different TMRCAs estimations are shown. (A) TMRCA for the subtype F epidemic. The first distribution represents the divergence between subtypes F1 and F2. The second distribution represents the TMRCA for the American subtype F strains. (B) TMRCA for the subtype B epidemics. The main TMRCA estimated for this study is represented by the third distribution, which is an overlap of two almost identical distributions that represent the TMRCAs for the Argentinean subtype B clade and for the Brazilian subtype B clade. This information is shown in the context of two other TMRCAs also estimated with the same dataset and previously found by Gilbert et al. that represent the divergence between subtypes B and D (first distribution) and the introduction of subtype B in the Americas (second distribution). (C) TMRCA for the CRF12_BF. The distribution shown in this figure represents the TMRCA for the subtype B strain that gave rise to the recombinants with a CRF12_BF-like structure. (*) Distributions marked with an asterisk are those obtained from the dataset under analysis but previously reported elsewhere.⁶

In the case of the subtype B sequences we estimated over dataset 2 the TMRCAs for both the Argentinean and Brazilian clades shown in Fig. 2B. Our analysis found identical results for both clades: located at 1972 (95%HPD: 1970–1974) for the Argentinean clade and 1972 (95%HPD: 1969–1974) for the Brazilian clade (Fig. 4B). The majority-rule consensus tree obtained from the posterior sample of this analysis confirms the monophyletic nature of these two clades with a posterior probability of 0.78 and 1.00, respectively (Supplementary Fig. S2). Regarding the substitution rates, we found a rate of 0.00419 (0.00374–0.00464) substitutions/site/year for the env gene of subtype B and 0.00167 (0.00115–0.00210) substitutions/site/year for the gag gene of subtype F.

Finally, we aimed to estimate the TMRCA for the CRF12_BF by analyzing dataset 3. In this case, the maximum parsimony tree differed, within the “CRF12_BF clade,” from the Bayesian tree in one sequence: as shown in Fig. 3 the recombinant BF sequence BREPM12609 (Accession number: DQ085873) was located outside the “CRF12_BF clade” in the maximum parsimony tree. Also, we found within this clade, a more resolved phylogeny when the parsimony method was applied compared to that obtained with the Bayesian method in which the majority of nodes were collapsed. In spite of this, the Bayesian method still resolves the “CRF12_BF clade” and therefore allows us to estimate the TMRCA. By carrying out this analysis we estimated that the maximum probability for the most recent common ancestor of the CRF12_BF was the early 1970s (TMRCA: 1969, 95%HPD: 1946–1981, Fig. 4C).

Discussion

In the present study we looked for an explanation of the particular HIV-1 epidemic that is present in South America where, for many years, several studies have consistently reported the cocirculation of subtype B strains and BF recombinants but a very low prevalence of pure subtype F strains.^{10,12

–16,30,31} Because we were analyzing the origin of the epidemic in South America and considering that the first detected cases of AIDS in the region occurred in Brazil and Argentina, we looked mainly for sequences derived from those two countries. In addition, the number of HIV sequences from other countries is very limited. The similarity in the recombination structures of the BF sequences suggested a common recombinant ancestor, but that hypothesis was never exhaustively evaluated before. The advantage of having, on the one hand, sequences of both subtype B and recombinants BF recovered from stored plasma of the past 20 years in our region and, on the other hand, an accurate phylogenetic reconstruction of subtype B in the Americas obtained by Gilbert et al. ⁶ allowed us to dissect the most ancestral phylogenetic relationships in the initiation of the recombinant BF epidemic.

Because TMRCA calculations critically depend on the monophyletic nature of regional epidemics and considering that monophyletic clades could be a consequence of the lack of sequences that, if included, could split the identified clusters, we included in this study the most representative set of sequences available. In this sense, the subtype F and recombinants BF sequences analyzed were obtained from plasma samples stored from 1987 through 2006 in five different hospitals in Argentina that function as voluntary and counseling testing (VCT) sites.

Because VCTs are considered one of the preferred sources of samples for epidemiological studies over the general population,³² we expect those sequences to be a representative sample of the local epidemic. On the other hand, sequences obtained from two cross-sectional prevalence studies performed over two at-risk cohorts of MSM³³ and female sex workers³¹ designed to obtain a representative sample of each group were also included. Sequences from other regions were limited to the availability in databases. However, in all the monophyletic clades analyzed, we found some level of heterogeneity in terms of the sequence origin, suggesting that we are likely analyzing unbiased samples of viral sequences. The same applies to subtype B sequences, with the exception that for this case we also added those sequences previously described by Gilbert et al.⁶ because previous epidemiological data suggested that the local epidemics in the American continent were related. In fact, it is known that one arm of the global dispersion of HIV began by the introduction in the Caribbean/U.S. region of subtype B strains that subsequently dispersed to the rest of America and to the rest of the world.^5,6,34 In the particular case of CRF12_BF, the 13 complete genomes sequences were obtained from three different sources: five of them were obtained from the cross-sectional study of prevalence in female sex workers mentioned above,³¹ three of them from a study of pregnant women seeking HIV testing,^15,35 and five of them from a study of drug resistance over primary and chronically infected patients who attended to a VCT site.^12,36 Therefore, although the number of complete CRF12_BF genomes available is limited, in the phylogenetic reconstruction it is likely that they represent distantly related sequences within a hypothetical group of all the CRF12_BF circulating.

The fact that all the South American subtype B sequences grouped in only two monophyletic clades with a geographic correlate, which is also the case for subtype F sequences, suggests that these epidemics are predominantly fed by local events of transmission even when they occur in populations of geographically close regions as is the case of Argentina and Brazil. The presence of sequences from Brazil in the Argentinean clade and sequences of Argentina in the Brazilian clade shows that circulation among regions actually occurs as recently suggested,²¹ although the monophyletic nature of both clades suggests an independent origin by one or a few related strains.

Previous studies have reported the initiation of HIV epidemics through a founder effect in San Francisco³⁷ and Trinidad and Tobago,³⁸ and other regions of the world,^39

–43 and in our study this may be explained by a divergence of the subtype F strains in both populations at the very early moment of the colonization, after their introduction in South America probably from Africa, as subtype F1 is not found in other regions of the world. For subtype B, considering our results, a likely explanation would be that one or a few related strains imported from the subtype B epidemic established in North America successfully spread and established these local epidemics. Taking into account that previous reports have shown that HIV transmission is primarily attributable to HIV-vulnerable groups,⁴⁴ i.e., concentrated epidemic,⁴⁵ we consider a possible explanation for the fact that a few strains successfully spreading over others that finally were limited. Also, the dates estimated for the TMRCA of these local epidemics in South America suggest that there had been a significant number of autochthonous infections by the time the first cases started to be identified.

Our estimations for the Brazilian clade subtype B are consistent with those previously reported by Bello et al. ⁴⁶ and almost identical to our estimations for the Argentinean clade B, suggesting that the introduction of subtype B in Argentina and in Brazil occurred almost at the same time and very soon after the introduction of subtype B in America. Based on these results, it is likely that HIV strains that circulate in different regions may have differentiated genetic characteristics even when they belong to the same viral subtypes and the differences may be attributable to founder events in the origin of the epidemics in different regions. In this sense, it would be interesting to perform phylogeography studies of HIV epidemics in other regions of the world to evaluate whether the majority of them were initiated by founder events and to better understand the genetic characteristics of strains currently circulating.

In addition, the fact that the subtype B regions of the different CRF12_BF-like recombinant genomes have a common ancestor closer to themselves than to any other subtype B sequences suggests that an event of recombination occurred in the very early years of the epidemic generating a recombinant form, subsequently identified as CRF12_BF. In fact, when we analyzed the subtype B region of the BF recombinants, we found them to have a TMRCA distribution similar to that estimated for the American subtype F strains (Fig. 4). This means that the subtype B strain that originated the CRF12_BF by recombination with a subtype F strain was circulating in the region at the same time as the subtype F strains arrived to America. Therefore, it is likely that the recombination events occurred in the early moments of the initiation of the F/BF epidemic. This is in accordance with the observation that subtype F is almost absent in the region while the majority of non-B strains are recombinant BF and also with the fact that the majority of BF recombinants share similar patterns of recombination highly related to those found in CRF12_BF.^{12

–16,19}

Regarding the estimated evolutionary rates, our estimation of around 0.00162 substitution per year per site for the gag dataset is similar to that obtained for the gag gene by Korber et al.⁵ (0.0009–0.0027) and our estimation of around 0.00419 for the env gene is also similar to that estimated for env by Leitner et al. ⁴⁷ (0.0046–0.0088) and Robbins et al. ³⁴ (0.00473), but higher than that obtained by Korber et al. ⁵ for env (0.0018–0.0028). The last difference is likely to be due to the fact that posterior studies demonstrated that a relaxed molecular clock better explains the phylogenetic reconstruction of HIV sequences. In this sense, we also found a significant variation in evolutionary rates among branches in the tree as shown by the standard deviation of the uncorrelated lognormal relaxed clock found in our analysis (ucld.std 95%HD: 0.195 ± 0.037 for env gene B and 0.456 ± 0.076 for gag gene F). Recently, Abecasis et al. ⁴⁸ estimated the evolutionary rates for the different HIV-1 subtypes, obtaining a rate for subtype B very similar to ours when analyzing the env gene [a mean rate of 0.0045 substitutions/site/year compared to the 0.00419 (0.00374–0.00464) estimated by us]. In addition, while Abecasis et al. ⁴⁸ found subtype G to have the higher evolutionary rate compared to subtypes B, C, A1, D, and recombinants AG and AE, contextualizing our estimations in these set of results, subtype F in the gag gene seems to have an evolutionary rate similar to that observed for the pol gene in the other HIV-1 subtypes [0.00167 (0.00115–0.00210) compared to a rate between 0.001 and 0.002 for the majority of the other subtypes].

Finally, it is important to remember that the HIV-1 epidemic in South America is also composed by a minority, although still an important proportion, of subtype C infections. In the present study we limited our analysis to both subtype B and subtype F/recombinant BF epidemics because those are the viral subtypes found in the oldest samples obtained from patients in the region. In fact, a later initiation of the subtype C epidemic in South America compared to the subtype B epidemic was previously reported⁴⁹ with an estimated initiation at the beginning of the 1990s. The higher prevalence of subtype B and recombinants BF compared to subtype C infections may be a consequence of an earlier initiation of the first ones, although additional factors influencing the evolutionary rate of circulating strains, such as the significantly higher rate estimated previously for subtype C strains in Brazil,⁴⁹ make the current scenario not necessarily stable and may lead to a shift in the future of the epidemic.

Conclusions

In conclusion, by applying parsimony and Bayesian methods to the reconstruction of the evolutionary history of HIV in America we were able to dissect the origin of regional HIV epidemics after their introduction in the Americas. According to our results subtype B strains of HIV coming from the north of the Americas arrived in South America in the early 1970s and established the South American subtype-B epidemic through a founder event.

The relative short period of time between the estimated origin of the U.S. subtype B epidemic and the South American U.S.-related subtype B epidemic suggests rapid spreading of the infection to the whole American continent after successful colonization through the Caribbean/U.S. region. At the same time a second wave of HIV strains arrived from Africa and established the subtype F epidemic of South America. Soon after the introduction of subtype F strains in this region, it is likely that a recombination with a subtype B strain had occurred giving rise to a recombinant form with a CRF12_BF-like recombination pattern. This founder recombinant event established a BF epidemic, giving rise to other unique recombinant forms by further recombination with subtype B strains. Finally, our results suggest that by the time the first cases of AIDS were identified, not only the Caribbean/North American HIV epidemic but also several regional HIV epidemics in South America were already established.

Footnotes

Acknowledgments

We thank Dr. Michael Worobey, Department of Ecology and Evolutionary Biology University of Arizona, for kindly providing the dataset for env analysis previously published by his group. DAD acknowledges the support from the Fogarty International Center/NIH grant through the AIDS International Training and Research Program at Mount Sinai School of Medicine-Argentina Program: Grant D43 TW001037.

Author Disclosure Statement

No institutional of commercial affiliations that might pose a conflict of interest regarding the publication of this manuscript exist.

References

UNAIDS: AIDS epidemic update. Geneva, UNAIDS, Geneva, 2009.

Epidemiologic aspects of the current outbreak of Kaposi's sarcoma and opportunistic infections. N Engl J Med, 1982; 306:248–252.

CDC: Pneumocystis Pneumonia––Los Angeles. Morb Mortal Weekly Report, 1981; 30:1–3.

Barre-Sinoussi

, Chermann

, Rey

, Nugeyre

, Chamaret

et al.

Isolation of a T-lymphotropic retrovirus from a patient at risk for acquired immune deficiency syndrome (AIDS)

Science, 1983; 220:868–871.

Korber

, Muldoon

, Theiler

, Gao

, Gupta

et al. Timing the ancestor of the HIV-1 pandemic strains. Science, 2000; 288:1789–1796.

Gilbert

, Rambaut

, Wlasiuk

, Spira

, Pitchenik

et al. The emergence of HIV/AIDS in the Americas and beyond. Proc Natl Acad Sci USA, 2007; 104:18566–18570.

Gomez-Carrillo

, Pampuro

, Duran

, Losso

, Harris

et al. Analysis of HIV type 1 diversity in pregnant women from four Latin American and Caribbean countries. AIDS Res Hum Retroviruses, 2006; 22:1186–1191.

Monteiro

, Alcantara

, de Oliveira

, Oliveira

, Melo

et al. Genetic variability of human immunodeficiency virus-1 in Bahia state, Northeast, Brazil: High diversity of HIV genotypes. J Med Virol, 2009; 81:391–399.

de Souza

, de Oliveira

, Rodrigues

, Silva

, Levi

. Molecular characterization of HIV type 1 BF pol recombinants from Sao Paulo, Brazil. AIDS Res Hum Retroviruses, 2008; 24:1521–1525.

10.

Brennan

, Brites

, Bodelle

, Golden

, Hackett

et al. HIV-1 strains identified in Brazilian blood donors: Significant prevalence of B/F1 recombinants. AIDS Res Hum Retroviruses, 2007; 23:1434–1441.

11.

Pereira

, Stefani

, Araujo Filho

, Souza

, Stefani

et al. Human immunodeficiency virus type 1 (HIV-1) and Mycobacterium leprae co-infection: HIV-1 subtypes and clinical, immunologic, and histopathologic profiles in a Brazilian cohort. Am J Trop Med Hyg, 2004; 71:679–684.

12.

Thomson

, Villahermosa

, Vazquez-de-Parga

, Cuevas

, Delgado

et al. Widespread circulation of a B/F intersubtype recombinant form among HIV-1-infected individuals in Buenos Aires, Argentina. AIDS, 2000; 14:897–899.

13.

Quarleri

, Rubio

, Carobene

, Turk

, Vignoles

et al. HIV type 1 BF recombinant strains exhibit different pol gene mosaic patterns: Descriptive analysis from 284 patients under treatment failure. AIDS Res Hum Retroviruses, 2004; 20:1100–1107.

14.

Dilernia

, Gomez

, Lourtau

, Marone

, Losso

et al. HIV type 1 genetic diversity surveillance among newly diagnosed individuals from 2003 to 2005 in Buenos Aires, Argentina. AIDS Res Hum Retroviruses, 2007; 23:1201–1207.

15.

Carr

, Avila

, Gomez Carrillo

, Salomon

, Hierholzer

et al. Diverse BF recombinants have spread widely since the introduction of HIV-1 into South America. AIDS, 2001; 15:F41–47.

16.

Avila

, Pando

, Carrion

, Peralta

, Salomon

et al. Two HIV-1 epidemics in Argentina: Different genetic subtypes associated with different risk groups. J Acquir Immune Defic Syndr, 2002; 29:422–426.

17.

Rios

, Delgado

, Perez-Alvarez

, Fernandez

, Galvez

et al. Antiretroviral drug resistance and phylogenetic diversity of HIV-1 in Chile. J Med Virol, 2007; 79:647–656.

18.

Carrion

, Laguna-Torres

, Soto-Castellares

, Castillo

, Salazar

et al. Molecular characterization of the human immunodeficiency virus type 1 among children in Lima, Peru. AIDS Res Hum Retroviruses, 2009; 25:833–835.

19.

Thomson

, Delgado

, Herrero

, Villahermosa

, Vazquez-de Parga

et al. Diversity of mosaic structures and common ancestry of human immunodeficiency virus type 1 BF intersubtype recombinant viruses from Argentina revealed by analysis of near full-length genome sequences. J Gen Virol, 2002; 83:107–119.

20.

De Sa Filho

, Sucupira

, Caseiro

, Sabino

, Diaz

et al. Identification of two HIV type 1 circulating recombinant forms in Brazil. AIDS Res Hum Retroviruses, 2006; 22:1–13.

21.

Zhang

, Foley

, Schultz

, Macke

, Bulla

et al. The role of recombination in the emergence of a complex and dynamic HIV epidemic. Retrovirology, 2010; 7:25.

22.

Dilernia

, Jones

, Rodriguez

, Turk

, Rubio

et al. HLA-driven convergence of HIV-1 viral subtypes B and F toward the adaptation to immune responses in human populations. PLoS ONE, 2008; 3:e3429.

23.

Hall

. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser, 1999; 41:95–98.

24.

Tamura

, Dudley

, Nei

, Kumar

. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol, 2007; 24:1596–1599.

25.

Goloboff

, Farris

, Nixon

. TNT, a free program for phylogenetic analysis. Cladistics, 2008; 24:774–786.

26.

Goloboff

. Analyzing large data sets in reasonable times: Solutions for composite optima. Cladistics, 1999; 15:415–428.

27.

Nixon

. The parsimony ratchet, a new method for rapid parsimony analysis. Cladistics, 1999; 15:407–414.

28.

Farris

. Methods for computing Wagner trees. Syst Zool, 1970; 19:83–92.

29.

Posada

, Crandall

. MODELTEST: Testing the model of DNA substitution. Bioinformatics, 1998; 14:817–818.

30.

Thomson

, Sierra

, Tanuri

, May

, Casado

et al. Analysis of near full-length genome sequences of HIV type 1 BF intersubtype recombinant viruses from Brazil reveals their independent origins and their lack of relationship to CRF12_BF. AIDS Res Hum Retroviruses, 2004; 20:1126–1133.

31.

Pando

, Berini

, Bibini

, Fernandez

, Reinaga

et al. Prevalence of HIV and other sexually transmitted infections among female commercial sex workers in Argentina. Am J Trop Med Hyg, 2006; 74:233–238.

32.

Guidelines for surveillance of HIV drug resistance. World Health Organization, 2003.

33.

Pando

, Eyzaguirre

, Segura

, Bautista

, Marone

et al. First report of an HIV-1 triple recombinant of subtypes B, C and F in Buenos Aires, Argentina. Retrovirology, 2006; 3:59.

34.

Robbins

, Lemey

, Pybus

, Jaffe

, Youngpairoj

et al. U.S. human immunodeficiency virus type 1 epidemic: Date of origin, population history, and characterization of early strains. J Virol, 2003; 77:6359–6366.

35.

de los Angeles Pando

, Biglione

, Toscano

, Rey

, Russell

et al. Human immunodeficiency virus type 1 and other viral co-infections among young heterosexual men and women in Argentina. Am J Trop Med Hyg, 2004; 71:153–159.

36.

Kijak

, Pampuro

, Avila

, Zala

, Cahn

et al. Resistance profiles to antiretroviral drugs in HIV-1 drug-naive patients in Argentina. Antivir Ther, 2001; 6:71–77.

37.

Foley

, Pan

, Buchbinder

, Delwart

. Apparent founder effect during the early years of the San Francisco HIV type 1 epidemic (1978–1979) AIDS Res Hum Retroviruses, 2000; 16:1463–1469.

38.

Cleghorn

, Jack

, Carr

, Edwards

, Mahabir

et al. A distinctive clade B HIV type 1 is heterosexually transmitted in Trinidad and Tobago. Proc Natl Acad Sci USA, 2000; 97:10532–10537.

39.

Bobkov

, Cheingsong-Popov

, Selimova

, Ladnaya

, Kazennova

et al. An HIV type 1 epidemic among injecting drug users in the former Soviet Union caused by a homogeneous subtype A strain. AIDS Res Hum Retroviruses, 1997; 13:1195–1201.

40.

Delwart

, Shpaer

, Louwagie

, McCutchan

, Grez

et al. Genetic relationships determined by a DNA heteroduplex mobility assay: Analysis of HIV-1 env genes. Science, 1993; 262:1257–1261.

41.

Grez

, Dietrich

, Balfe

, von Briesen

, Maniar

et al. Genetic analysis of human immunodeficiency virus type 1 and 2 (HIV-1 and HIV-2) mixed infections in India reveals a recent spread of HIV-1 and HIV-2 from a single ancestor for each of these viruses. J Virol, 1994; 68:2161–2168.

42.

Lukashov

, Karamov

, Eremin

, Titov

, Goudsmit

. Extreme founder effect in an HIV type 1 subtype A epidemic among drug users in Svetlogorsk, Belarus. AIDS Res Hum Retroviruses, 1998; 14:1299–1303.

43.

, Takebe

, Weniger

, Luo

, Kalish

et al. Independent introduction of two major HIV-1 genotypes into distinct high-risk populations in Thailand. Lancet, 1993; 341:1171–1174.

44.

Vignoles

, Avila

, Osimani

, de Los Angeles Pando

, Rossi

et al. HIV seroincidence estimates among at-risk populations in Buenos Aires and Montevideo: Use of the serologic testing algorithm for recent HIV seroconversion. J Acquir Immune Defic Syndr, 2006; 42:494–500.

45.

UNAIDS. 2008. Report on the global AIDS epidemic.

46.

Bello

, Eyer-Silva

, Couto-Fernandez

, Guimaraes

, Chequer-Fernandez

et al. Demographic history of HIV-1 subtypes B and F in Brazil. Infect Genet Evol, 2007; 7:263–270.

47.

Leitner

, Albert

. The molecular clock of HIV-1 unveiled through analysis of a known transmission history. Proc Natl Acad Sci USA, 1999; 96:10752–10757.

48.

Abecasis

, Vandamme

, Lemey

. Quantifying differences in the tempo of human immunodeficiency virus type 1 subtype evolution. J Virol, 2009; 83:12917–12924.

49.

Salemi

, de Oliveira

, Soares

, Pybus

, Dumans

et al. Different epidemic potentials of the HIV-1B and C subtypes. J Mol Evol, 2005; 60:598–605.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

2.23 MB

Analysis of HIV Type 1 BF Recombinant Sequences from South America Dates the Origin of CRF12_BF to a Recombination Event in the 1970s

Abstract

Introduction

Materials and Methods

Sequences analyzed

Viral subtype characterization

The parsimony analysis of datasets

Estimation of the time to the most recent common ancestor

Statistical analysis

Results

The phylogeography of subtype F

The origin of the subtype B epidemic in South America

Phylogenetic reconstruction of subtype B regions in CRF12_BF resembles the monophyletic origin of subtype F

Timing the origin of HIV epidemics in the region

Discussion

Conclusions

Footnotes

Acknowledgments

Author Disclosure Statement

References

Supplementary Material