Analysis of the evolutionary relationships of HIV-1 and SIVcpz sequences using bayesian inference: implications for the origin of HIV-1

The most plausible origin of HIV-1 group M is an SIV lineage currently represented by SIVcpz isolated from the chimpanzee subspecies Pan troglodytes troglodytes. The origin of HIV-1 group O is less clear. Putative recombination between any of the HIV-1 and SIVcpz sequences was tested using bootscanning and Bayesian-scanning plots, as well as a new method using a Bayesian multiple change-point (BMCP) model to infer parental sequences and crossing-over points. We found that in the case of highly divergent sequences, such as HIV-1/SIVcpz, Bayesian scanning and BMCP methods are more appropriate than bootscanning analysis to investigate spatial phylogenetic variation, including estimating the boundaries of the regions with discordant evolutionary relationships and the levels of support of the phylogenetic clusters under study. According to the Bayesian scanning plots and BMCP method, there was strong evidence for discordant phylogenetic clustering throughout the genome: (1) HIV-1 group O clustered with SIVcpzANT/TAN in middle pol, and partial vif/env; (2) SIVcpzGab1 clustered with SIVcpzANT/TAN in 3'pol/vif, and middle env; (3) HIV-1 group O grouped with SIVcpzCamUS and SIVcpzGab1 in p17/p24; (4) HIV-1 group M was more closely related to SIVcpzCamUS in 3'gag/pol and in middle pol, whereas in partial gp120 group M clustered with group O. Conditionally independent phylogenetic analysis inferred by maximum likelihood (ML) and Bayesian methods further confirmed these findings. The discordant phylogenetic relationships between the HIV-1/SIVcpz sequences may have been caused by ancient recombination events, but they are also due, at least in part, to altered rates of evolution between parental SIVcpz lineages.     

Inconsistencies in estimating the age of HIV-1 subtypes due to heterotachy

Rate heterogeneity among lineages is a common feature of molecular evolution, and it has long impeded our ability to accurately estimate the age of evolutionary divergence events. The development of relaxed molecular clocks, which model variable substitution rates among lineages, was intended to rectify this problem. Major subtypes of pandemic HIV-1 group M are thought to exemplify closely related lineages with different substitution rates. Here, we report that inferring the time of most recent common ancestor of all these subtypes in a single phylogeny under a single (relaxed) molecular clock produces significantly different dates for many of the subtypes than does analysis of each subtype on its own. We explore various methods to ameliorate this problem. We conclude that current molecular dating methods are inadequate for dealing with this type of substitution rate variation in HIV-1. Through simulation, we show that heterotachy causes root ages to be overestimated.     

Dating the Age of the SIV Lineages That Gave Rise to HIV-1 and HIV-2

Great strides have been made in understanding the evolutionary history of simian immunodeficiency virus (SIV) and the zoonoses that gave rise to HIV-1 and HIV-2. What remains unknown is how long these SIVs had been circulating in non-human primates before the transmissions to humans. Here, we use relaxed molecular clock dating techniques to estimate the time of most recent common ancestor for the SIVs infecting chimpanzees and sooty mangabeys, the reservoirs of HIV-1 and HIV-2, respectively. The date of the most recent common ancestor of SIV in chimpanzees is estimated to be 1492 (1266–1685), and the date in sooty mangabeys is estimated to be 1809 (1729–1875). Notably, we demonstrate that SIV sequences sampled from sooty mangabeys possess sufficient clock-like signal to calibrate a molecular clock; despite the differences in host biology and viral dynamics, the rate of evolution of SIV in sooty mangabeys is indistinguishable from that of its human counterpart, HIV-2. We also estimate the ages of the HIV-2 human-to-human transmissible lineages and provide the first age estimate for HIV-1 group N at 1963 (1948–1977). Comparisons between the SIV most recent common ancestor dates and those of the HIV lineages suggest a difference on the order of only hundreds of years. Our results suggest either that SIV is a surprisingly young lentiviral lineage or that SIV and, perhaps, HIV dating estimates are seriously compromised by unaccounted-for biases.     

The molecular population genetics of HIV-1 group O

HIV-1 group O originated through cross-species transmission of SIV from chimpanzees to humans and has established a relatively low prevalence in Central Africa. Here, we infer the population genetics and epidemic history of HIV-1 group O from viral gene sequence data and evaluate the effect of variable evolutionary rates and recombination on our estimates. First, model selection tools were used to specify suitable evolutionary and coalescent models for HIV group O. Second, divergence times and population genetic parameters were estimated in a Bayesian framework using Markov chain Monte Carlo sampling, under both strict and relaxed molecular clock methods. Our results date the origin of the group O radiation to around 1920 (1890-1940), a time frame similar to that estimated for HIV-1 group M. However, group O infections, which remain almost wholly restricted to Cameroon, show a slower rate of exponential growth during the twentieth century, explaining their lower current prevalence. To explore the effect of recombination, the Bayesian framework is extended to incorporate multiple unlinked loci. Although recombination can bias estimates of the time to the most recent common ancestor, this effect does not appear to be important for HIV-1 group O. In addition, we show that evolutionary rate estimates for different HIV genes accurately reflect differential selective constraints along the HIV genome.     

Adaptation of HIV-1 to its human host

Human immunodeficiency virus type 1 (HIV-1) originated from three independent cross-species transmissions of simian immunodeficiency virus (SIVcpzPtt) infecting chimpanzees (Pan troglodytes troglodytes) in west central Africa, giving rise to pandemic (group M) and non-pandemic (groups N and O) clades of HIV-1. To identify host-specific adaptations in HIV-1 we compared the inferred ancestral sequences of HIV-1 groups M, N and O to 12 full length genome sequences of SIVcpzPtt and four of the outlying but closely related SIVcpzPts (from P. t. schweinfurthii). This analysis revealed a single site that was completely conserved among SIVcpzPtt strains but different (due to the same change) in all three groups of HIV-1. This site, Gag-30, lies within p17, the gag-encoded matrix protein. It is Met in SIVcpzPtt, underwent a conservative replacement by Leu in one lineage of SIVcpzPts but changed radically to Arg on all three lineages leading to HIV-1. During subsequent diversification this site has been conserved as a basic residue (Arg or Lys) in most lineages of HIV-1. Retrospective analysis revealed that Gag-30 had reverted to Met in a previous experiment in which HIV-1 was passaged through chimpanzees. To examine whether this substitution conferred a species specific growth advantage, we used site-directed mutagenesis to generate variants of these chimpanzee-adapted HIV-1 strains with Lys at Gag-30, and tested their replication in both human and chimpanzee CD4+ T lymphocytes. Remarkably, viruses encoding Met replicated to higher titers than viruses encoding Lys in chimpanzee T cells, but the opposite was found in human T cells. Taken together, these observations provide compelling evidence for host-specific adaptation during the emergence of HIV-1 and identify the viral matrix protein as a modulator of viral fitness following transmission to the new human host.     

Molecular epidemiology of simian immunodeficiency virus SIVsm in U.S. primate centers unravels the origin of SIVmac and SIVstm

Retrospective molecular epidemiology was performed on samples from four sooty mangabey (SM) colonies in the United States to characterize simian immunodeficiency virus SIVsm diversity in SMs and to trace virus circulation among different primate centers (PCs) over the past 30 years. The following SIVsm sequences were collected from different monkeys: 55 SIVsm isolates from the Tulane PC sampled between 1984 and 2004, 10 SIVsm isolates from the Yerkes PC sampled in 2002, 7 SIVsm isolates from the New Iberia PC sampled between 1979 and 1986, and 8 SIVsm isolates from the California PC sampled between 1975 and 1977. PCR and sequencing were done to characterize the gag, pol, and env gp36 genes. Phylogenetic analyses were correlated with the epidemiological data. Our analysis identified nine different divergent phylogenetic lineages that cocirculated in these four SM colonies in the Unites States in the past 30 years. Lineages 1 to 5 have been identified previously. Two of the newly identified SIVsm lineages found in SMs are ancestral to SIVmac251/SIVmac239/SIVmne and SIVstm. We further identified the origin of these two macaque viruses in SMs from the California National Primate Research Center. The diversity of SIVsm isolates in PCs in the United States mirrors that of human immunodeficiency virus type 1 (HIV-1) group M subtypes and offers a model for the molecular epidemiology of HIV and a new approach to vaccine testing. The cocirculation of divergent SIVsm strains in PCs resulted in founder effects, superinfections, and recombinations. This large array of SIVsm strains showing the same magnitude of diversity as HIV-1 group M subtypes should be extremely useful for modeling the efficacy of vaccination strategies under the real-world conditions of HIV-1 diversity. The genetic variability of SIVsm strains among PCs may influence the diagnosis and monitoring of SIVsm infection and, consequently, may bias the results of pathogenesis studies.     

Cyclophilin A is required for the replication of group M human immunodeficiency virus type 1 (HIV-1) and simian immunodeficiency virus SIV(CPZ)GAB but not group O HIV-1 or other primate immunodeficiency viruses.

The human immunodeficiency virus type 1 (HIV-1) Gag polyprotein binds to cyclophilin A and incorporates this cellular peptidyl prolyl-isomerase into virions. Disruption of cyclophilin A incorporation, either by gag mutations or by cyclosporine A, inhibits virion infectivity, indicating that cyclophilin A plays an essential role in the HIV-1 life cycle. Using assays for packaging of cyclophilin A into virions and for viral replication sensitivity to cyclosporine A, as well as information gleaned from the alignment of Gag residues encoded by representative viral isolates, we demonstrate that of the five lineages of primate immunodeficiency viruses, only HIV-1 requires cyclophilin A for replication. Cloned viral isolates from clades A, B, and D of HIV-1 group M, as well as a phylogenetically related isolate from chimpanzee, all require cyclophilin A for replication. In contrast, the replication of two outlier (group O) HIV-1 isolates is unaffected by concentrations of cyclosporine A which disrupt cyclophilin A incorporation into virions, indicating that these viruses are capable of replicating independently of cyclophilin A. These studies identify the first phenotypic difference between HIV-1 group M and group O and are consistent with phylogenetic studies suggesting that the two HIV-1 groups were introduced into human populations via separate zoonotic transmission events.     

Estimating the Rate of Intersubtype Recombination in Early HIV-1 Group M Strains

West Central Africa has been implicated as the epicenter of the HIV-1 epidemic, and almost all group M subtypes can be found there. Previous analysis of early HIV-1 group M sequences from Kinshasa in the Democratic Republic of Congo, formerly Zaire, revealed that isolates from a number of individuals fall in different positions in phylogenetic trees constructed from sequences from opposite ends of the genome as a result of recombination between viruses of different subtypes. Here, we use discrete ancestral trait mapping to develop a procedure for quantifying HIV-1 group M intersubtype recombination across phylogenies, using individuals'' gag (p17) and env (gp41) subtypes. The method was applied to previously described HIV-1 group M sequences from samples obtained in Kinshasa early in the global radiation of HIV. Nine different p17 and gp41 intersubtype recombinant combinations were present in the data set. The mean number of excess ancestral subtype transitions (NEST) required to map individuals'' p17 subtypes onto the gp14 phylogeny samples, compared to the number required to map them onto the p17 phylogenies, and vice versa, indicated that excess subtype transitions occurred at a rate of approximately 7 × 10⁻³ to 8 × 10⁻³ per lineage per year as a result of intersubtype recombination. Our results imply that intersubtype recombination may have occurred in approximately 20% of lineages evolving over a period of 30 years and confirm intersubtype recombination as a substantial force in generating HIV-1 group M diversity.     

Multifractal analysis of HIV-1 genomes

Pathogens like HIV-1, which evolve into many closely related variants displaying differential infectivity and evolutionary dynamics in a short time scale, require fast and accurate classification. Conventional whole genome sequence alignment-based methods are computationally expensive and involve complex analysis. Alignment-free methodologies are increasingly being used to effectively differentiate genomic variations between viral species. Multifractal analysis, which explores the self-similar nature of genomes, is an alignment-free methodology that has been applied to study such variations. However, whether multifractal analysis can quantify variations between closely related genomes, such as the HIV-1 subtypes, is an open question. Here we address the above by implementing the multifractal analysis on four retroviral genomes (HIV-1, HIV-2, SIVcpz, and HTLV-1), and demonstrate that individual multifractal properties can differentiate between different retrovirus types easily. However, the individual multifractal measures do not resolve within-group variations for different known subtypes of HIV-1 M group. We show here that these known subtypes can instead be classified correctly using a combination of the crucial multifractal measures. This method is simple and computationally fast in comparison to the conventional alignment-based methods for whole genome phylogenetic analysis.     

Evolutionary modeling of rate shifts reveals specificity determinants in HIV-1 subtypes

A hallmark of the human immunodeficiency virus 1 (HIV-1) is its rapid rate of evolution within and among its various subtypes. Two complementary hypotheses are suggested to explain the sequence variability among HIV-1 subtypes. The first suggests that the functional constraints at each site remain the same across all subtypes, and the differences among subtypes are a direct reflection of random substitutions, which have occurred during the time elapsed since their divergence. The alternative hypothesis suggests that the functional constraints themselves have evolved, and thus sequence differences among subtypes in some sites reflect shifts in function. To determine the contribution of each of these two alternatives to HIV-1 subtype evolution, we have developed a novel Bayesian method for testing and detecting site-specific rate shifts. The RAte Shift EstimatoR (RASER) method determines whether or not site-specific functional shifts characterize the evolution of a protein and, if so, points to the specific sites and lineages in which these shifts have most likely occurred. Applying RASER to a dataset composed of large samples of HIV-1 sequences from different group M subtypes, we reveal rampant evolutionary shifts throughout the HIV-1 proteome. Most of these rate shifts have occurred during the divergence of the major subtypes, establishing that subtype divergence occurred together with functional diversification. We report further evidence for the emergence of a new sub-subtype, characterized by abundant rate-shifting sites. When focusing on the rate-shifting sites detected, we find that many are associated with known function relating to viral life cycle and drug resistance. Finally, we discuss mechanisms of covariation of rate-shifting sites.     

Human Immunodeficiency Virus Type 1 Intergroup (M/O) Recombination in Cameroon

Here we describe, for the first time, recombinants between two highly divergent major groups of human immunodeficiency virus type 1 (HIV-1), M and O, within a Cameroonian woman infected with three different HIV-1 strains, a group O virus, a subtype D virus, and a recently reported IBNG (A/G)-like recombinant virus. Using nested extra-long PCR amplification, we sequenced from the pol region to the env region including accessory genes of the viral genome obtained from the patient's uncultured peripheral blood mononuclear cells and examined the phylogenetic position of each gene. Compared with sequential blood samples obtained in 1995 and 1996, there were multiple segmental exchanges between three HIV-1 strains (O, D, and IBNG) and all the recombinants appeared to be derived from a common M/O ancestor. Importantly, recombination between groups M and O occurred, even though the homology between these two groups is 69, 76, 68, and 55% in the gag, pol, vif-vpr, and env regions, respectively. Recombination between strains with such distant lineages may contribute substantially to generating new HIV-1 variants.     

Characterization of a Highly Replicative Intergroup M/O Human Immunodeficiency Virus Type 1 Recombinant Isolated from a Cameroonian Patient

A Cameroonian patient with antibodies reacting simultaneously to human immunodeficiency virus type 1 (HIV-1) group O- and group M-specific V3-loop peptides was identified. In order to confirm that this patient was coinfected with both viruses, PCRs with O- and M-specific discriminating primers corresponding to different regions of the genome were carried out with both primary lymphocyte DNA and the corresponding viral strains isolated from three consecutive patient samples. The PCR data suggested that this patient is coinfected with a group M virus and a recombinant M/O virus. Indeed, only type M gag sequences could be amplified, while for the env region, both type M and O sequences were amplified, from plasma or from DNA extracted from primary lymphocytes. Sequence analysis of a complete recombinant genome isolated from the second sample (97CA-MP645 virus isolate) revealed two intergroup breakpoints, one in the vpr gene and the second in the long terminal repeat region around the TATA box. Comparison of the type M sequences shared by the group M and the recombinant M/O viruses showed that these sequences were closely related, with only 3% genetic distance, suggesting that the M virus was one of the parental viruses. In this report we describe for the first time a recombination event in vivo between viruses belonging to two different groups, leading to a replicative virus. Recombination between strains with such distant lineages (65% overall homology) may contribute substantially to the emergence of new HIV-1 variants. We documented that this virus replicates well and became predominant in vitro. At this time, group O viruses represent a minority of the strains responsible for the HIV-1 pandemic. If such recombinant intergroup viruses gained better fitness, inducing changes in their biological properties compared to the parental group O virus, the prevalences of group O sequences could increase rapidly. This will have important implications for diagnosis of HIV-1 infections by serological and molecular tests, as well as for antiviral treatment.     

Multiple barriers to recombination between divergent HIV-1 variants revealed by a dual-marker recombination assay

Recombination is a major force for generating human immunodeficiency virus type 1 (HIV-1) diversity and produces numerous recombinants circulating in the human population. We previously established a cell-based system using green fluorescent protein gene (gfp) as a reporter to study the mechanisms of HIV-1 recombination. We now report an improved system capable of detecting recombination using authentic viral sequences. Frameshift mutations were introduced into the gag gene so that parental viruses do not express full-length Gag; however, recombination can generate a progeny virus that expresses a functional Gag. We demonstrate that this Gag reconstitution assay can be used to detect recombination between two group M HIV-1 variants of the same or of different subtypes. Using both gfp and gag assays, we found that, similar to group M viruses, group O viruses also recombine frequently. When recombination between a group M virus and a group O virus was examined, we found three distinct barriers for intergroup recombination. First, similar to recombination within group M viruses, intergroup recombination is affected by the identity of the dimerization initiation signal (DIS); variants with the same DIS recombined at a higher rate than those with different DIS. Second, using the gfp recombination assay, we showed that intergroup recombination occurs much less frequently than intragroup recombination, even though the gfp target sequence is identical in all viruses. Finally, Gag reconstitution between variants from different groups is further reduced compared with green fluorescent protein, indicating that sequence divergence interferes with recombination efficiency in the gag gene. Compared with identical sequences, we estimate that recombination rates are reduced by 3-fold and by 10- to 13-fold when the target regions in gag contain 91% and 72-73% sequence identities, respectively. These results show that there are at least three distinct mechanisms preventing exchange of genetic information between divergent HIV-1 variants from different groups.     

The evolution of HIV-1 and the origin of AIDS

The major cause of acquired immune deficiency syndrome (AIDS) is human immunodeficiency virus type 1 (HIV-1). We have been using evolutionary comparisons to trace (i) the origin(s) of HIV-1 and (ii) the origin(s) of AIDS. The closest relatives of HIV-1 are simian immunodeficiency viruses (SIVs) infecting wild-living chimpanzees (Pan troglodytes troglodytes) and gorillas (Gorilla gorilla gorilla) in west central Africa. Phylogenetic analyses have revealed the origins of HIV-1: chimpanzees were the original hosts of this clade of viruses; four lineages of HIV-1 have arisen by independent cross-species transmissions to humans and one or two of those transmissions may have been via gorillas. However, SIVs are primarily monkey viruses: more than 40 species of African monkeys are infected with their own, species-specific, SIV and in at least some host species, the infection seems non-pathogenic. Chimpanzees acquired from monkeys two distinct forms of SIVs that recombined to produce a virus with a unique genome structure. We have found that SIV infection causes CD4⁺ T-cell depletion and increases mortality in wild chimpanzees, and so the origin of AIDS is more ancient than the origin of HIV-1. Tracing the genetic changes that occurred as monkey viruses adapted to infect first chimpanzees and then humans may provide insights into the causes of the pathogenicity of these viruses.     

The replicative fitness of primary human immunodeficiency virus type 1 (HIV-1) group M, HIV-1 group O, and HIV-2 isolates

The main (M) group of human immunodeficiency virus type 1 (HIV-1) is responsible for the global AIDS epidemic while HIV-1 group O (outlier) and HIV type 2 are endemic only in west and central Africa. The failure of HIV-2 and especially HIV-1 group O to spread following the initial zoonotic jumps is not well understood. This study was designed to examine the relative replicative capacities between these human lentiviruses. A pairwise competition experiment was performed with peripheral blood mononuclear cells with eight HIV-2 isolates, 6 group O viruses, and 15 group M viruses of subtype A (2 viruses), B (5 viruses), C (4 viruses), D (2 viruses) and CRF01_AE (2 viruses). HIV-1 group M isolates of any subtype were typically 100-fold-more fit than group O or HIV-2 strains when competed in peripheral blood mononuclear cells from various humans. This order in replicative fitness was also observed when virus pairs were added to human dendritic cells and then cocultured with primary, quiescent T cells, which is the model for HIV-1 transmission. These results suggest that reduced replicative and transmission fitness may be contributing to the low prevalence and limited geographical spread of HIV-2 and group O HIV-1 in the human population.     

SIVcpz from a naturally infected Cameroonian chimpanzee: Biological and genetic comparison with HIV-1 N

Thus far, simian immunodeficiency virus from chimpanzees (SIVcpz) genomes have been characterized as Pan troglodytes troglodytes and show a strong relation with human immunodeficiency virus (HIV)-1 N in their env genes. We fully characterized another SIVcpz from P. t. troglodytes . This chimpanzee (Cam5) was, as was also the host of SIVcpz-cam3, wild born in Cameroon, a region where all three groups of HIV-1 (M, N and O) co-occur. In contrast to other SIVcpz, SIVcpz-cam5 was isolated immediately after the rescue of the animal. Our data demonstrate that SIVcpz-cam5, like SIVcpz-cam3, grows easily on human peripheral blood mononuclear cells (PBMCs) and uses CCR5 as a co-receptor similar to HIV-1 N YBF30. Phylogenetic analysis based on the entire env gene shows that SIVcpz-cam5 falls into the same unique subcluster as HIV-1 N YBF30, SIVcpz-cam3 and SIVcpz-US. A phylogenetic relationship was also found with the vif gene of HIV-1 N. This study provides proof that HIV-1 N related viruses circulate in wild P. t. troglodytes .     

SIVcpz from a naturally infected Cameroonian chimpanzee: biological and genetic comparison with HIV-1 N

Thus far, simian immunodeficiency virus from chimpanzees (SIVcpz) genomes have been characterized as Pan troglodytes troglodytes and show a strong relation with human immunodeficiency virus (HIV)-1 N in their env genes. We fully characterized another SIVcpz from P. t. troglodytes. This chimpanzee (Cam5) was, as was also the host of SIVcpz-cam3, wild born in Cameroon, a region where all three groups of HIV-1 (M, N and O) co-occur. In contrast to other SIVcpz, SIVcpz-cam5 was isolated immediately after the rescue of the animal. Our data demonstrate that SIVcpz-cam5, like SIVcpz-cam3, grows easily on human peripheral blood mononuclear cells (PBMCs) and uses CCR5 as a co-receptor similar to HIV-1 N YBF30. Phylogenetic analysis based on the entire env gene shows that SIVcpz-cam5 falls into the same unique subcluster as HIV-1 N YBF30, SIVcpz-cam3 and SIVcpz-US. A phylogenetic relationship was also found with the vif gene of HIV-1 N. This study provides proof that HIV-1 N related viruses circulate in wild P. t. troglodytes.     

Mitochondrial complex I activity is impaired during HIV-1-induced T-cell apoptosis

Studies carried out till date to elucidate the pathways involved in HIV-1-induced T-cell depletion has revealed that apoptosis underlie the etiology, however, a clear molecular understanding of HIV-1-induced apoptosis has remained elusive. Although evidences pointing towards the importance of mitochondrial energy generating system in apoptosis exist but it's exact role remains to be clearly understood. Here, we describe for the first time specific downregulation of a complex I subunit NDUFA6 with simultaneous impairment of mitochondrial complex I activity in HIV infection. We also show that NDUFA6 gene silencing induces apoptosis and its overexpression reduces apoptosis in HIV-infected cells. Finally, sensitivity to complex I inhibitor Rotenone is reduced in HIV-1-infected T cells indicating an important role for it in the death process. Our data provide a novel molecular basis as to how the virus might interfere with host cell energy generating system during apoptotic cell death.     

New insights into the evolutionary rate of HIV-1 at the within-host and epidemiological levels

Over calendar time, HIV-1 evolves considerably faster within individuals than it does at the epidemic level. This is a surprising observation since, from basic population genetic theory, we would expect the genetic substitution rate to be similar across different levels of biological organization. Three different mechanisms could potentially cause the observed mismatch in phylogenetic rates of divergence: temporal changes in selection pressure during the course of infection; frequent reversion of adaptive mutations after transmission; and the storage of the virus in the body followed by the preferential transmission of stored ancestral virus. We evaluate each of these mechanisms to determine whether they are likely to make a major contribution to the mismatch in phylogenetic rates. We conclude that the cycling of the virus through very long-lived memory CD4(+) T cells, a process that we call 'store and retrieve', is probably the major contributing factor to the rate mismatch. The preferential transmission of ancestral virus needs to be integrated into evolutionary models if we are to accurately predict the evolution of immune escape, drug resistance and virulence in HIV-1 at the population level. Moreover, early infection viruses should be the major target for vaccine design, because these are the viral strains primarily involved in transmission.     

Reviewing the history of HIV-1: spread of subtype B in the Americas

The dispersal of HIV-1 subtype B (HIV-1B) is a reflection of the movement of human populations in response to social, political, and geographical issues. The initial dissemination of HIV-1B outside Africa seems to have included the passive involvement of human populations from the Caribbean in spreading the virus to the United States. However, the exact pathways taken during the establishment of the pandemic in the Americas remain unclear. Here, we propose a geographical scenario for the dissemination of HIV-1B in the Americas, based on phylogenetic and genetic statistical analyses of 313 available sequences of the pol gene from 27 countries. Maximum likelihood and bayesian inference methods were used to explore the phylogenetic relationships between HIV-1B sequences, and molecular variance estimates were analyzed to infer the genetic structure of the viral population. We found that the initial dissemination and subsequent spread of subtype B in the Americas occurred via a single introduction event in the Caribbean around 1964 (1950-1967). Phylogenetic trees present evidence of several primary outbreaks in countries in South America, directly seeded by the Caribbean epidemic. Cuba is an exception insofar as its epidemic seems to have been introduced from South America. One clade comprising isolates from different countries emerged in the most-derived branches, reflecting the intense circulation of the virus throughout the American continents. Statistical analysis supports the genetic compartmentalization of the virus among the Americas, with a close relationship between the South American and Caribbean epidemics. These findings reflect the complex establishment of the HIV-1B pandemic and contribute to our understanding between the migration process of human populations and virus diffusion.