Next Article in Journal
HIV/AIDS Mathematical Model of Triangle Transmission
Previous Article in Journal
Changes following the Onset of the COVID-19 Pandemic in the Burden of Hospitalization for Respiratory Syncytial Virus Acute Lower Respiratory Infection in Children under Two Years: A Retrospective Study from Croatia
Previous Article in Special Issue
Prevalence of Pretreatment HIV-1 Drug Resistance in Armenia in 2017–2018 and 2020–2021 following a WHO Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Institute of Microbiology and Immunology, Faculty of Medicine, University of Belgrade, 11000 Belgrade, Serbia
2
Institute for Oncology and Radiology of Serbia, 11000 Belgrade, Serbia
3
Laboratory of T-Lymphotropic Viruses, N.F. Gamaleya National Research Center of Epidemiology and Microbiology, 123098 Moscow, Russia
4
St. Petersburg City AIDS Center, 190103 St. Petersburg, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Viruses 2022, 14(12), 2748; https://doi.org/10.3390/v14122748
Submission received: 29 July 2022 / Revised: 29 November 2022 / Accepted: 29 November 2022 / Published: 9 December 2022

Abstract

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

1. Introduction

Russia is home to over 140 million people, and the site of the fastest growing HIV/AIDS epidemic in Europe. Newly diagnosed infections in 2020 in the Russian Federation contributed 57% of all new cases in the World Health Organization (WHO) European Region, with the highest rate of 40.8 new HIV diagnoses per 100,000 population [1]. According to the latest official data, since 1987 until December 2020 nearly 1.5 million cases of HIV infection have been registered, whereas over 1.1 million people in Russia currently live with HIV [2]. In contrast to the global epidemic pattern, the HIV epidemic in Russia continues to expand significantly [3].
Nevertheless, HIV epidemic burden is non-uniform in different Russian regions, with the major metropolitan cities, Ural and Siberia, disproportionately bearing the burden: the most affected regions are Urals (Yekaterinburg-107398, Tyumen-62414), the Volga region (Samara-79442) and Siberia (Irkutsk-64554, Kemerovo-77971), as well as the largest cities, Moscow (67913) and St. Petersburg (69123) cases. An explosive epidemic has been documented among people who inject drugs (PWID), starting from the mid-1990s, whereas presently, the majority of new infections are linked to sexual transmission [1,2]. Nationwide, among the cases with known infection risk, about 65% reported heterosexual contact, about 2% are men who have sex with men (MSM), and about 32% are PWID, while male to female ratio among newly diagnosed cases is approximately 60:40, pointing to heterosexual risk, MSM, and PWID as the key populations in the context of the Russian HIV epidemic.
Regarding HIV subtype distribution, nationwide epidemiological surveys and a number of local studies indicated a large predominance of the A6 sub-subtype virus (previously called AFSU), with the significant presence of other subtypes, namely subtype B and CRF063_02A [4,5,6,7,8,9,10,11,12]. The abovementioned big HIV epidemic among PWID has been shown to be the result of two single introductions of sub-subtype A6 and subtype B (labelled also BFSU or, IDU-B) [13,14].
The spread of subtype B seems to have gone through several independent viral lines and has circulated through sexual contacts, mainly among MSM at relatively low transmission rates, prior to its more significant presence in the overwhelming HIV epidemic in Russia [15,16,17].
The occurrence of HIV subtype B was considered rare in Russia; nevertheless, a recent analysis showed its significant presence in some regions, e.g., in the Russian Far East, probably linked to labor migrants [16]. Indeed, the population of HIV-1 subtype B strains in Russia and countries of the former Soviet Union has a complex structure. Besides numerous strains related to subtype B sequences circulating in Western Europe and seemingly originating from multiple independent introductions, most probably from Europe and the USA (“Western B”), there are at least two circulating subtype B viruses that are specific for this region. One of them is associated with the outbreak among PWID in Nikolaev, Ukraine, in 1996, namely IDU-B (= BFSU = FSU-B = BFSU) [14]. The initial origin of this strain has not been elucidated; it might be a recombinant of two or more subtype B variants. The circulation of subtype B variant related to the Brazilian subtype B has also been described, along with Thai subtype B’ virus [13].
In St. Petersburg, the first HIV-1 infections among PWIDs were diagnosed in the mid-1990s [18]. A number of studies so far have addressed subtype distribution and HIV molecular epidemiology in Russian’s second largest city, however, in depth phylogenetic studies are still scarce, especially regarding non-subtype A [19,20].
This study explores HIV subtype B sequences from St. Petersburg, collected from 2006 to 2020, in order to phylogenetically investigate and characterize transmission clusters, focusing on their evolutionary dynamics and potential for further growth, along with socio-demographic analysis of the available metadata.

2. Methods

2.1. Study Population

The study was approved by the “Health committee of the government of Saint Petersburg. Saint Petersburg state budget health institution Center for prevention and control of AIDS and infectious diseases”, Protocol No199, 2019. The study included HIV-1 subtype B infected, consenting patients, and referred to the St. Petersburg municipal AIDS center for genotypic drug resistance testing by DNA sequencing of partial pol gene. Blood samples were collected in the period 2006–2020, in the scope of routine laboratory monitoring. The eligibility criteria for participant selection included: HIV-1subtype B infected; seropositive adults (men and non-pregnant women) aged 18 years or above; with a place of residence in St. Petersburg; and viral load higher than 1000 copies per ml of blood.
Epidemiological, clinical, and behavioral data were collected using a standardized questionnaire. Transmission risk was categorized as heterosexual contact, MSM, PWID, transfusion, vertical transmission, or unknown.

2.2. RNA Isolation, RT PCR and DNA Sequencing

Plasma samples were genotyped using the ViroSeq HIV-1 Genotyping System (Abbott, Chicago, IL, USA) or AmpliSens® HIV-Resist-Seq (InterLabService Ltd., Moscow, Russia) in accordance with the manufacturer’s instructions.

2.3. Phylogenetic Analysis

Before starting this phylogenetic investigation, dataset was carefully checked for the presence of multiple sequences from the same patent, longitudinally followed during the study period. In the final dataset only single sequence per patient was included. Sequence alignment was performed by the Clustal W algorithm implemented in the MEGA 10.2 software, followed by maximum likelihood (ML) phylogenetic reconstruction. (https://www.megasoftware.net/, accessed on 29 November 2022). Prior to further phylogenetic analysis, sequences were codon-stripped for major drug resistance-associated positions, to avoid differential selection pressure bias between treated and naïve patients. [21]. The best fitting nucleotide-substitution model for the representative dataset, as determined using model selection employed in the jModeltest software, was the general time reversible model of nucleotide-substitution with a proportion of invariant sites (i) and gamma distribution of rates (G) (GTR + I + G), selected according to the Akaike Information Criterion (AIC), by [22,23]. Bootstrap support for the tree nodes of the reconstructed phylogenetic trees were calculated with 1000 replicates. Genetic distance calculations and phylogenetic tree topologies was also confirmed under the ML algorithm in the Paup software. Phylogenetic trees were visualized using the FigTree program v 1.4.2 [24]. Bayesian phylogenetic analyses, estimating the posterior distributions of phylogenetic tree topologies, were also performed, using Markov chain Monte Carlo (MCMC) methods in the MrBayes software v 3.2.7 [25].
Transmission clusters investigation was performed through an analysis of a series of criteria sets as described previously, aimed to maintain specificity in identifying true transmission clusters, yet to avoid an underestimation of the incidence of transmission chains, in the view of an extended 15 years timeframe of sampling in the present study [26,27,28]. According to the first set of criteria phylogenetic clades consisting of three or more sequences, with the genetic distance of 4.5% or less, minimal bootstrap support of 90% and the statistical support of a posterior probability higher than 0.9 in the Bayesian maximum clade credibility tree (MCCT), were considered as transmission clusters. A second criteria set was based on bootstrap support over 75% and genetic distance of less than 4.5%. These phylogenetic clades were further confirmed by means of BLAST analysis (http://blast.ncbi.nlm.nih.gov accessed on 10 May 2022), with the purpose to identify highly similar sequences, based on score significance, not included in the initial dataset. In particular, for every sequence belonging to a cluster identified by the second criteria set, five most similar sequences identified using the BLAST tool were included in tree reconstruction as previously described [26,28]. In the posterior analyses, only those clades with no intermixing of foreign sequences downloaded using BLAST, were considered as transmission clusters. In total, 220 BLAST identified sequences were retrieved from the NCBI database and included in the tree reconstruction through phylogenetic analyses as described.

2.4. Estimation of the Evolutionary History and the Effective Number of Infections

Bayesian Markov Chain Monte Carlo (MCMC) inference was applied to estimate the effective population size and tMRCA of HIV subtypes B, using the BEAST v1.10.4 software [29,30]. Bayesian analysis was performed using the following parameters: uncorrelated lognormal relaxed molecular clock; GTR nucleotide substitution model; estimated base frequencies, and gamma distribution model + invariant sites for heterogeneity among nucleotide sites. The analysis was performed using Bayesian skyline plot as tree prior. The MCMC chain length was set at 2 × 108, with a burn-in of 10%, which gave an effective sample size (ESS) of >200. The output was analyzed using TRACER v1.7.1, with 5% uncertainty in parameter estimates reflected as the 95% highest posterior density (HPD) [31]. The maximum clade credibility (MCC) tree was obtained with the TreeAnnotator software v1.8.3 and visualized in FigTree v1.4.2. [24].

2.5. Birth-Death Skyline (BDSKY] Analysis

For this investigation, aimed to explore effective reproductive number (Re) over time, the BEAST 2 software package v 2.6.5 was used [32,33,34]. The effective reproduction number (Re), denotes the average number of secondary infections caused by an infected person at a given time during the epidemic, where Re above 1 means that the number of cases are increasing and Re below 1 means that the epidemic will die out. As included sequences are sampled at different time points, Birth Death Skyline Serial model and lognormal prior for Re was chosen. Other priors were set based on prior knowledge of the trends of HIV epidemic and analyzed dataset, as follows: the rate of becoming uninfectious was set as a log-normal prior with M = 1.8 and S = 0.2 (95% CI 1.7–10.2 corresponding to the maximal natural infective period of 10 years in HIV patients), as reported previously [35], the sampling probability was estimated using a prior beta (1.0, 9999) corresponding to a minority of cases sampled, while the number of dimension for the reproductive number was set to 10 for the whole dataset and to five for the four analyzed transmission clusters. The number of dimensions was set according the space between the origin of the epidemic and the present time of the investigated clade, reflecting the number of changes of Re during the investigated period [32]. In order to plot the output of the BDSKY analysis, the R software with the bdskytools software package, available on GitHub, was used. (https://github.com/cran, accessed on 30 May 2022).

2.6. Latent Class Analysis

Latent class analysis (LCA) was performed in order to reveal hidden subgroups or ‘discrete classes’ of participants’ epidemiological, demographic, and clinical profiles, especially to identify how these characteristics of the patient’s affect transmission clustering patterns of the viral sequences. The R software with polytomous latent class analysis software package “poLCA” distributed through the Comprehensive R Archive Network, CRAN, was used [36]. Herein, we examined participants risk profiles regarding eight types of latent class indicators (categorical latent variable). Latent class indicators included: gender, age, transmission risk, date of first reactive HIV test (time-period of diagnosis), CDC stage at the time of diagnosis, CD4 count, city/country where HIV infection was most probably acquired, and the inclusion of sequence within clusters/network. The analysis was begun with a 1-class model and increased the number of classes in each subsequent model seeking to minimize both the BIC and the AIC value before these values increased with the addition of another class. With nreps = 10, the latent class with the lowest value of both BIC and AIC was chosen.

2.7. Statisticall Analysis

To illustrate the overall epidemiological context of HIV epidemics in St. Petersburg, surveillance data on HIV incidence for the period 1999–2020 for Russia and St. Petersburg were obtained from the publically available national data (Rospotrebnadzor Reports (HIV Infection in the Russian Federation). Available online: http://www.hivrussia.info/dannyepo-vich-infektsii-v-rossii/ accessed on 24 May 2022) and reports by St. Petersburg City AIDS Center, St. Petersburg, Russia (personal communication).
Linear regression and descriptive statistics were conducted using standard statistical methods.

3. Results

3.1. Study Population and HIV Incidence Trends

In total, 247 patients were included in the study. After exclusion of mulitple sequences per patient, the final sequence dataset contained partial pol gene sequences from 198 patients, single sequence per patient. GenBank accession numbers of the included sequences are given in the Supplementary File S1. The majority of patients included in the study were on highly active antiretroviral therapy (HAART) 75.7% (150/198), at the time of study. The overall population was 93% (185/198) men, with a great majority, 69.7%, reporting sex with men as risk for infection. Regarding the overall study population, MSM was the most frequently reported risk for infection, in 65.1% of all study patients (129/198), followed by heterosexual contact in 25.7% (51/198). Intravenous drug use (IVDU) was reported by a very small proportion of the study population, 9% (18/198), and vertical transmission in one patient.
A linear regression analysis of HIV incidence data for Russia and St. Petersburg revealed differing epidemiological trends (Figure 1). The overall downward trend in HIV incidence was observed in St. Petersburg, whereas an upward trend in HIV incidence is seen on the national level. Furthermore of note, is the fact that up to 2012 the epidemiological situation in St. Petersburg was worse than in Russia, illustrated by higher levels of HIV incidence in the city.

3.2. Phylogenetic Investigation of Transmission Clusters

A phylogenetic analysis of the dataset encompassing 198 HIV subtype B sequences from Saint Petersburg, sampled from 2006 to 2020, revealed the presence of 17 transmission clusters and 16 transmission pairs (Figure 2).
In total, 54% (107/198) of analyzed subtype B sequences were found grouped in 17 clusters that accomplished predefined sets of criteria. Clusters ranged from three to 14 sequences, with collection dates within the span of 1–12 years, and median pairwise genetic distance of 2.3% (range 0.1–4.5%). All identified transmission clusters fulfilled a second, less stringent criteria set, whereas 10 clusters, containing 58.8% (63/107) of the total number of clustering sequences, were defined by the first, most stringent criteria set. Four transmission clusters were found with over 10 sequences. Furthermore, 15.6% (32/198) sequences were identified within transmission pairs. Importantly, confirmatory cluster analysis with sequences obtained through BLAST search showed very few sequences identified by BLAST analysis to intermix within the identified transmission clusters. Of note, except for a single sequence reported from Cyprus, virtually no foreign sequence was BLAST identified, all of them had been sampled in Russia.

3.3. Timing the Origin and the Effective Population Size

Using Bayesian MCMC inference, tMRCA of HIV-1 subtype B was estimated to be around 1986 (95% HPD 1984–1991), Furthermore, temporal origin was also estimated for four local transmission clusters with the number of sequences above 10. Considering estimated temporal origin for the clusters 1, 2, 3 and 4, it was found to be in the early two-thousands. The median tMRCA for clusters 1, 2, 3 and 4 were estimated at 2002 (95% HPD 1997–2006), 2003 (95% HPD 1997–2006), 2001 (95% HPD 1995–2007) and 2005 (95% HPD 2000–2007), respectively.
A BSP coalescent tree prior enables the estimation of the effective population growth from the analyzed data of infected subjects as it progresses from the origin of the epidemic through time. The BSP analysis of the whole subtype B dataset identified exponential growth in viral effective population size (correlating with effective number of infections and/or transmission opportunities) from the mid-nineties to the mid-noughties (2005–2006) followed by constant growth phase onwards (Figure 3).
Effective population growth in four analyzed transmission clusters showed similar patterns as for the whole dataset with exponential growth phase starting from the beginning of the 2000s followed by an asymptotic phase approaching the present time.

3.4. Estimation of the Effective Reproductive Number

An estimation of the Re in BEAST 2.4.6 software package was firstly performed on the whole dataset, including 198 subtype B sequences (Figure 3). Results revealed the increase in Re value above 1 in the late nineties with the highest Re reaching the value of 1.5 from the beginning of 2000 to 2005, followed by a sharp decrease, however still reaming the value of around 1 till the end of the analyzed period. Secondly, Re was estimated for four largest transmission clusters revealing the value above 1 in all investigated clusters since their origination, followed gradually by decline. Of note, all four transmission clusters had a value above 1 approaching the present time (Figure 4).

3.5. Latent Class Analysis

LCA, performed aimed to correlate results of transmission clustering patterns obtained by phylogenetic analysis with available public health and clinical determinants, revealed the presence of four latent classes, based on AIC and BIC values (Figure 5).
Large transmission clusters, containing > 10 sequences, were dominantly present in the first latent class, which encompassed male subjects with sexually acquired HIV infection, predominantly MSM, late presenters aged between 30–50, with low CD4 cell count (below 100), and HIV diagnosed in the last 10 years of the study period. Notably, around one third of patients in this latent class reported HIV infection most probably acquired outside Saint Petersburg. The second latent class encompassed predominantly younger (under 30) MSM subjects, with HIV diagnosis in the preceding five years and presenting in the early stage of infection. The second latent class was the one encompassing smaller transmission clusters with the number of sequences up to 10. The third latent class encompassed an equal proportion of men and women, PWIDs and reporting heterosexual transmission risk for HIV infection, with the time of diagnosis dispersed across the whole study period (2006 to 2020) and mostly presenting at a late stage of HIV disease. The fourth latent class encompassed male subjects aged over five, reporting MSM risk and presenting in the late stage of HIV infection, with diagnosis of HIV infection established during the first 10 years of the study period.

4. Discussion

Here, we present the study that integrates an analysis of viral phylogenetics and analysis of socio-demographic metadata regarding HIV subtype B spread in St. Petersburg, Russia. Phylogenetic methods can provide unique insights in the transmission networks and the spread of the virus and have been used worldwide to map local HIV epidemics in correlation with transmission pathway, drug resistance, risk behavior, and cluster size [26,28,37]. Many studies were based on partial pol gene phylogeny, which has shown to be adequate to infer transmission events and to characterize epidemiological patterns of public health relevance [38]. When coupled to metadata analysis, they may give a quantitative description of transmission networks, by identifying socio-demographic correlates of clustering [28,35,39]. Furthermore, detailed analysis of HIV burden distribution may prompt targeted epidemiological analyses and preventive measures to reduce epidemic spread [40,41,42]. Insight into overall HIV incidence trend in St. Petersburg, and at the national level, reveals an opposing trend—a downward trend in HIV incidence in St. Petersburg vs. growing HIV epidemic in Russia. Historically, the first cases of HIV infection in Russia were registered in St. Petersburg long before the start of a large-scale epidemic. This early epidemic phase was linked to African students enrolled at the universities of the city at the time, that were since mostly lost to follow-up (either left the country or died), and the underlying HIV subtypes were not determined. Over the next years, certain cases were identified related to the introduction from Europe and the United States, all associated to MSM transmission and caused by HIV-1 subtype B [20,43]. Later, a large-scale outbreak of sub-subtype A6 occurred, mainly among PWIDs, overshadowing the continuing “hidden”, predominantly subtype B spread among MSM [43]. In the last decade, a changing pattern of HIV subtype distribution has been observed, involving all key populations [9,17]. This phenomenon required a deeper study, since the knowledge of the molecular structure of the epidemic as a whole and in separate groups gives grounds for the development of aimed preventive measures.
Previous studies have depicted transmission clusters in HIV spread in different regions of Russia [9,10]. To the best of our knowledge, the present study is the first one to address transmission clusters of HIV subtype B in St. Petersburg. In addition to 15.6% of analyzed subtype B sequences forming transmission pairs, our results imply that 54% of analyzed subtype B sequences, collected within the time span of 15 years, contributed to 17 clusters (of three sequences at minimum), similar as reported in comparable studies elsewhere [26,28,44,45,46,47]. Notably, confirmatory cluster analysis implies a completely autonomous Russian epidemic, without any significant intermixing of foreign sequences obtained through a BLAST search into the designated transmission clusters. The estimation of the effective reproductive number (Re) for the whole subtype B dataset, as well as for the four largest transmission clusters (containing 10 or more sequences), revealed a steady value above 1 in all investigated clades, implying a permanent increase in population size. In spite of the more or less pronounced gradual decline after the mid-noughties (2005–2006) all the four transmission clusters retained Re value above 1 until the end of the study period.
An estimation of the time of origin of HIV-1 subtype B in St. Petersburg, within the presented study, was to be around 1986 (95% HPD 1984–1991), in line with historical data and similar to the estimate regarding subtype B epidemic in Moscow [9]. This estimate coincides with major events of sociopolitical changes and turbulence in the former Soviet Union and Russia, when travel to Europe and USA, as well as contacts with the population of these countries, increased significantly. At the same time, other social events took place during “perestroika”, in particular an increase in illicit drug use, which, apparently, contributed to the spread of subtype B virus from Ukraine to Russia, and also to the subsequent formation of the CRF03_AB recombinant. Regarding the four large (over 10 sequences) subtype B St. Petersburg clusters, their origin has been estimated to be more recent, ranging from 2001 to 2005. This can be attributed to the changing pattern of HIV epidemic—notably, the explosive spread of HIV in Russia has been linked to illicit drug use, an epidemic predominantly caused by HIV sub-subtype A6. However, recent data imply an increasing dominance of sexual transmission, especially among newly diagnosed patients [1,2]. Notably, St. Petersburg is the first urban metropolis in the Russian Federation where a steady decline in new HIV infections has been recorded [48]. Nationwide, a concentrated HIV epidemic among people who inject drugs tends to dissolve to the generalized one, with prevalence in some regions reaching alarming heights [2]. However, through the City AIDS Center, St. Petersburg is among the few cities in the Russian Federation that has provided harm reduction programs, albeit on a limited scale [48]. With over 900 people living with HIV (PLWH) per 100,000 population, St. Petersburg is still among the Russian regions with the highest HIV prevalence, nevertheless, the incidence rate in the city is decreasing. Changing practices have been described among PWIDs in St. Petersburg, diminishing the risk of ongoing HIV transmission in this key population [49]. A recent study has shown the attribution of the transmission pathway with different HIV subtypes in St. Petersburg, with subtype B epidemic being linked to MSM transmission [43]. In the presented study, the overall population was 93% men, with over two thirds reporting sex with men as the risk for infection, while intravenous drug use was reported by a very small proportion of the population, 9%. At present, in Russian society, the disclosure of MSM infection risk is linked to substantial social stigma, even personal risk, making the studied dataset rather unique. Hence, complex metadata analyses may provide valuable insight into social networks implicated in HIV transmission and give grounds to planed preventive interventions. We performed latent class analysis, as an approach to reveal non-evident subgroups or ‘discrete classes’ of participants’ epidemiological, demographic, and clinical profiles within the study population, and correlate those with clustering pattern. Through this approach, four latent classes have been identified, with clearly distinguished social and epidemiological features. Two of the four identified latent classes encompassed mostly men with MSM transmission risk and were linked to transmission clusters. Larger clusters were associated to middle aged subjects with epidemiological links outside the city and more advanced HIV disease, whereas smaller clusters were associated to younger men with early stage HIV disease. In conjunction with the obtained results of phylodynamic and temporal analysis of transmission clusters, this points to the rising local spread of HIV subtype B among younger MSMs, who tend to present to care at earlier disease stage, along with maintaining continual spread within an older age group and more advanced HIV disease.
Of note, in view of the registered number of HIV cases in St. Petersburg of over 60,000, the sample size of the studied dataset may seem limited, and thus present as possible study limitation. However, in this regard, the documented low prevalence of subtype B in Russian epidemic should be considered, notably about 4% in the northwestern federal distric, where St. Petesburg is the main urban center [50]. Further, in Russia, the standard of care for HIV patients does not include genotyping prior to therapy initiation. Even in the case of therapy failure, access to genotyping is available for around 10% of patients in need of this analysis. Thus, the present study encompasses the total of all HIV subtype B sequence data, which have been obtained since genotypic resistance testing has been implemented in clinical care in St. Petersburg.
In conclusion, the results of our study imply a complex pattern of the epidemic spread of HIV subtype B in St. Petersburg, Russia, still in the exponential growth phase and in connection to the MSM transmission, highlighting the key population of young MSM as drivers of HIV subtype B epidemic, and providing a useful insight needed for the design of public health priorities and interventions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14122748/s1, Siljic et al_S1 supplementary file

Author Contributions

M.B. and M.S. (Maja Stanojevic) designed the study; M.S. (Marina Siljic), V.C., M.S. (Maja Stanojevic) and L.J. conducted the analyses; M.S. (Marina Siljic), V.C. and M.S. (Maja Stanojevic) drafted the manuscript; A.A., A.L., E.O., A.K., A.E. and N.M. contributed to database validation, curation and analysis; T.V. and M.B. critically edited the manuscript. 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 23 October 2022) and by the Ministry of Education, Science and Technological Development of the Republic of Serbia, grant No. 200110.

Institutional Review Board Statement

This study was carried out in accordance with the Declaration of Helsinki. The study was approved by “Health committee of the government of Saint Petersburg. Saint Petersburg state budget health institution Center for prevention and control of AIDS and infectious diseases”, Protocol No199.

Informed Consent Statement

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

Data Availability Statement

NCBI accession numbers of the analyzed sequences are available as Supplementary material.

Conflicts of Interest

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

References

  1. European Centre for Disease Prevention and Control/WHO Regional Office for Europe. HIV/AIDS Surveillance in Europe 2021–2020 Data; ECDC: Stockholm, Sweden, 2021.
  2. Federal Russian AIDS Center. Report 2021. Available online: http://www.hivrussia.info/wp-content/uploads/2022/05/Byulleten-46-VICH-infektsiya-za-2020-g.-.pdf (accessed on 30 May 2022). (In Russian).
  3. Beyrer, C.; Wirtz, A.L.; O’Hara, G.; Léon, N.; Kazatchkine, M. The expanding epidemic of HIV-1 in the Russian Federation. PLoS Med. 2017, 14, e1002462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Bobkov, A.; Kazennova, E.; Selimova, L.; Bobkova, M.; Khanina, T.; Ladnaya, N.; Kravchenko, A.; Pokrovsky, V.; Cheingsong-Popov, R.; Weber, J. A sudden epidemic of HIV type 1 among injecting drug users in the former Soviet Union: Identification of subtype A.; subtype B.; and novel gagA/envB recombinants. AIDS Res. Hum. Retrovir. 1998, 14, 669–676. [Google Scholar] [CrossRef] [PubMed]
  5. Bobkov, A.; Kazennova, E.; Khanina, T.; Bobkova, M.; Selimova, L.; Kravchenko, A.; Pokrovsky, V.; Weber, J. An HIV type 1 subtype A strain of low genetic diversity continues to spread among injecting drug users in Russia: Study of the new local outbreaks in Moscow and Irkutsk. AIDS Res. Hum. Retrovir. 2001, 17, 257–261. [Google Scholar] [CrossRef] [PubMed]
  6. Gashnikova, N.M.; Bogachev, V.V.; Baryshev, P.B.; Totmenin, A.V.; Gashnikova, M.P.; Kazachinskaya, A.G.; Ismailova, T.N.; Stepanova, S.A.; Chernov, A.S.; Mikheev, V.N. A rapid expansion of HIV-1 CRF63_02A1 among newly diagnosed HIV-infected individuals in the Tomsk Region.; Russia. AIDS Res. Hum. Retrovir. 2015, 31, 456–460. [Google Scholar] [CrossRef]
  7. Gashnikova, N.M.; Zyryanova, D.P.; Astakhova, E.M.; Ivlev, V.V.; Gashnikova, M.P.; Moskaleva, N.V.; Aikin, S.S.; Bulatova, T.N.; Pustylnikov, S.V.; Bocharov, E.F.; et al. 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. 2017, 162, 379–390. [Google Scholar] [CrossRef] [PubMed]
  8. Karamov, E.; Epremyan, K.; Siniavin, A.; Zhernov, Y.; Cuevas, M.T.; Delgado, E.; Sánchez-Martínez, M.; Carrera, C.; Kornilaeva, G.; Turgiev, A.; et al. HIV-1 Genetic Diversity in Recently Diagnosed Infections in Moscow: Predominance of AFSU.; Frequent Branching in Clusters.; and Circulation of the Iberian Subtype G Variant. AIDS Res. Hum. Retrovir. 2018, 34, 629–634. [Google Scholar] [CrossRef]
  9. 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]
  10. Neshumaev, D.; Lebedev, A.; Malysheva, M.; Boyko, A.; Skudarnov, S.; Ozhmegova, E.; Antonova, A.; Kazennova, E.; Bobkova, M. Molecular Surveillance of HIV-1 Infection in Krasnoyarsk Region, Russia: Epidemiology, Phylodynamics and Phylogeography. Curr. HIV Res. 2019, 17, 114–125. [Google Scholar] [CrossRef]
  11. Maksimenko, L.V.; Totmenin, A.V.; Gashnikova, M.P.; Astakhova, E.M.; Skudarnov, S.E.; Ostapova, T.S.; Yaschenko, S.V.; Meshkov, I.O.; Bocharov, E.F.; Maksyutov, R.A.; et al. Genetic Diversity of HIV-1 in Krasnoyarsk Krai: Area with High Levels of HIV-1 Recombination in Russia. BioMed Res. Int. 2020, 2020, 9057541. [Google Scholar] [CrossRef]
  12. 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]
  13. 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]
  14. Nabatov, A.A.; Kravchenko, O.N.; Lyulchuk, M.G.; Shcherbinskaya, A.M.; Lukashov, V.V. Simultaneous introduction of HIV type 1 subtype A and B viruses into injecting drug users in southern Ukraine at the beginning of the epidemic in the former Soviet Union. AIDS Res. Hum. Retrovir. 2002, 18, 891–895. [Google Scholar] [CrossRef] [PubMed]
  15. Lukashov, V.V.; Cornelissen, M.T.; Goudsmit, J.; Papuashvilli, M.N.; Rytik, P.G.; Khaitov, R.M.; Karamov, E.V.; de Wolf, F. Simultaneous introduction of distinct HIV-1 subtypes into different risk groups in Russia, Byelorussia and Lithuania. AIDS 1995, 9, 435–439. [Google Scholar] [CrossRef]
  16. Kazennova, E.; Laga, V.; Lapovok, I.; Glushchenko, N.; Neshumaev, D.; Vasilyev, A.; Bobkova, M. HIV-1 genetic variants in the Russian Far East. AIDS Res. Hum. Retrovir. 2014, 30, 742–752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Kazennova, E.; Laga, V.; Gromov, K.; Lebedeva, N.; Zhukova, E.; Pronin, A.; Grezina, L.; Dement’eva, N.; Shemshura, A.; Bobkova, M. Genetic Variants of HIV Type 1 in Men Who Have Sex with Men in Russia. AIDS Res. Hum. Retrovir. 2017, 33, 1061–1064. [Google Scholar] [CrossRef]
  18. Kozlov, A.P.; Volkova, G.V.; Malykh, A.G.; Stepanova, G.S.; Glebov, A.V. Epidemiology of HIV infection in St. Petersburg.; Russia. J. Acquir. Immune. Defic. Syndr. 1993, 6, 208–212. [Google Scholar]
  19. Lukashov, V.V.; Huismans, R.; Rakhmanova, A.G.; Lisitsina, Z.N.; Akhtyrskaya, N.A.; Vlasov, N.N.; Melnick, O.B.; Goudsmit, J. Circulation of subtype A and gagA/envB recombinant HIV type 1 strains among injecting drug users in St. Petersburg, Russia, correlates with geographical origin of infections. AIDS Res. Hum. Retrovir. 1999, 15, 1577–1583. [Google Scholar] [CrossRef]
  20. Thomson, M.M.; Vinogradova, A.; Delgado, E.; Rakhmanova, A.; Yakovlev, A.; Cuevas, M.T.; Muñoz, M.; Pinilla, M.; Vega, Y.; Pérez-Alvarez, L.; et al. Molecular epidemiology of HIV-1 in St. Petersburg, Russia: Predominance of subtype A.; former Soviet Union variant, and identification of intrasubtype subclusters. J. Acquir. Immune Defic. Syndr. 2009, 51, 332–339. [Google Scholar] [CrossRef]
  21. Wensing, A.M.; Calvez, V.; Ceccherini-Silberstein, F.; Charpentier, C.; Günthard, H.F.; Paredes, R.; Shafer, R.W.; Richman, D.D. Update of the drug resistance mutations in HIV-1. Top Antivir. Med. 2019, 27, 111–121. [Google Scholar]
  22. Posada, D. JModelTest: Phylogenetic model averaging. Mol. Biol. Evol. 2008, 25, 1253–1256. [Google Scholar] [CrossRef]
  23. Shapiro, B.; Rambaut, A.; Drummond, A.J. Choosing appropriate substitution models for the phylogenetic analysis of protein-coding sequences. Mol. Biol. Evol. 2006, 23, 7–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Rambaut, A. FigTree v1.3.1: Tree Figure Drawing Tool. 2009. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 30 May 2022).
  25. Ronquist, F.; Huelsenbeck, J.P. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19, 1572–1574. [Google Scholar] [CrossRef] [PubMed]
  26. Siljic, M.; Salemovic, D.; Jevtovic, D.; Pesic-Pavlovic, I.; Zerjav, S.; Nikolic, V.; Ranin, J.; Stanojevic, M. Molecular typing of the local HIV-1 epidemic in Serbia. Infect. Genet. Evol. 2013, 19, 378–385. [Google Scholar] [CrossRef] [PubMed]
  27. Siljic, M.; Salemovic, D.; Cirkovic, V.; Pesic-Pavlovic, I.; Ranin, J.; Todorovic, M.; Nikolic, S.; Jevtovic, D.; Stanojevic, M. Forensic application of phylogenetic analyses—Exploration of suspected HIV-1 transmission case. Forensic Sci. Int. Genet. 2017, 27, 100–105. [Google Scholar] [CrossRef]
  28. Jovanovic, L.; Siljic, M.; Cirkovic, V.; Salemovic, D.; Pesic-Pavlovic, I.; Todorovic, M.; Ranin, J.; Jevtović, D.; Stanojević, M. Exploring evolutionary and transmission dynamics of HIV epidemic in Serbia: Bridging socio-demographic with phylogenetic approach. Front. Microbiol. 2019, 10, 287. [Google Scholar] [CrossRef] [Green Version]
  29. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 2012, 29, 1969–1973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Suchard, M.A.; Lemey, P.; Baele, G.; Ayres, D.L.; Drummond, A.J.; Rambaut, A. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 2018, 4, vey016. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Rambaut, A.; Drummond, A.J.; Xie, D.; Baele, G.; Suchard, M.A. Posterior summarisation in Bayesian phylogenetics using Tracer 1.7. Syst. Biol. 2018, 67, 901–904. [Google Scholar] [CrossRef] [Green Version]
  32. Stadler, T.; Kühnert, D.; Bonhoeffer, S.; Drummond, A.J. Birth-death skyline plot reveals temporal changes of epidemic spread in HIV and hepatitis C virus (HCV). Proc. Natl. Acad. Sci. USA 2013, 110, 228–233. [Google Scholar] [CrossRef] [Green Version]
  33. Bouckaert, R.; Heled, J.; Kühnert, D.; Vaughan, T.; Wu, C.H.; Xie, D.; Suchard, M.A.; Rambaut, A.; Drummond, A.J. Beast 2: A software platform for bayesian evolutionary analysis. PLoS Comput. Biol. 2014, 10, e1003537. [Google Scholar] [CrossRef] [Green Version]
  34. Bouckaert, R.; Vaughan, T.G.; Barido-Sottani, J.; Duchêne, S.; Fourment, M.; Gavryushkina, A.; Heled, J.; Jones, G.; Kühnert, D.; De Maio, N.; et al. Beast 2.5: An advanced software platform for bayesian evolutionary analysis. PLoS Comput. Biol. 2019, 15, 1006650. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Paraskevis, D.; Kostaki, E.; Gargalianos, P.; Xylomenos, G.; Lazanas, M.; Chini, M.; Skoutelis, A.; Papastamopoulos, V.; Paraskeva, D.; Antoniadou, A. Transmission Dynamics of HIV-1 Drug Resistance among Treatment-Naïve Individuals in Greece: The Added Value of Molecular Epidemiology to Public Health. Genes 2017, 8, 322. [Google Scholar] [CrossRef] [PubMed]
  36. Linzer, D.A.; Lewis, J.B. PoLCA: An r package for polytomous variable latent class analysis. J. Stat. Softw. 2011, 42, 1–29. [Google Scholar] [CrossRef] [Green Version]
  37. Van de Vijver, D.A.M.C.; Boucher, C.A.B. Insights on transmission of HIV from phylogenetic analysis to locally optimize HIV prevention strategies. Curr. Opin. HIV AIDS 2018, 13, 95–101. [Google Scholar] [CrossRef]
  38. Hué, S.; Clewley, J.P.; Cane, P.A.; Pillay, D. HIV-1 pol gene variation is sufficient for reconstruction of transmissions in the era of antiretroviral therapy. AIDS 2004, 18, 719–728. [Google Scholar] [CrossRef]
  39. Kostaki, E.; Magiorkinis, G.; Psichogiou, M.; Flampouris, A.; Iliopoulos, P.; Papachristou, E.; Daikos, G.L.; Bonovas, S.; Otelea, D.; Friedman, S.R.; et al. Detailed Molecular Surveillance of the HIV-1 Outbreak Among People who Inject Drugs (PWID) in Athens During a Period of Four Years. Curr. HIV. Res. 2017, 15, 396–404. [Google Scholar] [CrossRef]
  40. Shiino, T.; Hachiya, A.; Hattori, J.; Sugiura, W.; Yoshimura, K. Nation-Wide Viral Sequence Analysis of HIV-1 Subtype B Epidemic in 2003–2012 Revealed a Contribution of Men Who Have Sex with Men to the Transmission Cluster Formation and Growth in Japan. Front. Reprod. Health 2020, 2, 531212. [Google Scholar] [CrossRef]
  41. 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]
  42. Li, X.; Li, Y.; Liu, H.; Trovão, N.S.; Foley, B.T. The emergence and transmission dynamics of HIV-1 CRF07_BC in Mainland China. Virus Evol. 2022, 8, veac014. [Google Scholar] [CrossRef]
  43. Leigh Brown, A.J.; Lycett, S.J.; Weinert, L.; Hughes, G.J.; Fearnhill, E.; Dunn, D.T.; UK HIV Drug Resistance Collaboration. Transmission network parameters estimated from HIV sequences for a nationwide epidemic. J. Infect. Dis. 2011, 204, 1463–1469. [Google Scholar] [CrossRef] [Green Version]
  44. Ragonnet-Cronin, M.; Lycett, S.J.; Hodcroft, E.B.; Hué, S.; Fearnhill, E.; Brown, A.E.; Delpech, V.; Dunn, D.; Leigh Brown, A.J. United Kingdom HIV Drug Resistance Database. Transmission of Non-B HIV Subtypes in the United Kingdom Is Increasingly Driven by Large Non-Heterosexual Transmission Clusters. J. Infect. Dis. 2016, 213, 1410–1418. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Chaillon, A.; Essat, A.; Frange, P.; Smith, D.M.; Delaugerre, C.; Barin, F.; Ghosn, J.; Pialoux, G.; Robineau, O.; Rouzioux, C.; et al. Spatiotemporal dynamics of HIV-1 transmission in France (1999–2014) and impact of targeted prevention strategies. Retrovirology 2017, 14, 15. [Google Scholar] [CrossRef] [PubMed]
  46. Parczewski, M.; Leszczyszyn-Pynka, M.; Witak-Jędra, M.; Szetela, B.; Gąsiorowski, J.; Knysz, B.; Bociąga-Jasik, M.; Skwara, P.; Grzeszczuk, A.; Jankowska, M.; et al. Expanding HIV-1 subtype B transmission networks among men who have sex with men in Poland. PLoS ONE 2017, 12, e0172473. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Preventing and Treating HIV in Saint Petersburg. Available online: https://www.unaids.org/en/resources/presscentre/featurestories/2018/february/20180220_saint_petersburg (accessed on 30 May 2022).
  48. Meylakhs, P.; Friedman, S.R.; Meylakhs, A.; Mateu-Gelabert, P.; Ompad, D.C.; Alieva, A.; Dmitrieva, A. A New Generation of Drug Users in St. Petersburg, Russia: HIV, HCV, and Overdose Risks in a Mixed-Methods Pilot Study of Young Hard Drug Users. AIDS Behav. 2019, 23, 3350–3365. [Google Scholar] [CrossRef] [PubMed]
  49. Dukhovlinova, E.; Masharsky, A.; Toussova, O.; Verevochkin, S.; Solovyeva, T.; Meringof, M.; Paintsil, E.; White, E.; Barbour, R.; Heimer, R.; et al. Two Independent HIV Epidemics in Saint Petersburg, Russia Revealed by Molecular Epidemiology. AIDS Res. Hum. Retrovir. 2015, 31, 608–614. [Google Scholar] [CrossRef] [Green Version]
  50. Shchemelev, A.N.; Ostankova, Y.V.; Zueva, E.B.; Semenov, A.V.; Totolian, A.A. Detection of Patient HIV-1 Drug Resistance Mutations in Russia’s Northwestern Federal District in Patients with Treatment Failure. Diagnostics 2022, 12, 1821. [Google Scholar] [CrossRef]
Figure 1. Linear regression trend of HIV incidence on 100,000 inhabitans in Russia (left) and St. Petesburg (right).
Figure 1. Linear regression trend of HIV incidence on 100,000 inhabitans in Russia (left) and St. Petesburg (right).
Viruses 14 02748 g001
Figure 2. ML phylogenetic tree based on 198 HIV subtype B pol sequences sampled from 2006 to 2020 in St. Petersburg; tree was constructed in MEGA software under the GTR + G + I nucleotide substitution model. Transmission clusters are labeled in red. Transmission pairs are mark in green.
Figure 2. ML phylogenetic tree based on 198 HIV subtype B pol sequences sampled from 2006 to 2020 in St. Petersburg; tree was constructed in MEGA software under the GTR + G + I nucleotide substitution model. Transmission clusters are labeled in red. Transmission pairs are mark in green.
Viruses 14 02748 g002
Figure 3. Bayesian skyline plot, representing estimation of the effective population growth of HIV subtype B in St. Petersburg (based on the ‘relaxed clock’ coalescent framework analysis on 198 sequences included in the study, in BEAST 1.10. software) (left); Effective reproductive number of the studied dataset of HIV subtype B in St. Petersburg, based on birth date skyline serial prior and 10 dimensions of Re in BEAST 2 software. Re plot was obtained in R software under the bdskytool software package (right).
Figure 3. Bayesian skyline plot, representing estimation of the effective population growth of HIV subtype B in St. Petersburg (based on the ‘relaxed clock’ coalescent framework analysis on 198 sequences included in the study, in BEAST 1.10. software) (left); Effective reproductive number of the studied dataset of HIV subtype B in St. Petersburg, based on birth date skyline serial prior and 10 dimensions of Re in BEAST 2 software. Re plot was obtained in R software under the bdskytool software package (right).
Viruses 14 02748 g003
Figure 4. Effective reproductive number of four large transmission clusters (>10 sequences) obtained by birth death skyline analysis performed in BEAST 2 software under the following parameters; relaxed lognormal clock, birth death skyline serial prior, reproductive number with five dimensions. Re plot was obtained in R software under the bdskytool software package.
Figure 4. Effective reproductive number of four large transmission clusters (>10 sequences) obtained by birth death skyline analysis performed in BEAST 2 software under the following parameters; relaxed lognormal clock, birth death skyline serial prior, reproductive number with five dimensions. Re plot was obtained in R software under the bdskytool software package.
Viruses 14 02748 g004
Figure 5. LCA result highlighted existence of four latent classes (indicated as I, II, III, IV), LCA was performed with poLCA software package in R, with study metadata and results of transmission clusters analysis as input. The distribution of participant characteristics within the latent class and predicted proportions in each cluster size by latent class are shown. Darker colors represent higher proportion as shown in the legend.
Figure 5. LCA result highlighted existence of four latent classes (indicated as I, II, III, IV), LCA was performed with poLCA software package in R, with study metadata and results of transmission clusters analysis as input. The distribution of participant characteristics within the latent class and predicted proportions in each cluster size by latent class are shown. Darker colors represent higher proportion as shown in the legend.
Viruses 14 02748 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Siljic, M.; Cirkovic, V.; Jovanovic, L.; Antonova, A.; Lebedev, A.; Ozhmegova, E.; Kuznetsova, A.; Vinogradova, T.; Ermakov, A.; Monakhov, N.; et al. Reconstructing the Temporal Origin and the Transmission Dynamics of the HIV Subtype B Epidemic in St. Petersburg, Russia. Viruses 2022, 14, 2748. https://doi.org/10.3390/v14122748

AMA Style

Siljic M, Cirkovic V, Jovanovic L, Antonova A, Lebedev A, Ozhmegova E, Kuznetsova A, Vinogradova T, Ermakov A, Monakhov N, et al. Reconstructing the Temporal Origin and the Transmission Dynamics of the HIV Subtype B Epidemic in St. Petersburg, Russia. Viruses. 2022; 14(12):2748. https://doi.org/10.3390/v14122748

Chicago/Turabian Style

Siljic, Marina, Valentina Cirkovic, Luka Jovanovic, Anastasiia Antonova, Aleksey Lebedev, Ekaterina Ozhmegova, Anna Kuznetsova, Tatiyana Vinogradova, Aleksei Ermakov, Nikita Monakhov, and et al. 2022. "Reconstructing the Temporal Origin and the Transmission Dynamics of the HIV Subtype B Epidemic in St. Petersburg, Russia" Viruses 14, no. 12: 2748. https://doi.org/10.3390/v14122748

APA Style

Siljic, M., Cirkovic, V., Jovanovic, L., Antonova, A., Lebedev, A., Ozhmegova, E., Kuznetsova, A., Vinogradova, T., Ermakov, A., Monakhov, N., Bobkova, M., & Stanojevic, M. (2022). Reconstructing the Temporal Origin and the Transmission Dynamics of the HIV Subtype B Epidemic in St. Petersburg, Russia. Viruses, 14(12), 2748. https://doi.org/10.3390/v14122748

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