Molecular epidemiology of HIV-1 in St Petersburg, Russia: predominance of subtype A, former Soviet Union variant, and identification of intrasubtype subclusters

Virally Suppressed People Living with HIV Who Use Opioids Have Diminished Latency Reversal

Of the 12 million people who inject drugs worldwide, 13% live with HIV. Whether opioid use impacts HIV pathogenesis and latency is an outstanding question. To gain insight into whether opioid use influences the proviral landscape and latent HIV reservoir, we performed intact proviral DNA assays (IPDA) on peripheral blood mononuclear cells (PBMCs) from antiretroviral therapy (ART)-suppressed people living with HIV (PWH) with or without current opioid use. No differences were observed between PWH with and without opioid use in the frequency of HIV intact and defective proviral genomes. To evaluate the latent reservoir, we activated PBMCs from ART-suppressed PWH with or without opioid use and assessed the induction of HIV RNA. PWH using opioids had diminished responses to ex vivo HIV reactivation, suggesting a smaller reversible reservoir of HIV-1 latently infected cells. However, in vitro studies using primary CD4+ T cells treated with morphine showed no effect of opioids on HIV-1 infection, replication or latency establishment. The discrepancy in our results from in vitro and clinical samples suggests that while opioids may not directly impact HIV replication, latency and reactivation in CD4+ T cells, opioid use may indirectly shape the HIV reservoir in vivo by modulating general immune functions.

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

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 A_FSU) 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.

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.

Molecular Epidemiology of HIV-1 in Oryol Oblast, Russia

The HIV/AIDS epidemic in Russia is growing, with approximately 100,000 people infected annually. Molecular epidemiology can provide insight into the structure and dynamics of the epidemic. However, its applicability in Russia is limited by the weakness of genetic surveillance, as viral genetic data are only available for <1 per cent of cases. Here, we provide a detailed description of the HIV-1 epidemic for one geographic region of Russia, Oryol Oblast, by collecting and sequencing viral samples from about a third of its known HIV-positive population (768 out of 2,157 patients). We identify multiple introductions of HIV-1 into Oryol Oblast, resulting in eighty-two transmission lineages that together comprise 66 per cent of the samples. Most introductions are of subtype A (315/332), the predominant HIV-1 subtype in Russia, followed by CRF63 and subtype B. Bayesian analysis estimates the effective reproduction number R_e for subtype A at 2.8 [1.7–4.4], in line with a growing epidemic. The frequency of CRF63 has been growing more rapidly, with the median R_e of 11.8 [4.6–28.7], in agreement with recent reports of this variant rising in frequency in some regions of Russia. In contrast to the patterns described previously in European and North American countries, we see no overrepresentation of males in transmission lineages; meanwhile, injecting drug users are overrepresented in transmission lineages. This likely reflects the structure of the HIV-1 epidemic in Russia dominated by heterosexual and, to a smaller extent, people who inject drugs transmission. Samples attributed to men who have sex with men (MSM) transmission are associated with subtype B and are less prevalent than expected from the male-to-female ratio for this subtype, suggesting underreporting of the MSM transmission route. Together, our results provide a high-resolution description of the HIV-1 epidemic in Oryol Oblast, Russia, characterized by frequent interregional transmission, rapid growth of the epidemic, and rapid displacement of subtype A with the recombinant CRF63 variant.

Dating the Origin and Estimating the Transmission Rates of the Major HIV-1 Clusters in Greece: Evidence about the Earliest Subtype A1 Epidemic in Europe

Our aim was to estimate the date of the origin and the transmission rates of the major local clusters of subtypes A1 and B in Greece. Phylodynamic analyses were conducted in 14 subtype A1 and 31 subtype B clusters. The earliest dates of origin for subtypes A1 and B were in 1982.6 and in 1985.5, respectively. The transmission rate for the subtype A1 clusters ranged between 7.54 and 39.61 infections/100 person years (IQR: 9.39, 15.88), and for subtype B clusters between 4.42 and 36.44 infections/100 person years (IQR: 7.38, 15.04). Statistical analysis revealed that the average difference in the transmission rate between the PWID and the MSM clusters was 6.73 (95% CI: 0.86 to 12.60; p = 0.026). Our study provides evidence that the date of introduction of subtype A1 in Greece was the earliest in Europe. Transmission rates were significantly higher for PWID than MSM clusters due to the conditions that gave rise to an extensive PWID HIV-1 outbreak ten years ago in Athens, Greece. Transmission rate can be considered as a valuable measure for public health since it provides a proxy of the rate of epidemic growth within a cluster and, therefore, it can be useful for targeted HIV prevention programs.

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.

Genetic Diversity of HIV-1 in Krasnoyarsk Krai: Area with High Levels of HIV-1 Recombination in Russia

More than a quarter of HIV-infected individuals registered in Russia live in Siberia. Unlike Central Russia where HIV-1 subtype A6 is predominant, in most Siberian regions since 2012, a new HIV-1 CRF63_02A1 genetic variant has spread, with the share of this variant attaining 75-85% among newly identified HIV cases. Krasnoyarsk Krai is considered to be a high-risk territory according to morbidity rate and HIV infection incidence among the population. The current paper aims to study the molecular epidemiologic characteristics of HIV-1 spreading in Krasnoyarsk Krai. Phylogenetic and recombination analyses of pol (PR-RT, IN) and env regions of the virus were used for genotyping 159 HIV-1 isolated in Krasnoyarsk Krai. 57.2% of the isolates belonged to subtype A (A6) specific to Russia, 12.6% to CRF63_02A1, and 0.6% to CRF02_AGСА, and in 29.6% HIV-1 URFs were detected, including URF63/А (23.9%), URFА/В (4.4%), and URF02/А (1.3%). In 6 of 7, HIV-1 URFА/В identical recombination model was detected; the origin of 38 URF63/А was proven to be the result of individual recombination events. Since 2015, a share of the population with newly diagnosed HIV who were infected with HIV-1 URF reached an exceptionally high rate of 38.6%. As distinct from adjacent Siberian regions, the HIV-1 CRF63_02A1 prevalence rate in Krasnoyarsk Krai is within 16%; however, the increased contribution of new HIV-1 into the regional epidemic development was observed due to the recombination of viruses of subtypes А, В, and CRF63_02A1. The difference between the described molecular epidemiologic picture in Krasnoyarsk Krai and in adjacent areas is likely caused by differences in predominant routes of HIV transmission and by more recent HIV-1 CRF63_02A1 transmission in the PWID group, which had a high prevalence of HIV-1 subtype A by the time of the new virus transmission, resulting in increased possibility of coinfection with various HIV-1 genetic variants.

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.

A cross-sectional study to characterize local HIV-1 dynamics in Washington, DC using next-generation sequencing

Washington, DC continues to experience a generalized HIV-1 epidemic. We characterized the local phylodynamics of HIV-1 in DC using next-generation sequencing (NGS) data. Viral samples from 68 participants from 2016 through 2017 were sequenced and paired with epidemiological data. Phylogenetic and network inferences, drug resistant mutations (DRMs), subtypes and HIV-1 diversity estimations were completed. Haplotypes were reconstructed to infer transmission clusters. Phylodynamic inferences based on the HIV-1 polymerase (pol) and envelope genes (env) were compared. Higher HIV-1 diversity (n.s.) was seen in men who have sex with men, heterosexual, and male participants in DC. 54.0% of the participants contained at least one DRM. The 40-49 year-olds showed the highest prevalence of DRMs (22.9%). Phylogenetic analysis of pol and env sequences grouped 31.9-33.8% of the participants into clusters. HIV-TRACE grouped 2.9-12.8% of participants when using consensus sequences and 9.0-64.2% when using haplotypes. NGS allowed us to characterize the local phylodynamics of HIV-1 in DC more broadly and accurately, given a better representation of its diversity and dynamics. Reconstructed haplotypes provided novel and deeper phylodynamic insights, which led to networks linking a higher number of participants. Our understanding of the HIV-1 epidemic was expanded with the powerful coupling of HIV-1 NGS data with epidemiological data.

The phylogenetic approach for viral infectious disease evolution and epidemiology: An updating review

In the last decade, the phylogenetic approach is recurrent in molecular evolutionary analysis. On 12 May, 2019, about 2 296 213 papers are found, but typing "phylogeny" or "epidemiology AND phylogeny" only 199 804 and 20 133 are retrieved, respectively. Molecular epidemiology in infectious diseases is widely used to define the source of infection as so as the ancestral relationships of individuals sampled from a population. Coalescent theory and phylogeographic analysis have had scientific application in several, recent pandemic events, and nosocomial outbreaks. Hepatitis viruses and immunodeficiency virus (human immunodeficiency virus) have been largely studied. Phylogenetic analysis has been recently applied on Polyomaviruses so as in the more recent outbreaks due to different arboviruses type as Zika and chikungunya viruses discovering the source of infection and the geographic spread. Data on sequences isolated by the microorganism are essential to apply the phylogenetic tools and research in the field of infectious disease phylodinamics is growing up. There is the need to apply molecular phylogenetic and evolutionary methods in areas out of infectious diseases, as translational genomics and personalized medicine. Lastly, the application of these tools in vaccine strategy so as in antibiotic and antiviral researchers are encouraged.

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.

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

HIV-1 protease-reverse transcriptase sequences from 62 HIV-1-infected individuals recently diagnosed in Moscow were analyzed. Subtype A former Soviet Union (FSU) (A_FSU) variant was the predominant clade (62.9%), followed by subtype B (22.6%), unique recombinants (6.5%), subtype G (6.5%), and CRF01_AE (1.6%). A_FSU predominated among people who inject drugs (88.9%) and heterosexually acquired infections (77.8%), while subtype B was the most prevalent genetic form among men who have sex with men (44%), although A_FSU was also frequent in this population (36%). Forty-eight (77.4%) viruses branched within intrasubtype clusters, three of which, of subtype B, had a majority of viruses collected outside of FSU. The four subtype G viruses identified in this study belonged to the Portuguese-Spanish (Iberian) variant and, together with three from databases, formed a Russian cluster closely related to viruses from Denmark. This is the first report of the circulation of the Iberian subtype G variant in Russia.

Molecular epidemiology of HIV-1 subtype A in former Soviet Union countries

While in other parts of the world it is on decline, incidence of HIV infection continues to rise in the former Soviet Union (FSU) countries. The present study was conducted to investigate the patterns and modes of HIV transmission in FSU countries. We performed phylogenetic analysis of publicly available 2705 HIV-1 subtype A pol sequences from thirteen FSU countries: Armenia, Azerbaijan, Belarus, Estonia, Georgia, Kazakhstan, Kyrgyzstan, Latvia, Lithuania, Moldova, Russia, Ukraine and Uzbekistan. Our analysis showed that the clusters from FSU countries were intermixed, indicating a possible role of transmigration in HIV transmission. Injection drug use was found to be the most frequent mode of transmission, while the clusters from PWID and heterosexual transmission were intermixed, indicating bridging of HIV infection across populations. To control the expanding HIV epidemic in this region, harm reduction strategies should be focused on three modes of transmission, namely, cross-border migration, injection drug use and heterosexual.

Characterization of HIV diversity, phylodynamics and drug resistance in Washington, DC

BACKGROUND

Washington DC has a high burden of HIV with a 2.0% HIV prevalence. The city is a national and international hub potentially containing a broad diversity of HIV variants; yet few sequences from DC are available on GenBank to assess the evolutionary history of HIV in the US capital. Towards this general goal, here we analyze extensive sequence data and investigate HIV diversity, phylodynamics, and drug resistant mutations (DRM) in DC.

METHODS

Molecular HIV-1 sequences were collected from participants infected through 2015 as part of the DC Cohort, a longitudinal observational study of HIV+ patients receiving care at 13 DC clinics. Sequences were paired with Cohort demographic, risk, and clinical data and analyzed using maximum likelihood, Bayesian and coalescent approaches of phylogenetic, network and population genetic inference. We analyzed 601 sequences from 223 participants for int (~864 bp) and 2,810 sequences from 1,659 participants for PR/RT (~1497 bp).

RESULTS

Ninety-nine and 94% of the int and PR/RT sequences, respectively, were identified as subtype B, with 14 non-B subtypes also detected. Phylodynamic analyses of US born infected individuals showed that HIV population size varied little over time with no significant decline in diversity. Phylogenetic analyses grouped 13.5% of the int sequences into 14 clusters of 2 or 3 sequences, and 39.0% of the PR/RT sequences into 203 clusters of 2-32 sequences. Network analyses grouped 3.6% of the int sequences into 4 clusters of 2 sequences, and 10.6% of the PR/RT sequences into 76 clusters of 2-7 sequences. All network clusters were detected in our phylogenetic analyses. Higher proportions of clustered sequences were found in zip codes where HIV prevalence is highest (r = 0.607; P<0.00001). We detected a high prevalence of DRM for both int (17.1%) and PR/RT (39.1%), but only 8 int and 12 PR/RT amino acids were identified as under adaptive selection. We observed a significant (P<0.0001) association between main risk factors (men who have sex with men and heterosexuals) and genotypes in the five well-supported clusters with sufficient sample size for testing.

DISCUSSION

Pairing molecular data with clinical and demographic data provided novel insights into HIV population dynamics in Washington, DC. Identification of populations and geographic locations where clustering occurs can inform and complement active surveillance efforts to interrupt HIV transmission.

The global spread of HIV-1 subtype B epidemic

Human immunodeficiency virus type 1 (HIV-1) was discovered in the early 1980s when the virus had already established a pandemic. For at least three decades the epidemic in the Western World has been dominated by subtype B infections, as part of a sub-epidemic that traveled from Africa through Haiti to United States. However, the pattern of the subsequent spread still remains poorly understood. Here we analyze a large dataset of globally representative HIV-1 subtype B strains to map their spread around the world over the last 50years and describe significant spread patterns. We show that subtype B travelled from North America to Western Europe in different occasions, while Central/Eastern Europe remained isolated for the most part of the early epidemic. Looking with more detail in European countries we see that the United Kingdom, France and Switzerland exchanged viral isolates with non-European countries than with European ones. The observed pattern is likely to mirror geopolitical landmarks in the post-World War II era, namely the rise and the fall of the Iron Curtain and the European colonialism. In conclusion, HIV-1 spread through specific migration routes which are consistent with geopolitical factors that affected human activities during the last 50years, such as migration, tourism and trade. Our findings support the argument that epidemic control policies should be global and incorporate political and socioeconomic factors.

Molecular epidemiology of HIV-1 infection in Europe: An overview

Human Immunodeficiency Virus type 1 (HIV-1) is characterised by vast genetic diversity. Globally circulating HIV-1 viruses are classified into distinct phylogenetic strains (subtypes, sub-subtypes) and several recombinant forms. Here we describe the characteristics and evolution of European HIV-1 epidemic over time through a review of published literature and updated queries of existing HIV-1 sequence databases. HIV-1 in Western and Central Europe was introduced in the early-1980s in the form of subtype B, which is still the predominant clade. However, in Eastern Europe (Former Soviet Union (FSU) countries and Russia) the predominant strain, introduced into Ukraine in the mid-1990s, is subtype A (A_FSU) with transmission mostly occurring in People Who Inject Drugs (PWID). In recent years, the epidemic is evolving towards a complex tapestry with an increase in the prevalence of non-B subtypes and recombinants in Western and Central Europe. Non-B epidemics are mainly associated with immigrants, heterosexuals and females but more recently, non-B clades have also spread amongst groups where non-B strains were previously absent - non-immigrant European populations and amongst men having sex with men (MSM). In some countries, non-B clades have spread amongst the native population, for example subtype G in Portugal and subtype A in Greece, Albania and Cyprus. Romania provides a unique case where sub-subtype F1 has predominated throughout the epidemic. In contrast, HIV-1 epidemic in FSU countries remains more homogeneous with A_FSU clade predominating in all countries. The differences between the evolution of the Western epidemic and the Eastern epidemic may be attributable to differences in transmission risk behaviours, lifestyle and the patterns of human mobility. The study of HIV-1 epidemic diversity provides a useful tool by which we can understand the history of the pandemic in addition to allowing us to monitor the spread and growth of the epidemic over time.

HIV-1 A1 Subtype Epidemic in Italy Originated from Africa and Eastern Europe and Shows a High Frequency of Transmission Chains Involving Intravenous Drug Users

Background

Subtype A accounts for only 12% of HIV-1 infections worldwide but predominates in Russia and Former Soviet Union countries of Eastern Europe. After an early propagation via heterosexual contacts, this variant spread explosively among intravenous drug users. A distinct A1 variant predominates in Greece and Albania, which penetrated directly from Africa. Clade A1 accounts for 12.5% of non-B subtypes in Italy, being the most frequent after F1 subtype.

Aim

Aim of this study was to investigate the circulation of A1 subtype in Italy and trace its origin and diffusion through phylogenetic and phylodynamic approaches.

Results

The phylogenetic analysis of 113 A1 pol sequences included in the Italian ARCA database, indicated that 71 patients (62.8%) clustered within 5 clades. A higher probability to be detected in clusters was found for patients from Eastern Europe and Italy (88.9% and 60.4%, respectively) compared to those from Africa (20%) (p < .001). Higher proportions of clustering sequences were found in intravenous drug users with respect to heterosexuals (85.7% vs. 59.3%, p = .056) and in women with respect to men (81.4% vs. 53.2%, p < .006). Subtype A1 dated phylogeny indicated an East African origin around 1961. Phylogeographical reconstruction highlighted 3 significant groups. One involved East European and some Italian variants, the second encompassed some Italian and African strains, the latter included the majority of viruses carried by African and Italian subjects and all viral sequences from Albania and Greece.

Conclusions

Subtype A1 originated in Central Africa and spread among East European countries in 1982. It entered Italy through three introduction events: directly from East Africa, from Albania and Greece, and from the area encompassing Moldavia and Ukraine. As in previously documented A1 epidemics of East European countries, HIV-1 A1 subtype spread in Italy in part through intravenous drug users. However, Eastern European women contributed to the penetration of such variant, probably through sex work.

High recombination potential of subtype A HIV-1

Recombination can assort polymorphic alleles to increase diversity in the HIV-1 population. To better understand the recombination potential of subtype A HIV-1, we generated viruses containing sequences from two variants circulating in Russia and analyzed the polymerase gene (pol) of the recombinants after one round of HIV-1 replication using single-genome sequencing. We observed that recombination occurred throughout pol and could easily assort alleles containing mutations that conferred resistance to currently approved antivirals. We measured the recombination rate in various regions of pol including a G-rich region that has been previously proposed to be a recombination hot spot. Our study does not support a recombination hot spot in this G-rich region. Importantly, of the 58 proviral sequences containing crossover event(s) in pol, we found that each sequence was a unique genotype indicating that recombination is a powerful genetic mechanism in assorting the genomes of subtype A HIV-1 variants.

Epidemiological Surveillance of HIV-1 Transmitted Drug Resistance in Spain in 2004-2012: Relevance of Transmission Clusters in the Propagation of Resistance Mutations

Our objectives were to carry out an epidemiological surveillance study on transmitted drug resistance (TDR) among individuals newly diagnosed of HIV-1 infection during a nine year period in Spain and to assess the role of transmission clusters (TC) in the propagation of resistant strains. An overall of 1614 newly diagnosed individuals were included in the study from January 2004 through December 2012. Individuals come from two different Spanish regions: Galicia and the Basque Country. Resistance mutations to reverse transcriptase inhibitors (RTI) and protease inhibitors (PI) were analyzed according to mutations included in the surveillance drug-resistance mutations list updated in 2009. TC were defined as those comprising viruses from five or more individuals whose sequences clustered in maximum likelihood phylogenetic trees with a bootstrap value ≥90%. The overall prevalence of TDR to any drug was 9.9%: 4.9% to nucleoside RTIs (NRTIs), 3.6% to non-nucleoside RTIs (NNRTIs), and 2.7% to PIs. A significant decrease of TDR to NRTIs over time was observed [from 10% in 2004 to 2% in 2012 (p=0.01)]. Sixty eight (42.2%) of 161 sequences with TDR were included in 25 TC composed of 5 or more individuals. Of them, 9 clusters harbored TDR associated with high level resistance to antiretroviral drugs. T215D revertant mutation was transmitted in a large cluster comprising 25 individuals. The impact of epidemiological networks on TDR frequency may explain its persistence in newly diagnosed individuals. The knowledge of the populations involved in TC would facilitate the design of prevention programs and public health interventions.

Short communication: Molecular epidemiology, phylogeny, and phylodynamics of CRF63_02A1, a recently originated HIV-1 circulating recombinant form spreading in Siberia

The HIV-1 epidemic in Russia is dominated by the former Soviet Union subtype A (A(FSU)) variant, but other genetic forms are circulating in the country. One is the recently described CRF63_02A1, derived from recombination between a CRF02_AG variant circulating in Central Asia and A(FSU), which has spread in the Novosibirsk region, Siberia. Here we phylogenetically analyze pol and env segments from 24 HIV-1 samples from the Novosibirsk region collected in 2013, with characterization of three new near full-length genome CRF63_02A1 sequences, and estimate the time of the most recent common ancestor (tMRCA) and the demographic growth of CRF63_02A1 using a Bayesian method. The analyses revealed that CRF63_02A1 is highly predominant in the Novosibirsk region (81.2% in pol sequences) and is transmitted both among injecting drug users and by heterosexual contact. Similarity searches with database sequences combined with phylogenetic analyses show that CRF63_02A1 is circulating in East Kazakhstan and the Eastern area of Russia bordering China. The analyses of near full-length genome sequences show that its mosaic structure is more complex than reported, with 18 breakpoints. The tMRCA of CRF63_02A1 was estimated around 2006, with exponential growth in 2008-2009 and subsequent stabilization. These results provide new insights into the molecular epidemiology, phylogeny, and phylodynamics of CRF63_02A1.

Phylogeographic analysis on the travel-related introduction of HIV-1 non-B subtypes to Northern Poland

Phylodynamic, sequence data based reconstructions for the surveillance of the geographic spatial spread are a powerful tool in molecular epidemiology. In this study region of origin for the set of 57 partial pol sequences derived from the patients the history of travel-related HIV transmission was analyzed using phylogeographic approach. Maximum likelihood trees based on the sets of country-annotated reference sequences were inferred for identified non-B variants. Region of sequence import was assigned using on the highest approximate likelihood ratios. Import of the A1 clades was traced to the Eastern Europe and associated with immigration from this region. Subtype C infections clustered most frequently with sequences of the South African origin while majority of subtype Ds were similar to the European clades. Subtype G sequences clustered with Portuguese lineage, CRF01_AE with Eastern or South-Eastern Asian. Eastern European, Middle African or Western African lineage was assigned for the CFR02_AG. Rare circulating recombinants originated either from Central Africa (CRF11_cpx - Democratic Republic of Congo, CRF13_cpx - Central African Republic, CRF37_cpx - Cameroon) or South America (CRF28_BF and CRF46_BF - Brazil). Import of the HIV-1 non-B variants, including recombinant forms previously rarely found in Poland and Europe is frequent among travelers. Observed founder events result in the heterosexually-driven introduction of the novel HIV-1 variants into the population.

Consensus HIV-1 subtype A integrase and its raltegravir-resistant variants: design and characterization of the enzymatic properties

Model studies of the subtype B and non-subtype B integrases are still required to compare their susceptibility to antiretroviral drugs, evaluate the significance of resistance mutations and identify the impact of natural polymorphisms on the level of enzymatic reactivity. We have therefore designed the consensus integrase of the HIV-1 subtype A strain circulating in the former Soviet Union territory (FSU-A) and two of its variants with mutations of resistance to the strand transfer inhibitor raltegravir. Their genes were synthesized, and expressed in E coli; corresponding His-tagged proteins were purified using the affinity chromatography. The enzymatic properties of the consensus integrases and their sensitivity to raltegravir were examined in a series of standard in vitro reactions and compared to the properties of the integrase of HIV-1 subtype B strain HXB2. The consensus enzyme demonstrated similar DNA-binding properties, but was significantly more active than HXB-2 integrase in the reactions of DNA cleavage and integration. All integrases were equally susceptible to inhibition by raltegravir and elvitegravir, indicating that the sporadic polymorphisms inherent to the HXB-2 enzyme have little effect on its susceptibility to drugs. Insensitivity of the mutated enzymes to the inhibitors of strand transfer occurred at a cost of a 30-90% loss of the efficacies of both 3'-processing and strand transfer. This is the first study to describe the enzymatic properties of the consensus integrase of HIV-1 clade A and the effects of the resistance mutations when the complex actions of sporadic sequence polymorphisms are excluded.

Consensus HIV-1 FSU-A integrase gene variants electroporated into mice induce polyfunctional antigen-specific CD4+ and CD8+ T cells

Our objective is to create gene immunogens targeted against drug-resistant HIV-1, focusing on HIV-1 enzymes as critical components in viral replication and drug resistance. Consensus-based gene vaccines are specifically fit for variable pathogens such as HIV-1 and have many advantages over viral genes and their expression-optimized variants. With this in mind, we designed the consensus integrase (IN) of the HIV-1 clade A strain predominant in the territory of the former Soviet Union and its inactivated derivative with and without mutations conferring resistance to elvitegravir. Humanized IN gene was synthesized; and inactivated derivatives (with 64D in the active site mutated to V) with and without elvitegravir-resistance mutations were generated by site-mutagenesis. Activity tests of IN variants expressed in E coli showed the consensus IN to be active, while both D64V-variants were devoid of specific activities. IN genes cloned in the DNA-immunization vector pVax1 (pVaxIN plasmids) were highly expressed in human and murine cell lines (>0.7 ng/cell). Injection of BALB/c mice with pVaxIN plasmids followed by electroporation generated potent IFN-γ and IL-2 responses registered in PBMC by day 15 and in splenocytes by day 23 after immunization. Multiparametric FACS demonstrated that CD8+ and CD4+ T cells of gene-immunized mice stimulated with IN-derived peptides secreted IFN-γ, IL-2, and TNF-α. The multi-cytokine responses of CD8+ and CD4+ T-cells correlated with the loss of in vivo activity of the luciferase reporter gene co-delivered with pVaxIN plasmids. This indicated the capacity of IN-specific CD4+ and CD8+ T-cells to clear IN/reporter co-expressing cells from the injection sites. Thus, the synthetic HIV-1 clade A integrase genes acted as potent immunogens generating polyfunctional Th1-type CD4+ and CD8+ T cells. Generation of such response is highly desirable for an effective HIV-1 vaccine as it offers a possibility to attack virus-infected cells via both MHC class I and II pathways.

The analysis of near full-length genome sequences of HIV type 1 subtype A viruses from Russia supports the monophyly of major intrasubtype clusters

The HIV-1 epidemic in Russia has been insufficiently studied, with only 11 complete genome sequences from this country currently available, only three of which are of the locally predominant genetic form, the former Soviet Union (FSU) subtype A variant (A(FSU)). Here we analyze 10 newly derived A(FSU) near full-length genome sequences from Russia. Samples were selected based on phylogenetic clustering in protease-reverse transcriptase in two of the major A(FSU) clusters, V77I(PR) (n=6), widely circulating in Russia and other FSU countries, and A(SP1) (n=4), predominant in St. Petersburg. The phylogenetic analysis shows that the V77I(PR) genomes group in a monophyletic cluster together with 10 previously obtained A(FSU) genome sequences from Uzbekistan, Kazakhstan, Russia, and Cyprus, all bearing the V77I substitution in protease. Similarly, the four A(SP1) genomes group in a monophyletic cluster. These results therefore show that the monophyly of V77I(PR) and A(SP1) A(FSU) clusters is supported in near complete genomes.

Isolation, characterization, and sequence analysis of a novel HIV type 1 subtype A1 strain from Russia

A new syncytia-inducing human immunodeficiency virus type 1 (HIV-1) strain designated PokA79 was first isolated from a peripheral blood mononuclear cells (PBMC) pool from 24 HIV-1-infected patients from Russia and adapted to growth in MT2, MT4, and U937 cell cultures. The full-length genome of this isolate was amplified by RT-two-round PCR and sequenced. The obtained sequence belongs to the A1 subtype and clustered with other A1 former Soviet Union (FSU) strains. Full-length genome analysis showed that some of the PokA79 genes contain insertions and frame shifts, and the gp120 V3 loop has the three amino acid insertion inside its apex, which is rare for subtype A1 sequences. Nonidentical insertions of the same size are also present in the V3 apex region of several HIV-1 sequences from Central Africa.

Phylogenetic structure in African HIV-1 subtype C revealed by selective sequential pruning

Subtype C is the most prevalent clade in the HIV-1 pandemic. Previous studies suggested that African HIV-1 subtype C lacks a well-defined phylogenetic structure. Here we show that, by sequential pruning of ambiguously positioned taxa, a well-defined intrasubtype C phylogenetic structure becomes apparent, with 52% African HIV-1 subtype C isolates analyzed in envelope sequences branching within 11 clusters, also supported in a tree of full-length genomes, and all with geographical associations. Among 46 viruses recently transmitted in South Africa, 70% branched within 7 clusters (41% in the largest one) and 15% additional isolates were intercluster recombinants. Choice of the outgroup sequence and inclusion of intrasubtype recombinant viruses in the analyses could greatly affect support of clusters. The identification of clusters comprising a large proportion of African HIV-1 subtype C viruses may have implications for the design of vaccines intended for use in areas where subtype C is prevalent.

Phylodynamics of HIV-1 from a phase III AIDS vaccine trial in Bangkok, Thailand

BACKGROUND

In 2003, a phase III placebo-controlled trial (VAX003) was completed in Bangkok, Thailand. Of the 2,546 individuals enrolled in the trial based on high risk for infection through injection drug use (IDU), we obtained clinical samples and HIV-1 sequence data (envelope glycoprotein gene gp120) from 215 individuals who became infected during the trial. Here, we used these data in combination with other publicly available gp120 sequences to perform a molecular surveillance and phylodynamic analysis of HIV-1 in Thailand.

METHODOLOGY AND FINDINGS

Phylogenetic and population genetic estimators were used to assess HIV-1 gp120 diversity as a function of vaccination treatment, viral load (VL) and CD4(+) counts, to identify transmission clusters and to investigate the timescale and demographics of HIV-1 in Thailand. Three HIV-1 subtypes were identified: CRF01_AE (85% of the infections), subtype B (13%) and CRF15_AE (2%). The Bangkok IDU cohort showed more gp120 diversity than other Asian IDU cohorts and similar diversity to that observed in sexually infected individuals. Moreover, significant differences (P<0.02) in genetic diversity were observed in CRF01_AE IDU with different VL and CD4(+) counts. No phylogenetic structure was detected regarding any of the epidemiological and clinical factors tested, although high proportions (35% to 50%) of early infections fell into clusters, which suggests that transmission chains associated with acute infection play a key role on HIV-1 spread among IDU. CRF01_AE was estimated to have emerged in Thailand in 1984.5 (1983-1986), 3-6 years before the first recognition of symptomatic patients (1989). The relative genetic diversity of the HIV-1 population has remained high despite decreasing prevalence rates since the mid 1990s.

CONCLUSIONS

Our study and recent epidemiological reports indicate that HIV-1 is still a major threat in Thailand and suggest that HIV awareness and prevention needs to be strengthened to avoid AIDS resurgence.

Epidemiological study of phylogenetic transmission clusters in a local HIV-1 epidemic reveals distinct differences between subtype B and non-B infections

Background

The number of HIV-1 infected individuals in the Western world continues to rise. More in-depth understanding of regional HIV-1 epidemics is necessary for the optimal design and adequate use of future prevention strategies. The use of a combination of phylogenetic analysis of HIV sequences, with data on patients' demographics, infection route, clinical information and laboratory results, will allow a better characterization of individuals responsible for local transmission.

Methods

Baseline HIV-1 pol sequences, obtained through routine drug-resistance testing, from 506 patients, newly diagnosed between 2001 and 2009, were used to construct phylogenetic trees and identify transmission-clusters. Patients' demographics, laboratory and clinical data, were retrieved anonymously. Statistical analysis was performed to identify subtype-specific and transmission-cluster-specific characteristics.

Results

Multivariate analysis showed significant differences between the 59.7% of individuals with subtype B infection and the 40.3% non-B infected individuals, with regard to route of transmission, origin, infection with Chlamydia (p = 0.01) and infection with Hepatitis C virus (p = 0.017). More and larger transmission-clusters were identified among the subtype B infections (p < 0.001). Overall, in multivariate analysis, clustering was significantly associated with Caucasian origin, infection through homosexual contact and younger age (all p < 0.001). Bivariate analysis additionally showed a correlation between clustering and syphilis (p < 0.001), higher CD4 counts (p = 0.002), Chlamydia infection (p = 0.013) and primary HIV (p = 0.017).

Conclusions

Combination of phylogenetics with demographic information, laboratory and clinical data, revealed that HIV-1 subtype B infected Caucasian men-who-have-sex-with-men with high prevalence of sexually transmitted diseases, account for the majority of local HIV-transmissions. This finding elucidates observed epidemiological trends through molecular analysis, and justifies sustained focus in prevention on this high risk group.

Short communication: Molecular epidemiology of HIV type 1 in the Republic of Dagestan, Russian Federation: virtually uniform circulation of subtype A, former Soviet Union variant, with predominance of the V77I(PR) subvariant

We examine the distribution of viral genetic forms and the presence of antiretroviral drug resistance mutations in HIV-1 infections in the Republic of Dagestan, in the North Caucasus area of Russia, where a recent large increase in HIV-1 infections has been documented. Samples were collected from 41 HIV-1-infected individuals from Dagestan, most of them from the cities of Derbent (n = 21) and Mahachkala (n = 18). Thirty six were injecting drug users and five were infected by heterosexual contact. None was on antiretroviral drug treatment. HIV-1 protease and a segment of reverse transcriptase were amplified by RT-PCR from plasma RNA and sequenced, and phylogenetic trees were constructed via maximum likelihood. Forty (97.6%) of 41 samples were of subtype A, former Soviet Union variant (A(FSU)), of which 27 (67.5%) clustered with the subvariant containing the V77I substitution in protease (V77I(PR)). Within this cluster, 13 viruses formed a local subcluster, 10 of which were from Derbent. Four viruses clustered with the A(SP2) subcluster, recently identified in St. Petersburg, two with a virus from Georgia and one with a virus from Azerbaijan. No mutations associated with antiretroviral drug resistance were detected. The results, therefore, show the relationship of the HIV-1 epidemic in Dagestan with that of other areas of Russia and of neighboring countries, and reveal the spread of the A(FSU) V77I(PR) variant in the North Caucasus area.

Near full-length genome characterization of a newly identified HIV type 1 subtype F variant circulating in St. Petersburg, Russia

We report the near full-length genome characterization of an HIV-1 subtype F virus (D88_845) collected in St. Petersburg, Russia, from a 25-year-old Russian woman perinatally infected in 1982. In a Bayesian phylogenetic analysis, the genome sequence branched basally to the subsubtype F1 clade. In partial sequences, D88_845 clustered with 13 other subtype F sequences from Russia, corresponding to gag (n = 2), pol (n = 3), and env (n = 8) segments. At least 11 of these sequences are from samples collected in St. Petersburg from heterosexually infected Russian individuals. In each of these segments, the Russian viruses formed a monophyletic cluster that branched as a sister clade of the F1 subsubtype. One sequence from Belgium branched with D88_845 with a posterior probability of 0.99. This is the first report on the identification and near full-length genome characterization of the subtype F variant circulating in St. Petersburg, which is closely related to, but distinct from, the F1 subsubtype.