Molecular epidemiology of an HIV-1 subtype A subcluster among injection drug users in the Southern Ukraine

Phylodynamic evolution of HIV-1 A6 sub-subtype epidemics in Poland

The human immunodeficiency virus type 1 (HIV-1) A6 sub-subtype is highly prevalent in Eastern Europe. Over the past decade, the dissemination of the A6 lineage has been expanding in Poland. The recent Russian invasion of Ukraine may further escalate the spread of this sub-subtype. While evolutionary studies using viral sequences have been instrumental in identifying the HIV epidemic patterns, the origins, and dynamics of the A6 sub-subtype in Poland remain to be explored. We analyzed 1185 HIV-1 A6 pol sequences from Poland, along with 8318 publicly available sequences from other countries. For analyses, phylogenetic tree construction, population dynamics inference, Bayesian analysis, and discrete phylogeographic modeling were employed. Of the introduction events to Poland, 69.94% originated from Ukraine, followed by 29.17% from Russia. Most A6 sequences in Poland (53.16%) formed four large clades, with their introductions spanning 1993-2008. Central and Southern Polish regions significantly influenced migration events. Transmissions among men who have sex with men (MSM) emerged as the dominant risk group for virus circulation, representing 72.92% of migration events. Sequences from migrants were found primarily outside the large clades. Past migration from Ukraine has fueled the spread of the A6 sub-subtype and the current influx of war-displaced people maintains the growing national epidemic.

HIV-1 Drug Resistance among Treatment-Naïve Patients in Russia: Analysis of the National Database, 2006-2022

In Russia, antiretroviral therapy (ART) coverage has significantly increased, which, in the absence of routine genotyping testing, could lead to an increase in HIV drug resistance (DR). The aim of this study was to investigate the patterns and temporal trends in HIV DR as well as the prevalence of genetic variants in treatment-naïve patients from 2006 to 2022, using data from the Russian database (4481 protease and reverse transcriptase and 844 integrase gene sequences). HIV genetic variants, and DR and DR mutations (DRMs) were determined using the Stanford Database. The analysis showed high viral diversity, with the predominance of A6 (78.4%), which was the most common in all transmission risk groups. The overall prevalence of surveillance DRMs (SDRMs) was 5.4%, and it reached 10.0% in 2022. Most patients harbored NNRTI SDRMs (3.3%). The prevalence of SDRMs was highest in the Ural (7.9%). Male gender and the CRF63_02A6 variant were association factors with SDRMs. The overall prevalence of DR was 12.7% and increased over time, primarily due to NNRTIs. Because baseline HIV genotyping is unavailable in Russia, it is necessary to conduct surveillance of HIV DR due to the increased ART coverage and DR prevalence. Centralized collection and unified analysis of all received genotypes in the national database can help in understanding the patterns and trends in DR to improve treatment protocols and increase the effectiveness of ART. Moreover, using the national database can help identify regions or transmission risk groups with a high prevalence of HIV DR for epidemiological measures to prevent the spread of HIV DR in the country.

Molecular Epidemiology of HIV-1 in Eastern Europe and Russia

The HIV epidemic in Eastern Europe and Russia is large and not well-controlled. To describe the more recent molecular epidemiology of HIV-1, transmitted drug resistance, and the relationship between the epidemics in this region, we sequenced the protease and reverse transcriptase genes of HIV-1 from 812 people living with HIV from Ukraine (n = 191), Georgia (n = 201), and Russia (n = 420) before the initiation of antiretroviral therapy. In 190 Ukrainian patients, the integrase gene sequence was also determined. The most reported route of transmission was heterosexual contact, followed by intravenous drug use, and men having sex with men (MSM). Several pre-existing drug resistance mutations were found against non-nucleoside reverse transcriptase inhibitors (RTIs) (n = 103), protease inhibitors (n = 11), and nucleoside analogue RTIs (n = 12), mostly polymorphic mutations or revertants. In the integrase gene, four strains with accessory integrase strand transfer inhibitor mutations were identified. Sub-subtype A6 caused most of the infections (713/812; 87.8%) in all three countries, including in MSM. In contrast to earlier studies, no clear clusters related to the route of transmission were identified, indicating that, within the region, the exchange of viruses among the different risk groups may occur more often than earlier reported.

Identifying HIV-1 Transmission Clusters in Uzbekistan through Analysis of Molecular Surveillance Data

The CRF02_AG and sub-subtype A6 are currently the predominant HIV-1 variants in the Republic of Uzbekistan, but little is known about their time-spatial clustering patterns in high-risk populations. We have applied molecular evolution methods and network analyses to better understand the transmission patterns of these subtypes by analyzing 316 pol sequences obtained during the surveillance study of HIV drug resistance. Network analysis showed that about one third of the HIV infected persons were organized into clusters, including large clusters with more than 35 members. These clusters were composed mostly of injecting drug users and/or heterosexuals, with women having mainly high centrality within networks identified in both subtypes. Phylogenetic analyses of the 'Uzbek' sequences, including those publicly available, show that Russia and Ukraine played a role as the main sources of the current subtype A6 epidemic in the Republic. At the same time, Uzbekistan has been a local center of the CRF02_AG epidemic spread in the former USSR since the early 2000s. Both of these HIV-1 variants continue to spread in Uzbekistan, highlighting the importance of identifying transmission networks and transmission clusters to prevent further HIV spread, and the need for HIV prevention and education campaigns in high-risk groups.

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.

HIV-1 Sub-Subtype A6: Settings for Normalised Identification and Molecular Epidemiology in the Southern Federal District, Russia

Russia has one of the largest and fastest growing HIV epidemics. However, epidemiological data are scarce. Sub-subtype A6 is most prevalent in Russia but its identification is challenging. We analysed protease/reverse transcriptase-, integrase-sequences, and epidemiological data from 303 patients to develop a methodology for the systematisation of A6 identification and to describe the HIV epidemiology in the Russian Southern Federal District. Drug consumption (32.0%) and heterosexual contact (27.1%) were the major reported transmission risks. This study successfully established the settings for systematic identification of A6 samples. Low frequency of subtype B (3.3%) and large prevalence of sub-subtype A6 (69.6%) and subtype G (23.4%) were detected. Transmitted PI- (8.8%) and NRTI-resistance (6.4%) were detected in therapy-naive patients. In therapy-experienced patients, 17.3% of the isolates showed resistance to PIs, 50.0% to NRTI, 39.2% to NNRTIs, and 9.5% to INSTIs. Multiresistance was identified in 52 isolates, 40 corresponding to two-class resistance and seven to three-class resistance. Two resistance-associated-mutations significantly associated to sub-subtype A6 samples: A62V_RT and G190S_RT. This study establishes the conditions for a systematic annotation of sub-subtype A6 to normalise epidemiological studies. Accurate knowledge on South Russian epidemiology will allow for the development of efficient regional frameworks for HIV-1 infection management.

Molecular Surveillance of HIV-1 Infection in Krasnoyarsk Region, Russia: Epidemiology, Phylodynamics and Phylogeography

BACKGROUND

The information about the dynamics of the viral population and migration events that affect the epidemic in different parts of the Russia is insufficient. Possibly, the huge size of the country and limited transport accessibility to certain territories may determine unique traits of the HIV-1 evolutionary history in different regions.

OBJECTIVE

The aim of this study was to explore the genetic diversity of HIV-1 in the Krasnoyarsk region and reconstruct spatial-temporal dynamics of the infection in the region.

METHODS

The demographic and virologic data from 281 HIV-infected individuals in Krasnoyarsk region collected during 2011-2016 were analyzed. The time to the most recent common ancestor, evolutionary rates, population growth, and ancestral geographic movements was estimated using Bayesian coalescent-based methods.

RESULTS

The study revealed moderate diversity of the HIV-1 subtypes found in the region, which included A6 (92.3%), CRF063_02A (4.3%), B (1.1%), and unique recombinants (2.5%). Phylogenetic reconstruction revealed that the A6 subtype was introduced into Krasnoyarsk region by one viral lineage, which arose around 1996.9 (1994.5-1999.5). The phylogeography analysis pointed to Krasnoyarsk city as the geographical center of the epidemic, which further spread to central neighboring districts of the region. At least two epidemic growth phases of subtype A6 were identified which included exponential growth in early-2000s followed by the decline in the mid/late 2010s.

CONCLUSION

This study demonstrates a change in the genetic diversity of HIV-1 in the Krasnoyarsk region. At the beginning of the epidemic, subtype A6 prevailed, subtypes B and CRF063_02A appeared in the region later.

A Phylogenetic Analysis of Human Immunodeficiency Virus Type 1 Sequences in Kiev: Findings Among Key Populations

Background

The human immunodeficiency virus (HIV) epidemic in Ukraine has been driven by a rapid rise among people who inject drugs, but recent studies have shown an increase through sexual transmission.

Methods

Protease and reverse transcriptase sequences from 876 new HIV diagnoses (April 2013-March 2015) in Kiev were linked to demographic data. We constructed phylogenetic trees for 794 subtype A1 and 64 subtype B sequences and identified factors associated with transmission clustering. Clusters were defined as ≥2 sequences, ≥80% local branch support, and maximum genetic distance of all sequence pairs in the cluster ≤2.5%. Recent infection was determined through the limiting antigen avidity enzyme immunoassay. Sequences were analyzed for transmitted drug resistance mutations.

Results

Thirty percent of subtype A1 and 66% of subtype B sequences clustered. Large clusters (maximum 11 sequences) contained mixed risk groups. In univariate analysis, clustering was significantly associated with subtype B compared to A1 (odds ratio [OR], 4.38 [95% confidence interval {CI}, 2.56-7.50]); risk group (OR, 5.65 [95% CI, 3.27-9.75]) for men who have sex with men compared to heterosexual males; recent, compared to long-standing, infection (OR, 2.72 [95% CI, 1.64-4.52]); reported sex work contact (OR, 1.93 [95% CI, 1.07-3.47]); and younger age groups compared with age ≥36 years (OR, 1.83 [95% CI, 1.10-3.05] for age ≤25 years). Females were associated with lower odds of clustering than heterosexual males (OR, 0.49 [95% CI, .31-.77]). In multivariate analysis, risk group, subtype, and age group were independently associated with clustering (P < .001, P = .007, and P = .033, respectively). Eighteen sequences (2.1%) indicated evidence of transmitted drug resistance.

Conclusions

Our findings suggest high levels of transmission and bridging between risk groups.

Reconstruction of the Genetic History and the Current Spread of HIV-1 Subtype A in Germany

HIV-1 non-B infections have been increasing in Europe for several years. In Germany, subtype A belongs to the most abundant non-B subtypes showing an increasing prevalence of 8.3% among new infections in 2016. Here we trace the origin and examine the current spread of the German HIV-1 subtype A epidemic. Bayesian coalescence and birth-death analyses were performed with 180 German HIV-1 pol sequences and 528 related and publicly available sequences to reconstruct the population dynamics and fluctuations for each of the transmission groups. Our reconstructions indicate two distinct sources of the German subtype A epidemic, with an Eastern European and an Eastern African lineage both cocirculating in the country. A total of 13 German-origin clusters were identified; among these, 6 clusters showed recent activity. Introductions leading to further countrywide spread originated predominantly from Eastern Africa when introduced before 2005. Since 2005, however, spreading introductions have occurred exclusively within the Eastern European clade. Moreover, we observed changes in the main route of subtype A transmission. The beginning of the German epidemic (1985 to 1995) was dominated by heterosexual transmission of the Eastern African lineage. Since 2005, transmissions among German men who have sex with men (MSM) have been increasing and have been associated with the Eastern European lineage. Infections among people who inject drugs dominated between 1998 and 2005. Our findings on HIV-1 subtype A infections provide new insights into the spread of this virus and extend the understanding of the HIV epidemic in Germany.IMPORTANCE HIV-1 subtype A is the second most prevalent subtype worldwide, with a high prevalence in Eastern Africa and Eastern Europe. However, an increase of non-B infections, including subtype A infections, has been observed in Germany and other European countries. There has simultaneously been an increased flow of refugees into Europe and especially into Germany, raising the question of whether the surge in non-B infections resulted from this increased immigration or whether German transmission chains are mainly involved. This study is the first comprehensive subtype A study from a western European country analyzing in detail its phylogenetic origin, the impact of various transmission routes, and its current spread. The results and conclusions presented provide new and substantial insights for virologists, epidemiologists, and the general public health sector. In this regard, they should be useful to those authorities responsible for developing public health intervention strategies to combat the further spread of HIV/AIDS.

The Changing Epidemiological Profile of HIV-1 Subtype B Epidemic in Ukraine

While HIV-1 subtype B has caused a large epidemic in the Western world, its transmission in Ukraine remains poorly understood. We assessed the genetic diversity of HIV-1 subtype B viruses circulating in Ukraine, characterized the transmission group structure, and estimated key evolutionary and epidemiological parameters. We analyzed 120 HIV-1 subtype B pol sequences (including 46 newly generated) sampled from patients residing in 11 regions of Ukraine between 2002 and 2017. Phylogenies were estimated using maximum likelihood and Bayesian phylogenetic methods. A Bayesian molecular clock coalescent analysis was used to estimate effective population size dynamics and date the most recent common ancestors of identified clades. A phylodynamic birth-death model was used to estimate the effective reproductive number (R_e) of these clades. We identified two phylogenetically distinct predominantly Ukrainian (≥75%) clades of HIV-1 subtype B. We found no significant transmission group structure for either clade, suggesting frequent mixing among transmission groups. The estimated dates of origin of both subtype B clades were around early 1970s, similar to the introduction of HIV-1 subtype A into Ukraine. R_e was estimated to be 1.42 [95% highest posterior density (HPD) 1.26-1.56] for Clade 1 and 1.69 (95% HPD 1.49-1.84) for Clade 2. Evidently, the subtype B epidemic in the country is no longer concentrated in specific geographical regions or transmission groups. The study results highlight the necessity for strengthening preventive and monitoring efforts to reduce the further spread of HIV-1 subtype B.

HIV-1 Genetic Diversity in Recently Diagnosed Infections in Moscow: Predominance of AFSU, Frequent Branching in Clusters, and Circulation of the Iberian Subtype G Variant

HIV-1 protease-reverse transcriptase sequences from 62 HIV-1-infected individuals recently diagnosed in Moscow were analyzed. Subtype A former Soviet Union (FSU) (A_FSU) variant was the predominant clade (62.9%), followed by subtype B (22.6%), unique recombinants (6.5%), subtype G (6.5%), and CRF01_AE (1.6%). A_FSU predominated among people who inject drugs (88.9%) and heterosexually acquired infections (77.8%), while subtype B was the most prevalent genetic form among men who have sex with men (44%), although A_FSU was also frequent in this population (36%). Forty-eight (77.4%) viruses branched within intrasubtype clusters, three of which, of subtype B, had a majority of viruses collected outside of FSU. The four subtype G viruses identified in this study belonged to the Portuguese-Spanish (Iberian) variant and, together with three from databases, formed a Russian cluster closely related to viruses from Denmark. This is the first report of the circulation of the Iberian subtype G variant in Russia.

Molecular epidemiology of HIV-1 infection in Europe: An overview

Human Immunodeficiency Virus type 1 (HIV-1) is characterised by vast genetic diversity. Globally circulating HIV-1 viruses are classified into distinct phylogenetic strains (subtypes, sub-subtypes) and several recombinant forms. Here we describe the characteristics and evolution of European HIV-1 epidemic over time through a review of published literature and updated queries of existing HIV-1 sequence databases. HIV-1 in Western and Central Europe was introduced in the early-1980s in the form of subtype B, which is still the predominant clade. However, in Eastern Europe (Former Soviet Union (FSU) countries and Russia) the predominant strain, introduced into Ukraine in the mid-1990s, is subtype A (A_FSU) with transmission mostly occurring in People Who Inject Drugs (PWID). In recent years, the epidemic is evolving towards a complex tapestry with an increase in the prevalence of non-B subtypes and recombinants in Western and Central Europe. Non-B epidemics are mainly associated with immigrants, heterosexuals and females but more recently, non-B clades have also spread amongst groups where non-B strains were previously absent - non-immigrant European populations and amongst men having sex with men (MSM). In some countries, non-B clades have spread amongst the native population, for example subtype G in Portugal and subtype A in Greece, Albania and Cyprus. Romania provides a unique case where sub-subtype F1 has predominated throughout the epidemic. In contrast, HIV-1 epidemic in FSU countries remains more homogeneous with A_FSU clade predominating in all countries. The differences between the evolution of the Western epidemic and the Eastern epidemic may be attributable to differences in transmission risk behaviours, lifestyle and the patterns of human mobility. The study of HIV-1 epidemic diversity provides a useful tool by which we can understand the history of the pandemic in addition to allowing us to monitor the spread and growth of the epidemic over time.

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.

Performance of an Early Infant Diagnostic Test, AmpliSens DNA-HIV-FRT, Using Dried Blood Spots Collected from Children Born to Human Immunodeficiency Virus-Infected Mothers in Ukraine

An accurate accessible test for early infant diagnosis (EID) is crucial for identifying HIV-infected infants and linking them to treatment. To improve EID services in Ukraine, dried blood spot (DBS) samples obtained from 237 HIV-exposed children (≤18 months of age) in six regions in Ukraine in 2012 to 2013 were tested with the AmpliSens DNA-HIV-FRT assay, the Roche COBAS AmpliPrep/COBAS TaqMan (CAP/CTM) HIV-1 Qual test, and the Abbott RealTime HIV-1 Qualitative assay. In comparison with the paired whole-blood results generated from AmpliSens testing at the oblast HIV reference laboratories in Ukraine, the sensitivity was 0.99 (95% confidence interval [CI], 0.95 to 1.00) for the AmpliSens and Roche CAP/CTM Qual assays and 0.96 (95% CI, 0.90 to 0.98) for the Abbott Qualitative assay. The specificity was 1.00 (95% CI, 0.97 to 1.00) for the AmpliSens and Abbott Qualitative assays and 0.99 (95% CI, 0.96 to 1.00) for the Roche CAP/CTM Qual assay. McNemar analysis indicated that the proportions of positive results for the tests were not significantly different (P > 0.05). Cohen's kappa (0.97 to 0.99) indicated almost perfect agreement among the three tests. These results indicated that the AmpliSens DBS and whole-blood tests performed equally well and were comparable to the two commercially available EID tests. More importantly, the performance characteristics of the AmpliSens DBS test meets the World Health Organization EID test requirements; implementing AmpliSens DBS testing might improve EID services in resource-limited settings.

Enhanced HIV-1 surveillance using molecular epidemiology to study and monitor HIV-1 outbreaks among intravenous drug users (IDUs) in Athens and Bucharest

BACKGROUND

A significant increase in HIV-1 diagnoses was reported among Injecting Drug Users (IDUs) in the Athens (17-fold) and Bucharest (9-fold) metropolitan areas starting 2011.

METHODS

Molecular analyses were conducted on HIV-1 sequences from IDUs comprising 51% and 20% of the diagnosed cases among IDUs during 2011-2013 for Greece and Romania, respectively. Phylodynamic analyses were performed using the newly developed birth-death serial skyline model which allows estimating of important epidemiological parameters, as implemented in BEAST programme.

RESULTS

Most infections (>90%) occurred within four and three IDU local transmission networks in Athens and Bucharest, respectively. For all Romanian clusters, the viral strains originated from local circulating strains, whereas in Athens, the local strains seeded only two of the four sub-outbreaks. Birth-death skyline plots suggest a more explosive nature for sub-outbreaks in Bucharest than in Athens. In Athens, two sub-outbreaks had been controlled (Re<1.0) by 2013 and two appeared to be endemic (Re∼1). In Bucharest one outbreak continued to expand (Re>1.0) and two had been controlled (Re<1.0). The lead times were shorter for the outbreak in Athens than in Bucharest.

CONCLUSIONS

Enhanced molecular surveillance proved useful to gain information about the origin, causal pathways, dispersal patterns and transmission dynamics of the outbreaks that can be useful in a public health setting.

Bayesian phylogeographic analyses clarify the origin of the HIV-1 subtype A variant circulating in former Soviet Union's countries

The HIV-1 subtype A variant dominating the HIV-1 epidemics in former Soviet Union (FSU) countries (A(FSU)) represents one of the major clades of the HIV-1 pandemic. This variant was reported to have begun spreading among injecting drug users (IDUs) in the Ukrainian city of Odessa in late 1994. Two competing hypotheses have been proposed on the ancestral origin of the A(FSU) variant, locating it either in the Democratic Republic of Congo (DRC) or in the Republic of Guinea (RG). The studies supporting these hypotheses employed phylogenetic analyses to identify HIV-1 sequences collected outside FSU countries ancestrally related to A(FSU). A different approach, based on Bayesian phylogenetic inference and coalescent-based population genetics, has been employed here to elucidate the ancestry of this HIV-1 variant and to improve our knowledge on its spread in FSU countries. The analyses were carried out using env (C2-V3-C3) and p24(gag) fragments of the HIV-1 genome. The inferred migration for the HIV-1 A(FSU) variant revealed only one significantly supported migration pathway from Africa to Eastern Europe, supporting the hypothesis of its origin in the DRC and estimating the upper limit of the migration of the ancestral virus from Africa around 1970. The support for an origin in the RG was negligible. The results supported the main role of Odessa as the epicenter of the A(FSU) epidemic, dating the tMRCA of the A(FSU) variant around 1984, ten years before its explosive expansion among IDUs. The estimated origin of the AFSU subcluster responsible for the IDU outbreak was also located in Odessa, with the estimated tMRCA around 1993. Statistically supported migration routes from Odessa to other cities of Ukraine, Russia, Kazakhstan, Uzbekistan and Belarus were also inferred by the Bayesian phylogeographic analysis. These results shed new light on the origin and spread of the HIV-1 A(FSU) variant.

Short communication: Molecular epidemiology, phylogeny, and phylodynamics of CRF63_02A1, a recently originated HIV-1 circulating recombinant form spreading in Siberia

The HIV-1 epidemic in Russia is dominated by the former Soviet Union subtype A (A(FSU)) variant, but other genetic forms are circulating in the country. One is the recently described CRF63_02A1, derived from recombination between a CRF02_AG variant circulating in Central Asia and A(FSU), which has spread in the Novosibirsk region, Siberia. Here we phylogenetically analyze pol and env segments from 24 HIV-1 samples from the Novosibirsk region collected in 2013, with characterization of three new near full-length genome CRF63_02A1 sequences, and estimate the time of the most recent common ancestor (tMRCA) and the demographic growth of CRF63_02A1 using a Bayesian method. The analyses revealed that CRF63_02A1 is highly predominant in the Novosibirsk region (81.2% in pol sequences) and is transmitted both among injecting drug users and by heterosexual contact. Similarity searches with database sequences combined with phylogenetic analyses show that CRF63_02A1 is circulating in East Kazakhstan and the Eastern area of Russia bordering China. The analyses of near full-length genome sequences show that its mosaic structure is more complex than reported, with 18 breakpoints. The tMRCA of CRF63_02A1 was estimated around 2006, with exponential growth in 2008-2009 and subsequent stabilization. These results provide new insights into the molecular epidemiology, phylogeny, and phylodynamics of CRF63_02A1.

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.

Genotypic and phenotypic heterogeneity in the U3R region of HIV type 1 subtype C

Approximately 20% of all HIV-1 mother-to-child transmission (MTCT) occurs in utero (IU). In a chronic HIV infection, HIV-1 exists as a complex swarm of genetic variants, and following IU MTCT, viral genomic diversity is restricted through a mechanism that remains to be described. The 5' U3R region of the HIV-1 long terminal repeat (LTR) contains multiple transcription factor (TF) binding sites and regulates viral transcription. In this study, we tested the hypothesis that sequence polymorphisms in the U3R region of LTR are associated with IU MTCT. To this end, we used single template amplification to isolate 517 U3R sequences from maternal, placental, and infant plasma derived from 17 HIV-infected Malawian women: eight whose infants remained HIV uninfected (NT) and nine whose infants became HIV infected IU. U3R sequences show pairwise diversities ranging from 0.2% to 2.3%. U3R sequences from one participant contained two, three, or four putative NF-κB binding sites. Phylogenetic reconstructions indicated that U3R sequences from eight of nine IU participants were consistent with placental compartmentalization of HIV-1 while only one of eight NT cases was consistent with such compartmentalization. Specific TF sequence polymorphisms were not significantly associated with IU MTCT. To determine if replication efficiency of the U3R sequences was associated with IU MTCT, we cloned 90 U3R sequences and assayed promoter activity in multiple cell lines. Although we observed significant, yet highly variable promoter activity and TAT induction of promoter activity in the cell lines tested, there was no association between measured promoter activity and MTCT status. Thus, we were unable to detect a promoter genotype or phenotype associated with IU MTCT.

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.

The analysis of near full-length genome sequences of HIV type 1 subtype A viruses from Russia supports the monophyly of major intrasubtype clusters

The HIV-1 epidemic in Russia has been insufficiently studied, with only 11 complete genome sequences from this country currently available, only three of which are of the locally predominant genetic form, the former Soviet Union (FSU) subtype A variant (A(FSU)). Here we analyze 10 newly derived A(FSU) near full-length genome sequences from Russia. Samples were selected based on phylogenetic clustering in protease-reverse transcriptase in two of the major A(FSU) clusters, V77I(PR) (n=6), widely circulating in Russia and other FSU countries, and A(SP1) (n=4), predominant in St. Petersburg. The phylogenetic analysis shows that the V77I(PR) genomes group in a monophyletic cluster together with 10 previously obtained A(FSU) genome sequences from Uzbekistan, Kazakhstan, Russia, and Cyprus, all bearing the V77I substitution in protease. Similarly, the four A(SP1) genomes group in a monophyletic cluster. These results therefore show that the monophyly of V77I(PR) and A(SP1) A(FSU) clusters is supported in near complete genomes.

Short communication: Molecular epidemiology of HIV type 1 in the Republic of Dagestan, Russian Federation: virtually uniform circulation of subtype A, former Soviet Union variant, with predominance of the V77I(PR) subvariant

We examine the distribution of viral genetic forms and the presence of antiretroviral drug resistance mutations in HIV-1 infections in the Republic of Dagestan, in the North Caucasus area of Russia, where a recent large increase in HIV-1 infections has been documented. Samples were collected from 41 HIV-1-infected individuals from Dagestan, most of them from the cities of Derbent (n = 21) and Mahachkala (n = 18). Thirty six were injecting drug users and five were infected by heterosexual contact. None was on antiretroviral drug treatment. HIV-1 protease and a segment of reverse transcriptase were amplified by RT-PCR from plasma RNA and sequenced, and phylogenetic trees were constructed via maximum likelihood. Forty (97.6%) of 41 samples were of subtype A, former Soviet Union variant (A(FSU)), of which 27 (67.5%) clustered with the subvariant containing the V77I substitution in protease (V77I(PR)). Within this cluster, 13 viruses formed a local subcluster, 10 of which were from Derbent. Four viruses clustered with the A(SP2) subcluster, recently identified in St. Petersburg, two with a virus from Georgia and one with a virus from Azerbaijan. No mutations associated with antiretroviral drug resistance were detected. The results, therefore, show the relationship of the HIV-1 epidemic in Dagestan with that of other areas of Russia and of neighboring countries, and reveal the spread of the A(FSU) V77I(PR) variant in the North Caucasus area.

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.

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.

New insights into the origin of the HIV type 1 subtype A epidemic in former Soviet Union's countries derived from sequence analyses of preepidemically transmitted viruses

The HIV-1 subtype A epidemic affecting injecting drug users (IDU) in former Soviet Union (FSU) countries started dramatically in Odessa, southern Ukraine, in 1995, and is caused by a variant of monophyletic origin, often designated IDU-A. We phylogenetically analyzed one near full-length genome and two partial sequences of three HIV-1 subtype A viruses collected in St. Petersburg, Russia, heterosexually transmitted in 1992-1994. The sequences branched basally to the IDU-A clade, together with eight viruses from Odessa collected in 1993, all presumably acquired heterosexually, and two viruses from the Democratic Republic of Congo. Of all other FSU sequences in databases, only those from three recently collected viruses, one from Ukraine and two from northwestern Russia, at least one of them acquired heterosexually, branched basally to the IDU-A cluster. The results indicate that the FSU IDU-A variant derives from a strain that initially propagated heterosexually in Ukraine and originated in central Africa.

Unequal evolutionary rates in the human immunodeficiency virus type 1 (HIV-1) pandemic: the evolutionary rate of HIV-1 slows down when the epidemic rate increases

HIV-1 sequences in intravenous drug user (IDU) networks are highly homogenous even after several years, while this is not observed in most sexual epidemics. To address this disparity, we examined the human immunodeficiency virus type 1 (HIV-1) evolutionary rate on the population level for IDU and heterosexual transmissions. All available HIV-1 env V3 sequences from IDU outbreaks and heterosexual epidemics with known sampling dates were collected from the Los Alamos HIV sequence database. Evolutionary rates were calculated using phylogenetic trees with a t test root optimization of dated samples. The evolutionary rate of HIV-1 subtype A1 was found to be 8.4 times lower in fast spread among IDUs in the former Soviet Union (FSU) than in slow spread among heterosexual individuals in Africa. Mixed epidemics (IDU and heterosexual) showed intermediate evolutionary rates, indicating a combination of fast- and slow-spread patterns. Hence, if transmissions occur repeatedly during the initial stage of host infection, before selective pressures of the immune system have much impact, the rate of HIV-1 evolution on the population level will decrease. Conversely, in slow spread, where HIV-1 evolves under the pressure of the immune system before a donor infects a recipient, the virus evolution at the population level will increase. Epidemiological modeling confirmed that the evolutionary rate of HIV-1 depends on the rate of spread and predicted that the HIV-1 evolutionary rate in a fast-spreading epidemic, e.g., for IDUs in the FSU, will increase as the population becomes saturated with infections and the virus starts to spread to other risk groups.

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.

HIV epidemiology in the Northwestern Federal District of Russia: dominance of HIV type 1 subtype A

A rapidly advancing epidemic of HIV-1 infection has been documented among injecting drug users (IDUs) in Russia. The Northwestern Federal District was the first of the seven Russian Federal Districts involved in a drug-related HIV epidemic through an outbreak in Kaliningrad in 1996. The Northwestern Federal District has a high HIV prevalence rate having reached 252 per 100,000 by the end of 2003, exceeding the Russian average (180) by 1.4 times. The epidemic peaked in 2001. Since then the annual number of new cases has decreased, probably reflecting saturation among at least some IDU populations. However, at the same time, the heterosexual spread of HIV has become more prominent. To study the genetic epidemiology of HIV-1, samples were collected from 150 individuals covering a wide geographical area and different transmission groups in the Northwestern Federal District. Phylogenetic analysis revealed that an Eastern European subtype A HIV-1 strain similar to those reported earlier among IDUs in other regions of Russia accounted for 80% of HIV-1 infections and was the predominant subtype in six out of seven administrative territories studied both among IDUs and heterosexually infected persons. As an exception to the dominant role of the Eastern European subtype A strain, the CRF03-AB strain was found to be dominant in the city of Cherepovets located in the north central European Russian territory of Vologda Oblast. This is the first report of the CRF03-AB strain causing an outbreak outside the Kaliningrad region.

Molecular analysis and phylogenetic characterization of HIV in Iran

The rate of human immunodeficiency virus type 1 (HIV-1) infection in Iran has increased dramatically in the last few years. While the earliest cases were found in hemophiliacs, intravenous drug users are now fueling the outbreak. In this study, both the 122 clones of HIV-1 gag p17 and the 131 clones of env V1-V5 region were obtained from 61 HIV-1 seropositives belonging to these two groups in Iran. HIV-1 subtyping and phylogenetic analysis was done by heteroduplex mobility assays (HMA) and multiple clone sequencing. The result indicated all hemophiliacs are infected with HIV-1 subtype B and all intravenous drug users are infected with HIV-1 subtype A. Since intravenous drug abuse is the major transmission route in Iran, HIV-1 subtype A is likely to be the dominant viral subtype circulating in the country. The analysis of genetic distances showed subtype B viruses in Iran to be twice as heterogeneous as the subtype A viruses. In conclusion, this first molecular study of HIV-1 genotypes in Iran suggests two parallel outbreaks in distinct high-risk populations and may offer clues to the origin and spread of infection in Iran.

Molecular epidemiology of HIV Type 1 in Ukraine: birthplace of an epidemic

During the 1990s, HIV-1 spread rapidly through drug networks in Ukraine and from there throughout the former Soviet Union. To examine the origins of this epidemic, the genetics of HIV-1 in Ukraine were studied. Proviral DNA from PBMC was extracted and PCR amplified. Part of pol and nearly full genomes of HIV-1 were sequenced and characterized. The predominant genetic form in 163 strains was subtype A (66%), followed by subtypes B (30%), C (2%), D (1%), and a new AB recombinant form (1%). HIV strains from Kiev were diverse having subtypes A, B, C, and D. In Crimea, Donetsk, Poltava, and Odessa, however, the strains were overwhelmingly subtype A, while in Nikolaev subtype B predominated. After the near simultaneous introduction of subtypes A and B in Ukraine, subtype B remained where it was introduced while subtype A spread widely, creating the fastest growing epidemic in the world.

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.

Diversity of human immunodeficiency virus type 1 subtype A and CRF03_AB protease in Eastern Europe: selection of the V77I variant and its rapid spread in injecting drug user populations

To characterize polymorphisms of the subtype A protease in the former Soviet Union, proviral DNA samples were obtained, with informed consent, from 119 human immunodeficiency virus type 1 (HIV-1)-positive untreated injecting drug users (IDUs) from 16 regions. All individuals studied have never been treated with antiretroviral drugs. The isolates were defined as IDU-A (n = 115) and CRF03_AB (n = 4) by using gag/env HMA/sequencing. The pro region was analyzed by using sequencing and original HIV-ProteaseChip hybridization technology. The mean of pairwise nucleotide distance between 27 pro sequences (23 IDU-A and 4 CRF03_AB) was low (1.38 +/- 0.79; range, 0.00 to 3.23). All sequences contained no primary resistance mutations. However, 13 of 23 (56.5%) subtype A isolates bore the V77I substitution known as the secondary protease mutation. V77I was associated with two synonymous substitutions in triplets 31 and 78, suggesting that all V77I-bearing viruses evolved from a single source in 1997. Hybridization analysis showed that 55 of 115 (47.8%) HIV-1 isolates contained V77I, but this variant was not found in any of 31 DNA samples taken from regions, where the HIV-1 epidemic among IDUs started earlier 1997, as well as in any of four CRF03_AB isolates. The results of analysis of 12 additional samples derived from epidemiologically linked subjects showed that in all four epidemiological clusters the genotype of the donor and the recipients was the same irrespective of the route of transmission. This finding demonstrates the transmission of the V77I mutant variant, which is spreading rapidly within the circulating viral pool in Russia and Kazakhstan. The continued molecular epidemiological and virological monitoring of HIV-1 worldwide thus remains of great importance.

Two viral strains and a possible novel recombinant are responsible for the explosive injecting drug use-associated HIV type 1 epidemic in Estonia

HIV-1 infection has been rare in Estonia. In 2000, an explosive epidemic among injecting drug users was detected in the Eastern border region, resulting in 3603 newly reported cases by the end of 2003. The molecular epidemiology of the outbreak was studied to establish whether the Estonian epidemic is linked to the epidemics in Eastern Europe. Over 200 newly infected individuals were prospectively sampled from June 2000 to March 2002 in a geographically representative way, with known dates of diagnosis and information of probable route of transmission. Viral regions coding for two viral gene regions were directly sequenced from plasma viral RNA and phylogenetically analyzed. In addition, a larger region coding for the entire env gene was sequenced from one sample and studied for indications of possible recombinant structure. The Estonian HIV outbreak was found to be caused by simultaneous introduction of two strains: a minor subtype A strain very similar to the Eastern European subtype A strain (approximately 8% of cases), and a second major strain (77%) found to be most closely related to the CRF06-cpx strain, previously described only from African countries. The variability in the two clusters was very low, suggesting point source introductions. Ten percent of cases seemed to be newly generated recombinants of the A and CRF06-cpx strains. Analysis of viral diversification over time revealed a rate of change within the V3 region of 0.83%/year for the CRF06-cpx strain, consistent with findings from other subtypes. Due to the relatively frequently found novel recombinant forms, the Estonian HIV-1 epidemic may allow studies of coinfection and intersubtype recombination in detail.

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.

Rapid epidemic spread of HIV type 1 subtype A1 among intravenous drug users in Latvia and slower spread of subtype B among other risk groups

To investigate the rapid HIV epidemic in Latvia, 97 newly detected individuals were sampled in 2000-2001. To establish the molecular epidemiology we sequenced the env V3 and gag p17 regions of the HIV genome and compared them with reference sequences using phylogenetic analyses. As expected, the vast majority (n = 88; 91%) were intravenous drug users (IDUs) from the Riga region. Also, the majority of the investigated individuals (n = 93; 96%) were found to carry a subtype A1 virus that may have entered the Latvian IDU population several times. In addition, one IDU was infected with CRF03_AB and three other individuals, who had been infected through sexual contacts, carried subtype B virus. Thus, subtype A1 dominates the Latvian epidemic and is strongly associated with the IDU risk group. Although some spread of subtype A1 has occurred in the heterosexual group, subtype B dominates among homosexually and heterosexually infected individuals.

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.

First report of CRF03_AB recombinant HIV type 1 in injecting drug users in Ukraine

It has been reported that subtypes A throughout H of HIV-1 are circulating in the former Soviet Union. In this sequence note, we analyzed the genetic prevalence of HIV-1 among injecting drug users (IDUs) in Ukraine. The subjects studied included two individuals from Kiev and six individuals from Simferopol', the latter located in the Crimean Peninsula. We found that one of eight IDUs was infected with a CRF03 gagA/envB recombinant HIV-1 and was from Simferopol', whereas the others were infected with HIV-1 subtype A. There combinant was closely related to other A/B recombinants reported previously, and had silent mutations Inthe V3 region, the same as other envB strains of HIV-1 circulating among mDUs in the former Soviet Union. The data supported reports that the Russian AIB recombinant HIV-1 was probably from Ukraine. This is the first report of a CRF03 gagA/envB recombinant HIV-1 found in Ukraine.

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.

HIV type 1 genetic diversity and genotypic drug susceptibility in the Republic of Moldova

HIV-1 genetic diversity and, for the first time, genotypic drug susceptibility was investigated for strains circulating in the Republic of Moldova (of the former Soviet Union). Eighty-three samples from adults recently infected by intravenous drug use (IDU) (n = 60), heterosexual contact (n = 8), and from blood donors (n = 15) that tested positive from 1997 to 1998, and originating from different regions of Moldova were serotyped. By group-specific and subtype-specific peptide ELISA, patients were infected by serotype A (n = 65), serotype B (n = 1), or were nontypable (n = 17). Heteroduplex mobility assay (HMA) confirmed 11 subtype A and the one subtype B infection. Analyses of pol and env sequences for six of the IDUs confirmed that they were infected with subtype A strain. These strains clustered tightly with subtype A strains isolated from the former Soviet Union in phylogenetic analysis. No mutations associated with drug resistance were detected. The Republic of Moldova is culturally more closely related to Romania (where subtype F dominates the epidemic), but depends economically on Russia (where subtype A is established among IDUs). Thus, our results suggest that the spread of HIV in this region is driven by drug networks rather than being due to cultural similarities.

Genetic subtypes of HIV type 1 viruses circulating in the Czech Republic

We analyzed the genetic diversity of HIV-1 strains circulating in the Czech Republic. Phylogenetic analysis of the env and gag gene sequence fragments from 39 isolates revealed that the majority of these strains (32 of 39, 82%) were of subtype B; other genetic subtypes identified were A, C, F, and recombinant circulating form CRF01_AE. The isolates that did not cluster with subtype B originated almost exclusively from a heterosexual route of transmission. The molecular epidemiological data are suggestive of multiple entry of HIV-1 infection into the Czech Republic and show that the genetic pattern of the HIV-1 strains circulating in this country corresponds to that found in other European countries.

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.

The molecular epidemiology of HIV type 1 of men in Mexico

Genotypic characteristics of human immunodeficiency virus type 1 (HIV-1) in Mexico were investigated in a multicenter study that involved centers in five geographic regions of the country. Study samples (n = 65) collected from male patients in 1998-1999 were sequenced within the C2-V5 region of the gp120 env gene. Phylogenetic analysis revealed that subtype B predominates in Mexico. The level of interpatient nucleotide diversity (mean value of 8.9%) was congruent with multiple introductions of the virus and the "aging" epidemic in Mexico. One-third of samples (30.8% of cases) showed polymorphism within the crown of the V3 loop demonstrating non-GPGR motifs. Two new motifs in the V3 loop crown - HPGG and GPEG - were observed. The evolution of the AIDS epidemic in Mexico should be closely monitored since non-B HIV-1 subtypes might be introduced. The nucleotide sequences were deposited in the GenBank under accession numbers AF200855-AF200869, AF200871-AF200892, and AF200894-AF200921.

Serological approaches to subtyping of HIV-1 in injecting drug users in Russia: evidence of subtype homogeneity at the main sites of the epidemic

The aim of this study was to develop and evaluate a simple V3 peptide-based enzyme immunoassay (PEIA) for large-scale serotyping of HIV-1 specimens derived from injecting drug users (IDUs) in the Russian Federation. Two synthetic peptides were evaluated, named P1 (RKSIHIGPGRAFYATGD) and P2 (RTSVRIGPGQVFYKTGD), in an PEIA on 63 HIV-1 IDUs sera for which genotypes had been determined by heteroduplex mobility assay (HMA) and sequencing. The sensitivities of P1 (subtype B) and P2 (subtype A) were 87% and 75% respectively. Specificity of the assay was 100% for both peptides, with 100% predictive values of a monoreactive positive test for both peptides. Using the PEIA with peptides P1 and P2, we have serotyped 375 of 477 serum samples derived from IDUs in 4 main sites of the HIV-1 epidemic in Russia. The results demonstrated a high level of subtype homogeneity in all regions studied. In 3 of 4 territories, Tver' (n=345), and Rostov-on-Don (n=61) regions, and Krasnodar Kray (n=27), 100% of typable sera were found to belong to env subtype A. On the other hand, all specimens serotyped in the Kaliningrad region (n=38) belonged to env subtype B, and there is strong evidence that the recombinant gagAenvB virus which has caused the largest outbreak of HIV-1 in Russia is located in this region. At the present time another parental strain with gagBenvB genotype is of minor importance in the IDUs HIV-1 epidemic in Russia.

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.

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. UNAIDS Virus Isolation Network

In the former Soviet Union (SU) increasing numbers of HIV-1 infections among injecting drug users (IDU) have been reported, especially in the Ukraine. The main subtype transmitted among the IDUs seems to be subtype A, but limited numbers of subtype B cases have also been reported. In Kaliningrad, Russia, an AB recombinant strain was earlier shown to be responsible for the local outbreak. Here we describe the genetic relationship of HIV-1 strains circulating among IDUs in the former SU. For subtype A and the AB recombinant strains nearly full-length genomes were sequenced, and for one subtype B strain the entire envelope gene was cloned. The relationship between the AB recombinant strain and the subtype A and subtype B strains and the mosaic structure of the recombinant was studied by phylogenetic analysis. Ukrainian A and B strains were shown to be the probable parental viruses of the Kaliningrad AB recombinant strain. In the envelope gene the recombination breakpoint could also be precisely mapped to a region of similarity of only 14 base pairs. This suggests that only short stretches of absolute sequence identity may be needed for efficient RNA recombination between HIV-1 subtypes.

Silent mutation in the V3 region characteristic of HIV type 1 env subtype B strains from injecting drug users in the former Soviet Union

New independent states of the former Soviet Union are facing a rapidly growing epidemic of HIV-1 among injecting drug users (IDUs). This epidemic is caused by three HIV-1 populations, one belonging to HIV-1 subtype A (IDU-A), another to subtype B (IDU-B), and the third being a recombinant of the IDU-A and IDU-B viruses (IDU-A/B, gagA/envB). Each of these populations is characterized by a high level of genetic homogeneity. We identified a unique synonymous nucleotide substitution in the first isoleucine codon at the IHIGPGR motif (ATT), which was observed in the env subtype B V3 sequences derived from IDUs in Russia and the Ukraine. This substitution was observed in none of 179 sequences obtained from IDUs in western Europe, northern America, and Asia. Molecular epidemiological analysis of HIV-1 strains based on this sequence pattern could be useful for tracing the origin and spread of the IDU-B viruses to other countries and risk groups.

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

Countries of the former Soviet Union are experiencing an emerging HIV-1 epidemic due to a rapid expansion of HIV-1 among injecting drug users (IDUs). To study the molecular epidemiology of HIV-1 among IDUs in St. Petersburg, Russia, virus sequences were obtained from 22 individuals. Phylogenetic analysis of the env and gag regions revealed circulation of two major HIV-1 populations, one belonging to HIV-1 subtype A, and another being a recombinant of subtype A and B viruses (gagA/envB). Both virus populations were highly homogeneous, with a mean pairwise genetic distance of <2%, and similar to viruses obtained earlier from IDUs in other regions of the former Soviet Union. Distribution of the two major HIV-1 genotypes in St. Petersburg correlated with geographical origin of infections. In one individual, a virus type previously unseen among IDUs was found, which demonstrates the possibility that new viruses are entering this risk group.