Trends in Antibiotic Consumption and Antimicrobial Resistance in China: An Ecological Analysis from 2016 to 2022

Fei Zhao Zinan Zhao Di Chen Bolin Zhu Tianqi Zhang Yuanchao Zhu Xuelin Sun Liang Liang Nan Zheng Lili Zou Wenfeng Xu Sirui Guo Yue Wang Ming Zhao Xin Hu Pengfei Jin

Fei Zhao, Zinan Zhao, Di Chen, Bolin Zhu, Tianqi Zhang, Yuanchao Zhu, Xuelin Sun, Liang Liang, Nan Zheng, Lili Zou, Wenfeng Xu, Sirui Guo, Yue Wang, Ming Zhao, Xin Hu, Pengfei Jin. Trends in Antibiotic Consumption and Antimicrobial Resistance in China: An Ecological Analysis from 2016 to 2022[J]. Biomedical and Environmental Sciences. doi: 10.3967/bes2025.169
Citation: Fei Zhao, Zinan Zhao, Di Chen, Bolin Zhu, Tianqi Zhang, Yuanchao Zhu, Xuelin Sun, Liang Liang, Nan Zheng, Lili Zou, Wenfeng Xu, Sirui Guo, Yue Wang, Ming Zhao, Xin Hu, Pengfei Jin. Trends in Antibiotic Consumption and Antimicrobial Resistance in China: An Ecological Analysis from 2016 to 2022[J]. Biomedical and Environmental Sciences. doi: 10.3967/bes2025.169

doi: 10.3967/bes2025.169

Trends in Antibiotic Consumption and Antimicrobial Resistance in China: An Ecological Analysis from 2016 to 2022

More Information
    Author Bio:

    Fei Zhao, Associate Researcher, PhD, Research Direction: Pharmacoepidemiology, Tel: 86-10-85133648, E-mail: zhaofei21cn@126.com

    Zinan Zhao, Chief Pharmacist, Master's Degree, Research Direction: Clinical Pharmacy, Tel: 86-10-85133648, E-mail: bellazhao2015@163.com

    Corresponding author: Pengfei Jin, Prof, Tel: 86-10-85133620, E-mail: j790101@163.com
  • PFJ, XH, ZNZ, and FZ conceptualized the study. FZ and ZNZ performed the analyses. DC, BLZ, TQZ, YCZ, XLS, LL, NZ, LLZ, MZ, WFX, YW and SRG were involved in data curation and database management. FZ and TQZ designed this study. ZNZ, DC, YCZ, LL and BLZ accessed and verified the data. FZ, ZNZ, XH and PFJ interpreted the data. FZ wrote the first draft of the manuscript, which was edited and reviewed by all the authors. All the authors have read and approved the final version of the manuscript. All authors had full access to all data in the study and had the final responsibility for the decision to submit the manuscript for publication.
  • All authors declare no competing interests.
  • Not applicable.
  • &These authors contributed equally to this work.
    • 关键词:
    •  / 
    •  / 
    •  / 
    •  
    PFJ, XH, ZNZ, and FZ conceptualized the study. FZ and ZNZ performed the analyses. DC, BLZ, TQZ, YCZ, XLS, LL, NZ, LLZ, MZ, WFX, YW and SRG were involved in data curation and database management. FZ and TQZ designed this study. ZNZ, DC, YCZ, LL and BLZ accessed and verified the data. FZ, ZNZ, XH and PFJ interpreted the data. FZ wrote the first draft of the manuscript, which was edited and reviewed by all the authors. All the authors have read and approved the final version of the manuscript. All authors had full access to all data in the study and had the final responsibility for the decision to submit the manuscript for publication.
    All authors declare no competing interests.
    Not applicable.
    &These authors contributed equally to this work.
    注释:
    1) Authors’ Contributions: 2) Competing Interests: 3) Ethic:
  • Figure  1.  Drug-resistant proportions of CRO and VRE in China from 2016 to 2022. Decreasing (CRAB and CRPA) and increasing (CRKP) trends.

    Figure  2.  Drug-resistant proportions of CRO and VRE across the provinces in China from 2016 to 2022. A: Drug resistance proportions of CRAB, CRPA, CRKP, and CREC. B: Proportions of drug-resistant VR E. faecium and VR E. faecalis. Interprovincial heterogeneity in the resistance proportions has been observed across China.

    Figure  3.  The total consumption of 15 antibiotic classes in China from 2016 to 2022.

    Figure  4.  The consumption of 15 antibiotic classes across provinces in China from 2016 to 2022.

    Table  1.   Association analysis of CRO drug resistance with single−factor inclusion in hierarchical model in China from 2016 to 2022

    Variable CRAB CRPA CRKP CREC
    β (95% CI) p β (95% CI) p β (95% CI) p β (95% CI) p
    Year −0.04
    (−0.06, −0.02)
    <0.001 −0.05
    (−0.06, −0.04)
    <0.001 0.06
    (0.03, 0.09)
    <0.001 9.87×10−4
    (−0.01, 0.02)
    0.896
    Population 6.96×10−5
    (2.33×10−5, 11.59×10−5)
    0.003 4.52×10−5
    (1.05×10−6, 8.94×10−5)
    0.045 1.55×10−4
    (5.70×10−5, 2.53×10−4)
    0.002 7.83×10−5
    (2.64×10−5, 1.30×10−4)
    0.003
    Aminoglycosides 3.34×10−7
    (−2.35×10−7, 9.03×10−7)
    0.250 1.90×10−8
    (−3.67×10−7, 4.05×10−7)
    0.923 −6.71×10−7
    (−1.29×10−6, −5.60×10−8)
    0.033 1.54×10−7
    (−3.98×10−7, 7.06×10−7)
    0.584
    Amphenicols 2.92×10−6
    (−5.34×10−6, 1.12×10−5)
    0.489 5.51×10−6
    (4.51×10−7, 1.06×10−5)
    0.033 −7.23×10−6
    (−1.53×10−5, 8.70×10−7)
    0.080 −6.51×10−6
    (−1.31×10−5, 5.68×10−8)
    0.052
    β−lactamase inhibitors 2.63×10−8
    (−4.59×10−8, 9.85×10−8)
    0.476 −4.19×10−8
    (−9.28×10−8, 8.97×10−9)
    0.106 −8.63×10−8
    (−1.66×10−7, −6.45×10−9)
    0.034 4.15×10−9
    (−6.94×10−8, 7.77×10−8)
    0.912
    Carbapenems 1.07×10−7
    (−1.67×10−7, 3.82×10−7)
    0.444 −2.42×10−8
    (−2.32×10−7, 1.84×10−7)
    0.819 −2.34×10−7
    (−5.87×10−7, 1.18×10−7)
    0.192 2.02×10−7
    (−6.57×10−8, 4.69×10−7)
    0.139
    Cephalosporins 1.00×10−8
    (−4.93×10−9, 2.50×10−8)
    0.189 −1.37×10−9
    (−1.47×10−8, 1.20×10−8)
    0.841 −2.13×10−8
    (−4.16×10−8, −1.04×10−9)
    0.039 9.73×10−10
    (−1.72×10−8, 1.92×10−8)
    0.916
    Vancomycin & norvancomycin 1.24×10−7
    (−1.73×10−6, 1.92×10−6)
    0.896 −4.58×10−7
    (−1.92×10−6, 1.00×10−6)
    0.538 −1.71×10−6
    (−4.15×10−6, 7.42×10−7)
    0.172 1.55×10−6
    (−8.48×10−8, 3.18×10−6)
    0.063
    Teicoplanin 2.31×10−6
    (−1.86×10−6, 6.48×10−6)
    0.278 5.23×10−7
    (−2.50×10−6, 3.54×10−6)
    0.734 −2.21×10−6
    (−7.11×10−6, 2.69×10−6)
    0.376 3.21×10−6
    (−1.33×10−6, 7.75×10−6)
    0.166
    Lincomycin−Clindamycin 2.07×10−8
    (−1.43×10−7, 1.85×10−7)
    0.804 −2.76×10−8
    (−1.52×10−7, 9.63×10−8)
    0.662 1.70×10−8
    (−1.95×10−7, 2.29×10−7)
    0.875 7.49×10−8
    (−1.06×10−7, 2.55×10−7)
    0.416
    Macrolides 1.27×10−8
    (−7.27×10−9, 3.26×10−8)
    0.213 6.07×10−9
    (−8.00×10−9, 2.01×10−8)
    0.398 −3.07×10−8
    (−5.28×10−8, −8.57×10−9)
    0.007 −1.09×10−9
    (−2.45×10−8, 2.23×10−8)
    0.928
    Oxazolidinones −3.21×10−8
    (−9.66×10−7, 9.01×10−7)
    0.946 −7.48×10−7
    (−1.61×10−6, 1.10×10−7)
    0.088 6.44×10−8
    (−1.29×10−6, 1.41×10−6)
    0.925 2.09×10−7
    (−5.26×10−7, 9.44×10−7)
    0.577
    Penicillin 2.74×10−9
    (−2.17×10−8, 2.72×10−8)
    0.826 −1.46×10−8
    (−3.17×10−8, 2.60×10−9)
    0.096 −4.56×10−8
    (−6.96×10−8, −2.16×10−8)
    <0.001 7.94×10−10
    (−1.97×10−8, 2.13×10−8)
    0.940
    Polymyxin −1.67×10−6
    (−6.49×10−6, 3.15×106)
    0.498 −5.26×10−6
    (−9.97×10−6, −5.47×10−7)
    0.029 −6.42×10−6
    (−1.54×10−5, 2.55×10−6)
    0.162 −1.59×10−6
    (−4.98×10−6, 1.80×10−6)
    0.357
    Quinolones 8.73×10−9
    (−2.02×10−8, 3.77×10−8)
    0.554 −7.44×10−9
    (−2.87×10−8, 1.38×10−8)
    0.492 −2.51×10−8
    (−5.74×10−8, 7.26×10−9)
    0.129 8.73×10−9
    (−1.76×10−8, 3.51×10−8)
    0.516
    Sulfonamides −2.00×10−8
    (−4.31×10−7, 3.91×10−7)
    0.924 −6.80×10−8
    (−3.52×10−7, 2.16×10−7)
    0.638 −3.24×10−7
    (−8.61×10−7, 2.12×10−7)
    0.236 2.08×10−7
    (−2.26×10−7, 6.41×10−7)
    0.348
    Tetracyclines −5.21×10−9
    (−4.12×10−8, 3.08×10−8)
    0.777 −2.68×10−8
    (−6.28×10−8, 9.29×10−9)
    0.146 3.06×10−8
    (−1.53×10−8, 7.64×10−8)
    0.191 1.73×10−8
    (−1.22×10−8, 4.68×10−8)
    0.251
      Note. CRAB: Acinetobacter baumannii; CRPA: Pseudomonas aeruginosa; CRKP: Klebsiella pneumoniae; CREC: Escherichia coli. 95% CI: 95% confidence interval.
    下载: 导出CSV

    Table  2.   Association analysis of CRO drug resistance with multifactor inclusion in hierarchical model in China from 2016 to 2022

    Variabl×10 CRAB CRPA CRKP CR×10C
    β (95% CI) p β (95% CI) p β (95% CI) p β (95% CI) p
    Year −2.31×10−2 (−4.28×10−2, −3.32×10−3) 0.022 −5.70×10−2 (−7.35×10−2, −4.05×10−2) <0.001 8.12×10−2 (5.13×10−2, 1.11×10−1) <0.001 −4.81×10−3 (−3.08×10−2, 2.12×10−2) 0.717
    Population 8.34×10−5 (2.41×10−5, 1.43×10−4) 0.006 5.03×10−5 (5.05×10−6, 9.55×10−5) 0.029 1.20×10−4 (3.63×10−6, 2.37×10−4) 0.043 6.54×10−5 (1.03×10−5, 1.21×10−4) 0.020
    Aminoglycosides −2.29×10−7 (−8.23×10−7, 3.65×10−7) 0.449 −1.82×10−7 (−6.41×10−7, 2.77×10−7) 0.437 −1.94×10−7 (−9.81×10−7, 5.92×10−7) 0.628 5.41×10−7 (−1.44×10−7, 1.23×10−6) 0.122
    Amphenicols −1.37×10−6 (−6.68×10−6, 3.94×10−6) 0.613 2.03×10−7 (−4.00×10−6, 4.41×10−6) 0.924 2.66×10−6 (−4.30×10−6, 9.61×10−6) 0.454 −4.40×10−6 (−1.09×10−5, 2.06×10−6) 0.182
    β−lactamase inhibitors −4.96×10−8 (−1.63×10−7, 6.38×10−8) 0.391 −1.26×10−7 (−2.20×10−7, −3.16×10−8) 0.009 −1.65×10−8 (−1.74×10−7, 1.41×10−7) 0.838 −1.14×10−7 (−2.45×10−7, 1.75×10−8) 0.089
    Carbapenems 5.75×10−7 (1.60×10−7, 9.89×10−7) 0.007 5.26×10−7 (2.08×10−7, 8.45×10−7) 0.001 4.61×10−7 (−8.85×10−8, 1.01×10−6) 0.100 2.91×10−7 (−1.73×10−7, 7.55×10−7) 0.219
    Cephalosporins 3.09×10−8 (−9.83×10−9, 7.16×10−8) 0.137 −4.29×10−9 (−3.66×10−8, 2.80×10−8) 0.794 2.12×10−8 (−3.32×10−8, 7.56×10−8) 0.445 −3.11×10−8 (−8.08×10−8, 1.87×10−8) 0.221
    Vancomycin & norvancomycin −2.21×10−6 (−5.56×10−6, 1.15×10−6) 0.197 −6.03×10−7 (−3.28×10−6, 2.08×10−6) 0.660 1.61×10−6 (−2.94×10−6, 6.16×10−6) 0.488 3.64×10−6 (−1.46×10−7, 7.43×10−6) 0.060
    Teicoplanin 7.82×10−10 (−4.11×10−6, 4.11×10−6) 1.000 −1.34×10−6 (−4.56×10−6, 1.87×10−6) 0.413 1.80×10−7 (−5.21×10−6, 5.57×10−6) 0.948 −7.90×10−7 (−5.37×10−6, 3.79×10−6) 0.735
    Lincomycin−Clindamycin −7.22×10−8 (−2.89×10−7, 1.45×10−7) 0.514 −3.76×10−8 (−2.20×10−7, 1.45×10−7) 0.686 1.25×10−7 (−1.89×10−7, 4.39×10−7) 0.435 6.01×10−9 (−2.62×10−7, 2.74×10−7) 0.965
    Macrolides −1.90×10−8 (−6.32×10−8, 2.52×10−8) 0.400 1.15×10−8 (−2.41×10−8, 4.72×10−8) 0.527 −7.35×10−8 (−1.30×10−7, −1.68×10−8) 0.011 −1.68×10−8 (−7.07×10−8, 3.70×10−8) 0.540
    Oxazolidinones −1.28×10−6 (−3.40×10−6, 8.32×10−7) 0.234 −1.50×10−7 (−1.78×10−6, 1.48×10−6) 0.857 −4.01×10−7 (−2.83×10−6, 2.03×10−6) 0.747 7.93×10−7 (−1.56×10−6, 3.15×10−6) 0.509
    Penicillin −2.13×10−8 (−5.97×10−8, 1.72×10−8) 0.278 −3.57×10−8 (−6.75×10−8, −3.87×10−9) 0.028 −5.81×10−8 (−1.09×10−7, −6.90×10−9) 0.026 −2.50×10−8 (−7.17×10−8, 2.16×10−8) 0.293
    Polymyxin 2.45×10−6 (−3.91×10−6, 8.80×10−6) 0.451 3.84×10−6 (−1.29×10−6, 8.97×10−6) 0.142 −1.93×10−5 (−2.78×10−5, −1.08×10−5) <0.001 −1.13×10−5 (−1.92×10−5, −3.42×10−6) 0.005
    Quinolones 3.79×10−8 (−4.18×10−8, 1.18×10−7) 0.351 3.85×10−8 (−2.38×10−8, 1.01×10−7) 0.226 5.86×10−8 (−4.04×10−8, 1.57×10−7) 0.246 5.70×10−8 (−3.32×10−8, 1.47×10−7) 0.215
    Sulfonamides 1.76×10−8 (−3.67×10−7, 4.02×10−7) 0.929 2.51×10−7 (−8.00×10−8, 5.82×10−7) 0.137 −4.24×10−7 (−1.13×10−6, 2.86×10−7) 0.242 2.06×10−7 (−3.41×10−7, 7.52×10−7) 0.461
    Tetracyclines −4.39×10−8 (−1.03×10−7, 1.49×10−8) 0.143 −1.47×10−8 (−6.14×10−8, 3.21×10−8) 0.539 5.17×10−8 (−2.46×10−8, 1.28×10−7) 0.184 3.96×10−11 (−6.79×10−8, 6.80×10−8) 0.999
      Note. CRAB: Acinetobacter baumannii; CRPA: Pseudomonas aeruginosa; CRKP: Klebsiella pneumoniae; CREC: Escherichia coli. 95% CI: 95% confidence interval.
    下载: 导出CSV

    Table  3.   Association analysis of VRE drug resistance with single−factor inclusion in hierarchical model in China from 2016 to 2022

    Variable VR E. faecium VR E. faecalis
    β (95% CI) p β (95% CI) p
    Year −0.13 (−0.19, −0.06) < 0.001 −0.14 (−0.19, −0.10) < 0.001
    Population 1.47×10−5 (−7.00×10−5, 9.94×10−5) 0.733 −4.39×10−5 (−1.05×10−4, 1.73×10−5) 0.160
    Aminoglycosides 2.18×10−6 (1.14×10−6, 3.23×10−6) < 0.001 −3.34×10−8 (−1.03×10−6, 9.60×10−7) 0.947
    Amphenicols 1.72×10−5 (1.05×10−6, 3.33×10−5) 0.037 1.65×10−5 (1.70×10−6, 3.14×10−5) 0.029
    β−lactamase inhibitors 3.39×10−7 (1.96×10−7, 4.83×10−7) < 0.001 −5.53×10−8 (−1.99×10−7, 8.82×10−8) 0.450
    Carbapenems 1.42×10−6 (9.41×10−7, 1.90×10−6) < 0.001 −1.81×10−7 (−7.24×10−7, 3.62×10−7) 0.514
    Cephalosporins 7.50×10−8 (5.00×10−8, 9.99×10−8) < 0.001 1.81×10−9 (−2.52×10−8, 2.88×10−8) 0.896
    Vancomycin & norvancomycin 9.14×10−6 (6.38×10−6, 1.19×10−5) < 0.001 8.72×10−7 (−2.58×10−6, 4.33×10−6) 0.621
    Teicoplanin 3.94×10−7 (−7.47×10−6, 8.26×10−6) 0.922 −1.41×10−6 (−8.05×10−6, 5.24×10−6) 0.678
    Lincomycin−Clindamycin 6.69×10−7 (3.79×10−7, 9.59×10−7) < 0.001 −1.64×10−7 (−4.82×10−7, 1.55×10−7) 0.314
    Macrolides 7.33×10−8 (4.07×10−8, 1.06×10−7) < 0.001 −3.30×10−9 (−3.89×10−8, 3.23×10−8) 0.856
    Oxazolidinones 7.42×10−6 (5.03×10−6, 9.81×10−6) < 0.001 −2.06×10−6 (−5.66×10−6, 1.53×10−6) 0.261
    Penicillin 1.53×10−7 (1.02×10−7, 2.03×10−7) < 0.001 −2.10×10−8 (−8.23×10−8, 4.03×10−8) 0.503
    Polymyxin 3.41×10−5 (1.59×10−5, 5.23×10−5) < 0.001 −1.19×10−5 (−3.13×10−5, 7.59×10−6) 0.232
    Quinolones 1.40×10−7 (9.07×10−8, 1.90×10−7) < 0.001 −1.66×10−8 (−7.87×10−8, 4.55×10−8) 0.600
    Sulfonamides 1.88×10−6 (1.09×10−6, 2.66×10−6) < 0.001 1.61×10−7 (−6.58×10−7, 9.79×10−7) 0.701
    Tetracyclines 2.48×10−7 (1.37×10−7, 3.59×10−7) < 0.001 −6.69×10−8 (−1.65×10−7, 3.16×10−8) 0.183
      Note. VR E. faecium: Enterococcus faecium; VR E. faecalis: Enterococcus faecalis. 95% CI: 95% confidence interval.
    下载: 导出CSV

    Table  4.   Association analysis of VRE drug resistance with multifactor inclusion in hierarchical model in China from 2016 to 2022

    Variable VR E. faecium VR E. faecalis
    β (95% CI) p β (95% CI) p
    Year −0.31 (−0.37, −0.26) < 0.001 −0.20 (−0.27, −0.14) < 0.001
    Population −7.77×10−6 (−9.94×10−5, 8.39×10−5) 0.868 −5.48×10−5 (−1.39×10−4, 2.89×10−5) 0.199
    Aminoglycosides −4.88×10−7 (−1.82×10−6, 8.44×10−7) 0.473 7.89×10−8 (−1.37×10−6, 1.52×10−6) 0.915
    Amphenicols −1.34×10−5 (−2.66×10−5, −1.62×10−7) 0.047 1.48×10−5 (8.88×10−7, 2.87×10−5) 0.037
    β−lactamase inhibitors −1.23×10−7 (−4.00×10−7, 1.54×10−7) 0.383 2.70×10−8 (−2.41×10−7, 2.95×10−7) 0.844
    Carbapenems 6.26×10−7 (−5.15×10−7, 1.77×10−6) 0.282 −1.18×10−6 (−2.63×10−6, 2.62×10−7) 0.109
    Cephalosporins −1.03×10−9 (−9.58×10−8, 9.37×10−8) 0.983 −3.42×10−8 (−1.48×10−7, 7.91×10−8) 0.554
    Vancomycin & norvancomycin 1.87×10−5 (1.06×10−5, 2.69×10−5) < 0.001 1.36×10−5 (5.30×10−6, 2.20×10−5) 0.001
    Teicoplanin −9.71×10−6 (−1.89×10−5, −5.50×10−7) 0.038 7.62×10−7 (−9.30×10−6, 1.08×10−5) 0.882
    Lincomycin−Clindamycin 1.20×10−6 (6.82×10−7, 1.71×10−6) < 0.001 −7.13×10−8 (−7.11×10−7, 5.68×10−7) 0.827
    Macrolides −5.38×10−8 (−1.66×10−7, 5.79×10−8) 0.345 2.96×10−8 (−1.08×10−7, 1.67×10−7) 0.673
    Oxazolidinones 8.08×10−6 (2.58×10−6, 1.36×10−5) 0.004 3.85×10−6 (−3.43×10−6, 1.11×10−5) 0.300
    Penicillin 7.72×10−8 (−2.84×10−8, 1.83×10−7) 0.152 1.03×10−9 (−1.12×10−7, 1.14×10−7) 0.986
    Polymyxin −2.46×10−6 (−2.10×10−5, 1.61×10−5) 0.795 3.63×10−6 (−1.93×10−5, 2.66×10−5) 0.757
    Quinolones −3.53×10−7 (−5.62×10−7, −1.44×10−7) 0.001 −1.74×10−7 (−4.11×10−7, 6.22×10−8) 0.149
    Sulfonamides 7.85×10−7 (−2.10×10−7, 1.78×10−6) 0.122 4.91×10−7 (−6.32×10−7, 1.62×10−6) 0.391
    Tetracyclines −1.76×10−7 (−3.27×10−7, −2.54×10−8) 0.022 3.01×10−8 (−1.41×10−7, 2.01×10−7) 0.731
      Note. VR E. faecium: Enterococcus faecium; VR E. faecalis: Enterococcus faecalis. 95% CI: 95% confidence interval.
    下载: 导出CSV
  • [1] WHO. Antimicrobial resistance. 2023. https://www.who.int/news-room/fact-sheets/detail/antimicrobial-resistance. [2023-11-21]
    [2] Antimicrobial Resistance Collaborators. Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. Lancet, 2022; 399, 629−55. doi:  10.1016/S0140-6736(21)02724-0
    [3] Temkin E, Fallach N, Almagor J, et al. Estimating the number of infections caused by antibiotic-resistant Escherichia coli and Klebsiella pneumoniae in 2014: a modelling study. Lancet Glob Health, 2018; 6, e969−79. doi:  10.1016/S2214-109X(18)30278-X
    [4] Cimen C, Berends MS, Bathoorn E, et al. Vancomycin-resistant enterococci (VRE) in hospital settings across European borders: a scoping review comparing the epidemiology in the Netherlands and Germany. Antimicrob Resist Infect Control, 2023; 12, 78. doi:  10.1186/s13756-023-01278-0
    [5] Tamma PD, Aitken SL, Bonomo RA, et al. Infectious diseases society of America 2023 guidance on the treatment of antimicrobial resistant gram-negative infections. Clin Infect Dis, 2023; 18, ciad428.
    [6] Wang Q, Chen MY, Ou Q, et al. Carbapenem-resistant hypermucoviscous Klebsiella pneumoniae clinical isolates from a tertiary hospital in China: antimicrobial susceptibility, resistance phenotype, epidemiological characteristics, microbial virulence, and risk factors. Front Cell Infect Microbiol, 2022; 12, 1083009. doi:  10.3389/fcimb.2022.1083009
    [7] Huang L, Zhang R, Hu YY, et al. Epidemiology and risk factors of methicillin-resistant Staphylococcus aureus and vancomycin-resistant enterococci infections in Zhejiang China from 2015 to 2017. Antimicrob Resist Infect Control, 2019; 8, 90. doi:  10.1186/s13756-019-0539-x
    [8] Zhou WQ, Zhou H, Sun YH, et al. Characterization of clinical enterococci isolates, focusing on the vancomycin-resistant enterococci in a tertiary hospital in China: based on the data from 2013 to 2018. BMC Infect Dis, 2020; 20, 356. doi:  10.1186/s12879-020-05078-4
    [9] CDC. Antibiotic resistance threats in the United States, 2019. Atlanta, GA: U. S. Department of Health and Human Services, CDC. 2019.
    [10] Todorić Z, Majdandžić I, Keretić Kregar T, et al. Increasing trend in enterococcal bacteraemia and vancomycin resistance in a tertiary care hospital in Croatia, 2017-2021. Infect Dis, 2023; 55, 9−16. doi:  10.1080/23744235.2022.2131901
    [11] European Centre for Disease Prevention and Control. Antimicrobial resistance in the EU/EEA (EARS-Net) - Annual epidemiological report 2021. Stockholm: European Centre for Disease Prevention and Control. 2022. (
    [12] Roberts SC, Zembower TR. Global increases in antibiotic consumption: a concerning trend for WHO targets. Lancet Infect Dis, 2021; 21, 10−1. doi:  10.1016/S1473-3099(20)30456-4
    [13] Aslam B, Wang W, Arshad MI, et al. Antibiotic resistance: a rundown of a global crisis. Infect Drug Resist, 2018; 11, 1645−58. doi:  10.2147/IDR.S173867
    [14] Abejew AA, Wubetu GY, Fenta TG. Relationship between antibiotic consumption and resistance: a systematic review. Can J Infect Dis Med Microbiol, 2024; 2024, 9958678.
    [15] Gong W, Tang W, Luo L, et al. Trends and correlation between antimicrobial resistance and antibiotics consumption in a specialist children's hospital from 2016 to 2021. Infect Drug Resist, 2022; 15, 5679−89. doi:  10.2147/IDR.S381604
    [16] Pérez-Lazo G, Abarca-Salazar S, Lovón R, et al. Antibiotic consumption and its relationship with bacterial resistance profiles in ESKAPE pathogens in a Peruvian hospital. Antibiotics (Basel), 2021; 10, 1221. doi:  10.3390/antibiotics10101221
    [17] Boutzoukas AE, Komarow L, Chen L, et al. International epidemiology of carbapenemase-producing Escherichia coli. Clin Infect Dis, 2023; 77, 499−509. doi:  10.1093/cid/ciad288
    [18] Requirements for the quality and capacity of hospital microbiological laboratories in the network of national bacterial drug resistance monitoring network. http://www.carss.cn/Download/Details/656. [2020-01-03]. (In Chinese)
    [19] National surveillance report on bacterial resistance in 2022. https://www.carss.cn/Report/Details?aId=917. [2023-11-20]. (In Chinese)
    [20] Standard specification for data format of national bacterial resistance monitoring network. http://www.carss.cn/Download/Details/812. [2021-12-28]
    [21] Technical scheme of national bacterial resistance monitoring network (version 2020). http://www.carss.cn/Download/Details/657. [2020-01-03]. (in Chinese)
    [22] Chen X, Yu Y, Tang G. Pharamacology. Beijing: People's Medical Publishing House. 2011.
    [23] National Bureau of Statistics. China statistical yearbook 2023. China Statistics Press. 2023. (In Chinese)
    [24] Aliabadi S, Anyanwu P, Beech E, et al. Effect of antibiotic stewardship interventions in primary care on antimicrobial resistance of Escherichia coli bacteraemia in England (2013-18): a quasi-experimental, ecological, data linkage study. Lancet Infect Dis, 2021; 21, 1689−700. doi:  10.1016/S1473-3099(21)00069-4
    [25] López-Lozano JM, Lawes T, Nebot C, et al. A nonlinear time-series analysis approach to identify thresholds in associations between population antibiotic use and rates of resistance. Nat Microbiol, 2019; 4, 1160−72. doi:  10.1038/s41564-019-0410-0
    [26] Maugeri A, Barchitta M, Magnano San Lio R, et al. Socioeconomic and governance factors disentangle the relationship between temperature and antimicrobial resistance: a 10-year ecological analysis of European countries. Antibiotics (Basel), 2023; 12, 777. doi:  10.3390/antibiotics12040777
    [27] Zhang D, Hu SS, Sun JY, et al. Antibiotic consumption versus the prevalence of carbapenem-resistant Gram-negative bacteria at a tertiary hospital in China from 2011 to 2017. J Infect Public Health, 2019; 12, 195−9. doi:  10.1016/j.jiph.2018.10.003
    [28] Terahara F, Nishiura H. Carbapenem-resistant Pseudomonas aeruginosa and carbapenem use in Japan: an ecological study. J Int Med Res, 2019; 47, 4711−22. doi:  10.1177/0300060519864181
    [29] Botelho J, Grosso F, Peixe L. Antibiotic resistance in Pseudomonas aeruginosa - Mechanisms, epidemiology and evolution. Drug Resist Updat, 2019; 44, 100640. doi:  10.1016/j.drup.2019.07.002
    [30] Jung H, Pitout JDD, Matsumura Y, et al. Genomic epidemiology and molecular characteristics of blaNDM-1-positive carbapenem-resistant Pseudomonas aeruginosa belonging to international high-risk clone ST773 in the Gauteng region, South Africa. Eur J Clin Microbiol Infect Dis, 2024; 43, 627−40. doi:  10.1007/s10096-024-04763-5
    [31] Liu S, Li Y, He ZQ, et al. A molecular study regarding the spread of vanA vancomycin-resistant Enterococcus faecium in a tertiary hospital in China. J Glob Antimicrob Resist, 2022; 31, 270−8. doi:  10.1016/j.jgar.2022.10.010
    [32] Yan MY, He YH, Ruan GJ, et al. The prevalence and molecular epidemiology of vancomycin-resistant Enterococcus (VRE) carriage in patients admitted to intensive care units in Beijing, China. J Microbiol Immunol Infect, 2023; 56, 351−7. doi:  10.1016/j.jmii.2022.07.001
    [33] Zhao YC, Sun ZH, Li JK, et al. Exploring the causes of the prevalence of vancomycin-resistant Enterococcus faecalis. Environ Sci Eur, 2024; 36, 92. doi:  10.1186/s12302-024-00923-8
    [34] Ayobami O, Willrich N, Reuss A, et al. The ongoing challenge of vancomycin-resistant Enterococcus faecium and Enterococcus faecalis in Europe: an epidemiological analysis of bloodstream infections. Emerg Microbes Infect, 2020; 9, 1180−93. doi:  10.1080/22221751.2020.1769500
  • [1] Shengfen Wang, Xichao Ou, Yang Zhou, Bing Zhao, Hui Xia, Yuanyuan Song, Ruida Xing, Yang Zheng, Yanlin Zhao.  Drug-Resistant Pulmonary Tuberculosis Among Adolescents and Young Adults in China . Biomedical and Environmental Sciences, doi: 10.3967/bes2025.159
    [2] Lulu Zhang, Junping Peng.  High-Throughput Sequencing for Tuberculosis Diagnosis and Antimicrobial Resistance Detection: Progress, Challenges, and Future Perspectives . Biomedical and Environmental Sciences, doi: 10.3967/bes2025.164
    [3] DING Rui, MA Rui Rui, LIU Ya Li, ZHAO Ying, GUO Li Na, DOU Hong Tao, SUN Hong Li, LIU Wen Jing, ZHANG Li, WANG Yao, LI Ding Ding, YI Qiao Lian, XU Ying Chun.  Causative Microorganisms Isolated from Patients with Intra-Abdominal Infections and Their Drug Resistance Profiles: An 11-Year (2011–2021) Single-Center Retrospective Study . Biomedical and Environmental Sciences, doi: 10.3967/bes2023.072
    [4] HUANG Yin Yan, XIE Li, WU Yi Fei, JIA Qing Jun, CHENG Qing Lin, LI Qing Chun, AI Li Yun, BAI Xue Xin.  Characterization of Mutations in Genes Related to Rifampicin and Isoniazid Resistance in Multidrug-resistant Mycobacterium tuberculosis Strains from Hangzhou, China . Biomedical and Environmental Sciences, doi: 10.3967/bes2023.112
    [5] QU Feng Tian, WANG Wen Qing, LIU Qian, ZHOU Hai Jian, HU Jin Rui, DU Xiao Li, WANG Yue, XUE Jia Qi, CUI Zhi Gang, XIE Gui Lin, MENG Shuang.  Genetic Diversity, Antibiotic Resistance, and Pathogenicity of Aeromonas Species from Food Products in Shanghai, China . Biomedical and Environmental Sciences, doi: 10.3967/bes2022.108
    [6] MENG Shuang, DU Xiao Li, WANG Yong Lu, QU Feng Tian, XIE Gui Lin, ZHOU Hai Jian, HU Jin Rui, QIN Zheng, WANG Yue, KAN Biao, CUI Zhi Gang.  Comparative Study of the Genetic Diversity, Antimicrobial Resistance, and Pathogenicity of Aeromonas Isolates from Clinical Patients and Healthy Individuals . Biomedical and Environmental Sciences, doi: 10.3967/bes2021.062
    [7] LIU Chun Fa, SONG Yi Meng, HE Ping, LIU Dong Xin, HE Wen Cong, LI Yan Ming, ZHAO Yan Lin.  Evaluation of Multidrug Resistant Loop-mediated Isothermal Amplification Assay for Detecting the Drug Resistance of Mycobacterium tuberculosis . Biomedical and Environmental Sciences, doi: 10.3967/bes2021.085
    [8] MENG Shuang, WANG Yong Lu, LIU Chen Geng, YANG Jing, YUAN Min, BAI Xiang Ning, JIN Dong, LIANG Jun Rong, CUI Zhi Gang, LI Juan.  Genetic Diversity, Antimicrobial Resistance, and Virulence Genes of Aeromonas Isolates from Clinical Patients, Tap Water Systems, and Food . Biomedical and Environmental Sciences, doi: 10.3967/bes2020.053
    [9] CHE Jie, LU Jin Xing, LI Wen Ge, ZHANG Yun Fei, ZHAO Xiao Fei, YUAN Min, BAI Xue Mei, CHEN Xia, LI Juan.  A New High-throughput Real-time PCR Assay for the Screening of Multiple Antimicrobial Resistance Genes in Broiler Fecal Samples from China . Biomedical and Environmental Sciences, doi: 10.3967/bes2019.111
    [10] JU Chang Yan, ZHANG Mao Jun, MA Yan Ping, LU Jing Rang, YU Mu Hua, CHEN Hui, LIU Chu Yun, GU Yi Xin, FU Yan Yan, DUAN Yong Xiang.  Genetic and Antibiotic Resistance Characteristics of Campylobacter jejuni Isolated from Diarrheal Patients, Poultry and Cattle in Shenzhen . Biomedical and Environmental Sciences, doi: 10.3967/bes2018.079
    [11] CUI Jing Hua, YU Bo, XIANG Yun, ZHANG Zhen, ZHANG Ting, ZENG Ying Chun, CUI Zhi Gang, HUO Xi Xiang.  Two Cases of Multi-antibiotic Resistant Cronobacter spp.Infections of Infants in China . Biomedical and Environmental Sciences, doi: 10.3967/bes2017.079
    [12] WANG Sheng Fen, ZHOU Yang, PANG Yu, ZHENG Hui Wen, ZHAO Yan Lin.  Prevalence and Risk Factors of Primary Drug-Resistant Tuberculosis in China . Biomedical and Environmental Sciences, doi: 10.3967/bes2016.010
    [13] YAN Shao Fei, WANG Wei, BAI Li, HU Yu Jie, DONG Yin Ping, XU Jin, LI Feng Qin.  Antimicrobial Resistance, Virulence Profile, and Molecular Characterization of Listeria monocytogenes Isolated from Ready-to-eat Food in China, 2013-2014 . Biomedical and Environmental Sciences, doi: 10.3967/bes2016.058
    [14] Ashraf M.Ahmed, Akito Maruyama, Hazim O.Khalifa, Toshi Shimamoto, Tadashi Shimamoto.  Seafood as a Reservoir of Gram-negative Bacteria Carrying Integrons and Antimicrobial Resistance Genes in Japan . Biomedical and Environmental Sciences, doi: 10.3967/bes2015.128
    [15] HUANG Jing Jing, XI Jing Ying, HU Hong Ying, TANG Fang, PANG Yu Chen.  Inactivation and Regrowth of Antibiotic-resistant Bacteria by PAA Disinfection in the Secondary Effluent of a Municipal Wastewater Treatment Plant . Biomedical and Environmental Sciences, doi: 10.3967/bes2013.012
    [16] CHEN Yu Qiao, CHEN Chao, ZHANG Xiao Jian, ZHENG Qi, LIU Yuan Yuan.  Inactivation of Resistant Mycobacteria mucogenicum in Water: Chlorine Resistance and Mechanism Analysis* . Biomedical and Environmental Sciences, doi: 10.3967/0895-3988.2012.02.016
    [17] QIAO Rong, HUANG ChengYu, DU HuiZhang, ZENG Guo, LI Ling, YE Sheng.  Milk Consumption and Lactose Intolerance in Adults . Biomedical and Environmental Sciences, doi: 10.3967/0895-3988.2011.05.009
    [18] CHANG LIU, ZHUO-YANG ZHANG, KE DONG, JIAN-PING YUAN, XIAO-KUI GUO.  Antibiotic Resistance of Probiotic Strains of Lactic Acid Bacteria Isolated from Marketed Foods and Drugs . Biomedical and Environmental Sciences,
    [19] ZAHOOR QADIR SAMRA, MARIAM NASEEM, SUMARIA JAVED KHAN, NADIA DAR, MUHAMMAD AMIN ATHAR.  PCR Targeting of Antibiotic Resistant Bacteria in Public Drinking Water of Lahore Metropolitan, Pakistan . Biomedical and Environmental Sciences,
    [20] JAGDISH S PATEL, PRERNA C PATEL, KIRAN KALIA.  Isolation and Characterization of Nickel Uptake by Nickel Resistant Bacterial Isolate (NiRBI) . Biomedical and Environmental Sciences,
  • 加载中
图(4) / 表ll (4)
计量
  • 文章访问数:  32
  • HTML全文浏览量:  12
  • PDF下载量:  0
  • 被引次数: 0
出版历程
  • 收稿日期:  2025-06-30
  • 录用日期:  2025-10-17

Trends in Antibiotic Consumption and Antimicrobial Resistance in China: An Ecological Analysis from 2016 to 2022

doi: 10.3967/bes2025.169
    作者简介:

    Fei Zhao, Associate Researcher, PhD, Research Direction: Pharmacoepidemiology, Tel: 86-10-85133648, E-mail: zhaofei21cn@126.com

    Zinan Zhao, Chief Pharmacist, Master's Degree, Research Direction: Clinical Pharmacy, Tel: 86-10-85133648, E-mail: bellazhao2015@163.com

    通讯作者: Pengfei Jin, Prof, Tel: 86-10-85133620, E-mail: j790101@163.com
注释:
1) Authors’ Contributions: 2) Competing Interests: 3) Ethic:

English Abstract

Fei Zhao, Zinan Zhao, Di Chen, Bolin Zhu, Tianqi Zhang, Yuanchao Zhu, Xuelin Sun, Liang Liang, Nan Zheng, Lili Zou, Wenfeng Xu, Sirui Guo, Yue Wang, Ming Zhao, Xin Hu, Pengfei Jin. Trends in Antibiotic Consumption and Antimicrobial Resistance in China: An Ecological Analysis from 2016 to 2022[J]. Biomedical and Environmental Sciences. doi: 10.3967/bes2025.169
Citation: Fei Zhao, Zinan Zhao, Di Chen, Bolin Zhu, Tianqi Zhang, Yuanchao Zhu, Xuelin Sun, Liang Liang, Nan Zheng, Lili Zou, Wenfeng Xu, Sirui Guo, Yue Wang, Ming Zhao, Xin Hu, Pengfei Jin. Trends in Antibiotic Consumption and Antimicrobial Resistance in China: An Ecological Analysis from 2016 to 2022[J]. Biomedical and Environmental Sciences. doi: 10.3967/bes2025.169
    • Antimicrobial resistance (AMR) poses a major threat to human health globally[1]. AMR killed at least 1.27 million people worldwide in 2019 and may result in 10 million deaths annually by 2050 if the current trend in AMR is not suppressed[2]. Among the resistant pathogens, carbapenem-resistant organisms (CRO) and vancomycin-resistant enterococci (VRE) are of particular concern, especially in healthcare settings, where they can cause severe infections in people with weakened immune systems[3-5].

      Global control of these resistant bacteria is challenging. The epidemiology of AMR in China shows both similarities and differences with global patterns. According to nationwide bacterial resistance surveillance data in China, the resistance rates of carbapenem-resistant Klebsiella pneumoniae (CRKP) to meropenem (MEM) and imipenem (IPM) were 2.9% and 3.0% in 2005, respectively, but increased to 26.3% and 25.0%, respectively, by 2020, paralleling global trends[6]. In contrast, vancomycin resistance in Enterococcus faecium (E. faecium) and Enterococcus faecalis (E. faecalis) has remained at relatively low levels in China[7,8]. This differs from the United States and Europe. In a recent report on the threat of antibiotic resistance in the United States 2019, approximately 30% of all healthcare-associated enterococcal infections were resistant to vancomycin. Notably, more than 70% of E. faecium isolates were resistant to vancomycin in solid organ transplant units[9]. VRE resistance rates show a consistently increasing trend across Europe[4,10,11].

      Antimicrobial drugs reduce the global burden of common infectious diseases. Globally, an increase in the annual consumption of various types of antibiotics is apparent[12]. However, the widespread and indiscriminate use of antimicrobial drugs is a major driving force of AMR[13]. Several studies have shown that antibiotic drug use is closely associated with antibiotic resistance[14]. Most of these studies provided snapshots of the correlation between antibiotic consumption and AMR. However, the correlation between antibiotic consumption and AMR has not been consistent. Gong et al. reported a strong correlation between the consumption of carbapenems and carbapenem-resistant Klebsiella pneumoniae[15]. Conversely, Pérez-Lazo reported a strong negative correlation between imipenem and imipenem- and meropenem-resistant K. pneumoniae[16]. Additionally, in the context of the global increase in antibiotic consumption, the resistance rates of VRE in China have remained lower than those reported in other countries, whereas the resistance rates of CRO have shown an upward trend, similar to that observed in other countries[17]. Therefore, additional studies are required to provide comprehensive insights into the antimicrobial resistance of CRO and VRE with respect to antibiotic consumption in China.

      In this study, we examined the association between resistance to six WHO-designated priority pathogens and the consumption of 15 related classes of antibiotics to provide a clearer understanding of how improved medical practices and informed public health policies can help reduce the spread of resistant bacteria.

    • We conducted an ecological study of data pertaining to January 2016 and December 2022.

    • We collected the drug-resistant proportions of carbapenem-resistant organisms (CRO) and vancomycin-resistant enterococci (VRE) from 31 provinces from the China Antibacterial Resistance Surveillance System (CARSS), and at the national level on the Chinese mainland between 2016–2022. Resistant bacteria included carbapenem-resistant Escherichia coli (CREC), Acinetobacter baumannii (CRAB), Klebsiella pneumoniae (CRKP), Pseudomonas aeruginosa (CRPA), vancomycin-resistant Enterococcus faecium (VRE. faecium) and Enterococcus faecalis (VR E. faecalis)).

      The CARSS is a nationwide, voluntary bacterial resistance surveillance network covering secondary and tertiary hospitals that submits standardized laboratory results according to unified operating procedures and quality control requirements[18]. By 2022, 1997, hospitals had submitted drug resistance data to the CARSS, with 29% being secondary hospitals and 71% being tertiary hospitals[19]. Participation in the CARSS is dynamic, but all submitted data follow the same format and laboratory quality standards, which ensure comparability across sites[20]. Hierarchical information on bacterial drug resistance was obtained to understand the epidemic trend and new threats of bacterial drug resistance nationwide[21].

      Antibiotic consumption data were drawn from a subset of hospitals reported to the National Hospital Information Network of the Chinese Pharmaceutical Association. We selected 829 hospitals for antibiotic use analyses, choosing sites such that the distribution of hospital levels (secondary vs. tertiary) and provincial coverage broadly matched that of CARSS contributors. Hospitals with continuous reporting from 2016 to 2022 were prioritized to enhance temporal comparability. Consumption data comprised annual hospital-level totals for 161 antibiotics, which were extracted and classified into 15 pharmacological classes, including 1) Aminoglycosides, 2) Amphenicols, 3) Carbapenem, 4) Cephalosporins, 5) Lincomycin and Clindamycin, 6) Macrolides, 7) Oxazolidinones, 8) Penicillin, 9) Polymyxin, 10) Quinolone, 11) Sulfonamides, 12) Teicoplanin, 13) Tetracyclines, 14) Vancomycin and norvancomycin, and 15) β-Amidase inhibitor (including β-Amidase inhibitor, Cephalosporins/β-Amidase inhibitor, and Penicillin/β-Amidase inhibitor)[22] (Supplementary Table S1). The selection of antibiotics was guided by two criteria: their clinical relevance for treating infections caused by the six WHO-designated priority pathogens, and prior evidence of association with antimicrobial resistance in published literature (Supplementary Table S2). Antibiotic consumption was standardized to defined daily doses (DDDs) and reported annually as DDDs per patient-day at both provincial and national levels. We collected data on all antibacterial agents used systemically. Data regarding antifungals or antivirals were not considered.

      To ensure comparability between AMR and antibiotic consumption datasets, we harmonized both sources at the provincial–year level. Only hospitals that submitted data within a given year were included without the imputation of missing values. In addition, the population of different provinces was collected as a covariate variable from the China Statistics Yearbook from 2016 to 2022[23].

    • The median and interquartile range (25th and 75th percentiles) were reported for DDDs as continuous variables with a non-normal distribution, and the resistance proportions were summarized for each bacterial species. Temporal trends in nationwide resistance were analyzed using beta regression. The resistance proportion of the six bacteria was considered the dependent variable, whereas antibiotic consumption and provincial population by year were considered independent variables. To account for possible changes in hospital participation and clustering of data within provinces over seven consecutive years, we aggregated data at the provincial-year level and applied hierarchical models with province-level random effects. The association between resistance proportions as dependent and independent variables was examined using hierarchical models with a beta-link function. Before multivariate model fitting, we excluded variables with VIF > 10 to avoid multicollinearity. Continuous independent variables (e.g., year, population, and antibiotic consumption) were mean-centered prior to modeling to improve numerical stability and interpretability of regression coefficients. This approach reduces the bias from uneven hospital counts across provinces and partially mitigates effects of inconsistent hospital-level participation. Statistical significance was set than 0.05. All statistical analyses were performed using the R software (R Core Team 2024).

    • The funder played no role in the study design, data collection, data analysis, data interpretation, or manuscript writing. The corresponding author had full access to all data in the study and had the final responsibility for the decision to submit the manuscript for publication.

    • This study included 1,302 drug-resistant proportions at the provincial level, spanning seven years, across six bacterial strains in 31 provinces of mainland China. After removing the random effects of provinces, there was no statistically significant trend in the population of each province over time (β = 10.577, z = 1.765, P = 0.088).

    • The drug-resistant proportions of CRAB, CRPA, VR E. faecium, and VR E. faecalis at the national level showed a significant decline from 2016 to 2022 (P < 0.01, Figure 1, Supplementary Table S2). However, the drug-resistant proportions of CRKP exhibited a significant increase during the same period (β = 0.061, z = 4.442, P < 0.01, Figure 1). The provincial-level median resistance proportion of CRKP increased from 4.6% in 2016 to 8.2% in 2022, representing a nearly a two-fold increase The drug-resistant proportions of each bacterial strain varied significantly across the provinces in each respective year (P < 0.01, Figure 2).

      Figure 1.  Drug-resistant proportions of CRO and VRE in China from 2016 to 2022. Decreasing (CRAB and CRPA) and increasing (CRKP) trends.

      Figure 2.  Drug-resistant proportions of CRO and VRE across the provinces in China from 2016 to 2022. A: Drug resistance proportions of CRAB, CRPA, CRKP, and CREC. B: Proportions of drug-resistant VR E. faecium and VR E. faecalis. Interprovincial heterogeneity in the resistance proportions has been observed across China.

    • The annual DDDs sums for the 15 antibiotic classes nationwide are presented in Figure 3. After adjusting for the clustering effect of provinces, there was a statistically significant increase in DDDs for the following antibiotic classes from 2016 to 2022: 1) β-lactamase inhibitors, 2) Carbapenems, 3) Vancomycin and norvancomycin, 4) Lincomycin Clindamycin, 5) Oxazolidinones, 6) Penicillin, 7) Polymyxin, 8) Quinolones, 9) Sulfonamides, and 10) Tetracyclines (β > 0, P < 0.05, Supplementary Table S3). At the provincial level, the median DDDs of Penicillin increased from 550.87 × 105 in 2016 to 1490.21 × 105 in 2022, representing a 2.7-fold rise. The median DDDs of Polymyxin rose sharply from < 0.01 × 102 in 2016 to 4176.32 × 102 in 2022. Conversely, the median DDDs of Amphenicols at the provincial level declined substantially, from 1322.06×102 in 2016 to 155.48 × 102 in 2022 (β = −767.23, P = 0.019, Supplementary Table S3). Furthermore, for each antibiotic class, there were significant differences in DDDs among provinces (P < 0.01, Figure 4).

      Figure 3.  The total consumption of 15 antibiotic classes in China from 2016 to 2022.

      Figure 4.  The consumption of 15 antibiotic classes across provinces in China from 2016 to 2022.

    • After adjusting for clustering effects in provinces using a hierarchical model, the proportions of CRAB and CRPA decreased from 2016 to 2022 (P < 0.01). However, there was an increase in the proportion of CRKP (β = 0.06, P < 0.001). Meanwhile, no statistically significant change was observed in the proportion of CREC over time (β < 0.01, P = 0.896). See Table 1.

      Table 1.  Association analysis of CRO drug resistance with single−factor inclusion in hierarchical model in China from 2016 to 2022

      Variable CRAB CRPA CRKP CREC
      β (95% CI) p β (95% CI) p β (95% CI) p β (95% CI) p
      Year −0.04
      (−0.06, −0.02)
      <0.001 −0.05
      (−0.06, −0.04)
      <0.001 0.06
      (0.03, 0.09)
      <0.001 9.87×10−4
      (−0.01, 0.02)
      0.896
      Population 6.96×10−5
      (2.33×10−5, 11.59×10−5)
      0.003 4.52×10−5
      (1.05×10−6, 8.94×10−5)
      0.045 1.55×10−4
      (5.70×10−5, 2.53×10−4)
      0.002 7.83×10−5
      (2.64×10−5, 1.30×10−4)
      0.003
      Aminoglycosides 3.34×10−7
      (−2.35×10−7, 9.03×10−7)
      0.250 1.90×10−8
      (−3.67×10−7, 4.05×10−7)
      0.923 −6.71×10−7
      (−1.29×10−6, −5.60×10−8)
      0.033 1.54×10−7
      (−3.98×10−7, 7.06×10−7)
      0.584
      Amphenicols 2.92×10−6
      (−5.34×10−6, 1.12×10−5)
      0.489 5.51×10−6
      (4.51×10−7, 1.06×10−5)
      0.033 −7.23×10−6
      (−1.53×10−5, 8.70×10−7)
      0.080 −6.51×10−6
      (−1.31×10−5, 5.68×10−8)
      0.052
      β−lactamase inhibitors 2.63×10−8
      (−4.59×10−8, 9.85×10−8)
      0.476 −4.19×10−8
      (−9.28×10−8, 8.97×10−9)
      0.106 −8.63×10−8
      (−1.66×10−7, −6.45×10−9)
      0.034 4.15×10−9
      (−6.94×10−8, 7.77×10−8)
      0.912
      Carbapenems 1.07×10−7
      (−1.67×10−7, 3.82×10−7)
      0.444 −2.42×10−8
      (−2.32×10−7, 1.84×10−7)
      0.819 −2.34×10−7
      (−5.87×10−7, 1.18×10−7)
      0.192 2.02×10−7
      (−6.57×10−8, 4.69×10−7)
      0.139
      Cephalosporins 1.00×10−8
      (−4.93×10−9, 2.50×10−8)
      0.189 −1.37×10−9
      (−1.47×10−8, 1.20×10−8)
      0.841 −2.13×10−8
      (−4.16×10−8, −1.04×10−9)
      0.039 9.73×10−10
      (−1.72×10−8, 1.92×10−8)
      0.916
      Vancomycin & norvancomycin 1.24×10−7
      (−1.73×10−6, 1.92×10−6)
      0.896 −4.58×10−7
      (−1.92×10−6, 1.00×10−6)
      0.538 −1.71×10−6
      (−4.15×10−6, 7.42×10−7)
      0.172 1.55×10−6
      (−8.48×10−8, 3.18×10−6)
      0.063
      Teicoplanin 2.31×10−6
      (−1.86×10−6, 6.48×10−6)
      0.278 5.23×10−7
      (−2.50×10−6, 3.54×10−6)
      0.734 −2.21×10−6
      (−7.11×10−6, 2.69×10−6)
      0.376 3.21×10−6
      (−1.33×10−6, 7.75×10−6)
      0.166
      Lincomycin−Clindamycin 2.07×10−8
      (−1.43×10−7, 1.85×10−7)
      0.804 −2.76×10−8
      (−1.52×10−7, 9.63×10−8)
      0.662 1.70×10−8
      (−1.95×10−7, 2.29×10−7)
      0.875 7.49×10−8
      (−1.06×10−7, 2.55×10−7)
      0.416
      Macrolides 1.27×10−8
      (−7.27×10−9, 3.26×10−8)
      0.213 6.07×10−9
      (−8.00×10−9, 2.01×10−8)
      0.398 −3.07×10−8
      (−5.28×10−8, −8.57×10−9)
      0.007 −1.09×10−9
      (−2.45×10−8, 2.23×10−8)
      0.928
      Oxazolidinones −3.21×10−8
      (−9.66×10−7, 9.01×10−7)
      0.946 −7.48×10−7
      (−1.61×10−6, 1.10×10−7)
      0.088 6.44×10−8
      (−1.29×10−6, 1.41×10−6)
      0.925 2.09×10−7
      (−5.26×10−7, 9.44×10−7)
      0.577
      Penicillin 2.74×10−9
      (−2.17×10−8, 2.72×10−8)
      0.826 −1.46×10−8
      (−3.17×10−8, 2.60×10−9)
      0.096 −4.56×10−8
      (−6.96×10−8, −2.16×10−8)
      <0.001 7.94×10−10
      (−1.97×10−8, 2.13×10−8)
      0.940
      Polymyxin −1.67×10−6
      (−6.49×10−6, 3.15×106)
      0.498 −5.26×10−6
      (−9.97×10−6, −5.47×10−7)
      0.029 −6.42×10−6
      (−1.54×10−5, 2.55×10−6)
      0.162 −1.59×10−6
      (−4.98×10−6, 1.80×10−6)
      0.357
      Quinolones 8.73×10−9
      (−2.02×10−8, 3.77×10−8)
      0.554 −7.44×10−9
      (−2.87×10−8, 1.38×10−8)
      0.492 −2.51×10−8
      (−5.74×10−8, 7.26×10−9)
      0.129 8.73×10−9
      (−1.76×10−8, 3.51×10−8)
      0.516
      Sulfonamides −2.00×10−8
      (−4.31×10−7, 3.91×10−7)
      0.924 −6.80×10−8
      (−3.52×10−7, 2.16×10−7)
      0.638 −3.24×10−7
      (−8.61×10−7, 2.12×10−7)
      0.236 2.08×10−7
      (−2.26×10−7, 6.41×10−7)
      0.348
      Tetracyclines −5.21×10−9
      (−4.12×10−8, 3.08×10−8)
      0.777 −2.68×10−8
      (−6.28×10−8, 9.29×10−9)
      0.146 3.06×10−8
      (−1.53×10−8, 7.64×10−8)
      0.191 1.73×10−8
      (−1.22×10−8, 4.68×10−8)
      0.251
        Note. CRAB: Acinetobacter baumannii; CRPA: Pseudomonas aeruginosa; CRKP: Klebsiella pneumoniae; CREC: Escherichia coli. 95% CI: 95% confidence interval.

      When considering only the association between DDDs of a single drug class and resistance proportions, the DDDs of the carbapenem class did not correlate significantly with the resistance proportions of four CRO bacteria (P > 0.05, Table 1), while some other antibiotics had associations with the resistance proportions, including Polymyxin with CRPA, β-lactamase inhibitors, Macrolides, Penicillin, respectively, with CRKP (Table 1). However, when adjusted for all antibiotics in the hierarchical model, the DDDs of the carbapenem class were positively associated with the proportions of CRAB (β = 5.75×10−7, P = 0.007) and CRPA (β = 5.26×10−7, P = 0.001). Table 2.

      Table 2.  Association analysis of CRO drug resistance with multifactor inclusion in hierarchical model in China from 2016 to 2022

      Variabl×10 CRAB CRPA CRKP CR×10C
      β (95% CI) p β (95% CI) p β (95% CI) p β (95% CI) p
      Year −2.31×10−2 (−4.28×10−2, −3.32×10−3) 0.022 −5.70×10−2 (−7.35×10−2, −4.05×10−2) <0.001 8.12×10−2 (5.13×10−2, 1.11×10−1) <0.001 −4.81×10−3 (−3.08×10−2, 2.12×10−2) 0.717
      Population 8.34×10−5 (2.41×10−5, 1.43×10−4) 0.006 5.03×10−5 (5.05×10−6, 9.55×10−5) 0.029 1.20×10−4 (3.63×10−6, 2.37×10−4) 0.043 6.54×10−5 (1.03×10−5, 1.21×10−4) 0.020
      Aminoglycosides −2.29×10−7 (−8.23×10−7, 3.65×10−7) 0.449 −1.82×10−7 (−6.41×10−7, 2.77×10−7) 0.437 −1.94×10−7 (−9.81×10−7, 5.92×10−7) 0.628 5.41×10−7 (−1.44×10−7, 1.23×10−6) 0.122
      Amphenicols −1.37×10−6 (−6.68×10−6, 3.94×10−6) 0.613 2.03×10−7 (−4.00×10−6, 4.41×10−6) 0.924 2.66×10−6 (−4.30×10−6, 9.61×10−6) 0.454 −4.40×10−6 (−1.09×10−5, 2.06×10−6) 0.182
      β−lactamase inhibitors −4.96×10−8 (−1.63×10−7, 6.38×10−8) 0.391 −1.26×10−7 (−2.20×10−7, −3.16×10−8) 0.009 −1.65×10−8 (−1.74×10−7, 1.41×10−7) 0.838 −1.14×10−7 (−2.45×10−7, 1.75×10−8) 0.089
      Carbapenems 5.75×10−7 (1.60×10−7, 9.89×10−7) 0.007 5.26×10−7 (2.08×10−7, 8.45×10−7) 0.001 4.61×10−7 (−8.85×10−8, 1.01×10−6) 0.100 2.91×10−7 (−1.73×10−7, 7.55×10−7) 0.219
      Cephalosporins 3.09×10−8 (−9.83×10−9, 7.16×10−8) 0.137 −4.29×10−9 (−3.66×10−8, 2.80×10−8) 0.794 2.12×10−8 (−3.32×10−8, 7.56×10−8) 0.445 −3.11×10−8 (−8.08×10−8, 1.87×10−8) 0.221
      Vancomycin & norvancomycin −2.21×10−6 (−5.56×10−6, 1.15×10−6) 0.197 −6.03×10−7 (−3.28×10−6, 2.08×10−6) 0.660 1.61×10−6 (−2.94×10−6, 6.16×10−6) 0.488 3.64×10−6 (−1.46×10−7, 7.43×10−6) 0.060
      Teicoplanin 7.82×10−10 (−4.11×10−6, 4.11×10−6) 1.000 −1.34×10−6 (−4.56×10−6, 1.87×10−6) 0.413 1.80×10−7 (−5.21×10−6, 5.57×10−6) 0.948 −7.90×10−7 (−5.37×10−6, 3.79×10−6) 0.735
      Lincomycin−Clindamycin −7.22×10−8 (−2.89×10−7, 1.45×10−7) 0.514 −3.76×10−8 (−2.20×10−7, 1.45×10−7) 0.686 1.25×10−7 (−1.89×10−7, 4.39×10−7) 0.435 6.01×10−9 (−2.62×10−7, 2.74×10−7) 0.965
      Macrolides −1.90×10−8 (−6.32×10−8, 2.52×10−8) 0.400 1.15×10−8 (−2.41×10−8, 4.72×10−8) 0.527 −7.35×10−8 (−1.30×10−7, −1.68×10−8) 0.011 −1.68×10−8 (−7.07×10−8, 3.70×10−8) 0.540
      Oxazolidinones −1.28×10−6 (−3.40×10−6, 8.32×10−7) 0.234 −1.50×10−7 (−1.78×10−6, 1.48×10−6) 0.857 −4.01×10−7 (−2.83×10−6, 2.03×10−6) 0.747 7.93×10−7 (−1.56×10−6, 3.15×10−6) 0.509
      Penicillin −2.13×10−8 (−5.97×10−8, 1.72×10−8) 0.278 −3.57×10−8 (−6.75×10−8, −3.87×10−9) 0.028 −5.81×10−8 (−1.09×10−7, −6.90×10−9) 0.026 −2.50×10−8 (−7.17×10−8, 2.16×10−8) 0.293
      Polymyxin 2.45×10−6 (−3.91×10−6, 8.80×10−6) 0.451 3.84×10−6 (−1.29×10−6, 8.97×10−6) 0.142 −1.93×10−5 (−2.78×10−5, −1.08×10−5) <0.001 −1.13×10−5 (−1.92×10−5, −3.42×10−6) 0.005
      Quinolones 3.79×10−8 (−4.18×10−8, 1.18×10−7) 0.351 3.85×10−8 (−2.38×10−8, 1.01×10−7) 0.226 5.86×10−8 (−4.04×10−8, 1.57×10−7) 0.246 5.70×10−8 (−3.32×10−8, 1.47×10−7) 0.215
      Sulfonamides 1.76×10−8 (−3.67×10−7, 4.02×10−7) 0.929 2.51×10−7 (−8.00×10−8, 5.82×10−7) 0.137 −4.24×10−7 (−1.13×10−6, 2.86×10−7) 0.242 2.06×10−7 (−3.41×10−7, 7.52×10−7) 0.461
      Tetracyclines −4.39×10−8 (−1.03×10−7, 1.49×10−8) 0.143 −1.47×10−8 (−6.14×10−8, 3.21×10−8) 0.539 5.17×10−8 (−2.46×10−8, 1.28×10−7) 0.184 3.96×10−11 (−6.79×10−8, 6.80×10−8) 0.999
        Note. CRAB: Acinetobacter baumannii; CRPA: Pseudomonas aeruginosa; CRKP: Klebsiella pneumoniae; CREC: Escherichia coli. 95% CI: 95% confidence interval.

      There was a decline in the proportions of VR E. faecium (β = −0.13, P < 0.001, Table 3) and VR E. faecalis (β = −0.14, P < 0.001, Table 3) from 2016 to 2022. The proportion of VR E. faecium was positively correlated with the DDDs of vancomycin and norvancomycin, irrespective of the usage of other antibiotics (β > 0, P < 0.001, Table 3 and Table 4). However, with consideration for the usage of other antibiotics, the results indicated a positive correlation between the DDDs of vancomycin and norvancomycin and the proportion of VR E. faecalis (β = 1.36×10−5, P = 0.001, Table 4), despite the non-significant correlation when only considering the association between the DDDs of vancomycin and norvancomycin and resistance proportions (β = 8.72×10−7, P = 0.621, Table 3).

      Table 3.  Association analysis of VRE drug resistance with single−factor inclusion in hierarchical model in China from 2016 to 2022

      Variable VR E. faecium VR E. faecalis
      β (95% CI) p β (95% CI) p
      Year −0.13 (−0.19, −0.06) < 0.001 −0.14 (−0.19, −0.10) < 0.001
      Population 1.47×10−5 (−7.00×10−5, 9.94×10−5) 0.733 −4.39×10−5 (−1.05×10−4, 1.73×10−5) 0.160
      Aminoglycosides 2.18×10−6 (1.14×10−6, 3.23×10−6) < 0.001 −3.34×10−8 (−1.03×10−6, 9.60×10−7) 0.947
      Amphenicols 1.72×10−5 (1.05×10−6, 3.33×10−5) 0.037 1.65×10−5 (1.70×10−6, 3.14×10−5) 0.029
      β−lactamase inhibitors 3.39×10−7 (1.96×10−7, 4.83×10−7) < 0.001 −5.53×10−8 (−1.99×10−7, 8.82×10−8) 0.450
      Carbapenems 1.42×10−6 (9.41×10−7, 1.90×10−6) < 0.001 −1.81×10−7 (−7.24×10−7, 3.62×10−7) 0.514
      Cephalosporins 7.50×10−8 (5.00×10−8, 9.99×10−8) < 0.001 1.81×10−9 (−2.52×10−8, 2.88×10−8) 0.896
      Vancomycin & norvancomycin 9.14×10−6 (6.38×10−6, 1.19×10−5) < 0.001 8.72×10−7 (−2.58×10−6, 4.33×10−6) 0.621
      Teicoplanin 3.94×10−7 (−7.47×10−6, 8.26×10−6) 0.922 −1.41×10−6 (−8.05×10−6, 5.24×10−6) 0.678
      Lincomycin−Clindamycin 6.69×10−7 (3.79×10−7, 9.59×10−7) < 0.001 −1.64×10−7 (−4.82×10−7, 1.55×10−7) 0.314
      Macrolides 7.33×10−8 (4.07×10−8, 1.06×10−7) < 0.001 −3.30×10−9 (−3.89×10−8, 3.23×10−8) 0.856
      Oxazolidinones 7.42×10−6 (5.03×10−6, 9.81×10−6) < 0.001 −2.06×10−6 (−5.66×10−6, 1.53×10−6) 0.261
      Penicillin 1.53×10−7 (1.02×10−7, 2.03×10−7) < 0.001 −2.10×10−8 (−8.23×10−8, 4.03×10−8) 0.503
      Polymyxin 3.41×10−5 (1.59×10−5, 5.23×10−5) < 0.001 −1.19×10−5 (−3.13×10−5, 7.59×10−6) 0.232
      Quinolones 1.40×10−7 (9.07×10−8, 1.90×10−7) < 0.001 −1.66×10−8 (−7.87×10−8, 4.55×10−8) 0.600
      Sulfonamides 1.88×10−6 (1.09×10−6, 2.66×10−6) < 0.001 1.61×10−7 (−6.58×10−7, 9.79×10−7) 0.701
      Tetracyclines 2.48×10−7 (1.37×10−7, 3.59×10−7) < 0.001 −6.69×10−8 (−1.65×10−7, 3.16×10−8) 0.183
        Note. VR E. faecium: Enterococcus faecium; VR E. faecalis: Enterococcus faecalis. 95% CI: 95% confidence interval.

      Table 4.  Association analysis of VRE drug resistance with multifactor inclusion in hierarchical model in China from 2016 to 2022

      Variable VR E. faecium VR E. faecalis
      β (95% CI) p β (95% CI) p
      Year −0.31 (−0.37, −0.26) < 0.001 −0.20 (−0.27, −0.14) < 0.001
      Population −7.77×10−6 (−9.94×10−5, 8.39×10−5) 0.868 −5.48×10−5 (−1.39×10−4, 2.89×10−5) 0.199
      Aminoglycosides −4.88×10−7 (−1.82×10−6, 8.44×10−7) 0.473 7.89×10−8 (−1.37×10−6, 1.52×10−6) 0.915
      Amphenicols −1.34×10−5 (−2.66×10−5, −1.62×10−7) 0.047 1.48×10−5 (8.88×10−7, 2.87×10−5) 0.037
      β−lactamase inhibitors −1.23×10−7 (−4.00×10−7, 1.54×10−7) 0.383 2.70×10−8 (−2.41×10−7, 2.95×10−7) 0.844
      Carbapenems 6.26×10−7 (−5.15×10−7, 1.77×10−6) 0.282 −1.18×10−6 (−2.63×10−6, 2.62×10−7) 0.109
      Cephalosporins −1.03×10−9 (−9.58×10−8, 9.37×10−8) 0.983 −3.42×10−8 (−1.48×10−7, 7.91×10−8) 0.554
      Vancomycin & norvancomycin 1.87×10−5 (1.06×10−5, 2.69×10−5) < 0.001 1.36×10−5 (5.30×10−6, 2.20×10−5) 0.001
      Teicoplanin −9.71×10−6 (−1.89×10−5, −5.50×10−7) 0.038 7.62×10−7 (−9.30×10−6, 1.08×10−5) 0.882
      Lincomycin−Clindamycin 1.20×10−6 (6.82×10−7, 1.71×10−6) < 0.001 −7.13×10−8 (−7.11×10−7, 5.68×10−7) 0.827
      Macrolides −5.38×10−8 (−1.66×10−7, 5.79×10−8) 0.345 2.96×10−8 (−1.08×10−7, 1.67×10−7) 0.673
      Oxazolidinones 8.08×10−6 (2.58×10−6, 1.36×10−5) 0.004 3.85×10−6 (−3.43×10−6, 1.11×10−5) 0.300
      Penicillin 7.72×10−8 (−2.84×10−8, 1.83×10−7) 0.152 1.03×10−9 (−1.12×10−7, 1.14×10−7) 0.986
      Polymyxin −2.46×10−6 (−2.10×10−5, 1.61×10−5) 0.795 3.63×10−6 (−1.93×10−5, 2.66×10−5) 0.757
      Quinolones −3.53×10−7 (−5.62×10−7, −1.44×10−7) 0.001 −1.74×10−7 (−4.11×10−7, 6.22×10−8) 0.149
      Sulfonamides 7.85×10−7 (−2.10×10−7, 1.78×10−6) 0.122 4.91×10−7 (−6.32×10−7, 1.62×10−6) 0.391
      Tetracyclines −1.76×10−7 (−3.27×10−7, −2.54×10−8) 0.022 3.01×10−8 (−1.41×10−7, 2.01×10−7) 0.731
        Note. VR E. faecium: Enterococcus faecium; VR E. faecalis: Enterococcus faecalis. 95% CI: 95% confidence interval.
    • This nationwide ecological study provided an overview of trends in carbapenem-resistant organisms (CRO) and vancomycin-resistant enterococci (VRE), along with antibiotic consumption in mainland China from 2016 to 2022. Although the observed declines in CRAB, CRPA, and VRE were encouraging, the continued increase in CRKP underscores the persistent threat of carbapenem resistance. Furthermore, significant changes in the DDDs of various antibiotic classes were observed from 2016 to 2022, providing evidence for the global trend of antibiotic overuse.

      The association between drug resistance and antibiotic use was also investigated. Several studies have indicated a strong correlation between antibiotic drug use and antibiotic resistance[24]. Lopez-Lozano et al. found empirical evidence of nonlinear relationships between population antibiotic consumption and resistance rates in five European settings[25]. Abejew et al. systematically reviewed the correlation between antibiotic resistance and global antibiotic consumption globally[14]. These results inspired the present study and provided a theoretical basis. Maugeri et al. found that antibiotic consumption is one of the main predictors of antimicrobial resistance in European countries[26]. Similarly, in our study, the consumption of carbapenem antibiotics was a predictor of CRAB, CRPA, and CRKP resistance, whereas the use of Vancomycin and Norvancomycin predicted VR E. faecium and VR E. faecalis resistance, consistent with most global research results[27,28].

      Many studies have shown a comprehensive association between antibiotic consumption and resistance to CRPA and CRKP due to their complex drug resistance mechanisms and cross-resistance[29,30]. Similarly, in our study, exploring the correlation between a single antibiotic and the resistance proportions of CRO and VRE was insufficient. We did not find an association between the consumption of carbapenem antibiotics and resistance to CRAB, CRPA, CRKP, and CREC, nor between the consumption of vancomycin and vancomycin and resistance to VR E. faecalis when considering only the association with a single drug class. However, this association was evident after adjusting for consumption of all other antibiotics. Therefore, it is necessary to consider the effects of other antibiotics on the development of drug resistance. Furthermore, the use of polymyxin was negatively associated with the proportions of CRKP and CREC, and penicillin usage showed a negative correlation with the resistance of CRPA and CRKP. Similar results have been previously reported[14]. These findings underscore the complex relationship between antibiotic use and emergence of drug-resistant bacterial strains, highlighting the critical need for judicious antibiotic management.

      Previous studies and our study indicate that the resistance rate to VRE in China differs from that in Europe and the United States, remaining low[31]. Although the use of various antibiotics has significantly increased in China, there is a clear association between vancomycin and VRE consumption. Therefore, compared with the high drug resistance of VRE in Europe and the United States, it is necessary to strengthen surveillance by conducting more sampling in China, especially in high-risk hospital units, such as hematology, ICU, and renal units[32-34]. However, the observed changes in statistical significance between the unadjusted and adjusted models likely reflect confounding from concurrent antibiotic consumption and the complex interplay among antibiotic classes. Although our ecological design could not confirm causal interactions, these findings are consistent with the possibility of synergistic or indirect effects, particularly in the case of vancomycin use and VRE resistance.

      Our findings have important implications for antimicrobial stewardship and public health policies. Continued efforts are needed to sustain the downward trends in CRAB, CRPA, and VRE, likely supported by strengthened infection prevention and rational prescriptions. Simultaneously, targeted interventions, including stricter carbapenem stewardship, improved diagnostic capacity, and expanded access to effective alternative therapies, are urgently required to curb the alarming increase in the incidence of CRKP. Finally, pronounced interprovincial variability highlights the need for region-specific strategies tailored to local epidemiology and resource capacities. Collectively, these measures could contribute to a more comprehensive national response to antimicrobial resistance in China.

      Compared with previous studies, our analysis presents novelties that also represent some of its strengths. First, we evaluated the association between antibiotic consumption and AMR at the provincial level in mainland China. These results provide antibiotics as indicators of better AMR control. Second, our study included a majority of drug-resistant CRO and VRE bacteria. These results provide a comprehensive overview of bacterial resistance and antibiotic use in mainland China. Finally, the introduction of a hierarchical model enabled us to disentangle the relationship between antibiotic consumption and AMR.

      This study has some limitations. First, as an ecological study, we could not infer causality from observed relationships. Based on previous studies and biological plausibility, we selected a set of antibiotics that may be effective in treating bacterial infections that are likely to induce drug resistance. However, these findings should be interpreted with caution, as the ecological design precludes causal inference and cannot disentangle the influence of concurrent infection control measures, regional prescribing practices, or socioeconomic determinants. Second, differences in reporting practices between the two surveillance networks and incomplete hospital participation across years may have introduced bias, although our harmonization procedures and modeling approach helped minimize these effects. Third, all data were available at provincial and yearly levels. Therefore, it was impossible to evaluate the observed relationships within a smaller geographic unit or season. However, this remains a direction for future research. Future research should account for additional confounders and consider use of spatial statistical approaches to explore potential geographic clustering in antimicrobial resistance and antibiotic consumption.

    • The proportions of CRAB, CRPA, VR E. faecium, and VR E. faecalis declined across the study period, whereas the proportion of CRKP increased in mainland China. The DDDs of 10 antibiotic classes exhibited a statistically significant nationwide increase. Our study corroborates previous findings, suggesting that the consumption of carbapenems could be the main indicator of carbapenem-resistant strains of Acinetobacter baumannii and Pseudomonas aeruginosa, whereas the consumption of vancomycin serves as a key indicator of VRE. However, it is crucial to consider the effect of other antibiotics on drug resistance. Further studies at a more granular level are necessary to thoroughly explore the relationship between AMR and antibiotic use.

参考文献 (34)

目录

    /

    返回文章
    返回