Phylodynamics of the major HIV-1 CRF02_AG African lineages and its global dissemination

Mapping Transmission Dynamics and Drug Resistance Surveillance in the Cyprus HIV-1 Epidemic (2017-2021)

The human immunodeficiency virus type 1 (HIV-1) epidemic has been a major public health threat on a global scale since the early 1980s. Despite the introduction of combination antiretroviral therapy (cART), the incidence of new HIV-1 infections continues to rise in some regions around the world. Thus, with the continuous transmission of HIV-1 and the lack of a cure, it is imperative for molecular epidemiological studies to be performed, to monitor the infection and ultimately be able to control the spread of this virus. This work provides a comprehensive molecular epidemiological analysis of the HIV-1 infection in Cyprus, through examining 305 HIV-1 sequences collected between 9 March 2017 and 14 October 2021. Employing advanced statistical and bioinformatic techniques, the research delved deeply into understanding the transmission dynamics of the HIV-1 epidemic in Cyprus, as well as the monitoring of HIV-1's genetic diversity and the surveillance of transmitted drug resistance. The characterization of Cyprus's HIV-1 epidemic revealed a diverse landscape, comprising 21 HIV-1 group M pure subtypes and circulating recombinant forms (CRFs), alongside numerous uncharacterized recombinant strains. Subtypes A1 and B emerged as the most prevalent strains, followed by CRF02_AG. The findings of this study also revealed high levels of transmitted drug resistance (TDR) patterns, raising concerns for the efficacy of cART. The demographic profiles of individuals involved in HIV-1 transmission underscored the disproportionate burden borne by young to middle-aged Cypriot males, particularly those in the MSM community, who reported contracting the virus in Cyprus. An assessment of the spatiotemporal evolutionary dynamics illustrated the global interconnectedness of HIV-1 transmission networks, implicating five continents in the dissemination of strains within Cyprus: Europe, Africa, Asia, North America, and Oceania. Overall, this study advances the comprehension of the HIV-1 epidemic in Cyprus and highlights the importance of understanding HIV-1's transmission dynamics through continuous surveillance efforts. Furthermore, this work emphasizes the critical role of state-of-the-art bioinformatics analyses in addressing the challenges posed by HIV-1 transmission globally, laying the groundwork for public health interventions aimed at curbing its spread and improving patient outcomes.

Phylogenetic evidence of extensive spatial mixing of diverse HIV-1 group M lineages within Cameroon but not between its neighbours

From the perspective of developing relevant interventions for treating HIV and controlling its spread, it is particularly important to comprehensively understand the underlying diversity of the virus, especially in countries where the virus has been present and evolving since the cross-species transmission event that triggered the global pandemic. Here, we generate and phylogenetically analyse sequences derived from the gag-protease (2010 bp; n = 115), partial integrase (345 bp; n = 36), and nef (719 bp; n = 321) genes of HIV-1 group M (HIV-1M) isolates sampled between 2000 and 2022 from two cosmopolitan cities and 40 remote villages of Cameroon. While 52.4% of all sequenced viruses belonged to circulating recombinant form (CRF) 02_AG (CRF02_AG), the remainder were highly diverse, collectively representing seven subtypes and sub-subtypes, eight CRFs, and 36 highly divergent lineages that fall outside the established HIV-1M classification. Additionally, in 77 samples for which at least two genes were typed, 31% of the studied viruses apparently had fragments from viruses belonging to different clades. Furthermore, we found that the distribution of HIV-1M populations is similar between different regions of Cameroon. In contrast, HIV-1M demographics in Cameroon differ significantly from those in its neighbouring countries in the Congo Basin (CB). In phylogenetic trees, viral sequences cluster according to the countries where they were sampled, suggesting that while there are minimal geographical or social barriers to viral dissemination throughout Cameroon, there is strongly impeded dispersal of HIV-1M lineages between Cameroon and other locations of the CB. This suggests that the apparent stability of highly diverse Cameroonian HIV-1M populations may be attributable to the extensive mixing of human populations within the country and the concomitant trans-national movements of major lineages with very similar degrees of fitness; coupled with the relatively infrequent inter-national transmission of these lineages from neighbouring countries in the CB.

Differences in neutralizing antibody sensitivities and envelope characteristics indicate distinct antigenic properties of Nigerian HIV-1 subtype G and CRF02_AG

The magnitude of the HIV-1 epidemic in Nigeria is second only to the subtype C epidemic in South Africa, yet the subtypes prevalent in Nigeria require further characterization. A panel of 50 subtype G and 18 CRF02_AG Nigerian HIV-1 pseudoviruses (PSV) was developed and envelope coreceptor usage, neutralization sensitivity and cross-clade reactivity were characterized. These PSV were neutralized by some antibodies targeting major neutralizing determinants, but potentially important differences were observed in specific sensitivities (eg. to sCD4, MPER and V2/V3 monoclonal antibodies), as well as in properties such as variable loop lengths, number of potential N-linked glycans and charge, demonstrating distinct antigenic characteristics of CRF02_AG and subtype G. There was preferential neutralization of the matched CRF/subtype when PSV from subtype G or CRF02_AG were tested using pooled plasma. These novel Nigerian PSV will be useful to study HIV-1 CRF- or subtype-specific humoral immune responses for subtype G and CRF02_AG.

Criteria for classification, nomenclature, and reference sequence selection for HIV sub-subtypes of CRF01_AE and CRF07_BC strains in China

The available knowledge regarding classification, nomenclature, and reference sequence selection for the various sub-subtypes of circulating recombinant forms (CRFs) is inadequate to fulfill the growing demands of research focused on HIV prevention. We analyzed the spread of CRF01_AE and CRF07_BC strains, mainly in China, to complement and update the existing nomenclature and to propose a reference sequence selection criteria for sub-subtypes of CRFs.

How to report and make sense of a new HIV-1 circulating recombinant form?

Co-circulation of multiple HIV-1 subtypes in the same high-risk groups leads to the on-going generation of various inter-subtype recombinants, including unique (URFs) and circulating (CRFs) recombinant forms, which brings a new challenge for the prevention and eradication of HIV/AIDS. Identification and prompt reporting of new CRFs will provide not only new insights into the understanding of genetic diversity and evolution of HIV-1, but also an early warning of potential prevalence of these variants. Currently, 140 HIV-1 CRFs have been described; however, their prevalence and clinical importance are less concerned. Apart from the mosaic genomic maps, less other valuable information, including the clinical and demographic data, genomic sequence characteristics, origin and evolutionary dynamics, as well as representative genomic fragments for determining the variants, are available for most of these CRFs. Accompanied with the growing increase of HIV-1 full-length genomic sequences, more and more CRFs will be identified in the near future due to the high recombination potential of HIV-1. Here, we discuss the prevalence and clinical importance of various HIV-1 CRFs and propose how to report and make sense of a new HIV-1 CRF.

Evaluation of the ABL NGS assay for HIV-1 drug resistance testing

HIV evolution and variability around the world requires special monitoring of the viral strains in infected people. High-throughput HIV sequencing and drug resistance testing techniques have become routinely available over the last few years. We conducted a study to assess the new CE-marked ABL NGS HIV genotyping assay on an Illumina® platform, to compare the results (the detection of resistance associated mutations (RAMs) detected in the three main targets: reverse transcriptase, protease, and integrase) with those produced by three Sanger-based assays, and to compare the assays' respective costs. For the 10 samples and a 20 % sensitivity threshold for the NGS technology, the percent agreement between the four assays ranged from 99.5 % to 100 %. We detected 4 more and 10 more RAMs of interest when we lowered the NGS assay's threshold to 10 % and 3 %, respectively. At a threshold of 3 %, the antiretroviral sensitivity interpretation algorithm (for protease inhibitors) was modified for only two patients. The NGS assay's unit cost fell rapidly as the number of samples per run increased. Compared with Sanger sequencing, the ABL NGS HIV genotyping assay is just as robust and somewhat more expensive but opens up interesting multiplexing perspectives for virology laboratories.

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.

Characterization of HIV-1 Epidemic in Kyrgyzstan

Kyrgyzstan has one of the highest rates of HIV-1 spread in Central Asia. In this study, we used molecular–epidemiological approaches to examine the HIV-1 epidemic in Kyrgyzstan. Samples were obtained from HIV-positive individuals who visited HIV/AIDS clinics. Partial pol gene sequences were used to identify HIV-1 subtypes and drug resistance mutations (DRMs) and to perform phylogenetic analysis. Genetic diversity and history reconstruction of the major HIV-1 subtypes were explored using BEAST. This study includes an analysis of 555 HIV-positive individuals. The study population was equally represented by men and women aged 1–72 years. Heterosexual transmission was the most frequent, followed by nosocomial infection. Men were more likely to acquire HIV-1 during injection drug use and while getting clinical services, while women were more likely to be infected through sexual contacts (p < 0.01). Heterosexual transmission was the more prevalent among individuals 25–49 years old; individuals over 49 years old were more likely to be persons who inject drugs (PWID). The major HIV-1 variants were CRF02_AG, CRF63_02A, and sub-subtype A6. Major DRMs were detected in 26.9% of the study individuals; 62.2% of those had DRMs to at least two antiretroviral (ARV) drug classes. Phylogenetic analysis revealed a well-defined structure of CRF02_AG, indicating locally evolving sub-epidemics. The lack of well-defined phylogenetic structure was observed for sub-subtype A6. The estimated origin date of CRF02_AG was January 1997; CRF63_02A, April 2004; and A6, June 1995. A rapid evolutionary dynamic of CRF02_AG and A6 among Kyrgyz population since the mid-1990s was observed. We observed the high levels of HIV-1 genetic diversity and drug resistance in the study population. Complex patterns of HIV-1 phylogenetics in Kyrgyzstan were found. This study highlights the importance of molecular–epidemiological analysis for HIV-1 surveillance and treatment implementation to reduce new HIV-1 infections.

The Role of Phylogenetics in Discerning HIV-1 Mixing among Vulnerable Populations and Geographic Regions in Sub-Saharan Africa: A Systematic Review

To reduce global HIV-1 incidence, there is a need to understand and disentangle HIV-1 transmission dynamics and to determine the geographic areas and populations that act as hubs or drivers of HIV-1 spread. In Sub-Saharan Africa (sSA), the region with the highest HIV-1 burden, information about such transmission dynamics is sparse. Phylogenetic inference is a powerful method for the study of HIV-1 transmission networks and source attribution. In this review, we assessed available phylogenetic data on mixing between HIV-1 hotspots (geographic areas and populations with high HIV-1 incidence and prevalence) and areas or populations with lower HIV-1 burden in sSA. We searched PubMed and identified and reviewed 64 studies on HIV-1 transmission dynamics within and between risk groups and geographic locations in sSA (published 1995-2021). We describe HIV-1 transmission from both a geographic and a risk group perspective in sSA. Finally, we discuss the challenges facing phylogenetic inference in mixed epidemics in sSA and offer our perspectives and potential solutions to the identified challenges.

Elucidation of Early Evolution of HIV-1 Group M in the Congo Basin Using Computational Methods

The Congo Basin region is believed to be the site of the cross-species transmission event that yielded HIV-1 group M (HIV-1M). It is thus likely that the virus has been present and evolving in the region since that cross-species transmission. As HIV-1M was only discovered in the early 1980s, our directly observed record of the epidemic is largely limited to the past four decades. Nevertheless, by exploiting the genetic relatedness of contemporary HIV-1M sequences, phylogenetic methods provide a powerful framework for investigating simultaneously the evolutionary and epidemiologic history of the virus. Such an approach has been taken to find that the currently classified HIV-1 M subtypes and Circulating Recombinant Forms (CRFs) do not give a complete view of HIV-1 diversity. In addition, the currently identified major HIV-1M subtypes were likely genetically predisposed to becoming a major component of the present epidemic, even before the events that resulted in the global epidemic. Further efforts have identified statistically significant hot- and cold-spots of HIV-1M subtypes sequence inheritance in genomic regions of recombinant forms. In this review we provide ours and others recent findings on the emergence and spread of HIV-1M variants in the region, which have provided insights into the early evolution of this virus.

Forensic genomics of a novel Klebsiella quasipneumoniae type from a neonatal intensive care unit in China reveals patterns of colonization, evolution and epidemiology

During March 2017, a neonatal patient with severe diarrhoea subsequently developed septicaemia and died, with Klebsiella isolated as the causative microorganism. In keeping with infection control protocols, the coincident illness of an attending staff member and three other neonates with Klebsiella infection triggered an outbreak response, leading to microbiological assessment of isolates collected from the staff member and all 21 co-housed neonates. Multilocus sequence typing and genomic sequencing identified that the isolates from the 21 neonates were of a new Klebsiella sequence type, ST2727, and taxonomically belonged to K. quasipneumoniae subsp. similipneumoniae (formerly referred to as KpIIB). Genomic characterization showed that the isolated ST2727 strains had diverged from other K. quasipneumoniae subsp. similipneumoniae strains at least 90 years ago, whereas the neonatal samples were highly similar with a genomic divergence of 3.6 months. There was no relationship to the Klebsiella isolate from the staff member. This demonstrates that no transmission occurred from staff to patient or between patients. Rather, the data suggest that ST2727 colonized each neonate from a common hospital source. Sequence-based analysis of the genomes revealed several genes for antimicrobial resistance and some virulence features, but suggest that ST2727 is neither extremely-drug resistant nor hypervirulent. Our results highlight the clinical significance and genomic properties of ST2727 and urge genome-based measures be implemented for diagnostics and surveillance within hospital environments. Additionally, the present study demonstrates the need to scale the power of genomic analysis in retrospective studies where relatively few samples are available.

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.

High HIV burden and recent transmission chains in rural forest areas in southern Cameroon, where ancestors of HIV-1 have been identified in ape populations

We studied HIV prevalence and genetic diversity in rural forest areas in Cameroon, where chimpanzee and gorilla populations infected with the ancestors of the different HIV-1 groups have been identified and transmitted to humans during the 20th century. A total of 2812 individuals were studied, 924 from south-central, 1116 from south-east and 772 from south-west Cameroon. Of 208 (7.4%) samples that were confirmed for HIV-1 infection all belong to HIV-1 group M. In all sites and in all age categories, HIV-1 prevalence was higher in women (160/1599 (10.0%)) as compared to men (48/1213 (4.0%)) with the highest prevalence in women aged between 25 and 34 years (>17%). For 188/208 (92.3%) HIV-1 positive individuals, a fragment of the pol gene was successfully amplified and sequenced. Phylogenetic analysis showed predominance of CRF02_AG (58%), a large diversity of subtypes (A, D, F2 and G), nine different CRFs and more than 12% URFs. Interestingly, 35/188 (18.6%) HIV-1 strains form 12 recent transmission chains. The majority of the clusters are composed of two (n = 8) or three (n = 3) sequences but one cluster included ten HIV-1 strains from women living in four different villages on a major road for logging concessions in the south-east (60 km distance). In the three regions of Cameroon where the ancestors of the four HIV-1 groups have been transmitted to humans, we observed a high HIV prevalence, especially in the southeast where HIV-1 M originated. Many factors allowing rapid establishment in the human population and subsequent rapid spread to urban areas of a new retrovirus or other pathogens of zoonotic origin are now present. Our study shows clearly that some rural areas should also be considered as hot-spots for HIV infection. Prevention efforts together with growing access to HIV diagnosis and antiretroviral treatment are urgently needed in these remote areas.

Characterisation of HIV-1 Molecular Epidemiology in Nigeria: Origin, Diversity, Demography and Geographic Spread

Nigeria has the highest number of AIDS-related deaths in the world. In this study, we characterised the HIV-1 molecular epidemiology by analysing 1442 HIV-1 pol sequences collected 1999-2014 from four geopolitical zones in Nigeria using state-of-the-art maximum-likelihood and Bayesian phylogenetic analyses. The main circulating forms were the circulating recombinant form (CRF) 02_AG (44% of the analysed sequences), CRF43_02G (16%), and subtype G (8%). Twenty-three percent of the sequences represented unique recombinant forms (URFs), whereof 37 (11%) could be grouped into seven potentially novel CRFs. Bayesian phylodynamic analysis suggested that five major Nigerian HIV-1 sub-epidemics were introduced in the 1960s and 1970s, close to the Nigerian Civil War. The analysis also indicated that the number of effective infections decreased in Nigeria after the introduction of free antiretroviral treatment in 2006. Finally, Bayesian phylogeographic analysis suggested gravity-like dynamics in which virus lineages first emerge and expand within large urban centers such as Abuja and Lagos, before migrating towards smaller rural areas. This study provides novel insight into the Nigerian HIV-1 epidemic and may have implications for future HIV-1 prevention strategies in Nigeria and other severely affected countries.

HIV-1 drug resistance testing is essential for heavily-treated patients switching from first- to second-line regimens in resource-limited settings: evidence from routine clinical practice in Cameroon

Background

With the phase-out of stavudine (d4T), change to first-line regimens with zidovudine (AZT) or tenofovir (TDF) in resource-limited settings (RLS) might increase risks of cross-resistance to nucleos(t) ide reverse transcriptase inhibitors (NRTI). This would restrict the scope of switching to the World Health Organisation (WHO)-recommended standard second-line combinations (SLC) without HIV drug resistance (HIVDR)-testing in routine clinical practice.

Methods

An observational study was conducted among 101 Cameroonian patients (55.4% male, median [IQR] age 34 [10–41] years) failing first-line antiretroviral therapy (ART) in 2016, and stratified into three groups according to NRTIs exposure: exposure to both thymidine analogues AZT “and” D4T (group-A, n = 55); exposure to both TDF and AZT “or” D4T (group-B, n = 22); exposure solely to D4T (group-C, n = 24). Protease-reverse transcriptase HIVDR was interpreted using the HIVdb penalty scores (≥60: high-resistance; 20–59: intermediate-resistance; < 20: susceptible). The acceptable threshold for potential-efficacy was set at 80%.

Results

The median [IQR] CD4, viral RNA, and time on ART, were respectively 129 [29–466] cells/μl, 71,630 [19,041-368,000] copies/ml, and 4 [2–5] years. Overall HIVDR-level was 89.11% (90/101), with 83.2% harbouring M184 V (high-level 3TC/FTC-resistance) and only 1.98% (2/101) major HIVDR-mutations to ritonavir-boosted protease-inhibitors (PI/r). Thymidine-analogue mutations (TAMs)-1 [T215FY (46.53%), M41 L (22.77%), L210 W (8.91%)], with cross-resistance to AZT and TDF, were higher compared to TAMs-2 [D67N (21.78%), K70R (19.80%), K219QE (18.81%)]. As expected, K65R was related with TDF-exposure: 0% (0/55) in group-A, 22.72% (5/22) group-B, 4.17% (1/24) group-C (p = 0.0013). The potential-efficacy of AZT vs. TDF was respectively 43.64% (24/55) vs. 70.91% (39/55) in group-A (p = 0.0038); 63.64% (14/22) vs. 68.28% (15/22) in group-B (p = 1.0000); and 37.50% (9/24) vs. 83.33% (20/24) in group-C (p = 0.0032). CRF02_AG was the prevailing subtype (63.40%), followed by CRF11.cpx (8.91%), A₁ (7.92%), G (5.94%); without any significant effect of the subtype-distribution on HIVDR (92.2% in CRF02_AG vs. 83.8% in non-AG; p = 0.204).

Conclusion

First-line ART-failure exhibits high-level NRTI-resistance, with potential lower-efficacy of AZT compared to TDF. Significantly, using our 80% efficacy-threshold, only patients without NRTI-substitution on first-line could effectively switch to SLC following the WHO-approach. Patients with multiple NRTI-substitutions (exposed to both thymidine-analogues and TDF) on first-line ART would require HIVDR-testing to select active NRTIs for SLC.

Electronic supplementary material

The online version of this article (10.1186/s12879-019-3871-0) contains supplementary material, which is available to authorized users.

Spatiotemporal Characteristics of the Largest HIV-1 CRF02_AG Outbreak in Spain: Evidence for Onward Transmissions

Background and Aim: The circulating recombinant form 02_AG (CRF02_AG) is the predominant clade among the human immunodeficiency virus type-1 (HIV-1) non-Bs with a prevalence of 5.97% (95% Confidence Interval-CI: 5.41–6.57%) across Spain. Our aim was to estimate the levels of regional clustering for CRF02_AG and the spatiotemporal characteristics of the largest CRF02_AG subepidemic in Spain.

Methods: We studied 396 CRF02_AG sequences obtained from HIV-1 diagnosed patients during 2000–2014 from 10 autonomous communities of Spain. Phylogenetic analysis was performed on the 391 CRF02_AG sequences along with all globally sampled CRF02_AG sequences (N = 3,302) as references. Phylodynamic and phylogeographic analysis was performed to the largest CRF02_AG monophyletic cluster by a Bayesian method in BEAST v1.8.0 and by reconstructing ancestral states using the criterion of parsimony in Mesquite v3.4, respectively.

Results: The HIV-1 CRF02_AG prevalence differed across Spanish autonomous communities we sampled from (p < 0.001). Phylogenetic analysis revealed that 52.7% of the CRF02_AG sequences formed 56 monophyletic clusters, with a range of 2–79 sequences. The CRF02_AG regional dispersal differed across Spain (p = 0.003), as suggested by monophyletic clustering. For the largest monophyletic cluster (subepidemic) (N = 79), 49.4% of the clustered sequences originated from Madrid, while most sequences (51.9%) had been obtained from men having sex with men (MSM). Molecular clock analysis suggested that the origin (t_MRCA) of the CRF02_AG subepidemic was in 2002 (median estimate; 95% Highest Posterior Density-HPD interval: 1999–2004). Additionally, we found significant clustering within the CRF02_AG subepidemic according to the ethnic origin.

Conclusion: CRF02_AG has been introduced as a result of multiple introductions in Spain, following regional dispersal in several cases. We showed that CRF02_AG transmissions were mostly due to regional dispersal in Spain. The hot-spot for the largest CRF02_AG regional subepidemic in Spain was in Madrid associated with MSM transmission risk group. The existence of subepidemics suggest that several spillovers occurred from Madrid to other areas. CRF02_AG sequences from Hispanics were clustered in a separate subclade suggesting no linkage between the local and Hispanic subepidemics.

Combining Phylogenetic and Network Approaches to Identify HIV-1 Transmission Links in San Mateo County, California

The HIV epidemic in San Mateo County is sustained by multiple overlapping risk groups and is an important hub for HIV transmission in northern California. Limited access to care has led historically to delayed clinical presentation, higher rates of opportunistic infections, and an increased prevalence of antiretroviral drug resistance. The virologic and clinical consequences of treatment within these multiple ethnic and behavioral groups are poorly understood, highlighting the need for efficient surveillance strategies that are able to elucidate transmission networks and drug resistance patterns. We obtained sequence data from a group of 316 HIV-positive individuals in the San Mateo AIDS Program over a 14-year period and integrated epidemiologic, phylogenetic, and network approaches to characterize transmission clusters, risk factors and drug resistance. Drug resistance mutations were identified using the Stanford HIV Drug Resistance Database. A maximum likelihood tree was inferred in RAxML and subjected to clustering analysis in Cluster Picker. Network analysis using pairwise genetic distances was performed in HIV-TRACE. Participants were primarily male (60%), white Hispanics and non-Hispanics (32%) and African American (20.6%). The most frequent behavior risk factor was male-male sex (33.5%), followed by heterosexual (23.4%) and injection drug use (9.5%). Nearly all sequences were subtype B (96%) with subtypes A, C, and CRF01_AE also observed. Sequences from 65% of participants had at least one drug resistance mutation. Clustered transmissions included a higher number of women when compared to non-clustered individuals and were more likely to include heterosexual or people who inject drugs (PWID). Detailed analysis of the largest network (N = 47) suggested that PWID played a central role in overall transmission of HIV-1 as well as bridging men who have sex with men (MSM) transmission with heterosexual/PWID among primarily African American men. Combined phylogenetic and network analysis of HIV sequence data identified several overlapping risk factors in the epidemic, including MSM, heterosexual and PWID transmission with a disproportionate impact on African Americans and a high prevalence of drug resistance.

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

Eastern European countries, including Russia, Ukraine, and other former Soviet Union (FSU) countries, have experienced a human immunodeficiency virus (HIV) epidemic spreading mostly among people who inject drugs (PWID). We aimed to investigate the origin and the dispersal patterns of HIV-1 CRF02_AG in Russia and other FSU countries. We studied 136 CRF02_AG sequences originating from Russia, and FSU countries along with a globally sampled dataset of 3,580 CRF02_AG sequences. Maximum-likelihood phylogeny reconstruction with bootstrap evaluation was conducted in RAxML. Bayesian phylogeographic analysis was performed in BEAST v1.8 using the discrete trait model. We found that all CRF02_AG sequences from Russia and other FSU countries formed a single monophyletic cluster within CRF02_AG radiation. The Russian/FSU clade was classified as CRF63_02A1. Sequences from the FSU countries clustered further within distinct subclades (two from Russia, three from Uzbekistan, and one Kazakhstan) according to the geographic origin of sampling. Molecular clock analysis revealed that the time to the most recent common ancestor (tMRCA) of the CRF63_02A1 epidemic was in 1996 [95% higher posterior density (95% HPD): 1992-1999], while for the two Russian subclades, tMRCA was estimated in 2003 (95% HPD: 2001-2004) and in 2007 (95% HPD: 2005-2008). Phylogeographic analysis suggested that the potential origin of the epidemic was in Uzbekistan. Early dispersal of CRF63_02A1 occurred in Kazakhstan and Uzbekistan and thereafter the epidemic spread to Russia. Notably, spillover transmissions to Russia kept occurring from both countries. Previous studies have shown that Russia and Ukraine have provided the source for the PWID-driven, HIV-1 subtype-A epidemic, spreading across the FSU countries (AFSU). In great contrast, CRF63_02A1 established an epidemic in central Asia (Uzbekistan and Kazakhstan), from where it subsequently disseminated to Russia. Our study suggests that cross-border transmissions among PWID occur bidirectionally between Russian and other FSU populations. These results are of public health importance and suggest that prevention actions have to be reinforced in this area to assist the management of high-risk practices.

Divergent HIV-1 strains (CRF92_C2U and CRF93_cpx) co-circulating in the Democratic Republic of the Congo: Phylogenetic insights on the early evolutionary history of subtype C

Molecular epidemiological studies revealed that the epicenter of the HIV pandemic was Kinshasa, the capital city of the Democratic Republic of the Congo (DRC) in Central Africa. All known subtypes and numerous complex recombinant strains co-circulate in the DRC. Moreover, high intra-subtype diversity has been also documented. During two previous surveys on HIV-1 antiretroviral drug resistance in the DRC, we identified two divergent subtype C lineages in the protease and partial reverse transcriptase gene regions. We sequenced eight near full-length genomes and classified them using bootscanning and likelihood-based phylogenetic analyses. Four strains are more closely related to subtype C although within the range of inter sub-subtype distances. However, these strains also have small unclassified fragments and thus were named CRF92_C2U. Another strain is a unique recombinant of CRF92_C2U with an additional small unclassified fragment and a small divergent subtype A fragment. The three remaining strains represent a complex mosaic named CRF93_cpx. CRF93_cpx have two fragments of divergent subtype C sequences, which are not conventional subtype C nor the above described C2, and multiple divergent subtype A-like fragments. We then inferred the time-scaled evolutionary history of subtype C following a Bayesian approach and a partitioned analysis using major genomic regions. CRF92_C2U and CRF93_cpx had the most recent common ancestor with conventional subtype C around 1932 and 1928, respectively. A Bayesian demographic reconstruction corroborated that the subtype C transition to a faster phase of exponential growth occurred during the 1950s. Our analysis showed considerable differences between the newly discovered early-divergent strains and the conventional subtype C and therefore suggested that this virus has been diverging in humans for several decades before the HIV/M diversity boom in the 1950s.

Assessing the danger of self-sustained HIV epidemics in heterosexuals by population based phylogenetic cluster analysis

Assessing the danger of transition of HIV transmission from a concentrated to a generalized epidemic is of major importance for public health. In this study, we develop a phylogeny-based statistical approach to address this question. As a case study, we use this to investigate the trends and determinants of HIV transmission among Swiss heterosexuals. We extract the corresponding transmission clusters from a phylogenetic tree. To capture the incomplete sampling, the delayed introduction of imported infections to Switzerland, and potential factors associated with basic reproductive number R0, we extend the branching process model to infer transmission parameters. Overall, the R0 is estimated to be 0.44 (95%-confidence interval 0.42—0.46) and it is decreasing by 11% per 10 years (4%—17%). Our findings indicate rather diminishing HIV transmission among Swiss heterosexuals far below the epidemic threshold. Generally, our approach allows to assess the danger of self-sustained epidemics from any viral sequence data.

In epidemiology, the “basic reproductive number” describes how efficiently a disease is transmitted, and represents the average number of new infections that an infected individual causes. If this number is less than one, many people do not infect anybody and hence the transmission chains die out. On the other hand, if the basic reproductive number is larger than one, an infected person infects on average more than one new individual, which leads to the virus or bacteria spreading in a self-sustained way.

Turk et al. have now developed a method to estimate the basic reproductive number using the genetic sequences of the virus or bacteria, and have used it to investigate how efficiently HIV spreads among Swiss heterosexuals. The results show that the basic reproductive number of HIV in this group is far below the critical value of one and that over the last years this number has been decreasing. Furthermore, the basic reproductive number differs for different subtypes of the HIV virus, indicating that the geographical region where the infection was acquired may play a role in transmission. Turk et al. also found that people who are diagnosed later or who often have sex with occasional partners spread the virus more efficiently.

These findings might be helpful for policy makers as they indicate that the risk of self-sustained transmission in this group in Switzerland is small. Furthermore the method allows HIV epidemics to be monitored at high resolution using sequence data, assesses the success of currently implemented preventive measures, and helps to target subgroups who are at higher risk of an infection – for instance, by supporting frequent HIV testing of these people. The method developed by Turk et al. could also prove useful for assessing the danger of other epidemics.

Biased phylodynamic inferences from analysing clusters of viral sequences

Phylogenetic methods are being increasingly used to help understand the transmission dynamics of measurably evolving viruses, including HIV. Clusters of highly similar sequences are often observed, which appear to follow a ‘power law’ behaviour, with a small number of very large clusters. These clusters may help to identify subpopulations in an epidemic, and inform where intervention strategies should be implemented. However, clustering of samples does not necessarily imply the presence of a subpopulation with high transmission rates, as groups of closely related viruses can also occur due to non-epidemiological effects such as over-sampling. It is important to ensure that observed phylogenetic clustering reflects true heterogeneity in the transmitting population, and is not being driven by non-epidemiological effects. We qualify the effect of using a falsely identified ‘transmission cluster’ of sequences to estimate phylodynamic parameters including the effective population size and exponential growth rate under several demographic scenarios. Our simulation studies show that taking the maximum size cluster to re-estimate parameters from trees simulated under a randomly mixing, constant population size coalescent process systematically underestimates the overall effective population size. In addition, the transmission cluster wrongly resembles an exponential or logistic growth model 99% of the time. We also illustrate the consequences of false clusters in exponentially growing coalescent and birth-death trees, where again, the growth rate is skewed upwards. This has clear implications for identifying clusters in large viral databases, where a false cluster could result in wasted intervention resources.

Time-scale of minor HIV-1 complex circulating recombinant forms from Central and West Africa

Background

Several HIV-1 circulating recombinant forms with a complex mosaic structure (CRFs_cpx) circulate in central and western African regions. Here we reconstruct the evolutionary history of some of these complex CRFs (09_cpx, 11_cpx, 13_cpx and 45_cpx) and further investigate the dissemination dynamic of the CRF11_cpx clade by using a Bayesian coalescent-based method.

Results

The analysis of two HIV-1 datasets comprising 181 pol (36 CRF09_cpx, 116 CRF11_cpx, 20 CRF13_cpx and 9 CRF45_cpx) and 125 env (12 CRF09_cpx, 67 CRF11_cpx, 17 CRF13_cpx and 29 CRF45_cpx) sequences pointed to quite consistent onset dates for CRF09_cpx (~1966: 1958–1979), CRF11_cpx (~1957: 1950–1966) and CRF13_cpx (~1965: 1958–1973) clades; while some divergence was found for the estimated date of origin of CRF45_cpx clade [pol = 1970 (1964–1976); env = 1960 (1952–1969)]. Phylogeographic reconstructions indicate that the HIV-1 CRF11_cpx clade most probably emerged in Cameroon and from there it was first disseminated to the Central Africa Republic and Chad in the early 1970s and to other central and western African countries from the early 1980s onwards. Demographic reconstructions suggest that the CRF11_cpx epidemic grew between 1960 and 1990 with a median exponential growth rate of 0.27 year⁻¹, and stabilized after.

Conclusions

These results reveal that HIV-1 CRFs_cpx clades have been circulating in Central Africa for a period comparable to other much more prevalent HIV-1 group M lineages. Cameroon was probably the epicenter of dissemination of the CRF11_cpx clade that seems to have experienced a long epidemic growth phase before stabilization. The epidemic growth of the CRF11_cpx clade was roughly comparable to other HIV-1 group M lineages circulating in Central Africa.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0824-8) contains supplementary material, which is available to authorized users.

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

Kemerovo Oblast (KO) has had the highest rate of HIV spread in Russia since 2011. The aim of this work was to study the genetic variation of HIV-1 in Kemerovo Oblast. Blood was sampled from a total of 91 HIV-positive antiretroviral-therapy-naïve individuals in 2013 (38) and 2015 (53). HIV-1 subtypes, pol gene drug resistance mutations, and viral tropism were analyzed. In 2013-2015, the prevalence of HIV-1 subtype A decreased in KO from 60.5 to 7.5 %. The samples collected in 2015 from the patients with newly diagnosed HIV demonstrate the current dominance of HIV-1 CRF63_02A1 (71.7 %) and HIV-1 URF63_A1 (20.8 %), their parental viruses being CRF63_02A1 and subtype A. The initially predominant genetic variant, HIV-1 subtype A, was replaced in KO. An unusually high incidence of HIV-1 unique recombinant forms is probably the result of HIV-1 CRF63_02A1 introduction in the group of injection drug users with the initial HIV-1 subtype A infection and the practice of risky behavior that promotes reinfection. HIV-1 CRF63_02A1, which recently emerged in Siberia, and its recombinant forms have an ever-increasing impact on the current HIV epidemic in Russia, making urgent the need for in-depth study of this HIV-1 genetic variant.