Background: Emerging HIV drug resistance caused by increased usage of antiretroviral drugs (ARV) could jeopardize the success of standardized HIV management protocols in resource-limited settings.
Objective: We aimed to characterize pretreatment HIV drug resistance (PDR) among HIV-positive individuals and risk factors in China in 2022.
Methods: This cross-sectional study was conducted using 2-stage systematic sampling according to the World Health Organization’s surveillance guidelines in 8 provincial-level administrative divisions in 2022. Demographic information and plasma samples were obtained from study participants. PDR was analyzed using the Stanford HIV drug resistance database, and the Tamura-Nei 93 model in HIV-TRACE was used to calculate pairwise matches with a genetic distance of 0.01 substitutions per site. Logistic regression was used to identify and estimate factors associated with PDR.
Results: PDR testing was conducted on 2568 participants in 2022. Of the participants, 34.8% (n=893) were aged 30-49 years, 81.4% (n=2091) were male, and 3.2% (n=81) had prior ARV exposure. The prevalence of PDR to protease and reverse transcriptase regions, nonnucleoside reverse transcriptase inhibitors, nucleoside reverse transcriptase inhibitors, and protease inhibitors were 7.4% (n=190), 6.3% (n=163), 1.2% (n=32), and 0.2% (n=5), respectively. Yunnan, Jilin, and Zhejiang had much higher PDR incidence than did Sichuan. The prevalence of nonnucleoside reverse transcriptase inhibitor–related drug resistance was 6.1% (n=157) for efavirenz and 6.3% (n=163) for nevirapine. Multivariable logistic regression models indicated that participants who had prior ARV exposure (odds ratio [OR] 7.45, 95% CI 4.50-12.34) and the CRF55_01B HIV subtype (OR 2.61, 95% CI 1.41-4.83) were significantly associated with PDR. Among 618 (24.2%) sequences (nodes) associated with 253 molecular transmission clusters (size range 2-13), drug resistance mutation sites included K103, E138, V179, P225, V106, V108, L210, T215, P225, K238, and A98.
Conclusions: The overall prevalence of PDR in China in 2022 was modest. Targeted genotypic PDR testing and medication compliance interventions must be urgently expanded to address PDR among newly diagnosed people living with HIV in China.
With the implementation of the “Treat all” program in 2016 in China, all people living with HIV are eligible to be treated, with their consent, regardless of their CD4 cell count and viral load level. Early treatment has been shown to effectively prolong the life span of people living with HIV and prevent secondary transmission. However, treatment of all individuals diagnosed with HIV has the potential to increase the prevalence of HIV drug resistance (DR) [- ]. Pretreatment HIV drug resistance (PDR) refers to resistance that is detected among antiretroviral drug (ARV)–naive people initiating antiretroviral therapy (ART) or people with prior ARV drug exposure [ ]. PDR may compromise the effectiveness of ART, increase mortality, and lead to secondary transmission of HIV [ ]. Thus, assessing PDR levels is important for informing national policies on ARV and for improving ART programs and service delivery.
In China, individual ART plans are determined primarily by national standards for first- and second-line regimens, and potential DR is not a major consideration among individuals preparing to initiate ART. Due partly to resource limitations, individual HIV genotype resistance testing is not routinely conducted prior to the initiation of ART in China. This approach to developing individual ART regimens may be associated with DR . A national meta-analysis indicated that PDR rapidly rose from 4.05% in 2011 to 5% in 2017. Specifically, resistance to nonnucleoside reverse transcriptase inhibitors (NNRTIs) increased from 2.15% to 3.81% from 2012 to 2017 [ ]. However, it is important to note that PDR testing and prevalence varies widely by geographic region. Other surveys conducted on PDR in certain provinces of China showed that the prevalence of PDR was 10.5% in Chongqing, 9.1% in Shenyang, 4.1% in Beijing, 6% in Shenzhen, and 10.8% in Pu’er County of Yunnan [ , - ]. The rates of PDR among HIV-positive individuals in China are comparatively lower than those in other countries, such as Cameroon (10.4%), Namibia (12.7%), Israel (12.1%), Mozambique (14.4%), Mexico (14.8%), the United States (22.5%), and Japan (12.5%) [ - ]. Yet China still lacks PDR surveillance methods that meet World Health Organization (WHO) standards. Assessing PDR prevalence with WHO standardized protocols will facilitate more accurate estimates of PDR prevalence in China, which will in turn better inform evidence-based ART policies in China.
Molecular transmission network analysis can be used to identify the maximum number of clusters and relatively active clusters based on genetic distance thresholds. Such methods have been used to evaluate the efficacy of ART on preventing secondary HIV transmission . Moreover, our understanding of DR has been greatly enhanced by studies integrating analysis of transmission characteristics, DR strains, and molecular transmission networks [ , ].
The WHO has identified DR strategic surveillance of PDR as a very high priority in countries using efavirenz (EFV)/nevirapine (NVP) for first-line ART and recommends implementation or updates every 3 years . Based on WHO PDR surveillance guidelines, the aim of this novel multicenter observational study is to investigate the prevalence of PDR and drug resistance mutations (DRMs) to partial pol regions, including protease (PR) and reverse transcriptase regions (RT), in diverse provinces in China.
Study Design and Study Population
We selected 8 provincial-level administrative divisions (PLADs) in China as the study area to conduct a large multicenter observational cross-sectional study. Two-stage systematic sampling in each PLAD was conducted according to the WHO PDR surveillance guidelines . Provinces were stratified into the following 3 categories based on their HIV incidence rates: high (>2.00 cases per 10,000 individuals; included Sichuan Province, Chongqing Municipality, and Yunnan Province), moderate (0.52-1.99 cases per 10,000 individuals; included Jiangsu Province and Zhejiang Province), and low (<0.52 cases per 10,000 individuals; included Jilin Province, Hebei Province, and Hubei Province).
Guided by WHO PDR surveillance guidelines, we calculated the minimum necessary sample size by assuming the prevalence of PDR was 10% and the CI range was ±5%. We sampled 11-20 clinics in each PLAD according to the number of HIV-positive individuals starting or restarting ART in the past 12 months. The minimum sample size in each PLAD was 158-179 samples. To ensure a minimum necessary sample size, we assumed that 40% of collected specimens would be lost to factors such as laboratory failure and a design effect. Expected and actual sample sizes for each province are shown in. Recruitment and specimen collection of HIV-positive individuals occurred in 2022 as individuals initiated ART.
The study eligibility criteria for individuals included the following: (1) tested HIV-positive, (2) ≥18 years old, (3) started or restarted first-line ART between January and June 2022 (includes women exposed to ARV drugs for preventing mother-to-child transmission of HIV, people who have received pre-exposure prophylaxis, and individuals reinitiating first-line ART), and (4) provided informed consent.
Each study participant was assigned a confidential and unique patient identifier used for linking basic epidemiological information and specimens. Basic sociodemographic data were collected, including age, sex, ethnicity, education level, marital status, occupation, route of infection, prior ARV exposure, time elapsed since HIV diagnosis, and PLADs. Clinical information included patients’ CD4 cell count before ART and HIV subtype.
Viral RNA was extracted from 200 µL plasma samples using the QIAamp viral RNA mini kit (Qiagen, Hilton, Germany) according to the operational instructions. Extracted RNA was used for first-round polymerase chain reaction, and the products were used for second-round polymerase chain reaction amplification; both rounds were performed in 25 µL volume reactions to amplify the HIV pol gene region (HXB2: 2253-3553 nt). This covered the PR region (4-99 amino acids) and the partial RT region (1-251 amino acids) . The amplified products were sequenced by Sanger sequencing.
Subtype and Drug Resistance Analysis
Sequencher (version 4.10.1, Gene Codes Corp) software was used to recombine and edit the original sequence fragments and nucleosides. BioEdit (version 7.0.9, Informer Technologies) software was used to realign the HXB2 reference sequence. The sequences were realigned, and some common reference sequences (including A to K) were used to identify the subtype by constructing neighbor joining phylogenetic trees in MEGA (version 7.0.26, MEGA Software). Subtypes with a bootstrap value higher than 70% were determined using HXB2 as the reference sequence. Sequences were realigned using the HIV database website . Realigned sequences were uploaded to the Stanford Drug Resistance Database [ ] to identify the degree of resistance of the 12 listed drugs that the WHO recommends for monitoring, namely, EFV, NVP, abacavir, zidovudine, lamivudine, tenofovir, emtricitabine, stavudine, didanosine, lopinavir/r, atazanavir/r, and darunavir/r. Corresponding DR mutant genes with levels of 1, 2, 3, 4, and 5 were classified as sensitive, potential resistance, low resistance, moderate resistance, and high resistance, respectively. Only mutations characterized as low, moderate, or high resistance were considered resistant. All possible sample mutual contamination, laboratory contamination, and other sequence quality controls were monitored by constructing neighbor joining phylogenetic trees and using the WHO HIVDR QC tool [ ].
Construction of the Molecular Network
Realigned sequences were conducted to calculate pairwise genetic distances though the Tamura-Nei 93 model in the transmission cluster engine, HIV-TRACE . Sequences no longer than 1000 bp or containing ≥5% ambiguities were excluded. A threshold genetic distance of 0.01 substitutions per site was used to observe transmission events that may occur within 5-6 years [ ]. Each individual in the molecular network was represented by a node, and each node in the molecular transmission network represented HIV-positive individuals who were matched with epidemiological information. We then linked nodes to each other if their pairwise genetic distance was up to 0.01 substitutions per site. Web pages were used to visualize molecular networks, which were based on technical guides for HIV transmission network monitoring and intervention [ ]. The results and subsequent molecular network diagram are accessible on the web [ ]. Our previous studies used other software, such as HYPHY to calculate the pairwise gene distance and Cytoscape for visualization [ , , ]. Clusters were defined as containing 2 or more nodes with HIV-positive individuals in the same cluster having potential for transmission.
All quantitative data analyses were performed using the Statistical Analysis System version 9.4 (SAS Institute). Two-sided P values <.05 were considered statistically significant. All variables were assessed using univariate binary logistic regression. Variables significantly associated (P<.05) with the outcome were entered into a stepwise regression model to remove collinear variables. The remaining variables were entered into a multivariate binary logistic regression model to identify correlates of PDR. The P value of the final multivariate logistic regression model was <.05.
This study was approved by the ethics committee of the National Center for AIDS/STD Control and Prevention, China Centers for Disease Control and Prevention (approval number X140617334).
Characteristics of the Study Participants
A total of 2869 HIV-positive individuals were initially recruited into the study. After removal of 16 HIV-positive individuals who were younger than 18 years and 9 individuals who had not started ART, the penultimate sample size was 2844. PDR test results were available for 90.3% (n=2568) of participants in RT/PR regions, thus yielding a final sample size of 2568. Of the 2568 study participants with valid PDR results, 34.8% (n=893) were aged 30-49 years, 81.4% (n=2091) were male, 86.8% (n=2229) were of Han ethnicity, 30.5% (n=783) had a junior high school education, 37% (n=950) were single, 36% (n=924) were farmers, and 57.6% (n=1480) and 37% (n=951) were infected through heterosexual and homosexual intercourse, respectively. In addition, 3.2% (n=81) of participants had prior ARV drug exposure, 84.5% (n=2170) were diagnosed in 2022, and 37.5% (n=962) had CD4 cell counts of 0-199 cells/mm3 prior to ART initiation ().
|Characteristic||Participants (n=2568), n (%)|
|Primary and below||797 (31)|
|Junior high school||783 (30.5)|
|Senior high school||360 (14)|
|Married or cohabiting||1144 (44.6)|
|Divorced or widowed||457 (17.8)|
|Route of infection|
|Heterosexual intercourse||1480 (57.6)|
|Homosexual intercourse||951 (37)|
|Intravenous drug use||21 (0.8)|
|Prior ARVa exposure|
|CD4 cell count before ARTb (cell/mm3)|
aARV: antiretroviral drugs.
bART: antiretroviral therapy.
PDR and Mutations
Among the 2568 included study participants, PDR was detected in 7.4% (n=190).shows the prevalence of overall PDR in Hebei (13/175; 7.4%), Jilin (19/206; 9.2%), Jiangsu (24/335; 7.2%), Zhejiang (34/385; 8.8%), Hubei (20/362; 5.5%), Chongqing (20/378; 5.3%), Sichuan (17/350; 4.9%) and Yunnan (43/377; 11.4%) in 2022. The prevalence of PDR in Yunnan (P=.002), Jilin (P=.047), and Zhejiang (P=.04) was much higher than that in Sichuan. Additionally, the prevalence of PDR in Yunnan was significantly higher than those in Chongqing (P=.002) and Hubei (P=.004).
Of the 2568 cases studied, the prevalence of PDR to NNRTIs, nucleoside reverse transcriptase inhibitors (NRTIs), and protease inhibitors (PIs) were 6.3% (n=163), 1.2% (n=32), and 0.2% (n=5), respectively. NNRTI-related DR to EFV and NVP was 6.1% (n=157) and 6.3% (n=163), respectively. There were 14 NNRTI-related resistance mutations; the most common mutations were K103N/S (n=81, 3.2%), followed by V179D/E/IL (n=39, 1.5%), and E138A/G/K (n=31, 1.2%). Among NRTI-related DR, the most frequent PDR drug was abacavir (n=19, 0.7%), followed by stavudine (n=17, 0.7%), emtricitabine (n=15, 0.6%), lamivudine (n=15, 0.6%), didanosine (n=13, 0.5%), zidovudine (n=10, 0.4%), and tenofovir (n=8, 0.3%). We detected 10 NRTI-related resistance mutations, of which M184V/I (n=13, 0.5%) was the most common. Even in the PR region, 0.2% (n=5) PI-related PDR to atazanavir/r and lopinavir/r was detected, including the PI-related DRMs L90M/IM (n=3, 0.1%), N88T/D (n=1, 0.1%), and V82IM (n=1, 0.1%) ().
|Province||Prevalence, n (%)||Participants (n=2568), n|
|PR/RTa (n=190)||NNRTIsb (n=163)||NRTIsc (n=32)||PIsd (n=5)|
|Hebei||13 (7.4)||9 (5.1)||4 (2.3)||0 (0)||175|
|Jilin||19 (9.2)||17 (8.2)||4 (1.9)||0 (0)||206|
|Jiangsu||24 (7.2)||22 (6.5)||1 (0.2)||1 (0.2)||335|
|Zhejiang||34 (8.8)||27 (7)||7 (1.8)||0 (0)||385|
|Hubei||20 (5.5)||19 (5.2)||1 (0.2)||0 (0)||362|
|Chongqing||20 (5.3)||14 (3.7)||7 (1.8)||1 (0.2)||378|
|Sichuan||17 (4.9)||16 (4.5)||0 (0)||1 (0.2)||350|
|Yunnan||43 (11.4)||39 (10.3)||8 (2.1)||2 (0.5)||377|
|All provinces||190 (7.4)||163 (6.3)||32 (1.2)||5 (0.2)||2568|
aPR/RT: protease and reverse transcriptase.
bNNRTIs: nonnucleoside reverse transcriptase inhibitors.
cNRTIs: nucleoside reverse transcriptase inhibitors.
dPIs: protease inhibitors.
|Antiretroviral drug||Prevalence, n||% (95% CI)||HIV drug resistance mutations and combination of mutations, n (%)|
|PR/RTa region (NNRTIsb, NRTIsc, PIsd)||190||7.4 (6.4-8.4)|
|Efavirenz||157||6.1 (5.2-7)||K103N/S: 81 (3.2)|
|Nevirapine||163||6.3 (5.4-7.3)||V179D/E/IL: 39 (1.5)|
V106A/M/I: 15 (0.6)
K101E/P: 9 (0.4)
Y181C: 9 (0.4)
A98G: 8 (0.5)
P225H: 8 (0.3)
Y188F/L/C: 6 (0.2)
V108I: 5 (0.2)
H221Y: 5 (0.2)
K238T: 3 (0.1)
|Abacavir||19||0.7 (0.4-1.1)||M184V/I: 13 (0.5)|
|Emtricitabine||15||0.6 (0.3-0.9)||T215A/S/I: 5 (0.2)|
|Lamivudine||15||0.6 (0.3-0.9)||T69D/DN/ADN/del: 5 (0.2)|
|Didanosine||13||0.5 (0.2-0.8)||D67N: 5 (0.2)|
|Stavudine||17||0.7 (0.3-1)||K70E/T/R: 5 (0.2)|
|Tenofovir||8||0.3 (0.1-0.5)||K65R: 3 (0.1)|
|Zidovudine||10||0.4 (0.1-0.6)||L210W/MRW: 2 (0.1)|
L74I: 2 (0.1)
K219Q: 1 (0.1)
Y115F: 1 (0.1)
|Atazanavir/r||4||0.1 (0-0.3)||L90M/IM: 3 (0.1)|
|Darunavir/r||0||N88T: 1 (0.1)|
|Lopinavir/r||4||0.1 (0-0.3)||V82IM: 1 (0.1)|
|Mutual drug resistance to NNRTIs and NRTIs||10||0.4 (0.1-0.6)|
aPR/RT: protease and reverse transcriptase.
bNNRTIs: nonnucleoside reverse transcriptase inhibitors.
cNRTIs: nucleoside reverse transcriptase inhibitors.
dPIs: protease inhibitors.
Factors Associated With PDR
lists the factors associated with PDR in PR/RT regions. Compared with participants without prior ARV drug exposure, those who had prior ARV drug exposure had 7.45 times greater odds of PDR (odds ratio [OR] 7.45, 95% CI 4.50-12.34). Compared to participants with the CRF07_BC subtype, those with the CRF55_01B subtype had 2.61 times greater odds of PDR (OR 2.61, 95% CI 1.41-4.83).
|Variable||Number (n=2568)||PDR (n=190), n (%)||Odds ratio (95% CI)||P value||Adjusted odds ratio (95% CI)||P value|
|30-49||893||68 (7.6)||0.74 (0.64-1.37)||.74||—||—|
|50-69||894||60 (6.7)||0.82 (0.55-1.21)||.32||—||—|
|≥70||175||13 (7.4)||0.91 (0.48-1.72)||.78||—||—|
|Female||477||41 (8.6)||1.23 (0.85-1.76)||.27||—||—|
|Other||339||32 (9.4)||1.37 (0.92-2.04)||.12||—||—|
|Primary and below||797||63 (7.9)||1.00||N/A||—||—|
|Junior high school||783||57 (7.3)||0.92 (0.631.33)||.64||—||—|
|Senior high school||360||21 (5.8)||0.72 (0.43-1.20)||.21||—||—|
|College||607||49 (8.1)||1.02 (0.69-1.51)||.91||—||—|
|Married or cohabiting||1144||72 (6.3)||0.69 (0.50-0.96)||.03||—||—|
|Divorced or widowed||457||34 (7.4)||0.83 (0.55-1.26)||.38||—||—|
|Other||1644||132 (8)||0.27 (0.16-2.64)||.10||—||—|
|Route of infectionc|
|Heterosexual intercourse||1480||102 (6.9)||1.00||N/A||—||—|
|Homosexual intercourse||951||72 (7.6)||1.11 (0.81-1.51)||.53||—||—|
|Intravenous drug use||21||6 (28.6)||5.41 (2.05-14.23)||<.001||—||—|
|Other||116||10 (8.6)||1.28 (0.65-2.51)||.48||—||—|
|Prior ARVd exposurec,e|
|Yes||81||27 (33.3)||7.13 (4.37-11.62)||<.001||7.45 (4.50-12.34)||.001|
|CD4 cell count before ART (cell/mm3)|
|≥200||1554||122 (7.9)||1.18 (0.86-1.61)||.31||—||—|
|Missing||52||3 (5.8)||0.85 (0.26-2.79)||.78||—||—|
|CRF08_BC||376||31 (8.2)||1.39 (0.89-2.16)||.14||1.09 (0.68-1.73)||.73|
|CRF01_AE||786||60 (7.6)||1.28 (0.89-1.83)||.18||1.27 (0.89-1.83)||.19|
|CRF55_01B||105||14 (13.3)||2.38 (1.29-4.39)||.006||2.61 (1.41-4.83)||.002|
|B||77||8 (10.4)||1.79 (0.83-3.88)||.14||1.61 (0.73-3.57)||.24|
|Other||87||8 (9.2)||1.57 (0.73-3.38)||.25||1.43 (0.65-3.14)||.37|
|Provincial-level administrative divisionsc|
|Hebei||175||13 (7.4)||1.57 (0.75-3.32)||.24||—||—|
|Jilin||206||19 (9.2)||1.99 (1.01-3.92)||.047||—||—|
|Jiangsu||335||24 (7.2)||1.51 (0.80-2.87)||.21||—||—|
|Zhejiang||385||34 (8.8)||1.90 (1.04-3.46)||.04||—||—|
|Hubei||362||20 (5.5)||1.15 (0.59-2.23)||.69||—||—|
|Chongqing||378||20 (5.3)||1.09 (0.56-2.13)||.79||—||—|
|Yunnan||377||43 (11.4)||2.52 (1.41-4.51)||.002||—||—|
aN/A: not applicable.
cVariables screened by stepwise regression.
dARV: antiretroviral drugs.
eVariables selected into multivariable logistic regression models (without prior ARV exposure serving as the reference group and with HIV subtype represented as 5 dummy variables with CRF07_BC serving as the reference group).
Molecular Transmission Network
Of the 2568 sequences in the PR/RT regions, 13 sequences were removed because the sequences in the PR/RT regions were shorter than 1000 bp or the number of mixed bases was >5%. Hence, a total of 2555 sequences in the PR/RT regions were used to construct a molecular transmission network. Under the optimal genetic distance threshold of 0.01 substitutions per site, 618 (24.2%) sequences (nodes) were linked to a total of 253 clusters (size range 2-13). We identified 38 sequences with PDR strains among 19 clusters. Most PDR strains linked to clusters had resistance to NNRTIs (36 of 38 sequences); DRM sites included K103, E138, V179, P225, V106, V108, L210, T215, P225, K238, and A98. Individual PDR status was not statistically significantly associated with belonging to a cluster. Cluster 01 with PDR consisted of 6 nodes in Jiangsu. Cluster 02 included HIV-positive individuals from 5 nodes in the same clinic sites in Zhejiang and 1 node in Hubei ().
In 2022, a cross-sectional survey was conducted to examine the prevalence of PDR among HIV-positive individuals in 8 PLADs in China. The results showed that the overall prevalence of PDR was 7.4%. According to the WHO’s qualitative classification of PDR prevalence levels, a prevalence below 5% is considered low, a prevalence between 5% and 15% is moderate, and a prevalence ≥15% is high . Therefore, the overall prevalence of PDR in 2022 in China and those of 7 PLADs (all but Sichuan) were considered moderate. The prevalence of PDR in 13 PLADs and 31 PLADs in China was 6.8% and 4.4% in 2017 and 2018, respectively [ , ]. Monitoring and surveillance are needed to assess the extent to which resistant strains have spread, and new technologies and strategies should be applied to effectively curb the spread of PDR.
Our study revealed that the PDR rate of NNRTIs was 6.3%. However, the PDR NNRTI prevalence was 10.3% for Yunnan. We observed higher PDR rates in Zhejiang and Yunnan in 2022 compared to surveys conducted in 2018 (3.7% in Zhejiang and 3.8% in Yunnan) and this increase in PDR prevalence was driven mainly by NNRTI resistance. According to WHO ART guidelines developed in 2016, it is recommended to consider implementing supplementary first-line antiviral treatment plans when the prevalence of PDR is equal to or greater than 10% for NNRTIs in HIV-positive individuals who are starting first-line ART . Further research is needed to identify the underlying reasons for the high NNRTI PDR rates in Yunnan. From 2014 to 2020, 21 of 30 countries that reported investigations to the WHO had an NVP or EFV PDR rate that exceeded 10% among HIV-positive individuals starting first-line ART [ ]. In our study, we found that the prevalence of EFV and NVP DR was 6.1% and 6.3%, respectively. The main DRM sites identified were K103, E138, V179, and V106. These results are consistent with earlier surveys conducted in China in 2015 and 2018 [ , ]. Studies from Zimbabwe, Kenya, and Argentina have identified K103 as the most common DRM site for NNRTIs used as first-line antiviral drugs [ - ]. National HIV treatment strategies and regimens should be reconsidered given the high prevalence of NNRTI DR. In 2023, the National Free HIV Antiviral Treatment Drug Treatment Manual was released in China, and rilpivirine was added as the first-line treatment regimen [ ]. The RPV resistance rate in our study was 5.3% (137/2568). The PDR after rilpivirine is applied to first-line treatment alternatives should be continuously monitored.
Our 2022 study indicated that PDR rates have increased in China over the past 7 years for NRTI DR (from 1.1% in 2015 and 0.8% in 2018 to 1.2% in 2022) and PI DR (from 0.02% in 2015 and 0.07% in 2018 to 0.2% in 2022) . Nonetheless, PDR rates for NRTIs and PIs in China generally appear to be lower than the WHO’s estimated regional PDR rates in Africa (NRTI PDR: 6.1%; PI PDR: 0.3%), the Americas (NRTI PDR: 6.4%; PI PDR: 0.8%), Southeast Asia (NRTI PDR: 3.1%; PI PDR: 0.4%), and the Western Pacific (NRTI PDR: 4%; PI PDR: 0%) [ ].
Results of the multivariable analysis indicated that individuals who had prior ARV drug exposure had 7.45 times greater odds of PDR compared to those without prior ARV drug exposure. Such findings are consistent with research from Mexico, African countries, and Northeast China [- ]. There was an increasing trend of the proportion of CRF55_01B samples, rising from 2.3% in 2015 to 3.9% in 2018 and 4.1% in 2022 in China [ ]. Additionally, participants with the CRF55_01B subtype had 2.61 times greater odds of PDR (OR 2.61, 95% CI 1.41-4.83) compared to participants whose subtype was CRF07_BC. This may be due to association with V179 as a characteristic DR site of CRF55_01B, and V179 being frequently associated with E138G and causing low resistance to EFV and NVP [ ]. Controlling the spread of CRF55_01B may be helpful for mitigating the proliferation of drug resistance. Furthermore, it is imperative to enhance the quality and effectiveness of intervention services for HIV-positive individuals.
In this study, we did not detect any statistically significant associations between having a PDR strain and belonging to clusters. Such results are similar with findings from Shijiazhuang city, Hebei Province and Liangshan region, Sichuan Province [, ]. One possible explanation for these findings is that the transmission potential of DR strains is weaker than that of nonresistant wild strains [ , ]. The predominant subtype among the sequences in the molecular transmission network was CRF07_BC (representing 46.6% of strains), which is consistent with findings from a 2018 national molecular transmission network survey in China [ ]. Clustered transmission of K103N DRMs in our study was observed in the molecular transmission network and is similar to findings from molecular transmission network studies in coastal Kenya [ ]. It is worth noting that 6 DR strains in cluster 01 were from HIV-positive individuals at identical sites in Jiangsu Province. Epidemiological data indicated that all 6 were adult men aged younger than 30 years, of whom 5 reported homosexual transmission and 1 reported heterosexual transmission. Cluster 02 contained 6 samples with DR strains from HIV-positive individuals in the same clinic sites in Zhejiang Province, of whom 5 reported heterosexual transmission and 1 reported homosexual transmission, suggesting a possible clustered transmission of DR strains in that area.
Our study has several limitations. Due to sampling limitations, clinic sites were established in different counties and districts, and only 20-80 patients were enrolled in each site during a certain period of time. This may have resulted in smaller observable clusters (containing 2-13 nodes) in the molecular network and a lower network access rate. In addition, when detecting drug-resistant strains using genotyping, Sanger sequencing may not be able to detect up to 20% of less prevalent PDR strains . This may have resulted in an underestimation of the PDR rate in HIV-positive individuals.
The findings from our 2022 PDR survey in China indicate a need to review and reconsider the standard ART protocol in some PLADs. Additional resources and strategies are needed in regions within China experiencing a higher prevalence of PDR. Possible interventions include extending pilot genotypic PDR testing to rapidly adjust medication regimens for HIV-positive individuals identified with PDR and targeted measures to improve medication compliance and prevention.
This work was supported by the Ministry of Science and Technology of China (grants 2022YFC2305205 and 2022YFC2305201), National Natural Science Foundation of China (grant 82260670), China Capital’s Funds for Health Improvement and Research (grant 2022-1G-3011), Beijing Municipal Science & Technology Commission (grant D161100000416002), Guangxi Bagui Honor Scholarship, and Chinese State Key Laboratory of Infectious Disease Prevention and Control. The authors thank the study participants for their support and the staff at the clinical sites and provincial Chinese Center for Disease Control and Prevention for their help in collecting blood samples and data.
HC, YR, HX, LL, YF, and YS were responsible for the study design and planning. HC, CS, YZ, XX, and ML were responsible for the statistics and figures. HC, J Hao, DL, XL, PX, DW, CS, J Hu, and JZ conducted the experiments and collected the data. HC, YR, and SWP drafted the manuscript. YR and HX guided the whole study and revised the article. All authors read and approved the final version of the manuscript.
Conflicts of Interest
Sample size of HIV-positive individuals in 8 provincial-level administrative divisions in 2022 in China.DOCX File , 24 KB
- Zhang T, Liao L, Shao Y, Feng Y, Ruan Y, Xing H. Relationship between drug resistance and death in HIV-infected patients receiving antiretroviral therapy - 7 PLADs, China, 2010-2019. China CDC Wkly 2021 May 02;3(14):291-297 [https://europepmc.org/abstract/MED/34594868] [CrossRef] [Medline]
- Volk JE, Nguyen DP, Hare CB, Marcus JL. HIV infection and drug resistance with unsupervised use of HIV pre-exposure prophylaxis. AIDS Res Hum Retroviruses 2018 May;34(4):329-330 [https://europepmc.org/abstract/MED/29262689] [CrossRef] [Medline]
- Zhao Y, Zhang M, Shi CX, Huang J, Zhang F, Rou K, et al. Mortality and virological failure among HIV-infected people who inject drugs on antiretroviral treatment in China: an observational cohort study. Drug Alcohol Depend 2017 Jan 01;170:189-197 [https://europepmc.org/abstract/MED/27987476] [CrossRef] [Medline]
- Yu F, Li Q, Wang L, Zhao H, Wu H, Yang S, et al. Drug resistance to HIV-1 integrase inhibitors among treatment-naive patients in Beijing, China. Pharmgenomics Pers Med 2022;15:195-203 [https://europepmc.org/abstract/MED/35300056] [CrossRef] [Medline]
- HIV drug resistance report 2019. World Health Organization. 2019. URL: https://apps.who.int/iris/handle/10665/325891 [accessed 2023-10-31]
- Beck IA, Levine M, McGrath CJ, Bii S, Milne RS, Kingoo JM, et al. Pre-treatment HIV-drug resistance associated with virologic outcome of first-line NNRTI-antiretroviral therapy: a cohort study in Kenya. EClinicalMedicine 2020 Jan;18:100239 [https://linkinghub.elsevier.com/retrieve/pii/S2589-5370(19)30244-5] [CrossRef] [Medline]
- National Center for AIDS/STD Control and Prevention. National Handbook of Free AIDS Antiviral Drug Treatment (4th Edition). Beijing: People's Medical Publishing House; 2016.
- 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 Jan;18:100238 [https://linkinghub.elsevier.com/retrieve/pii/S2589-5370(19)30243-3] [CrossRef] [Medline]
- Lu D, Wu H, Yarla NS, Xu B, Ding J, Lu T. HAART in HIV/AIDS treatments: future trends. Infect Disord Drug Targets 2018;18(1):15-22 [CrossRef] [Medline]
- Zhao B, Song W, Kang M, Dong X, Li X, Wang L, et al. Molecular network analysis reveals transmission of HIV-1 drug-resistant strains among newly diagnosed HIV-1 infections in a moderately HIV endemic city in China. Front Microbiol 2021;12:797771 [https://europepmc.org/abstract/MED/35069498] [CrossRef] [Medline]
- Ye J, Hao M, Xing H, Zhang F, Wu H, Lv W, et al. Transmitted HIV drug resistance among individuals with newly diagnosed HIV infection: a multicenter observational study. AIDS 2020 Mar 15;34(4):609-619 [CrossRef] [Medline]
- Liu M, He X, Deng R, Tang S, Harypursat V, Lu Y, et al. Pretreatment drug resistance in people living with HIV: a large retrospective cohort study in Chongqing, China. HIV Med 2022 Mar;23 Suppl 1(Suppl 1):95-105 [https://europepmc.org/abstract/MED/35293098] [CrossRef] [Medline]
- Hattori J, Shiino T, Gatanaga H, Mori H, Minami R, Uchida K, et al. Characteristics of transmitted drug-resistant HIV-1 in recently infected treatment-naive patients in Japan. J Acquir Immune Defic Syndr 2016 May 01;71(4):367-373 [CrossRef] [Medline]
- Taffa N, Roscoe C, Sawadogo S, De Klerk M, Baughman AL, Wolkon A, et al. Pretreatment HIV drug resistance among adults initiating ART in Namibia. J Antimicrob Chemother 2018 Dec 01;73(11):3137-3142 [https://europepmc.org/abstract/MED/30137412] [CrossRef] [Medline]
- Kassaye SG, Grossman Z, Balamane M, Johnston-White B, Liu C, Kumar P, et al. Transmitted HIV drug resistance is high and longstanding in metropolitan Washington, DC. Clin Infect Dis 2016 Oct 15;63(6):836-843 [https://europepmc.org/abstract/MED/27307507] [CrossRef] [Medline]
- Tchouwa GF, Eymard-Duvernay S, Cournil A, Lamare N, Serrano L, Butel C, et al. Prevalence of pretreatment HIV drug resistance in Cameroon following a nationally representative WHO survey. J Antimicrob Chemother 2018 Oct 01;73(9):2468-2474 [CrossRef] [Medline]
- Wagner T, Zuckerman NS, Halperin T, Chemtob D, Levy I, Elbirt D, et al. Epidemiology and transmitted HIV-1 drug resistance among treatment-naïve individuals in Israel, 2010-2018. Viruses 2021 Dec 31;14(1):71 [https://www.mdpi.com/resolver?pii=v14010071] [CrossRef] [Medline]
- Ismael N, Wilkinson E, Mahumane I, Gemusse H, Giandhari J, Bauhofer A, et al. Molecular epidemiology and trends in HIV-1 transmitted drug resistance in Mozambique 1999-2018. Viruses 2022 Oct 09;14(9):1992 [https://www.mdpi.com/resolver?pii=v14091992] [CrossRef] [Medline]
- 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 Mar 01;75(3):656-667 [https://europepmc.org/abstract/MED/31819984] [CrossRef] [Medline]
- Chen J, Chen H, Li J, Luo L, Kang R, Liang S, et al. Genetic network analysis of human immunodeficiency virus sexual transmission in rural Southwest China after the expansion of antiretroviral therapy: a population-based study. Front Microbiol 2022;13:962477 [https://europepmc.org/abstract/MED/36060743] [CrossRef] [Medline]
- Yuan D, Du Z, Zhou J, Ye L, Su L, Yang H, et al. HIV-1 subtype diversity, drug resistance, and genetic transmission networks in men who have sex with men with virologic failure in antiretroviral therapy in Sichuan, China, 2011 to 2017. Medicine (Baltimore) 2019 Oct;98(43):e17585 [https://europepmc.org/abstract/MED/31651864] [CrossRef] [Medline]
- HIV drug resistance strategy, 2021 update. World Health Organization. 2021. URL: https://apps.who.int/iris/handle/10665/343175 [accessed 2023-10-31]
- Surveillance of HIV drug resistance in adults initiating antiretroviral therapy (pre-treatment HIV drug resistance). World Health Organization. 2014. URL: http://apps.who.int/iris/bitstream/handle/10665/112802/9789241507196_eng.pdf?sequence=1 [accessed 2023-10-31]
- Li M, Liang S, Zhou C, Chen M, Liang S, Liu C, et al. HIV drug resistance mutations detection by next-generation sequencing during antiretroviral therapy interruption in China. Pathogens 2021 Mar 25;10(3):264 [https://www.mdpi.com/resolver?pii=pathogens10030264] [CrossRef] [Medline]
- HIVAlign. HIV Sequence Database. URL: https://www.hiv.lanl.gov/content/sequence/VIRALIGN/viralign.html [accessed 2023-10-31]
- HIV drug resistance database. Stanford University. URL: https://hivdb.stanford.edu/ [accessed 2023-10-31]
- WHO HIVDR QC Tool. URL: https://sequenceqc-dev.bccfe.ca/who_qc [accessed 2023-10-31]
- Kosakovsky Pond SL, Weaver S, Leigh Brown AJ, Wertheim JO. HIV-TRACE (TRAnsmission Cluster Engine): a tool for large scale molecular epidemiology of HIV-1 and other rapidly evolving pathogens. Mol Biol Evol 2018 Jul 01;35(7):1812-1819 [https://europepmc.org/abstract/MED/29401317] [CrossRef] [Medline]
- PS18-1802 list of attachments. Centers for Disease Control and Prevention. URL: https://www.cdc.gov/hiv/funding/announcements/ps18-1802/attachments.html [accessed 2023-10-31]
- WHO, UNODC, UNAIDS technical guide for countries to set targets for universal access to HIV prevention, treatment and care for injecting drug users. World Health Organization. URL: https://www.unodc.org/documents/hiv-aids/idu_target_setting_guide.pdf [accessed 2023-11-06]
- HIV-TRACE. URL: https://veg.github.io/hivtrace-viz/ [accessed 2023-10-31]
- Chen H, Hu J, Song C, Li M, Zhou Y, Dong A, et al. Molecular transmission network of pretreatment drug resistance among human immunodeficiency virus-positive individuals and the impact of virological failure on those who received antiretroviral therapy in China. Front Med (Lausanne) 2022;9:965836 [https://europepmc.org/abstract/MED/36106325] [CrossRef] [Medline]
- 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 May 01;35(6):947-955 [CrossRef] [Medline]
- Hao J, Zheng S, Gan M, Dong A, Kang R, Li M, et al. Changing proportions of HIV-1 subtypes and transmitted drug resistance among newly diagnosed HIV/AIDS individuals - China, 2015 and 2018. China CDC Wkly 2021 Dec 31;3(53):1133-1138 [https://europepmc.org/abstract/MED/35036036] [CrossRef] [Medline]
- Kang R, Liang S, Ma Y, Liang S, Xiao L, Zhang X, et al. Pretreatment HIV drug resistance in adults initiating antiretroviral therapy in China, 2017. Infect Dis Poverty 2020 May 24;9(1):54 [https://idpjournal.biomedcentral.com/articles/10.1186/s40249-020-00668-5] [CrossRef] [Medline]
- Consolidated guidelines on the use of antiretroviral drugs for treating and preventing HIV infection: recommendations for a public health approach, 2nd ed. World Health Organization. 2016. URL: https://www.who.int/publications/i/item/9789241549684 [accessed 2023-10-31]
- HIV drug resistance report 2021. World Health Organization. 2021. URL: https://apps.who.int/iris/handle/10665/349340 [accessed 2023-10-31]
- Mungati M, Mhangara M, Gonese E, Mugurungi O, Dzangare J, Ngwende S, et al. Pre-treatment drug resistance among patients initiating antiretroviral therapy (ART) in Zimbabwe: 2008-2010. BMC Res Notes 2016 Jul 10;9:302 [https://bmcresnotes.biomedcentral.com/articles/10.1186/s13104-016-2101-8] [CrossRef] [Medline]
- Bissio E, Barbás MG, Bouzas MB, Cudolá A, Salomón H, Espínola L, et al. Pretreatment HIV-1 drug resistance in Argentina: results from a surveillance study performed according to WHO-proposed new methodology in 2014-15. J Antimicrob Chemother 2017 Mar;72(2):504-510 [CrossRef] [Medline]
- Silverman RA, Beck IA, Kiptinness C, Levine M, Milne R, McGrath CJ, et al. Prevalence of pre-antiretroviral-treatment drug resistance by gender, age, and other factors in HIV-infected individuals initiating therapy in Kenya, 2013-2014. J Infect Dis 2017 Dec 19;216(12):1569-1578 [https://europepmc.org/abstract/MED/29040633] [CrossRef] [Medline]
- "National Free HIV Antiviral Drug Treatment Manual (2023 Edition)" officially released. National Center for AIDS/STD Control and Prevention, China CDC. URL: https://ncaids.chinacdc.cn/zxzx/zxdteff/202306/t20230616_266839.htm [accessed 2023-10-31]
- García-Morales C, Tapia-Trejo D, Matías-Florentino M, Quiroz-Morales VS, Dávila-Conn V, Beristain-Barreda Á, et al. HIV pretreatment drug resistance trends in Mexico City, 2017-2020. Pathogens 2021 Dec 08;10(12):1587 [https://www.mdpi.com/resolver?pii=pathogens10121587] [CrossRef] [Medline]
- Crowell TA, Danboise B, Parikh A, Esber A, Dear N, Coakley P, et al. Pretreatment and acquired antiretroviral drug resistance among persons living with HIV in four African countries. Clin Infect Dis 2021 Oct 05;73(7):e2311-e2322 [https://europepmc.org/abstract/MED/32785695] [CrossRef] [Medline]
- Gao L, Xia H, Zeng R, Wu Y, Zaongo SD, Hu Y, et al. Pre-treatment and acquired antiretroviral drug resistance among people living with HIV in Tianjin, China. HIV Med 2022 Mar;23 Suppl 1:84-94 [CrossRef] [Medline]
- Liu Y, Li H, Wang X, Han J, Jia L, Li T, et al. Natural presence of V179E and rising prevalence of E138G in HIV-1 reverse transcriptase in CRF55_01B viruses. Infect Genet Evol 2020 Jan;77:104098 [CrossRef] [Medline]
- Wang X, Liu X, Li F, Zhou H, Li J, Wang Y, et al. Epidemiological surveillance of HIV-1 transmitted drug resistance among newly diagnosed individuals in Shijiazhuang, northern China, 2014-2015. PLoS One 2018;13(6):e0198005 [https://dx.plos.org/10.1371/journal.pone.0198005] [CrossRef] [Medline]
- Liu L, Dong A, Liao L, Feng Y, Shao Y, Liang S, et al. Survey of pretreatment HIV drug resistance and genetic transmission network analysis among HIV patients in a high drug-use area of southwest China. Curr HIV Res 2019;17(6):441-451 [https://europepmc.org/abstract/MED/31778107] [CrossRef] [Medline]
- Hofstra LM, Nijhuis M, Pingen M, Mudrikova T, Riezebos-Brilman A, Simoons-Smit AM, et al. Evolution and viral characteristics of a long-term circulating resistant HIV-1 strain in a cluster of treatment-naive patients. J Antimicrob Chemother 2013 Jul;68(6):1246-1250 [CrossRef] [Medline]
- Hassan AS, Esbjörnsson J, Wahome E, Thiong'o A, Makau GN, Price MA, et al. HIV-1 subtype diversity, transmission networks and transmitted drug resistance amongst acute and early infected MSM populations from coastal Kenya. PLoS One 2018;13(12):e0206177 [https://dx.plos.org/10.1371/journal.pone.0206177] [CrossRef] [Medline]
|ART: antiretroviral therapy|
|ARV: antiretroviral drugs|
|DR: drug resistance|
|DRM: drug resistance mutation|
|NNRTIs: nonnucleoside reverse transcriptase inhibitors|
|NRTIs: nucleoside reverse transcriptase inhibitors|
|OR: odds ratio|
|PDR: pretreatment HIV drug resistance|
|PIs: protease inhibitors|
|PLADs: provincial-level administrative divisions|
|RT: reverse transcriptase|
|WHO: World Health Organization|
Edited by A Mavragani, T Sanchez; submitted 17.07.23; peer-reviewed by Y Tao, X Han; comments to author 23.08.23; revised version received 10.09.23; accepted 06.10.23; published 17.11.23Copyright
©Hongli Chen, Jingjing Hao, Jing Hu, Chang Song, Yesheng Zhou, Miaomiao Li, Jin Chen, Xiu Liu, Dong Wang, Xiaoshan Xu, Peixian Xin, Jiaxin Zhang, Lingjie Liao, Yi Feng, Dan Li, Stephen W Pan, Yiming Shao, Yuhua Ruan, Hui Xing. Originally published in JMIR Public Health and Surveillance (https://publichealth.jmir.org), 17.11.2023.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR Public Health and Surveillance, is properly cited. The complete bibliographic information, a link to the original publication on https://publichealth.jmir.org, as well as this copyright and license information must be included.