Reconstructing the Timing and Dispersion Routes of HIV-1 Subtype B Epidemics in The Caribbean and Central America: A Phylogenetic Story

The Caribbean and Central America are among the regions with highest HIV-1B prevalence worldwide. Despite of this high virus burden, little is known about the timing and the migration patterns of HIV-1B in these regions. Migration is one of the major processes shaping the genetic structure of virus populations. Thus, reconstruction of epidemiological network may contribute to understand HIV-1B evolution and reduce virus prevalence. We have investigated the spatio-temporal dynamics of the HIV-1B epidemic in The Caribbean and Central America using 1,610 HIV-1B partial pol sequences from 13 Caribbean and 5 Central American countries. Timing of HIV-1B introduction and virus evolutionary rates, as well as the spatial genetic structure of the HIV-1B populations and the virus migration patterns were inferred. Results revealed that in The Caribbean and Central America most of the HIV-1B variability was generated since the 80 s. At odds with previous data suggesting that Haiti was the origin of the epidemic in The Caribbean, our reconstruction indicated that the virus could have been disseminated from Puerto Rico and Antigua. These two countries connected two distinguishable migration areas corresponding to the (mainly Spanish-colonized) Easter and (mainly British-colonized) Western islands, which indicates that virus migration patterns are determined by geographical barriers and by the movement of human populations among culturally related countries. Similar factors shaped the migration of HIV-1B in Central America. The HIV-1B population was significantly structured according to the country of origin, and the genetic diversity in each country was associated with the virus prevalence in both regions, which suggests that virus populations evolve mainly through genetic drift. Thus, our work contributes to the understanding of HIV-1B evolution and dispersion pattern in the Americas, and its relationship with the geography of the area and the movements of human populations.     

Human Immunodeficiency Virus Type 1 (HIV-1) Subtype B Epidemic in Panama Is Mainly Driven by Dissemination of Country-Specific Clades

The Human immunodeficiency virus type-1 (HIV-1) subtype B is the most predominant clade in Central America; but information about the evolutionary history of this virus in this geographic region is scarce. In this study, we reconstructed the spatiotemporal and population dynamics of the HIV-1 subtype B epidemic in Panama. A total of 761 HIV-1 subtype B pol sequences obtained in Panama between 2004 and 2013 were combined with subtype B pol sequences from the Americas and Europe. Maximum Likelihood phylogenetic analyses revealed that HIV-1 subtype B infections in Panama derived from the dissemination of multiple founder viruses. Most Panamanian subtype B viruses (94.5%) belong to the pandemic viral strain proposed as originated in the US, whereas others (5.5%) were intermixed among non-pandemic Caribbean strains. The bulk (76.6%) of subtype B sequences from Panama grouped within 12 country-specific clades that were not detected in other Central American countries. Bayesian coalescent-based analyses suggest that most Panamanian clades probably originated between the early 1970s and the early 1980s. The root location of major Panamanian clades was traced to the most densely populated districts of Panama province. Major Panamanian clades appear to have experienced one or two periods of exponential growth of variable duration between the 1970s and the 2000s, with median growth rates from 0.2 to 0.4 year^{− 1}. Thus, the HIV-1 subtype B epidemic in Panama is driven by the expansion of local viral strains that were introduced from the Caribbean and other American countries at an early stage of the AIDS pandemic.     

Spatiotemporal Dynamics of Dissemination of Non-Pandemic HIV-1 Subtype B Clades in the Caribbean Region

The Human immunodeficiency virus type-1 (HIV-1) epidemic in the Caribbean region is mostly driven by subtype B; but information about the pattern of viral spread in this geographic region is scarce and different studies point to quite divergent models of viral dissemination. In this study, we reconstructed the spatiotemporal and population dynamics of the HIV-1 subtype B epidemic in the Caribbean. A total of 1,806 HIV-1 subtype B pol sequences collected from 17 different Caribbean islands between 1996 and 2011 were analyzed together with sequences from the United States (n = 525) and France (n = 340) included as control. Maximum Likelihood phylogenetic analyses revealed that HIV-1 subtype B infections in the Caribbean are driven by dissemination of the pandemic clade (B_PANDEMIC) responsible for most subtype B infections across the world, and older non-pandemic lineages (B_CAR) characteristics of the Caribbean region. The non-pandemic B_CAR strains account for >40% of HIV-1 infections in most Caribbean islands; with exception of Cuba and Puerto Rico. Bayesian phylogeographic analyses indicate that B_CAR strains probably arose in the island of Hispaniola (Haiti/Dominican Republic) around the middle 1960s and were later disseminated to Trinidad and Tobago and to Jamaica between the late 1960s and the early 1970s. In the following years, the B_CAR strains were also disseminated from Hispaniola and Trinidad and Tobago to other Lesser Antilles islands at multiple times. The B_CAR clades circulating in Hispaniola, Jamaica and Trinidad and Tobago appear to have experienced an initial phase of exponential growth, with mean estimated growth rates of 0.35–0.45 year⁻¹, followed by a more recent stabilization since the middle 1990s. These results demonstrate that non-pandemic subtype B lineages have been widely disseminated through the Caribbean since the late 1960s and account for an important fraction of current HIV-1 infections in the region.     

A Single Early Introduction of HIV-1 Subtype B into Central America Accounts for Most Current Cases

Human immunodeficiency virus type 1 (HIV-1) variants show considerable geographical separation across the world, but there is limited information from Central America. We provide the first detailed investigation of the genetic diversity and molecular epidemiology of HIV-1 in six Central American countries. Phylogenetic analysis was performed on 625 HIV-1 pol gene sequences collected between 2002 and 2010 in Honduras, El Salvador, Nicaragua, Costa Rica, Panama, and Belize. Published sequences from neighboring countries (n = 57) and the rest of the world (n = 740) were included as controls. Maximum likelihood methods were used to explore phylogenetic relationships. Bayesian coalescence-based methods were used to time HIV-1 introductions. Nearly all (98.9%) Central American sequences were of subtype B. Phylogenetic analysis revealed that 437 (70%) sequences clustered within five significantly supported monophyletic clades formed essentially by Central American sequences. One clade contained 386 (62%) sequences from all six countries; the other four clades were smaller and more country specific, suggesting discrete subepidemics. The existence of one large well-supported Central American clade provides evidence that a single introduction of HIV-1 subtype B in Central America accounts for most current cases. An introduction during the early phase of the HIV-1 pandemic may explain its epidemiological success. Moreover, the smaller clades suggest a subsequent regional spread related to specific transmission networks within each country.     

Phylodynamics of HIV-1 subtype B among the men-having-sex-with-men (MSM) population in Hong Kong

The men-having-sex-with-men (MSM) population has become one of the major risk groups for HIV-1 infection in the Asia Pacific countries. Hong Kong is located in the centre of Asia and the transmission history of HIV-1 subtype B transmission among MSM remained unclear. The aim of this study was to investigate the transmission dynamics of HIV-1 subtype B virus in the Hong Kong MSM population. Samples of 125 HIV-1 subtype B infected MSM patients were recruited in this study. Through this study, the subtype B epidemic in the Hong Kong MSM population was identified spreading mainly among local Chinese who caught infection locally. On the other hand, HIV-1 subtype B infected Caucasian MSM caught infection mainly outside Hong Kong. The Bayesian phylogenetic analysis also indicated that 3 separate subtype B epidemics with divergence dates in the 1990s had occurred. The first and latest epidemics were comparatively small-scaled; spreading among the local Chinese MSM while sauna-visiting was found to be the major sex partner sourcing reservoir for the first subtype B epidemic. However, the second epidemic was spread in a large-scale among local Chinese MSM with a number of them having sourced their sex partners through the internet. The epidemic virus was estimated to have a divergence date in 1987 and the infected population in Hong Kong had a logistic growth throughout the past 20 years. Our study elucidated the evolutionary and demographic history of HIV-1 subtype B virus in Hong Kong MSM population. The understanding of transmission and growth model of the subtype B epidemic provides more information on the HIV-1 transmission among MSM population in other Asia Pacific high-income countries.     

Evolution of the human immunodeficiency virus type 1 subtype-specific V3 domain is confined to a sequence space with a fixed distance to the subtype consensus.

Human immunodeficiency virus type 1 (HIV-1) strains can be separated into genetic subtypes based on phylogenetic analysis of the envelope gene. Once it had been shown that population-wide intrasubtype genetic variation of HIV-1 strains increases in the course of the AIDS epidemic, it remained uncertain whether HIV-1 subtypes are phenotypic entities spreading as distinct virus populations. To examine this, we applied Eigen's concepts of sequence geometry and fitness topography to the analysis of intrasubtype evolution of the gp120 V3 domain of HIV-1 subtypes A, B, C, and D in the course of the global AIDS epidemic. We observed that despite the high evolution rate of HIV-1, the nonsynonymous distances to the subtype consensus of sequences obtained early in the epidemic are similar to those obtained more than 10 years later, in contrast to the synonymous distances, which increased steadily over time. For HIV-1 subtype B, we observed that the evolution rate of the individual sequences is independent of their distance from the subtype B consensus, but for the individual sequences most distant from the consensus evolution away from the consensus is constrained. As a result, individual HIV-1 genomes fluctuate within a sequence space with fixed distance to the subtype consensus. Our findings suggest that the evolution of the V3 domain of HIV-1 subtypes A, B, C, and D is confined to an area in sequence space within a fixed distance to the consensus of a respective subtype. This in turn indicates that each HIV-1 subtype is a distinct viral quasispecies that is well adapted to the present environment, able to maintain its identity in the V3 region over time, and unlikely to merge during progression of the AIDS epidemic.     

HIV-1 genetic diversity in Western Brittany,France

We describe human immunodeficiency virus type 1 (HIV-1) diversity in Western Brittany, France, and trace the dissemination of HIV-1 non-B subtype infection. The strategy for HIV-1 subtyping used involved subtype specific enzyme immunoassays, heteroduplex mobility assays and phylogenetic analysis of the sequences of env encoding the V3 loop region. Samples were obtained from 567 patients: 465 (82%) were of subtype B and 66 (11.6%) were not (20 were subtype A, 11 subtype C, four subtype D, seven subtype F, five subtype G and 19 others with circulating recombinant forms: 4CRF01_AE, 11CRF02_AG, 1H, 3CRF11_cpx). These findings are consistent with other studies of French populations. There is an epidemiological correlation between subtype B and homosexual or heterosexual contamination in France and between non-B subtype and heterosexual contamination in Africa.     

Identification of reassortant pandemic H1N1 influenza virus in Korean pigs

Since the 2009 pandemic human H1N1 influenza A virus emerged in April 2009, novel reassortant strains have been identified throughout the world. This paper describes the detection and isolation of reassortant strains associated with human pandemic influenza H1N1 and swine influenza H1N2 (SIV) viruses in swine populations in South Korea. Two influenza H1N2 reassortants were detected, and subtyped by PCR. The strains were isolated using Madin- Darby canine kidney (MDCK) cells, and genetically characterized by phylogenetic analysis for genetic diversity. They consisted of human, avian, and swine virus genes that were originated from the 2009 pandemic H1N1 virus and a neuraminidase (NA) gene from H1N2 SIV previously isolated in North America. This identification of reassortment events in swine farms raises concern that reassortant strains may continuously circulate within swine populations, calling for the further study and surveillance of pandemic H1N1 among swine.     

Phylogeographic analysis of HIV-1 subtype C dissemination in Southern Brazil

The HIV-1 subtype C has spread efficiently in the southern states of Brazil (Rio Grande do Sul, Santa Catarina and Paraná). Phylogeographic studies indicate that the subtype C epidemic in southern Brazil was initiated by the introduction of a single founder virus population at some time point between 1960 and 1980, but little is known about the spatial dynamics of viral spread. A total of 135 Brazilian HIV-1 subtype C pol sequences collected from 1992 to 2009 at the three southern state capitals (Porto Alegre, Florianópolis and Curitiba) were analyzed. Maximum-likelihood and Bayesian methods were used to explore the degree of phylogenetic mixing of subtype C sequences from different cities and to reconstruct the geographical pattern of viral spread in this country region. Phylogeographic analyses supported the monophyletic origin of the HIV-1 subtype C clade circulating in southern Brazil and placed the root of that clade in Curitiba (Paraná state). This analysis further suggested that Florianópolis (Santa Catarina state) is an important staging post in the subtype C dissemination displaying high viral migration rates from and to the other cities, while viral flux between Curitiba and Porto Alegre (Rio Grande do Sul state) is very low. We found a positive correlation (r(2) = 0.64) between routine travel and viral migration rates among localities. Despite the intense viral movement, phylogenetic intermixing of subtype C sequences from different Brazilian cities is lower than expected by chance. Notably, a high proportion (67%) of subtype C sequences from Porto Alegre branched within a single local monophyletic sub-cluster. These results suggest that the HIV-1 subtype C epidemic in southern Brazil has been shaped by both frequent viral migration among states and in situ dissemination of local clades.     

HIV-1 subtype B epidemic and transmission patterns in Slovenia

In the present study the epidemic of human immunodeficiency virus type 1 (HIV-1) subtype B in Slovenia during the 10-year period was investigated using phylogenetic analysis of pol gene sequences. 119 pol sequences generated on samples dated from January 1996 to December 2005 were retrieved from the database of Slovenian HIV/AIDS Reference Laboratory. The phylogenetic analysis revealed 14 potentially significant transmission clusters (bootstrap value > or = 98%), comprising 34 HIV-1 strains. The vast majority of clustered individuals were men (91%), and of them, 79% were men who have sex with men. Factors significantly associated with clustering were: recent infection (HIV-1 infection during or after year 2003), diagnosis of primary HIV-1 infection, higher CD4 cell count and acquiring HIV-1 infection in Slovenia. Recent subtype B HIV-1 infections are the important driving force of current HIV-1 epidemic in Slovenia.     

Complete genomic sequence of a dengue type 2 virus from the French West Indies

Severe forms of dengue fever, dengue haemorrhagic fever, and dengue shock syndrome, were not prominent in the Americas until the epidemic of Cuba in 1981. Since that time, they have spread to other countries in Central and South America, correlating with the spread of dengue type 2 viruses related to Southeast Asian strains. We report here the complete genomic sequence of a dengue type 2 virus isolated during the epidemic in La Martinique in 1998. This constitutes the first complete genetic characterization of a dengue virus strain from French West Indies, and also the first molecular identification in this region of a dengue 2 strain phylogenetically related to the emerging American type 2 dengue viruses.     

Low Genetic Diversity in Melanaphis sacchari Aphid Populations at the Worldwide Scale

Numerous studies have examined the genetic diversity and genetic structure of invading species, with contrasting results concerning the relative roles of genetic diversity and phenotypic plasticity in the success of introduced populations. Increasing evidence shows that asexual lineages of aphids are able to occupy a wide geographical and ecological range of habitats despite low genetic diversity. The anholocyclic aphid Melanaphis sacchari is a pest of sugarcane and sorghum which originated in the old world, was introduced into the Americas, and is now distributed worldwide. Our purpose was to assess the genetic diversity and structuring of populations of this species according to host and locality. We used 10 microsatellite markers to genotype 1333 individuals (57 samples, 42 localities, 15 countries) collected mainly on sugarcane or sorghum. Five multilocus lineages (MLL) were defined, grouping multilocus genotypes (MLG) differing by only a few mutations or scoring errors. Analysis of a 658 bp sequence of mitochondrial COI gene on 96 individuals revealed five haplotypes, with a mean divergence of only 0.19 %. The distribution of MLL appeared to be strongly influenced by geography but not by host plant. Each of the five MLL grouped individuals from (A) Africa, (B) Australia, (C) South America, the Caribbean and the Indian Ocean including East Africa, (D) USA, and (E) China. The MLL A and C, with a wide geographic distribution, matched the definition of superclone. Among aphids, M. sacchari has one of the lowest known rates of genetic diversity for such a wide geographical distribution.     

Invasion origin,rapid population expansion,and the lack of genetic structure of cotton bollworm (Helicoverpa armigera) in the Americas

In 2013, Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) was officially declared as present in Brazil and, after two years, the species was detected in the Caribbean and North America. Information on genetic features and accurate distribution of pests is the basis for agricultural protection policies. Furthermore, such knowledge is imperative to develop control strategies, understand the geographical range, and genetic patterns of this species in the Americas. Here, we carried out the widest sampling of H. armigera in the South American continent and Puerto Rico, after we estimated the diversity, demographic parameters, and genetic structure. The Internal Transcribed Spacer 1 (ITS1) nuclear marker was used to investigate the presence of putative hybrids between H. armigera and H. zea, and they were observed at a frequency of 1.5%. An ABC analysis, based in COI gene fragment, suggested Europe as the origin of South America specimens of H. armigeraand following a movement northward through the Caribbean. Three mtDNA genes and three nDNA markers revealed high genetic diversity distributed without the defined population structure of H. armigera in South America. Most of the genetic variation is within populations with a multidirectional expansion of H. armigera among morphoclimatic regions. High genetic diversity, rapid population expansion, and hybridization have implications for pest management since they suggest that adaptive alleles are spread through wide areas in South America that favor rapid local adaptation of H. armigera to new and disturbed environments (e.g., in agricultural areas).     

Molecular evolution of Zika virus as it crossed the Pacific to the Americas

Zika virus was previously considered to cause only a benign infection in humans. Studies of recent outbreaks of Zika virus in the Pacific, South America, Mexico and the Caribbean have associated the virus with severe neuropathology. Viral evolution may be one factor contributing to an apparent change in Zika disease as it spread from Southeast Asia across the Pacific to the Americas. To address this possibility, we have employed computational tools to compare the phylogeny, geography, immunology and RNA structure of Zika virus isolates from Africa, Asia, the Pacific and the Americas. In doing so, we compare and contrast methods and results for tree search and rooting of Zika virus phylogenies. In some phylogenetic analyses we find support for the hypothesis that there is a deep common ancestor between African and Asian clades (the “Asia/Africa” hypothesis). In other phylogenetic analyses, we find that Asian lineages are descendent from African lineages (the “out of Africa” hypothesis). In addition, we identify and evaluate key mutations in viral envelope protein coding and untranslated terminal RNA regions. We find stepwise mutations that have altered both immunological motif sets and regulatory sequence elements. Both of these sets of changes distinguish viruses found in Africa from those in the emergent Asia–Pacific–Americas lineage. These findings support the working hypothesis that mutations acquired by Zika virus in the Pacific and Americas contribute to changes in pathology. These results can inform experiments required to elucidate the role of viral genetic evolution in changes in neuropathology, including microcephaly and other neurological and skeletomuscular issues in infants, and Guillain‐Barré syndrome in adults.     

Evolutionary analyses on the HA gene ofpandemic H1N1/09: early findings

The HA protein is responsible for influenza virus attachment and the subsequent fusion of viral and cellular membranes. Antigenic drift is driven by an accumulation of point mutations in the HA. And, the receptor-binding specificity of HA is responsible for the host range restriction of the virus. In April 2009, large outbreaks of novel H1N1 influenza in human population were reported from North America. The pandemic H1N1 virus originated from swine influenza virus. Evolutionary process of the pandemic virus after its introduction to human population remains to be clarified. We conducted phylogenetic analyses constructing a phylogenetic tree for and calculating site-by-site selective pressures in the HA gene. Phylogenetic tree showed that pandemic viruses were not clustered clearly by their geographical location or isolation time in the phylogenetic tree. The virus has been circulating the globe extensively with multiple introductions into most geographical areas. We found 3 sites positively selected in the HA gene for pandemic H1N1 virus. Among them, position 206 is located in an antigenic site. We did not find significant negative selection on any of the receptor binding sites. The virus has been evolving under unique selective pressure.     

The HIV-1 Subtype C Epidemic in South America Is Linked to the United Kingdom

Background

The global spread of HIV-1 has been accompanied by the emergence of genetically distinct viral strains. Over the past two decades subtype C viruses, which predominate in Southern and Eastern Africa, have spread rapidly throughout parts of South America. Phylogenetic studies indicate that subtype C viruses were introduced to South America through a single founder event that occurred in Southern Brazil. However, the external route via which subtype C viruses spread to the South American continent has remained unclear.

Methodology/Principal Findings

We used automated genotyping to screen 8,309 HIV-1 subtype C pol gene sequences sampled within the UK for isolates genetically linked to the subtype C epidemic in South America. Maximum likelihood and Bayesian approaches were used to explore the phylogenetic relationships between 54 sequences identified in this screen, and a set of globally sampled subtype C reference sequences. Phylogenetic trees disclosed a robustly supported relationship between sequences from Brazil, the UK and East Africa. A monophyletic cluster comprised exclusively of sequences from the UK and Brazil was identified and dated to approximately the early 1980s using a Bayesian coalescent-based method. A sub-cluster of 27 sequences isolated from homosexual men of UK origin was also identified and dated to the early 1990s.

Conclusions

Phylogenetic, demographic and temporal data support the conclusion that the UK was a crucial staging post in the spread of subtype C from East Africa to South America. This unexpected finding demonstrates the role of diffuse international networks in the global spread of HIV-1 infection, and the utility of globally sampled viral sequence data in revealing these networks. Additionally, we show that subtype C viruses are spreading within the UK amongst men who have sex with men.     

Molecular Evolution of the HIV-1 Thai Epidemic between the Time of RV144 Immunogen Selection to the Execution of the Vaccine Efficacy Trial

The RV144 HIV-1 vaccine trial (Thailand, 2003 to 2009), using immunogens genetically matched to the regional epidemic, demonstrated the first evidence of efficacy for an HIV-1 vaccine. Here we studied the molecular evolution of the HIV-1 epidemic from the time of immunogen selection to the execution of the efficacy trial. We studied HIV-1 genetic diversity among 390 volunteers who were deferred from enrollment in RV144 due to preexisting HIV-1 infection using a multiregion hybridization assay, full-genome sequencing, and phylogenetic analyses. The subtype distribution was 91.7% CRF01_AE, 3.5% subtype B, 4.3% B/CRF01_AE recombinants, and 0.5% dual infections. CRF01_AE strains were 31% more diverse than the ones from the 1990s Thai epidemic. Sixty-nine percent of subtype B strains clustered with the cosmopolitan Western B strains. Ninety-three percent of B/CRF01_AE recombinants were unique; recombination breakpoint analysis showed that these strains were highly embedded within the larger network that integrates recombinants from East/Southeast Asia. Compared to Thai sequences from the early 1990s, the distance to the RV144 immunogens increased 52% to 68% for CRF01_AE Env immunogens and 12% to 29% for subtype B immunogens. Forty-three percent to 48% of CRF01_AE sequences differed from the sequence of the vaccine insert in Env variable region 2 positions 169 and 181, which were implicated in vaccine sieve effects in RV144. In conclusion, compared to the molecular picture at the early stages of vaccine development, our results show an overall increase in the genetic complexity of viruses in the Thai epidemic and in the distance to vaccine immunogens, which should be considered at the time of the analysis of the trial results.     

Evolution of the HIV-1 V3 region in the Italian epidemic

The epidemic of Human Immunodeficiency Virus type 1 infection in Italy is mostly ascribed to the B subtype, which represents the prevalent subtype in Western Countries. The virus isolates of the B subtype, moreover, show an increasing nucleotide heterogeneity over time, indicating a continuous intra-subtype dynamic evolution, typical of long-lasting epidemics. In recent years, however, the progressive decrease in the transmission rate among the historically defined risk groups (i.e. homosexuals and IDUs) and the parallel increase in heterosexual transmission are slowly introducing variants of non-B subtype into the Italian HIV-1 epidemic. This appears to be strictly linked to the growing number of immigrants from non-Western Countries, where non-B clades and CRFs are prevalent, and consequent inter-racial blending. The distribution of these novel genetic forms needs to be evaluated by continuous molecular monitoring nationwide to verify whether they will overcome the pre-existing B-clade epidemic, which could have significant implications for diagnosis, treatment and vaccine development. Here we review the genetic evolution of HIV-1 spreading within the Italian epidemic.     

Molecular epidemiology of 58 new African human T-cell leukemia virus type 1 (HTLV-1) strains: identification of a new and distinct HTLV-1 molecular subtype in Central Africa and in Pygmies.

To gain new insights on the origin, evolution, and modes of dissemination of human T-cell leukemia virus type I (HTLV-1), we performed a molecular analysis of 58 new African HTLV-1 strains (18 from West Africa, 36 from Central Africa, and 4 from South Africa) originating from 13 countries. Of particular interest were eight strains from Pygmies of remote areas of Cameroon and the Central African Republic (CAR), considered to be the oldest inhabitants of these regions. Eight long-term activated T-cell lines producing HTLV-1 gag and env antigens were established from peripheral blood mononuclear cell cultures of HTLV-1 seropositive individuals, including three from Pygmies. A fragment of the env gene encompassing most of the gp21 transmembrane region was sequenced for the 58 new strains, while the complete long terminal repeat (LTR) region was sequenced for 9 strains, including 4 from Pygmies. Comparative sequence analyses and phylogenetic studies performed on both the env and LTR regions by the neighbor-joining and DNA parsimony methods demonstrated that all 22 strains from West and South Africa belong to the widespread cosmopolitan subtype (also called HTLV-1 subtype A). Within or alongside the previously described Zairian cluster (HTLV-1 subtype B), we discovered a number of new HTLV-1 variants forming different subgroups corresponding mainly to the geographical origins of the infected persons, Cameroon, Gabon, and Zaire. Six of the eight Pygmy strains clustered together within this Central African subtype, suggesting a common origin. Furthermore, three new strains (two originating from Pygmies from Cameroon and the CAR, respectively, and one from a Gabonese individual) were particularly divergent and formed a distinct new phylogenetic cluster, characterized by specific mutations and occupying in most analyses a unique phylogenetic position between the large Central African genotype (HTLV-1 subtype B) and the Melanesian subtype (HTLV-1 subtype C). We have tentatively named this new HTLV-1 genotype HTLV-1 subtype D. While the HTLV-1 subtype D strains were not closely related to any known African strain of simian T-cell leukemia virus type 1 (STLV-1), other Pygmy strains and some of the new Cameroonian and Gabonese HTLV-1 strains were very similar (>98% nucleotide identity) to chimpanzee STLV-1 strains, reinforcing the hypothesis of interspecies transmission between humans and monkeys in Central Africa.