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

Epidemic trend, genetic characteristics, and transmission networks of HIV-1 among treatment-naive men who have sex with men in Hebei province, China

Introduction

Homosexual transmission has contributed greatly to the current HIV-1 epidemic in Hebei province, China. Dolutegravir (DTG) will be conditionally used as a component of free antiretroviral therapy (ART) according to manual for national free anti-AIDS treatment drugs (2023 edition) issued by China in June 2023. However, current genetic characteristics and pretreatment drug resistance (PDR) to proteinase inhibitors (PIs), reverse transcriptase inhibitors (RTs) and integrase strand transfer inhibitors (INSTIs) of HIV-1 in this population have remained unclear.

Methods

Serial consecutive cross-sectional analyses for HIV- 1 infection trend, genetic characteristics, PDR and molecular transmission networks were conducted from 2018 to 2022. All of participants were HIV-1- infected MSM newly diagnosed at the HIV surveillance points (HSPs) in Hebei, China. Evidence of PDR was confirmed using the world health organization (WHO) list for surveillance of drug resistance mutations.

Results

In this study, a total of 14 HIV-1 subtypes were circulating in the HSPs of Hebei province, China. CRF01_ AE (51.9%, 350/675), CRF07_BC (30.4%, 205/675), B (6.2%, 42/675) and URFs (5.8%, 39/675) were the four most predominant subtypes among MSM. And, CRF07_BC (r > 0) and URFs (r > 0) indicated an increasing trend, respectively; however, CRF01_AE (r < 0) showed a decline trend. The overall prevalence of HIV-1 PDR showed a substantial increase from 6.3% in 2018 to 7.9% in 2022. The prevalence of NNRTI-PDR was the highest (5.8%, 39/675), followed by INSTIs (2.4%, 16/675), NRTIs (0.6%, 4/675) and PIs (0.3%, 2/675). Furthermore, extensive HIV-1 strains bearing PDR were circulating in the MSM population via molecular transmission networks for major HIV-1 subtypes, especially CRF01_AE and CRF07_BC.

Discussion

Our findings reflect that HIV-1 epidemic in the MSM population is complex and severe in Hebei, China. Therefore, it is urgent for us to implement more effective intervention measures to limit the further dissemination of HIV-1, especially the spread of HIV-1 INSTI-PDR strains.

Rev Protein Diversity in HIV-1 Group M Clades

The HIV-1 Rev protein expressed in the early stage of virus replication is involved in the nuclear export of some forms of virus RNA. Naturally occurring polymorphisms in the Rev protein could influence its activity. The association between the genetic features of different virus variants and HIV infection pathogenesis has been discussed for many years. In this study, Rev diversity among HIV-1 group M clades was analyzed to note the signatures that could influence Rev activity and, subsequently, clinical characteristics. From the Los Alamos HIV Sequence Database, 4962 Rev sequences were downloaded and 26 clades in HIV-1 group M were analyzed for amino acid changes, conservation in consensus sequences, and the presence of clade-specific amino acid substitutions (CSSs) and the Wu-Kabat protein variability coefficient (WK). Subtypes G, CRF 02_AG, B, and A1 showed the largest amino acid changes and diversity. The mean conservation of the Rev protein was 80.8%. In consensus sequences, signatures that could influence Rev activity were detected. In 15 out of 26 consensus sequences, an insertion associated with the reduced export activity of the Rev protein, 95QSQGTET96, was identified. A total of 32 CSSs were found in 16 clades, wherein A6 had the 41Q substitution in the functionally significant region of Rev. The high values of WK coefficient in sites 51 and 82, located on the Rev interaction surface, indicate the susceptibility of these positions to evolutionary replacements. Thus, the noted signatures require further investigation.

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.

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.

A critical analysis of the current state of virus taxonomy

Taxonomical classification has preceded evolutionary understanding. For that reason, taxonomy has become a battleground fueled by knowledge gaps, technical limitations, and a priorism. Here we assess the current state of the challenging field, focusing on fallacies that are common in viral classification. We emphasize that viruses are crucial contributors to the genomic and functional makeup of holobionts, organismal communities that behave as units of biological organization. Consequently, viruses cannot be considered taxonomic units because they challenge crucial concepts of organismality and individuality. Instead, they should be considered processes that integrate virions and their hosts into life cycles. Viruses harbor phylogenetic signatures of genetic transfer that compromise monophyly and the validity of deep taxonomic ranks. A focus on building phylogenetic networks using alignment-free methodologies and molecular structure can help mitigate the impasse, at least in part. Finally, structural phylogenomic analysis challenges the polyphyletic scenario of multiple viral origins adopted by virus taxonomy, defeating a polyphyletic origin and supporting instead an ancient cellular origin of viruses. We therefore, prompt abandoning deep ranks and urgently reevaluating the validity of taxonomic units and principles of virus classification.

Recombinant characteristics of three novel HIV-1 second-generation recombinant forms composed of CRF01_AE and CRF07_BC isolated in Hebei province, China

Human immunodeficiency virus 1 (HIV-1) is a fast-evolving, genetically diverse virus. The HIV-1 evolution rate is also significantly influenced by the frequency of HIV-1 spread in a population. Transmission via homosexual contact has become the predominant transmission route, leading to an increase in the HIV-1 epidemic in Hebei province, China. In this study, we report three novel HIV-1 CRF01_AE/CRF07_BC recombinant forms isolated from three men who have sex with men (MSM) in the cities of Shijiazhuang (20747) and Langfang (20809 and 20820). Phylogenetic analysis based on HIV-1 near-full-length genome (NFLG) sequences indicated that the three novel recombinant forms formed a distinct monophyletic branch that was separate from all known HIV-1 subtypes and circulating recombinant forms (CRFs). Breakpoint analysis showed that the three NFLGs displayed different recombinant patterns. NFLGs 20747 and 20809 had a recombinant pattern with subtype CRF01_AE gene fragments inserted into a CRF07_BC backbone, spanning from the gag to env gene regions, whereas NFLG 20820 had a recombinant pattern with subtype CRF07_BC gene fragments inserted into a CRF01_AE backbone. Subregion phylogenetic analysis confirmed that these three NFLGs comprised CRF01_AE and CRF07_BC. Our findings confirm the emergence of novel recombinant forms and highlight the need for continuous monitoring of the diversity of HIV-1 among sexually active populations, especially MSM, to better control the HIV-1 epidemic.

HIV and Drug-Resistant Subtypes

Acquired Immunodeficiency Syndrome (AIDS) is a human viral infectious disease caused by the positive-sense single-stranded (ss) RNA Human Immunodeficiency Virus (HIV) (Retroviridae family, Ortervirales order). HIV-1 can be distinguished into various worldwide spread groups and subtypes. HIV-2 also causes human immunodeficiency, which develops slowly and tends to be less aggressive. HIV-2 only partially homologates to HIV-1 despite the similar derivation. Antiretroviral therapy (ART) is the treatment approved to control HIV infection, based on multiple antiretroviral drugs that belong to different classes: (i) NNRTIs, (ii) NRTIs, (iii) PIs, (iv) INSTIs, and (v) entry inhibitors. These drugs, acting on different stages of the HIV life cycle, decrease the patient's total burden of HIV, maintain the function of the immune system, and prevent opportunistic infections. The appearance of several strains resistant to these drugs, however, represents a problem today that needs to be addressed as best as we can. New outbreaks of strains show a widespread geographic distribution and a highly variable mortality rate, even affecting treated patients significantly. Therefore, novel treatment approaches should be explored. The present review discusses updated information on HIV-1- and HIV-2-resistant strains, including details on different mutations responsible for drug resistance.

Confidence intervals in molecular dating by maximum likelihood

Molecular dating has been widely used to infer the times of past evolutionary events using molecular sequences. This paper describes three bootstrap methods to infer confidence intervals under a penalized likelihood framework. The basic idea is to use data pseudoreplicates to infer uncertainty in the branch lengths of a phylogeny reconstructed with molecular sequences. The three specific bootstrap methods are nonparametric (direct tree bootstrapping), semiparametric (rate smoothing), and parametric (Poisson simulation). Our extensive simulation study showed that the three methods perform generally well under a simple strict clock model of molecular evolution; however, the results were less positive with data simulated using an uncorrelated or a correlated relaxed clock model. Several factors impacted, possibly in interaction, the performance of the confidence intervals. Increasing the number of calibration points had a positive effect, as well as increasing the sequence length or the number of sequences although both latter effects depended on the model of evolution. A case study is presented with a molecular phylogeny of the Felidae (Mammalia: Carnivora). A comparison was made with a Bayesian analysis: the results were very close in terms of confidence intervals and there was no marked tendency for an approach to produce younger or older bounds compared to the other.

Combining biomarker and virus phylogenetic models improves HIV-1 epidemiological source identification

To identify and stop active HIV transmission chains new epidemiological techniques are needed. Here, we describe the development of a multi-biomarker augmentation to phylogenetic inference of the underlying transmission history in a local population. HIV biomarkers are measurable biological quantities that have some relationship to the amount of time someone has been infected with HIV. To train our model, we used five biomarkers based on real data from serological assays, HIV sequence data, and target cell counts in longitudinally followed, untreated patients with known infection times. The biomarkers were modeled with a mixed effects framework to allow for patient specific variation and general trends, and fit to patient data using Markov Chain Monte Carlo (MCMC) methods. Subsequently, the density of the unobserved infection time conditional on observed biomarkers were obtained by integrating out the random effects from the model fit. This probabilistic information about infection times was incorporated into the likelihood function for the transmission history and phylogenetic tree reconstruction, informed by the HIV sequence data. To critically test our methodology, we developed a coalescent-based simulation framework that generates phylogenies and biomarkers given a specific or general transmission history. Testing on many epidemiological scenarios showed that biomarker augmented phylogenetics can reach 90% accuracy under idealized situations. Under realistic within-host HIV-1 evolution, involving substantial within-host diversification and frequent transmission of multiple lineages, the average accuracy was at about 50% in transmission clusters involving 5-50 hosts. Realistic biomarker data added on average 16 percentage points over using the phylogeny alone. Using more biomarkers improved the performance. Shorter temporal spacing between transmission events and increased transmission heterogeneity reduced reconstruction accuracy, but larger clusters were not harder to get right. More sequence data per infected host also improved accuracy. We show that the method is robust to incomplete sampling and that adding biomarkers improves reconstructions of real HIV-1 transmission histories. The technology presented here could allow for better prevention programs by providing data for locally informed and tailored strategies.