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

Reconstructing the Temporal Origin and the Transmission Dynamics of the HIV Subtype B Epidemic in St. Petersburg, Russia

The HIV/AIDS epidemic in Russia is among the fastest growing in the world. HIV epidemic burden is non-uniform in different Russian regions and diverse key populations. An explosive epidemic has been documented among people who inject drugs (PWID) starting from the mid-1990s, whereas presently, the majority of new infections are linked to sexual transmission. Nationwide, HIV sub-subtype A6 (previously called A_FSU) predominates, with the increasing presence of other subtypes, namely subtype B and CRF063_02A. This study explores HIV subtype B sequences from St. Petersburg, collected from 2006 to 2020, in order to phylogenetically investigate and characterize transmission clusters, focusing on their evolutionary dynamics and potential for further growth, along with a socio-demographic analysis of the available metadata. In total, 54% (107/198) of analyzed subtype B sequences were found grouped in 17 clusters, with four transmission clusters with the number of sequences above 10. Using Bayesian MCMC inference, tMRCA of HIV-1 subtype B was estimated to be around 1986 (95% HPD 1984-1991), whereas the estimated temporal origin for the four large clusters was found to be more recent, between 2001 and 2005. The results of our study imply a complex pattern of the epidemic spread of HIV subtype B in St. Petersburg, Russia, still in the exponential growth phase, and in connection to the men who have sex with men (MSM) transmission, providing a useful insight needed for the design of public health priorities and interventions.

Prevalence and spatiotemporal dynamics of HIV-1 Circulating Recombinant Form 03_AB (CRF03_AB) in the Former Soviet Union countries

Background

HIV-1 circulating recombinant forms (CRFs) infections has been increasing in Former Soviet Union (FSU) countries in the recent decade. One is the CRF03_AB, which circulated in the region since late 1990s and probably became widespread in northwestern FSU countries. However, there is not much information provided about the dissemination of this recombinant. Here, we examine the prevalence, evolutionary dynamics and dispersion pattern of HIV-1 CRF03_AB recombinant.

Methods

We analyzed 32 independent studies and 151 HIV-1 CRF03_AB pol sequences isolated from different FSU countries over a period of 22 years. Pooled prevalence was estimated using a random effects model. Bayesian coalescent-based method was used to estimate the evolutionary, phylogeographic and demographic parameters.

Results

Our meta-analysis showed that the pooled prevalence of CRF03_AB infection in northwestern FSU region was 5.9% [95%CI: 4.1–7.8]. Lithuania (11.6%), Russia (5.9%) and Belarus (2.9%) were the most affected by CRF03_AB. We found that early region wide spread of HIV-1 CRF03_AB originated from one viral clade that arose in the city of Kaliningrad in 1992 [95%HPD: 1990–1995]. Fourteen migration route of this variant were found. The city of Kaliningrad is involved in most of these, confirming its leading role in CRF03_AB spread within FSU. Demographic reconstruction point to this is that CRF03_AB clade seems to have experienced an exponential growth until the mid-2000s and a decrease in recent years.

Conclusion

These data provide new insights into the molecular epidemiology of CRF03_AB as well as contributing to the fundamental understanding of HIV epidemic in FSU.

Local Epidemics Gone Viral: Evolution and Diffusion of the Italian HIV-1 Recombinant Form CRF60_BC

The molecular epidemiology of HIV-1 in Italy is becoming increasingly complex, mainly due to the spread of non-B subtypes and the emergence of new recombinant forms. We previously characterized the outbreak of the first Italian circulating recombinant form (CRF60_BC), occurring among young MSM living in Apulia between the years 2009 and 2011. Here we show a 5-year follow-up surveillance to trace the evolution of CRF60_BC and to investigate its further spread in Italy. We collected additional sequences and clinical data from patients harboring CRF60_BC, enrolled at the Infectious Diseases Clinic of the University of Bari. In addition to the 24 previously identified sequences, we retrieved 27 CRF60_BC sequences from patients residing in Apulia, whose epidemiological and clinical features did not differ from those of the initial outbreak, i.e., the Italian origin, young age at HIV diagnosis (median: 24 years; range: 18–37), MSM risk factor (23/25, 92%) and recent infection (from 2008 to 2017). Sequence analysis revealed a growing overall nucleotide diversity, with few nucleotide changes that were fixed over time. Twenty-seven additional sequences were detected across Italy, spanning multiple distant regions. Using a BLAST search, we also identified a CRF60_BC sequence isolated in United Kingdom in 2013. Three patients harbored a unique second generation recombinant form in which CRF60_BC was one of the parental strains. Our data show that CRF60_BC gained epidemic importance, spreading among young MSM in multiple Italian regions and increasing its population size in few years, as the number of sequences identified so far has triplicated since our first report. The observed further divergence of CRF60_BC is likely due to evolutionary bottlenecks and host adaptation during transmission chains. Of note, we detected three second-generation recombinants, further supporting a widespread circulation of CRF60_BC and the increasing complexity of the HIV-1 epidemic in Italy.

A New HIV-1 Circulating Recombinant Form (CRF98_cpx) Between CRF06_cpx and Subtype B Identified in Southwestern France

During a recent study on the sequencing data of our database between 2012 and 2016 in Southwestern France, we observed that eight patients harbored what seemed to be the same virus. Indeed, routine genotyping at the time of HIV diagnosis showed that protease and reverse transcriptase were related to CRF06_cpx and subtype B, respectively. The integrase sequences (available for three patients) were clustering with CRF06_cpx and envelope (Env) gp120 sequences (available for two patients) with subtype B. Since such a recombinant has not been recorded in the Los Alamos database, we decided to characterize the full-length genome of this virus. The data suggest the identification of a new circulating recombinant form (CRF) between CRF06_cpx and subtype B, the structure of which is very complex with multiple breakpoints. We will refer this CRF as CRF98_cpx.

Overview of HIV molecular epidemiology among people who inject drugs in Europe and Asia

HIV strains continuously evolve, tend to recombine, and new circulating variants are being discovered. Novel strains complicate efforts to develop a vaccine against HIV and may exhibit higher transmission efficiency and virulence, and elevated resistance to antiretroviral agents. The United Nations Joint Programme on HIV/AIDS (UNAIDS) set an ambitious goal to end HIV as a public health threat by 2030 through comprehensive strategies that include epidemiological input as the first step of the process. In this context, molecular epidemiology becomes invaluable as it captures trends in HIV evolution rates that shape epidemiological pictures across several geographical areas. This review briefly summarizes the molecular epidemiology of HIV among people who inject drugs (PWID) in Europe and Asia. Following high transmission rates of subtype G and CRF14_BG among PWID in Portugal and Spain, two European countries, Greece and Romania, experienced recent HIV outbreaks in PWID that consisted of multiple transmission clusters including subtypes B, A, F1, and recombinants CRF14_BG and CRF35_AD. The latter was first identified in Afghanistan. Russia, Ukraine, and other Former Soviet Union (FSU) states are still facing the devastating effects of epidemics in PWID produced by A_FSU (also known as IDU-A), B_FSU (known as IDU-B), and CRF03_AB. In Asia, CRF01_AE and subtype B (Western B and Thai B) travelled from PWID in Thailand to neighboring countries. Recombination hotspots in South China, Northern Myanmar, and Malaysia have been generating several intersubtype and inter-CRF recombinants (e.g. CRF07_BC, CRF08_BC, CRF33_01B etc.), increasing the complexity of HIV molecular patterns.

HIV-1 A1 Subtype Epidemic in Italy Originated from Africa and Eastern Europe and Shows a High Frequency of Transmission Chains Involving Intravenous Drug Users

Background

Subtype A accounts for only 12% of HIV-1 infections worldwide but predominates in Russia and Former Soviet Union countries of Eastern Europe. After an early propagation via heterosexual contacts, this variant spread explosively among intravenous drug users. A distinct A1 variant predominates in Greece and Albania, which penetrated directly from Africa. Clade A1 accounts for 12.5% of non-B subtypes in Italy, being the most frequent after F1 subtype.

Aim

Aim of this study was to investigate the circulation of A1 subtype in Italy and trace its origin and diffusion through phylogenetic and phylodynamic approaches.

Results

The phylogenetic analysis of 113 A1 pol sequences included in the Italian ARCA database, indicated that 71 patients (62.8%) clustered within 5 clades. A higher probability to be detected in clusters was found for patients from Eastern Europe and Italy (88.9% and 60.4%, respectively) compared to those from Africa (20%) (p < .001). Higher proportions of clustering sequences were found in intravenous drug users with respect to heterosexuals (85.7% vs. 59.3%, p = .056) and in women with respect to men (81.4% vs. 53.2%, p < .006). Subtype A1 dated phylogeny indicated an East African origin around 1961. Phylogeographical reconstruction highlighted 3 significant groups. One involved East European and some Italian variants, the second encompassed some Italian and African strains, the latter included the majority of viruses carried by African and Italian subjects and all viral sequences from Albania and Greece.

Conclusions

Subtype A1 originated in Central Africa and spread among East European countries in 1982. It entered Italy through three introduction events: directly from East Africa, from Albania and Greece, and from the area encompassing Moldavia and Ukraine. As in previously documented A1 epidemics of East European countries, HIV-1 A1 subtype spread in Italy in part through intravenous drug users. However, Eastern European women contributed to the penetration of such variant, probably through sex work.

Phylogeographic analysis on the travel-related introduction of HIV-1 non-B subtypes to Northern Poland

Phylodynamic, sequence data based reconstructions for the surveillance of the geographic spatial spread are a powerful tool in molecular epidemiology. In this study region of origin for the set of 57 partial pol sequences derived from the patients the history of travel-related HIV transmission was analyzed using phylogeographic approach. Maximum likelihood trees based on the sets of country-annotated reference sequences were inferred for identified non-B variants. Region of sequence import was assigned using on the highest approximate likelihood ratios. Import of the A1 clades was traced to the Eastern Europe and associated with immigration from this region. Subtype C infections clustered most frequently with sequences of the South African origin while majority of subtype Ds were similar to the European clades. Subtype G sequences clustered with Portuguese lineage, CRF01_AE with Eastern or South-Eastern Asian. Eastern European, Middle African or Western African lineage was assigned for the CFR02_AG. Rare circulating recombinants originated either from Central Africa (CRF11_cpx - Democratic Republic of Congo, CRF13_cpx - Central African Republic, CRF37_cpx - Cameroon) or South America (CRF28_BF and CRF46_BF - Brazil). Import of the HIV-1 non-B variants, including recombinant forms previously rarely found in Poland and Europe is frequent among travelers. Observed founder events result in the heterosexually-driven introduction of the novel HIV-1 variants into the population.

Distinct HIV type 1 strains in different risk groups and the absence of new infections by drug-resistant strains in Lithuania

To analyze HIV-1 genotypes in Lithuania and the transmission of drug-resistant viruses, HIV-1 sequences were obtained from 138 individuals, who were diagnosed as HIV-1 infected in 1990-2008 and represented all major risk groups. Subtype A strains, dominating in the former Soviet Union (90% of cases), were found in 60% of individuals, followed by subtype B (22%) and CRF03_AB (12%) strains. The remaining 7% of the strains included variants belonging to subtype C, CRF01_AE, CRF02_AG, more complex recombinant forms, and strains that could not be reliably genotyped. Analysis of virus genotypes per risk group revealed the circulation of distinct HIV-1 strains in different risk groups: subtype A viruses were present in 82% of injecting drug users (IDUs), but less than a half of heterosexually infected individuals and cases with unknown transmission route, and none of men having sex with men (MSM). We observed no mutations causing drug resistance among 27 newly diagnosed HIV-1 cases.

Transmission of risk-group specific HIV-1 strains among Dutch drug users for more than 20 years and their replacement by nonspecific strains after switching to low-harm drug practices

OBJECTIVES

To characterize HIV-1 epidemiological networks of men having sex with men (MSM) and drug users (DUs) in the Netherlands for >30 years.

DESIGN AND METHODS

Previously, we demonstrated different origin of the HIV-1 epidemics in Dutch MSM and DUs. To achieve the study objectives, risk group-specific genetic markers in the pol gene were examined in 315 participants of the Amsterdam Cohort Studies on HIV/AIDS who were registered as HIV-1 infected in 1981-2011.

RESULTS

Phylogenetic analysis demonstrated circulation of distinct virus strains in the 2 networks, with 98% of viruses of MSM clustering together and apart from strains of 73% DUs. Nine genetic markers that significantly distinguished virus strains specific for DUs were identified, of which 3 were ≥90% conserved. Over the total observation period, only 6% of viruses (4 of MSM and 14 of DUs) clustered with those of the other risk group. Among these sequences, the 3 most conserved genetic markers of that other risk group were 87% conserved.All 4 cases of DU-specific viruses among MSM occurred in 1980s-early 1990s. Viruses nonspecific for DUs were causing new infections among DUs at the rate of 20% till 2002 and replaced DU-specific strains among new infections thereafter, coinciding with switching of DUs to low-harm drug practices.

CONCLUSIONS

Dutch MSM and DUs have remained separate epidemiological networks for decades, despite their geographical and behavioral overlap. Switching to low-harm drug practices among DUs resulted in new infections caused by HIV-1 strains originating from other risk groups.

Molecular typing of the recently expanding subtype B HIV-1 epidemic in Romania: evidence for local spread among MSMs in Bucharest area

HIV-1 subtype B is predominant in Europe except in some countries from Eastern Europe which are characterized by a high prevalence of non-B subtypes and circulating recombinant forms (CRFs). Romania is a particular case: the HIV-1 epidemic started with subtype F1 which is still the most prevalent. Previous studies have shown an increasing prevalence of subtype B which is the second most frequent one among the newly diagnosed individuals, followed by subtype C and several CRFs as well as unique recombinant forms (URFs). Our objective was to analyze in detail the characteristics (way of dispersal, association with transmission risk groups) of the subtype B infections in Romania by means of phylogenetic analysis. Among all the individuals sampled during 2003-2010, 71 out of 1127 patients (6.3%) have been identified to be infected with subtype B strains. The most frequent route of infection identified in HIV-1 subtype B patients in Romania was MSM transmission (39.6%), followed by the heterosexual route (35.2%). Many of the patients acquired the infection abroad, mainly in Western European countries. Phylogenetic analysis indicated the existence of a local transmission network (monophyletic clade) including 14 patients, mainly MSM living in the Bucharest area. We estimate the origin of the local transmission network that dates at the beginning of the 90s; the introduction of the F1 and C subtypes occurred earlier. The rest of the sequences were intermixed with reference strains sampled across Europe suggesting that single infection were not followed by subsequent dispersal within the local population. Although HIV-1 subtype B epidemic in Romania is recent, there is evidence for local spread among the MSMs, in addition to multiple introductions.

A new unique recombinant HIV type 1 isolated from a child born to an HIV-infected mother

The new HIV-1 recombinant, with a B(gag) A(pol)A(env) structure, is described. This recombinant virus differs from the classical "Kaliningrad" (AF193276.1) virus with an A(gag) B(pol)B(env) structure. The number of new HIV cases in Belarus has been increasing in the past few years. Within the 12 months of 2010, 1069 new cases of HIV infection were registered. Molecular epidemiological investigations have shown that though HIV-1 subtype A (84.5%) still dominates in HIV/AIDS patients, the quantity of CRFs has also increased to 7.1%. Although cases with the CRF03_AB virus were previously described in patients from Belarus, CRF06_cpx and CRF02_AG are described in Belarus for the first time.

Subtype diversity associated with the development of HIV-1 resistance to integrase inhibitors

We used genotypic and phylogenetic analysis to determine integrase diversity among subtypes, and studied natural polymorphisms and mutations implicated in resistance to integrase inhibitors (INI) in treatment-naïve persons (n = 220) and -experienced individuals (n = 24). Phylogenetics revealed 7 and 10% inter-subtype diversity in the integrase and reverse transcriptase (RT)/protease regions, respectively. Integrase sequencing identified a novel A/B recombinant in which all viruses in a male-sex-male (MSM) transmission cluster (n = 12) appeared to possess subtype B in integrase and subtype A in the remainder of the pol region. Natural variations and signature polymorphisms were observed at codon positions 140, 148, 151, 157, and 160 among HIV subtypes. These variations predicted higher genetic barriers to G140S and G140C in subtypes C, CRF02_AG, and A/CRF01_AE, as well as higher genetic barriers toward acquisition of V151I in subtypes CRF02_AG and A/CRF01_AE. The E157Q and E160Q mutational motif was observed in 35% of INI-naïve patients harboring subtype C infections, indicating intra-subtype variations. Thirteen patients failed raltegravir (RAL)-containing regimens within 8 ± 1 months, in association with the major Q148K/R/H and G140A/S (n = 8/24) or N155H (n = 5/24) mutational pathways. Of note, the remaining patients on RAL regimens for 14 ± 3 months harbored no or only minor integrase mutations/polymorphisms (T66I, T97A, H114P, S119P, A124S, G163R, I203M, R263K). These results demonstrate the importance of understanding subtype variability in the development of resistance to INIs.

Epidemiological networks and drug resistance of HIV type 1 in Krasnoyarsk region, Russia

To study the molecular epidemiology of HIV-1 in Krasnoyarsk region, Russia, where HIV-1 has spread rapidly since 2000, we obtained pol sequences from individuals living in this region (n = 67) as well as in the geographically closely related Altay region (n = 13). In both regions, subtype A viruses specific for the former Soviet Union (IDU-A strains) were dominant (92.5%). Virus sequences clustered according to the geographic origin of the infected individuals rather than to their risk group, demonstrating the role of geographically defined epidemiological networks in the propagation of the HIV-1 epidemic in the region. Six viruses belonged to subtype B. Three of them were phylogenetically (and therefore epidemiologically) closely related to each other, demonstrating that even though IDU-A viruses dominate the epidemic, the spread of other virus strains does occur. Most viruses (75%) had an A62V mutation in reverse transcriptase, specific for HIV-1 strains in Russia. Remarkably, 26 of 47 (55%) patients under HAART with detectable virus loads did not have any known drug-resistant mutation, indicating the need to increase compliance to therapy.

Molecular epidemiology of HIV-1 in St Petersburg, Russia: predominance of subtype A, former Soviet Union variant, and identification of intrasubtype subclusters

OBJECTIVES

To examine HIV-1 genetic diversity in St. Petersburg.

METHODS

Partial HIV-1 pol sequences from 102 plasma samples collected in 2006 were analyzed with a Bayesian phylogeny inference method.

RESULTS

Subtype A, former Soviet Union (FSU) variant (AFSU), was the predominant clade (89.3%); other clades were subtypes B (9.7%) and F1 (1%). AFSU was predominant both among injecting drug users (98.2%) and heterosexually infected individuals (91.4%), whereas subtype B was more prevalent among homosexual men (75%). Within the AFSU variant, most sequences (93.5%) branched within 1 of 4 strongly supported subclusters. The largest comprised 63% AFSU viruses and was uncommon outside St Petersburg. A second subcluster (17.4% AFSU viruses) corresponds to the variant with the V77I substitution in protease, which is widely circulating in different FSU countries. Two minor subclusters comprised 8.7% and 6.5% AFSU viruses, respectively. There was no correlation between risk exposure and AFSU subclusters. Six of 8 subtype B sequences, 4 of them from homosexual men, grouped in a monophyletic subcluster.

CONCLUSIONS

The results of this study show a great predominance of AFSU viruses in St Petersburg and point to a few phylogenetically identifiable introductions as the origin of most current HIV-1 AFSU infections in the city.

Identifying recombinants in human and primate immunodeficiency virus sequence alignments using quartet scanning

Background

Recombination has a profound impact on the evolution of viruses, but characterizing recombination patterns in molecular sequences remains a challenging endeavor. Despite its importance in molecular evolutionary studies, identifying the sequences that exhibit such patterns has received comparatively less attention in the recombination detection framework. Here, we extend a quartet-mapping based recombination detection method to enable identification of recombinant sequences without prior specifications of either query and reference sequences. Through simulations we evaluate different recombinant identification statistics and significance tests. We compare the quartet approach with triplet-based methods that employ additional heuristic tests to identify parental and recombinant sequences.

Results

Analysis of phylogenetic simulations reveal that identifying the descendents of relatively old recombination events is a challenging task for all methods available, and that quartet scanning performs relatively well compared to the triplet based methods. The use of quartet scanning is further demonstrated by analyzing both well-established and putative HIV-1 recombinant strains. In agreement with recent findings, we provide evidence that the presumed circulating recombinant CRF02_AG is a 'pure' lineage, whereas the presumed parental lineage subtype G has a recombinant origin. We also demonstrate HIV-1 intrasubtype recombination, confirm the hybrid origin of SIV in chimpanzees and further disentangle the recombinant history of SIV lineages in a primate immunodeficiency virus data set.

Conclusion

Quartet scanning makes a valuable addition to triplet-based methods for identifying recombinant sequences without prior specifications of either query and reference sequences. The new method is available in the VisRD v.3.0 package .

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.

The AIDS epidemic, inflicted by the human immunodeficiency virus (HIV), has already claimed 25 million lives, thus posing a global threat. Since its discovery, several HIV subtypes have emerged, characterized by distinct genomic sequences and variable geographic locations. Here, we investigate the nature of the genetic differences among the subtypes. The neutral theory of evolution suggests that most genetic differences marginally affect the function of the encoded proteins (hence neutral) and thus occur randomly. Alternatively, changes in protein function are reflected by a pattern of nonrandom genetic differences. To address this issue, we developed a computational method, which studies the differences between sequences of different HIV subtypes, and estimates which of the explanations is more likely. Using a large sample of HIV protein sequences, we discovered that part of the variability among the subtypes is not random and possibly reflects different functional constraints imposed on the subtypes during the course of their evolution. An in-depth inspection of these nonrandom changes revealed a correlation with biological traits, such as drug resistance and mechanisms facilitating viral entry into the host cell. Interestingly, nonrandom changes are also characteristic of a viral strain that recently emerged in the former Soviet Union.

Identification of a possible ancestor of the subtype A1 HIV Type 1 variant circulating in the former Soviet Union

The HIV-1 subsubtype A1 epidemic of the former Soviet Union (FSU) was caused by an A1 variant apparently derived from a single introduction event. We identified an A1 variant highly related to the A1 lineage circulating in the FSU in a patient from Conakry, Republic of Guinea, who was diagnosed with HIV-1 when he moved to Italy in 2002. The most probable route of infection was two blood transfusions received in his country of origin in 1998. Full-length (9781 bp) molecular characterization revealed that this strain was evolutionarily parental to, but distinct from, the A1 lineage circulating in the FSU. Similar results were obtained analyzing partial genome sequences. A full-length sequence similarity plot and bootscanning analysis supported this evidence and excluded any potential recombination events with other HIV-1 variants. This viral strain represents the first evidence of an African patient infected by an A1 subtype related to the A1 variants spreading in FSU countries. The identification of this distinct A1 variant may support the origin of the Russian A epidemic from West Africa.

The HIV type 1 epidemic in Bulgaria involves multiple subtypes and is sustained by continuous viral inflow from West and East European countries

Little is known about the HIV-1 epidemic in Balkan countries. To fill the gap, we investigated the viral genetic diversity in Bulgaria, by sequencing and phylogenetic characterization of 86 plasma samples collected between 2002 and 2006 from seropositive individuals diagnosed within 1986-2006. Analysis of pol gene sequences assigned 51% of the samples to HIV-1 subtype B and 27% to subtype A1. HIV-1 subtype C, F, G, H, and a few putative recombinant forms were also found. Phylogenetic and molecular clock analysis showed a continuous exchange of subtype A and B between Bulgaria and Western as well as other Eastern European countries. At least three separate introductions of HIV-1 subtype A and four of HIV-1 subtype B have occurred within the past 25 years in Bulgaria. The central geographic location of Bulgaria, the substantial genetic heterogeneity of the epidemic with multiple subtypes, and the significant viral flow observed to and from the Balkan countries have the potential to modify the current HIV-1 epidemiological structure in Europe and highlight the importance of more extensive and continuous monitoring of the epidemic in the Balkans.

High-resolution molecular epidemiology and evolutionary history of HIV-1 subtypes in Albania

Background

HIV-1 epidemic in Western Europe is largely due to subtype B. Little is known about the HIV-1 in Eastern Europe, but a few studies have shown that non-B subtypes are quite common. In Albania, where a recent study estimated a ten-fold increase of AIDS incidence during the last six years, subtype A and B account for 90% of the know infections.

Methodology/Principal Findings

We investigated the demographic history of HIV-1 subtype A and B in Albania by using a statistical framework based on coalescent theory and phylogeography. High-resolution phylogenetic and molecular clock analysis showed a limited introduction to the Balkan country of subtype A during the late 1980s followed by an epidemic outburst in the early 1990s. In contrast, subtype B was apparently introduced multiple times between the mid-1970s and mid-1980s. Both subtypes are growing exponentially, although the HIV-1A epidemic displays a faster growth rate, and a significantly higher basic reproductive number R₀. HIV-1A gene flow occurs primarily from the capital Tirane, in the center of the country, to the periphery, while HIV-1B flow is characterized by a balanced exchange between center and periphery. Finally, we calculated that the actual number of infections in Albania is at least two orders of magnitude higher than previously thought.

Conclusions/Significance

Our analysis demonstrates the power of recently developed computational tools to investigate molecular epidemiology of pathogens, and emphasize the complex factors involved in the establishment of HIV-1 epidemics. We suggest that a significant correlation exists between HIV-1 exponential spread and the socio-political changes occurred during the Balkan wars. The fast growth of a relatively new non-B epidemic in the Balkans may have significant consequences for the evolution of HIV-1 epidemiology in neighboring countries in Eastern and Western Europe.

The feasibility of HIV vaccine efficacy trials among Russian injection drug users

IDU exposure remains a primary driver of the Russian HIV epidemic, and recent incidence data provide little evidence that this epidemic is slowing. While there are multiple important challenges that need to be further explored before starting vaccine trials, most importantly access to evidence-based drug treatment services for trial participants, the current context of high HIV incidence and low genetic diversity of HIV strains, suggests the need for intensified prevention strategies and supports the feasibility of mounting efficacy trials of HIV vaccines among IDUs in the Russian Federation.

Identification of a novel circulating recombinant form (CRF) 36_cpx in Cameroon that combines two CRFs (01_AE and 02_AG) with ancestral lineages of subtypes A and G

An array of CRFs have been identified in Cameroon, the most notable being CRF02_AG. HIV-1 in the East Province of Cameroon is particularly diverse: in a recent study, we found a high proportion of unique recombinant forms (URFs). Herein we describe the analysis of the full-length sequences of two of these URFs, which, after preliminary analysis of gag, pol, and env fragments, appeared to be a novel CRF. This novel strain, CRF36_cpx, contains fragments that can be assigned to the CRF01_AE, CRF02_AG, and subtype A and G radiations. Forty percent of the genome can be classified as CRF02_AG, including regions in gag, pol, env, and the accessory genes. Twenty-seven percent is CRF01_AE, comprising the majority of gag, the beginning of env, and the end of env into the 3' LTR. Twenty percent of the genome can be assigned to subtype A, with segments in pol and env. The remaining 13% of the sequence is classifiable as subtype G, in pol and vpu. The subtype A and G lineages formed by the CRF36_cpx sequences are unique and appear ancestral in nature. CRF36_cpx is both the first to combine more than one CRF and the first to include fragments of CRF02_AG. The ancestral sequences present in CRF36_cpx represent a link to extinct strains, and, potentially, insight into the evolution of HIV-1.

Molecular determinants of HIV-1 intersubtype recombination potential

Sequence differences in the dimerization initiation signal (DIS) affect the rate of recombination between subtype B and subtype C HIV-1. To test the hypothesis that DIS sequences can be used to predict intersubtype recombination potentials, we measured the recombination rate between CRF01_A/E (AE) and B, which contain mismatches in the DIS, and between AE and C, which have an identical DIS. Compared with the intrasubtype recombination rate, the recombination rate between AE and subtype B virus was 9-fold lower, and the rate between AE and subtype C virus was 2-fold lower. Thus, DIS sequences can be used to predict the recombination potential between HIV-1 subtypes. Further analyses revealed that the 2-fold lower recombination rate between AE and C viruses can be restored to the intrasubtype recombination rate by matching a part of the LTR and a portion of the viral genome. Therefore, the lower intersubtype recombination rate between AE and C is not caused by a given region but is a cumulative effect by more than one region.

The rate of epidemiological and virological changes during the transition from nascent to concentrated HIV epidemic stage in the former Soviet Union countries

The rate of processes accompanying the transition of the HIV-1 epidemic from nascent stage to concentrated one in the Former Soviet Union (FSU) during intravenous drug user (IDU)-associated HIV infection outbreaks in 1994-1999 has not been analyzed. To define the rates, we studied susceptible populations and circulating viruses before, during, and after the outbreaks. Our findings included the following: (1) the pattern of high HIV-1 genetic diversity characteristic of the nascent epidemic changed to a concentrated one within 1 year in St. Petersburg and in Moscow; (2) different FSU regions were at different stages of the HIV-1 epidemic in 1994-1996; (3) the change of serotypic patterns characteristic of different stages of the HIV/AIDS epidemic for the non-IDU risk group occurred within 1 year in Moscow, suggesting an extremely high rate of IDU-associated epidemic pattern distributions in regions and susceptible populations in the FSU.

A diversidade do HIV-1: uma ferramenta para o estudo da pandemia

Uma das características mais marcantes do HIV-1 é a imensa diversidade observada entre as cepas que compõem a pandemia de HIV/AIDS. Na última década, a classificação das variantes do vírus em grupos, subtipos e formas recombinantes circulantes (CRF) e a observação de padrões específicos de mutação têm provado serem ferramentas poderosas para os estudos da dinâmica molecular do vírus. O acompanhamento da distribuição mundial da diversidade do HIV-1 tem sido empregado, por exemplo, em programas de vigilância epidemiológica, bem como na reconstrução da história de epidemias regionais. Além disto, a observação de padrões específicos de distribuição espacial do vírus sugere a existência de diferenças na patogenia e transmissibilidade entre os diversos subtipos. A análise molecular das seqüências do vírus também permite a estimativa do tempo de divergência entre as variantes e das forças dinâmicas que modelam as árvores filogenéticas.

Analysis of drug resistance-associated mutations in treatment-naïve individuals infected with different genetic forms of HIV-1 circulating in countries of the former Soviet Union

There are few data on drug resistance-associated mutations in the former Soviet Union since, studies have usually been focused on the env or gag genes for subtype information. This study examines the prevalence and patterns of resistance-associated mutations to reverse transcriptase and protease inhibitors (RTI, PRI) in 278 HIV-1-infected treatment-naïve subjects from countries of Eastern Europe, and defines characteristic polymorphisms of RT and PR sequences in HIV-1 subtype A viruses. Blood samples were collected between 1997 and 2004. Plasma RNA was used for PR-RT amplification by reverse transcription coupled with nested PCR and sequencing. Phylogenetic analysis was done with neighbor-joining trees and bootscanning. Analysis of drug resistance mutations, with Stanford University HIV Drug Resistance Database's algorithm, resulted in an overall prevalence of 12.9% resistance to RTI and 3.9% to PRI. The most frequent substitutions in the RT region were at positions 62 and 236. V77I substitution in PR was found in 47.8% of samples. Polymorphisms in subtype A sequences were identified. This is the first study reporting the prevalence and patterns of both PRI and RTI resistance-associated mutations in naïve HIV-1 infected patients from the former Soviet Union. These data underline the importance of genotypic resistance testing of chronically HIV-1-infected patients before initiating treatment, in order to select the most suitable drug regimen.

The HIV type 1 epidemic in Belarus: predominance of Eastern European subtype A strains and circulation of subtype B viruses

To study the molecular epidemiology of HIV-1 in Belarus, where the rapid spread of HIV-1 has been registered since 1996, we obtained HIV-1 sequences from 30 individuals living in five cities in both the main geographic areas of the epidemic (Gomel and Minsk regions) and territories where spreading of the epidemic remains limited (Grodno region). Analysis of env V3 and gag p17/p24 sequences demonstrated that infections in all 12 injecting drug users and 14 of 18 individuals infected through sexual contacts were caused by subtype A viruses that are specific for the epidemic in the former Soviet Union (IDU-A viruses), while the remaining four infections were caused by phylogenetically unrelated to each other subtype B viruses. Extrapolation of these results to the total population of HIV-1-infected individuals in Belarus allowed us to estimate that IDU-A viruses account for nearly 95% of HIV-1 infections in Belarus.

Temporal trends in the HIV-1 epidemic in Russia: predominance of subtype A

During the period 1996-1997, three highly homogeneous variants of HIV-1 were identified, circulating among injecting drug users (IDUs) in the former Soviet Union republics. One of these belonged to HIV-1 genetic subtype A (IDU-A), another belonged to HIV-1 genetic subtype B (IDU-B) and the third was a recombinant between the first two variants (CRF03_AB). However, since 1997, the HIV-1 epidemic has affected an increasing number of geographic regions in Russia. This study was undertaken to survey the prevailing genetic variants and to estimate the current proportions of these three HIV-1 genetic subtypes in Russia. Blood samples were taken in 1999-2003 from 1090 HIV-infected individuals and analysed by gag/env HMA. The IDU-A variant was found to be the majority variant (89.7-100%) in 44 of 45 regions of the Russian Federation studied. The IDU-A variant was also found to spreading rapidly through heterosexual transmission in 1999-2003 (30/34, 88%). CRF03_AB predominates in the Kaliningrad region only (28/29, 96.6%). The IDU-B variant is currently of minor importance in the IDU epidemic but other European subtype B variants predominate among men having sex with men (18/18, 100%). Sequence analysis of the env V3 encoding regions derived from HIV-1 infected individuals in Yekaterinburg (the main centre of the HIV-1 epidemic in Russia in 2002-2003) showed that the IDU-A variant is still highly homogeneous. The mean pairwise nucleotide distance (n = 9) was 2.89 +/- 1.14 (range 1.36-6.14). However, the mean genetic distance between each sequence within the samples collected from the Yekaterinburg IDU-A variant subset and the IDU-A consensus is 2.51 +/- 1.06 (range 1.36-4.66) and considerably higher than in South Russia in 1996 (0.79 +/- 0.51, range 0.38-1.90). The current HIV-1 epidemic in Russia is almost entirely caused by a highly homogeneous A-subtype strain, which will influence vaccine development strategies and must be taken into account in the quality control of molecular tests for the diagnosis of HIV-1.

Molecular cloning and analysis of full-length genome of HIV type 1 strains prevalent in countries of the former Soviet Union

The HIV-1 epidemic among injecting drug users (IDUs) in countries of the former Soviet Union (FSU) was caused mainly by two HIV-1 variants: subtype A and CRF03-AB. To date only three full-length HIV-I genomes from the FSU have been sequenced: one subtype A from Byelorussia and two CRF03-AB from Russia. We report the full-length genome cloning and analysis of two more HIV-1 strains from the FSU countries (98UA0116 of subtype A and 98BY10443 of CRF03-AB). Isolate 98UA0116 is the second cloned and sequenced full-length HIV-1 genome of subtype A lineage from the FSU, which may be a novel subsubtype within sub-type A. Isolate 98BY10443 is the third full-length HIV-1 genome of CRF03-AB in the world to be cloned and sequenced. Additionally, it is the first CRF03-AB strain discovered in Byelorussia. Cloned genomic sequences of the FSU HIV-1 isolates are being used for the development of a region-specific HIV-1 vaccine.

Human immunodeficiency virus in Uzbekistan: epidemiological and genetic analyses

This study investigates the molecular epidemiology of HIV in Uzbekistan--a former Soviet Union (FSU) country located in central Asia. A total of 18,910,370 subjects were involved in an HIV serological examination through a population survey conducted in 1987-2002. Rapid changes in epidemiological dynamics and transmission modes have been observed since 1999: incidence rose from 25 newly HIV-infected subjects per year to more than 100 new cases per month within the first half of 2002, and the rate of intravenous drug use (IVDU)-associated HIV infection increased to 75% per year during the same period. Thirty HIV-1 strains, isolated from specimens obtained in 1999-2000, were directly sequenced in the env region. Phylogenetic analysis revealed a relationship to genotype A in 56.7%, and to 03_CRFAB in 13.3%; both variants have been previously reported in the FSU. The majority (85.7%) of these strains were isolated from IVDUs. The demographic history of the most prevalent HIV strain in Uzbekistan was inferred from reconstructed molecular phylogenies; exponential growth of the viral population size was thus observed to occur after the mid-1990s. In summary, detectable HIV seroprevalence remains low in the general population of Uzbekistan. However, the current study demonstrates a substantially increasing number of new infections, in association with IVDU, and an exponentially growing effective population size of the IVDU-associated HIV strain.

Spreading of HIV-1 subtype G and envB/gagG recombinant strains among injecting drug users in Lisbon, Portugal

We have evaluated the genetic diversity of HIV-1 strains infecting injecting drug users (IDUs) in Lisbon, Portugal. Heteroduplex mobility assay and/or phylogenetic analysis revealed that env (C2V3C3 or gp41) subtype B is present in 63.7% of the 135 viral samples studied, followed by subtypes G (23.7%), A (6.7%), F (5.2%), and D (0.7%). Similar analysis of gag (p24/p7) performed on 91 of the specimens demonstrated that 49.5% of the infections were caused by subtype G viruses; other gag subtypes identified were B (39.5%), F (3.3%), A and D (1.1.% each), and the recombinant circulating form CRF02_AG (5.5%). Discordant env/gag sub-types were detected in 34.1% of the strains and may reflect the presence of dual infections and/or recombinant viruses. The presumptive B/G recombinant form was highly predominant (21 of 31). The genetic pattern of HIV-1 subtype B and G strains is suggestive of multiple introductions and recombination episodes and of a longstanding presence of both subtypes in the country. C2V3C3 amino acid sequences from IDU-derived subtype G viruses presented highly significant signatures, which distinguish the variants from this transmission group. The unusually high prevalence of subtype G sequences (34.1%), independent of the geographic origin of the infected individuals, makes this IDU HIV-1 epidemic unique.

Precise mapping of recombination breakpoints suggests a common parent of two BC recombinant HIV type 1 strains circulating in China

Two different BC recombinant HIV-1 strains have arisen and begun to circulate among intravenous drug users in China. The recombinants are mostly subtype C with a few small subtype B segments. Additional full-genome sequences of the two recombinants, termed CRF07_BC and CRFO&_BC, are now available for analysis. Four CRF07_BC strains, including c54, 97CNU01, 98CN009, and a new strain CNGL-179, described here, and four CRF08_BC strains, including 97CNGX-6, 97CNGX-7, 97CNGX-9, and 98CN006, were compared for their recombination breakpoints by bootscanning and software for fine mapping of recombinants. The four CRF07_BC strains shared an identical recombination structure and the four CRF08_BC strains shared an identical, but different, recombination structure. The two CRFs share five precise subtype B/C boundaries, although although other segments differ between them, suggesting that they shared a common ancestor, itself a BC recombinant that separately "back-crossed" onto different subtype C strains. Both CRFs are broadly distributed from north to south in western China and have maintained low interpatient diversity.

Simultaneous introduction of HIV type 1 subtype A and B viruses into injecting drug users in southern Ukraine at the beginning of the epidemic in the former Soviet Union

The vast majority of HIV-1 strains from the epidemic in the former Soviet Union (FSU) belong to subtype A (IDU-A) and CRF03_AB (IDU-A/B), for which IDU-A is one of parental strains; no epidemic by another parental virus, belonging to subtype B (IDU-B), has yet been identified. To characterize viruses present during the early stage of the epidemic in southern Ukraine, where the first outbreaks in the FSU were registered, we obtained partial env and pol sequences from IDUs from Odessa and Nikolaev and compared them with viruses from other outbreaks. All viruses from Odessa belonged to the IDU-A type, which is in accord with previous studies. At the same time, we found that the outbreak in Nikolaev was caused by IDU-B viruses, indicating that this outbreak is the result of an independent virus introduction. Phylogenetic analysis of viruses from the FSU supported the epidemiological data suggesting that the HIV-1 epidemic in the FSU started in southern Ukraine.

On-going generation of multiple forms of HIV-1 intersubtype recombinants in the Yunnan Province of China

OBJECTIVES

To investigate the molecular epidemiology of HIV in China's Yunnan Province, where the initial HIV-1 outbreak among injecting drug users (IDU) occurred in 1989, and to analyse the genesis and interrelationship of the epidemic with that in surrounding areas.

DESIGN

A molecular epidemiological investigation was conducted among IDU in three prefectures in Yunnan Province, including Wenshan (east), Honghe (southeast) and Dehong (west).

METHODS

Thirty-nine specimens were collected from consenting IDU in 2000-2001. The nucleotide sequences of 2.6 kb gag-RT and 340 base pair (bp) env (C2/V3) regions were determined. Phylogenetic tree and recombination breakpoint analyses were performed.

RESULTS

The circulating recombinant form (CRF), CRF08_BC, predominated in east Yunnan near Guangxi Province (89% in Wenshan and 81% in Honghe), whereas it was not detected in Dehong (0/14) in the west. In contrast, 71% (10/14) of the Dehong isolates were unique recombinant forms (URF), mostly between subtypes B' (Thailand variant of subtype B) and C, with distinct profiles of recombination breakpoints. The subtype B' accounts for the remaining 29% (4/14) of Dehong isolates. Interestingly, two Honghe isolates (2/16) shared some of the precise B'/C recombination breakpoints with CRF07_BC.

CONCLUSION

New recombinant strains are arising continually in west Yunnan near the Myanmar border. Some appeared to be secondary recombinants derived from CRF07_BC that had further recombined with other strains. The uneven distribution of subtypes, CRF and URF, suggests the presence of independent transmission networks and clusters among IDU in Yunnan.

Prevalence and origin of HIV-1 group M subtypes among patients attending a Belgian hospital in 1999

HIV-1 group M strains are usually subtyped based on gag and/or env gene sequences. In our lab, part of the pol gene sequence was available in order to determine the genotypic anti-HIV drug resistance profile. To estimate the prevalence of the different HIV-1 subtypes in patients visiting the University Hospitals in Leuven in 1999 and for whom a genotypic drug resistance test was needed, we tried to use the pol sequence for subtyping. Recombination was investigated by similarity plots and bootscanning and subtyping was performed by phylogenetic analysis. The overall region spanning the entire protease and 747 nucleotides of the reverse transcriptase proved very suitable for subtyping, although there was a low phylogenetic signal at the beginning of the reverse transcriptase (nucleotides 0-250), as we demonstrated by likelihood mapping. Of the 41 samples analyzed, 21 belonged to subtype B. Of the other 20 non-B strains, 9 belonged to subtype C, 2 to subtype D and 1 to subtype A, G, H and J, respectively, 3 were CRF_02 (Circulating Recombinant Form), 1 was recombinant with a novel breakpoint and 1 sample was untypable. Although subtype B is still the most prevalent subtype in Belgium, it seems to be responsible for only half of the infections in this study. We could also document that the prevalence of subtype C is high in the Belgian native patients, especially among the heterosexually infected population. This could possibly be an indication for an epidemic spread of HIV-1 subtype C in Belgium, as for one third of these patients, no link to an endemic region could be found. The other non-B subtypes and the recombinants are mainly introduced by immigrants or by Belgian citizens traveling abroad.

CRF06-cpx: a new circulating recombinant form of HIV-1 in West Africa involving subtypes A, G, K, and J

Phylogenetic analysis of numerous strains of HIV-1 isolated from diverse geographic origins has revealed three distinct groups of HIV-1: groups M, N, and O. Within group M, subtypes, sub-subtypes and circulating recombinant forms (CRFs) exist. Recently, two near-full-length genomes of similar complex mosaic viruses containing fragments of subtypes A, G, I, and J were described in patients from Burkina Faso (BFP-90) and Mali (95ML-84). Here, we report on the characterization of two additional full-length genome sequences with similar mosaic structure in epidemiologically unlinked individuals from Senegal (97SE-1078) and Mali (95ML-127). Phylogenetic and recombinant analysis confirmed that the previously described strains, BFP-90 and 95ML-84, were indeed a new CRF of HIV-1, which we can now designate as CRF06-cpx. This new CRF fits the complex (cpx) designation, because four different subtypes (A, G, K, and J) were involved in the mosaic genome structure. The fragment in the pol gene, which was initially characterized as unknown in the BFP-90 strain and subsequently as subtype I in the 95ML-84 strain, is now, with the recent description of the new K subtype, clearly identified as subtype K. CRF06-cpx circulates in Senegal, Mali, Burkina Faso, Ivory Coast, and Nigeria, although the exact prevalence remains to be determined. Importantly, this new variant has also been documented on other continents (Europe [France] and Australia), showing that these viruses are spreading not only locally but globally.

Molecular epidemiology of HIV type 1 infection in Portugal: high prevalence of non-B subtypes

In this study, we have investigated the diversity of current HIV-1 strains circulating in the metropolitan area of Lisbon, Portugal. A total of 217 HIV-1-positive blood samples, collected between October 1998 and December 2000, was genetically characterized in the gp120 C2V3C3 region (n = 205) or part of the gp41 N-terminal segment (n = 12) by heteroduplex mobility assay (HMA) and/or DNA sequencing. The HMA subtyping efficiency (number of samples unambiguously subtyped by HMA divided by the total number of samples subtyped) was 65.9% (143 of 217), with indeterminate migration patterns of subtype A and G strains contributing significantly to this value. On the overall, subtype B was the most prevalent (50.2%), followed by subtypes G (21.7%), A (17.5%), and F (5.5%), whereas subtypes C, D, H, and J accounted altogether for 5.1% of the infections. Non-B subtypes were responsible for 77.4 and 33.1% of the infections among African immigrants and Portuguese subjects, respectively. Angolan individuals (n = 25) were the only ones infected with all the HIV-1 subtypes documented, probably reflecting a high degree of viral genetic diversification in their country of origin. Phylogenetic analysis showed a predominance of IbNG-like viruses among subtype A sequences and two new major subclusters within subtype G (G(P) and G(P)'). The majority of the Portuguese G sequences described formed a well-defined subcluster (G(P)), supported by bootstrap values >90%, phylogenetically distant from clade G sequences in databases. gag (p24/p7) sequence analysis of these variants confirmed the maintenance of the subtype G subclusters. The multiple subclustering observed for the major clades A, B, D, and G, as well as the variety of subtypes found, indicate a high diversity of HIV-1 variants circulating in Portugal and suggest a need for continuous epidemiologic surveillance.

Human immunodeficiency virus type 1 reverse-transcriptase and protease subtypes: classification, amino acid mutation patterns, and prevalence in a northern California clinic-based population

Phylogenetic analysis of the reverse transcriptase (RT) and protease of 117 published complete human immunodeficiency virus (HIV) type 1 genome sequences demonstrated that these genes cluster into distinct subtypes. There was a slightly higher proportion of informative sites in the RT (40.4%) than in the protease (34.8%; P= .03). Although most variation between subtypes was due to synonymous nucleotide substitutions, several subtype-specific amino acid patterns were observed. In the protease, the subtype-specific variants included 7 positions associated with drug resistance. Variants at positions 10, 20, 36, and 82 were more common in non-B isolates, whereas variants at positions 63, 77, and 93 were more common in subtype B isolates. In the RT, the subtype-specific mutations did not include positions associated with anti-retroviral drug resistance. RT and protease sequences from 2246 HIV-infected persons in northern California were also examined: 99.4% of the sequences clustered with subtype B, whereas 0.6% clustered with subtype A, C, or D.

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

An explosive epidemic of human immunodeficiency virus type 1 (HIV-1) has been documented among injecting drug users (IDUs) in the former Soviet Union republics. In 1999, the two largest local IDU outbreaks of HIV-1 infection in the Russian Federation were registered in the Moscow and Irkutsk regions, where 13,004 HIV-1 cases were identified (44% of the total number of HIV-1 infections in Russia in 1999). To study the prevailing genetic variants and to estimate the genetic diversity of HIV-1 in these outbreaks, 60 samples from Moscow (n = 36) and from Irkutsk (n = 24) were analyzed using the gag/env heteroduplex mobility assay, and the env gp120 V3 encoding regions obtained from 23 individuals were sequenced. Both virus populations were highly homogeneous (the means of pairwise nucleotide distance were 1.75 +/- 0.83 and 2.35 +/- 1.59 for Irkutsk and Moscow, respectively), and similar to the subtype A viruses obtained earlier from IDUs in the former Soviet Union. The subtype A HIV-1 variant thus dominates in the largest HIV-1 outbreaks among IDUs in the Russian Federation.

A recent outbreak of human immunodeficiency virus type 1 infection in southern China was initiated by two highly homogeneous, geographically separated strains, circulating recombinant form AE and a novel BC recombinant

New outbreaks of human immunodeficiency virus type 1 (HIV-1) among injecting drug users (IDUs) are spreading in China along heroin trafficking routes. Recently, two separate HIV-1 epidemics among IDUs were reported in Guangxi, Southern China, where partial sequencing of the env gene showed subtype C and circulating recombinant form (CRF) AE. We evaluated five virtually full-length HIV-1 genome sequences from IDUs in Guangxi to determine the genetic diversity and the presence of intersubtype recombinants. Sequence analysis showed two geographically separated, highly homogeneous HIV-1 strains. B/C intersubtype recombinants were found in three IDUs from Baise City, in a mountainous region near the Yunnan-Guangxi border. These were mostly subtype C, with portions of the capsid and reverse transcriptase (RT) genes from subtype B. The subtype B portion of the capsid was located in the N-terminal domain, which has been shown to influence virus core maturation, virus infectivity, and binding to cyclophilin A, whereas the subtype B portion of RT was located in the palm subdomain, which is the active site of the enzyme. These BC recombinants differed from a BC recombinant found in Xinjiang Province in northwestern China. CRF AE strains were found in IDUs from Nanning, the capital of Guangxi, and in IDUs from Pingxiang City near the China-Vietnam border. The AE and BC recombinants were both remarkable for their low interpatient diversity, less than 1% for the full genome. Rapid spread of HIV-1 among IDUs may foster the emergence of highly homogeneous strains, including novel recombinants in regions with multiple subtypes.