Spatiotemporal Characteristics of the HIV-1 CRF02_AG/CRF63_02A1 Epidemic in Russia and Central Asia

Medical Interns' Knowledge, Attitude, and Practice Toward People Living with HIV: Multicenter Experience from Saudi Arabia

Background

Discrimination by some healthcare providers toward people living with HIV/AIDS has been documented. Differences in cultural backgrounds make it harder for future doctors, who need a lot of knowledge and a positive attitude to treat patients. In conservative countries like Saudi Arabia, not enough is known about how much medical interns know about HIV and how they feel about people living with HIV/AIDS.

Methods

From April to September 2021, this cross-sectional study use non-probability random sampling and utilized a self-administered questionnaire to collected the data from 346 medical interns who had graduated from five different medical schools.

Results

Most of the subjects correctly identified the main transmission routes, such as unprotected sex (94.57%), blood and body fluid exchange (94.19%), and sharing needles or syringes (91.47%). But they did not know what the most common co-infections were for HIV patients or how to protect themselves after exposure. This paper showed that medical interns have some stigmatizing behaviors toward patients living with HIV, as 31.1% and 22.9% agreed, respectively, that they would feel more sympathetic toward people who get AIDS from blood transfusions compared to IV drug users (IDU).

Conclusion

Medical interns also showed some positive attitudes, as more than half of the sample (56.2%) would not isolate beds for people living with HIV/AIDS. The study's conclusion is that HIV education and training programs should be added for medical interns, which might have a significant positive impact on their attitude.

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

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

Prevalence of Pretreatment HIV-1 Drug Resistance in Armenia in 2017-2018 and 2020-2021 following a WHO Survey

The increased antiretroviral therapy (ART) coverage of patients in the absence of routine genotyping tests and in the context of active labor migration highlight the importance of HIV-1 drug resistance (DR) surveillance in Armenia. We conducted a two-phase pretreatment DR (PDR) study in 2017-2018 (phase I; 120 patients) and 2020-2021 (phase II; 133 patients) according to the WHO-approved protocol. The analysis of HIV-1 genetic variants showed high degrees of viral diversity, with the predominance of A6. The prevalence of any PDR was 9.2% in phase I and 7.5% in phase II. PDR to protease inhibitors was found only in 0.8% in phase II. PDR to efavirenz and nevirapine was found among 5.0% and 6.7% of patients in phase I, and 6.0% and 6.8% of patients in phase II, respectively. The prevalence of PDR to nucleoside reverse-transcriptase inhibitors decreased from 5.0% in phase I to 0.8% in phase II. In addition, we identified risk factors associated with the emergence of DR-male, MSM, subtype B, and residence in or around the capital of Armenia-and showed the active spread of HIV-1 among MSM in transmission clusters, i.e., harboring DR, which requires the immediate attention of public health policymakers for the prevention of HIV-1 DR spread in the country.

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

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

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

Prevalence of HIV-1 drug resistance in Eastern European and Central Asian countries

BACKGROUND

Eastern Europe and Central Asia (EECA) is one of the regions where the HIV epidemic continues to grow at a concerning rate. Antiretroviral therapy (ART) coverage in EECA countries has significantly increased during the last decade, which can lead to an increase in the risk of emergence, transmission, and spread of HIV variants with drug resistance (DR) that cannot be controlled. Because HIV genotyping cannot be performed in these countries, data about HIV DR are limited or unavailable.

OBJECTIVES

To monitor circulating HIV-1 genetic variants, assess the prevalence of HIV DR among patients starting antiretroviral therapy, and reveal potential transmission clusters among patients in six EECA countries: Armenia, Azerbaijan, Belarus, Russia, Tajikistan, and Uzbekistan.

MATERIALS AND METHODS

We analyzed 1071 HIV-1 pol-gene fragment sequences (2253-3369 bp) from patients who were initiating or reinitiating first-line ART in six EECA counties, i.e., Armenia (n = 120), Azerbaijan (n = 96), Belarus (n = 158), Russia (n = 465), Tajikistan (n = 54), and Uzbekistan (n = 178), between 2017 and 2019. HIV Pretreatment DR (PDR) and drug resistance mutation (DRM) prevalence was estimated using the Stanford HIV Resistance Database. The PDR level was interpreted according to the WHO standard PDR survey protocols. HIV-1 subtypes were determined using the Stanford HIV Resistance Database and subsequently confirmed by phylogenetic analysis. Transmission clusters were determined using Cluster Picker.

RESULTS

Analyses of HIV subtypes showed that EECA, in general, has the same HIV genetic variants of sub-subtype A6, CRF63_02A1, and subtype B, with different frequencies and representation for each country. The prevalence of PDR to any drug class was 2.8% in Uzbekistan, 4.2% in Azerbaijan, 4.5% in Russia, 9.2% in Armenia, 13.9% in Belarus, and 16.7% in Tajikistan. PDR to protease inhibitors (PIs) was not detected in any country. PDR to nucleoside reverse-transcriptase inhibitors (NRTIs) was not detected among patients in Azerbaijan, and was relatively low in other countries, with the highest prevalence in Tajikistan (5.6%). The prevalence of PDR to nonnucleoside reverse-transcriptase inhibitors (NNRTIs) was the lowest in Uzbekistan (2.8%) and reached 11.1% and 11.4% in Tajikistan and Belarus, respectively. Genetic transmission network analyses identified 226/1071 (21.1%) linked individuals, forming 93 transmission clusters mainly containing two or three sequences. We found that the time since HIV diagnosis in clustered patients was significantly shorter than that in unclustered patients (1.26 years vs 2.74 years). Additionally, the K103N/S mutation was mainly observed in clustered sequences (6.2% vs 2.8%).

CONCLUSIONS

Our study demonstrated different PDR prevalence rates and DR dynamics in six EECA countries, with worrying levels of PDR in Tajikistan and Belarus, where prevalence exceeded the 10% threshold recommended by the WHO for immediate public health action. Because DR testing for clinical purposes is not common in EECA, it is currently extremely important to conduct surveillance of HIV DR in EECA due to the increased ART coverage in this region.

Genetic Diversity and Drug Resistance Mutations in Reverse Transcriptase and Protease Genes of HIV-1 Isolates from Southwestern Siberia

The analysis of a pol gene fragment encoding protease and part of reverse transcriptase was carried out for 55 sera collected in 2016 and 2018 from HIV-1-infected patients diagnosed in 2014-2018 living in the south of Western Siberia, Russia: Altai Territory (n = 11), Republic of Altai (n = 15), Kemerovo region (n = 18), and Novosibirsk region (n = 11). CRF63_02A was the dominant genetic form (>70%) in the Altai Territory and Kemerovo and Novosibirsk regions, with subsubtype A6 comprising <30% of samples. In the Altai Republic, subsubtype A6 was predominant (53%), with 33% of viruses belonging to CRF63_02A. Four CRF63_02A/A6 unique recombinant forms were identified in the Altai Territory, Kemerovo Region, and the Altai Republic. A majority (11 of 15) of CRF63_02A viruses from Kemerovo were grouped in a cluster. Antiretroviral (ARV) drug resistance mutations were found in 6 (14%) of 43 drug-naive patients. This study provides new insights in HIV-1 molecular epidemiology and prevalence of transmitted ARV drug resistance mutations in Southwestern Siberia.

Prevalence of Antiretroviral Drug Resistance Mutations Among Pretreatment and Antiretroviral Therapy-Failure HIV Patients in Uzbekistan

To evaluate the national prevalence of antiretroviral therapy (ART)-resistant HIV-1 viruses among both ART-initiators (pretreatment drug resistance, PDR) and ART-failure HIV patients in Uzbekistan. A nation-wide, cross-sectional active HIV-1 PDR surveillance was conducted in Uzbekistan from 2015 to 2016. In total, 713 blood plasma samples from adults were collected, including samples from ART-naive patients initiating ART and ART-failure HIV patients. HIV-1 genome polregion viral sequences were obtained from 309 patients, of those 106 on ART and 203 on ART-initiators. Analysis of HIV-1 subtypes and drug resistance mutations (DRMs) to HIV protease and reverse transcriptase inhibitors was performed. Among all the viruses studied, HIV-1 CRF 02_AG recombinant was the most common-57% (176/309). The second major group was represented by A1-40.5% (125/309). Two viruses were found to be recombinants formed by subtypes A1 and CRF02_AG sequences. ART-naive cohort I (PDR) included six samples that contained at least one surveillance drug resistance mutation (SDRM) (2.96%), with the most common being K103N mutation (4/6). In ART-experienced patients, cohort II, 77.4% (82/106) of viruses contained at least one mutation against PIs, NRTIs, or NNRTIs, with the most common mutations of M184V/I (49.1%; 52/106), K65R (18.9%; 20/106), K103N (23.6%; 25/106), and G190S (22.6%; 24/106). The significant difference in frequency of mutations was found between two dominant subtypes, A1 and CRF02_AG. The molecular epidemiological profile of HIV infection in Uzbekistan has changed toward a predominance of CRF02_AG viruses. In the first national-scale study of the PDR prevalence, it was found to be relatively low (2.96%). The DR mutations in failure patients correspond to the main therapy regimens (NRTI/NNRTI) adopted in the country. The observations provide new evidence for differences in ART efficacy and resistance profiles for different subtypes.

Global and regional epidemiology of HIV-1 recombinants in 1990-2015: a systematic review and global survey

BACKGROUND

Global HIV-1 genetic diversity and evolution form a major challenge to treatment and prevention efforts. An increasing number of distinct HIV-1 recombinants have been identified worldwide, but their contribution to the global epidemic is unknown. We aimed to estimate the global and regional distribution of HIV-1 recombinant forms during 1990-2015.

METHODS

We assembled a global HIV-1 molecular epidemiology database through a systematic literature review and a global survey. We searched the PubMed, Embase (Ovid), CINAHL (Ebscohost), and Global Health (Ovid) databases for HIV-1 subtyping studies published from Jan 1, 1990, to Dec 31, 2015. Unpublished original HIV-1 subtyping data were collected through a survey among experts in the field who were members of the WHO-UNAIDS Network for HIV Isolation and Characterisation. We included prevalence studies with HIV-1 subtyping data collected during 1990-2015. Countries were grouped into 14 regions and analyses were done for four time periods (1990-99, 2000-04, 2005-09, and 2010-15). The distribution of circulating recombinant forms (CRFs) and unique recombinant forms (URFs) in individual countries was weighted according to the UNAIDS estimates of the number of people living with HIV in each country to generate regional and global estimates of numbers and proportions of HIV-1 recombinants in each time period. The systematic review is registered with PROSPERO, CRD42017067164.

FINDINGS

Our global data collection yielded an HIV-1 molecular epidemiology database of 383 519 samples from 116 countries in 1990-2015. We found that the proportion of recombinants increased over time, both globally and in most regions, reaching 22·8% (7 978 517 of 34 921 639) of global HIV-1 infections in 2010-15. Both the proportion and the number of distinct CRFs detected increased over time to 16·7% and 57 CRFs in 2010-15. The global and regional distribution of HIV-1 recombinants was diverse and evolved over time, and we found large regional variation in the numbers (0-44 CRFs), types (58 distinct CRFs), and proportions (0-80·5%) of HIV-1 recombinants. Globally, CRF02_AG was the most prevalent recombinant, accounting for 33·9% (2 701 364 of 7 978 517) of all recombinant infections in 2010-15. URFs accounted for 26·7% (2 131 450 of 7 978 517), CRF01_AE for 23·0% (1 838 433), and other CRFs for 16·4% (1 307 270) of all recombinant infections in 2010-15. Although other CRFs accounted for small proportions of infections globally (<1% each), they were prominent in regional epidemics, including in east and southeast Asia, west and central Africa, Middle East and north Africa, and eastern Europe and central Asia. In addition, in 2010-15, central Africa (21·3% [243 041 of 1 143 531]), west Africa (15·5% [838 476 of 5 419 010]), east Africa (12·6% [591 140 of 4 704 986]), and Latin America (9·6% [153 069 of 1 586 605]) had high proportions of URFs.

INTERPRETATION

HIV-1 recombinants are increasingly prominent in global and regional HIV epidemics, which has important implications for the development of an HIV vaccine and the design of diagnostic, resistance, and viral load assays. Continued and improved surveillance of the global molecular epidemiology of HIV is crucial.

FUNDING

None.

Prevalence and spatiotemporal dynamics of HIV-1 Circulating Recombinant Form 03_AB (CRF03_AB) in the Former Soviet Union countries

Background

HIV-1 circulating recombinant forms (CRFs) infections has been increasing in Former Soviet Union (FSU) countries in the recent decade. One is the CRF03_AB, which circulated in the region since late 1990s and probably became widespread in northwestern FSU countries. However, there is not much information provided about the dissemination of this recombinant. Here, we examine the prevalence, evolutionary dynamics and dispersion pattern of HIV-1 CRF03_AB recombinant.

Methods

We analyzed 32 independent studies and 151 HIV-1 CRF03_AB pol sequences isolated from different FSU countries over a period of 22 years. Pooled prevalence was estimated using a random effects model. Bayesian coalescent-based method was used to estimate the evolutionary, phylogeographic and demographic parameters.

Results

Our meta-analysis showed that the pooled prevalence of CRF03_AB infection in northwestern FSU region was 5.9% [95%CI: 4.1–7.8]. Lithuania (11.6%), Russia (5.9%) and Belarus (2.9%) were the most affected by CRF03_AB. We found that early region wide spread of HIV-1 CRF03_AB originated from one viral clade that arose in the city of Kaliningrad in 1992 [95%HPD: 1990–1995]. Fourteen migration route of this variant were found. The city of Kaliningrad is involved in most of these, confirming its leading role in CRF03_AB spread within FSU. Demographic reconstruction point to this is that CRF03_AB clade seems to have experienced an exponential growth until the mid-2000s and a decrease in recent years.

Conclusion

These data provide new insights into the molecular epidemiology of CRF03_AB as well as contributing to the fundamental understanding of HIV epidemic in FSU.

SARS-CoV-2 exhibits intra-host genomic plasticity and low-frequency polymorphic quasispecies

In December 2019, an outbreak of atypical pneumonia (Coronavirus disease 2019 -COVID-19) associated with a novel coronavirus (SARS-CoV-2) was reported in Wuhan city, Hubei province, China. The outbreak was traced to a seafood wholesale market and human to human transmission was confirmed. The rapid spread and the death toll of the new epidemic warrants immediate intervention. The intra-host genomic variability of SARS-CoV-2 plays a pivotal role in the development of effective antiviral agents and vaccines, as well as in the design of accurate diagnostics. We analyzed NGS data derived from clinical samples of three Chinese patients infected with SARS-CoV-2, in order to identify small- and large-scale intra-host variations in the viral genome. We identified tens of low- or higher- frequency single nucleotide variations (SNVs) with variable density across the viral genome, affecting 7 out of 10 protein-coding viral genes. The majority of these SNVs (72/104) corresponded to missense changes. The annotation of the identified SNVs but also of all currently circulating strain variations revealed colocalization of intra-host as well as strain specific SNVs with primers and probes currently used in molecular diagnostics assays. Moreover, we de-novo assembled the viral genome, in order to isolate and validate intra-host structural variations and recombination breakpoints. The bioinformatics analysis disclosed genomic rearrangements over poly-A / poly-U regions located in ORF1ab and spike (S) gene, including a potential recombination hot-spot within S gene. Our results highlight the intra-host genomic diversity and plasticity of SARS-CoV-2, pointing out genomic regions that are prone to alterations. The isolated SNVs and genomic rearrangements reflect the intra-patient capacity of the polymorphic quasispecies, which may arise rapidly during the outbreak, allowing immunological escape of the virus, offering resistance to anti-viral drugs and affecting the sensitivity of the molecular diagnostics assays.

Why COVID-19 Transmission Is More Efficient and Aggressive Than Viral Transmission in Previous Coronavirus Epidemics?

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is causing a pandemic of coronavirus disease 2019 (COVID-19). The worldwide transmission of COVID-19 from human to human is spreading like wildfire, affecting almost every country in the world. In the past 100 years, the globe did not face a microbial pandemic similar in scale to COVID-19. Taken together, both previous outbreaks of other members of the coronavirus family (severe acute respiratory syndrome (SARS-CoV) and middle east respiratory syndrome (MERS-CoV)) did not produce even 1% of the global harm already inflicted by COVID-19. There are also four other CoVs capable of infecting humans (HCoVs), which circulate continuously in the human population, but their phenotypes are generally mild, and these HCoVs received relatively little attention. These dramatic differences between infection with HCoVs, SARS-CoV, MERS-CoV, and SARS-CoV-2 raise many questions, such as: Why is COVID-19 transmitted so quickly? Is it due to some specific features of the viral structure? Are there some specific human (host) factors? Are there some environmental factors? The aim of this review is to collect and concisely summarize the possible and logical answers to these questions.

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

BACKGROUND

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Near Full-length Genomic Sequencing and Molecular Analysis of HIV-Infected Individuals in a Network-based Intervention (TRIP) in Athens, Greece: Evidence that Transmissions Occur More Frequently from those with High HIV-RNA

Background:

TRIP (Transmission Reduction Intervention Project) was a network-based, contact tracing approach to locate and link to care, mostly people who inject drugs (PWID) with recent HIV infection.

Objective:

We investigated whether sequences from HIV-infected participants with high viral load cluster together more frequently than what is expected by chance.

Methods:

Paired end reads were generated for 104 samples using Illumina MiSeq next-generation se-quencing.

Results:

63 sequences belonged to previously identified local transmission networks of PWID (LTNs) of an HIV outbreak in Athens, Greece. For two HIV-RNA cut-offs (105 and 106 IU/mL), HIV transmissions were more likely between PWID with similar levels of HIV-RNA (p<0.001). 10 of the 14 sequences (71.4%) from PWID with HIV-RNA >106 IU/mL were clustered in 5 pairs. For 4 of these clusters (80%), there was in each one of them at least one sequence from a recently HIV-infected PWID.

Conclusion:

We showed that transmissions are more likely among PWID with high viremia.

Human Immunodeficiency Virus-1 Diversity in the Moscow Region, Russia: Phylodynamics of the Most Common Subtypes

This study analyzes the HIV-1 subtype diversity and its phylodynamics in Moscow region, which is the most densely populated area of Russia characterized by high rates of internal and external migration. The demographic and viral data from 896 HIV-infected individuals collected during 2011–2016 were analyzed. The study revealed broad diversity in the HIV-1 subtypes found in Moscow, which included A6 (85.1%), B (7.6%), CRF02_AG (1.2%) and URF_A6/B recombinants (4.2%). Other HIV-1 subtypes were detected as single cases. While A6 was most prevalent (>86.0%) among heterosexuals, injecting drug users and cases of mother-to-child transmission of HIV, subtype B (76.3%) was more common in men who have sex with men. Phylogenetic reconstruction revealed that the A6 sequences were introduced into the epidemic cluster that arose approximately around 1998. Within the subtype B, six major epidemic clusters were identified, each of which contained strains associated with only one or two dominant transmission routes. The date of origin of these clusters varied between 1980 and 1993, indicating that the HIV-1 B epidemic began much earlier than the HIV-1 A6 epidemic. Reconstruction of the demographic history of subtypes A6 and B identified at least two epidemic growth phases, which included an initial phase of exponential growth followed by a decline in the mid/late 2010s. Thus, our results indicate an increase in HIV-1 genetic diversity in Moscow region. They also help in understanding the HIV-1 temporal dynamics as well as the genetic relationships between its circulating strains.