First description of two new HIV-1 recombinant forms CRF82_cpx and CRF83_cpx among drug users in Northern Myanmar

Optimizing ancestral trait reconstruction of large HIV Subtype C datasets through multiple-trait subsampling

Large datasets along with sampling bias represent a challenge for phylodynamic reconstructions, particularly when the study data are obtained from various heterogeneous sources and/or through convenience sampling. In this study, we evaluate the presence of unbalanced sampled distribution by collection date, location, and risk group of human immunodeficiency virus Type 1 Subtype C using a comprehensive subsampling strategy and assess their impact on the reconstruction of the viral spatial and risk group dynamics using phylogenetic comparative methods. Our study shows that a most suitable dataset for ancestral trait reconstruction can be obtained through subsampling by all available traits, particularly using multigene datasets. We also demonstrate that sampling bias is inflated when considerable information for a given trait is unavailable or of poor quality, as we observed for the trait risk group. In conclusion, we suggest that, even if traits are not well recorded, including them deliberately optimizes the representativeness of the original dataset rather than completely excluding them. Therefore, we advise the inclusion of as many traits as possible with the aid of subsampling approaches in order to optimize the dataset for phylodynamic analysis while reducing the computational burden. This will benefit research communities investigating the evolutionary and spatio-temporal patterns of infectious diseases.

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.

The Emergence and Transmission Dynamics of HIV-1 CRF07_BC in Mainland China

A total of 1155 partial pol gene sequences of human immunodeficiency virus (HIV)-1 CRF07_BC were sampled between 1997 and 2015, spanning 13 provinces in Mainland China and risk groups [heterosexual, injecting drug users (IDU), and men who have sex with men (MSM)] to investigate the evolution, adaptation, spatiotemporal and risk group dynamics, migration patterns, and protein structure of HIV-1 CRF07_BC. Due to the unequal distribution of sequences across time, location, and risk group in the complete dataset (‘full1155’), subsampling methods were used. Maximum-likelihood and Bayesian phylogenetic analysis as well as discrete trait analysis of geographical location and risk group were carried out. To study mutations of a cluster of HIV-1 CRF07_BC (CRF07-1), we performed a comparative analysis of this cluster to the other CRF07_BC sequences (‘backbone_295’) and mapped the mutations observed in the respective protein structure. Our findings showed that HIV-1 CRF07_BC most likely originated among IDU in Yunnan Province between October 1992 to July 1993 [95 per cent hightest posterior density (HPD): May 1989–August 1995] and that IDU in Yunnan Province and MSM in Guangdong Province likely served as the viral sources during the early and more recent spread in Mainland China. We also revealed that HIV-1 CRF07-1 has been spreading for roughly 20 years and continues to cause local transmission in Mainland China and worldwide. Overall, our study sheds light on the dynamics of HIV-1 CRF07_BC distribution patterns in Mainland China. Our research may also be useful in formulating public health policies aimed at controlling acquired immune deficiency syndrome in Mainland China and globally.

Genetic Characterization of a Novel HIV-1 Second-Generation Recombinant Form Originating from CRF86_BC and a Unique Recombinant Form in Yunnan, China

Yunnan is the first place where human immunodeficiency virus type 1 (HIV-1) became prevalent in China, and it is also the place with the most complicated HIV-1 genetic diversity in China. On October 23, 2019, a patient newly diagnosed with acquired immunodeficiency syndrome from a hospital in Baoshan, Yunnan, was recruited for genetic analysis. Near full-length genome of HIV-1 was amplified from the plasma sample. Phylogenetic analysis revealed that this sequence (BS6F24) has a close relationship with CRF86_BC and a unique recombinant form (URF) (KY406739), which was formed by recombination of subtypes B and C. Bootscan analysis confirmed that the first part (HXB2:1022-5832) and last part (HXB2:5833-9120) genomes of BS6F24 had the same recombinant structures as KY406739 and CRF86_BC, respectively. A second-generation recombinant form that originated from CRF86_BC and a URF were reported for the first time. This indicates the need for continuous monitoring of the genetic diversity of HIV-1 in Yunnan, China.

Phylogenetic Analysis of Sequences in the HIV Database Revealed Multiple Potential Circulating Recombinant Forms in China

HIV recombination contributes greatly to its diversity and produces many circulating recombinant forms (CRFs) and unique recombinant forms (URFs). In China, 24 CRFs have been reported to date, and CRFs cause more than 80% of HIV infections. However, the prevalence of CRFs might still be underestimated, as a high level of onward transmission of URFs has been reported. In this study, we analyzed all Chinese pol region (2,253-3,252) sequences in the HIV Database to evaluate potential new CRFs in China. HIV-1 genotypes were verified by the Context-based Modeling for Expeditious Typing (COMET) tool. Maximum-likelihood (ML) trees were constructed based on sequences with unassigned genotypes. Cluster Picker 1.2.1 was used to identify transmission clusters. Meanwhile, a jumping-profile hidden Markov model (jpHMM) was used to perform recombination breakpoint analysis. Beast 1.7.5 was used to estimate the time of the most recent common ancestor of new CRFs. In the HIV databases, CRF01_AE was the most prevalent genetic form in China, accounting for 39.69% of all national infections, followed by CRF07_BC (20.47%), subtype B (17.50%), CRF08_BC (6.60%), subtype C (6.28%), CRF55_01B (2.06%), and other CRFs (1.77%). The URFs were responsible for 5.31% of all infections nationwide. Among URFs, genomes comprising BC, 01BC, 01B, and 01C were dominant. Finally, 17 potential CRFs and 1 novel CRF were identified. BEAST analysis indicates that novel CRF originated around in 2009. The data highlight that more CRFs have been spreading in China. HIV-1 pol sequences that are commonly used to explore drug resistance are helpful for the surveillance of epidemics of different HIV-1 genotypes.

Identification of a newly emerging HIV-1 CRF111_01C comprising CRF01_AE and C in Yunnan, China: Genetic characterization and recombinant history

New human immunodeficiency virus (HIV)-1 circulating recombinant forms (CRFs) have recently emerged and disseminated rapidly in China; in total 38 CRFs have been identified thus far. Yunnan province shares its border with Myanmar, and is regarded as a "hotspot" for the occurrence of new HIV-1 recombinations; more than half of novel CRFs reported in China have been first documented in Yunnan province. In the present study, based on the information available on four existing near-full-length genome (NFLG) sequences, combined with data on four other closely related sequences obtained via Basic Local Alignment Search Tool (BLAST) analysis, NFLG/subregion phylogenetic, bootscanning, and time to the most recent common ancestor (TMRCA) analyses were performed. Phylogenetic analysis revealed that the eight strains demonstrated the formation of a distinct monophyletic branch with a bootstrap value of 100%. Strains in this branch were distantly related to all known HIV-1 CRFs; it was temporarily named CRF111_01C. Bootscanning analysis revealed that CRF111_01C consisted of a CRF01_AE backbone and four inserted subtype C segments. Remarkably, CRF111_01C shared six mosaic fragment identities with the previously identified CRF100_01C. Furthermore, CRF111_01C may be deemed a potential second-generation CRF consisting of CRF100_01C and C. Coalescent Bayesian analyses revealed that the TMRCA of CRF111_01C was approximately the period 1999-2002. The emergence of such second-generation recombinants highlights that continuous molecular screening is necessary to carefully monitor the evolutionary dynamics of HIV-1 epidemics.

Tree Shrew Cells Transduced with Human CD4 and CCR5 Support Early Steps of HIV-1 Replication, but Viral Infectivity Is Restricted by APOBEC3

The host range of human immunodeficiency virus type 1 (HIV-1) is narrow. Therefore, using ordinary animal models to study HIV-1 replication, pathogenesis, and therapy is impractical. The lack of applicable animal models for HIV-1 research spurred our investigation on whether tree shrews (Tupaia belangeri chinensis), which are susceptible to many types of human viruses, can act as an animal model for HIV-1. Here, we report that tree shrew primary cells are refractory to wild-type HIV-1 but support the early replication steps of HIV-1 pseudotyped with the vesicular stomatitis virus glycoprotein envelope (VSV-G), which can bypass entry receptors. The exogenous expression of human CD4 renders the tree shrew cell line infectible to X4-tropic HIV-1IIIB, suggesting that tree shrew CXCR4 is a functional HIV-1 coreceptor. However, tree shrew cells did not produce infectious HIV-1 progeny virions, even with the human CD4 receptor. Subsequently, we identified tree shrew (ts) apolipoprotein B editing catalytic polypeptide 3 (tsAPOBEC3) proteins as active inhibitors of HIV-1 particle infectivity, with virus infectivity reduced 10- to 1,000-fold. Unlike human APOBEC3G, the tsA3Z2c-Z1b protein was not degraded by the HIV-1 viral infectivity factor (Vif) but markedly restricted HIV-1 replication through mutagenicity and reverse transcription inhibition. The pooled knockout of tsA3Z2c-Z1b partially restored the infectivity of the HIV-1 progeny. This work suggests that tsAPOBEC3 proteins serve as an additional barrier to the development of HIV-1 tree shrew models, even when virus entry is overcome by exogenous expression of human CD4. IMPORTANCE The development of animal models is critical for studying human diseases and their pathogenesis and for evaluating drug and vaccine efficacy. For improved AIDS research, the ideal animal model of HIV-1 infection should be a small laboratory mammal that closely mimics virus replication in humans. Tree shrews exhibit considerable potential as animal models for the study of human diseases and therapeutic responses. Here, we report that human CD4-expressing tree shrew cells support the early steps of HIV-1 replication and that tree shrew CXCR4 is a functional coreceptor of HIV-1. However, tree shrew cells harbor additional restrictions that lead to the production of HIV-1 virions with low infectivity. Thus, the tsAPOBEC3 proteins are partial barriers to developing tree shrews as an HIV-1 model. Our results provide insight into the genetic basis of HIV inhibition in tree shrews and build a foundation for the establishment of gene-edited tree shrew HIV-1-infected models.

Genetic Characteristics of Three Unique Recombinant Forms of HIV-1 in Yunnan, China

To explore the molecular epidemiological status of human immunodeficiency virus type 1 (HIV-1) in Yunnan, China, three HIV-1 near full-length genomes were amplified and sequenced from plasma samples that were collected from Burmese patients newly diagnosed with HIV-1 in Dehong Prefecture in Yunnan Province in 2017. Phylogenetic and bootscanning analyses revealed that all the sequences might be HIV-1 second-generation recombinant forms of circulating recombinant forms (CRF07_BC and CRF83_cpx) and unique recombinant forms. One of the sequences contained six CRF01_AE fragments, five subtype C fragments, and two subtype B fragments, which were separated by 12 breakpoints. These results revealed that the second-generation recombination of HIV-1 within different strains is still ongoing in Dehong, China. Systematic surveys and immediate interventions are urgently needed to prevent the formation of increasingly complex HIV-1 recombinant forms.

The Genetic Diversity of HIV-1 Within Antiretroviral-Naive Outpatients in Ganzhou, China

To explore the molecular epidemiological status of human immunodeficiency virus type 1 (HIV-1) in Ganzhou, China, eight HIV-1-positive outpatients were recruited from July 5 to 21, 2018. Six HIV-1 near-full-length sequences were amplified and sequenced from the plasma samples that were collected before the patients' antiretroviral treatments. Phylogenetic and bootscan analyses revealed that one of the sequences was CRF01_AE, one was CRF55_01B, and two were CRF07_BC. Notably, one of the sequences was a unique recombinant form, and one of them was a second-generation recombinant form of CRF07_BC and subtype C. These results revealed that multiple HIV-1 subtypes are circulating in Ganzhou, China. Systematic surveys with large sample sizes are urgently needed to explore the exact molecular epidemiological characteristics and to trace the origins of HIV-1 in Ganzhou, China.

HIV-1 gp41 genetic diversity and enfuvirtide resistance-associated mutations among enfuvirtide-naïve patients in southern China

BACKGROUND

Human immunodeficiency virus (HIV) increasing molecular diversity and emergence of drug resistant mutants remain a major concern in China. Enfuvirtide (ENF/T-20) is the first entry inhibitor used in patients failing highly active antiretroviral therapy (HAART). However, data on HIV-1 gp41genetic diversity and primary ENF resistance-associated mutations among treatment-naïve patients in China is limited. The objective was to identify molecular diversity and ENF resistance patterns of HIV-1 in southern China, using envelope (env) gp41 sequences and bioinformatics tools, which may help optimize antiretroviral therapy.

METHODS

From November 2018 to January 2019, 439 blood plasma samples from ENF-naïve patients were collected from Shenzhen, Wuhan and Chongqing, of which 396 HIV env regions were sequenced and subtyped, and were performed the analysis of drug resistance-associated mutations (DRMs).

RESULTS

The main subtypes were circulating recombinant form (CRF) 01_AE (30.6 %) and CRF07_BC (48.7 %). CRF55_01B had been the fourth subtype in the study, and many rare CRFs were observed. Notably, CRF02_AG and CRF_BF strains typically found in Africa and US respectively were identified amongst Chinese patients. Known DRMs were detected in 27.5 % (109/396) of ENF treatment-naïve patients. One major DRM (L44 M), many secondary DRMs (including N126 K, E137 K, S138A), and lots of polymorphisms were found in the study, which have been proved to elevate resistance to ENF.

CONCLUSIONS

HIV-1 molecular diversity was observed in the study, which indicating that HIV-1 variability is becoming increasingly complex in southern China. A diverse set of primary DRMs discovered in this study described the serious threat to ART, which reminds us the urgent need of timely surveillance of HIV-1 viral diversity and drug resistance in China.

HIV-1 Drug Resistance in ART-Naïve Individuals in Myanmar

Background

Estimating the prevalence and characterizing the transmission of HIV-1 drug resistance in treatment-naïve individuals are very important in the prevention and control of HIV/AIDS. As one of the areas most affected by HIV/AIDS, few data are currently available for HIV-1 drug resistance in antiretroviral therapy (ART)-naïve individuals in Myanmar, which borders Yunnan, China.

Methods

HIV-1 pol sequences from ART-naïve HIV-1-infected individuals during 2008 and 2014 in Myanmar were retrieved from our previous studies. HIV-1 transmitted drug resistance (TDR) and susceptibility to antiretroviral drugs were predicted using the Stanford HIVdb program. HIV-1 transmission cluster (TC) was determined by Cluster Picker.

Results

A total of 169 partial pol sequences from ART-naïve HIV-1 positive Burmese were analyzed. The prevalence of TDR was 20.1%. CRF01_AE and BC recombinants appeared to have a higher prevalence of TDR than other subtypes. The V179D/T was found to be very common in the China–Myanmar border region and was involved in half of the transmission clusters formed by HIV-1 drug-resistance strains in this region. Comparison showed that drug-resistance mutation profile in Myanmar was very similar to that in Dehong prefecture of Yunnan. By further phylogenetic analysis with all available sequences from the China–Myanmar border region, four HIV-1 drug-resistance-related TCs were identified. Three of them were formed by Burmese long-distance truck drivers and the Burmese staying in Yunnan, and another was formed by Burmese injection drug users staying in Myanmar and Yunnan. These results suggest a potential transmission link of HIV-1 drug resistance between Myanmar and Yunnan.

Conclusion

Given the high prevalence of TDR in Myanmar, and the potential risk of cross-border transmission of HIV-1 drug-resistant strains between Myanmar and Yunnan, China, ongoing monitoring of HIV-1 drug resistance in ART-naïve individuals will provide a guideline for clinical antiretroviral treatment and benefit the prevention and control of HIV/AIDS in this border region.

Laos is Affected by HIV CRF01_AE and the Newly Identified CRF97_01B

Laos is the only landlocked country in Southeast Asia and borders Thailand, Myanmar and Cambodia, the three countries in this region that have been hardest hit by human immunodeficiency virus (HIV). Laos has been regarded as a low-HIV-prevalence country for decades. To understand the status of HIV in Laos in recent years, a retrospective study was performed among 2851 patients visiting a hospital in Vientiane, the capital of Laos, from November 2011 to May 2012. Whole blood samples were obtained from the patients, and DNA was extracted. HIV status was determined by HIV gag fragment-specific PCR assay. Sixty-nine samples were detected as HIV proviral DNA positive with a positive rate of 2.4% (69/2851). Sixty-one near full-length genomic sequences were obtained from the positive samples. The results of phylogenetic analysis showed that the vast majority (91.8%) of the HIV strains belonged to CRF01_AE, and the other five (8.2%) strains were identified as a new HIV circulating recombinant form CRF97_01B, which had a CRF01_AE backbone with an insertion of subtype B in the gag-pol region. Phylogeographic analysis revealed that HIV CRF01_AE circulating in Laos were multiply introduced from Thailand. These results indicated that Laos might be suffering a considerably more serious impact of HIV than previously believed. To keep this country from undergoing the same increase in HIV prevalence observed in its neighbors, immediate intervention measures and sufficient epidemiological research are urgently needed.

Spatial clusters of HIV-1 genotypes in a recently infected population in Yunnan, China

Background

As a gateway for HIV-1 in China, Yunnan has experienced dramatic changes in HIV-1 epidemics, during which HIV-1 genotypes have become complex. To track dynamic changes in HIV-1 genotypes, an HIV-1 molecular epidemiological study was implemented in the recently infected population in Yunnan.

Methods

From 6,357 HIV-1-positive samples diagnosed during the first half of 2015 in Yunnan, 586 samples were identified as recent infections with BED-capture enzyme immunoassay (CEIA) and were subjected to phylogenetic analyses. Spatial scanning analyses for the main HIV-1 genotypes were also performed.

Results

Among the 439 specimens successfully genotyped, more than ten genotypes were detected, including CRF08_BC (45.3%), CRF07_BC (19.4%), unique recombinant forms (URFs) (18.2%), CRF01_AE (11.4%), subtype C (2.1%), CRF85_BC (1.1%), CRF55_01B (0.9%), subtype B (0.5%), CRF64_BC (0.5%), CRF59_01B (0.2%), CRF83_cpx (0.2%) and CRF87_cpx (0.2%). Females, Chinese, heterosexual contact and intravenous drug injection were significantly associated with CRF08_BC infection; homosexual contact was significantly associated with CRF01_AE and CRF07_BC infection; males and non-Chinese had a higher risk of URF infection than females. Among all HIV-1 genotypes, the geographic coverage of CRF08_BC was the largest. For CRF08_BC, CRF07_BC, URFs and CRF01_AE, spatial clusters were detected. The two CRF08_BC clusters and one URF cluster were associated with heterosexual transmission, and two of CRF01_AE clusters were associated with homosexual transmission. Transmitted drug resistance (TDR)-associated mutations were detected in 2.4% of individuals.

Conclusions

The diversity of HIV-1 genotypes increased in recent infections because of a long-term HIV-1 epidemic in Yunnan. The predominant HIV-1 strains showed distinct demographic characteristics and formed spatial clusters. These findings improved our understanding of the evolution of HIV-1 in Yunnan and provided information for further HIV-1 control and prevention.

Electronic supplementary material

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

Burmese injecting drug users in Yunnan play a pivotal role in the cross-border transmission of HIV-1 in the China-Myanmar border region

Injecting drug users (IDUs) are the major risk group for HIV-1 infection in the China-Myanmar border area. There are a large number of Burmese IDUs living in Yunnan (Yunnan-mIDUs) who might be associated with the cross-border transmission of HIV-1. From 2010 to 2013, 617 Yunnan-mIDUs were recruited from three counties of Yunnan, 19.0% of whom were detected to be HIV-1 positive by serological testing. Partial HIV-1 p17, pol, vif-env, and env genes were amplified from the positive samples and were sequenced. Phylogenetic and HIV-1 subtyping analyses revealed that HIV-1 recombinant forms (RFs), including RF_BC (36.4%), RF_01BC (26.1%), RF_01C (9.1%) and RF_01B (1.1%), were predominant among this cohort. Of the identified HIV-1 strains, 14.8%, 9.1% and 3.4% belonged to subtype C, CRF01_AE and subtype B, respectively. Transmission cluster analysis showed that sequences from the Yunnan-mIDUs formed transmission clusters not only with those from Burmese IDUs but also with those from Chinese IDUs, indicating that Yunnan-mIDUs might acquire HIV-1 infection from or spread HIV-1 to both Burmese and Chinese IDUs. Phylogeographic analyses revealed three cross-border transmission patterns associated with Yunnan-mIDUs, in which Yunnan-mIDUs served as the crucial nodes linking the Burmese and Chinese IDUs. These results suggest that Yunnan-mIDUs are a potential viral reservoir for the diffusion of HIV-1 in Yunnan and play a pivotal role in the bidirectional cross-border transmission of HIV-1 in the China-Myanmar border region. More intervention efforts that focus on Yunnan-mIDUs are recommended in Yunnan’s campaign against HIV/AIDS.

Characterization of a Novel HIV-1 Circulating Recombinant Form, CRF01_AE/B'/C (CRF96_cpx), in Yunnan, China

Currently, complex HIV-1 recombinations among the B', C, and CRF01_AE genotypes have frequently arisen in Yunnan, China. A novel HIV-1 complex circulating recombinant form (CRF) consisting of B', C, and CRF01_AE (CRF96_cpx) was recently characterized from three epidemiologically unlinked individuals. Two strains of them were isolated from the injecting drug users in this study, the remaining one strain (JL. RL01) was obtained from a previous report in Jilin province. Phylogenetic analysis based on near full-length genome revealed that CRF96_cpx formed a distinct monophyletic cluster supported by a high bootstrap value of 100%, distantly related to all known HIV-1 subtypes/CRFs. CRF96_cpx had a CRF01_AE backbone with three subtype B' and C segments inserted, respectively, in the gag and pol region. Furthermore, subregion tree analysis showed that CRF01_AE backbone and subtype B segment inserted originated from a Thai-CRF01_AE lineage, whereas subtype C fragment inserted was from an India C lineage. They are different from previously documented CRF01_AE/B/C forms in its distinct backbone, inserted fragment size, and breakpoints. This highlighted the importance of continual monitoring of genetic diversity and complexity of HIV-1 strains in Yunnan, China.

Identification of Two New HIV-1 Circulating Recombinant Forms (CRF87_cpx and CRF88_BC) from Reported Unique Recombinant Forms in Asia

The on-going generation of HIV-1 intersubtype recombination has led to new circulating recombinant forms (CRFs) and unique recombinant forms (URFs) in Asia. In this study, we evaluated whether previously reported URFs were actually CRFs. All available complete or near full-length HIV-1 URF sequences from Asia were retrieved from the HIV Los Alamos National Laboratory Sequence database, and phylogenetic, transmission cluster, and bootscan analyses were performed using MEGA 6.0, Cluster Picker 1.2.1, and SimPlot3.5.1. According to the criterion of new CRFs, two new HIV-1 CRFs (CRF87_cpx and CRF88_BC) were identified from these available URFs. CRF87_cpx comprised HIV-1 subtypes B, C, and CRF01_AE, and CRF88_BC comprised subtypes B and C. HIV Blast and bootscan analysis revealed that besides the three representative strains, there were two additional CRF87_cpx strains. Furthermore, we defined seven dominant URFs (dURF01-dURF07), each of which contained two strains sharing same recombination map and can be used as sequence references to facilitate the finding of new potential CRFs in future. These results will benefit the molecular epidemiological investigation of HIV-1 in Asia.