Objectives: This research describes the genetic diversity of HIV-1 in Uzbekistan.
Methods: During 2002 and 2003, blood from HIV-positive patients in Uzbekistan was collected, and part of the proviral pol gene and nearly full-length genomes were sequenced and analyzed.
Results: Among 142 Uzbek strains, most clustered genetically with the subtype A strain common in the former Soviet Union. Most of these subtype A-infected drug-naive subjects (65.6%) had an accessory drug resistance mutation, A62V, in the reverse transcriptase gene. Thirteen of the strains (9.2%) clustered with CRF02_AG, an HIV strain common in West Africa. People infected with CRF02_AG were all residents of Tashkent and sampled in 2002. The CRF02_AG strains were monophyletic and probably descended from a single ancestor. Two strains were recombinant between CRF02_AG and subtype A, with each having a different subtype structure. The CRF02_AG and the subtype A elements of the recombinants were monophyletic with Uzbek CRF02_AG and subtype A. New full-length genomes of 12 Uzbek strains suggested that neither the subtype A and nor the CRF02_AG strains in this epidemic were mosaics with other subtypes or circulating recombinant forms.
Conclusion: A genetic analysis of Uzbek HIV strains demonstrated the predominance of subtype A in the epidemic. An outbreak of a West African strain of HIV-1, CRF02_AG, occurred in Tashkent, Uzbekistan in 2002, however. The cocirculation of the 2 strains has resulted in new recombinants that are apparently unique to Uzbekistan.
From the *US Military HIV Research Program, Henry M. Jackson Foundation, Rockville, MD; †Navy Medical Research Unit-3, Cairo, Egypt; ‡National AIDS Center, Tashkent, Uzbekistan; §US Military HIV Research Program, Walter Reed Army Institute of Research, Rockville, MD; and ∥Bloomberg School of Public Health, Johns Hopkins University, Baltimore, MD.
Received for publication May 19, 2004; accepted December 14, 2004
Supported by funds provided by the US Naval Medical Research and Development Command, Silver Spring, MD, Work Unit No. 62787A 873 H B0002, and by a cooperative agreement between the Henry M. Jackson Foundation and the US Department of Defense.
The views and opinions expressed herein do not necessarily reflect those of the US Army, the US Navy, or the Department of Defense.
Reprints: Jean K. Carr, Henry M. Jackson Foundation, 1 Taft Court, Suite 250, Rockville, MD 20850 (e-mail: firstname.lastname@example.org).
The HIV epidemic in former Soviet Union (FSU) countries is growing faster than it is in any other geographic region of the world, although it began there less than a decade ago.1-3 Between 1995 and 2001, the region reported one of the highest annual growth rates in HIV-infected cases; current estimates are that 1.2 to 1.8 million infected persons live in this region.2,3 Most transmission occurs among injecting drug users (IDUs) over a broad geographic area from Riga, Latvia in the west to Irkutsk, Russia in the east.4,5 An analysis of the molecular epidemiology demonstrates, however, that this is a single interconnected epidemic.6,7 Although HIV was initially concentrated in IDU networks, it is now spreading into heterosexual populations in several locations, such as Ukraine and the Russian Federation.7-9
The genetic diversity of HIV-1 is extensive because of error-prone reverse transcription and recombination.10 Worldwide, HIV-1 can be classified into 9 genetic subtypes and 2 sub-subtypes.11 Recombination between subtypes has generated many circulating recombinant forms (CRFs), with recombinants of defined structure spreading in populations; 17 different CRFs have already been identified.11 The 2 most important in the pandemic are the Southeast Asian CRF01_AE and the West African CRF02_AG.12,13
The HIV-1 genetic subtype infecting people in the FSU is, somewhat surprisingly, subtype A, otherwise common only in Africa.7,14,15 This distinguishes the epidemic from that in Europe, where subtype B is most common. Additionally, the strain of subtype A in the FSU is monophyletic, suggesting a common ancestor for all strains; this molecular form of subtype A is called FSU-A.16 Furthermore, the interstrain genetic diversity in the epidemic is unusually low even though the FSU epidemic has been growing for nearly a decade and covers an enormous geographic area.4,6
The epidemic in Uzbekistan, still nascent, is concentrated among IDU networks, where there was an exponential rise in prevalence between 2001 and 2002.7 As such, it is a prototypic “concentrated” epidemic. To date, the genetic subtypes in Uzbekistan have been predominantly FSU-A, with a minority of subtypes B and D.7 The current report documents an outbreak of CRF02_AG, a West African virus, in Tashkent, Uzbekistan in 2002.
After receiving informed consent, blood was drawn from HIV-positive patients at the National AIDS Center in Tashkent, Uzbekistan in 2002 and 2003. Blood was collected in CPT™ vacutainer tubes (Becton Dickinson), and the peripheral blood mononuclear cells (PBMCs) were purified, washed, and then frozen in 95% ethanol at −20°C. Cells were sent to the Navy Medical Research Unit No. 3 (NAMRU-3) in Cairo, Egypt for DNA extraction using QIAamp (QIAgen, Valencia, CA). Under a different study protocol, samples were also collected from residents of rural villages in Cameroon and processed using the same procedures as were used for Uzbek samples.
Partial pol polymerase chain reaction (PCR) amplification from PBMC DNA was performed on a convenience sample of 142 using a nested strategy. The first-round amplification was done with 2 primers: Pro5F (5′-AGAAATTGCAGGGCCCCTAGGAA) and RT3474R (5′-GAATCTCTCTGTTTTCTGCCAG), using Amplitaq Gold (Applied Biosystems, Foster City, CA) and 2 mM of MgCl++ in a total volume of 50 μL. The second-round amplification was completed using the following 2 primers: Pro3F (5′- AGAICAGAGCCAACAGCCCCACCA) and ProRT (5′- TTTCCCCACTAACTTCTGTATGTCATTGACA). Both amplification reactions began with a “hot start” to activate the polymerase, followed by 45 cycles for 30 seconds at 94°C, 30 seconds at 55°C, and 1.5 minutes at 55°C. The amplicon contained the coding sequences for the protease protein and part of the reverse transcriptase (RT) protein, corresponding to nts 2253 through 3266 on HXB-2 (Genbank Accession No. K03455). Nearly full-length genome sequences were also amplified, and all were sequenced in the Applied Biosystems 3100 automated sequencer using Big Dye terminators (Applied Biosystems).17,18
A multiple alignment of the newly derived protease/RT sequences and full-length genome sequences with selected reference sequences was constructed, consisting of 1069 nts and 9914 nts, respectively. Phylogenetic trees were generated, and the consistency of branching order was evaluated using SEQBOOT, DNADIST, NEIGHBOR, and CONSENSE modules of the Phylogeny Inference Package (V3.52c).19 Recombinant analysis used Simplot, version 3.4,20 and alignment examination determined precise breakpoints.13,21 After breakpoint identification, each segment was extracted and analyzed phylogenetically to confirm the assignment of subtype; breakpoint locations were designated relative to HXB-2 (Genbank Accession No. K03455). The Kimura 2-parameter method was used to calculate pairwise genetic distances after removal of hypermutated strains. Analysis of the protease-RT sequences for mutations that might lead to resistance to antiretroviral drugs was also performed.22
Sequences from Uzbekistan were submitted to Genbank under accession numbers AY845716 through AY845857 (partial pol) and AY829203 through AY829214 (nearly full length) and those from Cameroon under AY847362 through AY847453 (partial pol).
There were 142 partial pol sequences of HIV-1 from Uzbekistan and 12 nearly full-length genome sequences. Most of the samples were collected in 2002 (72.5%) and were from male subjects (90.6%), IDUs (84.8%), and residents of the capital, Tashkent (64.6%) (Table 1). Phylogenetic analysis of the sequences revealed that the strains were mostly subtype A (88.0%) and clustered with 97BL006 from Belarus (Fig. 1).23 This phylogenetic cluster of subtype A (FSU-A) has only been found in the FSU and not yet in any part of Africa (see Fig. 1).18 In addition to the subtype A strains, a significant minority were CRF02_AG (9.2%), a CRF common in West and West Central Africa and unexpected in Uzbekistan.7,12
There were also 2 subtype B strains and 2 strains (02UZ813, 02UZ678) that clustered neither with other strains nor with each other. These strains proved to be recombinants between CRF02_AG and FSU-A. 02UZ813 was predominantly CRF02_AG, whereas UZ678 was predominantly subtype A (Fig. 2). The subtype analysis showed that the subtype A portions of the recombinants were monophyletic with FSU-A, supporting the hypothesis that the recombinants were formed locally.
Isolated secondary antiretroviral resistance mutations were found in some patient sequences but none that would confer resistance. One, A62V, is a secondary RT mutation that does not confer drug resistance by itself but is common among the subtype A-infected subjects; 82 (65.6%) had this mutation.24 These strains are also monophyletic within subtype A in the phylogenetic analysis and thus are descended from a single ancestor (see Fig. 1).
The distribution of samples by subtype as well as demographic and risk factor information is shown in Table 1. The only male subject who reported having had sex with another man (MSM) was infected with subtype B. Neither sex nor transmission risk factors distinguished those having subtype A from those with CRF02_AG, but both populations were primarily male (90.4%) and IDUs (85.3%) (see Table 1). Among the 15 patients with CRF02_AG or its recombinants, only 1 reported travel (to Russia) and none reported contact with foreigners. CRF02_AG was significantly concentrated in samples collected in 2002 (P < 0.05), however, and from patients reporting Tashkent as their residence (P < 0.025). The data suggest that the CRF02_AG outbreak began in 2002 among IDUs in Tashkent.
A comparison between CRF02_AG strains from this outbreak and CRF02_AG strains from rural Cameroon, with predominantly heterosexual transmission, showed significantly less variation among the Uzbek strains (0.8%) than among the Cameroon strains (5.1%) (P < 0.001; Fig. 3) and no statistical difference from the diversity that exists among the subtype A strains (1.4%). The intersubject genetic distances between CRF02_AG strains from Uzbekistan, where injecting drug use is the predominant risk factor, were significantly lower than those between strains from sub-Saharan Africa.
A phylogenetic analysis of 12 selected full-length genomes from Uzbekistan consisted of 8 that were subtype A, 3 that were CRF02_AG, and 1 that was hypermutated subtype B (Fig. 4). Careful recombinant analysis of 11 of the strains demonstrated that all of the subtype A strains were entirely subtype A and that all the CRF02_AG strains were entirely CRF02_AG. Up to this time, 3 full-length examples of FSU-A existed, 1 each from Belarus, Russia, and Ukraine; now, there are 11. In addition, the 3 full-length sequences of the CRF02_AG from Uzbekistan should make the epidemic easier to track.
The Uzbek HIV epidemic is dominated by the same subtype A found in other countries of the FSU-FSU-A. A similar finding based on partial env sequences has been previously reported, strongly suggesting that the most common genetic form in Uzbekistan is FSU-A.7 A significant proportion of the subtype A strains in this analysis (65.6%) carried an accessory resistance mutation which, in conjunction with other mutations, would confer resistance to nucleoside RT inhibitors.24 The strains with this mutation form a monophyletic cluster, evidence of a point source introduction.
In this study of HIV-infected patients, a significant minority (9.2%) were infected with CRF02_AG, a recombinant common in West and West Central Africa.12 CRF02_AG has also spread to a small number of people in Estonia, Korea, Taiwan, Ecuador, and Israel, but most of these epidemics have been associated with self-contained transmission clusters.25,26 The outbreak of CRF02_AG in Uzbekistan may have been a single contained transmission cluster at the time of this report, but recent evidence suggests further spread (M. D. Saah, personal communication, 2004). It is possible that the strain is, in essence, a molecular marker for 1 particular drug-using network. The interconnectedness of drug-using networks in different countries of the FSU is not well understood, and the spread of this particular CRF02_AG strain might reveal some of those connections.
The presence of 2 unique recombinants of CRF02_AG and FSU-A suggests that 2 drug networks, 1 with CRF02_AG and 1 with FSU-A, may have overlapped to some extent. Unique recombinants are usually found in people who have been dually infected with both genetic forms; sampling from PBMC DNA in 1 such person over time has revealed either or both parental forms and a haphazard distribution of recombinants.27 The close genetic relationship that each of the Uzbek recombinants has to existing Uzbek strains argues for a local origin of the recombinants in Tashkent, Uzbekistan.
Low intersubject genetic diversity characterizes HIV-1 from IDU epidemics in many countries and subtypes.28 The low genetic diversity of the A epidemic in the FSU is, by far, the largest and the most enduring of these epidemics.6 It has been spreading since 1996 and now covers an immense territory. The Tashkent CRF02_AG outbreak in IDUs also has low diversity. HIV-1, like other retroviruses, is notoriously mutagenic, making approximately 1 error per replication cycle,10 with an estimated 1010 progeny virions produced each day in an untreated HIV patient29; the mechanism by which viral populations in IDUs have remained relatively unchanged is difficult to understand and deserves more thorough and systematic investigation.
The outbreak of CRF02_AG in Tashkent demonstrates the rapid spread of new viruses into drug-using networks of the FSU. When a new strain enters a drug network, it begins to recombine with extant strains and to create new recombinants unique to that drug network.30 Unfortunately, drug networks provide the ideal means for rapidly spreading HIV within core groups and represent an emerging threat for the introduction of novel or recombinant strains into larger heterosexually active communities.
The authors thank Robert Shafer of Stanford University for his helpful comments.
1. Hamers FF. HIV infection in Ukraine (1987-1996). Epidemiol Sante Publique
2. Joint United Nations Program on HIV/AIDS (UNAIDS). AIDS Epidemic Update
. UNAIDS; Geneva, Switzerland, 2003.
3. United Nations Development Programme, New York, NY, USA. HIV/AIDS in Eastern Europe and the Commonwealth of Independent States: Reversing the Epidemic, Facts and Policy Options
4. Ferdats A, Konicheva V, Dievberna I, et al. An HIV type 1 subtype A outbreak among injecting drug users in Latvia. AIDS Res Hum Retroviruses
5. Lukashov V, Huismans R, Rakhmanova A, et al. Circulation of subtype A and gagA/envB recombinant HIV type 1 strains among injecting drug users in St. Petersburg, Russia, correlates with geographical origin of infections. AIDS Res Hum Retroviruses
6. Bobkov A, Kazennova E, Khanina T, et al. An HIV type 1 subtype A strain of low genetic diversity continues to spread among injecting drug users in Russia: study of the new local outbreaks in Moscow and Irkutsk. AIDS Res Hum Retroviruses
7. Kurbanov F, Kondo M, Tanaka Y, et al. Human immunodeficiency virus in Uzbekistan: epidemiological and genetic analyses. AIDS Res Hum Retroviruses
8. Lowndes CM, Renton A, Alary M, et al. Conditions for widespread heterosexual spread of HIV in the Russian Federation: implications for research, monitoring and prevention. Int J Drug Policy
9. Mavrov GI, Bondarenko GM. The evolution of sexually transmitted infections in the Ukraine. Sex Transm Infect
10. Coffin JM. HIV population dynamics in vivo: implications for genetic variation, pathogenesis, and therapy. Science
11. Kuiken CL, Foley B, Freed E, et al. Sequence Compendium
, vol. IV. LA-UR 03-3564, Los Alamos, NM: Theoretical Biology and Biophysics Group, Los Alamos National Laboratory; 2002.
12. Carr JK, Salminen MO, Albert J, et al. Full genome sequences of human immunodeficiency virus type 1 subtypes G and A/G intersubtype recombinants. Virology
13. Carr JK, Salminen MO, Koch C, et al. Full-length sequence and mosaic structure of a human immunodeficiency virus type 1 isolate from Thailand. J Virol
14. Novitsky VA, Montano MA, Essex M. Molecular epidemiology of an HIV-1 subtype A subcluster among injection drug users in the southern Ukraine. AIDS Res Hum Retroviruses
15. Lukashov VV, Karamov EV, Eremin VF, et al. Extreme founder effect in an HIV type 1 subtype A epidemic among drug users in Svetlogorsk, Belarus. AIDS Res Hum Retroviruses
16. Masharsky A, Klimov N, Kozlov AP. Molecular cloning and analysis of full-length genome of HIV type 1 strains prevalent in countries of the former Soviet Union. AIDS Res Hum Retroviruses
17. Salminen M, Koch C, Sanders-Buell E, et al. Recovery of virtually full-length HIV-1 provirus of diverse subtypes from primary virus cultures using the polymerase chain reaction. Virology
18. Carr JK, Laukkanen T, Salminen MO, et al. Characterization of subtype A HIV-1 from Africa by full genome sequencing. AIDS
19. Felsenstein J. Phylip (Phylogeny Inference Package)
[computer program], V3.5c. Seattle: University of Washington Department of Genetics; 1992.
20. Ray SC, Lubaki N, Dhruva BR, et al. Autologous strain-specific cytolytic T lymphocyte responses directed against human immunodeficiency virus type 1 env. AIDS Res Hum Retroviruses
21. Salminen MO, Carr JK, Burke DS, et al. Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning. AIDS Res Hum Retroviruses
22. Shafer RW, Jung DR, Betts BJ, et al. Human immunodeficiency virus reverse transcriptase and protease sequence database. Nucleic Acids Res
23. Liitsola K, Holm K, Bobkov A, et al. An AB recombinant and its parental HIV type 1 strains in the area of the former Soviet Union: low requirements for sequence identity in recombination. AIDS Res Hum Retroviruses
24. Shafer RW, Jung DR, Betts BJ. Human immunodeficiency virus type 1 reverse transcriptase and protease mutation search engine for queries. Nat Med
25. Ustina V, Zilmer K, Tammai L, et al. Epidemiology of HIV in Estonia. AIDS Res Hum Retroviruses
26. Carrion G, Hierholzer J, Montano S, et al. Circulating recombinant form CRF02_AG in South America. AIDS Res Hum Retroviruses
27. Salminen MO, Carr JK, Robertson DL, et al. Evolution and probable transmission of intersubtype recombinant human immunodeficiency virus type 1 in a Zambian couple. J Virol
28. Oelrichs RB, Shrestha IL, Anderson DA, et al. The explosive human immunodeficiency virus type 1 epidemic among injecting drug users of Kathmandu, Nepal, is caused by a subtype C virus of restricted genetic diversity. J Virol
29. Perelson AS, Neumann AU, Markowitz M, et al. HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time. Science
30. Liitsola K, Tashkinova I, Laukkanen T, et al. HIV-1 genetic subtype A/B recombinant strain causing an explosive epidemic in injecting drug users in Kaliningrad. AIDS
Keywords:© 2005 Lippincott Williams & Wilkins, Inc.
HIV-1 molecular epidemiology; former Soviet Union; injecting drug user epidemic; genetic diversity