Data sources
A comprehensive literature search was conducted across both international and domestic databases to ensure global coverage of AMR research related to China. For international literature, we systematically searched the Web of Science Core Collection (including Science Citation Index Expanded, SCI-E) and the China National Knowledge Infrastructure (CNKI) from January 1, 2016, to December 31, 2025. Exact Boolean Query for SCI-E: TS = ((“antimicrobial resistance” OR “antibiotic resistance” OR “drug resistance”) AND (“China” OR “Chinese” OR “Mainland China”)). For CNKI, the Subject field was searched using the Chinese-language equivalents of “antibiotic resistance,” “antimicrobial resistance,” and “drug resistance,” combined with “China.” No formal keyword-synonym merging or thesaurus-based preprocessing was performed before the CiteSpace analysis. The final search was executed on April 28, 2026, to allow sufficient time for publications issued through December 31, 2025, to be indexed in the databases. The database publication-year filter was restricted to 2016–2025. Records labeled as Early Access were additionally checked during screening, and articles with a first online publication date in 2026 were excluded from the analysis. Keyword labels were checked against the original CiteSpace output to ensure consistent terminology in the text, tables, and figures.
For the CNKI search, 983 records were retrieved from the Peking University Core Journals (PKU Core), 802 from the Chinese Science Citation Database (CSCD), and 462 from the World Journal Clout Index (WJCI). Because a single article may be indexed in more than one of these categories, these numbers represent category-specific retrieval records and are not mutually exclusive publication counts. The three record sets were therefore merged before eligibility screening. Duplicate entries were identified first by exact DOI matching. For records without a DOI, normalized article titles were compared together with the first-author name and publication year. Records with incomplete or inconsistent metadata were reviewed manually. When the same article appeared in two or more CNKI indexing categories, it was retained only once in the combined CNKI dataset, irrespective of the number of categories in which it was indexed. After within-CNKI deduplication and eligibility screening, 1,096 unique Chinese-language articles were retained. No journal impact-factor threshold was applied, because the purpose of this study was to characterize the overall AMR research landscape rather than to restrict the corpus according to journal prestige. Formal methodological-quality or risk-of-bias assessment of individual publications was not performed because the objective of this study was to characterize the overall bibliometric landscape rather than to evaluate study-level methodological quality.
Across SCI-E and CNKI, 5,286 records were identified initially. The same DOI-, title-, author-, and year-based criteria were subsequently applied to identify potential overlaps between the CNKI and Web of Science datasets. A total of 2,190 records were removed during combined deduplication and eligibility screening, leaving 3,096 publications for analysis. The complete literature identification, screening, and inclusion process is presented in Supplementary Figure 1.
Inclusion criteria were: studies focusing on bacterial AMR in China; articles reporting resistance rates, molecular epidemiology, clinical outcomes, or intervention strategies; and publications between January 1, 2016, and December 31, 2025. Duplicate records, 2026 publications, reviews, or meta-analyses without original datasets, and studies limited to veterinary or environmental AMR without a human-health component were excluded. Reviews and meta-analyses without new primary datasets were excluded because the analysis was designed to characterize original research output. Restricting the corpus to original research improves the comparability of bibliometric indicators, including publication counts, co-authorship patterns, keyword frequencies, and citation-related measures, because evidence-synthesis articles aggregate prior studies and may have substantially different citation and keyword profiles. After removing duplicates across databases and screening titles/abstracts, a total of 3,096 publications were included for bibliometric analysis, comprising 1,096 Chinese-language articles (35.4%) and 2000 English-language publications (64.6%).
Analytical methods
CiteSpace (Drexel University, Philadelphia, PA, USA)17 was used for bibliometric analysis. The time span was 2016–2025 with a 1-year slice length. Node types were selected separately for author, institution, and keyword analyses. Node selection used the g-index with k = 10; the displayed configuration also used LRF = 2.5, L/N = 10, LBY = 5, and e = 1.0. Pathfinder pruning was applied to the author, keyword, and cluster networks to reduce redundant links while retaining the principal network structure. The 1-year time slice was selected to preserve annual changes within the 10-year study window.
Four main analytical modules were utilized: co-authorship analysis to map collaboration networks among authors and institutions; keyword co-occurrence analysis to identify research themes and their relationships; cluster analysis using the Log-Likelihood Ratio (LLR) algorithm to extract thematic clusters with automatic labeling; and burst detection using Kleinberg’s algorithm to identify emerging trends and research frontiers. Network metrics, including modularity Q and mean silhouette S, were used to evaluate cluster quality. A modularity Q value greater than 0.3 indicates a meaningful community structure, whereas a mean silhouette value greater than 0.7 indicates high within-cluster consistency and reliable cluster assignment18.
Statistical analysis of resistance trends
Resistance rate trends for major pathogens were derived from the China Antimicrobial Resistance Surveillance System (CARSS) annual reports (2016–2025), which aggregate data from participating sentinel hospitals nationwide. These surveillance data were analyzed independently of the bibliometric corpus and were not extracted from, pooled across, or included among the 3,096 Web of Science and CNKI publications used for the CiteSpace analyses.
For each pathogen–antimicrobial combination, we extracted the annual resistance rate (percentage of non-duplicate isolates resistant to the specified agent), the total number of isolates tested (denominator), and the number of sentinel hospitals contributing data. Susceptibility testing was performed using standardized broth microdilution or disk diffusion methods, with interpretations based on Clinical and Laboratory Standards Institute (CLSI) guidelines (M100 documents, editions 2016–2025).
To assess temporal trends, we performed linear regression analysis with the resistance rate as the dependent variable and year as the independent variable, covering the full 10-year period (2016–2025). The regression slope (β) and coefficient of determination (R2) were calculated for each pathogen–antimicrobial pair. To account for potential changes in surveillance coverage (e.g., expansion from 1,398 hospitals in 2016 to 2,349 hospitals in 2025), we also performed sensitivity analyses weighting each year's resistance rate by the inverse of the number of participating hospitals and by the geographic distribution of reporting sites. These weighted models yielded trend directions and significance levels consistent with the unweighted analyses; therefore, unweighted results are presented for simplicity and interpretability. All statistical analyses were performed using R version 4.3.2 (R Foundation for Statistical Computing, Vienna, Austria). A two-sided p-value <0.05 was considered statistically significant.