To examine intersubtype recombination, pol sequences were analysed by bootscanning, which revealed that 12 were recombinant (Fig. 2d–l), with subtypes (in 5′–3′ order) as follows (U denotes unknown subtype): three UK (CU33, CU34 and CU55, which grouped in the phylogenetic tree with the African MAL isolate [23–25]), two GBGB (CU100 and CU103), one BU (CU13), one UB (CU66), two GD (CU45 and CU57, which exhibit different crossover points), one DB (CU54), one BD (CU17) and one DUD (CU70). In all cases, grouping of different segments with different reference sequences was supported by significant bootstrap values, and results were consistent independently of reference sequences used. Bootscan analysis and phylogenetic trees of partial segments of the amplified sequences indicated that D subtype and U segments of recombinant pol sequences of CU13, CU66, CU45, CU57, CU54, CU17 and CU70 (but not unclassified segments of CU33, CU34 and CU55) grouped with viruses of the D subtype and U pol clusters of Cuba, respectively, described above. Similarly, the G subtype segments of CU57 GD recombinant sequence clustered with G subtype pol sequences of CU74, CU85 and CU87 (which form a monophyletic group). This indicates that these recombinant were probably generated by recombination of genetic forms circulating in Cuba.
Phylogenetic analysis of V3 sequences
In phylogenetic trees (Fig. 1b,d), sequences of the V3 region grouped with the following subtypes or CRF: A (n = 26), B (n = 54), C (n = 2), D (n = 1), F1 (n = 1), G (n = 6), H (n = 11), J (n = 1) and CRF01-AE (n = 1). In two samples, CU10 and CU72, which are of subtypes C and G, respectively, in pol, the V3 region could not be amplified.
Subtypes in V3 were coincident with those in pol for all viruses that were of non-recombinant subtypes B, C, F1, G, H and J in pol. Viruses of the Cuban D subtype pol cluster, except CU8, were of subtype A in V3, forming in this segment a monophyletic group, which also included viruses with recombinant pol sequences CU17 (BD in pol) and CU57 (GD in pol). Viruses forming the U cluster in pol, except CU46 and CU75, were of subtype H in V3, forming in this segment a monophyletic group, which also included CU8 (D in pol), CU45 (GD in pol) and CU70 (DUD in pol). CU46 grouped with subtype B, and CU75 with CRF01-AE in V3. The three viruses that were MAL-like in pol grouped with subtype A reference sequences in V3, similarly to the isolate X327 of Spain , but different from both MAL [23–25] and the partly MAL-like virus NOGIL3 of Norway , which are of subtypes D and H, respectively, in env. Subtypes in pol and V3 of all viruses are shown in Tables 1 and 2, with distribution of genetic forms in Cuba represented graphically in Fig. 3. Overall, 43 (78%) of the 55 non-B subtype viruses were recombinant. The number of genetic forms identified was 21, including seven non-recombinant subtypes and 14 recombinant forms.
Phylogenetic relations with database sequences
Homologies of non-B subtype sequences with sequences deposited in GenBank were searched using the BLAST algorithm. No sequences homologous to those of the U pol cluster grouping with these in phylogenetic trees were found in the database. There were 10 database H subtype sequences that grouped with the Cuban U/H viruses in V3, with significant bootstrap values (80%), seven of which were from the DRC (Fig. 4). One subtype A sequence from the Republic of Congo grouped with the Cuban D/A recombinants in V3. Pol sequences of these African viruses are not available; therefore, it is not known if they are phylogenetically related to the U/H and D/A recombinants of Cuba, respectively, in pol.
Pol sequences of CU33, CU34 and CU55 clustered with MAL isolate and with MAL-like pol sequences of NOGIL3 of Norway  and X327 of Spain . In V3, CU33 and CU34 grouped with X327 and with two other subtype A sequences of the Republic of Congo (Fig. 4). Pol sequences of these two African viruses are not available and, therefore, it is not known if they are also MAL-like in this segment.
Viruses of subtypes A and H in env, previously reported in Cuba , grouped with D/A and U/H recombinant viruses, respectively, of our study (Fig. 4).
There were differences in the prevalence of non-B subtype viruses between groups defined by gender, risk category and place of infection. Non-B subtype viruses (including recombinants) were more frequent among heterosexuals (72%) than among homo/bisexuals (37%) (P < 0.001), and among women (73%) than among men (44%) (P < 0.05). However, the prevalence of non-B subtype viruses among heterosexual men (69%) did not differ significantly from that in women.
Nine of ten infections (90%) reported to having been acquired in Africa were with non-B viruses. Notably, seven were non-recombinant. Subtypes of non-B viruses acquired in Africa were: three C (two acquired in Ethiopia and one in Zambia), one F1 (Angola), one G (DRC), one H (Angola), one J (Angola) and two MAL-like in pol and subtype A in env (both acquired in DRC). Among infections acquired in Cuba, 46 (50%) of 92 were with viruses of non-B subtypes.
Within subtype B and D/A recombinant clusters, there were some subclusters, supported by high bootstrap values, of viruses found in individuals sharing epidemiological features. One of these subclusters comprised 14 subtype B viruses (signalled with brackets in Fig. 1), all of which corresponded to homo- or bisexual men, 12 of them from Havana City, representing 46% of subtype B viruses among individuals of this risk category in this city. Within the D/A recombinant cluster, five of the six viruses found in homo- or bisexual men from Havana City formed a subcluster supported by high bootstrap values (signalled with brackets in Fig. 1). Genetic distances in phylogenetic trees (Fig. 1) and dates of diagnosis (1996 or later in all but one of the B subtype subcluster, and 1997 or 1998 in the D/A recombinant subcluster) were consistent with relatively recent outbreaks, each originating from a common source.
Here we report the most extensive survey to date on HIV genetic diversity in Cuba, with 105 samples analysed from an estimated HIV-infected population of approximately 2000 individuals. Previous studies were limited to 15 individuals in 1995  and a recently published study of 11 individuals . In both studies, only the V3 region was analysed. In the first study, it was stated that the predicted V3 amino acid sequences were similar to those of subtype B reference isolates. In the second study, V3 sequences were of subtypes B (n = 6), A (n = 2), C (n = 2) and H (n = 1). In our study, two segments of the viral genome, pol and env, were analysed phylogenetically using neighbour-joining trees and also by bootscanning to detect possible intersubtype recombination. The results show that approximately half (52%) of HIV-1 infections are caused by non-B subtype or recombinant viruses, with a high diversity of genetic forms, including non-recombinant viruses of seven subtypes (B, C, D, F1, G, H and J) and recombinant viruses of 14 different genetic forms in 78% of the non-B samples, containing segments of subtypes A, B, D, G, H and K and CRF01-AE, as well as unclassified segments. Only one of the recombinant forms, a virus related to a MAL-like isolate reported in Spain (but probably acquired in Africa) , has been reported previously. Nevertheless, in spite of this diversity, our results provide phylogenetic and epidemiological evidence of only four genetic forms (B and G subtypes and Dpol/Aenv and Upol/Henv recombinants) currently circulating in Cuba. The remaining viruses either were acquired in Africa or were detected in only one or two individuals. D/A recombinants (n = 21) and U/H viruses (n = 7) grouped in respective monophyletic clusters in both pol and env, suggesting that they are candidates for recognition as circulating recombinant forms, pending sequencing of full-length genomes . These two recombinant forms might be of African origin, since the parental viruses were not detected in Cuba. In GenBank, no sequences phylogenetically related to the pol sequences of the U/H or D/A recombinants of Cuba are found, but there are African viruses with env sequences related to these recombinants (Fig. 4). pol sequences of these African viruses are not available, which does not allow us to confirm if recombinant viruses related to those circulating in Cuba are present in Africa.
The HIV-1 genetic diversity found in Cuba has no parallel in the Americas. In other countries of the Caribbean area, infections are almost uniformly with subtype B, although occasional cases of infections with other subtypes have been reported [30–32]. In the American continent, non-B subtype infections are relatively common only in Brazil, Argentina and Uruguay, although these are limited to F and C subtypes and BF recombinants in Brazil [1,2,33,34], BF recombinants in Argentina  and F subtype (possibly BF recombinants) and imported CRF01-AE viruses in Uruguay [4,35].
A diversity of HIV-1 genetic forms of a degree comparable to that found in Cuba has only been reported in Central African countries, with the highest group M diversity found in the DRC [8,36]. However, considering that the estimated number of HIV infections in Cuba is only approximately 2000, the HIV genetic diversity found in Cuba, in relative terms, has no parallel in any other country.
Several factors may have contributed to the high HIV genetic diversity in Cuba. Large numbers of Cuban soldiers and civilian personnel were stationed in the 1970s and 1980s in Angola [6,14]. Although there are no published studies on HIV-1 genetic diversity in Angola, it would not be unexpected to find a high HIV-1 genetic diversity in this country, since it is bordering the DRC (in fact, the three infections that were acquired in Angola were of three different subtypes, F1, H and J). Additionally, large numbers of ‘internationalist’ Cuban aid workers, serving in several sub-Saharan African countries, may have also contributed to the introduction of non-B genetic forms [6,10,14]. In our study, only 10 of 105 infections were acquired in Africa, indicating that non-B viruses introduced from Africa have started to circulate in Cuba. A superimposed factor contributing to the generation of HIV genetic diversity in Cuba is recombination . Phylogenetic analyses indicate that 12 infections are with 11 different viruses generated by recombination between the four genetic forms circulating in Cuba. Considering the low prevalence of HIV infections in Cuba, the diversity of recombinants generated locally is disproportionately high. Although there is no direct evidence of this, it could be suspected that recombination of HIV genetic forms would have been facilitated by the former policy of prolonged quarantine of all HIV-infected individuals in selected sanatoria [10,14]. This, together with the scarcity of condoms in Cuba in earlier years , could have set conditions favourable for coinfections with multiple genetic forms, with subsequent generation of recombinants.
The finding of a high HIV diversity in Cuba may have implications for vaccine design [39–42] and for the use of tests for viral load determination [43,44] and for detection of drug resistance mutations [45–47]. Although the prevalence of HIV infections in Cuba is the lowest in the Americas, and travel of Cubans outside the country is presently limited, the possibility of exporting some of the multiple genetic forms of Cuba to other countries may not be negligible, considering the expansion of tourism in Cuba in recent years, with a concomitant increase in casual sexual contacts, in relation with declining economic conditions [14,38,48,49].
A recently reported HIV-1 complex recombinant virus from Cameroon (CM53379) sequenced in the full-length genome  clusters with high boostrap values, both in pol and env, with the U (pol)/H (env) recombinants from Cuba. This supports the Central African ancestry of these Cuban viruses.
In summary, in contrast to the almost uniform predominance of subtype B in other Caribbean countries, a high diversity of HIV-1 genetic forms, unprecedented in the Americas, has been found in Cuba. Some of these genetic forms were imported from Africa and others were generated locally by recombination. The causes of this diversity are related to historical circumstances peculiar to Cuba, which have contributed to make HIV diversity in this country a mosaic of the West (represented by B subtype sequences, presumably of North American origin) and of Africa (represented by multiple non-B subtype and recombinant viruses) grafted into the Caribbean.
We thank José Esparza for his contribution to organizing the UNAIDS program under which this study was done, and Francisco Parras for his support of this study.
1. Louwagie J, Delwart EL, Mullins JI, McCutchan FE, Eddy G, Burke DS. Genetic analysis of HIV-1
isolates from Brazil reveals presence of two distinct genetic subtypes
. AIDS Res Hum Retroviruses 1994, 10: 561–567.
2. Morgado MG, Sabino EC, Shpaer EG. et al
. V3 region polymorphisms in HIV-1
from Brazil: prevalence of subtype B strains divergent from North American/European prototype and detection of subtype F. AIDS Res Hum Retroviruses 1994, 10: 569–576.
3. Thomson MM, Villahermosa ML, Vázquez-de-Parga E. 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.
4. Russell KL, Carcamo C, Watts DM. et al
. Emerging genetic diversity of HIV-1
in South America. AIDS 2000, 14: 1785–1791.
5. Gómez CE, Iglesias E, Perdomo W. et al
. Isolates from four different HIV type 1 clades circulating in Cuba
identified by DNA sequence of the C2–V3 region. AIDS Res Hum Retroviruses 2001, 17: 55–58.
6. Torres-Anjel MJ. Macroepidemiology of the HIVs-AIDS (HAIDS) pandemic. Insufficiently considered zoological and geopolitical aspects.
Ann N Y Acad Sci 1992, 653: 257–273.
7. Vidal N, Peeters M, Mulanga-Kabeya C. et al
. Unprecedented degree of human immunodeficiency virus type 1 (HIV-1
) group M genetic diversity in the Democratic Republic of Congo suggests that the HIV-1
pandemic originated in Central Africa. J Virol 2000, 74: 10498–10507.
8. Scheper-Hughes N. AIDS, public health, and human rights in Cuba
. Lancet 1993, 342: 965–967.
9. Ubell RN. High-tech medicine in the Caribbean. 25 years of Cuban health care.
N Engl J Med 1983, 309: 1468–1472.
10. Pérez-Stable EJ. Cuba
's response to the HIV epidemic. Am J Public Health 1991, 81: 563–567.
11. UNAIDS. Epidemiological fact sheet on HIV/AIDS and sexually transmitted infections. 2000 update. http://www.unaids.org
12. Voelker R. HIV/AIDS in the Caribbean: big problems among small islands. JAMA 2001, 285: 2961–2963.
13. Hsieh YH, Chen CW, Lee SM, de Arazoza H. On the recent sharp increase in HIV detections in Cuba
. AIDS 2001, 15: 426–428.
14. Santana S, Faas L, Wald K. Human immunodeficiency virus in Cuba
: the public health response of a Third World country. Int J Health Serv 1991, 21: 511–537.
15. Santana S. Article on HIV in Cuba
criticized. Am J Public Health 1992, 82: 899–900.
16. Ruibal-Brunet IJ, Cuevas MT, Díaz-Torres H. et al
. Genotypic resistance mutations to antiretroviral drugs in HIV-1
B and non-B subtypes
. Pan Am J Public Health 2001, 10: 174–179.
17. Tenorio A, Echevarría JE, Casas I, Echevarría JM, Tabarés E. Detection and typing of human herpesviruses by multiplex polymerase chain reaction. J Virol Meth 1993, 44: 261–269.
18. Villahermosa ML, Thomson M, Vázquez-de Parga E. et al
. Improved conditions for extraction and amplification of human immunodeficiency virus type 1 RNA from plasma samples with low viral load. J Hum Virol 2000, 3: 27–34.
19. Leitner T, Korovina G, Marquina S, Smolskaya T, Albert J. Molecular epidemiology
and MT-2 cell tropism of Russian HIV type 1 variant. AIDS Res Hum Retroviruses 1996, 12: 1595–1603.
20. Werle E, Schneider C, Renner M, Volker M, Fiehn W. Convenient single-step, one tube purification of PCR products for direct sequencing. Nucl Acids Res 1994, 22: 4354–4355.
21. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl Acids Res 1997, 25: 4876–4882.
22. Hillis DM, Bull JJ. An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis. Syst Biol 1993, 42: 182–192.
23. Alizon M, Wain-Hobson S, Montagnier L, Sonigo P. Genetic variability of the AIDS virus: nucleotide sequence analysis of two isolates from African patients. Cell 1986, 46: 63–74.
24. Robertson DL, Hahn BH, Sharp PM. Recombination in AIDS viruses. J Mol Evol 1995, 40: 249–259.
25. Gao F, Robertson DL, Carruthers CD. et al
. An isolate of human immunodeficiency virus type 1 originally classified as subtype I represents a complex mosaic comprising three different group M subtypes
(A, G, and I). J Virol 1998, 72: 10234–10241.
26. Thomson MM, Delgado E, Manjón N. et al
genetic diversity in Galicia, Spain: BG intersubtype recombinant viruses are circulating among injecting drug users. AIDS 2001, 15: 509–516.
27. Jonassen TO, Grinde B, Asjo B, Hasle G, Hungnes O. Intersubtype recombinant HIV type 1 involving HIV-MAL-like and subtype H-like sequence in four Norwegian cases. AIDS Res Hum Retroviruses 2000, 16: 49–58.
28. Rolo FM, Miranda L, Wainberg MA. et al
. Envelope V3 region sequences of Cuban HIV-1
isolates. J Acquir Immune Defic Syndr Hum Retrovirol 1995, 9: 123–125.
29. Robertson DL, Anderson JP, Bradac JA. et al
nomenclature proposal. Science 2000, 288: 55–56.
30. Desgranges C, Fillon S, Audoly G. et al
. Presence of HIV-1 subtypes
B and F and HTLV-I in HIV/HTLV coinfected individuals of Martinique. J Acquir Immune Defic Syndr Hum Retrovirol 1996, 13: 468–470.
31. Ouka M, Menu E, Barin F. et al
. Introduction of HIV-1
subtype A in the French West Indies. J Acquir Immune Defic Syndr Hum Retrovirol 1998, 19: 308–310.
32. Flores I, Pieniazek D, Moran N. et al
subtype F in single and dual infections in Puerto Rico: a potential sentinel site for monitoring novel genetic HIV variants in North America. Emerg Infect Dis 1999, 5: 481–483.
33. Csillag C. HIV-1
subtype C in Brazil. Lancet 1994, 344: 1354.1354.
34. Ramos A, Tanuri A, Schechter M. et al
. Dual and recombinant infections: an integral part of the HIV-1
epidemic in Brazil. Emerg Infect Dis 1999, 5: 65–74.
35. Artenstein AW, Coppola J, Brown AE. et al
. Multiple introductions of HIV-1
subtype E into the western hemisphere. Lancet 1995, 346: 1197–1198.
36. Rambaut A, Robertson DL, Pybus OG, Peeters M, Holmes EC. Human immunodeficiency virus. Phylogeny and the origin of HIV-1.
Nature 2001, 410: 1047–1048.
37. Robertson DL, Sharp PM, McCutchan FE, Hahn BH. Recombination in HIV-1
. Nature 1995, 374: 124–126.
38. Veeken H. Cuba
: plenty of care, few condoms, no corruption. BMJ 1995, 311: 935–937.
39. Mascola JR, Louwagie J, McCutchan FE. et al
. Two antigenically distinct subtypes
of human immunodeficiency virus type 1: viral genotype predicts neutralization serotype. J Infect Dis 1994, 169: 48–54.
40. Rowland-Jones SL, Dong T, Fowke KR. et al
. Cytotoxic T cell responses to multiple conserved HIV epitopes in HIV- resistant prostitutes in Nairobi. J Clin Invest 1998, 102: 1758–1765.
41. Dorrell L, Dong T, Ogg GS. et al
. Distinct recognition of non-clade B human immunodeficiency virus type 1 epitopes by cytotoxic T lymphocytes generated from donors infected in Africa. J Virol 1999, 73: 1708–1714.
42. Cao H, Mani I, Vincent R. et al
. Cellular immunity to human immunodeficiency virus type 1 (HIV-1
) clades: relevance to HIV- 1 vaccine trials in Uganda. J Infect Dis 2000, 182: 1350–1356.
43. Burgisser P, Vernazza P, Flepp M. et al
. Performance of five different assays for the quantification of viral load in persons infected with various subtypes of HIV-1. Swiss HIV Cohort Study.
J Acquir Immune Defic Syndr 2000, 23: 138–144.
44. Jagodzinski LL, Wiggins DL, McManis JL. et al
. Use of calibrated viral load standards for group M subtypes
of human immunodeficiency virus type 1 to assess the performance of viral RNA quantitation tests. J Clin Microbiol 2000, 38: 1247–1249.
45. Vahey M, Nau ME, Barrick S. et al
. Performance of the Affymetrix GeneChip HIV PRT 440 platform for antiretroviral drug resistance genotyping of human immunodeficiency virus type 1 clades and viral isolates with length polymorphisms. J Clin Microbiol 1999, 37: 2533–2537.
46. Lee CC, Thoe SY, Leo YS, Chan KP, Ling AE. Detection of drug-selected mutations in HIV-1
subtype E reverse transcriptase sequence using the line probe assay. AIDS 2000, 14: 1064–1065.
47. Debyser Z, van Wijngaerden E, van Laethem K. et al
. Failure to quantify viral load with two of the three commercial methods in a pregnant woman harboring an HIV type 1 subtype G strain. AIDS Res Hum Retroviruses 1998, 14: 453–459.
48. Burr C. Havana. Assessing Cuba's approach to contain AIDS and HIV.
Lancet 1997, 350: 647.647.
49. Barry M. Effect of the U.S. embargo and economic decline on health in Cuba
. Ann Intern Med 2000, 132: 151–154.
50. Carr JK, Torimiro JN, Wolfe ND. et al
. The AG recombinant 1bNG and novel strains of group M HIV-1
are common in Cameroon. Virology 2001, 268: 168–181.
Keywords:© 2002 Lippincott Williams & Wilkins, Inc.
HIV-1; Cuba; subtypes; recombinant forms; molecular epidemiology