Next Article in Journal
Changes in Body Mass Index after Initiation of Antiretroviral Treatment: Differences by Class of Core Drug
Next Article in Special Issue
Prevalence and Structure of HIV-1 Drug Resistance to Antiretrovirals in the Volga Federal District in 2008–2019
Previous Article in Journal
Evaluation of EPISEQ SARS-CoV-2 and a Fully Integrated Application to Identify SARS-CoV-2 Variants from Several Next-Generation Sequencing Approaches
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 14
Issue 8
10.3390/v14081675
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Aleksey Lebedev
1,*,
Anna Kuznetsova
1,
Kristina Kim
1,
Ekaterina Ozhmegova
1,
Anastasiia Antonova
1,
Elena Kazennova
1,
Aleksandr Tumanov
1,
Adkhamjon Mamatkulov
2,
Evgeniya Kazakova
2,
Nargiz Ibadullaeva
2,
Krestina Brigida
2,
Erkin Musabaev
2,
Dildora Mustafaeva
3,
Visola Rakhimova
4 and
Marina Bobkova
1
1
Gamaleya National Research Center for Epidemiology and Microbiology, 123098 Moscow, Russia
2
Research Institute of Virology, Tashkent 100194, Uzbekistan
3
Republican AIDS Center, The Ministry of Health, Tashkent 100135, Uzbekistan
4
Center for Development of Profession Qualification of Medical Workers, Tashkent 100007, Uzbekistan
*
Author to whom correspondence should be addressed.
Viruses 2022, 14(8), 1675; https://doi.org/10.3390/v14081675
Submission received: 17 June 2022 / Revised: 25 July 2022 / Accepted: 27 July 2022 / Published: 29 July 2022
(This article belongs to the Special Issue Molecular Epidemiology of HIV Infection in Eastern Europe and Central Asia)
Browse Figures
Review Reports Versions Notes

Abstract

:
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.

1. Introduction

Thirty years ago, in 1991, the Soviet Union broke up and its 15 republics gained independence. After the initial socio-economic turmoil, most of them started to recover, but they are still facing the enormous challenge of controlling the rapid spread of HIV infection, which reached epidemic proportions. Approximately 2.2 million people in Eastern Europe and Central Asia were living with HIV in 2020; about 58,000 of these were in the Republic of Uzbekistan [1,2]. During 2020, 3983 new HIV infections were diagnosed in the Republic. Of these new infections, 73.2% were associated with sexual transmission, which is considered to be the most important factor driving HIV spread in Uzbekistan [3]. Nevertheless, around 18.2% of new infections occur in injecting drug users (IDUs) [3], with HIV prevalence in this group reaching nearly 30.0% [4,5]. In Uzbekistan and many other Central Asian countries, sexual risks—primarily heterosexual (HET)—become conflated with risks associated with injection drug usage [6,7,8]. At least 20% of female sex workers (FSW) are believed to inject drugs, which increase the risk of infection by 5–6 times [6,7]. The prevalence of HIV infection is high among Uzbek migrants (4.1%) [4], who are regarded as a driving force of the HIV-1 spread in former Soviet Union (FSU) countries [9,10,11]. Given all this, the identification of transmission networks and clusters affecting key populations is of paramount importance for preventing further HIV transmission.
The HIV-1 subtype patterns in Uzbekistan are complex and continuously changing in response to human population migration and active transmission networks. In the Republic, as many as six known HIV-1 subtypes and circulating recombinant forms (CRFs) have been detected, with CRF02_AG (57.0%) and sub-subtype A6 (41.0%) accounting for 98% of total infections [12]. The latest National HIV Molecular Epidemiological Survey performed by our group revealed that CRF02_AG, initially identified in IDUs in the capital [13], has quickly overtaken subtype A6 in this vulnerable population during the past few years. Furthermore, CRF02_AG also became dominant among HETs, increasing from 6.0% in 2003 to 53.0% in 2016 [12,14]. According to our estimates, CRF02_AG today can account for 60–65% of all HIV-infections in Uzbekistan and is likely to increase further. For this reason, it is imperative to understand the trends of CRF02_AG and A6 epidemics in Uzbekistan in order to inform targeted public health interventions.
The aim of our current study was to examine phylogenetic clustering, cluster composition and factors associated with clustering in a cohort of ‘Uzbek’ HIV+ people infected with two dominating HIV-1 subtypes: CRF02_AG and A6.

2. Materials and Methods

2.1. Sequence Dataset Compilation

HIV-1 sequences included in this study were derived mostly (95.6%) from participants in the Nationwide Surveillance of HIV Drug Resistance (NS-HIVDR) program [12], which constitutes a well-defined and demographically diverse cohort of antiretroviral therapy (ART)-naive and ART-experienced HIV-positive people from 13 different Uzbekistan regions. A total of 316 HIV-1 sequences covering the complete protease (PR) and a part of the reverse transcriptase (RT) of pol (polymerase) gene (nucleotides 2253 to 3344 of reference strain HXB2; 1092 bp long) were used, of which 59.8% (188/316) were sampled from participants infected with HIV-1 CRF02_AG, and 40.2% (127/316) from patients infected with HIV-1 subtype A6.
The blood plasma samples were collected during 2010–2016 at the Research Institute of Virology of the Ministry of Health of the Republic of Uzbekistan, or at local hospitals in different regions of the Republic; in addition, 37 sequences (unpublished studies) were retrieved from the Los Alamos HIV Sequence Database (https://www.hiv.lanl.gov/content/index, accessed on 20 May 2021) (Figure 1; Table S1). The epidemiological metadata, including transmission route, age, sex, country of infection, sampling area, immunological (CD4+ T-cell count) and virological (viral load) information, were available for 210 of 316 sequences and used as characteristic data. A total of 739 CRF02_AG/CRF063 (N = 346) and subtype A6 (N = 393) pol sequences most closely related to each Uzbekistan sequence from the BLAST search at GenBank (with >95% sequence identity and de-duplicated) were also used; in particular, the CRF063_02A6 sequences were used to improve the clock-like structure in datasets, but were not used to reconstruct migration pathways. The GenBank accession numbers for sequences used in this study are provided as Supplementary Information (Table S2).

2.2. Sequence Alignment and Subtype Assignment

Nucleotide sequences were aligned using MAFFT online service (https://mafft.cbrc.jp/alignment/server/, accessed on 20 May 2021). HIV-1 subtype analysis was performed by our group earlier and described in the article by Mamatkulov et al. [12]. The RIP tools (https://www.hiv.lanl.gov/content/sequence/RIP/RIP.html, accessed on 20 May 2021) were also used for recombination analysis.

2.3. Identification of Transmission Clusters and Uzbek-Origin Clusters

The clusters of closely related sequences (transmission clusters) and potential origin of the HIV-1 clades were evaluated by phylogenetic approach, as described in the study by Alexiev et al. [15]. Briefly, we used the MicrobeTrace tool [16] to cluster selections in ML trees using the sequence alignment, including all 316 ‘Uzbek’ pol sequences, at Tamura–Nei pairwise genetic distance (d) threshold of 0.015 substitutions/site and having at least 3 sequences; the ML trees in Newick format were importing into MicrobeTrace after the reconstruction by Nextstrain tool (https://nextstrain.org/, accessed on 20 May 2021) under GTR (general time reversible) model. Visualizing clusters by node and total number of links in the transmission network was performed in MicrobeTrace. The degree of centrality, measured as the number of direct connections a node has to other nodes (that is, how many direct “one hop” connections each node has to other nodes in the network), was manually estimated based on the formula contained in [17]. We then sought epidemiological predictors of clustering and calculated assortativity for epidemiological variables such as region, sex, and transmission category.
To evaluate the potential origins of HIV-1 clusters in Uzbekistan, all sequences were combined with those from other countries and analyzed with the Nextstrain tool (https://nextstrain.org/) that uses ML analysis implemented in Tree-Time [18,19]. The total dataset for Nextstrain ML analysis contained 1055 HIV-1 sequences, comprising 316 sequences from Uzbekistan and 739 global HIV-1 pol sequences, plus one HIV-1 subtype C sequence (GenBank accession numbers: AF443091) used as the outgroup. The output JASON files from Nextstrain pipeline were viewed using Auspice (https://auspice.us, accessed on 20 May 2021).

2.4. Statistical Analysis

Epidemiological data, including sex, age, country of birth, region of Uzbekistan, and transmission categories, were analyzed by Pearson’s chi-squared or two-tailed Fisher’s exact test (where necessary). Multivariate logistic regression analysis was performed with the outcome variable as HIV-1 subtype: CRF02_AG or A6. Odds ratios were estimated with 95% confidence intervals (CIs). The differences were considered significant at a p-value less than 0.05. All analyses were performed in STATISTICA v.10.0 software (StatSoft, Tulsa, OK, USA). Assortativity coefficients for selected variables were calculated using the Python v3.9.5 (https://www.python.org/downloads/release/python-395/, accessed on 20 May 2021) script by Sergey-Knyazev (https://github.com/Sergey-Knyazev/attribute_assortativity, accessed on 20 May 2021) and NetworkX package (https://networkx.org, accessed on 20 May 2021), with thresholds of d = 0.015 substitutions/site (1.5%).

3. Results

3.1. Association of Demographic Characteristics with HIV-1 Infection in Uzbekistan

To perform multivariable analysis, we used a subset of 210 HIV-positive people for whom there was information about all the variables considered: sex, risk group, region of residence, and age. Of these, 124 (59.0%) sequences belonged to CRF02_AG and 86 (41.0%) sequences belonged to sub-subtype A6 HIV-1. It was also decided to combine the sequences from the north-eastern part of the Republic (Andijan, Fergana, Namangan and Tashkent regions) into a single group termed ‘North Eastern Region’ (NER). The multivariable analysis found an association between the region of residence and HIV-1 subtype (Table 1). So, HIV-positive people from NER were more likely to have the CRF02_AG variant compared to those infected in other regions (OR = 2.62, 95% CI = 1.44–4.75, p = 0.001). Moreover, this analysis produced the same results as when using all sequences from Uzbekistan (N = 316); CRF02_AG was still the most prevalent HIV variant in NER, although the proportion (137/316, 43.4%) was higher than in the subset of 210 people (80/210, 38.1%) (Table S1). There was no evidence of association of other characteristics with HIV-1 subtypes.

3.2. HIV-1 Subtype Transmission Clusters in Uzbekistan

The phylogeny of HIV-1 CRF02_AG (N = 189) and A6 strains (N = 127), followed by processing and analysis in MicrobeTrace, indicated the presence of four transmission clusters, including two large clusters containing over 35 members, one for each subtype (Figure 2; Table 2). In total, 28.8% (91/316) of people from the dataset were members of transmission clusters. The largest cluster (CRF02_AG; N = 49) mainly comprised women (67.3%) and IDUs (38.8%), with a majority of people living in Andijan (N = 25, 51.0%), Tashkent (N = 7, 14.3%), Fergana (N = 4, 8.2%), and Samarkand (N = 8, 16.3%) regions; that is, about 90% of cluster members lived in the NER. By contrast, the second largest cluster (sub-subtype A6; N = 36) consisted mostly of men (52.7%) and HETs (44.0%); about two-thirds of them were from Tashkent (N = 22, 61.1%), while the rest were mainly from Andijan (N = 4, 11.1%) and the Samarkand region (N = 3, 8.3%), i.e., more than 72% were from the NER. The members from remaining clusters (N = 3, in each cluster) were mostly HETs from a single home region (such as Xorazm). It is interesting to note that women had a significant presence within networks identified in both the CRF02_AG and A6 clusters. With regard to the characteristics of people involved in the clusters, all were natives of Uzbekistan; the overall proportion of labor migrants among people in our study could be as much as 30% (personal data).
We also calculated assortativity among HIV-infected people in major clusters and dyads by sex, region, and transmission category using assortativity (r) coefficients (Table 3). The r-value characterizes the degree of mixing (i.e., disassortativity) or non-mixing (i.e., assortativity) of characteristics between entities that are connected within a network. The r-value ranges from −1 to 1, where negative values represent dissasortativity or mixing between entities of different characteristics, and positive values represent assortativity or non-mixing; zero value represents random mixing with no tendency for assortativity [20]. Our analysis found that region was the only characteristic with high assortativity (r > 0.25), being variable among HIV subtypes and highest among CRF02_AG (r = 0.30; total, node size = 67) (Table 3). The set was modestly dissasortative (r > 0.25) by HIV-transmission risk for HIV-1 subtype A6-infected people.

3.3. Origin and Cross(in)-Country Transmission of HIV-1 Sub-Subtype A6 and CRF02_AG in Uzbekistan

Here, we combined NS-HIVDR subtype A6 and CRF02_AG sequences with similar sequences of the same subtypes identified by BLAST search in HIV Databases at LANL (Los-Alamos National Laboratory), followed by data analysis in Nextstrain. A visual inspection of our phylogenetic trees suggested a certain degree of genetic structure based on geographical location (Figure 3). Analysis of CRF02_AG viruses shows that the majority of transmission events took place in Central Asia and formed a large older clade, with ancestor root location most probably traced to Uzbekistan. The CRF02_AG strains from Uzbekistan, Kyrgyzstan, Tajikistan and Kazakhstan, and a small number from Russia, are intermingled within this clade to form several minor country-specific subclades; overall, apart from one subclade formed by the Russian CRF063_026A sequences (marked with an asterisk), no other country showed a clearly compartmentalized grouping of sequences (Figure 3A). The only more or less large ‘Uzbek’ subclade (scl-1, N = 94) consisted of 30.2% (59/189) of ‘Uzbek’ sequences from people living mostly in the Andijan (50.9%, 29/57) and Tashkent (24.6%, 14/57) regions and included all sequences of a large transmission cluster, described above; the overwhelming majority of them were related to IDUs (36.8%, 21/57). The minor Uzbek subclade 2 (scl-2, N = 36) and subclade 3 (scl-3, N = 27) included of 77.8% (28/36) and 66.6% (18/27) Uzbek sequences, respectively; most peoples in both subclade (≥50%) were from Tashkent and with a comparable proportion IDUs and HET (~30%).
The finding that almost all ‘Uzbek’ CRF02_AG sequences could be traced back to a unique most recent common ancestor (MRCA) suggests a single major introduction event, followed by its local spread. Using an evolutionary time scale, the Nextstrain analysis indicated that the CRF02_AG outbreak cluster in Uzbekistan evolved from a MRCA introduced around 25-05-2000 [12-02-1998, 14-04-2001]. The MRCA for one large and two minor ‘Uzbek’ subclades corresponded to 02-03-2003 [28-04-2002, 12-12-2004], 29-08-2001 [27-01-2001, 12-06-2002] and 27-01-2003 [10-03-2002, 21-08-2004], respectively.
The reconstruction of migration pathways between different geographic areas showed that the early dissemination of CRF02_AG most likely occurred in Uzbekistan and Kazakhstan, and thereafter the epidemic spread to Kyrgyzstan and Russia (Figure 3C). We also found that the most significant migration routes of HIV-1 CRF02_AG showed directionality from Uzbekistan to these countries and only very few transmissions went in the opposite direction. Regarding the intra-country migrations events, there were 15 HIV transmission pathways found, of which 11 involved the city of Tashkent (Figure 4A).
Unlike CRF02_AG, the majority of sub-subtype A6 transmission events took place in East European countries (Belarus, Russia and Ukraine), followed by some regional isolation after local introduction (as seen for Uzbekistan) and spread and mixing in other geographic locations (as for infections in Russia) (Figure 3B). We revealed at least three clades consisting almost exclusively of ‘Uzbek’ sequences. The large clade (cl-1, N = 75) contained about half of all ‘Uzbek’ sequences (48.0%, 61/127) and included all sequences of the large transmission cluster described above. The majority of cases (60.7%, 37/61) in this clade were diagnosed in the capital city of Tashkent, with most of them 47.5% (29/61) being IDUs, followed by HETs (32.8%, 20/61). The medium-sized clade (cl-2, N = 14) contained 78.6% (11/14) Uzbek HIV cases from Samarkand (81.8%, 9/11) with an unknown risk of infection (54.5%, 6/11). The smallest clade (cl-3, N = 10) comprised 90% Uzbek HIV-positive people, mainly HETs (55.6%, 5/9) from six different regions.
The most probable geographic origin of the large and smallest clades was Russia; there were ten A6 sequences from Ukraine near the root of medium-sized clade. Maximum-likelihood phylodynamic analysis inferred the MRCA [95% HPD] for large, medium-sized and smallest subclades at 08-02-2000 [11-02-1999, 31-07-2001], 09-08-2000 [29-12-1998, 29-01-2002] and 09-08-2000 [29-11-1998, 29-07-2002], respectively. The MRCA for the complete clade of 519 subtype A6 taxa corresponded to 03-12-1994 [95% HPD: 03-11-1991, 07-03-1996].
The analysis of migration pathways between different geographic areas supported at least five HIV-1 sub-subtype A6 migration events into Uzbekistan (Figure 3D). As expected, the original HIV-1 A6 source for at least some of the infections in Uzbekistan probably originated from Ukraine; the same was true for Russia and Belarus, which were affected by strong gene flow from Ukraine. However, the most significant migrations routes of HIV-1 subtype A6 showed directionality from Russia to Uzbekistan and back. We also found HIV transmissions from Uzbekistan to Tajikistan, Kazakhstan and Kyrgyzstan. Intra-country migration events (N = 12) mainly related to the city of Tashkent (N = 7) as in the previous case, and also Fergana/Samarkand (N = 3) (Figure 4B).

4. Discussion

The HIV epidemic has been going on for over 20 years in Uzbekistan, and transmission patterns have gradually evolved from injecting drug usage to heterosexual contacts [3,12,13]. Unlike FSU countries with a single HIV subtype (e.g., Caucasian republics with subtype A6 at more than 85.0% [21,22,23]), today there are two principal HIV-1 subtypes, CRF02_AG and A6 in about equal proportions, in Uzbekistan [12]. Our current data suggest that the majority of CRF02_AG infected people in Uzbekistan (43.4%) live in the NER, where the possibility of being infected with this HIV-1 variant and not A6 sub-subtype is twice as high as subject from other regions. The earlier dispersal of CRF02_AG to Tashkent and regions bordering the capital (Andijan, Fergana, Namangan and Syrdaryo), together with tight epidemiological links between these locations due to their geographical proximity and high level of people migration, as well as possible introduction of CRF02_AG into highly-connected transmission networks, appear to be the most likely causes of this phenomenon. In this study, we investigated the transmission networks of the most prevalent HIV subtypes by association of all available sequences with the demographic and epidemiologic data.
To date, there is no consensus on methodology to infer HIV-1 molecular networks, but it is obvious that threshold selection is a key factor in network construction [24,25]. Here, we used single optimal thresholds (1.5%) for both subtypes. The network analysis showed that just over a third of HIV+ individuals were organized into clusters (including large clusters of over 35 members), the majority of which were clusters that contained non-dyads. Among people involved in clusters, all were natives of Uzbekistan, with a majority (over 70%) living in the NER (Andijan, Fergana, Namangan, and Tashkent regions). We found that women were slightly more likely to participate in transmission clusters, a finding that aligned with the current reports describing an increase in infections among women in Uzbekistan [3,12]. Moreover, these data highlight that females may be engaged in riskier sexual and drug-related behavior and they could be better supported by the expansion of molecular-based surveillance strategies to reduce transmission.
It was previously established that 2–3 out of every 1000 women in Tashkent are involved in sex work [26]. In addition to this, about 20–24% of male IDUs have contacts with FSWs [27]. Female IDUs also exist and are highly susceptible to HIV infection through exchanging sex for drugs or money and risky injecting drug use, which increases the overlap between these risk groups [6,27,28]. All of this may explain the patterns of group risk mixing that we have identified in transmission clusters. Data in this report also indicate that IDU/HET-women, having high centrality within the highly connected network, can contribute more CRF02_AG infections than IDU-men, and a strategy to disrupt these contacts would be more effective in reducing infections than randomly targeted prevention approaches.
Although the high percentage of non-clustered sequences might be explained by low number of sequences included in the study or low sequence completeness (lack of sequences from HIV-infected persons with closely related viruses) [29], dead-end transmission events (non-transmitted infections) or infection abroad can also be seen as possible causes. The Uzbek labor migrants (or any other immigrants) are at an increased risk of HIV infection, since they often stay away from their hometowns for long time periods, and instead of having a monogamous sexual partner have multiple sexual contacts with FSWs. Although the high assortativity by regions within the country mostly supports local genetic connections, and the local populations are more likely to be infected inside the country (e.g., through regional networks), the migrant population needs urgent attention because HIV prevalence appears to be growing in this group.
Migration is one of many social factors contributing to the HIV/AIDS epidemic; the male/female migrants become infected while away from home and infect their spouses or regular partners upon return [30,31,32]. Uzbek migrants are mostly males going to Russia and Kazakhstan for work. They tend to have limited HIV knowledge and engage in sexual contacts involving sporadic use of condoms [33]. Our study supports this to be the case, but there is additional evidence for the bi-directionality of HIV transmission (between the country of birth to the country of residence): the results of phylogenetic study of CRF02_AG and A6 sub-subtypes in Uzbekistan are indirect indications (but not evidence) of this, in addition to the patient interviews.
It has been previously shown that after leaving Africa, CRF02_AG initially extended to IDUs and grew rapidly in Uzbekistan, followed by dissemination across the FSU countries, including Kazakhstan and Russia [34,35,36]. According to the latest data, this expansion originated from a single viral clade that arose in 1996, and the subsequent epidemic growth occurred approximately during 1998–2003 [34,37]. Our results are in full accordance with these findings, since we estimated that the time to MRCA (TMRCA) of this and the other FSU countries’ epidemics to be around 1996, with the epicenter in Tashkent.
In the late 1990s/early 2000s, Central Asian FSU countries were experiencing an explosive HIV-infection epidemic, which affected IDUs and their sexual partners [38,39,40,41,42]. To a large extent this was due to the location of the Central Asian drug trafficking routes (‘northern routes’) from Afghanistan towards their primary markets in Russia and Europe. As a result, 524 new HIV cases were diagnosed in Uzbekistan alone between 1999 and 2000 [14]. As with other IDU-associated HIV-infection epidemics in FSU during this period [43,44,45,46], the main mode of HIV transmission was syringe- and needle-sharing. These historical records support our phylogenetic findings. According to our data, CRF02_AG, after some regional isolation in Uzbekistan, as evidenced by several all-‘Uzbek’ clades, was introduced to the Central Asian countries and then Russia during a relatively short period via the IDU network; this may also explain the ‘similar’ HIV-1 epidemic on both sides of the Uzbekistan-Kazakhstan border (or Kyrgyzstan, and or Tajikistan [47,48]). Interestingly, we found only marginal cross-border CRF02_AG transmissions from FSU countries back to Uzbekistan, even in the case of Russia. Being a major hub for the Uzbek migrant populations, Russia is thus acting as a conductor zone, providing HIV transmission between the other FSU countries’. Even though in Russia other AG recombinants (CRF063) are widespread, CRF02_AG today only accounts for 3–5% of all HIV-infections in the country [49,50,51,52].
As part of the A6 epidemic, which was first described in the IDUs’ population, but spread among HET and even MSM (men who have sex with men) populations soon afterwards [7,8,9,10], we inferred at least two large ‘Uzbek’ A6 clades dated to the year 1999. The ancestral origin for these clades was traced with high probability to Russia and Ukraine, which supports previous findings [53] and reflects the primary role of these countries as the center of HIV epidemic spread to Uzbekistan during the late 1990s/early 2000s. Moreover, our findings also suggest the possibility of multiple introductions of A6 into Uzbekistan from these countries. However, we also report the dissemination of sub-subtype A6 from Uzbekistan back to Russia and Ukraine (as well as exportation to Tajikistan), with the participation of migrants infected through regional FSW-networks. We cannot rule out the possibility that migrant men were infected from their wives, who acquired the infection from their parallel sexual partners [33]. Overall, labor HET-migrants today may be a main factor fueling the HIV-1 bidirectional transmission between FSU countries. All of this points to the need for HIV prevention and education campaigns among high-risk groups, including a range of harm-reduction services such as syringe exchange programs, condom distribution, HIV testing, and other outreach activities.
Some comments should be made concerning the limitations of this study. First, limited coverage of HIV testing and small sampling depth in our study may have restricted the ability to detect transmission clusters. Hence, at the population level, the infection network could be wider than the one studied, with unknown transmissions occurring. Second, route of transmission data is missing for a large proportion of samples because of stigma and discrimination on the basis of real or perceived participation in key risk groups (IDUs, sexual workers or MSM) [54,55,56,57]; this missing data, to some extent, makes it difficult to understand the characteristics of HIV transmission clusters at a population level. Third, the small sample size in some countries is not fully representative of their general HIV-infected population. When combined with the large sequence number variability in each country, this may lead to bias in assessment of migration routes and detection of viral movements. Lastly, our results are limited to the pol region of HIV-1; the use of complete genomes may provide more accurate inference of transmission histories. Thus, some caution is advised regarding the interpretation of the findings from the present study.

5. Conclusions

In summary, this study based on genetic data identifies the transmission networks of two major HIV-1 subtypes in Uzbekistan, emphasizes that IDU-related and sexual exposures are involved as mechanisms of transmission in the overlapping social networks, and indicates the complexity and high burden of HIV-1 infection in the population. We also demonstrate that the existing migration flow probably contributes to an increase in the HIV infected population, and that cross-border transmissions among high-risk groups occur bi-directionally between Uzbekistan and other FSU populations. These results are of public health importance and suggest that preventive action should be reinforced in this area to improve the management of high-risk practices.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14081675/s1, Table S1: General characteristics of HIV-1 infected patients from Uzbekistan stratified by subtypes; Table S2: GenBank Accession numbers, origin information and sampling date for the HIV-1 sequences used in this study.

Author Contributions

Conceptualization, A.L. and M.B.; methodology, A.L., A.K. and K.K.; validation, A.L., E.O. and A.A.; formal analysis, A.L., A.K. and K.K.; investigation, A.L., A.A., A.T., E.K. (Elena Kazennova); resources, E.K. (Evgeniya Kazakova), N.I., K.B., D.M., V.R. and E.M.; data curation, A.L., E.K. (Evgeniya Kazakova) and A.M.; writing-original draft preparation, A.L.; writing-review and editing, A.L., E.K. (Evgeniya Kazakova) and M.B.; visualization, A.L.; supervision, E.K. (Evgeniya Kazakova), A.M., A.K. and M.B.; project administration, M.B.; funding acquisition, M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grant from the Russian Science Foundation, project No. 22-15-00117 (https://rscf.ru/project/22-15-00117, accessed on 28 July 2022).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by Research Ethics Committee of the Ministry of Health of the Republic of Uzbekistan (№ 1/1-640, 2015-01-22).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Bobkova, M. Current status of HIV-1 diversity and drug resistance monitoring in the former USSR. AIDS Rev. 2013, 15, 204–212. [Google Scholar] [PubMed]
  2. Joint United Nations Programme on HIV and AIDS (UNAIDS). UNAIDS: Region Factsheets Eastern Europe and Central Asia 2020. Available online: https://www.unaids.org/en/regionscountries/easterneuropeandcentralasia (accessed on 20 May 2021).
  3. Joint United Nations Programme on HIV and AIDS (UNAIDS). UNAIDS: Uzbekistan 2020 Progress Country Report. Available online: https://www.unaids.org/sites/default/files/country/documents/UZB_2020_countryreport.pdf (accessed on 20 May 2021).
  4. Mathers, B.M.; Degenhardt, L.; Phillips, B.; Wiessing, L.; Hickman, M.; Strathdee, S.A.; Wodak, A.; Panda, S.; Tyndall, M.; Toufik, A.; et al. Global epidemiology of injecting drug use and HIV among people who inject drugs: A systematic review. Lancet 2008, 372, 1733–1745. [Google Scholar] [CrossRef]
  5. Davlidova, S.; Haley-Johnson, Z.; Nyhan, K.; Farooq, A.; Vermund, S.H.; Ali, S. Prevalence of HIV, HCV and HBV in Central Asia and the Caucasus: A systematic review. Int. J. Infect. Dis. 2021, 104, 510–525. [Google Scholar] [CrossRef] [PubMed]
  6. Todd, C.S.; Khakimov, M.M.; Alibayeva, G.; Abdullaeva, M.; Giyasova, G.M.; Saad, M.D.; Botros, B.A.; Bautista, C.T.; Sanchez, J.L.; Carr, J.K.; et al. Prevalence and correlates of human immunodeficiency virus infection among female sex workers in Tashkent, Uzbekistan. Sex. Transm. Dis. 2006, 33, 496–501. [Google Scholar] [CrossRef] [PubMed]
  7. Mamaev, T. Results of HIV sentinel epidemiological surveillance among sex workers in Osh City of Kyrgyz Republic. Zh. Mikrobiol. Epidemiol. Immunobiol. 2007, 3, 72–74. [Google Scholar]
  8. Beyrer, C.; Patel, Z.; Stachowiak, J.A.; Tishkova, F.K.; Stibich, M.A.; Eyzaguirre, L.M.; Carr, J.K.; Mogilnii, V.; Peryshkina, A.; Latypov, A.; et al. Characterization of the emerging HIV type 1 and HCV epidemics among injecting drug users in Dushanbe, Tajikistan. AIDS Res. Hum. Retrovir. 2009, 25, 853–860. [Google Scholar] [CrossRef] [PubMed]
  9. Amirkhanian, Y.A.; Kuznetsova, A.V.; Kelly, J.A.; Difranceisco, W.J.; Musatov, V.B.; Avsukevich, N.A.; Chaika, N.A.; McAuliffe, T.L. Male labor migrants in Russia: HIV risk behavior levels, contextual factors, and prevention needs. J. Immigr. Minor. Health 2011, 13, 919–928. [Google Scholar] [CrossRef] [Green Version]
  10. Weine, S.M.; Kashuba, A.B. Labor migration and HIV risk: A systematic review of the literature. AIDS Behav. 2012, 16, 1605–1621. [Google Scholar] [CrossRef] [PubMed]
  11. Aibekova, L.; Foley, B.; Hortelano, G.; Raees, M.; Abdraimov, S.; Toichuev, R.; Ali, S. Molecular epidemiology of HIV-1 subtype A in former Soviet Union countries. PLoS ONE 2018, 13, e0191891. [Google Scholar] [CrossRef]
  12. Mamatkulov, A.; Kazakova, E.; Ibadullaeva, N.; Joldasova, E.; Bayjanov, A.; Musabaev, E.; Kan, N.; Mustafaeva, D.; Lebedev, A.; Bobkova, M.; et al. Prevalence of Antiretroviral Drug Resistance Mutations Among Pretreatment and Antiretroviral Therapy-Failure HIV Patients in Uzbekistan. AIDS Res. Hum. Retrovir. 2021, 37, 38–43. [Google Scholar] [CrossRef] [PubMed]
  13. Kurbanov, F.; Kondo, M.; Tanaka, Y.; Zalalieva, M.; Giasova, G.; Shima, T.; Jounai, N.; Yuldasheva, N.; Ruzibakiev, R.; Mizokami, M.; et al. Human immunodeficiency virus in Uzbekistan: Epidemiological and genetic analyses. AIDS Res. Hum. Retrovir. 2003, 19, 731–738. [Google Scholar] [CrossRef]
  14. Carr, J.K.; Nadai, Y.; Eyzaguirre, L.; Saad, M.D.; Khakimov, M.M.; Yakubov, S.K.; Birx, D.L.; Graham, R.R.; Wolfe, N.D.; Earhart, K.C.; et al. Outbreak of a West African recombinant of HIV-1 in Tashkent, Uzbekistan. J. Acquir. Immune Defic. Syndr. 2005, 39, 570–575. [Google Scholar]
  15. Alexiev, I.; Campbell, E.M.; Knyazev, S.; Pan, Y.; Grigorova, L.; Dimitrova, R.; Partsuneva, A.; Gancheva, A.; Kostadinova, A.; Seguin-Devaux, C.; et al. Molecular Epidemiological Analysis of the Origin and Transmission Dynamics of the HIV-1 CRF01_AE Sub-Epidemic in Bulgaria. Viruses 2021, 13, 116. [Google Scholar] [CrossRef] [PubMed]
  16. Campbell, E.M.; Boyles, A.; Shankar, A.; Kim, J.; Knyazev, S.; Cintron, R.; Switzer, W.M. MicrobeTrace: Retooling molecular epidemiology for rapid public health response. PLoS Comput. Biol. 2021, 17, e1009300. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, J.; Luo, Y. Degree Centrality, Betweenness Centrality, and Closeness Centrality in Social Network. Adv. Intell. Syst. Res. 2017, 132, 300–303. [Google Scholar] [CrossRef] [Green Version]
  18. Hadfield, J.; Megill, C.; Bell, S.M.; Huddleston, J.; Potter, B.; Callender, C.; Sagulenko, P.; Bedford, T.; Neher, R.A. Nextstrain: Real-Time tracking of pathogen evolution. Bioinformatics 2018, 34, 4121–4123. [Google Scholar] [CrossRef]
  19. Sagulenko, P.; Puller, V.; Neher, R.A. TreeTime: Maximum-likelihood phylodynamic analysis. Virus Evol. 2018, 4, vex042. [Google Scholar] [CrossRef]
  20. Newman, M.E. Mixing patterns in networks. Phys. Rev. E. Stat. Nonlin. Soft. Matter Phys. 2003, 67, 026126. [Google Scholar] [CrossRef] [Green Version]
  21. Dvali, N.; Parker, M.M.; Chkhartishvili, N.; Sharvadze, L.; Gochitashvili, N.; Abutidze, A.; Karchava, M.; DeHovitz, J.A.; Tsertsvadze, T. Characterization of HIV-1 subtypes and drug resistance mutations among individuals infected with HIV in Georgia. J. Med. Virol. 2012, 84, 1002–1008. [Google Scholar] [CrossRef]
  22. Laga, V.; Vasilyev, A.; Lapovok, I.; Grigoryan, S.; Papoyan, A.; Glushchenko, N.; Kazennova, E.; Bobkova, M. HIV Type 1 Subtype A1 Dominates in Armenia. Cur. HIV Res. 2015, 13, 219–225. [Google Scholar] [CrossRef] [PubMed]
  23. Kirichenko, A.; Kireev, D.; Lopatukhin, A.; Murzakova, A.; Lapovok, I.; Saleeva, D.; Ladnaya, N.; Gadirova, A.; Ibrahimova, S.; Safarova, A.; et al. Prevalence of HIV-1 drug resistance in Eastern European and Central Asian countries. PLoS ONE 2022, 17, e0257731. [Google Scholar] [CrossRef] [PubMed]
  24. Hassan, A.S.; Pybus, O.G.; Sanders, E.J.; Albert, J.; Esbjornsson, J. Defining HIV-1 transmission clusters based on sequence data. AIDS 2017, 31, 1211–1222. [Google Scholar] [CrossRef]
  25. Liu, M.; Han, X.; Zhao, B.; An, M.; He, W.; Wang, Z.; Qiu, Y.; Ding, H.; Shang, H. Dynamics of HIV-1 Molecular Networks Reveal Effective Control of Large Transmission Clusters in an Area Affected by an Epidemic of Multiple HIV Subtypes. Front. Microbiol. 2020, 11, 604993. [Google Scholar] [CrossRef] [PubMed]
  26. United Nations Office on Drugs and Crime (UNODC). UNODC: Sex Work, HIV/AIDS and Human Rights in Central and Eastern Europe and Central Asia 2010. Available online: https://www.unodc.org/documents/hivaids/publications/CEEAndCAsiaharm_05_sex_work_east_eur_0408.pdf (accessed on 20 May 2021).
  27. Thorne, C.; Ferencic, N.; Malyuta, R.; Mimica, J.; Niemiec, T. Central Asia: Hotspot in the worldwide HIV epidemic. Lancet Infect. Dis. 2010, 10, 479–488. [Google Scholar] [CrossRef] [Green Version]
  28. Novitsky, V.; Moyo, S.; Lei, Q.; De Gruttola, V.; Essex, M. Impact of sampling density on the extent of HIV clustering. AIDS Res. Hum. Retrovir. 2014, 12, 1226–1235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Platt, L.; Rhodes, T.; Judd, A.; Koshkina, E.; Maksimova, S.; Latishevskaya, N.; Renton, A.; McDonald, T.; Parry, J.V. Effects of sex work on the prevalence of syphilis among injection drug users in 3 Russian cities. Am. J. Public Health 2007, 97, 478–485. [Google Scholar] [CrossRef]
  30. Decosas, J.; Adrien, A. Migration and HIV. AIDS 1997, 11 (Suppl. A), S77–S84. [Google Scholar] [PubMed]
  31. Mabey, D.; Mayaud, P. Sexually transmitted diseases in mobile populations. Genitourin. Med. 1997, 73, 18–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Lurie, M.N.; Williams, B.G.; Zuma, K.; Mkaya-Mwamburi, D.; Garnett, G.; Sturm, A.W.; Sweat, M.D.; Gittelsohn, J.; Abdool Karim, S.S. The impact of migration on HIV-1 transmission in South Africa: A study of migrant and nonmigrant men and their partners. Sex. Transm. Dis. 2003, 30, 149–156. [Google Scholar] [CrossRef] [Green Version]
  33. Smolak, A.; El-Bassel, N.; Malin, A.; Terlikbayeva, A.; Samatova, S. Sex workers, condoms, and mobility among men in Uzbekistan: Implications for HIV transmission. Int. J. STD AIDS 2016, 27, 268–272. [Google Scholar] [CrossRef] [PubMed]
  34. Mir, D.; Jung, M.; Delatorre, E.; Vidal, N.; Peeters, M.; Bello, G. Phylodynamics of the major HIV-1 CRF02_AG African lineages and its global dissemination. Infect. Genet. Evol. 2016, 46, 190–199. [Google Scholar] [CrossRef] [Green Version]
  35. Eyzaguirre, L.M.; Erasilova, I.B.; Nadai, Y.; Saad, M.D.; Kovtunenko, N.G.; Gomatos, P.J.; Zeman, V.V.; Botros, B.A.; Sanchez, J.L.; Birx, D.L.; et al. Genetic characterization of HIV-1 strains circulating in Kazakhstan. J. Acquir. Immune. Defic. Syndr. 2007, 46, 19–23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Moskaleychik, F.F.; Laga, V.Y.; Delgado, E.; Vega, Y.; Fernandez-Garcia, A.; Perez, A.; Kornilaeva, G.V.; Pronin, A.Y.; Zhernov, Y.V.; Thomson, M.M.; et al. Rapid spread of the HIV-1 circular recombinant CRF02-AG in Russia and neighboring countries. Vopr. Virol. 2015, 60, 14–19. [Google Scholar]
  37. Kostaki, E.G.; Karamitros, T.; Bobkova, M.; Oikonomopoulou, M.; Magiorkinis, G.; Garcia, F.; Hatzakis, A.; Paraskevis, D. Spatiotemporal Characteristics of the HIV-1 CRF02_AG/CRF63_02A1 Epidemic in Russia and Central Asia. AIDS Res. Hum. Retrovir. 2018, 34, 415–420. [Google Scholar] [CrossRef]
  38. Dehne, K.L.; Khodakevich, L.; Hamers, F.F.; Schwartlander, B. The HIV/AIDS epidemic in eastern Europe: Recent patterns and trends and their implications for policy-making. Aids 1999, 13, 741–749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Dehne, K.L.; Pokrovskiy, V.; Kobyshcha, Y.; Schwartlander, B. Update on the epidemics of HIV and other sexually transmitted infections in the newly independent states of the former Soviet Union. Aids 2000, 14 (Suppl. 3), S75–S84. [Google Scholar]
  40. Giyassova, G.M. National program for HIV prevention in Republic of Uzbekistan, 2007–2011. In Proceedings of the International Conference on the HIV Epidemic in Central Asia and Neighbouring Countries, Bishkek, Kyrgyzstan, 11–13 March 2008. [Google Scholar]
  41. Mamatov, M.S. Status of the HIV epidemic: Infections and medical treatment of persons living with HIV in the Kyrgyz Republic. In Proceedings of the International Conference on the HIV Epidemic in Central Asia and Neighbouring Countries, Bishkek, Kyrgyzstan, 11–13 March 2008. [Google Scholar]
  42. Khasanova, M.A. Analysis of HIV and AIDS epidemiologic situation in Kazakhstan. In Proceedings of the International Conference on the HIV Epidemic in Central Asia and Neighbouring Countries, Bishkek, Kyrgyzstan, 11–13 March 2008. [Google Scholar]
  43. Bobkov, A.; Cheingsong-Popov, R.; Selimova, L.; Ladnaya, N.; Kazennova, E.; Kravchenko, A.; Fedotov, E.; Saukhat, S.; Zverev, S.; Pokrovsky, V.; et al. An HIV type 1 epidemic among injecting drug users in the former Soviet Union caused by a homogeneous subtype A strain. AIDS Res. Hum. Retrovir. 1997, 13, 1195–1201. [Google Scholar] [CrossRef] [PubMed]
  44. Lukashov, V.V.; Karamov, E.V.; Eremin, V.F.; Titov, L.P.; Goudsmit, J. Extreme founder effect in an HIV type 1 subtype A epidemic among drug users in Svetlogorsk, Belarus. AIDS Res. Hum. Retrovir. 1998, 14, 1299–1303. [Google Scholar] [CrossRef] [PubMed]
  45. Novitsky, V.A.; Montano, M.A.; Essex, M. Molecular epidemiology of an HIV-1 subtype A subcluster among injection drug users in the Southern Ukraine. AIDS Res. Hum. Retrovir. 1998, 14, 1079–1085. [Google Scholar] [CrossRef]
  46. Pandrea, I.; Descamps, D.; Collin, G.; Robertson, D.L.; Damond, F.; Dimitrienco, V.; Gheorghita, S.; Pecec, M.; Simon, F.; Brun-Vezinet, F.; et al. HIV type 1 genetic diversity and genotypic drug susceptibility in the Republic of Moldova. AIDS Res. Hum. Retrovir. 2001, 17, 1297–1304. [Google Scholar] [CrossRef]
  47. Lapovok, I.; Kazennova, E.; Laga, V.; Vasilyev, A.; Utegenova, A.; Abishev, A.; Dzissyuk, N.; Tukeev, M.; Bobkova, M. Short communication: Molecular epidemiology of HIV type 1 infection in Kazakhstan: CRF02_AG prevalence is increasing in the southeastern provinces. AIDS Res. Hum. Retrovir. 2014, 30, 769–774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Sivay, M.V.; Totmenin, A.V.; Zyryanova, D.P.; Osipova, I.P.; Nalimova, T.M.; Gashnikova, M.P.; Ivlev, V.V.; Meshkov, I.O.; Chokmorova, U.Z.; Narmatova, E.; et al. Characterization of HIV-1 Epidemic in Kyrgyzstan. Front. Microbiol. 2021, 12, 753675. [Google Scholar] [CrossRef]
  49. Lapovok, I.A.; Lopatukhin, A.E.; Kireev, D.E.; Kazennova, E.V.; Lebedev, A.V.; Bobkova, M.R.; Kolomeets, A.N.; Turbina, G.I.; Shipulin, G.A.; Ladnaya, N.N.; et al. Molecular epidemiological analysis of HIV-1 variants circulating in Russia in 1987–2015. Ter. Arkh. 2017, 89, 44–49. [Google Scholar] [CrossRef] [PubMed]
  50. Lebedev, A.; Lebedeva, N.; Moskaleychik, F.; Pronin, A.; Kazennova, E.; Bobkova, M. Human Immunodeficiency Virus-1 Diversity in the Moscow Region, Russia: Phylodynamics of the Most Common Subtypes. Front. Microbiol. 2019, 10, 320. [Google Scholar] [CrossRef]
  51. Schlosser, M.; Kartashev, V.V.; Mikkola, V.H.; Shemshura, A.; Saukhat, S.; Kolpakov, D.; Suladze, A.; Tverdokhlebova, T.; Hutt, K.; Heger, E.; et al. HIV-1 Sub-Subtype A6: Settings for Normalised Identification and Molecular Epidemiology in the Southern Federal District, Russia. Viruses 2020, 12, 475. [Google Scholar] [CrossRef] [Green Version]
  52. Kazennova, E.V.; Neshumaev, D.A.; Rukavitsin, D.V.; Lapovok, I.A.; Laga, V.; Rakcheeva, O.V.; Vasil’ev, A.V.; Tumanov, A.S.; Vorontsova, G.A.; Kuznetsova, A.V.; et al. Molecular epidemiological analysis of the HIV infection in the Blagoveshchensk and Khabarovsk Area (Russian Far East). Vopr. Virusol. 2014, 59, 31–36. [Google Scholar] [PubMed]
  53. Diez-Fuertes, F.; Cabello, M.; Thomson, M.M. Bayesian phylogeographic analyses clarify the origin of the HIV-1 subtype A variant circulating in former Soviet Union’s countries. Infect. Genet. Evol. 2015, 33, 197–205. [Google Scholar] [CrossRef] [PubMed]
  54. Joint United Nations Programme on HIV and AIDS (UNAIDS). UNAIDS: The Response to HIV in Eastern Europe and Central Asia—Global AIDS Update 2019. Available online: https://www.unaids.org/sites/default/files/media_asset/2019_Regional_GR_Eastern-Europe-and-central-Asia_en.pdf (accessed on 20 May 2021).
  55. Rechel, B. HIV/AIDS in the Countries of the Former Soviet Union: Societal and Attitudinal Challenges. Cent. Eur. J. Public Health 2010, 18, 110–115. [Google Scholar] [CrossRef]
  56. Spicer, N.; Bogdan, D.; Brugha, R.; Harmer, A.; Murzalieva, G.; Semigina, T. ‘It’s risky to walk in the city with syringes’: Understanding access to HIV/AIDS services for injecting drug users in the former Soviet Union countries of Ukraine and Kyrgyzstan. Glob. Public Health 2011, 7, 22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Ancker, S.; Rechel, B. Policy responses to HIV/AIDS in Central Asia. Glob. Health 2015, 10, 817–833. [Google Scholar] [CrossRef]
Viruses 14 01675 g001 550
Figure 1. Subtypes, transmission route and geographical composition of Uzbek HIV-1 sequences enrolled in the study. Diagram (top panel) showing the number of two principal HIV-1 subtypes sequences and transmission route. Map of Uzbekistan (bottom panel) with region names and the proportions of CRF02_AG and A6 subtype sequences listed on the bubble; the circle size is proportional to the number of sequences in corresponding regions.
Figure 1. Subtypes, transmission route and geographical composition of Uzbek HIV-1 sequences enrolled in the study. Diagram (top panel) showing the number of two principal HIV-1 subtypes sequences and transmission route. Map of Uzbekistan (bottom panel) with region names and the proportions of CRF02_AG and A6 subtype sequences listed on the bubble; the circle size is proportional to the number of sequences in corresponding regions.
Viruses 14 01675 g001
Viruses 14 01675 g002 550
Figure 2. Network clustering results for the HIV-infected persons sequences in Uzbekistan generated by MicrobeTrace. The number and size of individual clusters for persons with HIV-1 CRF02_AG (A) and subtype A6 infection (B) are shown. Gender and transmission risk factors have been mapped to node shape and color, respectively, as described in the legend. Each node corresponds to an individual person, and each line represents two persons having a genetic distance 1.5% (d = 0.015). Different clusters are displayed in different polygon using a grey color gradient. Abbreviation: HET, heterosexual contacts; IDU, injection drug users; UNK, unknown.
Figure 2. Network clustering results for the HIV-infected persons sequences in Uzbekistan generated by MicrobeTrace. The number and size of individual clusters for persons with HIV-1 CRF02_AG (A) and subtype A6 infection (B) are shown. Gender and transmission risk factors have been mapped to node shape and color, respectively, as described in the legend. Each node corresponds to an individual person, and each line represents two persons having a genetic distance 1.5% (d = 0.015). Different clusters are displayed in different polygon using a grey color gradient. Abbreviation: HET, heterosexual contacts; IDU, injection drug users; UNK, unknown.
Viruses 14 01675 g002
Viruses 14 01675 g003 550
Figure 3. Maximum-likelihood phylogeny and inferred geographical transmission routes of Uzbek HIV-1 sequences with publicly available sequences from FSU and European countries. Shown is the analysis of HIV-1 CRF02_AG (left panel; A,B) and sub-subtypes A6 (right panel; C,D) (A,B). ML-tree shows a time calibrated phylogeny and revealing similarities of most closely related to Uzbekistan global sequences from a BLAST search at GenBank, including sequence from Caucasian and Central Asian republics, and Europe. Branches are colored by country of sampling as indicated in the legends. Branch lengths of ML tree are drawn to scale with the concentric circles indicating calendar years. The trees were rooted through HIV-1 subtype C. (C,D). The migration patterns of the HIV-1 strains in the Uzbekistan and neighboring countries. The image shows the HIV-1 transmissions (migrations) to/from different geographic areas that are colored according to the legends. The circle size is proportional to the number of sequences in corresponding regions. The graphics were generated using the auspice (https://auspice.us, accessed on 20 May 2021).
Figure 3. Maximum-likelihood phylogeny and inferred geographical transmission routes of Uzbek HIV-1 sequences with publicly available sequences from FSU and European countries. Shown is the analysis of HIV-1 CRF02_AG (left panel; A,B) and sub-subtypes A6 (right panel; C,D) (A,B). ML-tree shows a time calibrated phylogeny and revealing similarities of most closely related to Uzbekistan global sequences from a BLAST search at GenBank, including sequence from Caucasian and Central Asian republics, and Europe. Branches are colored by country of sampling as indicated in the legends. Branch lengths of ML tree are drawn to scale with the concentric circles indicating calendar years. The trees were rooted through HIV-1 subtype C. (C,D). The migration patterns of the HIV-1 strains in the Uzbekistan and neighboring countries. The image shows the HIV-1 transmissions (migrations) to/from different geographic areas that are colored according to the legends. The circle size is proportional to the number of sequences in corresponding regions. The graphics were generated using the auspice (https://auspice.us, accessed on 20 May 2021).
Viruses 14 01675 g003
Viruses 14 01675 g004 550
Figure 4. The in-country migration patterns of the HIV-1 strains (Uzbekistan). The image shows the HIV-1 sub-subtype A6 (A) and CRF02_AG (B) transmissions (migrations) to/from different geographic areas that are represented a three-letter code/color code: And, Andijan; Buk, Bukhara; Fer, Fergana; Jiz, Jizzakh; Kar, Karakalpakstan; Xor, Xorazm; Nam, Namangan; Nav, Navoiy; Qas, Qashqadaryo; Sam, Samarkand; Sir, Sirdaryo; Sur, Surxondaryo; Tas, Tashkent.
Figure 4. The in-country migration patterns of the HIV-1 strains (Uzbekistan). The image shows the HIV-1 sub-subtype A6 (A) and CRF02_AG (B) transmissions (migrations) to/from different geographic areas that are represented a three-letter code/color code: And, Andijan; Buk, Bukhara; Fer, Fergana; Jiz, Jizzakh; Kar, Karakalpakstan; Xor, Xorazm; Nam, Namangan; Nav, Navoiy; Qas, Qashqadaryo; Sam, Samarkand; Sir, Sirdaryo; Sur, Surxondaryo; Tas, Tashkent.
Viruses 14 01675 g004
Table
Table 1. Multivariable analysis of selected demographic characteristics of persons infected with HIV-1 CRF02_AG and sub-subtype A6 in Uzbekistan.
Table 1. Multivariable analysis of selected demographic characteristics of persons infected with HIV-1 CRF02_AG and sub-subtype A6 in Uzbekistan.
CharacteristicsEstimateStandard Errort Valuep ValueOdds Ratio
[95% CI]
Gender
(Male vs. Female)		−0.1880.304−0.6180.5360.828 [0.454–1.509]
Age (in years)0.0210.0161.3100.1911.022 [0.989–1.056]
Region
(NER vs. Other) 1		0.9630.3023.1860.0012.620 [1.443–4.757]
Transmission category
(IDU vs. HET)		0.2540.2970.8530.3941.289 [0.716–2.320]
1 NER (north-eastern region): Andijan, Fergana, Namangan and Tashkent region. Abbreviations: HET, heterosexual contacts; IDU, injection drug users; CI, confidence interval.
Table
Table 2. Characterization of Uzbek HIV-1 CRF02_AG and sub-subtype A6 transmissions clusters.
Table 2. Characterization of Uzbek HIV-1 CRF02_AG and sub-subtype A6 transmissions clusters.
CRF02_AGSub-Subtype A6
Cluster12Dyads12Dyads
Cluster (dyads) size, N (%)49 (25.5)3 (13.5)10 (13.5)36 (28.3)3 (8.7)8 (8.7)
Gender, N (%)
  Male16 (32.7)2 (66.7)7 (35.0)19 (52.7)1 (33.3)4 (50.0)
  Female33 (67.3)1 (33.3)13 (65.0)17 (47.3)2 (66.7)4 (50.0)
Transmission route, N (%)
  HET15 (30.6)3 (100.0)11 (55.0)14 (38.9)1 (33.4)2 (25.0)
  IDU19 (38.8)06 (30.0)13 (36.1)1 (33.3)2 (25.0)
  UNK15 (30.6)03 (15.0)9 (25.0)1 (33.3)4 (50.0)
Residence or region of sampling Viruses 14 01675 i001 Viruses 14 01675 i002N/a Viruses 14 01675 i003 Viruses 14 01675 i004N/a
Diagnosis date2013–201620152002–20152002–201620152013–2015
Likely phylogenetic origin (country)UzbekistanUzbekistanN/aRussiaRussiaN/a
TMRCA (95% CI)02-03-2003
[28-04-2002, 12-12-2004]		27-12-2005
[20-03-2004, 30-06-2008]		N/a08-02-2000
[11-02-1999, 31-07-2001]		07-03-2009
[12-09-2007, 08-01-2011]		N/a
Abbreviations: HET, heterosexual contacts; IDU, injection drug users; UNK, unknown; TMRCA, time to the most recent common ancestor; CI, confidence interval; N/a, not assessed.
Table
Table 3. Assortativity for the HIV-1 CRF02_AG and sub-subtype A6 clusters in Uzbekistan by gender, region and transmission category.
Table 3. Assortativity for the HIV-1 CRF02_AG and sub-subtype A6 clusters in Uzbekistan by gender, region and transmission category.
CRF063_02ASub-Subtype A6
Cluster1Total 11Total
Number of Nodes, N49723647
Assortativity, r
Region0.290.300.270.21
Sex0.030.03−0.08−0.09
Transmission route−0.04<0.01−0.18−0.06
1 Estimates across the clusters and dyads in aggregate.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lebedev, A.; Kuznetsova, A.; Kim, K.; Ozhmegova, E.; Antonova, A.; Kazennova, E.; Tumanov, A.; Mamatkulov, A.; Kazakova, E.; Ibadullaeva, N.; et al. Identifying HIV-1 Transmission Clusters in Uzbekistan through Analysis of Molecular Surveillance Data. Viruses 2022, 14, 1675. https://doi.org/10.3390/v14081675

AMA Style

Lebedev A, Kuznetsova A, Kim K, Ozhmegova E, Antonova A, Kazennova E, Tumanov A, Mamatkulov A, Kazakova E, Ibadullaeva N, et al. Identifying HIV-1 Transmission Clusters in Uzbekistan through Analysis of Molecular Surveillance Data. Viruses. 2022; 14(8):1675. https://doi.org/10.3390/v14081675

Chicago/Turabian Style

Lebedev, Aleksey, Anna Kuznetsova, Kristina Kim, Ekaterina Ozhmegova, Anastasiia Antonova, Elena Kazennova, Aleksandr Tumanov, Adkhamjon Mamatkulov, Evgeniya Kazakova, Nargiz Ibadullaeva, and et al. 2022. "Identifying HIV-1 Transmission Clusters in Uzbekistan through Analysis of Molecular Surveillance Data" Viruses 14, no. 8: 1675. https://doi.org/10.3390/v14081675

APA Style

Lebedev, A., Kuznetsova, A., Kim, K., Ozhmegova, E., Antonova, A., Kazennova, E., Tumanov, A., Mamatkulov, A., Kazakova, E., Ibadullaeva, N., Brigida, K., Musabaev, E., Mustafaeva, D., Rakhimova, V., & Bobkova, M. (2022). Identifying HIV-1 Transmission Clusters in Uzbekistan through Analysis of Molecular Surveillance Data. Viruses, 14(8), 1675. https://doi.org/10.3390/v14081675

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop