Diese retrospektive Studie wurde am 9. Mai 2025 vom Medizinischen Ethikausschuss des Baoquanling Hospital of Beidahuang Group geprüft und genehmigt (Genehmigungsnummer BQH-BDHG-EC-2025-056), und die informierte Zustimmung wurde aufgehoben. Die Studie wurde gemäß der Helsinkier Erklärung und der Datenschutzrichtlinie des Krankenhauses durchgeführt, und alle personenbezogenen Daten wurden vor der Analyse deklariert. Die in dieser Studie verwendeten Forschungsinstrumente sind in der Materialtabelle aufgeführt.
1. Studiendesign
Dies war eine retrospektive Querschnittsstudie mit nur einem Zentrum, die in der Abteilung für Endokrinologie und Stoffwechsel des Baoquanling Hospital in Beidahuang Group durchgeführt wurde und 126 berechtigte Patienten mit Typ-2-Diabetes umfasste, die zwischen dem 1. Januar 2022 und dem 31. Dezember 2024 das Krankenhaus besuchten. Das Datum des ambulanten Besuchs oder der Krankenhausaufnahme wurde als Indexdatum definiert, und die an diesem Tag erhaltenen Laborergebnisse wurden als Ausgangspunkt verwendet. Expositions- und Ergebnismessungen wurden zum gleichen Zeitpunkt oder innerhalb von 7 Tagen vor oder nach dem Indexdatum ermittelt. Alle Daten stammen aus historischen medizinischen Unterlagen ohne jegliche Intervention.
Der Bericht wurde gemäß der Erklärung Strengthening the Reporting of Observational Studies in Epidemiology (STROBE)11 für Querschnittsstudien erstellt, und die Forschungsfragen, Variablen, statistische Methoden und Strategien zur Sensitivitätsanalyse wurden vorab festgelegt, wobei die Ergebnisse gemäß dem vorgegebenen Plan präsentiert wurden.
2. Untersuchung der Population
Die Studienpopulation wurde nacheinander aus dem elektronischen Patientenaktensystem gemäß den vorgegebenen Ein- und Ausschlusskriterien gescreent. Die Einschlusskriterien waren 18–80 Jahre, eine dokumentierte Diagnose von Typ-2-Diabetes in der Krankenakte, verfügbare Ergebnisse für Nüchternplasmaglukose, Nüchterninsulin, HS-CRP, Urinalbumin, Urinkreatinin und Serumkreatinin innerhalb von 7 Tagen vor oder nach dem Indexdatum, eGFR ≥ 60 mL/min/1,73 m2 (CKD-EPI 2021), eine frisch entnommene Urinprobe nach dem Standardtestverfahren, und vollständige klinische Daten sowie wichtige Kovariaten oder einen fehlenden Anteil, der die vorgegebenen Behandlungskriterien erfüllte.
Die Ausschlusskriterien umfassten Hinweise auf akute Entzündungen oder Infektionen, hs-CRP > 10 mg/L, Harnwegsinfektionen, Hämaturie oder Pyurie, Schwangerschaft oder Laktation, eine eindeutige Vorgeschichte nicht-diabetischer Nierenerkrankung oder bildgebende Hinweise auf strukturelle Nierenerkrankungen, UACR ≥ 300 mg/g, systemische Glukokortikoidtherapie innerhalb der vorangegangenen 3 Monate, akute kardiovaskuläre oder zerebrovaskuläre Ereignisse oder größere Operationen innerhalb der letzten 3 Monate. und eine Malignität, die Chemotherapie oder Immuntherapie erhält.
Nachdem Unterlagen, die den Kriterien nicht entsprachen, ausgeschlossen wurden, wurde die endgültige Probe entnommen. Unter einer zweiseitigen α = 0,05 und 80 % potenziell entsprach die endgültige Stichprobengröße von 126 in dieser Studie einem nachweisbaren absoluten Korrelationskoeffizienten von etwa 0,25. Diese Aussage war lediglich eine nachträgliche Beschreibung des nachweisbaren Effektbereichs und stellte keine a priori Stichprobengrößenschätzung dar. Analysen mit HOMA-IR-Berechnungen beschränkten sich auf Teilnehmer, die keine exogene Insulintherapie erhielten, und die Stichprobengröße dieser Untergruppe wurde als beobachtet angegeben.
3. Datenquellen und Erhebungsverfahren
Laut einem einheitlichen Datenwörterbuch extrahierten die Forscher demografische Informationen (Alter, Geschlecht), Diabetesdauer, Größe, Gewicht, Rauchstatus, Alkoholkonsum, systolischen Blutdruck (SBP) und diastolischen Blutdruck (DBP) (die zweite Messung nach zwei Messungen am Indexdatum), Komorbiditätsaufzeichnungen und Medikamenteninformationen aus den elektronischen Krankenakten.
Medikamentenvariablen umfassten Renin-Angiotensin-Systeminhibitoren (RASi), Natrium-Glukose-Cotransporter-2-Inhibitoren (SGLT2i), glukagonähnliche Peptid-1-Rezeptoragonisten (GLP-1RA) und Statine; Kontinuierliche Nutzung für ≥ 3 Monate bevor das Indexdatum als "Nutzung" vermerkt wurde. Die Dauer des Diabetes wurde als die Anzahl der Jahre von der Diagnose bis zum Indexdatum definiert. Der Body-Mass-Index (BMI) wurde als Gewicht (kg)/Körpergröße2 (m2) berechnet.
Labortests wurden auf einer einheitlichen Plattform durchgeführt und unterlagen interner Qualitätskontrolle sowie externer Qualitätsbewertung. hs-CRP wurde mittels hochsensitiver Immunoturbidimetrie gemessen, mit einer Nachweisgrenze von ≤ 0,1 mg/L; FINS wurde mittels chemilumineszenz mit In-Batch-Kalibrierung mit Kalibratoren gemessen; FPG wurde mit der Hexokinase-Methode gemessen; Glykiertes Hämoglobin (HbA1c) wurde mittels Hochleistungs-Flüssigkeitschromatographie gemessen; Urinalbumin wurde mittels Immunturbidimetrie gemessen, Urinkreatinin mit der enzymatischen Methode, und das Verhältnis wurde als mg/g angegeben; Serumkreatinin wurde mit der enzymatischen Methode mit IDMS-nachverfolgbarer Kalibrierung gemessen.
Alle Proben wurden innerhalb von 2 Stunden nach der Entnahme oder nach kurzzeitiger Lagerung bei 4 °C getestet. Wenn mehrere Ergebnisse für denselben Besuch verfügbar waren, wurden die am selben Tag wie das Indexdatum erhaltenen Ergebnisse priorisiert.
4. Variablendefinitionen und Messung
Zu den Expositionsvariablen gehörten hs-CRP und HOMA-IR. hs-CRP wurde in mg/L angegeben, als kontinuierliche Variable in die primären Analysen eingetragen und für Trendanalysen in Stichprobentertilen kategorisiert. FINS wurde in μU/mL angegeben, FPG in mmol/L, und die Berechnungsformel für HOMA-IR lautete wie folgt:6:
HOMA - IR = (FINS FPG)/22,5
Analysen mit HOMA-IR beschränkten sich auf Teilnehmer, die kein exogenes Insulin verwendeten und am selben Tag FPG- und FINS-Ergebnisse zeigten. In dieser Untergruppe wurden HOMA-IR und hs-CRP gemeinsam in das mehrdimensionale Modell aufgenommen. UACR wurde mit Urinalbumin und Urinkreatinin berechnet, gemessen in derselben Probe, und wurde als mg/g12 angegeben. Urinkreatinin wurde bei Bedarf auf Gramm harmonisiert, um das Verhältnis zu berechnen. In diesem Datensatz war kein UACR-Wert null. Bei Urinalbuminwerten unterhalb der unteren Nachweisgrenze (<2,0 mg/L) gab das Labor das Ergebnis mit <2,0 mg/L an, und diese unterdetektiveren Werte wurden vor der Datenanalyse durch die Hälfte des unteren Nachweislimits (1,0 mg/L) ersetzt. Um den Einfluss der Rechtsverrägheit zu verringern, wurde eine natürliche Logaritmumtransformation angewandt, und als primäres Ergebnis wurde LN(UACR) verwendet.
Wenn innerhalb von 7 Tagen vor oder nach dem Indexdatum mehrere Urintests für denselben Teilnehmer verfügbar waren, wurde nur die Probe, die dem Indexdatum am nächsten kam, bezeichnet als Morgenurin oder Ersturinurin, aufbewahrt. Wenn nur eine zufällige Urinprobe vorlag, wurde der Probentyp erfasst und die Sensitivitätsanalysen auf die Teilmenge des morgendlichen Urins beschränkt. Mikroalbuminurie wurde als UACR ≥ 30 mg/g definiert und als binäres Surrogatergebnis in der logistischen Regression verwendet.
eGFR wurde mit der CKD-EPI 2021 Kreatiningleichung13 berechnet. Als Serumkreatinin in μmol/L angegeben wurde, wurde es vor der Berechnung in mg/dL (μmol/L ÷ 88,4) umgerechnet. eGFR wurde als mL/min/1,73 m² angegeben und als kontinuierliches sekundäres Endergebnis in Sensitivitätsanalysen verwendet. Die Konsistenz der Ergebnisse wurde in der eGFR 60–89 mL/min/1,73 m² Subgruppe untersucht. Das primäre Modell wurde vordefiniert, um Alter, Geschlecht, Diabetesdauer, SBP und HbA1c zu berücksichtigen, und die Sensitivitätsmodelle umfassten außerdem BMI, Rauchstatus, Alkoholkonsum, Nutzung von RASi, Verwendung von SGLT2i, Verwendung von GLP-1RA und Statinkonsum.
5. Datenmanagement und Vorverarbeitung
Nach der Deidentifikation wurden die Daten als analytischer Datensatz exportiert, und ein kurzes Datenwörterbuch mit Variablennamen, Definitionen, Einheiten und Codierungsregeln für die analytischen Variablen wurde in ergänzender Tabelle 1 bereitgestellt. Duplikate Datensätze wurden nach Indexdatum zusammengeführt, und Schlüsselvariablen wurden auf logische Konsistenz überprüft.
Als sowohl EMR- als auch LIS-Zeitstempel verfügbar waren, wurde die Zeit der LIS-Probenentnahme als primärer Zeitstempel für die zeitliche Ausrichtung verwendet; wenn nicht verfügbar, wurde die LIS-Berichts-/Verifikationszeit verwendet, während das EMR-Besuchsdatum nur zur Festlegung des Indexdatums verwendet wurde. Wenn innerhalb des vorgegebenen Zeitfensters mehrere zulässige Ergebnisse verfügbar waren, wurde das Ergebnis, das dem Indexdatum am nächsten lag, beibehalten; Wenn zwei Ergebnisse gleich nah beieinander lagen, wurde das Ergebnis am selben Tag priorisiert, wobei die Urinproben weiter nach der vorgegebenen Probentypregel ausgewählt wurden.
Fehlende Werte wurden nach der vorgegebenen hierarchischen Strategie behandelt: Wenn der fehlende Anteil eines einzelnen Kovariaten ≤ 10 % betrug, verwendeten die primären Analysen einen vollständigen Fallansatz; Wenn sie 10 % überschritt, wurde eine mehrfache Imputation mit Kettengleichungen mit 20 Imputationen durchgeführt, wobei das Imputationsmodell die Expositionen, das Ergebnis und alle Kovariaten einschloss, und die imputierten Ergebnisse mit den Vollfallergebnissen verglichen wurden.
Teilnehmer mit hs-CRP > 10 mg/L oder Hinweisen auf akute Entzündung wurden von den Primäranalysen ausgeschlossen. Kontinuierliche Variablen wurden mit dem Shapiro-Wilk-Test und den Q-Q-Diagrammen für die Verteilung bewertet, und UACR war natürlich logaritmiert. Wenn hs-CRP oder HOMA-IR ausgeprägte Verzerrungen zeigten, wurden in Sensitivitätsanalysen Rang- oder Logaritmumtransformationen durchgeführt. Kategoriale Variablen wurden gemäß vorgegebenen Regeln als binäre oder ordinale Variablen codiert. Einflussreiche Beobachtungen wurden durch absolute, studentisierte Residuen > 3 oder Cooks Abstand > 4/n identifiziert, und die Primärmodelle wurden nach Ausschluss dieser Beobachtungen in Sensitivitätsanalysen wiederholt.
Vordefinierte Untergruppenanalysen umfassten Geschlechtsuntergruppen und Untergruppen basierend auf HbA1c < 7 % und ≥ 7 %. Vordefinierte Sensitivitätsanalysen umfassten Einschränkung auf morgendliche Urinproben, Hinzufügen von Medikamentenvariablen zu den Hauptmodellen, Ausschluss einflussreicher Beobachtungen, Verwendung robuster Standardfehler anstelle konventioneller Standardfehler und Einschränkung auf Teilnehmer mit eGFR 60–89 mL/min/1,73 m2.
6. Statistische Analyse
Alle statistischen Analysen wurden in R-Version 4.3.2 durchgeführt. Kontinuierliche Variablen wurden mit dem Shapiro-Wilk-Test und Q-Q-Diagrammen für die Verteilung bewertet. Normal verteilte Daten wurden als Mittelwert ± Standardabweichung angegeben, nicht-normalverteilte Daten als Median (Interquartilbereich) und kategorische Variablen als Frequenz und Prozentsatz.
Die Ausgangswerte wurden je nach Erreichen des Mikroalbuminurie-Schwellenwerts beschrieben, und Gruppenvergleiche wurden mit dem unabhängigen Proben-t-Test, dem Mann-Whitney-U-Test, dem Chi-Quadrat-Test oder dem Fisher-Exakt-Test je nach Variablentyp und -verteilung durchgeführt. Die Spearman-Rangkorrelationsanalyse wurde zwischen ln-UACR und hs-CRP sowie zwischen ln-UACR und HOMA-IR durchgeführt, wobei Korrelationskoeffizienten und ihre 95%-Konfidenzintervalle berechnet wurden, wobei Konfidenzintervalle durch Fishers z-Transformation geschätzt wurden. Korrelationsdiagramme wurden als Streudiagramme präsentiert, überlagert mit lokal gewichteten regressionsangepassten Linien, um den Trend darzustellen.
Mehrvariable lineare Regressionsmodelle wurden mit ln-UACR als abhängiger Variable konstruiert. In der Subgruppe, die kein exogenes Insulin verwendete, wurden hs-CRP und HOMA-IR gemeinsam in das Modell aufgenommen, mit Anpassung für Alter, Geschlecht, Diabetesdauer, SBP und HbA1c; In der Gesamtstichprobe wurde ein separates Modell eingebaut, das nur HS-CRP umfasste, um die Gesamtkorrelation zu untersuchen. Standardisierte Regressionskoeffizienten, 95%-Konfidenzintervalle und p-Werte wurden berichtet, und die Änderung des Bestimmungskoeffizienten vor und nach der Einbeziehung der Expositionsvariablen wurde berichtet. Die Kollinearität wurde anhand von Varianzinflationsfaktoren mit einem Schwellenwert von 5 bewertet.
Nachdem hs-CRP und HOMA-IR in Stichproben-Tertilen kategorisiert wurden, definiert durch die empirischen 33,3. und 66,7. Perzentile der entsprechenden analytischen Stichprobe (hs-CRP: die Gesamtstichprobe für Gesamtanalysen und die nicht-insulinbasierte Untergruppe für Untergruppen-/Gemeinschaftsanalysen; HOMA-IR: nur die Nicht-Insulin-Untergruppe), wobei T1 ≤ untere Grenzwert, T2 > untere Grenzwert ≤ oberen Grenzwert und T3 > oberer Grenzwert eingegeben wurde, wurde der Medianwert jedes Tertils als kontinuierliche Variable eingegeben, um den linearen Trend zu testen, und die angepassten marginalen Mittelwerte von ln-UACR über die Tertile sowie der p-Wert für den Trend wurden berichtet.
In der explorativen Analyse wurden hohe HS-CRP und hohe HOMA-IR durch die oberen tertilen Cutoffs definiert, die aus der Nicht-Insulin-Subgruppe für die 2×2-Gelenk-Analyse abgeleitet wurden, und eine 2×2-Gruppierung wurde erstellt, um die adjustierten marginalen Mittelwerte von Ln-UACR zwischen den Gruppen zu vergleichen. Ein Interaktionsterm wurde hinzugefügt, um die statistische Wechselwirkung zu testen, und der angepasste Mittelwert zwischen den Gruppen "high × high" und "low × low" wurde berichtet; Diese Analyse wurde nicht kausal interpretiert.
Für die Surrogat-Outcome-Analyse wurde eine multivariable logistische Regression mit dem Erreichen der Mikroalbuminurie-Schwelle als abhängige Variable konstruiert, und das Odds Ratio sowie das 95%-Konfidenzintervall entsprechend jeder 1-Standardabweichungs-Erhöhung in HS-CRP oder HOMA-IR wurden gemeldet.
Die multivariable lineare Regression wurde mit eGFR als abhängige Variable konstruiert, um die Zusammenhängungsrichtung zwischen den Expositionen und der glomerulären Filtrationsrate zu untersuchen.
Die Sensibilitätsanalysen umfassten Einschränkung auf morgendliche Urinproben, zusätzliche Anpassung für BMI, Rauchstatus, Alkoholkonsum und die vier Medikamentenkategorien in den Primärmodellen, Ausschluss einflussreicher Beobachtungen, Verwendung von Huber-White-robusten Standardfehlern sowie Wiederholung der Primärmodelle in der eingeschränkten Population mit UACR < 300 mg/g und eGFR ≥ 60 mL/min/1,73 m2. Die Modelldiagnostik verwendete Residual-versus-Fitting-Diagramme, Q-Q-Diagramme und den Shapiro-Wilk-Test zur Bewertung von Residualnormalität und Linearität sowie den Breusch-Pagan-Test zur Beurteilung der Homoskedastik. War Heteroscedastizität vorhanden, wurden robuste Standardfehler gemeldet; Wenn eine klare Nichtlinearität festgestellt wurde, wurden in Sensitivitätsanalysen tertile Indikatorvariablen anstelle kontinuierlicher Variablen verwendet.
Korrelations- und Regressionsanalysen wurden mit dem Statistikpaket durchgeführt, die Kollinearität wurde mit dem Autopaket bewertet, robuste Standardfehler wurden mit den Sandwich- und LMTEST-Paketen erzielt, mehrere Imputationen wurden mit dem Mäusepaket durchgeführt und Zahlen wurden mit dem ggplot2-Paket generiert. Alle Tests waren beidseitig, die Signifikanzschwelle wurde auf α = 0,05 gesetzt, und es wurden 95%-Konfidenzintervalle angegeben.