Low-level viremia episodes appear to affect the provirus composition of the circulating cellular HIV reservoir during antiretroviral therapy.

Low-level viremia (LLV) ranging from 50 to 1,000 copies/ml is common in most HIV-1-infected patients receiving antiretroviral therapy (ART). However, the source of LLV and the impact of LLV on the HIV-1 reservoir during ART remain uncertain. We hypothesized that LLV may arise from the HIV reservoir and its occurrence affect the composition of the reservoir after LLV episodes. Accordingly, we investigated the genetic linkage of sequences obtained from plasma at LLV and pre-ART time points and from peripheral blood mononuclear cells (PBMCs) at pre-ART, pre-LLV, LLV, and post-LLV time points. We found that LLV sequences were populated with a predominant viral quasispecies that accounted for 67.29%∼100% of all sequences. Two episodes of LLV in subject 1, spaced 6 months apart, appeared to have originated from the stochastic reactivation of latently HIV-1-infected cells. Moreover, 3.77% of pre-ART plasma sequences were identical to 67.29% of LLV-3 plasma sequences in subject 1, suggesting that LLV may have arisen from a subset of cells that were infected before ART was initiated. No direct evidence of sequence linkage was found between LLV viruses and circulating cellular reservoirs in all subjects. The reservoir size, diversity, and divergence of the PBMC DNA did not differ significantly between the pre- and post-LLV sampling points (P > 0.05), but the composition of viral reservoir quasispecies shifted markedly before and after LLV episodes. Indeed, subjects with LLV had a higher total PBMC DNA level, greater viral diversity, a lower proportion of variants with identical sequences detected at two or more time points, and a shorter variant duration during ART compared with subjects without LLV. Overall, our findings suggested that LLV viruses may stem from an unidentified source other than circulating cellular reservoirs. LLV episodes may introduce great complexity into the HIV reservoir, which brings challenges to the development of treatment strategies.

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.

Comparative miRNA transcriptomics of macaques and mice reveals MYOC is an inhibitor for Cryptococcus neoformans invasion into the brain.

Cryptococcal meningoencephalitis (CM) is emerging as an infection in HIV/AIDS patients shifted from primarily ART-naive to ART-experienced individuals, as well as patients with COVID-19 and immunocompetent hosts. This fungal infection is mainly caused by the opportunistic human pathogen Cryptococcus neoformans. Brain or central nervous system (CNS) dissemination is the deadliest process for this disease; however, mechanisms underlying this process have yet to be elucidated. Moreover, illustrations of clinically relevant responses in cryptococcosis are currently limited due to the low availability of clinical samples. In this study, to explore the clinically relevant responses during C. neoformans infection, macaque and mouse infection models were employed and miRNA-mRNA transcriptomes were performed and combined, which revealed cytoskeleton, a major feature of HIV/AIDS patients, was a centric pathway regulated in both infection models. Notably, assays of clinical immune cells confirmed an enhanced macrophage "Trojan Horse" in patients with HIV/AIDS, which could be shut down by cytoskeleton inhibitors. Furthermore, myocilin, encoded by MYOC, was found to be a novel enhancer for the macrophage "Trojan Horse," and an enhanced fungal burden was achieved in the brains of MYOC-transgenic mice. Taken together, the findings from this study reveal fundamental roles of the cytoskeleton and MYOC in fungal CNS dissemination, which not only helps to understand the high prevalence of CM in HIV/AIDS but also facilitates the development of novel therapeutics for meningoencephalitis caused by C. neoformans and other pathogenic microorganisms.

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.

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 (R0) of two lineages were similar and both of the effective production number (Re) fell below 1 at the most recent sampling time, the starts of the Re 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.

Prevention and control of HIV/AIDS in China: lessons from the past three decades.

ABSTRACT: In the past 37 years, human immunodeficiency virus/acquired immunodeficiency syndrome (HIV/AIDS) has undergone various major transmission routes in China, with the world most complex co-circulating HIV-1 subtypes, even the prevalence is still low. In response to the first epidemic outbreak of HIV in injecting drug users and the second one by illegal commercial blood collection, China issued the Anti-Drug Law and launched the Blood Donation Act and nationwide nucleic acid testing, which has avoided 98,232 to 211,200 estimated infections and almost ended the blood product-related infection. China has been providing free antiretroviral therapy (ART) since 2003, which covered >80% of the identified patients and achieved a viral suppression rate of 91%. To bend the curve of increasing the disease burden of HIV and finally end the epidemic, China should consider constraining HIV spread through sexual transmission, narrowing the gaps in identifying HIV cases, and the long-term effectiveness and safety of ART in the future.

Transmitted drug resistance to Tenofovir/Emtricitabine among persons with newly diagnosed HIV infection in Shenyang city, Northeast China from 2016 to 2018.

BACKGROUND: To assess transmitted drug resistance (TDR) to tenofovir (TDF)/emtricitabine (FTC), using as pre-exposure prophylaxis, among newly diagnosed human immunodeficiency virus-1 (HIV-1)-infected residents in Shenyang city, northeast China. METHODS: Demographic and epidemiological information of all newly diagnosed HIV-1 infected residents in Shenyang city from 2016 to 2018 were anonymously collected from the local HIV epidemic database. HIV-1 pol sequences were amplified from RNA in cryopreserved plasma samples and sequenced directly. Viral subtypes were inferred with phylogenetic analysis and drug resistance mutations (DRMs) were determined according to the Stanford HIVdb algorithm. Recent HIV infection was determined with HIV Limiting Antigen avidity electro immunoassay. RESULTS: A total of 2176 sequences (92.4%, 2176/2354) were obtained; 70.9% (1536/2167) were CRF01_AE, followed by CRF07_BC (18.0%, 391/2167), subtype B (4.7%, 102/2167), other subtypes (2.6%, 56/2167), and unique recombinant forms (3.8%, 82/2167). The prevalence of TDR was 4.9% (107/2167), among which, only 0.6% (13/2167) was resistance to TDF/FTC. Most of these subjects had CRF01_AE strains (76.9%, 10/13), were unmarried (76.9%, 10/13), infected through homosexual contact (92.3%, 12/13), and over 30 years old (median age: 33). The TDF/FTC DRMs included K65R (8/13), M184I/V (5/13), and Y115F (2/13). Recent HIV infection accounted for only 23.1% (3/13). Most cases were sporadic in the phylogenetic tree, except two CRF01_AE sequences with K65R (Bootstrap value: 99%). CONCLUSIONS: The prevalence of TDR to TDF/FTC is low among newly diagnosed HIV-infected cases in Shenyang, suggesting that TDR may have little impact on the protective effect of the ongoing CROPrEP project in Shenyang city.

Priority Intervention Targets Identified Using an In-Depth Sampling HIV Molecular Network in a Non-Subtype B Epidemics Area.

Molecular network analysis based on the genetic similarity of HIV-1 is increasingly used to guide targeted interventions. Nevertheless, there is a lack of experience regarding molecular network inferences and targeted interventions in combination with epidemiological information in areas with diverse epidemic strains of HIV-1.We collected 2,173 pol sequences covering 84% of the total newly diagnosed HIV-1 infections in Shenyang city, Northeast China, between 2016 and 2018. Molecular networks were constructed using the optimized genetic distance threshold for main subtypes obtained using sensitivity analysis of plausible threshold ranges. The transmission rates (TR) of each large cluster were assessed using Bayesian analyses. Molecular clusters with the characteristics of ≥5 newly diagnosed cases in 2018, high TR, injection drug users (IDUs), and transmitted drug resistance (TDR) were defined as priority clusters. Several HIV-1 subtypes were identified, with a predominance of CRF01_AE (71.0%, 1,542/2,173), followed by CRF07_BC (18.1%, 393/2,173), subtype B (4.5%, 97/2,173), other subtypes (2.6%, 56/2,173), and unique recombinant forms (3.9%, 85/2,173). The overall optimal genetic distance thresholds for CRF01_AE and CRF07_BC were both 0.007 subs/site. For subtype B, it was 0.013 subs/site. 861 (42.4%) sequences of the top three subtypes formed 239 clusters (size: 2-77 sequences), including eight large clusters (size ≥10 sequences). All the eight large clusters had higher TR (median TR = 52.4/100 person-years) than that of the general HIV infections in Shenyang (10.9/100 person-years). A total of ten clusters including 231 individuals were determined as priority clusters for targeted intervention, including eight large clusters (five clusters with≥5 newly diagnosed cases in 2018, one cluster with IDUs, and two clusters with TDR (K103N, Q58E/V179D), one cluster with≥5 newly diagnosed cases in 2018, and one IDUs cluster. In conclusion, a comprehensive analysis combining in-depth sampling HIV-1 molecular networks construction using subtype-specific optimal genetic distance thresholds, and baseline epidemiological information can help to identify the targets of priority intervention in an area epidemic for non-subtype B.

Integrated analysis of lncRNA, miRNA and mRNA profiles reveals potential lncRNA functions during early HIV infection.

BACKGROUND: Long noncoding RNAs (lncRNAs) can regulate gene expression in a cis-regulatory fashion or as "microRNA sponges". However, the expression and functions of lncRNAs during early human immunodeficiency virus (HIV) infection (EHI) remain unclear. METHODS: 3 HAART-naive EHI patients and 3 healthy controls (HCs) were recruited in this study to perform RNA sequencing and microRNA (miRNA) sequencing. The expression profiles of lncRNAs, mRNAs and miRNAs were obtained, and the potential roles of lncRNAs were analysed based on discovering lncRNA cis-regulatory target mRNAs and constructing lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) networks. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on 175 lncRNA-associated differentially expressed (DE) mRNAs to investigate the potential functions of DE lncRNAs in ceRNA networks. RESULTS: A total of 242 lncRNAs, 1240 mRNAs and 21 mature known miRNAs were determined as differentially expressed genes in HAART-naive EHI patients compared to HCs. Among DE lncRNAs, 44 lncRNAs were predicted to overlap with 41 target mRNAs, and 107 lncRNAs might regulate their nearby DE mRNAs. Two DE lncRNAs might regulate their cis-regulatory target mRNAs BTLA and ZAP70, respectively, which were associated with immune activation. In addition, the ceRNA networks comprised 160 DE lncRNAs, 21 DE miRNAs and 175 DE mRNAs. Seventeen DE lncRNAs were predicted to regulate HIF1A and TCF7L2, which are involved in the process of HIV-1 replication. Twenty DE lncRNAs might share miRNA response elements (MREs) with FOS, FOSB and JUN, which are associated with both immune activation and HIV-1 replication. CONCLUSIONS: This study revealed that lncRNAs might play a critical role in HIV-1 replication and immune activation during EHI. These novel findings are helpful for understanding of the pathogenesis of HIV infection and provide new insights into antiviral therapy.

First complete-genome documentation of HIV-1 intersubtype superinfection with transmissions of diverse recombinants over time to five recipients.

Human immunodeficiency virus type 1 (HIV-1) recombinants in the world are believed to be generated through recombination between distinct HIV-1 strains among coinfection or superinfection cases. However, direct evidence to support transmission of HIV-1 recombinants from a coinfected/superinfected donor to putative recipient is lacking. Here, we report on the origin and evolutionary relationship between a set of recombinants from a CRF01_AE/CRF07_BC superinfected putative donor and diverse CRF01_AE/CRF07_BC recombinants from five putative recipients. Interviews on sociodemographic characteristics and sexual behaviors for these six HIV-1-infected men who have sex with men showed that they had similar ways of partner seeking: online dating sites and social circles. Phylogenetic and recombination analyses demonstrated that the near-full-length genome sequences from six patients formed a monophyletic cluster different from known HIV-1 genotypes in maximum likelihood phylogenetic trees, were all composed of CRF01_AE and CRF07_BC fragments with two common breakpoints on env, and shared 4-7 breakpoints with each other. Moreover, 3' half-genomes of recombinant strains from five recipients had identical/similar recombinant structures with strains at longitudinal samples from the superinfected donor. Recombinants from the donor were paraphyletic, whereas five recipients were monophyletic or polyphyletic in the maximum clade credibility tree. Bayesian analyses confirmed that the estimated time to the most recent common ancestor (tMRCA) of CRF01_AE and CRF07_BC strains of the donor was 2009.2 and 2010.7, respectively, and all were earlier than the emergence of recombinants from five recipients. Our results demonstrated that the closely related unique recombinant forms of HIV-1 might be the descendent of a series of recombinants generated gradually in a superinfected patient. This finding highlights the importance of early initiation of antiretroviral therapy as well as tracing and testing of partners in patients with multiple HIV-1 infection.

Characterization of a Novel HIV-1 CRF01_AE/CRF07_BC Recombinant Strain Among Men Who Have Sex with Men in Liaoning, China.

CRF01_AE and CRF07_BC are two widespread human immunodeficiency virus type 1 (HIV-1) strains among different high-risk populations, including men who have sex with men (MSM), in China. This co-epidemic of various HIV strains enables the production of second-generation recombinants. In this study, we detected a novel HIV-1 CRF01_AE/CRF07_BC recombinant from LN321945, an MSM lived in Liaoning province, northeast China. The phylogenetic and recombination analyses indicated the near full-length genome (NFLG) sequence of LN321945 had six recombination breakpoints, with three CRF07_BC fragments inserted into a CRF01_AE backbone. Further subregion trees analysis revealed that both CRF01_AE and CRF07_BC fragments were derived from two predominant HIV-1 strains among MSM. In addition, the NFLG of LN321945 was revealed to be clustered closely to another CRF01_AE/CRF07_BC recombinant previously identified in Shaanxi province, northwest China, but these two recombinants had distinct recombination structure and origin of CRF01_AE fragments. Hence, this study identified a second-generation recombinant between the main strains circulating among MSM, indicating more complicated trend of HIV-1 epidemic in China.

Independent epidemic patterns of HIV-1 CRF01_AE lineages driven by mobile population in Shenzhen, an immigrant city of China.

Shenzhen, a city with >12 million migrant population, may play a key role in the spread of human immunodeficiency virus (HIV)-1 in China. The transmission dynamics of CRF01_AE, a predominant subtype in Shenzhen, is a good model to characterize the impact of human mobility on HIV-1 epidemic locally and nationally. We used phylodynamic and phylogeographic methods to estimate the viral transmission dynamics and migration trajectory of variable lineages based on 1,423 CRF01_AE sequences in Shenzhen sampled between 2006 and 2015. Eleven lineages of CRF01_AE were detected in Shenzhen. Of those, four main lineages originated during the 1990s. Their basic viral reproduction number (R 0) ranged 1.96-3.92. The effective viral reproduction number (R e ) of two lineages prevalent among heterosexuals/people who inject drugs had reduced <1 at the end of sampling, and the main sources were the intra-provincial immigrants (72 per cent) for one and local residents of Shenzhen (91 per cent) for another. Within two lineages among men who have sex with men (MSM), R e had been above or close to 1 at the end of sampling, and the immigrants from Jiangxi/Shaanxi and Hubei as sources accounted for 93 per cent and 68 per cent of all viral migration events, respectively. Moreover, no obvious recipients were found throughout the viral migration history for any lineage. Our findings demonstrate that HIV epidemic is declining in Shenzhen, which coincided with the initiation of the interventions during the 2000s. However, the obvious differences of the epidemic patterns between lineages emphasize the importance of further targeting interventions and continued molecular tracing, focusing on high-risk transmission sources among MSM.

Dynamics of HIV-1 Molecular Networks Reveal Effective Control of Large Transmission Clusters in an Area Affected by an Epidemic of Multiple HIV Subtypes.

This study reconstructed molecular networks of human immunodeficiency virus (HIV) transmission history in an area affected by an epidemic of multiple HIV-1 subtypes and assessed the efficacy of strengthened early antiretroviral therapy (ART) and regular interventions in preventing HIV spread. We collected demographic and clinical data of 2221 treatment-naïve HIV-1-infected patients in a long-term cohort in Shenyang, Northeast China, between 2008 and 2016. HIV pol gene sequencing was performed and molecular networks of CRF01_AE, CRF07_BC, and subtype B were inferred using HIV-TRACE with separate optimized genetic distance threshold. We identified 168 clusters containing ≥ 2 cases among CRF01_AE-, CRF07_BC-, and subtype B-infected cases, including 13 large clusters (≥ 10 cases). Individuals in large clusters were characterized by younger age, homosexual behavior, more recent infection, higher CD4 counts, and delayed/no ART (P < 0.001). The dynamics of large clusters were estimated by proportional detection rate (PDR), cluster growth predictor, and effective reproductive number (R e ). Most large clusters showed decreased or stable during the study period, indicating that expansion was slowing. The proportion of newly diagnosed cases in large clusters declined from 30 to 8% between 2008 and 2016, coinciding with an increase in early ART within 6 months after diagnosis from 24 to 79%, supporting the effectiveness of strengthened early ART and continuous regular interventions. In conclusion, molecular network analyses can thus be useful for evaluating the efficacy of interventions in epidemics with a complex HIV profile.

Time-Resolved Luminescence Properties of Laser-Fabricated Nano-diamonds.

In the study, well-crystallized nano-diamonds with an average size of 3.8 nm are obtained via femtosecond laser ablation. Both steady-state and transient luminescence are observed. The luminescence peaks of nano-diamonds shift from 380 to 495 nm when the excitation wavelength changes from 280 to 420 nm. After passivation by polyethylene glycol-400N, the surface of nano-diamonds is significantly oxidized, which is verified by Raman and UV-Vis absorption spectra. Furthermore, there is no change in all the luminescence wavelengths, although the maximum intensity increases by 10 times. Time-resolved luminescence spectra reveal that trapping states can be modified by surface passivation, and this leads to stronger luminescence with a longer lifetime.

Molecular network-based intervention brings us closer to ending the HIV pandemic.

Precise identification of HIV transmission among populations is a key step in public health responses. However, the HIV transmission network is usually difficult to determine. HIV molecular networks can be determined by phylogenetic approach, genetic distance-based approach, and a combination of both approaches. These approaches are increasingly used to identify transmission networks among populations, reconstruct the history of HIV spread, monitor the dynamics of HIV transmission, guide targeted intervention on key subpopulations, and assess the effects of interventions. Simulation and retrospective studies have demonstrated that these molecular network-based interventions are more cost-effective than random or traditional interventions. However, we still need to address several challenges to improve the practice of molecular network-guided targeting interventions to finally end the HIV epidemic. The data remain limited or difficult to obtain, and more automatic real-time tools are required. In addition, molecular and social networks must be combined, and technical parameters and ethnic issues warrant further studies.

Natural polymorphisms in HIV-1 CRF01_AE strain and profile of acquired drug resistance mutations in a long-term combination treatment cohort in northeastern China.

BACKGROUND: The impacts of genetic polymorphisms on drug resistance mutations (DRMs) among various HIV-1 subtypes have long been debated. In this study, we aimed to analyze the natural polymorphisms and acquired DRM profile in HIV-1 CRF01_AE-infected patients in a large first-line antiretroviral therapy (ART) cohort in northeastern China. METHODS: The natural polymorphisms of CRF01_AE were analyzed in 2034 patients from a long-term ART cohort in northeastern China. The polymorphisms in 105 treatment failure (TF) patients were compared with those in 1148 treatment success (TS) patients. The acquired DRM profile of 42 patients who experienced TF with tenofovir/lamivudine/efavirenz (TDF/3TC/EFV) treatment was analyzed by comparing the mutations at TF time point to those at baseline. The Stanford HIVdb algorithm was used to interpret the DRMs. Binomial distribution, McNemar test, Wilcoxon test and CorMut package were used to analyze the mutation rates and co-variation. Deep sequencing was used to analyze the evolutionary dynamics of co-variation. RESULTS: Before ART, there were significantly more natural polymorphisms of 31 sites on reverse transcriptase (RT) in CRF01_AE than subtype B HIV-1 (|Z value| ≥ 3), including five known drug resistance-associated sites (238, 118, 179, 103, and 40). However, only the polymorphism at site 75 was associated with TF (|Z value| ≥ 3). The mutation rate at 14 sites increased significantly at TF time point compared to baseline, with the most common DRMs comprising G190S/C, K65R, K101E/N/Q, M184 V/I, and V179D/I/A/T/E, ranging from 66.7 to 45.2%. Moreover, two unknown mutations (V75 L and L228R) increased by 19.0 and 11.9% respectively, and they were under positive selection (Ka/Ks > 1, log odds ratio [LOD] > 2) and were associated with several other DRMs (cKa/Ks > 1, LOD > 2). Deep sequencing of longitudinal plasma samples showed that L228R occurred simultaneously or followed the appearance of Y181C. CONCLUSION: The high levels of natural polymorphisms in CRF01_AE had little impact on treatment outcomes. The findings regarding potential new CRF01_AE-specific minor DRMs indicate the need for more studies on the drug resistance phenotype of CRF01_AE.

Cross-Continental Dispersal of Major HIV-1 CRF01_AE Clusters in China.

Since the 1990s, several distinct clusters of human immunodeficiency virus-type 1 (HIV-1) CRF01_AE related to a large epidemic in China have been identified, but it is yet poorly understood whether its transmission has dispersed globally. We aimed to characterize and quantify the genetic relationship of HIV-1 CRF01_AEs circulating in China and other countries. Using representative sequences of Chinese clusters as queries, all relevant CRF01_AE pol sequences in two large databases (the Los Alamos HIV sequence database and the UK HIV Drug Resistance Database) were selected with the online basic local alignment search (BLAST) tool. Phylogenetic and phylogeographic analyses were then carried out to characterize possible linkage of CRF01_AE strains between China and the rest of the world. We identified that 269 strains isolated in other parts of the world were associated with five major Chinese CRF01_AE clusters. 80.7% were located within CN.01AE.HST/IDU-2, most of which were born in Southeast Asia. 17.8% were clustered with CN.01AE.MSM-4 and -5. Two distinct sub-clusters associated with Chinese men who have sex with men (MSM) emerged in HK-United Kingdom and Japan after 2000. Our analysis suggests that HIV-1 CRF01_AE strains related to viral transmission in China were initially brought to the United Kingdom or other countries during the 1990s by Asian immigrants or returning international tourists from Southeast Asia, and then after having circulated among MSM in China for several years, these Chinese strains dispersed outside again, possibly through MSM network. This study provided evidence of regional and global dispersal of Chinese CRF01_AE strains. It would also help understand the global landscape of HIV epidemic associated with CRF01_AE transmission and highlight the need for further international collaborative study in this field.

Factors associated with high-risk low-level viremia leading to virologic failure: 16-year retrospective study of a Chinese antiretroviral therapy cohort.

BACKGROUND: Low level viremia (LLV) often occurs during antiretroviral therapy (ART) against HIV-1. However, whether LLV increases the risk of virologic failure (VF) is controversial because of the non-uniform definitions of LLV and VF. METHODS: A long-term first line regimen ART cohort from 2002 to 2018 from Shenyang, northeast China, was retrospectively studied. All participants were followed up every 3 to 6 months to evaluate the treatment effect. The high-risk LLV subgroups leading to VF (with strict standards) were explored with Cox proportional hazards model and linear mixed-effect model. The association factors of high-risk LLV were further explored using multivariate logistic regression analyses. RESULTS: A total of 2155 HIV-1 infected participants were included; of these, 38.7% showed LLV. Both high level LLV (HLLV) and any other level LLV coupled with high level blip (HLB) showed higher risk of VF (hazards ratios, HRHLLV = 5.93, and HRHLB = 2.84, p <  0.05 respectively). Moreover, HR increased with prolonged duration of LLV. Independent factors associated with high-risk LLV included the zenith baseline viral load (VL) above 6 log copies/ml (aOR = 3.49, p = 0.002), nadir baseline CD4 + T cell counts below 200 cells/mm3 (aOR = 1.78, p = 0.011), Manchu (aOR = 2.03, p = 0.003), ART over 60 months (aOR = 1.81, p = 0.004), AZT + 3TC + NVP (aOR = 2.26, p <  0.001) or DDI-based regimen (aOR = 9.96, p = 0.002), and subtype B' infection (aOR = 8.22, p = 0.001). CONCLUSIONS: In case of VF with strict standards, high-risk LLV leading to VF includes VL above 400 copies/ml, occurring at least once. Serious laboratory indicators or advanced stage of infection, long term ART and subtype B' infection might also predict the occurrence of high-risk LLV.

The S68G polymorphism is a compensatory mutation associated with the drug resistance mutation K65R in CRF01_AE strains.

BACKGROUND: The rate of S68G mutation in human immunodeficiency virus type 1 (HIV-1) reverse transcriptase has increased and is closely related to the K65R mutation among CRF01_AE-infected patients who failed treatment. We aimed to explore the temporal association of S68G and K65R mutations and disclose the role of the former on susceptibility to nucleotide/nucleoside reverse transcriptase inhibitor (NRTI) and viral replication with the K65R double mutations among CRF01_AE-infected patients who failed treatment. METHODS: The occurrence of S68G and K65R mutations was evaluated among HIV-1 of various subtypes in the global HIV Drug Resistance Database. The temporal association of S68G and K65R mutations was analyzed through next-generation sequencing in four CRF01_AE-infected patients who failed treatment with tenofovir/lamivudine/efavirenz. The impact of the S68G mutation on susceptibility to NRTI and replication fitness was analyzed using pseudovirus phenotypic resistance assays and growth competition assays, respectively. RESULTS: The frequency of the S68G mutation increased by 1.4-9.7% in almost all HIV subtypes and circulating recombinant forms in treatment-experienced patients, except subtype F. The S68G mutation often occurred in conjunction with the K65R mutation among RTI-treated patients, with frequencies ranging 21.1-61.7% in various subtypes. Next-generation sequencing revealed that the S68G mutation occurred following the K65R mutation in three of the four CRF01_AE-infected patients. In these three patients, there was no significant change detected in the half maximal inhibitory concentration for zidovudine, tenofovir, or lamivudine between the K65R and K65R/S68G mutations, as demonstrated by the phenotypic resistance assays. Virus stocks of the K65R and K65R/S68G mutations were mixed with 4:6, 1:1, and 9:1 and cultured for 13 days, the K65R/S68G mutants outgrew those of the K65R mutants irrespective of the input ratio. CONCLUSIONS: S68G may be a natural polymorphism and compensatory mutation of K65R selected by NRTIs in the CRF01_AE strain of HIV-1. This mutation does not affect susceptibility to NRTI; however, it improves the replication fitness of K65R mutants. This study deciphers the role of the S68G mutation in the HIV reverse transcriptase of the CRF01_AE strain and provides new evidence for the interpretation of drug-resistant mutations in non-B subtypes of HIV-1.