Źródła danych
Przeprowadzono kompleksowy przegląd literatury w międzynarodowych i krajowych bazach danych, aby zapewnić globalny zasięg badań nad AMR w odniesieniu do Chin. W przypadku literatury międzynarodowej systematycznie przeszukiwano Web of Science Core Collection (w tym Science Citation Index Expanded, SCI-E) oraz China National Knowledge Infrastructure (CNKI) w okresie od 1 stycznia 2016 r. do 31 grudnia 2025 r. Dokładne zapytanie Boole'a dla SCI-E: TS = ((“antimicrobial resistance” OR “antibiotic resistance” OR “drug resistance”) AND (“China” OR “Chinese” OR “Mainland China”)). W przypadku CNKI przeszukiwano pole Temat (Subject) przy użyciu chińskojęzycznych odpowiedników terminów „oporność na antybiotyki”, „oporność na środki przeciwdrobnoustrojowe” oraz „oporność na leki” w połączeniu z terminem „Chiny”. Przed analizą w programie CiteSpace nie przeprowadzono formalnego łączenia synonimów słów kluczowych ani wstępnego przetwarzania w oparciu o tezaurus. Ostateczne wyszukiwanie wykonano 28 kwietnia 2026 r., aby zapewnić wystarczającą ilość czasu na zaindeksowanie w bazach danych publikacji wydanych do 31 grudnia 2025 r. Filtr roku publikacji w bazie danych ograniczono do lat 2016–2025. Podczas przesiewu dodatkowo sprawdzono rekordy oznaczone jako Early Access, a z analizy wykluczono artykuły, których data pierwszej publikacji online przypadała na 2026 r. Etykiety słów kluczowych zweryfikowano z oryginalnymi wynikami CiteSpace, aby zapewnić spójność terminologiczną w tekście, tabelach i rysunkach.
W wyszukiwaniu w bazie CNKI pobrano 983 rekordy z Peking University Core Journals (PKU Core), 802 z Chinese Science Citation Database (CSCD) oraz 462 z World Journal Clout Index (WJCI). Ponieważ pojedynczy artykuł może być zaindeksowany w więcej niż jednej z tych kategorii, liczby te reprezentują rekordy wyszukiwania specyficzne dla danej kategorii i nie są wzajemnie wykluczającymi się liczbami publikacji. W związku z tym trzy zestawy rekordów zostały połączone przed przesiewem kwalifikowalności. Powtórzone wpisy zidentyfikowano w pierwszej kolejności poprzez dokładne dopasowanie DOI. W przypadku rekordów bez numeru DOI porównano znormalizowane tytuły artykułów wraz z nazwiskiem pierwszego autora i rokiem publikacji. Rekordy z niepełnymi lub niespójnymi metadanymi sprawdzono ręcznie. Gdy ten sam artykuł pojawiał się w dwóch lub więcej kategoriach indeksowania CNKI, był on zachowywany tylko raz w połączonym zbiorze danych CNKI, niezależnie od liczby kategorii, w których został zaindeksowany. Po usunięciu duplikatów w obrębie bazy CNKI i przesiewie kwalifikowalności zachowano 1 096 unikalnych artykułów w języku chińskim. Nie zastosowano progu współczynnika wpływu (impact factor) czasopism, ponieważ celem niniejszego badania była charakterystyka ogólnego krajobrazu badań nad AMR, a nie ograniczenie korpusu ze względu na prestiż czasopisma. Nie przeprowadzono formalnej oceny jakości metodologicznej ani analizy ryzyka błędu (risk of bias) poszczególnych publikacji, ponieważ celem niniejszego badania była charakterystyka ogólnego krajobrazu bibliometrycznego, a nie ocena jakości metodologicznej na poziomie poszczególnych badań.
W bazach SCI-E oraz CNKI początkowo zidentyfikowano 5 286 rekordów. Następnie zastosowano te same kryteria oparte na DOI, tytule, autorze i roku w celu wykrycia potencjalnych powtórzeń między zbiorami danych CNKI i Web of Science. W procesie łącznej deduplikacji oraz oceny kwalifikowalności usunięto łącznie 2 190 rekordów, co pozostawiło 3 096 publikacji do analizy. Cały proces identyfikacji, przesiewania i włączania literatury przedstawiono na Rysunku uzupełniającym 1.
Kryteria włączenia obejmowały: badania koncentrujące się na bakteryjnej AMR w Chinach; artykuły raportujące wskaźniki oporności, epidemiologię molekularną, wyniki kliniczne lub strategie interwencyjne oraz publikacje wydane między 1 stycznia 2016 r. a 31 grudnia 2025 r. Wykluczono duplikaty, publikacje z 2026 r., przeglądy lub metaanalizy bez oryginalnych zbiorów danych oraz badania ograniczone do AMR w weterynarii lub środowisku bez komponentu dotyczącego zdrowia ludzi. Przeglądy i metaanalizy bez nowych pierwotnych zbiorów danych zostały wykluczone, ponieważ analiza została zaprojektowana w celu charakterystyki wyników badań oryginalnych. Ograniczenie korpusu do badań oryginalnych poprawia porównywalność wskaźników bibliometrycznych, w tym liczby publikacji, wzorców współautorstwa, częstotliwości słów kluczowych i miar związanych z cytowaniami, ponieważ artykuły syntetyzujące dowody agregują wcześniejsze badania i mogą mieć znacząco inne profile cytowań i słów kluczowych. Po usunięciu duplikatów z baz danych i przesiewowym przeglądzie tytułów/streszczeń, do analizy bibliometrycznej włączono łącznie 3 096 publikacji, w tym 1 096 artykułów w języku chińskim (35,4%) i 2000 publikacji w języku angielskim (64,6%).
Metody analityczne
Do analizy bibliometrycznej wykorzystano oprogramowanie CiteSpace (Drexel University, Philadelphia, PA, USA)17. Przedział czasowy obejmował lata 2016–2025 przy długości wycinka czasowego wynoszącej 1 rok. Typy węzłów wybrano oddzielnie dla analiz autorów, instytucji oraz słów kluczowych. Do selekcji węzłów wykorzystano wskaźnik g-index z parametrem k = 10; w wyświetlonej konfiguracji zastosowano również wartości LRF = 2.5, L/N = 10, LBY = 5 oraz e = 1.0. W celu redukcji nadmiarowych połączeń przy jednoczesnym zachowaniu głównej struktury sieci, w sieciach autorów, słów kluczowych i klastrów zastosowano przycinanie typu Pathfinder. Roczny wycinek czasowy wybrano, aby zachować coroczne zmiany w ramach 10-letniego okna badania.
Wykorzystano cztery główne moduły analityczne: analizę współautorstwa w celu zmapowania sieci współpracy między autorami i instytucjami; analizę współwystępowania słów kluczowych w celu zidentyfikowania tematów badawczych i ich wzajemnych relacji; analizę skupień z wykorzystaniem algorytmu Log-Likelihood Ratio (LLR) w celu wyodrębnienia klastrów tematycznych z automatycznym etykietowaniem; oraz detekcję gwałtownych wzrostów (burst detection) z wykorzystaniem algorytmu Kleinberga w celu zidentyfikowania wyłaniających się trendów i nowych obszarów badań. Do oceny jakości klastrów wykorzystano metryki sieciowe, w tym modularność Q oraz średnią sylwetkę S. Wartość modularności Q większa niż 0,3 wskazuje na znaczącą strukturę wspólnoty, natomiast średnia wartość sylwetki większa niż 0,7 świadczy o wysokiej spójności wewnątrz klastra i niezawodności przypisania do klastra18.
Analiza statystyczna trendów oporności
Trendy w zakresie wskaźników oporności głównych patogenów zostały opracowane na podstawie raportów rocznych Chińskiego Systemu Nadzoru nad Opornością na Środki Przeciwmikrobowe (CARSS) z lat 2016–2025, które agregują dane z uczestniczących szpitali sentinelowych w całym kraju. Dane nadzorcze te analizowano niezależnie od korpusu bibliometrycznego i nie zostały one wyodrębnione, połączone ani włączone do zbioru 3 096 publikacji z baz Web of Science i CNKI wykorzystanych w analizach CiteSpace.
Dla każdej kombinacji patogen–środek przeciwdrobnoustrojowy wyodrębniono roczną stopę oporności (procent niepowtarzających się izolatów opornych na określony preparat), całkowitą liczbę przebadanych izolatów (mianownik) oraz liczbę szpitali sentinelowych dostarczających dane. Badania wrażliwości przeprowadzono przy użyciu standaryzowanych metod mikrorozcieńczeń w bulionie lub dyfuzyjnych metod krążkowych, a interpretację wyników oparto na wytycznych Clinical and Laboratory Standards Institute (CLSI) (dokumenty M100, wydania 2016–2025).
Aby ocenić trendy czasowe, przeprowadzono analizę regresji liniowej, w której zmienną zależną była stopa oporności, a zmienną niezależną rok, obejmującą pełny 10-letni okres (2016–2025). Dla każdej pary patogen–lek przeciwdrobnoustrojowy obliczono nachylenie regresji (β) oraz współczynnik determinacji (R2). Aby uwzględnić potencjalne zmiany w zakresie nadzoru (np. rozszerzenie z 1398 szpitali w 2016 r. do 2349 szpitali w 2025 r.), przeprowadzono również analizy wrażliwości, ważąc stopę oporności z każdego roku odwrotnością liczby uczestniczących szpitali oraz rozkładem geograficznym ośrodków raportujących. Modele ważone wykazały kierunki trendów i poziomy istotności zgodne z analizami nieważonymi; w związku z tym, dla uproszczenia i przejrzystości interpretacji przedstawiono wyniki nieważone. Wszystkie analizy statystyczne wykonano przy użyciu programu R w wersji 4.3.2 (R Foundation for Statistical Computing, Wiedeń, Austria). Dwustronną wartość p <0,05 uznano za statystycznie istotną.