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

The Viral Founder Effect and Economic-Driven Human Mobility Shaped the Distinct Epidemic Pattern of HIV-1 CRF01_AE in Northeast China

Background: In China, two distinct lineages shaped the epidemic of HIV-1 CRF01_AE among men who have sex with men (MSM), of which the uneven distributions were observed geographically. One lineage spread across China, while another dominated in Northeast China. Understanding the drivers of viral diffusion would provide guidelines for identifying the source and hotspots of HIV transmission among MSM to target interventions in China.

Methods: We collected the pol sequences between 2002–2017 to reconstruct the spatiotemporal history of CRF01_AE lineages in Shenyang, one economic center of Northeast China, using the Bayesian phylogeographic and phylodynamic approaches. Importantly, for the datasets with the high sample density, we did the down-sampling to avoid the sampling bias.

Results: Two lineages accounted for 97%, including 426 and 1516 sequences, and homosexuals and bisexuals were above 80%. One lineage appeared earlier 7 years than another (1993 vs. 2002) among homosexuals and bisexuals, whereas among heterosexuals, both lineages were observed firstly in 2002. 96% viral migrations within one lineage were from homosexuals toward bisexuals (49%) and male-heterosexuals (46%). Within another, except for homosexuals (72%), bisexuals (23%) served as the top second source, and female-heterosexuals (11%) were the third recipients following bisexuals (44%) and male-heterosexuals (39%). Although the basic reproduction number (R₀) of two lineages were similar and both of the effective production number (R_e) fell below 1 at the most recent sampling time, the starts of the R_e declining varied.

Conclusions: Our findings revealed that throughout the viral national spread chain, Shenyang is the source for the initial expanding of one lineage, where is only a sink of another, proving that the viral founder effect and regional human mobility contributed to the uneven distribution of two lineages, and emphasizing the important roles of the area where the virus originated and economy-driven migrants in HIV transmission.

Identification of CRF66_BF, a New HIV-1 Circulating Recombinant Form of South American Origin

Circulating recombinant forms (CRFs) are important components of the HIV-1 pandemic. Among 110 reported in the literature, 17 are BF1 intersubtype recombinant, most of which are of South American origin. Among these, all 5 identified in the Southern Cone and neighboring countries, except Brazil, derive from a common recombinant ancestor related to CRF12_BF, which circulates widely in Argentina, as deduced from coincident breakpoints and clustering in phylogenetic trees. In a HIV-1 molecular epidemiological study in Spain, we identified a phylogenetic cluster of 20 samples from 3 separate regions which were of F1 subsubtype, related to the Brazilian strain, in protease-reverse transcriptase (Pr-RT) and of subtype B in integrase. Remarkably, 14 individuals from this cluster (designated BF9) were Paraguayans and only 4 were native Spaniards. HIV-1 transmission was predominantly heterosexual, except for a subcluster of 6 individuals, 5 of which were men who have sex with men. Ten additional database sequences, from Argentina (n = 4), Spain (n = 3), Paraguay (n = 1), Brazil (n = 1), and Italy (n = 1), branched within the BF9 cluster. To determine whether it represents a new CRF, near full-length genome (NFLG) sequences were obtained for 6 viruses from 3 Spanish regions. Bootscan analyses showed a coincident BF1 recombinant structure, with 5 breakpoints, located in p17 gag , integrase, gp120, gp41-rev overlap, and nef, which was identical to that of two BF1 recombinant viruses from Paraguay previously sequenced in NFLGs. Interestingly, none of the breakpoints coincided with those of CRF12_BF. In a maximum likelihood phylogenetic tree, all 8 NFLG sequences grouped in a strongly supported clade segregating from previously identified CRFs and from the CRF12_BF "family" clade. These results allow us to identify a new HIV-1 CRF, designated CRF66_BF. Through a Bayesian coalescent analysis, the most recent common ancestor of CRF66_BF was estimated around 1984 in South America, either in Paraguay or Argentina. Among Pr-RT sequences obtained by us from HIV-1-infected Paraguayans living in Spain, 14 (20.9%) of 67 were of CRF66_BF, suggesting that CRF66_BF may be one of the major HIV-1 genetic forms circulating in Paraguay. CRF66_BF is the first reported non-Brazilian South American HIV-1 CRF_BF unrelated to CRF12_BF.

Increased Frequency of Inter-Subtype HIV-1 Recombinants Identified by Near Full-Length Virus Sequencing in Rwandan Acute Transmission Cohorts

Most studies of HIV-1 transmission have focused on subtypes B and C. In this study, we determined the genomic sequences of the transmitted founder (TF) viruses from acutely infected individuals enrolled between 2005 and 2011 into IAVI protocol C in Rwanda and have compared these isolates to viruses from more recent (2016-2019) acute/early infections in three at risk populations - MSM, high risk women (HRW), and discordant couples (DC). For the Protocol C samples, we utilized near full-length single genome (NFLG) amplification to generate 288 HIV-1 amplicons from 26 acutely infected seroconverters (SC), while for the 21 recent seroconverter samples (13 from HRW, two from DC, and six from MSM), we PCR amplified overlapping half-genomes. Using PacBio SMRT technology combined with the MDPseq workflow, we performed multiplex sequencing to obtain high accuracy sequences for each amplicon. Phylogenetic analyses indicated that the majority of recent transmitted viruses from DC and HRW clustered within those of the earlier Protocol C cohort. However, five of six sequences from the MSM cohort branched together and were greater than 97% identical. Recombination analyses revealed a high frequency (6/26; 23%) of unique inter-subtype recombination in Protocol C with 19% AC and 4% CD recombinant viruses, which contrasted with only 6.5% of recombinants defined by sequencing of the pol gene previously. The frequency of recombinants was significantly higher (12/21; 57%) in the more recent isolates, although, the five related viruses from the MSM cohort had identical recombination break points. While major drug resistance mutations were absent from Protocol C viruses, 4/21 of recent isolates exhibited transmitted nevirapine resistance. These results demonstrate the ongoing evolution and increased prevalence of recombinant and drug resistant transmitted viruses in Rwanda and highlight the importance of defining NFLG sequences to fully understand the nature of TF viruses and in particular the prevalence of unique recombinant forms (URFs) in transmission cohorts.

Large Evolutionary Rate Heterogeneity among and within HIV-1 Subtypes and CRFs

HIV-1 is a fast-evolving, genetically diverse virus presently classified into several groups and subtypes. The virus evolves rapidly because of an error-prone polymerase, high rates of recombination, and selection in response to the host immune system and clinical management of the infection. The rate of evolution is also influenced by the rate of virus spread in a population and nature of the outbreak, among other factors. HIV-1 evolution is thus driven by a range of complex genetic, social, and epidemiological factors that complicates disease management and prevention. Here, we quantify the evolutionary (substitution) rate heterogeneity among major HIV-1 subtypes and recombinants by analyzing the largest collection of HIV-1 genetic data spanning the widest possible geographical (100 countries) and temporal (1981–2019) spread. We show that HIV-1 substitution rates vary substantially, sometimes by several folds, both across the virus genome and between major subtypes and recombinants, but also within a subtype. Across subtypes, rates ranged 3.5-fold from 1.34 × 10⁻³ to 4.72 × 10⁻³ in env and 2.3-fold from 0.95 × 10⁻³ to 2.18 × 10⁻³ substitutions site⁻¹ year⁻¹ in pol. Within the subtype, 3-fold rate variation was observed in env in different human populations. It is possible that HIV-1 lineages in different parts of the world are operating under different selection pressures leading to substantial rate heterogeneity within and between subtypes. We further highlight how such rate heterogeneity can complicate HIV-1 phylodynamic studies, specifically, inferences on epidemiological linkage of transmission clusters based on genetic distance or phylogenetic data, and can mislead estimates about the timing of HIV-1 lineages.

Global and Regional Estimates for Subtype-Specific Therapeutic and Prophylactic HIV-1 Vaccines: A Modeling Study

Global HIV-1 genetic diversity forms a major obstacle to the development of an HIV vaccine. It may be necessary to employ subtype-specific HIV-1 vaccines in individual countries according to their HIV-1 subtype distribution. We estimated the global and regional need for subtype-specific HIV-1 vaccines. We took into account the proportions of different HIV-1 variants circulating in each country, the genetic composition of HIV-1 recombinants, and the different genome segments (gag, pol, env) that may be incorporated into vaccines. We modeled different scenarios according to whether countries would employ subtype-specific HIV-1 vaccines against (1) the most common subtype; (2) subtypes contributing more than 5% of HIV infections; or (3) all circulating subtypes. For therapeutic vaccines targeting the most common HIV-1 subtype in each country, 16.5 million doses of subtype C vaccine were estimated globally, followed by subtypes A (14.3 million) and B (4.2 million). A vaccine based on env required 2.6 million subtype E doses, and a vaccine based on pol required 4.8 million subtype G doses. For prophylactic vaccines targeting the most common HIV-1 subtype in each country, 1.9 billion doses of subtype A vaccine were estimated globally, followed by subtype C (1.1 billion) and subtype B (1.0 billion). A vaccine based on env required 1.2 billion subtype E doses, and a vaccine based on pol required 0.3 billion subtype G doses. If subtype-specific HIV-1 vaccines are also directed against less common subtypes in each country, vaccines targeting subtypes D, F, H, and K are also needed and would require up to five times more vaccine doses in total. We conclude that to provide global coverage, subtype-specific HIV-1 vaccines need to be directed against subtypes A, B, and C. Vaccines targeting env also need to include subtype E and those targeting pol need to include subtype G.

Genetic Characterization of a New HIV-1 Sub-Subtype A in Cabo Verde, Denominated A8

Previous molecular characterization of Human immunodeficiency virus (HIV-1) samples from Cabo Verde pointed out a vast HIV-1 pol diversity, with several subtypes and recombinant forms, being 5.2% classified as AU-pol. Thus, the aim of the present study was to improve the characterization of these AU sequences. The genomic DNA of seven HIV-1 AU pol-infected individuals were submitted to four overlapping nested-PCR fragments aiming to compose the full-length HIV-1 genome. The final classification was based on phylogenetic trees that were generated using the maximum likelihood and bootscan analysis. The genetic distances were calculated using Mega 7.0 software. Complete genome amplification was possible for two samples, and partial genomes were obtained for the other five. These two samples grouped together with a high support value, in a separate branch from the other sub-subtypes A and CRF26_A5U. No recombination was verified at bootscan, leading to the classification of a new sub-subtype A. The intragroup genetic distance from the new sub-subtype A at a complete genome was 5.2%, and the intergroup genetic varied from 8.1% to 19.0% in the analyzed fragments. Our study describes a new HIV-1 sub-subtype A and highlights the importance of continued molecular surveillance studies, mainly in countries with high HIV molecular diversity.

Identification of CRF89_BF, a new member of an HIV-1 circulating BF intersubtype recombinant form family widely spread in South America

Circulating recombinant forms (CRFs) contribute substantially to the HIV-1 pandemic. Among 105 CRFs described in the literature, 16 are BF intersubtype recombinants, most of South American origin, of which CRF12_BF is the most widely spread. A BF recombinant cluster identified in Bolivia was suggested to represent a new CRF_BF. Here we find that it belongs to a larger cluster incorporating 39 viruses collected in 7 countries from 3 continents, 22 of them in Spain, most from Bolivian or Peruvian individuals, and 12 in South America (Bolivia, Argentina, and Peru). This BF cluster comprises three major subclusters, two associated with Bolivian and one with Peruvian individuals. Near full-length genome sequence analyses of nine viruses, collected in Spain, Bolivia, and Peru, revealed coincident BF mosaic structures, with 13 breakpoints, 6 and 7 of which coincided with CRF12_BF and CRF17_BF, respectively. In a phylogenetic tree, they grouped in a clade closely related to these CRFs, and more distantly to CRF38_BF and CRF44_BF, all circulating in South America. These results allowed to identify a new HIV-1 CRF, designated CRF89_BF. Through phylodynamic analyses, CRF89_BF emergence was estimated in Bolivia around 1986. CRF89_BF is the fifth CRF member of the HIV-1 recombinant family related to CRF12_BF.

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.