Canadian Journal of Gastroenterology and Hepatology

Canadian Journal of Gastroenterology and Hepatology / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 9701920 | https://doi.org/10.1155/2020/9701920

Lazzat Aibekova, Aizada Bexeitova, Arailym Aldabergenova, Gonzalo Hortelano, Zhangwen Ge, Feng Yi, Yiming Shao, Jack DeHovitz, Sten H. Vermund, Syed Ali, "Transmission of HIV and HCV within Former Soviet Union Countries", Canadian Journal of Gastroenterology and Hepatology, vol. 2020, Article ID 9701920, 10 pages, 2020. https://doi.org/10.1155/2020/9701920

Transmission of HIV and HCV within Former Soviet Union Countries

Academic Editor: Kevork M. Peltekian
Received17 Jan 2020
Accepted30 Apr 2020
Published15 Jul 2020

Abstract

Background. Following the collapse of the Union of Soviet Socialist Republic (USSR) in 1991, trans-border mobility increased within the former Soviet Union (FSU) countries. In addition, drug-trafficking and injection drug use began to rise, leading to the propagation and transmission of blood-borne infections within and across the FSU countries. To examine the transmission of blood-borne infections within this region, we analyzed the phylogenetic relationship of publically available sequences of two blood-borne viruses, hepatitis C virus (HCV) and human immunodeficiency virus (HIV), from FSU countries. Methods. We analysed 614 and 295 NS5B sequences from HCV genotypes 1b and 3a, respectively, from 9 FSU countries. From 13 FSU countries, we analysed 347 HIV gag and 1282 HIV env sequences. To examine transmission networks and the origins of infection, respectively, phylogenetic and Bayesian analyses were performed. Results. Our analysis shows intermixing of HCV and HIV sequences, suggesting transmission of these viruses both within and across FSU countries. We show involvement of three major populations in transmission: injection drug user, heterosexual, and trans-border migrants. Conclusion. This study highlights the need to focus harm reduction efforts toward controlling transmission of blood-borne infections among the abovementioned high-risk populations in the FSU countries.

1. Introduction

Following the collapse of the Union of Soviet Socialist Republics (USSR) in 1991, the ensuing economic crisis led to poverty and unemployment in the former Soviet republics. Existing social and cultural ties among the former Soviet Union (FSU) countries and visa-free travel across borders facilitated massive movement of migrants in search of employment [1]. High migration rates in the setting of economic destabilization were accompanied by increased rates of injected drug use, facilitating the transmission of blood-borne viruses such as human immunodeficiency virus (HIV) and hepatitis C virus (HCV) in the region [2, 3]. Due to their geographic location along drug-trafficking routes from Afghanistan, the main hub of opium production and supply for Russia and Europe, there has been increased trafficking and use of injectable drugs in Central Asia [4]. According to the World Health Organization (WHO), in Eastern Europe, 6.8 million people were estimated in 2015 to be positive for antibodies to HCV (3.3% prevalence) and 4.7 million people were living with chronic HCV (2.3% prevalence; 69% viremia rate), while in Central Asia, these figures were 4.5 million (5.4% prevalence) and 1.9 million (2.3% prevalence; 43% viremia rate) people [5]. The number of people living with HIV in Eastern European and Central Asian (EECA) countries, representing the only region in the world with rising HIV incidence, reached 1.6 million by 2016 [6].

The objective of this study was to investigate the epidemiology of blood-borne viruses, namely, HIV and HCV, among FSU countries. Using viral sequences from public databases, we have performed phylogenetic analysis to assess common routes of transmission of these two viruses within FSU countries.

2. Methods

2.1. Selection and Downloading of Sequences

For HCV, a 234 bp fragment of NS5B gene (corresponding to H77 8322–8555 nt) was studied. We downloaded 614 sequences for genotype 1b and 295 for genotype 3a from the Los Alamos HCV database (http://www.hcv.lanl.gov). Sequences from 9 FSU countries, namely, Russia, Uzbekistan, Tajikistan, Azerbaijan, Belarus, Lithuania, Latvia, Estonia, and Georgia, were used for this analysis. Sequence retrieval and length decisions were based on selecting the most represented gene fragment (and genotype sequences) available for most FSU countries in the said database (Supplementary Figure 1). Genotypes for all the sequences were ascertained using the Oxford HCV Automated Subtyping Tool 2.0 (http://www.bioafrica.net/rega-genotype/html/subtypinghcv.html). For constructing phylogenetic tress, 10 known genotype sequences were used as reference.

For HIV, subtype A gag and env sequences from Los Alamos HIV Sequence Database (http://www.hiv.lanl.gov) were downloaded. Available sequences from 13 FSU countries, namely, Armenia, Azerbaijan, Belarus, Estonia, Georgia, Kazakhstan, Kyrgyzstan, Latvia, Lithuania, Moldova, Russia, Ukraine, and Uzbekistan, were retrieved (Supplementary Figure 1). Recombinant and duplicate sequences were identified using, respectively, Recombinant HIV-1 Drawing Tool v2.1.0 and ElimDupes (http://www.hiv.lanl.gov), and eliminated. After alignment and trimming, 371 bp long, 334 gag sequences (HXB2: 859–1230 nt), and 257 bp long, 1282 env sequences (HXB2: 7071–7336 nt) were used for this analysis. To root the tree, 12 reference sequences of nonsubtype A and group M were used from Los Alamos HIV Sequence Database (http://www.hiv.lanl.gov). To construct the tree with risk group information, data from the HIV Los Alamos database were used.

2.2. Phylogenetic Analysis

SeaView v4.6.1 was used to align sequences using the ClustalO algorithm, and the alignments were manually edited using MEGA v6.06. The phylogenetic trees were constructed using maximum likelihood method with 100 replicates bootstrap evaluation under the Tamura-Nei model as nucleotide substitution model, including a uniform rate of heterogeneity among sites, as implemented in MEGAv6.06. The final tree was refined and visualized using Figtree 1.4.2 software (http://tree.bio.ed.ac.uk/software/). The clusters were picked based on the predominance of sequences from a particular country or risk group.

2.3. Bayesian Analysis

To investigate temporal signal in the data set, we used the maximum-likelihood tree for the analysis of correlation between root-to-tip genetic distance and year of sampling in the program TempEst v1.5. A Bayesian phylogenetic approach [7] was used for joint estimation of the ages of each of clusters and the demographic history of all of the strains. This was done by analyzing these sequences using the general time reversible (GTR) model [8] plus a gamma distribution among site rate heterogeneities [9] on the basis of the standard Akaike information criterion in a hierarchical likelihood ratio test in our study. Bayesian Markov chain Monte Carlo (MCMC) analyses were performed with the selected nucleotide substitution model using an uncorrelated lognormal relaxed-clock model [10] as implemented in BEAST, version 1.8.4. The analyses were performed using a Bayesian SkyGrid coalescent tree prior [11, 12] in order to investigate the demographic histories of HCV or HIV and the degree to which dating estimates are affected by the demographic model chosen. Sequences for which sampling year information was unavailable were excluded from the analysis.

To ensure convergence of parameter estimates, the MCMC was run for 100 million steps and sampled every 10000 steps. The resulting MCMC samples under different demographic models were used to estimate the rates of evolution and tMRCAs with the first 10% removed as burn-in. Convergence of parameters and sufficient sampling were assessed by calculating the effective sample size (ESS) using TRACER v1.6 (http://beast.bio.ed.ac.uk/software/tracer). A final maximum clade credibility (MCC) tree was constructed from the posterior tree distribution using the program Tree Annotator (http://beast.bio.ed.ac.uk) and visualized in FigTree.

3. Results

3.1. HCV

The maximum likelihood trees generated using genotype 1b or 3a NS5B sequences from 7 countries revealed phylogenetic links of Russian clusters with sequences from other FSU countries. The sequences analysed were deposited during 1995–2015 (Tables 1 and 2). Overall, 10 clusters were picked from HCV 1b NS5B tree based on the predominance of sequences from a particular country or risk group (Table 1 and Figure 1): cluster A (53 sequences, representing 5 countries), cluster B (37 sequences, representing 7 countries), cluster C (18 sequences, representing 4 countries), cluster D (23 sequences, representing 6 countries), cluster E (133 sequences, representing 6 countries), cluster F (33 sequences, representing 6 countries), cluster G (30 sequences, representing 5 countries), cluster H (31 sequences, representing 4 countries), cluster I (39 sequences, representing 6 countries), and cluster J (32 sequences, representing 5 countries). Genotype 1b tree showed intermixing of Russian sequences mainly with those from Estonia (clusters A, B, and D–J), Lithuania (clusters B-J), Azerbaijan (clusters A-E, G, and I), Uzbekistan (clusters A-G and I), and Tajikistan (clusters A, B, D–F, and H–J) (Figure 1, and Table 1). Seven clusters were picked from HCV 3a NS5B tree based on the predominance of sequences from a particular country or risk group (Table 2 and Figure 2): cluster A (44 sequences, representing 6 countries), cluster B (11 sequences, representing 3 countries), cluster C (11 sequences, representing 4 countries), cluster D (18 sequences, representing 3 countries), cluster E (23 sequences, representing 3 countries), cluster F (80 sequences, representing 6 countries), and cluster G (33 sequences, representing 4 countries). Genotype 3a tree showed a similar trend with an intermixing of HCV sequences from Russia mainly with those from Estonia (clusters A and C–G), Uzbekistan (clusters A, C, and E–G), and Azerbaijan (clusters A, B, D, F, and G) (Figure 2 and Table 2).


ClusterNumber of sequencesNumber of countriesCountriesYears

A535Uzbekistan, Russia, Estonia, Azerbaijan, Tajikistan1997–2014
B377Uzbekistan, Russia, Estonia, Azerbaijan, Tajikistan, Georgia, Lithuania1997–2014
C184Uzbekistan, Russia, Azerbaijan, Lithuania1999–2014
D236Uzbekistan, Russia, Estonia, Azerbaijan, Tajikistan, Lithuania1997–2014
E1336Uzbekistan, Russia, Estonia, Azerbaijan, Tajikistan, Lithuania1995–2014
F335Uzbekistan, Russia, Estonia, Tajikistan, Lithuania1997–2006
G305Uzbekistan, Russia, Estonia, Azerbaijan, Lithuania1998–2014
H314Russia, Estonia, Tajikistan, Lithuania1998–2015
I396Uzbekistan, Russia, Estonia, Azerbaijan, Tajikistan, Lithuania1995–2014
J325Uzbekistan, Russia, Estonia, Tajikistan, Lithuania1997–2014


ClusterNumber of sequencesNumber of countriesCountriesYears

A446Uzbekistan, Russia, Estonia, Lithuania, Azerbaijan, Georgia1998–2012
B113Lithuania, Russia, Azerbaijan1998–2011
C114Uzbekistan, Russia, Lithuania, Estonia2001–2006
D183Russia, Estonia, Azerbaijan1998–2009
E233Uzbekistan, Russia, Estonia1998–2009
F806Uzbekistan, Russia, Estonia, Lithuania, Azerbaijan, Tajikistan1999–2014
G334Uzbekistan, Russia, Azerbaijan, Estonia1999–2008

3.2. HIV

A phylogenetic tree was constructed with env sequences from 10 counties, spanning years 1993–2015, as follows: Armenia (2009), Belarus (1997–2015), Georgia (1997–2013), Estonia (2001–2005), Kyrgyzstan (2009), Latvia (1998–2005), Moldova (1997), Russia (1994–2015), Ukraine (1993–2012), and Uzbekistan (1999–2002). In addition, 12 nonsubtype A references were retrieved from the database and were used as outliers to root the phylogenetic tree. Fifteen clusters were picked from the tree for further analysis (Table 3 and Figure 3): cluster A (76 sequences, representing 5 countries), cluster B (53 sequences, representing 4 countries), cluster C (25 sequences, representing 2 countries), cluster D (83 sequences, representing 6 countries), cluster E (13 sequences, representing 3 countries), cluster F (111 sequences, representing 4 countries), cluster G (143 sequences, representing 5 countries), cluster H (8 sequences, representing 4 countries), cluster I (17 sequences, representing 5 countries), cluster J (55 sequences, representing 4 countries), cluster K (130 sequences, representing 4 countries), cluster L (16 sequences, representing 2 countries), cluster M (84 sequences, representing 5 countries), cluster N (16 sequences, representing 3 countries), and cluster O (6 sequences, representing 1 countries). More than half of the sequences used to construct the tree were from Russia; therefore, the results show a predominance of Russian branches in most clusters (Figure 3 and Table 3). Interestingly, sequences from Russia, Uzbekistan, and Latvia were intermixed in almost all the clusters, whereas sequences from Armenia were also represented prominently in most clusters (Figure 3 and Table 3).


ClusterNumber of sequencesNumber of countriesCountriesYears

A765Russia, Uzbekistan, Estonia, Kyrgyzstan, Armenia1999–2012
B534Russia, Latvia, Uzbekistan, Estonia1999–2015
C252Russia, Latvia2001–2012
D836Russia, Latvia, Uzbekistan, Belarus, Ukraine, Moldova1996–2012
E133Russia, Uzbekistan, Ukraine1995–2012
F1114Russia, Latvia, Armenia, Estonia1998–2013
G1435Russia, Latvia, Uzbekistan, Armenia, Ukraine1993–2015
H84Russia, Latvia, Uzbekistan, Armenia1999–2009
I175Russia, Latvia, Uzbekistan, Moldova, Armenia1997–2011
J554Russia, Latvia, Armenia, Belarus2004–2015
K1304Russia, Latvia, Armenia, Ukraine2001–2015
L162Russia, Latvia2001–2015
M845Russia, Latvia, Armenia, Ukraine, Georgia1999–2015
N163Russia, Armenia, Ukraine2009–2015
O61Ukraine1993

The second tree was constructed with HIV gag sequences from 8 FSU countries spanning the years 1997–2013 as follows: Belarus (1997–2013), Georgia 1999 (1999), Kazakhstan (2002), Latvia (1998–2005), Russia (1997–2011), Ukraine (2000–2012), and Uzbekistan (2002). The majority of gag sequences were from Russia and Latvia, whereas only one sequence was from Georgia. Twelve clusters were picked from the tree for the further analysis (Table 4 and Figure 4): cluster A (10 sequences, representing 1 countries), cluster B (13 sequences, representing 3 countries), cluster C (9 sequences, representing 4 countries), cluster D (7 sequences, representing 3 countries), cluster E (6 sequences, representing 4 countries), cluster F (17 sequences, representing 1 countries), cluster G (33 sequences, representing 4 countries), cluster H (13 sequences, representing 1 countries), cluster I (12 sequences, representing 2 countries), cluster J (20 sequences, representing 3 countries), cluster K (31 sequences, representing 1 countries), and cluster L (7 sequences, representing 3 countries). This tree clearly showed intermixing of Russian strains with Latvia (Figure 4). In agreement with the observations from the HIV env tree (Figure 3), in this tree as well, sequences from Russia and Latvia were intermixed in almost all the clusters throughout the tree (Figure 4 and Table 4). Additionally, as was observed in the env tree, the ancestral sequences, representing the earliest infections, in this tree as well originated from Ukraine, with some intermixing of sequences from Russia (Figures 3 and 4). Both in the env (cluster F) and gag (clusters A, F, H and K) trees, monophyletic clusters of Latvian sequences originating from Russian sequences were observed (Figures 3 and 4), possibly indicating epidemics in Latvia that resulted from migrant-associated transmissions from Russia.


ClusterNumber of sequencesNumber of countriesCountriesYears

A101Latvia2001–2005
B133Russia, Latvia, Ukraine2000–2003
C94Russia, Latvia, Kazakhstan, Uzbekistan2002–2005
D73Russia, Latvia, Ukraine2001–2011
E64Russia, Latvia, Ukraine, Belarus1997–2005
F171Latvia1998–2002
G334Russia, Latvia, Ukraine, Belarus1997–2006
H131Latvia2001–2003
I122Latvia, Russia2001–2005
J203Russia, Latvia, Ukraine1998–2008
K311Latvia2001–2005
L73Russia, Kazakhstan, Latvia2002–2006

3.3. Risk Group Analysis

To analyse modes of transmission, high-risk data associated with HIV env and gag sequences were also analysed. Only sequences with recorded high-risk information were included in this analysis. Both the env (Figure 5(a)) and gag (Figure 5(b)) trees showed a predominance of PWID-associated sequences, followed by those linked with heterosexual transmission. Results represented here showed intermixing of PWID and heterosexual routes of transmission, indicating the possibility of bridging of the initial PWID-epidemic into heterosexual populations.

3.4. Bayesian Analysis

As summarized in Table 5, our Bayesian phylogenetic analysis showed substitution rate for NS5B HCV 1b and HCV 3a to be 6.8 (5.8–7.9) × 10−3 and 1.1 (0.7–1.4) × 10−3 substitutions per site per year, respectively, with a SkyGrid coalescent model, whereas the estimated tMRCAs for HCV 1b and HCV 3a were 1977 (95% highest posterior density, HPD: 1962–1985) and 1965 (95% HPD: 1941–1978), respectively. The estimated evolutionary rates of HIV env and HIV gag were 3.5 (2.5–4.0) × 10−2 and 13 (8.8–17) × 10−3 substitutions per site per year, respectively. The likely years of origins of HIV env and gag were estimated to be 1992 (95% HPD: 1990–1993) and 1996 (95% HPD: 1995–1997), respectively.


tMRCA (median and 95% HPD interval)Substitution/site/year (median and 95% HPD interval)

HCV 1B37.7 (30.0, 52.6)6.8E – 3 (5.8E –3, 7.9E – 3)
HCV 3A49.3 (36.0, 73.1)1.1E – 3 (7.0E – 4, 1.4E – 3)
HIV ENV22.8 (22.0, 24.8)3.5E – 2 (2.5E – 2, 4.0E – 2)
HIV GAG16.9 (16.3, 17.5)1.3E – 2 (8.8E – 3, 1.7E – 2)

We further investigated the past population dynamics of HCV and HIV using a Bayesian skyline plot, which depicts the changes in effective population size over time. The effective population size of HCV 1b was observed to be fluctuating and slowly rising from 1998 to 2011. The change in effective population size for HCV 3a rose slowly and steadily starting in 1965, characterized by a phase of exponential growth (1981–1997), and the effective population size surged approximately 100-fold during the period. The epidemic history of HIV env and HIV gag seems to have experienced complex dynamics, characterized by two phases of wave-like fluctuating growth, during 1997–2005 and 2009–2014 for HIV env, and during 1995–2005 and 2009–2015 for HIV gag (Figure 6).

4. Discussion

In this study, we show evidence of transmission of two blood-borne viruses, HCV and HIV, within the FSU countries. According to our analysis, common routes of transmission appear to be injection drug use and heterosexual sex, with migrant mobility facilitating the cross-border spread of infection.

There have been reports of high prevalence of blood-borne infections in the FSU region. Within FSU countries, Russia has the highest rate of HCV genotype 1b (50%), followed by HCV 3a (45%) [13]. These two HCV genotypes also appear to be predominant in other FSU countries; Uzbekistan (1b, 64.2%; 3a, 25%), Azerbaijan (1b, 71.1%; 3a, 17%), Tajikistan (1b, 84.6%; 3a, 7.6%), Lithuania (1b, 54.3%; 3a, 23.9%), Estonia (1b, 71%; 3a, 24%), and Georgia (1b, 62%; 3a, 27%) [1419].

In our phylogenetic tree constructed with HCV 1b sequences (Figure 1), a preponderance of Russian sequences was seen—a reflection of a higher representation of sequences from this country in the Los Alamos database. Russian sequences were intermixed with clusters from most other FSU countries, with a relatively higher proportion from Lithuania and Estonia. This might reflect a higher transmigration between Russia, Lithuania, and Estonia, due to their close geographic proximity. Conversely, in clusters C, D, and F, sequences from Azerbaijan, Uzbekistan, and Tajikistan were seen intermixed with those from Russia, indicating that Russia might have served as a hub for transmission across these countries. Phylogenetic tree constructed with HCV 3a NS5B sequences (Figure 2) showed an overall agreement with the 1b tree, showing trends of possible transmission within Russia, Lithuania, Uzbekistan, Estonia, and Azerbaijan.

Such phylogenetic patterns may be supported by historical and socioeconomic factors that came into play after the disintegration of the Soviet Union. It has been established that the first nationalism-related conflicts between Armenia and Azerbaijan (1987–1991), followed by Uzbekistan conflict (1989), and civil wars in Georgia (1991–1993) and Tajikistan (1992–1997) led to mass migration from these countries to Russia [1]. Sequences from these countries that show intermixing with those from Russia also date to roughly the same time (Tables 1 and 2). Our Bayesian analysis predicts origin of the HCV epidemics in FSU countries to be as early as 1965, for HCV3a, and 1977, for genotype 1b (Table 5 and Figure 6), indicating that these HCV genotypes were possibly circulating in the region before the conflicts started. It is possible that HCV was transmitted during this time through injection use, leading to a rise in transmission in the mid-1990s, as suggested by our Skyline plot analysis (Figure 6). In addition, significant reduction of employment opportunities and a dramatic decline in wages in Tajikistan and Georgia also facilitated a high flow of migrant labor into Russia. Further movement involved repatriates and internally displaced people of Russian-speaking origin into Russia owing to a resurgence of nationalism as a result of implementation of national language policies that discouraged use of Russian. It is thought that the reported numbers of migrants in this region is much lower than the actual numbers, since visa-free movement across borders does not account for illegally entering migrant workers [1].

In the initial stages of HIV outbreak in central Asian region, most people infected with HIV acquired a homogenic variant of subtype A, known as AFSU (A6). This was thought to be a large-scale founder effect, resulting from transmigration within FSU countries [20]. The first HIV epidemic in FSU was recorded in the mid-1990s in Ukraine. A year later, outbreaks were detected in the European part of Russia, thereafter reaching other parts of Russia, Central Asia, and Eastern Europe.

The HIV epidemic expanded rapidly after the collapse of the USSR. As noted previously, political transitions led to the socioeconomic decline of newly formed independent counties. Consequently, large-scale labor migration to Russia from other countries of FSU provided epidemiological bridging of HIV transmission between Russia and the rest of FSU countries [20]. The database of international migration demonstrated a significant increase in migration between 1991 and 2000 within the FSU. The number of migrants from Central Asia, Caucasus, and Eastern Europe to Russia was highest from 1992 to 1996 [21], the same period during which our Skyline plot shows a gradual rise in the infections (Figure 6). Additionally, our phylogenetic trees revealed a dense intermixing of Russian sequences with other countries, dating around the same time (Figures 3 and 4 and Tables 3 and 4), which might suggest a link to migrant-associated HIV transmission. The HIV epidemic in FSU is considered to have initiated by transmissions through sexual contact and injection drug use [20]. Our results of phylogenetic tree in the context of high-risk behavior (Figures 5(a) and 5(b)) support that fact, showing the branches close to the ancestral node associated with heterosexual mode of transmission, with sequences represented from Uzbekistan (1999), Ukraine (1993), and Russia (1994).

As shown in previous reports, the initial HIV outbreak in FSU countries occurred in the mid-1990s in the city of Odessa, Ukraine, and spread to other countries of FSU [22]. Our phylogenetic tree also shows the ancestral node for the HIV subtype A epidemic originating from Ukrainian sequences deposited in 1993 (Figures 3 and 4), whereas the Bayesian analysis also predicts the origins of the epidemic to be around the same time (Table 1 and Figure 6). These results might indicate that Odessa had a favourable location for an explosive spread of HIV among PWID, as the city is situated on the Black Sea shore and was a major seaport and transportation hub during that period [20]. In the first postcollapse decade, there was an increase in injection drug use in the FSU [23], linked to the rise of drug production in Afghanistan during the mid-1990s and a shift in trafficking to the “Northern route,” which runs through Central Asia to Russia and South Caucasus to Eastern Europe [24, 25]. This increase in drug trafficking may have played a role in the early spread of HIV among PWID communities. Our phylogenetic tree also shows a high number of transmissions associated with drug use, later shifting the mode of transmission toward heterosexual sex. Furthermore, our phylogenetic analysis shows heavy intermixing of Armenian PWID sequences from the South Caucasian region, which is within “Northern route” of drug trafficking, with other Eastern European sequences. During the period of 1989–2001, another important scenario that may have facilitated transmission of HIV and HCV was the displacement of a population of about one million Armenians in aftermath of the 1988 earthquake, coupled with deteriorating socioeconomic situation in the 1990s. These displaced labor populations immigrated mainly to the Eastern European countries of Ukraine, Belarus, and Eastern part of Russia [26].

5. Conclusion

Our analyses highlight the main routes of transmission for two important blood-borne viruses in FSU countries. This information may be exploited to refine Public Health policies to better manage the infected populations and prevent further spread of the infections. Although this study was limited by the availability of the number of sequences from FSU countries in the open-access databases, the results highlight the importance of focusing the harm reduction efforts on communities where blood-borne infections might be currently spreading, namely, persons who inject drugs, heterosexual populations, migrant workers, and cross-border travelers.

Data Availability

The data used to support the findings of this study are available from the corresponding author.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Lazzat Aibekova, Aizada Bexeitova, and Arailym Aldabergenova have contributed equally to the manuscript.

Acknowledgments

Funding for this study was partially provided by Nazarbayev University under grant: 110119FD4516. SHV is supported, in part, by NIH Grant No. P30MH062294. JD was partially supported by Fogarty International Center and the National Institute on Drug Abuse of the National Institutes of Health under Award Number D43 TW010046.

Supplementary Materials

Supplementary Figure 1: numbers of HIV and HCV sequences available in Los Alamos Database are shown for each country. HCV Genotype 1b NS5B, HCV Genotype 3a NS5B, HIV Subtype A env, and HIV Subtype A gag are color-coded in boxes as yellow, green, light blue, and dark blue, respectively. (Supplementary Materials)

References

  1. O. Choudinovskikh and M. Denissenko, Migration between CIS Countries: Trends and Policy, National Research University, Moscow, Russia, 2014.
  2. V. Laga, A. Vasilyev, I. Lapovok et al., “HIV type 1 subtype A1 dominates in Armenia,” Current HIV Research, vol. 13, no. 3, pp. 219–225, 2015. View at: Publisher Site | Google Scholar
  3. V. A. Manuylov, E. V. Chub, V. S. Kichatova et al., “Hepatitis C virus subtype 3a was introduced in the USSR in the early 1980s,” Journal of General Virology, vol. 98, no. 8, pp. 2079–2087, 2017. View at: Publisher Site | Google Scholar
  4. A. Renton, D. Gzirishvilli, G. Gotsadze, and J. Godinho, “Epidemics of HIV and sexually transmitted infections in Central Asia: trends, drivers and priorities for control,” International Journal of Drug Policy, vol. 17, no. 6, pp. 494–503, 2006. View at: Publisher Site | Google Scholar
  5. World Health Organization, Guidelines for the Screening, Care and Treatment of Persons with Chronic Hepatitis C Infection, World Health Organization, Geneva, Switzerland, 2016.
  6. UNAIDS, UNAIDS Data 2017, UNAIDS, Geneva, Switzerland, 2017, http://www.unaids.org/sites/default/files/media_asset/20170720_Data_book_2017_en.pdf.
  7. X. Li, H. Liu, L. Liu et al., “Tracing the epidemic history of HIV-1 CRF01_AE clusters using near-complete genome sequences,” Scientific Reports, vol. 7, no. 1, p. 4024, 2017. View at: Publisher Site | Google Scholar
  8. F. Rodríguez, J. L. Oliver, A. Marín, and J. R. Medina, “The general stochastic model of nucleotide substitution,” Journal of Theoretical Biology, vol. 142, no. 4, pp. 485–501, 1990. View at: Publisher Site | Google Scholar
  9. Z. Yang, “Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods,” Journal of Molecular Evolution, vol. 39, no. 3, pp. 306–314, 1994. View at: Publisher Site | Google Scholar
  10. A. J. Drummond, S. Y. W. Ho, M. J. Phillips, and A. Rambaut, “Relaxed phylogenetics and dating with confidence,” PLoS Biology, vol. 4, no. 5, p. e88, 2006. View at: Publisher Site | Google Scholar
  11. A. J. Drummond, A. Rambaut, B. Shapiro, and O. G. Pybus, “Bayesian coalescent inference of past population dynamics from molecular sequences,” Molecular Biology and Evolution, vol. 22, no. 5, pp. 1185–1192, 2005. View at: Publisher Site | Google Scholar
  12. N. R. Faria, A. Rambaut, M. A. Suchard et al., “The early spread and epidemic ignition of HIV-1 in human populations,” Science, vol. 346, no. 6205, pp. 56–61, 2014. View at: Publisher Site | Google Scholar
  13. O. V. Kalinina, “Hepatitis C virus: variability mechanisms, classification, evolution,” Voprosy Virusologii, vol. 60, no. 5, pp. 5–10, 2015. View at: Google Scholar
  14. F. Kurbanov, Y. Tanaka, F. Sugauchi et al., “Hepatitis C virus molecular epidemiology in Uzbekistan,” Journal of Medical Virology, vol. 69, no. 3, pp. 367–375, 2003. View at: Publisher Site | Google Scholar
  15. F. Bokharaei-Salim, H. Keyvani, S. H. Monavari, S. M. Alavian, S. Fakhim, and S. Nasseri, “Distribution of hepatitis C virus genotypes among Azerbaijani patients in capital city of Iran-tehran,” Hepatitis Monthly, vol. 13, no. 9, Article ID e13699, 2013. View at: Publisher Site | Google Scholar
  16. A. Khan, F. Kurbanov, Y. Tanaka et al., “Epidemiological and clinical evaluation of hepatitis B, hepatitis C, and delta hepatitis viruses in Tajikistan,” Journal of Medical Virology, vol. 80, no. 2, pp. 268–276, 2008. View at: Publisher Site | Google Scholar
  17. A. Ambrozaitis, K. Z̆agminas, G. Balc̆iūnaite, and A. Widell, “Hepatitis C in Lithuania: incidence, prevalence, risk factors and viral genotypes,” Clinical and Diagnostic Virology, vol. 4, no. 4, pp. 273–284, 1995. View at: Publisher Site | Google Scholar
  18. T. Tallo, H. Norder, V. Tefanova et al., “Genetic characterization of hepatitis C virus strains in Estonia: fluctuations in the predominating subtype with time,” Journal of Medical Virology, vol. 79, no. 4, pp. 374–382, 2007. View at: Publisher Site | Google Scholar
  19. A. Petruzziello, S. Marigliano, G. Loquercio, A. Cozzolino, and C. Cacciapuoti, “Global epidemiology of hepatitis C virus infection: an up-date of the distribution and circulation of hepatitis C virus genotypes,” World Journal of Gastroenterology, vol. 22, no. 34, pp. 7824–7840, 2016. View at: Publisher Site | Google Scholar
  20. M. Bobkova, “Current status of HIV-1 diversity and drug resistance monitoring in the former USSR,” AIDS Reviews, vol. 15, no. 4, pp. 204–212, 2013. View at: Google Scholar
  21. UN International Migration Flows to and from Selected Countries: The 2015 Revision, 2015.
  22. F. Díez-Fuertes, M. Cabello, and M. M. Thomson, “Bayesian phylogeographic analyses clarify the origin of the HIV-1 subtype A variant circulating in former Soviet Union’s countries,” Infection, Genetics and Evolution, vol. 33, pp. 197–205, 2015. View at: Publisher Site | Google Scholar
  23. T. Rhodes, A. Ball, G. V. Stimson et al., “HIV infection associated with drug injecting in the newly independent states, eastern Europe: the social and economic context of epidemics,” Addiction, vol. 94, no. 9, pp. 1323–1336, 1999. View at: Publisher Site | Google Scholar
  24. UNODC, Global Illicit Drug Trends: 2003, UNODC, New York, NY, USA, 2003.
  25. S. A. Farooq, M. H. Rasooly, S. H. Abidi, K. Modjarrad, and S. Ali, “Opium trade and the spread of HIV in the golden crescent,” Harm Reduction Journal, vol. 14, no. 1, p. 47, 2017. View at: Publisher Site | Google Scholar
  26. International Labour Organization, ILO Migration and Development: Armenia Country Study, International Labour Organization, Yerevan, Armenia, 2009.

Copyright © 2020 Lazzat Aibekova et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views357
Downloads303
Citations

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.