Spatiotemporal Patterns of CRF07_BC in China: A Population-Based Study of the HIV Strain With the Highest Infection Rates

The distinct spatiotemporal evolutionary landscape of HBV and HDV largely determines the unique epidemic features of HDV globally

Chronic infection of hepatitis B virus (HBV) and hepatitis D virus (HDV) causes the most severe form of viral hepatitis. Due to the dependence on HBV, HDV was deemed to co-evolve and co-migrate with HBV. However, we previously found that the naturally occurred HDV/HBV combinations do not always reflect the most efficient virological adaptation (Wang et al., 2021). Moreover, regions with heavy HBV burden do not always correlate with high HDV prevalence (e.g., East Asia), and vice versa (e.g., Central Asia). Herein, we systematically elucidated the spatiotemporal evolutionary landscape of HDV to understand the unique epidemic features of HDV. We found that the MRCA of HDV was from South America around the late 13th century, was globally dispersed mainly via Central Asia, and evolved into eight genotypes from the 19th to 20th century. In contrast, the MRCA of HBV was from Europe ∼23.7 thousand years ago (Kya), globally dispersed mainly via Africa and East Asia, and evolved into eight genotypes ∼1100 years ago. When HDV stepped in, all present-day HBV genotypes had already formed and its global genotypic distribution had stayed stable geographically. Nevertheless, regionalized HDV adapted to local HBV genotypes and human lineages, contributing to the global geographical separation of HDV genotypes. Additionally, a sharp increase in HDV infections was observed after the 20th century. In conclusion, HDV exhibited a distinct spatiotemporal distribution path compared with HBV. This unique evolutionary relationship largely fostered the unique epidemic features we observe nowadays. Moreover, HDV infections may continue to ramp up globally, thus more efforts are urgently needed to combat this disease.

Molecular network characteristics and drug resistance analysis of 392 newly reported MSM HIV/AIDS cases in Chongqing, China

To comprehensively investigate the molecular transmission patterns of HIV-1 genotypes among men who have sex with men (MSM) in Chongqing, we employed 392 pol sequences of MSM to construct a phylogenetic tree and gene transmission network. Among the viral subtypes, CRF07_BC accounted for 73.2% (287/392) and CRF01_AE accounted for 20.7% (81/392), emerging as the predominant subtypes in this investigation. Additionally, we observed the presence of CRF55_01B, subtype B, CRF08_BC and other circulating recombinant forms. The HIV-1 molecular network was constructed with a gene distance threshold of 1.5%, resulting in an entry rate of 61.4% (241/392). Within the network, we identified a total of 23 molecular clusters, with the largest cluster being the CRF07_BC molecular cluster comprising 148 node values. Transmitted drug-resistance (TDR) mutations were found in 4.34% of the cases, with 1.79% associated with protease inhibitors (PIs), 0.51% with nucleoside reverse transcriptase inhibitors (NRTIs), and 2.55% with non-nucleoside reverse transcriptase inhibitors (NNRTIs). Statistical analysis indicated a higher enrollment rate in the HIV-1 molecular network among infected individuals with the CRF07_BC subtype, those identifying with same-sex sexual roles as "vers," and individuals with higher education levels. This suggests the need for strengthened investigation and intervention in this population to prevent the formation of larger transmission clusters. Furthermore, continuous monitoring of the HIV-1 molecular dynamics network is necessary to promptly and accurately track changes in molecular epidemic characteristics.

Near-full-length genome analysis of two novel HIV second recombinant forms in Hebei, China

The number of individuals infected with HIV-1 among men who have sex with men (MSM) has risen rapidly in recent years in China, and the subtypes CRF01_AE, CRF07_BC, and B, as well as many novel unique recombinant forms (URFs) are prevalent among them. Co-circulation of strains among MSM populations allows the generation of circulating recombinant forms (CRFs) and URFs. In this study, we identified two new URFs from two HIV-1-positive subjects who were infected through homosexual contact in Hebei, China. Analysis of near-full-length genome sequences, using phylogenetic and recombination analysis showed that the two URFs originated from CRF01_AE, CRF07_BC, and B, and CRF01_AE segments in the backbone of the URFs were derived from cluster 4 of CRF01_AE. The CRF07_BC segments of two URFs were clustered with 07BC_N in a phylogenetic tree. The identification of novel URFs with complex genomic structures shows that it is necessary to strengthen surveillance of HIV-1 variants in MSM populations in this region.

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.

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.

Genetic Diversity and Characteristics of Drug Resistance Among Treatment-Naive People Living with HIV in Xi'an, China

Purpose

The genetic diversity and genetic predisposition for drug resistance mutations are the primary features of human immunodeficiency virus type 1 (HIV-1), which could cause the incidence of failure of antiretroviral therapy (ART). This study investigates the distribution of various HIV-1 genotypes and the incidence of pretreatment drug resistance (PDR) in the antiretroviral-naive HIV-1 infected participants in Xi'an, China.

Patients and Methods

In this study, a cross-sectional analysis was carried out at the Xi'an Eighth Hospital between January 2020 and December 2021 among newly-diagnosed ART-naive HIV-1 infected participants. A nested PCR technique was used for amplifying the target segment of 1.3 kb present in the pol gene that spanned the reverse transcriptase and the protease regions. HIV-1 genotypes and the PDR-associated mutations were identified using the Stanford HIV Drug Resistance Database.

Results

A total of 317 pol gene sequences were retrieved, amplified, and sequenced. The circulating recombinant form (CRF), CRF07_BC (51.7%) was seen to be the most prevalent HIV-1 genotype, followed by other genotypes like CRF01_AE (25.9%), B (14.2%), and CRF55_01B (4.7%). PDR was found in 18.3% of the population. The PDR mutation frequency in the non-nucleoside reverse transcriptase inhibitor (NNRTI) (16.1%) was significantly higher compared to that of the nucleoside reverse transcriptase inhibitor (NRTI) (4.4%) and the protease inhibitor (0.9%). V179D/E (both were 4.4%) was seen to be the most predominant type of NNRTI mutation. K65R and M184V (1.3%) were the most frequent NRTI-associated mutations. About half (48.3%) of the sequenced HIV-1 strains that had mutations could show a potential low-level NNRTI resistance owing to V179D/E. Multivariate regression analysis revealed one PDR mutation associated with subtype CRF01_AE (p=0.002) and CRF55_01B (p<0.001) as a higher risk mutation.

Conclusion

Diverse and complex HIV-1 genotypes are distributed in Xi'an, China. Considering new evidence, it is necessary to screen for baseline HIV-1 drug resistance among the newly-diagnosed HIV-1 individuals.

A New HIV-1 K28E32-Reverse Transcriptase Variant Associated with the Rapid Expansion of CRF07_BC among Men Who Have Sex with Men

HIV-1 CRF07_BC originated among injection drug users (IDUs) in China. After diffusing into men who have sex with men (MSM), CRF07_BC has shown a rapid expansion in this group; however, the mechanism remains unclear. Here, we identified a new K₂₈E₃₂ variant of CRF07_BC that was characterized by five specific mutations (E28K, K32E, E248V, K249Q, and T338S) in reverse transcriptase. This variant was mainly prevalent among MSM, and was overrepresented in transmission clusters, suggesting that it could have driven the rapid expansion of CRF07_BC in MSM, though founder effects cannot be ruled out. It was descended from an evolutionary intermediate accumulating four specific mutations and formed an independent phylogenetic node with an estimated origin time in 2003. The K₂₈E₃₂ variant was demonstrated to have significantly higher in vitro HIV-1 replication ability than the wild type. Mutations E28K and K32E play a critical role in the improvement of in vitro HIV-1 replication ability, reflected by improved reverse transcription activity. The results could allow public health officials to use this marker (especially E28K and K32E mutations in the reverse transcriptase (RT) coding region) to target prevention measures prioritizing MSM population and persons infected with this variant for test and treat initiatives. IMPORTANCE HIV-1 has very high mutation rate that is correlated with the survival and adaption of the virus. The variants with higher transmissibility may be more selective advantage than the strains with higher virulence. Several HIV-1 variants were previously demonstrated to be correlated with higher viral load and lower CD4 T cell count. Here, we first identified a new variant (the K₂₈E₃₂ variant) of HIV-1 CRF07_BC, described its origin and evolutionary dynamics, and demonstrated its higher in vitro HIV-1 replication ability than the wild type. We demonstrated that five RT mutations (especially E28K and K32E) significantly improve in vitro HIV-1 replication ability. The appearance of the new K₂₈E₃₂ variant was associated with the rapidly increasing prevalence of CRF07_BC among MSM.