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

Identification of a novel first-generation HIV-1 circulating recombinant form (CRF152_DG) among people living with HIV in Karachi, Pakistan

The objective of this study was to characterize a novel circulating recombinant form of human immunodeficiency virus type 1 (HIV-1) among people living with HIV in Karachi, Pakistan. We conducted near-full-length genome (NFLG) sequencing on eight samples exhibiting D/G recombination signals in the pol gene region. We successfully obtained NFLG sequences (790-9,614; with reference to the HXB2 genome) from four of the eight samples and then conducted phylogenetic and recombination analyses on them. The four NFLG sequences from our study and one DG unique recombinant form previously identified in the United Kingdom (GenBank accession: MF109700) formed a distinct monophyletic cluster with an Shimodaira-Hasegawa approximate likelihood ratio test node support value of 100%. Bootscan analyses of the five NFLG sequences of DG recombinants showed that all five NFLGs shared the same unique mosaic pattern of recombination breakpoints between D and G clades, with two D fragments in the pol and vif regions inserted into a G backbone. Subregion phylogenetic analyses confirmed these sequences to be a novel circulating recombinant form (CRF) composed of subtypes D and G. The DG recombinant sequences were eventually designated as CRF152_DG by the Los Alamos HIV Sequence Database staff.

IMPORTANCE

In Pakistan, the genetic diversity of human immunodeficiency virus type 1 (HIV-1) is becoming increasingly complex, compared to the early years of the epidemic that started after the detection of the first cases of HIV-1 in 1987 in Karachi. Based on the available molecular studies, two dominant HIV-1 clades, sub-subtype A1 and CRF02_AG, have been found to co-circulate with other clades, namely B, C, D, G, CRF01_AE, CRF35_A1D, and CRF56_cpx, in various urban areas of Pakistan. Several novel recombinant forms have also been detected. This first report of CRF152_DG highlights the complex nature of the HIV epidemic in Pakistan and emphasizes the importance of continual molecular surveillance (ideally based on whole-genome sequences) of HIV.

SEEPS: A Simulation Tool for Understanding Impacts of Contact Tracing on Epidemiological Inference from Phylogenetic Data

Motivation

Robust sampling methods are foundational to inferences using phylogenies. Yet the impact of using contact tracing, a type of non-uniform sampling used in public health applications such as infectious disease outbreak investigations, is not well understood. To understand how this non-uniform sampling method influences a recovered phylogeny, a new simulation tool is needed.

Results

We developed a new simulation tool called SEEPS (Sequence Evolution and Epidemiological Process Simulator) that allows for the simulation of contact tracing and the resulting transmission tree, pathogen phylogeny, and corresponding virus genetic sequences. Importantly, SEEPS takes within-host evolution into account when generating pathogen phylogenies and sequences from transmission histories. Using SEEPS, we demonstrate that contact tracing can significantly impact the structure of the resulting tree, as described by popular tree statistics. We also examined real data from a 2007-2008 Swedish HIV-1 outbreak and the broader 1998-2010 European HIV-1 epidemic to highlight the differences in contact tracing and expected phylogenies. Aided by SEEPS, we show that the data collection of the Swedish outbreak was strongly influenced by contact tracing even after downsampling, while the broader European Union epidemic showed little evidence of universal contact tracing, agreeing with the known epidemiological information about sampling and spread. Overall, our results highlight the importance of including possible non-uniform sampling schemes when examining phylogenetic trees. For that, SEEPS serves as a useful tool to evaluate such impacts, thereby facilitating better phylogenetic inferences of the characteristics of a disease outbreak.

Availability

SEEPS is available at github.com/MolEvolEpid/SEEPS.

Human Immunodeficiency Virus Type-1 Genetic Diversity and Drugs Resistance Mutations among People Living with HIV in Karachi, Pakistan

The human immunodeficiency virus type-1 epidemic in Pakistan has significantly increased over the last two decades. In Karachi, Pakistan, there is a lack of updated information on the complexity of HIV-1 genetic diversity and the burden of drug resistance mutations (DRMs) that can contribute to ART failure and poor treatment outcomes. This study aimed to determine HIV-1 genetic diversity and identify drug-resistance mutations among people living with HIV in Karachi. A total of 364 HIV-positive individuals, with a median age of 36 years, were enrolled in the study. The HIV-1 partial pol gene was successfully sequenced from 268 individuals. The sequences were used to generate phylogenetic trees to determine clade diversity and also to assess the burden of DRMs. Based on the partial pol sequences, 13 distinct HIV-1 subtypes and recombinant forms were identified. Subtype A1 was the most common clade (40%), followed by CRF02_AG (33.2%). Acquired DRMs were found in 30.6% of the ART-experienced patients, of whom 70.7%, 20.7%, and 8.5% were associated with resistance to NNRTIs, NRTIs, and PIs, respectively. Transmitted DRMs were found in 5.6% of the ART-naïve patients, of whom 93% were associated with resistance against NNRTIs and 7% to PIs. The high prevalence of DRMs in ART-experienced patients poses significant challenges to the long-term benefits and sustainability of the ART program. This study emphasizes the importance of continuous HIV molecular epidemiology and drug resistance surveillance to support evidence-based HIV prevention, precise ART, and targeted AIDS care.

Long-term exposure to ambient particulate matter is associated with prognosis in people living with HIV/AIDS: Evidence from a longitudinal study

BACKGROUND

Evidence on the association between particulate matter (PM) exposure and prognosis in people living with HIV/AIDS (PWHA) is scarce. We aim to investigate the associations of long-term exposure to PM with AIDS-related deaths and complications.

METHODS

We collected follow-up information on 7444 PWHAs from 2000 to 2021 from the HIV/AIDS Comprehensive Response Information Management System of the Wuhan Center for Disease Control and Prevention. The AIDS-related deaths and complications were assessed by physicians every 3 to 6 months, and the monthly average PM concentrations for each PWHA were extracted from the China High Air Pollutants dataset. We employed time-varying Cox regression models to evaluate the associations of the average cumulative PM exposure concentrations with AIDS-related deaths and complications, as well as the mediating effects of AIDS-related complications in PM-induced AIDS-related deaths.

RESULTS

For each 1 μg/m3 increase in PM1, PM2.5, and PM10, the adjusted hazard ratios (HRs) for AIDS-related deaths were 1.021 (1.009, 1.033), 1.012 (1.005, 1.020), and 1.010 (1.005, 1.015), respectively; and the HRs for AIDS-related complications were 1.049 (1.034, 1.064), 1.029 (1.020, 1.038), and 1.031 (1.024, 1.037), respectively. AIDS-related complications mediated 18.38 % and 18.68 % of the association of exposure to PM1 and PM2.5 with AIDS-related deaths, respectively. The association of PM exposure with AIDS-related deaths was more significant in older PWHA. Meanwhile, the association between PM exposure and AIDS-related complications was stronger in PWHA with a BMI ≥ 24 kg/m2.

CONCLUSION

Long-term exposure to PM is positively associated with AIDS-related deaths and complications, and AIDS-related complications have mediating effects in PM-induced AIDS-related deaths. Our evidence emphasizes that enhanced protection against PM exposure for PWHAs is an additional mitigation strategy to reduce AIDS-related deaths and complications.

Optimization of genetic distance threshold for inferring the CRF01_AE molecular network based on next-generation sequencing

Introduction

HIV molecular network based on genetic distance (GD) has been extensively utilized. However, the GD threshold for the non-B subtype differs from that of subtype B. This study aimed to optimize the GD threshold for inferring the CRF01_AE molecular network.

Methods

Next-generation sequencing data of partial CRF01_AE pol sequences were obtained for 59 samples from 12 transmission pairs enrolled from a high-risk cohort during 2009 and 2014. The paired GD was calculated using the Tamura-Nei 93 model to infer a GD threshold range for HIV molecular networks.

Results

2,019 CRF01_AE pol sequences and information on recent HIV infection (RHI) from newly diagnosed individuals in Shenyang from 2016 to 2019 were collected to construct molecular networks to assess the ability of the inferred GD thresholds to predict recent transmission events. When HIV transmission occurs within a span of 1-4 years, the mean paired GD between the sequences of the donor and recipient within the same transmission pair were as follow: 0.008, 0.011, 0.013, and 0.023 substitutions/site. Using these four GD thresholds, it was found that 98.9%, 96.0%, 88.2%, and 40.4% of all randomly paired GD values from 12 transmission pairs were correctly identified as originating from the same transmission pairs. In the real world, as the GD threshold increased from 0.001 to 0.02 substitutions/site, the proportion of RHI within the molecular network gradually increased from 16.6% to 92.3%. Meanwhile, the proportion of links with RHI gradually decreased from 87.0% to 48.2%. The two curves intersected at a GD of 0.008 substitutions/site.

Discussion

A suitable range of GD thresholds, 0.008-0.013 substitutions/site, was identified to infer the CRF01_AE molecular transmission network and identify HIV transmission events that occurred within the past three years. This finding provides valuable data for selecting an appropriate GD thresholds in constructing molecular networks for non-B subtypes.

Molecular genetic characterization analysis of a novel HIV-1 circulating recombinant form (CRF156_0755) in Guangdong, China

Introduction

The characteristic of human immunodeficiency virus type 1 (HIV-1) is its susceptibility to erroneous replication and recombination, which plays a crucial role in the diverse and dynamic variation of HIV-1. The spread of different subtypes in the same population often leads to the emergence of circulating recombination forms (CRFs). At present, the main recombinant subtypes of HIV-1 in China are CRF07_BC, CRF01_AE, CRF08_BC and B' subtypes, while CRF55_01B has become the fifth major epidemic strain in China after rapid growth in recent years since it was first reported in 2013. In this study, we obtained five nearly full-length genomes (NFLGs) and one half-length genome from five different cities in Guangdong. Here, we focused on analyzing their characteristics, parental origin and drug resistance.

Methods

Plasma samples were collected from six HIV-1 infected patients in Guangdong Province who had no epidemiological association with each other. The NFLGs of HIV-1 were amplified in two overlapping segments by the near-terminal dilution method. The positive products were sequenced directly to obtain genomic sequences. The recombinant patterns and breakpoints of the NFLGs were determined using the Simplot software and confirmed by the maximum likelihood trees for segments using the IQ-TREE and BEAST software. The genotypic resistance profiles of the protease reverse transcriptase and integrase were resolved by the Stanford HIV drug resistance database.

Results

The six genomes shared highly similar recombinant pattern, with the CRF55_01B backbone substituted by CRF07_BC segments, therefore assigned as CRF156_0755. The evolutionary analysis of the segments showed that CRF07_BC segments were not clustered with the Chinese MSM variants in the CRF07_BC lineage. All the five NFLGs were identified with the non-nucleoside reverse-transcription inhibitors (NNRTIs) resistance mutation V179E.

Discussion

With the accumulation and evolution of recombination between CRF55_01B and CRF 07_BC, the prevalence of more recombinant strains of CRF55_01B and CRF 07_BC may occur. Therefore, it is necessary to strengthen the identification and monitoring of the recombination of CRF55_01B and CRF 07_BC.

HIV-1 Gag, Pol, and Env diversified with limited adaptation since the 1980s

Knowledge of HIV-1 global sequence diversity is critical for developing an effective prophylactic against HIV-1 infection. We developed the Hervé platform to analyze and visualize trends in HIV-1 diversification. Using Hervé, we analyzed 4,830 Env, 4,407 Gag, and 3,002 Pol publicly available independent sequences corresponding to subtypes A1, A6, B, C, D, F1, and G and circulating recombinant forms (CRFs) 01_AE, 02_AG, and 07_BC; sequences were sampled between 1980 and 2020 from 82 countries. HIV-1 diversified with a median of 1.82 amino acid substitutions per year in Env, 0.297 in Gag, and 0.779 in Pol. Yet, Env subtype B diversification plateaued post-2000. Pairwise diversity within subtypes and CRFs increased by 41.82% (range = 24.85%-54.41%) in Env, 56.93% (15.38%-89.16%) in Gag, and 46.12% (11.70%-70.57%) in Pol. Consensus sequences based on sequences sampled in each decade remained relatively stable over time. Similarly, at antibody epitope sites, only 0-8 residues that were minority variants became consensus over time in any subtype/CRF and only one known drug resistance mutation site differed from the reference (subtype G). The apparent contradiction between the fast diversification of HIV-1 and its limited adaptation illustrates that HIV-1 evolution is not directional and its consensus is at the intersection of millions of within-host selective processes occurring in a star-like manner. While a consensus sequence is a better representation of HIV-1 diversity than any individual sequence, consensus sequences have progressively become more distant from the circulating sequences they represent.

IMPORTANCE

Global surveillance of HIV-1 sequences is critical for designing relevant prophylactic and therapeutic interventions to infection. We designed an open-source platform, Hervé, for analyzing and visualizing the diversification dynamics of HIV-1 protein sequences. We characterized the evolution of over 12,000 HIV-1 Env, Gag, and Pol protein sequences from 1980-2020 and found that, despite a steady increase in intra-subtype and circulating recombinant form diversity, the most frequent residue at each site, i.e., the consensus, has varied only moderately.

HIV-1 subtypes and latent reservoirs

PURPOSE OF REVIEW

We explore the current status of research on HIV-1 subtype-specific variations and their impact on HIV-1 latency. We also briefly address the controversy surrounding the decision-making process governing the ON/OFF states of HIV-1 transcription, specifically focusing on the regulatory elements, the long terminal repeat (LTR), and Tat. Understanding the decision-making process is crucial for developing effective intervention strategies, such as the 'shock-and-kill' approach, to reactivate latent HIV-1.

RECENT FINDINGS

Attention has been drawn to subtype-specific transcription factor binding site (TFBS) variations and the possible impact of these variations on viral latency. Further, diverse subtype-specific assays have been developed to quantify the latent viral reservoirs. One interesting observation is the relatively larger latent reservoirs in HIV-1B infection than those of other viral subtypes, which needs rigorous validation. The emergence of LTR-variant viral strains in HIV-1C demonstrating significantly higher levels of latency reversal has been reported.

SUMMARY

Despite persistent and substantial efforts, latent HIV-1 remains a formidable challenge to a functional cure. Determined and continued commitment is needed to understand the ON/OFF decision-making process of HIV-1 latency, develop rigorous assays for accurately quantifying the latent reservoirs, and identify potent latency-reversing agents and cocktails targeting multiple latency stages. The review emphasizes the importance of including diverse viral subtypes in future latency research.

The origin, dissemination, and molecular networks of HIV-1 CRF65_cpx strain in Hainan Island, China

BACKGROUND

HIV-1 CRF65_cpx strain carries drug-resistant mutations, which raises concerns about its potential for causing virologic failure. The CRF65_cpx ranks as the fourth most prevalent on Hainan Island, China. However, the origin and molecular epidemiology of CRF65_cpx strains in this area remain unclear. This study aims to estimate the spatial origins and dissemination patterns of HIV-1 CRF65_cpx in this specific region.

METHODS

Between 2018 and 2021, a total of 58 pol sequences of the CRF65_cpx were collected from HIV-positive patients on Hainan Island. The available CRF65_cpx pol sequences from public databases were compiled. The HIV-TRACE tool was used to construct transmission networks. The evolutionary history of the introduction and dissemination of HIV-1 CRF65_cpx on Hainan Island were analyzed using phylogenetic analysis and the Bayesian coalescent-based approach.

RESULTS

Among the 58 participants, 89.66% were men who have sex with men (MSM). The median age was 25 years, and 43.10% of the individuals had a college degree or above. The results indicated that 39 (67.24%) sequences were interconnected within a single transmission network. A consistent expansion was evident from 2019 to 2021, with an incremental annual addition of four sequences into the networks. Phylodynamic analyses showed that the CRF65_cpx on Hainan Island originated from Beijing (Bayes factor, BF = 17.4), with transmission among MSM on Hainan Island in 2013.2 (95%HPD: 2012.4, 2019.5), subsequently leading to an outbreak. Haikou was the local center of the CRF65_cpx epidemic. This strain propagated from Haikou to other locations, including Sanya (BF > 1000), Danzhou (BF = 299.3), Chengmai (BF = 27.0) and Tunchang (BF = 16.3). The analyses of the viral migration patterns between age subgroups and risk subgroups revealed that the viral migration directions were from "25-40 years old" to "17-24 years old" (BF = 14.6) and to "over 40 years old" (BF = 17.6), and from MSM to heterosexuals (BF > 1000) on Hainan Island.

CONCLUSION

Our analyses elucidate the transmission dynamics of CRF65_cpx strain on Hainan Island. Haikou is identified as the potential hotspot for CRF65_cpx transmission, with middle-aged MSM identified as the key population. These findings suggest that targeted interventions in hotspots and key populations may be more effective in controlling the HIV epidemic.

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.

Changes in Primary HIV-1 Drug Resistance Due to War Migration from Eastern Europe

In recent years, especially as a result of war in Ukraine, enormous movements of migration to Poland from eastern European countries have been reported, including people living with Human Immunodeficiency Virus (HIV). We have conducted multi-center, prospective study, which aimed to establish HIV-1 subtype and assess the presence of primary drug resistance mutations to nucleoside reverse transcriptase inhibitors, non-nucleoside reverse transcriptase inhibitors and protease inhibitors in antiretroviral treatment naïve patients. The clinical trial recruited 117 individuals during 2 years period (2020-2022). The prevalence of HIV-1 subtype A was statistically significantly more frequent in Ukrainian, and HIV-1 subtype B in Polish patients (p < 0.05). Drug resistance mutations were detected in 44% of all cases and the comparison of presence of mutations in the analyzed groups, as well as in the subgroups of subtype A and B HIV-1 has not revealed any significant differences (p > 0.05), nevertheless Polish patients had multidrug resistance mutations more frequent (p < 0.05). The results from our trial show no increased risk of transmission of multidrug resistant HIV strains in our cohort of Ukrainian migrants.Clinical trials. Gov number NCT04636736; date of registration: November 19, 2020.

Genomic Characteristics of the New HIV-1 CRF07_BC K28E32 Variant

Accompanied with the appearance and prevalence of the new K28E32 variant among men who have sex with men, HIV-1 circulating recombinant form 07_BC (CRF07_BC) was becoming the most predominant subtype circulating in China. The K28E32 variant with five specific mutations in reverse transcriptase coding region appears to have significantly higher in vitro HIV-1 replication ability than the wild-type strain. In this study, we characterized the special mutations/substitutions in the K28E32 variant at the genomic level. Ten specific mutations that rarely appeared in other six main HIV-1 subtypes/CRFs (A-D, CRF01_AE, and CRF02_AG) were identified in the coding genes/regions of the K28E32 variant, including S77L and a novel seven-amino acid detection (32DKELYPL38) (p6Δ7) in p6, I135L in integrase, T189S in Vif, H/Y15L/F in Vpr, I264V/A and LV/LI328-329VG in gp41, and H82C and S97P in Rev. The special locations of the novel p6Δ7, and gp41 mutations I264V/A and LV/LI328-329VG in crucial protein functional domains suggest that these mutations might be functionally important to the K28E32 variant. Furthermore, eight specific substitutions were identified in Rev responsive element (RRE) of the K28E32 variant, and were revealed to increase the stability of RRE structure with a lower minimum free energy. Whether these mutations/substitutions contribute to improved transmissibility of the CRF07_BC K28E32 variant needs to be further confirmed.

Comprehensive Genetic Characterization of Four Novel HIV-1 Circulating Recombinant Forms (CRF129_56G, CRF130_A1B, CRF131_A1B, and CRF138_cpx): Insights from Molecular Epidemiology in Cyprus

Molecular investigations of the HIV-1 pol region (2253-5250 in the HXB2 genome) were conducted on sequences obtained from 331 individuals infected with HIV-1 in Cyprus between 2017 and 2021. This study unveiled four distinct HIV-1 putative transmission clusters, encompassing 19 previously unidentified HIV-1 recombinants. These recombinants, each comprising eight, three, four, and four sequences, respectively, did not align with previously established Circulating Recombinant Forms (CRFs). To characterize these novel HIV-1 recombinants, near-full-length genome sequences were successfully obtained for 16 of the 19 recombinants (790-8795 in the HXB2 genome) using an in-house-developed RT-PCR assay. Phylogenetic analyses, employing MEGAX and Cluster-Picker, along with confirmatory neighbor-joining tree analyses of subregions, were conducted to identify distinct clusters and determine subtypes. The uniqueness of the HIV-1 recombinants was evident in their exclusive clustering within generated maximum likelihood trees. Recombination analyses highlighted the distinct chimeric nature of these recombinants, with consistent mosaic patterns observed across all sequences within each of the four putative transmission clusters. Conclusive genetic characterization identified four novel HIV-1 CRFs: CRF129_56G, CRF130_A1B, CRF131_A1B, and CRF138_cpx. CRF129_56G exhibited two recombination breakpoints and three fragments of subtypes CRF56_cpx and G. Both CRF130_A1B and CRF131_A1B featured seven recombination breakpoints and eight fragments of subtypes A1 and B. CRF138_cpx displayed five recombination breakpoints and six fragments of subtypes CRF22_01A1 and F2, along with an unclassified fragment. Additional BLAST analyses identified a Unique Recombinant Form (URF) of CRF138_cpx with three additional recombination sites, involving subtype F2, a fragment of unknown subtype origin, and CRF138_cpx. Post-identification, all putative transmission clusters remained active, with CRF130_A1B, CRF131_A1B, and CRF138_cpx clusters exhibiting further growth. Furthermore, international connections were identified through BLAST analyses, linking one sequence from the USA to the CRF130_A1B strain, and three sequences from Belgium and Cameroon to the CRF138_cpx strain. This study contributes valuable insights into the dynamic landscape of HIV-1 diversity and transmission patterns, emphasizing the need for ongoing molecular surveillance and global collaboration in tracking emerging viral variants.

Identification of a HIV-1 circulating BF1 recombinant form (CRF75_BF1) of Brazilian origin that also circulates in Southwestern Europe

Introduction

The high recombinogenic potential of HIV-1 has resulted in the generation of countless unique recombinant forms (URFs) and around 120 reported circulating recombinant forms (CRFs). Here we identify through analyses of near full-length genomes (NFLG) a new HIV-1 CRF derived from subtypes B and F1.

Methods

HIV-1 protease-reverse transcriptase (Pr-RT) sequences were obtained by RT-PCR amplification from plasma RNA. Near full-length genome sequences were obtained after amplification by RT-PCR in 5 overlapping fragments. Phylogenetic sequence analyses were performed via maximum likelihood. Mosaic structures were analyzed by bootscanning and phylogenetic analyses of genome segments. Temporal and geographical estimations of clade emergence were performed with a Bayesian coalescent method.

Results

Through phylogenetic analyses of HIV-1 Pr-RT sequences obtained by us from samples collected in Spain and downloaded from databases, we identified a BF1 recombinant cluster segregating from previously reported CRFs comprising 52 viruses, most from Brazil (n = 26), Spain (n = 11), and Italy (n = 9). The analyses of NFLG genomes of 4 viruses of the cluster, 2 from Spain and 2 from Italy, allowed to identify a new CRF, designated CRF75_BF1, which exhibits a complex mosaic structure with 20 breakpoints. All 4 patients harboring CRF75_BF1 viruses studied by us had CD4+ T-cell lymphocyte counts below 220/mm3 less than one year after diagnosis, a proportion significantly higher (p = 0.0074) than the 29% found in other patients studied in Spain by us during the same period. The origin of the clade comprising CRF75_BF1 and related viruses was estimated around 1984 in Brazil, with subsequent introduction of CRF75_BF1 in Italy around 1992, and migration from Italy to Spain around 1999.

Conclusion

A new HIV-1 CRF, designated CRF75_BF1, has been identified. CRF75_BF1 is the 6th CRF of South American origin initially identified in Western Europe, reflecting the increasing relationship of South American and European HIV-1 epidemics. The finding of low CD4+ T-cell lymphocyte counts early after diagnosis in patients harboring CRF75_BF1 viruses warrants further investigation on the virulence of this variant.

Characterization of HIV-1 CRF02_AG/A3/G unique recombinant forms identified among children in Larkana, Pakistan

Co-circulation of different human immunodeficiency virus type 1 HIV-1 subtypes among infected populations can lead to the generation of new recombinants. In Pakistan, subtype A1 and CRF02_AG are the dominant strains circulating among key populations. The high prevalence of new HIV infections among the key populations highlights the possibility of recombination between the dominant strains, which can lead to the generation of new recombinants. Here, we identified a recombinant cluster composed of CRF02_AG, sub-subtype A3, and subtype G among HIV-infected children in Larkana. For the study, 10 retrospectively collected samples, with recombination signals in the pol gene, were used to perform a near full-length genome NFLG sequencing. Of the 10 samples, NFLG was successfully sequenced from seven samples. Phylogenetic analysis of the seven NFLGs showed that all recombinants formed a distinct monophyletic cluster and were distinct from known HIV-1 circulating recombinant forms CRFs. Recombination analyses showed that all seven NFLGs shared a similar recombinant structure consisting of CRF02_AG, sub-subtype A3, and subtype G, with a sub-subtype A3 fragment inserted into pol and vif regions spanning from (HXB2: 4218-5518), and a subtype G fragment inserted into vpu, rev, tat and env regions spanning from (HXB2: 5957-8250) of the CRF02_AG backbone. The identification of unique recombinant forms may indicate the presence and transmission of several co-circulating lineages in Larkana, giving rise to newer CRFs. This study also highlights the importance of continuous molecular surveillance to fully understand HIV-1 genetic diversity in Pakistan, particularly in Larkana, which is the epicenter of HIV outbreaks.

Epidemiological and molecular characterization of HBV and HCV infections in HIV-1-infected inmate population in Italy: a 2017-2019 multicenter cross-sectional study

HBV/HCV co-infection is common in HIV-1-infected prisoners. To investigate the characteristics of HIV co-infections, and to evaluate the molecular heterogeneity of HIV, HBV and HCV in prisoners, we carried-out a multicenter cross-sectional study, including 65 HIV-1-infected inmates enrolled in 5 Italian detention centers during the period 2017-2019. HIV-1 subtyping showed that 77.1% of inmates were infected with B subtype and 22.9% with non-B subtypes. Italian nationals were all infected with subtype B (93.1%), except two individuals, one infected with the recombinant form CRF72_BF1, and the other with the HIV-1 sub-subtype A6, both previously not identified in inmates of Italian nationality. Non-Italian nationals were infected with subtype B (52.6%), CRFs (36.8%) and sub-subtypes A1 and A3 (5.2%). HIV variants carrying resistance mutations to NRTI, NNRTI, PI and InSTI were found in 7 inmates, 4 of which were never exposed to the relevant classes of drugs associated with these mutations. HBV and/or HCV co-infections markers were found in 49/65 (75.4%) inmates, while 27/65 (41.5%) showed markers of both HBV and HCV coinfection. Further, Italian nationals showed a significant higher presence of HCV markers as compared to non-Italian nationals (p = 0.0001). Finally, HCV phylogenetic analysis performed in 18 inmates revealed the presence of HCV subtypes 1a, 3a, 4d (66.6%, 16.7% and 16.7%, respectively). Our data suggest the need to monitor HIV, HBV and HCV infections in prisons in order to prevent spreading of these viruses both in jails and in the general population, and to implement effective public health programs that limit the circulation of different genetic forms as well as of viral variants with mutations conferring resistance to treatment.

Alarming rise in HIV cases in Pakistan: Challenges and future recommendations at hand

Background

Human immunodeficiency virus (HIV) is a retrovirus that suppresses the immune system by reducing the CD4+ T lymphocytes level. It has become a global challenge with fast prevalence ratio. Like other developing countries, Pakistan is also struggling for overcoming this viral disease since very first reported case in 1987.

Aim

To update the society on the alarming rise in HIV cases in Pakistan: challenges and future recommendations at hand.

Materials and Methods

The review paper utilized different search engines such pubmed central, scopus, web of science, google scholar etc. to conduct this review paper.

Results

Lack of awareness, low literacy rate, practice of unhygienic equipment in healthcare departments, unstable economy, and unsafe sexual practices are the major factors behind the increasing rate of AIDS in Pakistan.

Conclusion

By regulating healthcare practices and policies, promoting psychological counseling to HIV positive patients, educating the society and minimizing commercial sex practices, Pakistan can overcome this viral disease.

Global associations of key populations with HIV-1 recombinants: a systematic review, global survey, and individual participant data meta-analysis

Introduction

Global HIV infections due to HIV-1 recombinants are increasing and impede prevention and treatment efforts. Key populations suffer most new HIV infections, but their role in the spread of HIV-1 recombinants is unknown. We conducted a global analysis of the associations between key populations and HIV-1 recombinants.

Methods

We searched PubMed, EMBASE, CINAHL, and Global Health for HIV-1 subtyping studies published from 1/1/1990 to 31/12/2015. Unpublished data was collected through a global survey. We included studies with HIV-1 subtyping data of key populations collected during 1990-2015. Key populations assessed were heterosexual people (HET), men who have sex with men (MSM), people who inject drugs (PWID), vertical transmissions (VERT), commercial sex workers (CSW), and transfusion-associated infections (BLOOD). Logistic regression was used to determine associations of key populations with HIV-1 recombinants. Subgroup analyses were performed for circulating recombinant forms (CRFs), unique recombinant forms (URFs), regions, and time periods.

Results

Eight hundred and eighty five datasets including 77,284 participants from 83 countries were included. Globally, PWID were associated with the greatest odds of recombinants and CRFs (OR 2.6 [95% CI 2.46-2.74] and 2.99 [2.83-3.16]), compared to HET. CSW were associated with increased odds of recombinants and URFs (1.59 [1.44-1.75] and 3.61 [3.15-4.13]). VERT and BLOOD were associated with decreased odds of recombinants (0.58 [0.54-0.63] and 0.43 [0.33-0.56]). MSM were associated with increased odds of recombinants in 2010-2015 (1.43 [1.35-1.51]). Subgroup analyses supported our main findings.

Discussion

As PWID, CSW, and MSM are associated with HIV-1 recombinants, increased preventative measures and HIV-1 molecular surveillance are crucial within these key populations.

Systematic review registration

PROSPERO [CRD42017067164].

Dolutegravir-associated resistance mutations after first-line treatment failure in Brazil

BACKGROUND

Since January 2017, the recommended first-line antiretroviral regimen in Brazil is the fixed-dose combination of tenofovir plus lamivudine with dolutegravir (TL + D). According to the literature, integrase resistance-associated mutations (INRAMs) are rarely found upon virologic failure to first-line dolutegravir plus two nucleoside reverse transcriptase inhibitors. We evaluated the HIV antiretroviral genotypic resistance profile of patients referred for genotyping in the public health system who failed first-line TL + D after at least six months of therapy on or before December 31, 2018.

METHODS

HIV Sanger sequences of the pol gene were generated from plasma of patients with confirmed virologic failure to first-line TL + D in the Brazilian public health system before December 31, 2018.

RESULTS

One hundred thirteen individuals were included in the analysis. Major INRAMs were detected in seven patients (6.19%), four with R263K, one with G118R, one with E138A, and one with G140R. Four patients with major INRAMs also had the K70E and M184V mutations in the RT gene. Sixteen (14.2%) additional individuals presented minor INRAMs, and five (4,42%) patients had both major and minor INRAMS. Thirteen (11.5%) patients also presented mutations in the RT gene selected by tenofovir and lamivudine, including four with both the K70E and M184V mutations and four with only M184V. The integrase mutations L101I and T124A, which are in the in vitro pathway for integrase inhibitor resistance, were found in 48 and 19 patients, respectively. Mutations not related to TL + D, thus probable transmitted resistance mutations (TDR), were present in 28 patients (24.8%): 25 (22.1%) to nucleoside reverse transcriptase inhibitors, 19 (16.8%) to non-nucleoside reverse transcriptase inhibitors, and 6 (5.31%) to protease inhibitors.

CONCLUSIONS

In marked contrast to previous reports, we report a relatively high frequency of INRAMs among selected patients failing first-line TL + D in the public health system in Brazil. Possible reasons for this discrepancy include delays in detecting virologic failure, patients inadvertently on dolutegravir monotherapy, TDR, and/or infecting subtype.

Risk Factors in HIV-1 Positive Patients on the Intensive Care Unit: A Single Center Experience from a Tertiary Care Hospital

HIV-positive patients with acquired immunodeficiency syndrome (AIDS) often require treatment on intensive care units (ICUs). We aimed to present data from a German, low-incidence region cohort, and subsequently evaluate factors measured during the first 24 h of ICU stay to predict short- and long-term survival, and compare with data from high-incidence regions. We documented 62 patient courses between 2009 and 2019, treated on a non-operative ICU of a tertiary care hospital, mostly due to respiratory deterioration and co-infections. Of these, 54 patients required ventilatory support within the first 24 h with either nasal cannula/mask (n = 12), non-invasive ventilation (n = 16), or invasive ventilation (n = 26). Overall survival at day 30 was 77.4%. While ventilatory parameters (all p < 0.05), pH level (c/o 7.31, p = 0.001), and platelet count (c/o 164,000/µL, p = 0.002) were significant univariate predictors of 30-day and 60-day survival, different ICU scoring systems, such as SOFA score, APACHE II, and SAPS 2 predicted overall survival (all p < 0.001). Next to the presence or history of solid neoplasia (p = 0.026), platelet count (HR 6.7 for <164,000/µL, p = 0.020) and pH level (HR 5.8 for <7.31, p = 0.009) remained independently associated with 30-day and 60-day survival in multivariable Cox regression. However, ventilation parameters did not predict survival multivariably.

HIV-1 Tat amino acid residues that influence Tat-TAR binding affinity: a scoping review

HIV-1 remains a global health concern and to date, nearly 38 million people are living with HIV. The complexity of HIV-1 pathogenesis and its subsequent prevalence is influenced by several factors including the HIV-1 subtype. HIV-1 subtype variation extends to sequence variation in the amino acids of the HIV-1 viral proteins. Of particular interest is the transactivation of transcription (Tat) protein due to its key function in viral transcription. The Tat protein predominantly functions by binding to the transactivation response (TAR) RNA element to activate HIV-1 transcriptional elongation. Subtype-specific Tat protein sequence variation influences Tat-TAR binding affinity. Despite several studies investigating Tat-TAR binding, it is not clear which regions of the Tat protein and/or individual Tat amino acid residues may contribute to TAR binding affinity. We, therefore, conducted a scoping review on studies investigating Tat-TAR binding. We aimed to synthesize the published data to determine (1) the regions of the Tat protein that may be involved in TAR binding, (2) key Tat amino acids involved in TAR binding and (3) if Tat subtype-specific variation influences TAR binding. A total of thirteen studies met our inclusion criteria and the key findings were that (1) both N-terminal and C-terminal amino acids outside the basic domain (47-59) may be important in increasing Tat-TAR binding affinity, (2) substitution of the amino acids Lysine and Arginine (47-59) resulted in a reduction in binding affinity to TAR, and (3) none of the included studies have investigated Tat subtype-specific substitutions and therefore no commentary could be made regarding which subtype may have a higher Tat-TAR binding affinity. Future studies investigating Tat-TAR binding should therefore use full-length Tat proteins and compare subtype-specific variations. Studies of such a nature may help explain why we see differential pathogenesis and prevalence when comparing HIV-1 subtypes.

Needle-Free Devices and CpG-Adjuvanted DNA Improve Anti-HIV Antibody Responses of Both DNA and Modified Vaccinia Ankara-Vectored Candidate Vaccines

The combination of mosaic Gag and CAP256 envelope in an HIV vaccine regimen comprising DNA prime and modified vaccinia Ankara (MVA) boost followed by protein boost has previously been shown to generate robust autologous Tier 2 neutralizing antibodies (nAbs) in rabbits. Further refinements of this strategy have been investigated to improve antibody responses. The delivery of both DNA and recombinant MVA vaccines with a needle-free device was compared to delivery by injection, and the effect of formulating the DNA vaccine with adjuvant CpG ODN 1826 was determined. The Pharmajet Stratis^® needle-free injection device (PharmaJet, Golden, CO, USA) improved binding antibody responses to the DNA vaccine as well as both binding and neutralizing antibody responses to the MVA vaccines. Formulation of the DNA vaccines with CpG adjuvant further improved the antibody responses. A shortened vaccination regimen of a single DNA inoculation followed by a single MVA inoculation did not elicit Tier 1B nor Tier 2 neutralization responses as produced by the two DNA, followed by two MVA vaccination regimen. This study showed the immunogenicity of HIV DNA and MVA vaccines administered in a DDMM regimen could be improved using the PharmaJet Stratis needle-free injection device and formulation of the DNA vaccines with CpG adjuvant.

Tracing the relationship among HIV-1 sub-subtype F1 strains: a phylodynamic perspective

BACKGROUND

The human immunodeficiency virus type 1, F1 sub-subtype (HIV-1 F1) circulates in three continents: Africa, Europe, and South America. In Brazil, this sub-subtype co-circulates with subtypes B and C and several recombinant forms, mainly BF1 variants.

OBJECTIVES

This study aimed to reconstruct the dynamic history of HIV-1 F1 in Brazil.

METHODS

HIV-1 near full-length genome and pol gene nucleotide sequences available in public databases were assembled in two datasets (POL671 and NFLG53) to cover the largest number of F1 sub-subtype sequences. Phylodynamic and temporal analyses were performed.

FINDINGS

Two main strains of the F1 sub-subtype are circulating worldwide. The first (F1.I) was found among Brazilian samples (75%) and the second (F1.II) among Romanian (62%) and other European and African isolates. The F1 subtype epidemic in Brazil originated from a single entry into the country around 1970. This ancestral sample is related to samples isolated in European countries (France, Finland, and Belgium), which are possibly of African origin. Moreover, further migration (1998 CI: 1994-2003) of strains from Brazil to Europe (Spain and the UK) was observed. Interestingly, all different recombinant BF patterns found, even those from outside Brazil, present the same F1 lineage (F1.I) as an ancestor, which could be related to the acquisition of adaptive advantages for the recombinant progenies.

MAIN CONCLUSIONS

These findings are important for the understanding of the origin and dynamics of the F1 sub-subtype and a consequent better and greater understanding of the HIV-1 F1 and BF epidemic that still spreads from Brazil to other countries.

Short-term associations of PM2.5 and PM2.5 constituents with immune biomarkers: A panel study in people living with HIV/AIDS

Studies on associations of fine particulate matter (PM_2.5) with immunity in people living with HIV/AIDS (PLWHA) were absent. We aimed to explore whether changes of immune biomarkers were associated with short-term exposure to PM_2.5 in PLWHA. Based on a panel study in Wuhan, we selected 163 PLWHA as participants with up to 4 repeated visits from March 2020 to January 2021. Immune biomarkers, including CD4⁺T cell count, CD8⁺T cell count, HIV viral load (VL) and CD4⁺T/CD8⁺T ratio were tested for all participants at each visit. Residential exposures of PM_2.5 and PM_2.5 constituents for each participant were assessed using spatial-temporal models. Linear mixed-effect models and general linear mixed models were applied to evaluate the associations between PM_2.5 and immune biomarkers. To estimate the combined effect of PM_2.5 constituents, weighted quantile sum regression and Bayesian kernel machine regression were employed. Each 10 μg/m³ increase of 7-day average PM_2.5 concentrations was associated with an 8.75 cells/mm³ (95%CI: -15.55, -1.98) decrease in CD4⁺T cell count and a 92% (OR: 1.92, 95%CI: 1.43, 2.58) increased odds ratio of detectable HIV VL. However, the odds ratio of inverted CD4⁺T/CD8⁺T was only positively associated with PM_2.5 concentrations at lag2 day (OR:1.27, 95%CI:1.02, 1.57). CD4⁺T may be a potential mediator between PM_2.5 and detectable HIV VL with 3.83% mediated proportion. Besides, the combined effect of PM_2.5 chemical constituents indicated that NO₃^- and SO₄^2- were the main constituents in reducing CD4⁺T cell count and increasing odds ratio of detectable HIV VL. Our finding revealed that short-term exposure to PM_2.5 was negatively associated with CD4⁺T cell count but positively related to the odds ratio of detectable HIV VL in PLWHA. This research may provide new evidence in associations between PM_2.5 and immune biomarkers as well as improving prognosis of PLWHA.

Identification and genomic characterization of a novel HIV-1 unique recombinant form (CRF01_AE/CRF07_BC) in Zhejiang Province, China

Mutation and recombination are important mechanisms leading to the frequent evolution and genetic diversity of viruses as HIV-1. In this study, we identified the near full-length genomic characterization of a novel HIV-1 unique recombinant form (URF) strain (Sample ID: ZJ20202195/ZJ/CHN/2020, hereafter referred to as ZJ20202195) isolated during the HIV-1 molecular surveillance in 2020 in Zhejiang Province, China, through different recombination analysis tools and phylogenetic analysis. Our results amply proved that the near full-length genome (NFLG) sequence of ZJ20202195 was a novel HIV-1 unique recombinant form (URF) consisting of CRF01_AE and CRF07_BC subtype, and delimited three recombinant segments, of which the Segment I (HXB2:776-5559 nucleotide (nt)) and Segment III (HXB2:6224-9412 nt) were mainly originated from CRF01_AE cluster g4a strains prevalent in China and Segment II (HXB2:5560-6223 nt) was from CRF07_BC subtype. Overall, our findings provide insight and a scientific basis in the genetic diversity and accurate determination of HIV-1 recombinant strains in China.

Identification of two novel HIV-1 circulating recombinant forms of CRF111_01C and CRF116_0108 in southwestern Yunnan, China

Yunnan, the region hardest hit by HIV/AIDS in China, is also an area with the most abundant HIV-1 genetic diversity. A large number of novel HIV-1 circulating recombinant forms (CRFs) and unique recombinants were identified among injection drug users in Yunnan; however, few were found among sexual contacts. Here, we obtained 15 near full-length genome sequences (NFLGs) from HIV-1 seropositive heterosexual contacts in Yunnan who received antiretroviral therapy during the period from 2014 to 2016. Phylogenetic analysis showed that six NFLGs belonged to CRF01_AE (n = 3) and CRF106_cpx (n = 3), and the other nine sequences were novel inter-subtype recombinants. Of the recombinants, two novel CRFs (CRF111_01 C (n = 4) and CRF116_0108 (n = 4)) and one CRF106_cpx variant (n = 1) were identified. CRF111_01 C had a CRF01_AE backbone with seven subtype C fragments inserted into the gag, pol, vif, env, nef and 3'LTR regions. CRF116_0108 had a CRF08_BC backbone with a CRF01_AE fragment inserted into the pol, tat, rev, vif, vpr, vpu and env regions. Phylogeographic analyses estimated that CRF111_01 C and CRF116_0108 originated approximately 1995.7-1998.6 and 1991.7-1993.7, respectively. These identifications of two novel HIV-1 CRFs highlighted the importance of continuous surveillance in heterosexual contacts and other high-risk groups in this region and the surrounding regions.

HIV-1 genetic diversity and transmitted drug resistance to integrase strand transfer inhibitors among recently diagnosed adults in Porto Alegre, South Brazil

OBJECTIVES

The HIV-1 genetic diversity and the presence of transmitted drug resistance mutations (TDRMs) against integrase strand transfer inhibitors (INSTIs) were assessed sequencing samples of antiretroviral (ARV)-naive HIV-1-infected individuals from South Brazil.

METHODS

Viral RNA from 42 ART-naive individuals was submitted to complete HIV-1 integrase gene amplification by RT-PCR and sequencing.

RESULTS

Viral strains carrying TDRMs against INSTIs were not detected in the present study. However, the polymorphisms L74M and L74I were each observed in 4.8% of the individuals. These accessory mutations have been reported as putative causes of TDRMs in ART with raltegravir, but only when associated with additional major mutations. When submitted to HIV-1 subtyping, 50% were classified as subtype C, 21% as recombinant BC, 19% as subtype B, 4.8% as subtype F1 and 4.8% as recombinant CF1.

CONCLUSIONS

All 42 ARV-naive individuals were apparently susceptible to INSTIs, included in the Brazilian therapeutic guideline since 2009. To the best of our knowledge, this is the first study to evaluate TDRMs against INSTIs in Brazil. The most prevalent HIV-1 subtypes were subtype C, followed by the recombinant BC and subtype B, which is in agreement with previous studies. However, the presence of subtype F1 and recombinant CF1 reported herein was not observed in previous studies.

HIV Infection and Oral Manifestations: An Update

Human immunodeficiency virus (HIV) causes a complete depletion of the immune system; it has been a major health issue around the world since the 1980s, and due to the reduction of CD4+ T lymphocytes levels, it can trigger various opportunistic infections. Oral lesions are usually accurate indicators of immunosuppression because these oral manifestations may occur as a result of the compromised immune system caused by HIV infection; therefore, oral lesions might be initial and common clinical features in people living with HIV. So, it is necessary to evaluate and understand the mechanism, prevalence, and risk factors of oral lesions to avoid the increase morbidity among those with oral diseases.

Identification of a new circulating recombinant form of human immunodeficiency virus type 1, CRF124_cpx involving subtypes A, G, H, and CRF27_cpx in Angola

Angola, located in Central Africa, has around 320,000 (270,000–380,000) people living with human immunodeficiency virus (HIV)/AIDS, equivalent to 1% of the country’s population at the end of 2021. A previous study conducted in 2012, using Angolan samples collected between 2008 and 2010 revealed a high prevalence of HIV-1 recombinants, around 42% of sequences, with 21% showing the same UH profile in partial pol region which were grouped into a monophyletic cluster with high bootstrap support. Thus, the objective of the present work was to obtain complete genomes of those sequences and characterize them, aiming at a description of a new circulating recombinant form (CRF). Whole blood from nine HIV-1 UH pol-infected individuals had their genomic DNA extracted, and nested PCR was used to amplify seven overlapping fragments targeting the full-length HIV-1 genome. The final classification was based on maximum likelihood trees, and recombination analyses were performed using a bootscan from the Simplot program. BLAST and Los Alamos Database inspections were used to search other similar H-like pol sequences. Complete genome amplification was possible for three samples, partial genomes were obtained for the other three, and only pol was available for the remaining three sequences. Bootscan analysis of the two whole-genome and three partial genome sequences retrieved from people living with HIV/AIDS (PLHIVA) without epidemiological linkage showed the same complex recombination profile involving HIV-1 subtypes A/G/H/CRF27_cpx, with a total of six recombinant breakpoints, aiming to classify a new HIV-1 CRF124_cpx. We found no other full-length HIV-1 genomes with the same mosaic profile; however, we identified 33 partial pol sequences, mainly sampled from Angola between 2001 to 2019, with the same H-like profile. Bayesian analysis of H and H-like pol sequences indicates that CRF124_cpx probably originated in Angola at mid-1970s, indicating that this CRF has been circulating in the country for a long time. In summary, our study describes a new CRF circulating principally in Angola and highlights the importance of continuing molecular surveillance studies, especially in countries with high molecular diversity of HIV.

Control of HIV-1 Replication by CD8+ T Cells Specific for Two Novel Pol Protective Epitopes in HIV-1 Subtype A/E Infection

Although many HIV-1-specific CD8⁺ T cell epitopes have been identified and used in various HIV-1 studies, most of these epitopes were derived from HIV-1 subtypes B and C. Only 17 well-defined epitopes, none of which were protective, have been identified for subtype A/E infection. The roles of HIV-1-specific T cells have been rarely analyzed for subtype A/E infection. In this study, we identified six novel HLA-B*15:02-restricted optimal HIV-1 subtype A/E epitopes and then analyzed the presentation of these epitopes by HIV-1 subtype A/E virus-infected cells and the T cell responses to these epitopes in treatment-naive HIV-1 subtype A/E-infected HLA-B*15:02⁺ Vietnamese individuals. Responders to the PolTY9 or PolLF10 epitope had a significantly lower plasma viral load (pVL) than nonresponders among HLA-B*15:02⁺ individuals, whereas no significant difference in pVL was found between responders to four other epitopes and nonresponders. The breadth of T cell responses to these two Pol epitopes correlated inversely with pVL. These findings suggest that HLA-B*15:02-restricted T cells specific for PolTY9 and PolLF10 contribute to the suppression of HIV-1 replication in HLA-B*15:02⁺ individuals. The HLA-B*15:02-associated mutation Pol266I reduced the recognition of PolTY9-specific T cells in vitro but did not affect HIV-1 replication by PolTY9-specific T cells in Pol266I mutant virus-infected individuals. These findings indicate that PolTY9-specific T cells suppress replication of the Pol266I mutant virus even though the T cells selected this mutant. This study demonstrates the effective role of T cells specific for these Pol epitopes to control circulating viruses in HIV-1 subtype A/E infection. IMPORTANCE It is expected that HIV-1-specific CD8⁺ T cells that effectively suppress HIV-1 replication will contribute to HIV-1 vaccine development and therapy to achieve an HIV cure. T cells specific for protective epitopes were identified in HIV-1 subtype B and C infections but not in subtype A/E infection, which is epidemic in Southeast Asia. In the present study, we identified six T cell epitopes derived from the subtype A/E virus and demonstrated that T cells specific for two Pol epitopes effectively suppressed HIV-1 replication in treatment-naive Vietnamese individuals infected with HIV-1 subtype A/E. One of these Pol protective epitopes was conserved among circulating viruses, and one escape mutation was accumulated in the other epitope. This mutation did not critically affect HIV-1 control by specific T cells in HIV-1 subtype A/E-infected individuals. This study identified two protective Pol epitopes and characterized them in cases of HIV-1 subtype A/E infection.

HIV-1 Circulating Recombinant Forms (CRFs) and Unique Recombinant Forms (URFs) in Israel, 2010-2018

Monitoring HIV-1 circulating recombinant forms (CRFs) and unique recombinant forms (URFs) is important for disease surveillance. Recombination may affect prevention efforts and interfere with the diagnosis and treatment of HIV-1 infection. Here, we characterized the epidemiology of HIV-1 CRFs and URFs in Israel. Partial pol sequences from treatment naïve patients diagnosed in 2010−2018 were assessed using the recombinant identification program (RIP), the recombinant detection program (RDP5), and using the maximum-likelihood phylogenetic method, using 410 reference sequences obtained from the Los Alamos database. CRFs and URFs were identified in 11% (213/1940) of all sequenced cases. The median age at diagnosis was 38 (30−47) years, 61% originated from Israel, and 82% were male. The most common were CRF02_AG (30.5%), CRF01_AE (16.9%), and the more complex forms CRF01_AE/CRF02_AG/A3 (10.8%) and B/F1 (7%). A significant increase in their overall proportion was observed in recent years (8.1% in 2010−2012, 20.3% in 2016−2018, p < 0.001). This increase was most prominent in individuals carrying CRF02_AG (2.5% in 2010−2015, 9.8% in 2016−2018, p < 0.001). Men who have sex with men (MSM) was the most common risk group; however, those infected with the secondary recombinant CRF02_AG/A6 were mainly injecting drug users (IDUs). The most common resistance mutations were K103N (5/213, 2.3%) and E138A (18/213, 8.5%) in the reverse transcriptase. Only E138A was more frequent in the recombinants compared with the classic subtypes and was significantly associated with a specific secondary CRF, CRF02_AG/A4. We concluded that CRFs and URFs were mainly detected in Israeli-born MSM and that an increase in the overall proportion of such HIV-1 sequences could be observed in more recent years.

Detection of HIV-1 Transmission Clusters from Dried Blood Spots within a Universal Test-and-Treat Trial in East Africa

The Sustainable East Africa Research in Community Health (SEARCH) trial was a universal test-and-treat (UTT) trial in rural Uganda and Kenya, aiming to lower regional HIV-1 incidence. Here, we quantify breakthrough HIV-1 transmissions occurring during the trial from population-based, dried blood spot samples. Between 2013 and 2017, we obtained 549 gag and 488 pol HIV-1 consensus sequences from 745 participants: 469 participants infected prior to trial commencement and 276 SEARCH-incident infections. Putative transmission clusters, with a 1.5% pairwise genetic distance threshold, were inferred from maximum likelihood phylogenies; clusters arising after the start of SEARCH were identified with Bayesian time-calibrated phylogenies. Our phylodynamic approach identified nine clusters arising after the SEARCH start date: eight pairs and one triplet, representing mostly opposite-gender linked (6/9), within-community transmissions (7/9). Two clusters contained individuals with non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance, both linked to intervention communities. The identification of SEARCH-incident, within-community transmissions reveals the role of unsuppressed individuals in sustaining the epidemic in both arms of a UTT trial setting. The presence of transmitted NNRTI resistance, implying treatment failure to the efavirenz-based antiretroviral therapy (ART) used during SEARCH, highlights the need to improve delivery and adherence to up-to-date ART recommendations, to halt HIV-1 transmission.

HIV-1 Subtype Shift in the Philippines is Associated with High Transmitted Drug Resistance, High Viral Loads and Fast Immunologic Decline

BACKGROUND

The Philippines has one of the fastest growing HIV epidemics. A subtype shift from B to CRF01_AE may have contributed to the rise in cases. We undertook a genotyping and transmitted drug resistance (TDR) study to determine if the dominant subtype has any advantages in resistance and transmission.

METHODS

Treatment-naive Filipinos living with HIV were recruited from two large government treatment hubs March 2016 to August 2018. HIV-1 viral load, CD4 count, genotyping and TDR testing were performed. Demographic and clinical data was collected and compared across subtypes.

RESULTS

Two hundred ninety-eight Filipino PLHIV were recruited. Median CD4 count was 143 cells/µL and HIV viral load was 2,345,431 copies/mL. Sanger-based sequencing showed 230/298 (77.2%) had subtype CRF01_AE, 41 (13.8%) subtype B, and the rest other subtypes or recombinants. Overall TDR was 11.7%. TDR was associated with lower viral loads and no previous HIV testing. CRF01_AE had a higher likelihood of a viral load >100,000 copies/mL and having a baseline CD4 count <50 cells/mm3.

CONCLUSIONS

TDR in the Philippines is high at 11.7%. CRF01_AE was observed to have a higher baseline viral load and lower CD4 counts compared to other co-circulating subtypes. Further research needs to confirm this observation since it suggests that CRF01_AE may have a survival advantage that led to replacement of subtype B as the dominant subtype. Drug-resistance testing is recommended in the Philippines when initiating NNRTI-based anti-retroviral therapy but may not be necessary for INSTI-based regimens.

The future of long-acting cabotegravir plus rilpivirine therapy: Deeds and misconceptions

HIV infection is currently managed as a chronic disease thanks to the improvement of antiretroviral therapy (ART). In this light, treatment is lifelong. Therefore, switch to new regimens is a natural event during long-term therapies to avoid problems related to toxicity, adherence, failure and potential selection of drug resistance. In this context, coformulations of multiple agents in one pill and the development of novel drug classes and drugs with high genetic barrier to resistance are now available. The recent approval of the long-acting once monthly or bimonthly injectable combination of the second-generation strand transfer integrase inhibitor (InSTI) cabotegravir (CAB) together with the non-nucleoside reverse transcriptase inhibitor (NNRTI) rilpivirine (RPV) represents the most recent achievement along the line of potent and convenient ART. Several pivotal trials (such as LATTE-2, ATLAS, FLAIR, and ATLAS-2M) largely demonstrated the high efficacy and safety of this long-acting formulation used as an induction-maintenance strategy. Few confirmed virological failures (CVF) have been observed. The combination of at least two baseline factors among HIV-1 subtype A6/A1, a body mass index ≥30 kg/m², and RPV resistance associated mutations, was associated with an increased risk of CVF at week 48. Available data indicate that this long-acting therapeutic strategy is attractive and potent, thus defining the most appropriate patient and how to handle practical issues is warranted.

Multiple CRF01_AE/CRF07_BC Recombinants Enhanced the HIV-1 Epidemic Complexity Among MSM in Shenyang City, Northeast China

The transmission of Unique Recombinant Forms (URFs) has complicated the molecular epidemic of HIV-1. This increasing genetic diversity has implications for prevention surveillance, diagnosis, and vaccine design. In this study, we characterized the HIV-1 URFs from 135 newly diagnosed HIV-1 infected cases between 2016 and 2020 in Shenyang, northeast China and analyzed the evolutionary relationship of them by phylogenetic and recombination approaches. Among 135 URFs, we found that the CRF01_AE/CRF07_BC recombinants were the most common (81.5%, 110/135), followed by CRF01_AE/B (11.9%, 16/135), B/C (3.7%, 5/135), and others (3.0%, 4/135). 94.8% (128/135) of patients infected by URFs were through homosexual contact. Among 110 URFs_0107, 60 (54.5%) formed 11 subclusters (branch support value = 1) and shared the consistent recombination structure, respectively. Four subclusters have caused small-scale spread among different high-risk populations. Although the recombination structures of URFs_0107 are various, the hotspots of recombinants gathered between position 2,508 and 2,627 (relative to the HXB2 position). Moreover, the CRF07_BC and CRF01AE fragments of URFs_0107 were mainly derived from the MSM population. In brief, our results reveal the complex recombinant modes and the high transmission risk of URFs_0107, which calls for more attention on the new URFs_0107 monitoring and strict control in the areas led by homosexual transmission route.

Antiretroviral Imprints and Genomic Plasticity of HIV-1 pol in Non-clade B: Implications for Treatment

Combinational antiretroviral therapy (cART) is the most effective tool to prevent and control HIV-1 infection without an effective vaccine. However, HIV-1 drug resistance mutations (DRMs) and naturally occurring polymorphisms (NOPs) can abrogate cART efficacy. Here, we aimed to characterize the HIV-1 pol mutation landscape in Cameroon, where highly diverse HIV clades circulate, and identify novel treatment-associated mutations that can potentially affect cART efficacy. More than 8,000 functional Cameroonian HIV-1 pol sequences from 1987 to 2020 were studied for DRMs and NOPs. Site-specific amino acid frequencies and quaternary structural features were determined and compared between periods before (≤2003) and after (2004-2020) regional implementation of cART. cART usage in Cameroon induced deep mutation imprints in reverse transcriptase (RT) and to a lower extent in protease (PR) and integrase (IN), according to their relative usage. In the predominant circulating recombinant form (CRF) 02_AG (CRF02_AG), 27 canonical DRMs and 29 NOPs significantly increased or decreased in RT during cART scale-up, whereas in IN, no DRM and only seven NOPs significantly changed. The profound genomic imprints and higher prevalence of DRMs in RT compared to PR and IN mirror the dominant use of reverse transcriptase inhibitors (RTIs) in sub-Saharan Africa and the predominantly integrase strand transfer inhibitor (InSTI)-naïve study population. Our results support the potential of InSTIs for antiretroviral treatment in Cameroon; however, close surveillance of IN mutations will be required to identify emerging resistance patterns, as observed in RT and PR. Population-wide genomic analyses help reveal the presence of selective pressures and viral adaptation processes to guide strategies to bypass resistance and reinstate effective treatment.

Updated HIV-1 Consensus Sequences Change but Stay Within Similar Distance From Worldwide Samples

HIV consensus sequences are used in various bioinformatic, evolutionary, and vaccine related research. Since the previous HIV-1 subtype and CRF consensus sequences were constructed in 2002, the number of publicly available HIV-1 sequences have grown exponentially, especially from non-EU and US countries. Here, we reconstruct 90 new HIV-1 subtype and CRF consensus sequences from 3,470 high-quality, representative, full genome sequences in the LANL HIV database. While subtypes and CRFs are unevenly spread across the world, in total 89 countries were represented. For consensus sequences that were based on at least 20 genomes, we found that on average 2.3% (range 0.8–10%) of the consensus genome site states changed from 2002 to 2021, of which about half were nucleotide state differences and the rest insertions and deletions. Interestingly, the 2021 consensus sequences were shorter than in 2002, and compared to 4,674 HIV-1 worldwide genome sequences, the 2021 consensuses were somewhat closer to the worldwide genome sequences, i.e., showing on average fewer nucleotide state differences. Some subtypes/CRFs have had limited geographical spread, and thus sampling of subtypes/CRFs is uneven, at least in part, due to the epidemiological dynamics. Thus, taken as a whole, the 2021 consensus sequences likely are good representations of the typical subtype/CRF genome nucleotide states. The new consensus sequences are available at the LANL HIV database.

Intersecting Pandemics of HIV and SARS-CoV-2: Commentary on the Special Issue

PURPOSE OF REVIEW

This editorial introduces this special Global Health Section on the interface of the HIV/AIDS and COVID-19 pandemics.

RECENT FINDINGS

Authors of articles in this special issue take on a variety of topics that capture how the acute COVID-19 pandemic affected global efforts towards HIV control, and how co-infection, stigma, and social determinants of disease have affected populations on multiple continents. Two historic pandemics -- HIV/AIDS and COVID-19 -- have affected the world in our lifetimes at a level reminiscent of the 1918-1919 H1N1 influenza pandemic. We have much to learn from both experiences to optimize pandemic disease control, prevention, and management.

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.