Spatiotemporal dynamics of HIV-1 CRF63_02A6 sub-epidemic

HIV-1 epidemic in Russia is one of the fastest growing in the world reaching 1.14 million people living with HIV-1 (PLWH) in 2021. Since mid-1990s, the HIV-1 epidemic in Russia has started to grow substantially due to the multiple HIV-1 outbreaks among persons who inject drugs (PWID) leading to expansion of the HIV-1 sub-subtype A6 (former Soviet Union (FSU) subtype A). In 2006, a local HIV-1 sub-epidemic caused by the distribution of novel genetic lineage CRF63_02A6 was identified in Siberia. In this study, we used a comprehensive dataset of CRF63_02A6 pol gene sequences to investigate the spatiotemporal dynamic of the HIV-1 CRF63_02A6 sub-epidemic. This study includes all the available CRF63_02A6 HIV-1 pol gene sequences from Los Alamos National Laboratory (LANL) HIV Sequence Database. The HIV-1 subtypes of those sequences were conferred using phylogenetic analysis, and two automated HIV-1 subtyping tools Stanford HIVdb Program and COMET. Ancestral state reconstruction and origin date were estimated using Nextstrain. Evolutionary rate and phylodynamic analysis were estimated using BEAST v 1.10.4. CRF63_02A6 was assigned for 872 pol gene sequences using phylogenetic analysis approach. Predominant number (n = 832; 95.4%) of those sequences were from Russia; the remaining 40 (4.6%) sequences were from countries of Central Asia. Out of 872 CRF63_02A6 sequences, the corresponding genetic variant was assigned for 75.7 and 79.8% of sequences by Stanford and COMET subtyping tools, respectively. Dated phylogenetic analysis of the CRF63_02A6 sequences showed that the virus most likely originated in Novosibirsk, Russia, in 2005. Over the last two decades CRF63_02A6 has been widely distributed across Russia and has been sporadically detected in countries of Central Asia. Introduction of new genetic variant into mature sub-subtype A6 and CRF02_AG_FSU epidemics could promote the increase of viral genetic diversity and emergence of new recombinant forms. Further HIV-1 studies are needed due to a continuing rapid virus distribution. Also, the implementation of HIV-1 prevention programs is required to reduce HIV-1 transmission. This study also highlights the discrepancies in HIV-1 subtyping approaches. The reference lists of HIV-1 sequences implemented in widely used HIV-1 automated subtyping tools need to be updated to provide reliable results.

Introduction

Russia, as well as other countries of Eastern Europe and Central Asia, has one of the fastest growing HIV-1 epidemics in the world (Frank et al., 2019). In 2021, 1.14 million people are living with HIV-1 (PLWH) in Russia (Central Research Institute of Epidemiology, Federal Service on Customers’ Rights Protection and Human Well-being Surveillance, 2021)¹. Russia has a diverse HIV-1 epidemic with unique patterns shaped by various socio-economical, geographical, and cultural factors influencing the development of local HIV-1 sub-epidemics within the country (Bobkova, 2013). Since mid-1990s, the HIV-1 epidemic in Russia has started to grow explosively due to the multiple HIV-1 outbreaks occurred among PWID leading to expansion of the HIV-1 sub-subtype A6, also known as former Soviet Union (FSU) subtype A (FSU-A). This genetic variant predominates in most of the Russian territories ever since (Bobkov et al., 2004; Thomson et al., 2009; Diez-Fuertes et al., 2015; Aibekova et al., 2018). In late-1990s, HIV-1 sub-epidemic among different high-risk groups was detected on the territories of the Russian Fast East (Kazennova et al., 2014; Kotova et al., 2019) and European part of Russia (Liitsola et al., 2000; Lebedev et al., 2020; Ozhmegova et al., 2020). Currently, along with the predominant sub-subtype A6, subtypes B, C, and CRF03_AB have been constantly identified (Karamov et al., 2018; Pasechnik and Blokh, 2018; Kotova et al., 2019; Safina et al., 2022). In 2006, a local HIV-1 sub-epidemic was caused by the distribution of the novel genetic lineage CRF63_02A6 in Siberia (Baryshev et al., 2012, 2014). The detailed analysis of CRF63_02A6 genome structure revealed that this viral variant originated due to recombination between Russian sub-subtype A6 and CRF02_AG_FSU (Central Asian lineage of the CRF02_AG; Carr et al., 2005; Baryshev et al., 2014). Several outbreaks of the CRF63_02A6 were detected in different Siberian regions between 2007 and 2014. This genetic variant has dominated in Siberia since that time (Gashnikova et al., 2015, 2017; Ulyanova et al., 2019; Kazennova et al., 2020; Maksimenko et al., 2020). The number of new HIV-1 cases has been gradually rising in Russia and has peaked in 2014, followed by the decline stage (Figure 1). At the same time, HIV-1 prevalence in Siberia was significantly higher than the country average.

FIGURE 1

Figure 1. Dynamics of HIV-1 incidence in Russia from 2006 to 2021. Data were adopted from the Central Scientific Research Institute of Epidemiology of Rospotrebnadzor, Russian Federal AIDS Center, Moscow, Russia. Incidence data collected in 2019–2021 was most likely biased due to the COVID-19 world pandemic and represented the under-registered new HIV-1 cases.

Initial HIV-1 epidemic in Russia occurred mainly among PWID and their sexual partners. But since 2016, the number of people infected through the heterosexual contacts exceeded the PWID cases and the epidemic has transmitted to the general population (Bobkov et al., 1998; Bobkova, 2013; Kazennova et al., 2014; Kotova et al., 2019; Schlösser et al., 2020).

Since 2008, our laboratory has extensively collected HIV-1 sequences in Siberian regions. Beginning in 2016, the HIV-1 surveillance studies have been conducted in the Central Asia countries as a part of the national HIV-1 drug resistance surveillance program. In this study, we used a comprehensive dataset of 872 pol gene HIV CRF63_02A6 sequences from Russia and countries of Central Asia to investigate the spatiotemporal dynamic of HIV-1 CRF63_02A6 sub-epidemic.

Materials and methods

Study population and HIV-1 sequence dataset

This study uses HIV-1partial pol gene sequences (1,047 base pairs, HXB2 #K03455 positions 2,252–3,299, n = 587) from our previous surveillance studies (Baryshev et al., 2012; Gashnikova et al., 2015, 2017; Ulyanova et al., 2019; Maksimenko et al., 2020; Sivay et al., 2021). Blood samples were collected from 2009 to 2020 from PWLH, both ART-naïve and therapy experienced, and were used for routine antiretroviral drug resistance testing. Samples were collected in six regions of Russia (Novosibirsk, Altai, Tomsk, Omsk, Krasnoyarsk, and Kemerovo). Further manipulations with the RNA/DNA samples as well as the population sequencing procedure were previously described (Sivay et al., 2021). Additional pol gene sequences were retrieved from Los Alamos National Laboratory (LANL) HIV Sequence Database (Los Alamos National Laboratory HIV Sequence Database, 2022a)²; all sequences classified as CRF63_02A6, CRF63_01A, CRF02A6, CRF02_AG, and “other” subtypes from FSU countries were included. A limited number of CRF02_AG sequences from African and European countries were also included in the dataset; those sequences were selected randomly.

HIV-1 subtyping

HIV-1 subtypes were identified by phylogenetic analysis. Study sequences were combined with the 2020 list of LANL HIV-1 subtype reference sequences (including all pure subtypes and CRFs) and the further phylogenetic analysis was performed with IQ-Tree (Nguyen et al., 2015) using maximum likelihood method and GRT + G4 + I substitution model and Shimodaira-Hasegawa approximate likelihood-ratio test with 1,000 replicates. HIV-1 subtype was assigned if study sequence falls into monophyletic clade with the corresponding subtype reference sequence with the branch support value ≥90. HIV-1 inter- and intra-subtype recombinants were identified using GARD (Kosakovsky Pond et al., 2006). In addition to the phylogenetic analysis, we applied the two widely used automated HIV-1 subtyping tools -- Stanford HIVdb program (Liu and Shafer, 2006) and COMET (Struck et al., 2014). The other widely used HIV-1 automated subtyping tool – REGA – was not used in this study; list of reference sequences that is used in REGA is outdated and does not include CRF63_02A6 sequences.

HIV-1 phylogenetic and transmission cluster analyses

Maximum likelihood phylogenetic tree was constructed using IQ-Tree under GTR + 4G + I substitution model as selected by jModelTest2 (Darriba et al., 2012). Transmission clusters were identified using Cluster Picker v1.2.3 (Ragonnet-Cronin et al., 2013) with the branch support threshold of 90 and maximum genetic distance threshold of 0.03.

Estimation of HIV-1 CRF63_02A6 origin time and location

To estimate ancestral state reconstruction and identify time to the most recent common ancestor (tMRCA) of the CRF63_02A6, a local instance of Nextstrain (Hadfield et al., 2018) was used. This is a recently developed software for real-time tracking of pathogens evolution such as COVID-19. To estimate the evolutionary rate and examine demographic history of the HIV-1 CRF63_02A6, a Bayesian Markov Chain Monte Carlo (MCMC) approach incorporated into BEAST v1.10.4 (Suchard et al., 2018) was used. To reduce computational complexity, a subsampling step was introduced, so that only a single sequence per transmission cluster was included in the further spatiotemporal analysis. The GTR + 4G substitution model, the lognormal uncorrelated relaxed clock model, and the Bayesian skyline coalescent model were used. Two independent MCMC runs of 15×10⁷ simulations were performed. Convergence of the MCMC results was examined in Tracer v1.7.1 (Rambaut et al., 2018) with effective sampling size (ESS) >200 for all parameters considered acceptable. The maximum clade credibility (MCC) tree was generated using TreeAnnotator v1.10.4 and visualized in FigTree v1.4.3 (Rambaut, 2009). We also estimated tMRCA and investigated geographic distribution of the CRF63_02A6 using BEAST v1.10.4 and compared it with the results obtained by the Nextstrain. Migration events of the HIV-1 CRF63_02A6 were reconstructed using the discrete Bayesian phylogeographic approach combined with the Bayesian Stochastic Search Variable Selection (BSSVS). Geographic location of sequences was assigned as a sampling city; if sampling city data was unavailable, then the country’s largest city location was used. Viral migration routes were summarized using SPREAD3 v0.9.7.1 (Bielejec et al., 2016).

Results

HIV-1 sequences dataset

A total of 2,637 HIV-1 pol gene sequences were included in the initial dataset using the described criteria and were used in phylogenetic analysis. Out of those sequences, 872 sequences were identified as HIV-1 CRF63_02A6, 1,157 sequences were identified as CRF02_AG_FSU, 412 sequences were from CRF02_AG_AF/EU clade, and 196 sequences were classified as unique recombinant form (URF) (Supplementary Figure S1). We compare HIV-1 subtyping results obtained from the phylogenetic analysis with those from Stanford and COMET automated HIV-1 subtyping tools. Out of the sequences identified as CRF63_02A6 by phylogenetic analysis, subtyping results by Stanford and COMET were concordant for 660 (75.7%) and 696 (79.8%) sequences, respectively. Moreover, out of the full sequence dataset, COMET classified 1,092 (41.4%) sequences as CRF63_02A6; 37.5% (409/1,092) of those sequences were classified as CRF_02AG_FSU by phylogenetic analysis. The detailed information about subtyping results is presented in Supplementary Figure S2.

Phylogenetic and transmission cluster analyses

Phylogenetic reconstruction was performed for 872 HIV-1 CRF63_02A6 pol gene sequences collected from 2009 to 2020 (Figure 2): 587 sequences were obtained by our research team and additional 285 sequences were retrieved form LANL HIV sequence database. Monophyletic structure of the CRF63_02A6 phylogenetic tree most likely indicates the initiation of the CRF63_02A6 sub-epidemic by a single ancestral virus. Predominant number (n = 832; 95.4%) of those sequences were from 25 regions of Russia; the remaining 40 (4.6%) sequences were from countries of Central Asia. A closer investigation ofHIV-1 CRF63_02A6 tree showed a total intermix of the sequences from different countries and regions across Russia supporting extensive virus distribution between those territories. Three well-supported (branch support value >95) location-specific monophyletic clades were detected: T/K clade was predominantly represented by CRF63_02A6 sequences from Tomsk and Krasnoyarsk regions (82% of the sequences within the clade), T clade was predominantly represented by sequences from Tomsk region (87% of the sequences within the clade), and JA/A clade was presented by sequences from Jewish Autonomous and Amur regions only. Those clades most likely represent separate introductions causing the super-spreading events of the virus. Recombination analysis revealed some evidence of recombination between CRF63_02A6 and sub-subtype A6 among sequences of the JA/A clade (n = 35, sequences dated 2018 and 2019); this could indicate the further local virus evolution.

FIGURE 2

Figure 2. IQ-Tree of HIV-1 CRF63_02A 872 pol gene sequences. Tree branches color-codded according to the country of isolation (top left chart). Transmission clusters were identified using Cluster Picker v1.2.3 using a branch support value ≥90 and a maximum genetic distance threshold of 0.03. Transmission clusters of ≥10 sequences are colored blue; transmission clusters of <10 sequences are colored yellow.

Out of 872 sequences, 519 (59.5%) sequences were found in 157 transmission clusters (using thresholds of 0.03 of genetic distance and 90 of branch support) of size 2–38 sequences per cluster; Seven large (≥10 sequences) clusters were identified. The remaining 91 clusters were dyads, 38 clusters were triplets, 8 clusters were the size of four, 5 cluster were the size of five, 3 clusters were the size of six, two clusters were the size of seven, one cluster was the size of eight, and two clusters were the size of nine. The detailed information about large transmission clusters is provided in Table 1.

TABLE 1

Table 1. Characteristics of the putative CRF63_02A6 transmission clusters.

Spatiotemporal dynamics of HIV-1 CRF63_02A6

Phylogenetic and evolutionary analysis of CRF63_02A6 872 pol gene sequences using Nextstrain identified that CRF63_02A6 originated in Novosibirsk around February 2005 (confidence interval: December, 2004 – July, 2005; Figure 3). The virus was further transmitted across other Siberian territories (Omsk, Altai, Tomsk, Krasnoyarsk, and Kemerovo). By 2015, CRF63_02A6 was identified in 13 regions of Russia, and by 2020 the virus was also detected in other nine Russian regions and six countries of Central Asia.

FIGURE 3

Figure 3. Nextstrain maximum likelihood analysis of HIV-1 CRF63_02A6 872 pol gene sequences. (A) RAx-ML tree was constructed using 832 sequences from Russia and 40 sequences from countries of Central Asia. Tree branches are colored according to the sampling city. (B) Geographical representation of CRF63_02A spatial dynamics. Lines between cities represent putative viral transitions; circle sizes are proportional to the square root of the number of sequences from the same location.

Phylodynamic reconstruction of HIV-1 CRF63_02A6

To investigate the evolutionary rate and population dynamics of the CRF63_02A6 using Bayesian phylogenetic inference, 510 pol gene sequences were included in the analysis. This number of sequences was defined using subsampling strategy described in Methods. The median estimated substitution rate was 2.3×10⁻³ substitutions/year (95% highest posterior density [HPD]: 1.9×10⁻³-2.7×10⁻³). The skyline plot shows an exponential growth of the ESS of the HIV-1 CRF63_02A6 from 2005 till 2010 followed by the stable phase (Figure 4). TMRCA was estimated as 2004 (95% highest posterior density: 2002–2006) using BEAST. Both results from the Nextstrain and Bayesian inference agreed and placed Novosibirsk as the most probable ancestral location of the HIV-1 CRF63_02A6 sub-epidemic. The tMRCA obtained by those two methods was also similar (2005 vs. 2004). Bayesian maximum clade credibility (MCC) tree and transmission history of the HIV-1 CRF63_02A6 constructed using SPREAD3 are presented in Supplementary Figure S3.

FIGURE 4

Figure 4. The Bayesian skyline plot of the 510 CRF63_02A6 pol gene sequences. The solid blue line indicates the median effective population size through time; the shaded blue area indicates the 95% highest posterior density (HPD) interval of the effective population size. The vertical solid black line indicates the median (2004) estimated tMRCA of the HIV-1 CRF63_02A6; the vertical dashed lines indicate the estimated tMRCA 95% HPD interval (2002–2006).

Discussion

This report represents the characterization of the spatiotemporal dynamics of HIV-1 CRF63_02A6 sub-epidemic. This analysis includes 872 CRF63_02A6 pol gene sequences sampled in seven Eurasian countries. By the time the study was conducted, CRF63_02A6 sequences were sampled predominantly in Russia, along with the sporadic cases detected in six countries of Central Asia. The current analysis revealed that the onset date of HIV-1 CRF63_02A6 sub-epidemic was 2005 and the most probable origin location was Novosibirsk; the virus was subsequently transmitted to other regions of Russia and was also detected in Kyrgyzstan, Uzbekistan, Tajikistan, Azerbaijan, Armenia, and Georgia.

The HIV-1 genetic diversity in Russia is broad and characterized by the region-specific sub-epidemics. The first HIV-1 infection outbreak was registered in 1988–1989 among nosocomially infected children in the South of Russia and was caused by the HIV-1 subtype G (Thomson and Najera, 2007; Serova et al., 2018), but this subtype is currently rarely detected in the country. In the mid-1990s, the multiple HIV-1 outbreaks among PWID were detected in Central Russia (Kazennova et al., 2014) and gave rise to an extensive expansion of the HIV-1 epidemic in the country caused by the sub-subtype A6. In the following years, the virus was detected in different Russian regions and neighboring countries (Thomson et al., 2009) and currently predominates across all FSU countries. Subtype B was identified among MSM population and presently accounts for up to 2–5% of all HIV-1 cases (Gashnikova et al., 2015; Pasechnik and Blokh, 2018; Lebedev et al., 2019). Recombinant form of subtype B and sub-subtype A6 - CRF03_AB caused several outbreaks among PWID in the Northwestern Russia (Lebedev et al., 2020) and predominates alongside with the sub-subtype A6 (Ozhmegova et al., 2020). HIV-1 genetic diversity in the Russian Far East is different from the rest of the country. Besides the predominant sub-subtype A6, subtype B and subtype C were found to account for 25% and 10% of the HIV-1 cases, respectively (Kazennova et al., 2008).

Up to mid-2000s, sub-subtype A6 predominated in Siberia. But since 2008, a rapid spread of novel HIV-1 genetic variant was detected in Novosibirsk (Gashnikova et al., 2011). Detailed full genome analysis identified this genetic variant as a recombinant of sub-subtype A6 and CRF02_AG_FSU; phylogenetic analysis revealed that the new recombinant form (subsequently named as CRF63_02A6) was genetically related to the CRF02_AG_FSU from Uzbekistan (Baryshev et al., 2012). A retrospective molecular-epidemiological analysis revealed that some of the PWLH in Novosibirsk were HIV-infected by this recombinant form back in 2004–2006, when the local outbreaks among PWID were registered (Gashnikova et al., 2011). The HIV-1 CRF63_02A6 sub-epidemic distribution in Russia continued and the number of CRF63_02A6 cases started to increase in other regions (Safina et al., 2022). CRF63_02A6 was also detected in the countries of Central Asia, however those cases are sporadic and most likely represent cross-border transmission events (Sivay et al., 2021). In the current study, a large-scale phylogenetic analysis of over 2,500 pol gene sequences showed that CRF63_02A6 forms a well-supported distinct monophyletic clade, which most likely indicates the initiation of the CRF63_02A6 epidemic by a single ancestral virus. By now, CRF63_02A6 was detected in over 20 Russian territories and six countries of Central Asia. The current phylogenetic analysis of the CRF63_02A6 revealed three distinct region-specific clades, most likely indicating the local super-spreading events.

Earlier studies (Kostaki et al., 2018) on demographic history and phylodynamic analysis of the CRF02_AG/CRF63_02A6 combined CRF63_02A6 and CRF02_AG_FSU lineages. The study identified the origin of the CRF63_02A6 in Kazakhstan and Uzbekistan in 1996; two Russian subclades dated 2003 and 2007 were detected. By the time that study was conducted, a very limited number (n = 136) of the CRF63_02A6 sequences were available. Our study and the study of the global dissemination of the CRF02_AG (Mir et al., 2016) showed that CRF63_02A6 and CRF02_AG_FSU represent two distinct HIV-1 genetic variants. Given the results of our previous study and earlier report by Baryshev et al. (2012), CRF63_02A6 most likely descended from CRF02_AG_FSU and is characterized by a significant genome change compared to CRF02_AG_FSU. Both of those CRFs have different geographic characteristics. Thus, CRF_02AG_FSU is endemic in Central Asia (Sivay et al., 2021), while CRF63_02A6 predominantly circulates in Russia. Current report includes all the available HIV-1 CRF63_02A6 and CRF02_AG_FSU sequences published in the LANL HIV sequence database from the Eurasian continent (Supplementary Figure S1). Our analysis shows that all the pol gene sequences from Kazakhstan available in LANL HIV Sequences Database belong to CRF02_AG_FSU. Our analysis also finds only a few CRF63_02A6 sequences from Uzbekistan dated 2015; most of the HIV-1 sequences from Uzbekistan are CRF02_AG_FSU.

HIV-1 epidemic in countries of Central Asia was driven by the drug trafficking through the so called “Northern route” from Afghanistan to Russia through Tajikistan, Uzbekistan, Kirgizia, and Kazakhstan (Vorobyev, 2014; Merz, 2018). Since earlier-2000s, the rapid increase of drug use was registered in Russia, and from 2002 to 2014 Siberia was leading by the number of PWID in the country (Andreev and Kokunova, 2017). Also, an intensive labor migration from Central Asian countries to Novosibirsk was registered. Those factors could create favorable conditions for HIV-1 transmission. According to Andreev and Kokunova (2017), a high level of drug use was found among HIV-infected illegal labor migrants. Rapid CRF63_02A6 sub-epidemic expansion in Siberia could also be linked to the high synthetic drug use starting in 2009. Since 2011, the decline in opioid drugs use together with an increase in synthetic drugs use were registered (Vorobyev, 2014; Andreev and Kokunova, 2017). Synthetic drugs induce both stimulant and sexual effects promoting high-risk sexual behavior (Luo et al., 2018). Active synthetic drug use could contribute to dissemination of distinct CRF63_02A6 lineages resulting to a fold-increase of new HIV-1 cases in Siberian regions (Gashnikova et al., 2015, 2017). The well-supported region-specific monophyletic CRF63_02A6 clades identified in the current study most likely indicate these super-spreading events (Figure 2).

Since 2005, the number of new HIV-1 cases in Novosibirsk has started to grow exponentially. The largest increase (from 19.6% in 2003 to 57% in 2007) in new HIV-1 cases in the region was recorded among prisoners (Novosibirsk Region Center on AIDS, 2014).³ Most likely such increase was driven by the high-risk behavior, such as intravenous drug use (IDU). Most people who entered the Russian penal system in mid-2000s were PWID (Severson, 2011). The penal system had transformed into a unique risk environment for PWID and HIV-1 transmission due to the lack of harm reduction programs. An earlier study among prisoners in Russia showed that almost a half of them had a history of drug injection; that 20% of prisoners had used drugs while imprisoned, 64% of those individuals had injected drugs with previously used injecting equipment, and that 13.5% had started injecting drugs while imprisoned (Jürgens and Betteridge, 2005). Several HIV-1 outbreaks associated with IDU were registered among prisoners in Russia (Bobrik et al., 2005; Severson, 2011). We admit that the emergence of the HIV-1 CRF63_02A6 could have happened within one of the prisons in the Novosibirsk.

In Russia, the guidelines for laboratory testing of HIV-1 infection requires a fourth-generation HIV-EIA test (4th Gen), which detects both antibody and antigen with the following confirmatory Western Blot testing and/or HIV-1 RNA testing (Ministry of Health of Russian Federation, 2020).⁴ However, HIV-1 RNA testing is not routinely used and is generally provided to those individuals who have previously confirmed HIV-1 infection, as well as children and pregnant women. Patients with the acute HIV-1 infection can be misdiagnosed during the HIV-1 seronegative window period and can easily transmit the virus. Earlier report described clinical characteristics of patients with acute HIV-1 infection in Novosibirsk. Most of the patients (76%) had CRF63_02A6 infection and 19.4% had sub-subtype A6 infection. Statistically significant differences were found between clinical manifestation in patients with the HIV-1 CRF63_02A6 and sub-subtype A6 infection. Thus, CRF63_02A6 infection is characterized by more severe primary symptoms, such as prolonged fever, generalized lymphadenopathy, higher (>10⁷) viral load, and lower CD4 cell counts compared to patients with the sub-subtype A6 infection (Ulyanova et al., 2019). These data indicate a faster disease progression and rapid virus transmission among patients with the CRF63_02A6 compared to sub-subtype A6 infection. Symptoms of acute HIV-infection are nonspecific, and the infection often can be misdiagnosed, especially within the ongoing COVID-19 pandemic and COVID-19 mass vaccination. Beginning in 2000, our research team initiated HIV-1 surveillance program in Novosibirsk, which was also launched in other regions of Siberia around 2010. Our studies detected the displacement of sub-subtype A6 by CRF63_02A6, and started in 2009 CRF63_02A6 has been detected in over 80% of new HIV-1 cases in Siberia. Also, various URFs of sub-subtype A6 and CRF63_02A6 has been persistently detected (Gashnikova et al., 2015, 2017; Ulyanova et al., 2019).

Current analysis identified a misclassification of CRF63_02A6 as CRF02_AG by widely used HIV-1 automated subtyping tools. HIV-1 has a high recombination rate promoting viral genetic diversification and escape from the host immune system (Levy et al., 2004). New CRFs are constantly detecting and 118 of them have been described by now (Los Alamos National Laboratory HIV Sequence Database, 2022b).⁵ Proper identification of HIV-1 subtypes is challenging due to the frequent recombination. Molecular phylogeny remains the most reliable HIV-1 subtyping approach (Fabeni et al., 2017), however this method requires the data manipulation and interpretation skills. In the clinical settings, HIV-1 subtyping is routinely based on fast and simple HIV-1 automated subtyping tools which have some limitations compared to the phylogenetic approach. Automated subtyping tools often have a very limited number of CRFs in the reference dataset. Most of the automated HIV-1 subtyping algorithms are designed mainly for the most prevalent viral subtypes such as B and C, and CRFs such as CRF01_AE and CRF02_AG (Fabeni et al., 2017). Also, those tools are not frequently updated. In the current report, we used two automated HIV-1 subtyping tools, Stanford and COMET, and maximum likelihood phylogenetic approach (combined with 2020 LANL HIV-1 subtype reference sequences) to identify CRF63_02A6 sequences. The HIV-1 CRFs were assigned based on high phylogenetic relatedness of the study sequences with the corresponding HIV-1 subtype reference sequences. HIV-1 subtyping results of CRF63_02A6 sequences obtained by Stanford and COMET automated subtyping tools were concordant with the phylogenetic analysis results in approximately 75 and 79% instances, respectively. Moreover, a third of CRF02_AG_FSU sequences were classified as CRF63_02A6 by COMET. Our data highlight a rather limited performance of Stanford and COMET subtyping tools in CRF63_02A6 assignment compared to the phylogenetic approach. Although the automated HIV-1 subtyping tools are widely used in clinical settings, the reference datasets of these tools need to be updated to improve reliability of the HIV-1 subtyping results.

This report provides a comprehensive analysis on the origin and dissemination of the recent HIV-1 CRF63_02A sub-epidemic in Russia and countries of Central Asia. The study describes the onset of the CRF63_02A6 in Siberia in the mid-2000s and its following spreading across Russia and the neighboring countries. The CRF63_02A6 distribution simultaneously occurred with the use of injectable and synthetic drugs. Since 2016, the HIV-1 heterosexual transmission in Siberia has prevailed over PWID, however HIV-1 incidence rate among PWID remains high. A high prevalence (from 2 to 30% across the regions) of HIV-1 URFs indicates a substantial level of co-infection with different HIV-1 subtypes. Further surveillance studies of the HIV-1 CRF63_02A6 are needed due to the continuing rapid distribution and high rate of recombination between CRF63_02A6 and other HIV-1 subtypes/CRFs. Studies of epidemiological characteristics and their relationships with the new emerging HIV-1 genetic variants are urgently needed. Also, the introduction of HIV-1 prevention and treatment programs are required to implement further HIV-1 transmission reduction strategies. This study also emphasizes importance of updating of the reference datasets of the widely used automated HIV-1 subtyping tools to provide reliable HIV-1 subtyping results.

Data availability statement

The data presented in this study are deposited in the GenBank repository, accession numbers are listed in Supplementary Data.

Ethics statement

The studies involving human participants were reviewed and approved by Local Ethical Committee of the State Budgetary Healthcare Institution of Novosibirsk Region “City Infectious Clinical Hospital #1”. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author contributions

MS, AT, and NG planned and designed the study. LM, IO, AN, MG, DZ, VE, TN, and VI performed HIV genotyping and collected sequences. MS performed the analysis of the epidemiological data, performed the phylogenetic and phylodynamic analysis, produced the illustrations, and wrote the manuscript. NG supervised the project and edited the manuscript. DK, LP, SS, TO, SY, ON, AC, and TI collected the epidemiological data and assisted with samples collection at the local study clinics. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the State Research and Development Program #6–21 State Research Center of Virology and Biotechnology “Vector.”

Acknowledgments

The authors want to thank all of the patient for their contributions to the study. We also thank the staff at the study HIV/AIDS clinics for laboratory testing.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.946787/full#supplementary-material

Footnotes

1. ^http://www.hivrussia.info/wp-content/uploads/2022/03/Spravka-VICH-v-Rossii-na-31.12.2021-g.pdf

2. ^http://www.hiv.lanl.gov/

3. ^https://oblmed.nsk.ru/upload/files/hiv.pdf

4. ^http://rushiv.ru/category/docs/national-recs/kr-20/

5. ^https://www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html

References

Aibekova, L., Foley, B., Hortelano, G., Raees, M., Abdraimov, S., Toichuev, R., et al. (2018). Molecular epidemiology of HIV-1 subtype A in former Soviet Union countries. PLoS One 13:e0191891. doi: 10.1371/journal.pone.0191891

PubMed Abstract | CrossRef Full Text | Google Scholar

Andreev, A. P., and Kokunova, S. D. (2017). Effect of illegal migration on drug use in Russian Federation. Counter-Terrorism 4, 19–27.

Google Scholar

Baryshev, P. B., Bogachev, V. V., and Gashnikova, N. M. (2012). Genetic characterization of an isolate of HIV type 1 AG recombinant form circulating in Siberia. Russia. Arch. Virol. 157, 2335–2341. doi: 10.1007/s00705-012-1442-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Baryshev, P. B., Bogachev, V. V., and Gashnikova, N. M. (2014). HIV-1 genetic diversity in Russia: CRF63_02A1, a new HIV type 1 genetic variant spreading in Siberia. AIDS Res. Hum. Retrovir. 30, 592–597. doi: 10.1089/aid.2013.0196

PubMed Abstract | CrossRef Full Text | Google Scholar

Bielejec, F., Baele, G., Vrancken, B., Suchard, M. A., Rambaut, A., and Lemey, P. (2016). SpreaD3: interactive visualization of spatiotemporal history and trait evolutionary processes. Mol. Biol. Evol. 33, 2167–2169. doi: 10.1093/molbev/msw082

PubMed Abstract | CrossRef Full Text | Google Scholar

Bobkov, A., Kazennova, E., Selimova, L., Bobkova, M., Khanina, T., Ladnaya, N., et al. (1998). A sudden epidemic of HIV type 1 among injecting drug users in the former Soviet Union: identification of subtype A, subtype B, and novel gag A/env B recombinants. AIDS Res. Hum. Retrovir. 14, 669–676. doi: 10.1089/aid.1998.14.669

PubMed Abstract | CrossRef Full Text | Google Scholar

Bobkov, A. F., Kazennova, E. V., Selimova, L. M., Khanina, T. A., Ryabov, G. S., Bobkova, M. R., et al. (2004). Temporal trends in the HIV-1 epidemic in Russia: predominance of subtype A. J. Med. Virol. 74, 191–196. doi: 10.1002/jmv.20177

PubMed Abstract | CrossRef Full Text | Google Scholar

Bobkova, M. (2013). Current status of HIV-1 diversity and drug resistance monitoring in the former USSR. AIDS Rev. 15, 204–212.

Google Scholar

Bobrik, A., Danishevski, K., Eroshina, K., and McKee, M. (2005). Prison health in Russia: the larger picture. J. Public Health Policy 26, 30–59. doi: 10.1057/palgrave.jphp.3200002

PubMed Abstract | CrossRef Full Text | Google Scholar

Carr, J. K., Nadai, Y., Eyzaguirre, L., Saad, M. D., Khakimov, M. M., Yakubov, S. K., et al. (2005). Outbreak of a west African recombinant of HIV-1 in Tashkent, Uzbekistan. JAIDS J. Acquired Immune Deficiency Syndromes 39, 570–575.

Google Scholar

Central Research Institute of Epidemiology, Federal Service on Customers' Rights Protection and Human Well-being Surveillance (2021). HIV statistics in Russian Federation. http://www.hivrussia.info/wp-content/uploads/2022/03/Spravka-VICH-v-Rossii-na-31.12.2021-g.pdf (Accessed April 5, 2022).

Google Scholar

Darriba, D., Taboada, G. L., Doallo, R., and Posada, D. (2012). jModelTest 2: more models, new heuristics and parallel computing. Nat. Methods 9:772. doi: 10.1038/nmeth.2109

PubMed Abstract | CrossRef Full Text | Google Scholar

Diez-Fuertes, F., Cabello, M., and Thomson, M. M. (2015). Bayesian phylogeographic analyses clarify the origin of the HIV-1 subtype A variant circulating in former Soviet Union’s countries. Infect. Genet. Evol. 33, 197–205. doi: 10.1016/j.meegid.2015.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Fabeni, L., Berno, G., Fokam, J., Bertoli, A., Alteri, C., Gori, C., et al. (2017). Comparative evaluation of subtyping tools for surveillance of newly emerging HIV-1 strains. J. Clin. Microbiol. 55, 2827–2837. doi: 10.1128/JCM.00656-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Frank, T. D., Carter, A., Jahagirdar, D., Biehl, M. H., Douwes-Schultz, D., Larson, S. L., et al. (2019). Global, regional, and national incidence, prevalence, and mortality of HIV, 1980–2017, and forecasts to 2030, for 195 countries and territories: a systematic analysis for the global burden of diseases, injuries, and risk factors study 2017. Lancet HIV 6, e831–e859. doi: 10.1016/S2352-3018(19)30196-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Gashnikova, N. M., Bogachev, V. V., Baryshev, P. B., Totmenin, A. V., Gashnikova, M. P., Kazachinskaya, A. G., et al. (2015). A rapid expansion of HIV-1 CRF63_02A1 among newly diagnosed HIV-infected individuals in the Tomsk region. Russia. AIDS Res. Hum. Retroviruses 31, 456–460. doi: 10.1089/aid.2014.0375

PubMed Abstract | CrossRef Full Text | Google Scholar

Gashnikova, N. M., Zyryanova, D. P., Astakhova, E. M., Ivlev, V. V., Gashnikova, M. P., Moskaleva, N. V., et al. (2017). Predominance of CRF63_02A1 and multiple patterns of unique recombinant forms of CRF63_A1 among individuals with newly diagnosed HIV-1 infection in Kemerovo oblast. Russia. Arch. Virol. 162, 379–390. doi: 10.1007/s00705-016-3120-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Gashnikova, N. M., Safronov, P. F., Nikonorova, Y. V., Unagaeva, N. V., Lapteva, T. A., Bogachev, V. V., et al. (2011). Properties of CRF02_AG HIV-1 isolates circulating in Novosibirsk region. Zh. Mikrobiol. Epidemiol. Immunobiol. 3, 38–43.

Google Scholar

Hadfield, J., Megill, C., Bell, S. M., Huddleston, J., Potter, B., Callender, C., et al. (2018). Nextstrain: real-time tracking of pathogen evolution. Bioinformatics 34, 4121–4123. doi: 10.1093/bioinformatics/bty407

PubMed Abstract | CrossRef Full Text | Google Scholar

Jürgens, R., and Betteridge, G. (2005). Prisoners who inject drugs: public health and human rights imperatives. Emerg. Issues in HIV/AIDS 8, 46–74. doi: 10.2307/4065334

CrossRef Full Text | Google Scholar

Karamov, E., Epremyan, K., Siniavin, A., Zhernov, Y., Cuevas, M. T., Delgado, E., et al. (2018). HIV-1 genetic diversity in recently diagnosed infections in Moscow: predominance of A_FSU, frequent branching in clusters, and circulation of the Iberian subtype G variant. AIDS Res. Hum. Retrovir. 34, 629–634. doi: 10.1089/aid.2018.0055

PubMed Abstract | CrossRef Full Text | Google Scholar

Kazennova, E. V., Bronnikova, A. V., Kuzin, S. N., Kirillova, I. L., Ershova, O. N., Afanaseva, L. R., et al. (2008). Molecular genetic characteristics of HIV-1 variants circulating in Cherepovets, Vologda region: the second case of the epidemic outbreak caused by the recombinant gagAenvB. Voprosi Virusologii. 53, 23–27.

PubMed Abstract | Google Scholar

Kazennova, E., Laga, V., Lapovok, I., Glushchenko, N., Neshumaev, D., Vasilyev, A., et al. (2014). HIV-1 genetic variants in the Russian Far East. AIDS Res. Hum. Retrovir. 30, 742–752. doi: 10.1089/aid.2013.0194

PubMed Abstract | CrossRef Full Text | Google Scholar

Kazennova, E. V., Antonova, A. A., Ozhmegova, E. N., Demyanenko, E. R., Minakova, M. V., Belousova, O. V., et al. (2020). Genetic analysis of HIV-1 in the Altai Kray: The further spread of the CRF63_02A1 variant in Western Siberia. HIV Infect. Immunosuppressive Disord. 12, 47–57. doi: 10.22328/2077-9828-2020-12-1-47-57

CrossRef Full Text | Google Scholar

Kotova, V. O., Trotsenko, O. E., Balakhontseva, L. A., and Bazykina, E. A. (2019). Molecular genetic characteristics of HIV-1 variants isolated in the subjects of the Russian Far East. Vopr. Virusol. 64, 79–89. doi: 10.18821/0507-4088-2019-64-2-79-89

PubMed Abstract | CrossRef Full Text | Google Scholar

Kosakovsky Pond, S. L., Posada, D., Gravenor, M. B., Woelk, C. H., and Frost, S. D. (2006). GARD: a genetic algorithm for recombination detection. Bioinformatics 22, 3096–3098. doi: 10.1093/bioinformatics/btl474

PubMed Abstract | CrossRef Full Text | Google Scholar

Kostaki, E. G., Karamitros, T., Bobkova, M., Oikonomopoulou, M., Magiorkinis, G., Garcia, F., et al. (2018). Spatiotemporal characteristics of the HIV-1 CRF02_AG/CRF63_02A1 epidemic in Russia and Central Asia. AIDS Res. Hum. Retrovir. 34, 415–420. doi: 10.1089/aid.2017.0233

PubMed Abstract | CrossRef Full Text | Google Scholar

Lebedev, A., Lebedeva, N., Moskaleychik, F., Pronin, A., Kazennova, E., and Bobkova, M. (2019). Human immunodeficiency virus-1 diversity in the Moscow region, Russia: Phylodynamics of the most common subtypes. Front. Microbiol. 10:320. doi: 10.3389/fmicb.2019.00320

PubMed Abstract | CrossRef Full Text | Google Scholar

Lebedev, A., Pasechnik, O., Ozhmegova, E., Antonova, A., Blokh, A., Grezina, L., et al. (2020). Prevalence and spatiotemporal dynamics of HIV-1 circulating recombinant form 03_AB (CRF03_AB) in the former Soviet Union countries. PLoS One 15:e0241269. doi: 10.1371/journal.pone.0241269

PubMed Abstract | CrossRef Full Text | Google Scholar

Levy, D. N., Aldrovandi, G. M., Kutsch, O., and Shaw, G. M. (2004). Dynamics of HIV-1 recombination in its natural target cells. Proc. Natl. Acad. Sci. 101, 4204–4209. doi: 10.1073/pnas.0306764101

PubMed Abstract | CrossRef Full Text | Google Scholar

Liitsola, K., Holm, K., Bobkov, A., Pokrovsky, V., Smolskaya, T., Leinikki, P., et al. (2000). An AB recombinant and its parental HIV type 1 strains in the area of the former Soviet Union: low requirements for sequence identity in recombination. AIDS Res. Hum. Retrovir. 16, 1047–1053. doi: 10.1089/08892220050075309

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, T. F., and Shafer, R. W. (2006). Web resources for HIV type 1 genotypic-resistance test interpretation. Clin. Infect. Dis. 42, 1608–1618. doi: 10.1086/503914

PubMed Abstract | CrossRef Full Text | Google Scholar

Los Alamos National Laboratory HIV Sequence Database. (2022a). http://www.hiv.lanl.gov/ (Accessed January 15, 2022).

Google Scholar

Los Alamos National Laboratory HIV Sequence Database (2022b). https://www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html (Accessed April 6, 20229).

Google Scholar

Luo, W., Hong, H., Wang, X., McGoogan, J. M., Rou, K., and Wu, Z. (2018). Synthetic drug use and HIV infection among men who have sex with men in China: a sixteen-city, cross-sectional survey. PLoS One 13:e0200816. doi: 10.1371/journal.pone.0200816

PubMed Abstract | CrossRef Full Text | Google Scholar

Maksimenko, L. V., Totmenin, A. V., Gashnikova, M. P., Astakhova, E. M., Skudarnov, S. E., Ostapova, T. S., et al. (2020). Genetic diversity of HIV-1 in Krasnoyarsk Krai: area with high levels of HIV-1 recombination in Russia. BioMed Research international, 2020. doi: 10.1155/2020/9057541

CrossRef Full Text | Google Scholar

Merz, F. (2018). United Nations Office on drugs and crime: world drug report 2017. SIRIUS-Zeitschrift für Strategische Analysen 2, 85–86. doi: 10.1515/sirius-2018-0016

CrossRef Full Text | Google Scholar

Ministry of Health of Russian Federation (2020). Clinical recommendations. HIV infection in adults and children. http://rushiv.ru/category/docs/national-recs/kr-20/

Google Scholar

Mir, D., Jung, M., Delatorre, E., Vidal, N., Peeters, M., and Bello, G. (2016). Phylodynamics of the major HIV-1 CRF02_AG African lineages and its global dissemination. Infect. Genet. Evol. 46, 190–199. doi: 10.1016/j.meegid.2016.05.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, L. T., Schmidt, H. A., Von Haeseler, A., and Minh, B. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274. doi: 10.1093/molbev/msu300

PubMed Abstract | CrossRef Full Text | Google Scholar

Novosibirsk Region Center on AIDS (2014). HIV-Infection. Incidence. HIV pre- and post-Testing Counseling. HIV Testing among high risk Groups. https://oblmed.nsk.ru/upload/files/hiv.pdf

Google Scholar

Ozhmegova, E. N., Antonova, A. A., Lebedev, A. V., Melnikova, T. N., Krylova, T. V., Kazachek, A. V., et al. (2020). Genetic profile of HIV-1 in the Vologda region: domination of CRF03_AB and rapid distribution of URFs. HIV Infect. Immunosuppressive Disord. (In Russ.) 12, 79–88. doi: 10.22328/2077-9828-2020-12-2-79-88

CrossRef Full Text | Google Scholar

Pasechnik, O. A., and Blokh, A. I. (2018). The prevalence of HIV recombinant forms in Russia and countries of the CIS: systematic review and metaanalysis. Russ. J. Infect. Immun. 8, 127–138. doi: 10.15789/2220-7619-2018-2-127-138

CrossRef Full Text | Google Scholar

Ragonnet-Cronin, M., Hodcroft, E., Hué, S., Fearnhill, E., Delpech, V., Brown, A. J. L., et al. (2013). Automated analysis of phylogenetic clusters. BMC Bioinform. 14, 1–10. doi: 10.1186/1471-2105-14-317

CrossRef Full Text | Google Scholar

Rambaut, A. (2009). FigTree. Tree Figure Drawing tool. http://tree.bio.ed.ac.uk/software/figtree/.

Google Scholar

Rambaut, A., Drummond, A. J., Xie, D., Baele, G., and Suchard, M. A. (2018). Posterior summarization in Bayesian phylogenetics using tracer 1.7. Syst. Biol. 67, 901–904. doi: 10.1093/sysbio/syy032

PubMed Abstract | CrossRef Full Text | Google Scholar

Safina, K. R., Sidorina, Y., Efendieva, N., Belonosova, E., Saleeva, D., Kirichenko, A., et al. (2022). Molecular epidemiology of HIV-1 in Oryol oblast, Russia. Virus Evolution. 8:1, veac044. doi: 10.1093/ve/veac044

CrossRef Full Text | Google Scholar

Schlösser, M., Kartashev, V. V., Mikkola, V. H., Shemshura, A., Saukhat, S., Kolpakov, D., et al. (2020). HIV-1 sub-subtype A6: settings for normalised identification and molecular epidemiology in the southern Federal District, Russia. Viruses 12:475. doi: 10.3390/v12040475

PubMed Abstract | CrossRef Full Text | Google Scholar

Serova, E., Murzakova, A., Nossik, M., Shemshura, A., Bazykin, G., and Kireev, D. (2018). The 1988 Elista outbreak and the epidemiology of HIV subtype G in Russia. Virus Evolution 4:vey010.010.

Google Scholar

Severson, J. L. (2011). HIV and AIDS in the Russian Federation: Prisons as a case Study of risk Environments and Agency [doctoral dissertation]. Austin (TX): University of Texas at Austin, 4. doi: 10.1093/ve/vey010.010

CrossRef Full Text | Google Scholar

Sivay, M. V., Totmenin, A. V., Zyryanova, D. P., Osipova, I. P., Nalimova, T. M., Gashnikova, M. P., et al. (2021). Characterization of HIV-1 epidemic in Kyrgyzstan. Front. Microbiol. 2985. doi: 10.3389/fmicb.2021.753675

CrossRef Full Text | Google Scholar

Struck, D., Lawyer, G., Ternes, A. M., Schmit, J. C., and Bercoff, D. P. (2014). COMET: adaptive context-based modeling for ultrafast HIV-1 subtype identification. Nucleic Acids Res. 42:e144. doi: 10.1093/nar/gku739

PubMed Abstract | CrossRef Full Text | Google Scholar

Suchard, M. A., Lemey, P., Baele, G., Ayres, D. L., Drummond, A. J., and Rambaut, A. (2018). Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 4:vey016. doi: 10.1093/ve/vey016

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomson, M. M., and Nájera, R. (2007). Increasing HIV-1 genetic diversity in Europe. J. Infect. Dis. 196, 1120–1124. doi: 10.1086/521683

CrossRef Full Text | Google Scholar

Thomson, M. M., Vinogradova, A., Delgado, E., Rakhmanova, A., Yakovlev, A., Cuevas, M. T., et al. (2009). Molecular epidemiology of HIV-1 in St Petersburg, Russia: predominance of subtype A, former Soviet Union variant, and identification of intrasubtype subclusters. JAIDS J. Acquired Immune Deficiency Syndromes 51, 332–339. doi: 10.1097/QAI.0b013e31819c1757

PubMed Abstract | CrossRef Full Text | Google Scholar

Ulyanova, Y. S., Gashnikova, N. M., Ivlev, V. V., Krasnova, E. I., Khokhlova, N. I., Totmenin, A. V., et al. (2019). Clinical and laboratory characterictic of acute HIV-infection in adult residents of Novosibirsk region. J. Infectol. 11, 40–44. doi: 10.22625/2072-6732-2019-11-2-40-44

CrossRef Full Text | Google Scholar

Vorobyev, M. V., Injection drug use in Russian Federation, 2005-2013 (2014). Current Problems of Science and Education, 5.

Google Scholar