Skip to main content
Download PDF

HIV-1 drug resistance and genetic transmission network among newly diagnosed people living with HIV/AIDS in Ningbo, China between 2018 and 2021

Virology Journal volume 20, Article number: 233 (2023) Cite this article

Abstract

Background

As the HIV epidemic continues to grow, transmitted drug resistance(TDR) and determining relationship of HIV transmission are major barriers to reduce the risk of HIV transmissions.This study aimed to examine the molecular epidemiology and TDR and evaluated the transmission pattern among newly diagnosed people living with HIV/AIDS(PLWHA) in Ningbo city, which could contribute to the development of targeted precision interventions.

Methods

Consecutive cross-sectional surveys were conducted in Ningbo City between January 2018 and December 2021. The HIV-1 pol gene region was amplified and sequenced for drug resistance and genetic transmission network analysis. TDR was determined using the Stanford University HIV Drug Resistance Database. Genetic transmission network was visualized using Cytoscape with the genetic distance threshold of 0.013.

Results

A total of 1006 sequences were sequenced successfully, of which 61 (6.1%) showed evidence of TDR. The most common mutations were K103N (2.3%), E138A/G/Q (1.7%) and V179D/E (1.2%). 12 HIV-1 genotypes were identified, with CRF07_BC being the major genotype (43.3%, 332/767), followed by CRF01_AE (33.7%, 339/1006). 444 (44.1%) pol sequences formed 856 links within 120 transmission clusters in the network. An increasing trend in clustering rate between 2018 and 2021(χ2 = 9.546, P = 0.023) was observed. The odds of older age (≥ 60 years:OR = 2.038, 95%CI = 1.072 ~ 3.872, compared to < 25 years), HIV-1 genotypes (CRF07_BC: OR = 2.147, 95%CI = 1.582 ~ 2.914; CRF55_01B:OR = 2.217, 95%CI = 1.201 ~ 4.091, compared to CRF01_AE) were significantly related to clustering. Compared with CRF01_AE, CRF07_BC were prone to form larger clusters. The largest cluster with CRF07_BC was increased from 15 cases in 2018 to 83 cases in 2021.

Conclusions

This study revealed distribution of HIV-1 genotypes, and genetic transmission network were diverse and complex in Ningbo city. The prevalence of TDR was moderate, and NVP and EFV were high-level NNRTI resistance. Individuals aged ≥ 60 years old were more easily detected in the networks and CRF07_BC were prone to form rapid growth and larger clusters. These date suggested that surveillance and comprehensive intervention should be designed for key rapid growth clusters to reduce the potential risk factors of HIV-1 transmission.

Background

Despite substantial efforts to control human immunodeficiency virus-1 (HIV-1), acquired immunodeficiency syndrome (AIDS) is one of the most serious public health problem worldwide [1]. At the end of 2021, there were 38.4 million people living with HIV/AIDS(PLWHA), and 28.7 million people were accessing antiretroviral therapy(ART) [2]. World Health Organization(WHO) revised global recommendations that ART should be initiated in everyone living with HIV, regardless of CD4 + T lymphocytes counts in 2015 [3]. Early ART can suppress viral replication, reduce opportunistic infections, and significantly reduce morbidity and mortality among PLWHA [4,5,6]. Since HPTN 052 study published in interim form in 2011 [7], the strategy of treatment as prevention had been generalized highly effective in stopping reducing the risk of HIV transmissions globally [8, 9]. However, with the widespread use of ART, HIV drug resistance would lead to treatment failure in HIV infected patients, and even increase HIV transmission among individuals with risk behaviors and transmitted drug resistance(TDR) [10, 11].

As the HIV epidemic continues to spread, it is necessary to locate the source and relationship of HIV transmission [12, 13]. Genetic transmission network, which constructed based on the genetic information of people infected with HIV through gene distance between sequences, widely used in the study of early HIV cases detection, long-term monitoring of drug resistance and targeted precision intervention [11, 14,15,16,17]. Real-time genetic transmission network have been improved and further combined with social transmission network to judge new transmission events, identify high-risk spreaders and potential infected individuals and evaluate the effect of intervention measures [18,19,20]. Clustering analyses based on HIV drug resistance surveillance would rapidly detect and respond to emerging clusters of HIV infection to further reduce new transmissions, which is one pillar of the ending the HIV Epidemic plan in the United States [21].

Ningbo is an eastern coastal city of China, nearby Shanghai, with an area of 9365 km2 and a population of approximately 9.54 million people. We previously conducted an HIV tracing epidemiological survey among MSM in Ningbo City during 2018–2020 [22]. This study aimed to examine the molecular epidemiology and TDR of HIV-1 and evaluated the transmission pattern among newly diagnosed PLWHA in Ningbo city during 2018–2021, which could contribute to the development of targeted precision interventions.

Methods

Study participants

A cross-sectional survey was conducted in Ningbo City between January, 2018 and December, 2021. The inclusion criteria were as follows: (1) aged 18 years and above; (2) newly diagnosed HIV-1 cases; (3) had no received ART before enrollment; (4) agreed to participate in the survey and signed an informed consent. After providing written informed consent, blood samples were collected for CD4 + T lymphocytes counts, drug resistance and HIV sequencing.

Laboratory tests

Blood specimens were sent to Ningbo CDC for below laboratory tests. CD4 + T lymphocytes counts determined in fresh whole blood by Flow cytometry (Becton Dickinson, NJ, USA). The remaining whole blood sample was centrifuged at 3000 rpm to produce blood plasma preserved in a -80℃ freezer. Viral RNA was extracted from blood plasma using Viral RNA Mini Kit (Tianlong, Suzhou, China) according to the manufacturer’s instructions. The obtained RNA samples were amplified using reverse transcription polymerase chain reaction (PCR) and nested PCR for the pol (HXB2: 2147–3462, encoding the protease gene and the first 300 codons of the reverse transcriptase gene) gene regions of HIV-1. After electrophoretic analysis, the amplified positive products were sent to Hangzhou Qingke Zixi Biotechnology Co. Ltd. for purification and gene sequencing.

Sequence analysis

The sequences were assembled and adjusted with Sequencher v5.0 software(Genecodes, Ann Arbor, MI). The assembled sequences were aligned, edited, and analyzed with Bio-Edit 7.2 software(Genecodes, Ann Arbor, MI). International reference sequences were selected and downloaded from the HIV databases of the Los Alamos National Laboratory (http://hiv.lanl.gov). To identify the HIV-1 subtypes, the phylogenetic tree of Neighbor-joining was constructed using Mega 11.0 software. Bayesian Information Criterion (BIC) scores was performed to determine the evolutionary model. Considering the models with the lowest scores, the Kimura 2-parameter model with 1000 bootstrap replicates was determined to be the best fitting model [23]. Sequences with possible intersubtype recombination were analyzed using the Recombination Identification Program (RIP) tool (https://www.hiv.lanl.gov/content/sequence/RIP/RIP.html).

Drug resistance analysis

The sequences were submitted to the Stanford University HIV Resistance Database (http://hivdb.stanford.edu/) to describe and interpret HIV-1 TDR. TDR level was classified according to the Stanford Penalty Score as high (60), intermediate (30–59), or low (15–29) to the following drugs: PI, Protease inhibitor; NRTI, Nucleoside reverse transcriptase inhibitor; NNRTI, Non-nucleoside reverse transcriptase inhibitor; INSTI, Integrase strand transfer inhibitors; NFV: Nelfinavir; ABC: Abacavir; AZT: Zidovudine; D4T: Stavudine; DDI: Didanosine; FTC: Emtricitabine; 3TC, Lamivudine; TDF: Tenofovir; DOR: Doravirine; EFV: Efavirenz; ETR: Etravirine; RPV: Rilpivirine. Drug resistance mutations were analyzed using the CRP tool (http://cpr.stanford.edu/cpr.cgi).

Genetic transmission network analysis

The genetic transmission network was inferred based on the nucleotide genetic distance (GD) between HIV-1 pol sequences from each participant. Hyphy2.2.4 software was used to calculate the GD between pairings based on the TN93 model. The GD threshold that could identify the maximum number of clusters in the genetic network was chosen for analysis [24]. Transmission partner was defined that the GD between two sequences was below the distance threshold. Cytoscape 3.6.0 software was used to construct HIV molecular transmission network diagram based on the optimal distance threshold of different gene distances and the number of molecular clusters. The GD threshold of 0.013 was chosen to identified the maximum number of clusters in the genetic network. A node in the network represents a case or sequence. The connection between two nodes was an edge or link. The number of node connections was the degree, indicating its importance in the network. The higher the degree, the higher the inferred communication relationship with more people, the higher the communication risk.The nucleotide sequences were submitted to GenBank under the accession numbers OR521366-OR522371.

Statistical analysis

Characteristics of all participants were described by categorical variables presented as absolute values and percentages. The demographic information and distribution of genotype were examined by chi-square tests. Univariate and multivariable forward stepwise logistic regression models were performed to examine risk factors associated with study participants within genetic transmission networks. The statistical significance was defined as P < 0.05. All statistical analyses were performed in SPSS (version 21.0, IBM, Armonk, NY, USA).

Results

Participant characteristics

A total of 1097 newly diagnosed PLWH were enrolled in this study, and we successfully sequenced and analyzed the samples collected from 1006 (91.7%) individuals. Among them, the median age was 40 years (interquartile range 28–54 years), ranging from 18 to 84 years. There were 86.3% (868/1006) and 13.7% (138/1006) males and females, respectively. Regarding marital status, 37.5% (377/1006) were single, 42.8% (431/1006) were married, and 19.7% (198/1006) were divorce or death. Most participants were middle school or less education level (61.8%), of Han ethnicity (95.8%), living non local (56.5%), heterosexual transmission (52.7%), having < 6 sexual partners (79.6%). The proportions of CD4 + T lymphocytes counts < 200 cells/uL, 200 ~ 499 cells/uL and ≥ 500 cells/uL were 37.7% (379/1006), 52.2% (525/1006) and 10.1% (102/1006), respectively (Table 1).

Table 1 Characteristics of study participants between three groups with different HIV-1 genotypes
Full size table

Drug resistance analysis

The prevalence of TDR among newly diagnosed PLWH was 6.1% (61/1006). TDR-associated mutations to NNRTI, NRTI and PI were detected in 7.9% (79/1006), 2.7% (27/1006) and 0.7% (7/1006). No INSTI-related mutations were identified. As for the NNRTI resistance, NVP (2.9%) was the highest proportion of high-level resistance, followed by EFV (2.7%). Most exhibited a low-level resistance to RPV (2.8%) and ETR (1.3%). The most frequent mutations were K103N (2.3%), E138A/G/Q (1.7%) and V179D/E (1.2%). As for the NRTI resistance, FTC and 3TC presented the highest levels (0.7%), which was related to the mutation of M184MV/V (0.7%). As for the PI resistance, the HIV-1 strains only exhibited resistance to NFV (0.7%), and mutations were M46L/I(0.7%). The above results presented in Fig. 1.

Fig. 1
figure 1

(A) Proportion of HIV-1 transmitted drug resistance (TDR) mutations to PIs, NRTIs, and NNRTIs. (B) Levels of HIV-1 TDR associated mutations to different ART drug among 1006 newly diagnosed HIV/AIDS patients. PI, Protease inhibitor;NRTI, Nucleoside reverse transcriptase inhibitor; NNRTI, Non-nucleoside reverse transcriptase inhibitor; NFV: Nelfinavir;ABC: Abacavir; AZT: Zidovudine; D4T: Stavudine; DDI: Didanosine; FTC: Emtricitabine;3TC, Lamivudine;TDF, Tenofovir; DOR,Doravirine; EFV,Efavirenz; ETR: Etravirine; RPV: Rilpivirine

Full size image

Genotype analysis

Among all participants, the most common genotype was CRF07_BC (40.7%, 409/1006), followed by CRF01_AE (33.7%, 339/1006). The other genotypes were CRF08_BC (6.9%, 69/1006), subtype C (5.4%, 54/1006), CRF55_01B (5.0%, 50/1006), subtype B (2.6%, 26/1006), unique recombinant forms (URFs) (2.6%, 26/1006), CRF85_BC (1.5%, 15/1006), CRF01_BC (0.6%, 6/1006), CRF59_01B (0.4%, 4/1006), CRF67_01B (0.4%, 4/1006) and CRF68_01B (0.4%, 4/1006). The URFs included 19 CRF01AE/C, four CRF01AE/B/C, two subtype B/C, one CRF01AE/B. As shown in Table 1, the distribution of genotype revealed significant differences according to gender (χ2 = 25.799, P < 0.001), age (χ2 = 51.906, P < 0.001), education level (χ2 = 22.299, P < 0.001), route of infection (χ2 = 42.681, P < 0.001), CD4 + T lymphocytes counts (χ2 = 19.400, P = 0.038), and TDR (χ2 = 25.686, P < 0.001).

Genetic transmission network analysis

Overall, 444 (44.1%) pol sequences formed 856 links within 120 transmission clusters in the network. Comparing the sequences included to those not included in the transmission networks, we observed significant differences in age (χ2 = 4.778, P = 0.029), number of sexual partners in a life time (χ2 = 5.253, P = 0.022), HIV-1 genotypes (χ2 = 27.539, P < 0.001) and CD4 + T lymphocytes counts (χ2 = 6.581, P = 0.037). After adjustment for confounders, older age (OR = 2.038, 95%CI = 1.072 ~ 3.872, P = 0.030; ≥60 vs. <25 years) was associated with higher adjusted odds of clustering. Additionally, compared with CRF01_AE, CRF07_BC (OR = 2.147, 95%CI = 1.582 ~ 2.914, P = 0.011) and CRF55_01B (OR = 2.217, 95%CI = 1.201 ~ 4.091, P < 0.001) were significantly related to clustering (Table 2).

Table 2 Factors associated with study participants within genetic transmission networks in Ningbo, China
Full size table

As seen in the set of Fig. 2, the clustering rates were 35.4% (67/189), 43.8% (49/112), 43.5% (137/315), 49.0% (191/390) between 2018 and 2021, respectively. The proportions of links with GD < 0.005 were 19.0% (36/189), 14.3% (16/112), 14.0% (44/315), 13.8% (54/390), respectively. We observed an increasing trend in clustering rate between 2018 and 2021 (χ2 = 9.546, P = 0.023). The proportion of links with GD < 0.005 did not vary significantly during the four years (χ2 = 3.153, P = 0.369).

Fig. 2
figure 2

Enrolment and clustering rate in Ningbo, China between 2018 and 2021.Dispersion measures indicate 95% confidence intervals

Full size image

Of the 444 sequences evaluated, 154 (34.7%) had one link, and 290 (65.3%) had more than two links. There were 105 (87.5%) clusters with size < 5, 10 (8.3%) clusters with 5 ≤ size < 10, and five (4.2%) clusters with size ≥ 10. The five large clusters had 83, 21, 14, 10 and 10 individuals, respectively. The median degrees were 4, 20, 11, 7.5 and 9 respectively. The proportions of links with GDs < 0.005 were 19.0% (84/213), 42.9% (85/198) and 22.5% (16/71), 2.8% (1/36), 39.5% (17/43) respectively.

As shown in Fig. 3, among 444 participants who were included in these networks, 375 (84.5%) were male and 69 (15.5%) were female. We analyzed these 120 clusters, 46 (38.3%) clusters were CRF07_BC, 42 (35.0%) were CRF01_AE, 8 (6.6%) were subtype C, 7 (5.8%) were CRF55_01B, 5 (4.2%) were CRF85_BC and CRF08_BC, 3 (2.5%) were URFs, 2 (1.6%) were subtype B, 1 (0.8%) were CRF01_BC and CRF59_01B. Of 211 CRF07_BC sequences, 137 (64.9%) were connected to at least one other sequence with genetics distance < 0.005. Compared with CRF01_AE, CRF07_BC were prone to form larger clusters. The largest cluster with CRF07_BC was increased from 15 cases in 2018 to 83 cases in 2021. Considering clusters with three or more nodes, the most frequent links included persons with 25–39 years (35.2%), followed 40–59 years (31.7%). The proportion of links with GD < 0.005 varied across age: 41.0%, 31.7%, 35.7% and 42.9% in persons < 25, 25–39, 40–59 and ≥ 60 years, respectively, and was significantly higher in persons ≥ 60 years compared to all other age teams (χ2 = 14.456, P = 0.002).Considering clusters included all the four years, there were two clusters with 83 and 21 nodes, respectively. The annual growth rate of participants who were included in these two clusters were 36.3% and 76.5%. We also observed that 34.4% (21/61) of TDR cases were included in the transmission. All these TDR cases were found in 11 clusters with size < 6. The possibility of TDR clustering was higher in gender being male (80.6%), infection route being heterosexual transmission (58.1%) and genotype being CRF01_AE (51.6%). Of 856 links in the genetic transmission networks, 13 (1.5%) links were TDR connected to TDR, 9 (1.1%) were TDR linked to No-TDR, and 834 (97.4%) were No-TDR connected to No-TDR.

Fig. 3
figure 3

Genetic transmission network among newly diagnosed people living with HIV/AIDS in Ningbo, China between 2018 and 2021.Clusters are shown ordered by size in each panel. Links are coloured by genetic distance.Shapes of the nodes represent different gender.Colour of the nodes represent different participant characteristics,including genotype,age,enrolment year and transmitted drug resistance

Full size image

Discussion

The AIDS epidemic has a serious impact on health, economy and society of people all over the world. Meanwhile, TDR has become the focus of AIDS prevention and treatment. Furthermore, traditional epidemiological studies focus on the risk of individual HIV infection, with less consideration given to the impact of individual interactions on HIV transmission. Therefore, this study described TDR and genotype in newly diagnosed PLWHA in Ningbo city between 2018 and 2021, and further analyzed the HIV genetic network according to the GD. The results showed that the prevalence of TDR was 6.1% and 12 HIV-1 genotypes and 120 transmission clusters were identified, which revealed the complexity and diversity of the AIDS epidemic in Ningbo City.

Since China had adopted Free ART Program in 2003 and ‘Treat for all’ policy in 2016, free virus load test were implemented for all of the PLWHA receiving ART once per year [25].However, free viral load test is usually offered to patients until a year after ART to assess the effectiveness of treatment. Due to lack of resistance monitoring and testing before ART, the prevalence of TDR rapidly rising in recent years [26]. In our study, the prevalence of TDR was 6.1% and TDR-associated mutations to NNRTI was 7.9%, which was moderate (5–15%) according to the WHO definition [27]. It was higher than the prevalence determined in Sichuan [28], Guangdong [29] and Jiangsu [30] in China. In this study, we found that NVP and EFV showed high-level NNRTI resistance, with the most frequent mutations being K103N, E138A/G/Q and V179D/E. Additionally, in case of NRTI, the mutation M184MV/V can cause the highest level resistance to FTC and 3TC, which were consistent to the national study in China [26]. Moreover, 0.7% of subjects carried PI resistance related the mutation M46L/I which conferred TDR to NFV. Currently, 3TC, EFV and NVP are widely-used free drugs in the first-line ART regimens in China [31]. Long-term use of limited number of drugs facilitated the generation and spread of TDR. Therefore, the emergence of these mutations to 3TC, EFV and NVP should be closely monitored. In particular, drug-resistance testing prior to ART initiation should be performed to reduce spread of TDR. Furthermore, the introduction and development of longer-acting ART should be considered in China.

As is reported, HIV-1 genotype played an important role in HIV transmission and responses to ART [32]. Our study showed that the major epidemic genotypes were CRF07_BC and CRF01_AE, which was consistent to the prior study [33].And the distribution of genotype revealed significant differences according to gender, age, education level route of infection, CD4 + T lymphocytes counts and TDR. Currently, CRF08_BC transmission population has changed from intravenous drug users to heterosexuals in the highest risk groups in China [34]. CRF08_BC appeared to be a distinctive strain in Ningbo city. CRF08_BC was also the predominant genotype among individuals aged more than 40 years, middle school or less educational level, heterosexually transmitted individuals in our study. The prevalence of TDR was higher than other genotypes among individuals infected with CRF08_BC. This finding revealed urgent measures should be taken to interrupt the spread of CRF08_BC, particularly among elderly individuals. CRF55_01B, which was first detected in MSM, had quickly spread to heterosexuals, according to our results and earlier study [35]. Besides, CRF55_01B may have a higher transmission risk than CRF01_AE and CRF07_BC [36]. Thus, we should continue to concern the prevention and control of CRF55_01B.

We constructed HIV-1 genetic transmission network analysis to understand the transmission characteristic among PLWHA in Ningbo city. Our study showed that 444 individuals could be segregated into 120 transmission clusters. After adjusting for confounders, we observed that age and genotypes were significantly related to clustering. GD < 0.005 might be more appropriate for distinguishing rapidly growing clusters [37]. In this study, individuals aged over 60 years old were more likely to enter the transmission network and the proportions of links with GD < 0.005 were higher than other age teams, indicated that they increased the risk of local HIV transmission and should be monitored as a priority [38]. Moreover, compared with CRF01_AE, CRF07_BC were prone to form larger clusters and the largest cluster with CRF07_BC was grown rapid, which homosexual transmission was the major route of transmission. Furthermore, most of CRF07_BC sequences were connected to smaller GD. We found clustering rate was significantly increased between 2018 and 2021. It suggested that comprehensive intervention should be designed for key rapid growth clusters, especially for those with multiple links [39]. In addition, we identified that 34.4% of TDR cases were included in 11 clusters, mostly composed of heterosexuals. However, heterosexual contacts present small clusters, indicating a low level of forward transmission, but still represent a high proportion of TDR transmission in Ningbo city. These results suggested that local government should continue monitoring TDR transmission to conduct targeted interventions and control epidemics in a timely fashion [40].

Limitation of this study

There are a few inherent limitations in this study. Firstly, the sample composition may not be representative all PLWAH in Ningbo city. Secondly, molecular transmission cluster only represented a group of highly associated infected individuals, and could not reflect direct transmission relationship. Thirdly, as integrase was not included in the amplified pol sequence, this study did not identify INSTI-related drug resistance mutations. Finally, with limited time and funds, it was based on sequence analysis of the pol gene, which may not holistically represent the epidemic.

Conclusions

This study revealed distribution of HIV-1 genotypes and genetic transmission network were diverse and complex in Ningbo city. The prevalence of TDR was moderate, and NVP and EFV were high-level NNRTI resistance. Individuals aged ≥ 60 years old were more easily detected in the networks and CRF07_BC were prone to form rapid growth and larger clusters. These date suggested that surveillance and comprehensive intervention should be designed for key rapid growth clusters to reduce the potential risk factors of HIV-1 transmission.

Data Availability

The datasets used and/or analyzed during this study is not publicly available, but may be available from the corresponding author upon reasonable request, and with permission from Ningbo Municipal Center for Disease Control and Prevention.

Abbreviations

HIV:

human immunodeficiency virus

AIDS:

acquired immunodeficiency syndrome

PLWHA:

people living with HIV/AIDS

ART:

antiretroviral therapy

WHO:

World Health Organization

TDR:

transmitted drug resistance

PCR:

polymerase chain reaction

RIP:

Recombination Identification Program

PI:

Protease inhibitor

NRTI:

Nucleoside reverse transcriptase inhibitor

NNRTI:

Non-nucleoside reverse transcriptase inhibitor

INSTI:

Integrase strand transfer inhibitors

NFV:

Nelfinavir

ABC:

Abacavir

AZT:

Zidovudine

D4T:

Stavudine

DDI:

Didanosine

FTC:

Emtricitabine

3TC:

Lamivudine

TDF:

Tenofovir

DOR:

Doravirine

EFV:

Efavirenz

ETR:

Etravirine

RPV:

Rilpivirine

GD:

genetic distance

URFs:

unique recombinant forms

References

  1. De Cock KM, Jaffe HW, Curran JW. Reflections on 40 years of AIDS. Emerg Infect Dis. 2021;27:1553–60.

    Article  PubMed  PubMed Central  Google Scholar 

  2. UNAIDS, Global HIV. & AIDS statistics — Fact sheet. https://www.unaids.org/en/resources/fact-sheet. Accessed 19 May 2023.

  3. WHO. Guideline on when to start antiretroviral therapy and on pre-exposure prophylaxis for HIV. https://www.who.int/publications-detail-redirect/9789241509565. Accessed 20 May 2023.

  4. Olson AD, Walker AS, Suthar AB, Sabin C, Bucher HC, Jarrin I, et al. Limiting cumulative HIV Viremia Copy-Years by Early Treatment reduces risk of AIDS and Death. J Acquir Immune Defic Syndr. 2016;73:100–8.

    Article  PubMed  PubMed Central  Google Scholar 

  5. TEMPRANO ANRS 12136 Study Group, Danel C, Moh R, Gabillard D, Badje A, Le Carrou J, et al. A trial of early antiretrovirals and Isoniazid Preventive Therapy in Africa. N Engl J Med. 2015;373:808–22.

    Article  Google Scholar 

  6. INSIGHT START Study Group, Lundgren JD, Babiker AG, Gordin F, Emery S, Grund B, et al. Initiation of antiretroviral therapy in early asymptomatic HIV infection. N Engl J Med. 2015;373:795–807.

    Article  Google Scholar 

  7. Cohen MS, Chen YQ, McCauley M, Gamble T, Hosseinipour MC, Kumarasamy N, et al. Prevention of HIV-1 infection with early antiretroviral therapy. N Engl J Med. 2011;365:493–505.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Rodger AJ, Cambiano V, Bruun T, Vernazza P, Collins S, Degen O, et al. Risk of HIV transmission through condomless sex in serodifferent gay couples with the HIV-positive partner taking suppressive antiretroviral therapy (PARTNER): final results of a multicentre, prospective, observational study. Lancet. 2019;393:2428–38.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Kroon EDMB, Phanuphak N, Shattock AJ, Fletcher JLK, Pinyakorn S, Chomchey N, et al. Acute HIV infection detection and immediate treatment estimated to reduce transmission by 89% among men who have sex with men in Bangkok. J Int AIDS Soc. 2017;20:21708.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Kiekens A, Dierckx de Casterlé B, Pellizzer G, Mosha IH, Mosha F, Rinke de Wit TF, et al. Exploring the mechanisms behind HIV drug resistance in sub-saharan Africa: conceptual mapping of a complex adaptive system based on multi-disciplinary expert insights. BMC Public Health. 2022;22:1–15.

    Article  Google Scholar 

  11. Pang X, Tang K, He Q, Huang J, Fang N, Zhou X, et al. HIV drug resistance and HIV transmission risk factors among newly diagnosed individuals in Southwest China. BMC Infect Dis. 2021;21:160.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Leitner T, Romero-Severson E. Phylogenetic patterns recover known HIV epidemiological relationships and reveal common transmission of multiple variants. Nat Microbiol. 2018;3:983–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Lundgren E, Romero-Severson E, Albert J, Leitner T. Combining biomarker and virus phylogenetic models improves HIV-1 epidemiological source identification. PLoS Comput Biol. 2022;18:e1009741.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Liu M, Han X, Zhao B, An M, He W, Wang Z, et al. Dynamics of HIV-1 Molecular Networks reveal Effective Control of large transmission clusters in an area affected by an epidemic of multiple HIV subtypes. Front Microbiol. 2020;11:604993.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Zhao B, Qiu Y, Song W, Kang M, Dong X, Li X, et al. Undiagnosed HIV infections may drive HIV Transmission in the era of treat all: a deep-sampling Molecular Network Study in Northeast China during 2016 to 2019. Viruses. 2022;14:1895.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. Matías-Florentino M, Chaillon A, Ávila-Ríos S, Mehta SR, Paz-Juárez HE, Becerril-Rodríguez MA, et al. Pretreatment HIV drug resistance spread within transmission clusters in Mexico City. J Antimicrob Chemother. 2020;75:656–67.

    Article  PubMed  Google Scholar 

  17. Han X, Zhao B, An M, Zhong P, Shang H. Molecular network-based intervention brings us closer to ending the HIV pandemic. Front Med. 2020;14:136–48.

    Article  PubMed  Google Scholar 

  18. Dávila-Conn V, García-Morales C, Matías-Florentino M, López-Ortiz E, Paz-Juárez HE, Beristain-Barreda Á, et al. Characteristics and growth of the genetic HIV transmission network of Mexico City during 2020. J Int AIDS Soc. 2021;24:e25836.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Pasquale DK, Doherty IA, Sampson LA, Hué S, Leone PA, Sebastian J, et al. Leveraging phylogenetics to understand HIV Transmission and Partner Notification Networks. J Acquir Immune Defic Syndr. 2018;78:367–75.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Kang R, Li J, Chen H, Tang Z, Pan SW, Luo L, et al. Using longitudinal genetic-network study to understand HIV treatment-as-prevention. AIDS. 2021;35:947–55.

    Article  PubMed  CAS  Google Scholar 

  21. Fauci AS, Redfield RR, Sigounas G, Weahkee MD, Giroir BP. Ending the HIV Epidemic: a plan for the United States. JAMA. 2019;321:844–5.

    Article  PubMed  Google Scholar 

  22. Hong H, Zhang DD, Jiang HB, Shi HB, Tan SW, Gu WZ, et al. [HIV infection and related factors of traceability efficiency among sex partners of HIV positive men who have sex with men]. Zhonghua Liu Xing Bing Xue Za Zhi. 2021;42:2100–5.

    PubMed  CAS  Google Scholar 

  23. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol. 2018;35:1547–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Hassan AS, Pybus OG, Sanders EJ, Albert J, Esbjörnsson J. Defining HIV-1 transmission clusters based on sequence data. AIDS. 2017;31:1211–22.

    Article  PubMed  Google Scholar 

  25. Cao W, Hsieh E, Li T. Optimizing treatment for adults with HIV/AIDS in China: successes over two decades and remaining Challenges. Curr HIV/AIDS Rep. 2020;17:26–34.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Zuo L, Liu K, Liu H, Hu Y, Zhang Z, Qin J, et al. Trend of HIV-1 drug resistance in China: a systematic review and meta-analysis of data accumulated over 17 years (2001–2017). eClinicalMedicine. 2020;18:100238.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Bennett DE, Camacho RJ, Otelea D, Kuritzkes DR, Fleury H, Kiuchi M, et al. Drug resistance mutations for surveillance of transmitted HIV-1 drug-resistance: 2009 update. PLoS ONE. 2009;4:e4724.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Zhou C, Liang S, Li Y, Zhang Y, Li L, Ye L, et al. Characterization of HIV-1 molecular epidemiology and transmitted drug-resistance in newly diagnosed HIV-infected patients in Sichuan, China. BMC Infect Dis. 2022;22:602.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Lan Y, Li L, He X, Hu F, Deng X, Cai W, et al. Transmitted drug resistance and transmission clusters among HIV-1 treatment-naïve patients in Guangdong, China: a cross-sectional study. Virol J. 2021;18:181.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Yin Y, Zhou Y, Lu J, Liu X, Shi L, Fu G, et al. Molecular Epidemiology of HIV-1 in Jiangsu Province, Southeast China: genotypes and HIV-1 transmission networks among newly diagnosed Men having sex with men in 2017. AIDS Res Hum Retroviruses. 2021;37:62–9.

    Article  PubMed  CAS  Google Scholar 

  31. Zhou F, Kominski GF, Qian H-Z, Wang J, Duan S, Guo Z, et al. Expenditures for the care of HIV-infected patients in rural areas in China’s antiretroviral therapy programs. BMC Med. 2011;9:6.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Cao Z, Li J, Chen H, Song C, Shen Z, Zhou X, et al. Effects of HIV-1 genotype on baseline CD4 + cell count and mortality before and after antiretroviral therapy. Sci Rep. 2020;10:15875.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Fan Q, Zhang J, Luo M, Feng Y, Ge R, Yan Y, et al. Molecular genetics and epidemiological characteristics of HIV-1 epidemic strains in various sexual risk behaviour groups in developed Eastern China, 2017–2020. Emerg Microbes Infect. 2022;11:2326–39.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Li K, Liu M, Chen H, Li J, Liang Y, Feng Y, et al. Using molecular transmission networks to understand the epidemic characteristics of HIV-1 CRF08_BC across China. Emerg Microbes Infect. 2021;10:497–506.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Gan M, Zheng S, Hao J, Ruan Y, Liao L, Shao Y, et al. The prevalence of CRF55_01B among HIV-1 strain and its connection with traffic development in China. Emerg Microbes Infect. 2021;10:256–65.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Wei L, Lu X, Li H, Zheng C, Li G, Yang Z, et al. Impact of HIV-1 CRF55_01B infection on CD4 counts and viral load in men who have sex with men naive to antiretroviral treatment. The Lancet. 2018;392:43.

    Article  Google Scholar 

  37. Ding X, Chaillon A, Pan X, Zhang J, Zhong P, He L, et al. Characterizing genetic transmission networks among newly diagnosed HIV-1 infected individuals in eastern China: 2012–2016. PLoS ONE. 2022;17:e0269973.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Zhang F, Yang Y, Liang N, Liang H, Chen Y, Lin Z, et al. Transmission network and phylogenetic analysis reveal older male-centered transmission of CRF01_AE and CRF07_BC in Guangxi, China. Emerg Microbes Infect. 2023;12:2147023.

    Article  PubMed  Google Scholar 

  39. Wertheim JO, Panneer N, France AM, Saduvala N, Oster AM. Incident infection in high-priority HIV molecular transmission clusters in the United States. AIDS. 2020;34:1187–93.

    Article  PubMed  Google Scholar 

  40. Beckwith CG, Min S, Manne A, Novitsky V, Howison M, Liu T, et al. HIV Drug Resistance and Transmission Networks among a justice-involved Population at the time of Community Reentry in Washington, D.C. AIDS Res Hum Retroviruses. 2021;37:903–12.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to express our gratitude to the participants, primary clinicians and nurses who participated in this study.

Funding

This study was supported by the Zhejiang Medical and Health Technology Discipling (2021KY1017), Natural Science Foundation of Zhejiang Province (LQ20H260005), Ningbo Health Branding Subject Fund (PPXK2018-10), and Zhejiang Medical Key Discipline (07–013). The funding organizations had no role in the design and conduct of this study.

Author information

Author notes

  1. Hang Hong, Chunlan Tang contributed equally to this work.

Authors and Affiliations

  1. School of Public health, Health Science Center, Ningbo University, Ningbo, Zhengjiang, 315211, China

    Hang Hong, Chunlan Tang, Ting Fang & Guozhang Xu

  2. Ningbo Center for Disease Control and Prevention, Ningbo, Zhengjiang, 315010, China

    Yuhui Liu & Haibo Jiang

Authors
  1. Hang Hong
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Chunlan Tang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yuhui Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Haibo Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Ting Fang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Guozhang Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

HH designed the study, HH and GX coordinated the study, HS, YL,HJ collected the data, HH and TF did the primary data analysis, HH and CT contributed to part of the analysis. HH and CT drafted the paper. All authors contributed to results interpretation interpreted the data, revised the article, and approved the final version and manuscript revision. GX had access to all the data and was responsible for the final decision to submit the manuscript for publication.

Corresponding author

Correspondence to Guozhang Xu.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

The study protocol was reviewed and approved by the Institutional Review Board of the Ningbo CDC in China (201706 and 201908). Participants provided written informed consent to use sequencing data both for HIV drug resistance and genetic transmission network studies.

Consent for publication

Not applicable.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hong, H., Tang, C., Liu, Y. et al. HIV-1 drug resistance and genetic transmission network among newly diagnosed people living with HIV/AIDS in Ningbo, China between 2018 and 2021. Virol J 20, 233 (2023). https://doi.org/10.1186/s12985-023-02193-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-023-02193-x

Keywords

Virology Journal

ISSN: 1743-422X