Deze retrospectieve studie werd beoordeeld en goedgekeurd door de Medische Ethische Commissie van het Baoquanling Ziekenhuis van Beidahuang Group op 9 mei 2025 (Goedkeuringsnummer BQH-BDHG-EC-2025-056), en informed consent werd opgegeven. De studie werd uitgevoerd in overeenstemming met de Verklaring van Helsinki en het privacybeleid van het ziekenhuis, en alle persoonlijk identificeerbare informatie werd vóór de analyse gedeïdentificeerd. De onderzoeksinstrumenten die in deze studie worden gebruikt, staan vermeld in de Materiaallijst.
1. Studieontwerp
Dit was een retrospectieve dwarsdoorsnedestudie met één centrum, uitgevoerd op de afdeling Endocrinologie en Metabolisme van het Baoquanling Hospital of Beidahuang Group, inclusief 126 in aanmerking komende patiënten met type 2 diabetes die tussen 1 januari 2022 en 31 december 2024 het ziekenhuis bezochten. De datum van het poliklinische bezoek of ziekenhuisopname werd gedefinieerd als de indexdatum, en de laboratoriumresultaten die op die datum werden verkregen, werden als uitgangspunt gebruikt. Blootstellings- en uitkomstmetingen werden op hetzelfde tijdstip of binnen 7 dagen voor of na de indexdatum verkregen. Alle gegevens zijn afkomstig uit historische medische dossiers zonder enige interventie.
Het rapport werd opgesteld in overeenstemming met de verklaring Strengthening the Reporting of Observational Studies in Epidemiology (STROBE)11 voor dwarsdoorsnedestudies, en de onderzoeksvraag, variabelen, statistische methoden en gevoeligheidsanalysestrategieën werden vooraf gespecificeerd, met resultaten gepresenteerd volgens het vooraf gespecificeerde plan.
2. Studiepopulatie
De studiepopulatie werd achtereenvolgens gescreend vanuit het elektronische medische dossiersysteem volgens de vooraf gespecificeerde inclusie- en exclusiecriteria. De inclusiecriteria waren leeftijd 18–80 jaar, een gedocumenteerde diagnose van type 2 diabetes in het medisch dossier, beschikbare resultaten voor nuchtere plasmaglucose, nuchtere insuline, hs-CRP, urinealbumine, urinecreatinine en serumcreatinine binnen 7 dagen voor of na de indexdatum, eGFR ≥ 60 mL/min/1,73m2 (CKD-EPI 2021), een vers afgenomen urinemonster verwerkt volgens de standaard testprocedure, en het voltooien van klinische gegevens en belangrijke covariaten of een ontbrekend aandeel dat voldeed aan de vooraf gespecificeerde behandelingscriteria.
De uitsluitingscriteria omvatten bewijs van acute ontsteking of infectie, hs-CRP > 10 mg/L, urineweginfectie, hematurie of pyurie, zwangerschap of lactatie, een duidelijke voorgeschiedenis van niet-diabetische nierziekte of beeldvormend bewijs van structurele nierziekte, UACR ≥ 300 mg/g, systemische glucocorticoïde therapie binnen de voorgaande 3 maanden, acute cardiovasculaire of cerebrovasculaire gebeurtenissen of grote operaties binnen de voorgaande 3 maanden. en maligniteit die chemotherapie of immunotherapie krijgt.
Nadat gegevens die niet aan de criteria voldeden waren uitgesloten, werd het definitieve monster verkregen. Onder een tweezijdige α = 0,05 en 80% power kwam de uiteindelijke steekproefgrootte van 126 in deze studie overeen met een detecteerbare absolute correlatiecoëfficiënt van ongeveer 0,25. Deze uitspraak was slechts een post-hoc beschrijving van het detecteerbare effectbereik en vormde geen a priori steekproefgrootteschatting. Analyses met HOMA-IR-berekeningen waren beperkt tot deelnemers die geen exogene insulinetherapie ontvingen, en de steekproefgrootte van deze subgroep werd als waargenomen gerapporteerd.
3. Gegevensbronnen en verzamelingsprocedures
Volgens een geïntegreerd datadictionary haalden de onderzoekers demografische informatie (leeftijd, geslacht), duur van diabetes, lengte, gewicht, rookstatus, alcoholgebruik, systolische bloeddruk (SBP) en diastolische bloeddruk (DBP) (de tweede meting na twee metingen op de indexdatum), comorbiditeitsdossiers en medicatiegegevens uit de elektronische medische dossiers.
Medicatievariabelen omvatten renin-angiotensine systeem remmers (RASi), natriumglucose cotransporter 2-remmers (SGLT2i), glucagonachtige peptide-1 receptoragonisten (GLP-1RA) en statines; Continu gebruik gedurende ≥ 3 maanden vóór de indexdatum werd als "gebruik" geregistreerd. De duur van diabetes werd gedefinieerd als het aantal jaren vanaf diagnose tot de indexdatum. De body mass index (BMI) werd berekend als gewicht (kg)/lengte2 (m2).
Laboratoriumtesten werden uitgevoerd op een uniform platform en stonden onder interne kwaliteitscontrole en externe kwaliteitsbeoordeling. hs-CRP werd gemeten met hooggevoelige immunoturbidimetrie, met een detectielimiet van ≤ 0,1 mg/L; FINS werd gemeten door chemiluminescentie met kalibratie binnen de batch met kalibratoren; FPG werd gemeten met de hexokinasemethode; geglyceerde hemoglobine (HbA1c) werd gemeten met high-performance vloeistofchromatografie; Urinealbumine werd gemeten met immunoturbidimetrie, urinecreatinine met de enzymatische methode, en de verhouding werd uitgedrukt als mg/g; serumcreatinine werd gemeten met de enzymatische methode en IDMS-traceerbare kalibratie.
Alle monsters werden binnen 2 uur na verzameling of na kortdurende opslag bij 4 °C getest. Als er meerdere resultaten beschikbaar waren voor hetzelfde bezoek, werden de resultaten verkregen op dezelfde dag als de indexdatum geprioriteerd.
4. Variabelendefinities en meting
Blootstellingsvariabelen omvatten hs-CRP en HOMA-IR. hs-CRP werd uitgedrukt in mg/L, opgenomen in de primaire analyses als een continue variabele en gecategoriseerd in steekproeftertielen voor trendanalyses. FINS werd uitgedrukt in μU/mL, FPG in mmol/L, en de rekenformule voor HOMA-IR was als volgt:6:
HOMA - IR = (FINS FPG)/22,5
Analyses met HOMA-IR waren beperkt tot deelnemers die geen exogene insuline gebruikten en dezelfde dag FPG- en FINS-resultaten hadden. In deze subgroep werden HOMA-IR en hs-CRP samen in het multivariabele model opgenomen. UACR werd berekend met urinealbumine en urinecreatinine gemeten in hetzelfde monster en werd uitgedrukt als mg/g12. Urinecreatinine werd geharmoniseerd tot gram voor de verhouding wanneer nodig. In deze dataset was geen UACR-waarde nul. Voor urinealbuminewaarden onder de onderste detectiegrens (<2,0 mg/L) rapporteerde het laboratorium het resultaat als <2,0 mg/L, en deze waarden onder de detectie werden vóór de data-analyse vervangen door de helft van de onderste detectielimiet (1,0 mg/L). Om de invloed van rechtsscheefheid te verminderen, werd een natuurlijke logaritme-transformatie toegepast en werd ln(UACR) als primaire uitkomst gebruikt.
Als er binnen 7 dagen voor of na de indexdatum meerdere urinetests beschikbaar waren voor dezelfde deelnemer, werd alleen het monster het dichtst bij de indexdatum bewaard, aangeduid als ochtendurine of eerste urine van de plasing. Als er slechts een willekeurig urinemonster beschikbaar was, werd het type monster geregistreerd en werden gevoeligheidsanalyses beperkt tot de ochtendurine-subset. Microalbuminurie werd gedefinieerd als UACR ≥ 30 mg/g en werd gebruikt als een binaire surrogaatuitkomst bij logistische regressie.
eGFR werd berekend met behulp van de CKD-EPI 2021 creatininevergelijking13. Wanneer serumcreatinine werd gerapporteerd in μmol/L, werd dit omgerekend naar mg/dL (μmol/L ÷ 88,4) vóór de berekening. eGFR werd uitgedrukt als mL/min/1,73 m2 en gebruikt als een continue secundaire uitkomst in gevoeligheidsanalyses. De consistentie van de resultaten werd onderzocht in de eGFR 60–89 mL/min/1,73 m2 subgroep. Het primaire model was vooraf gespecificeerd om te corrigeren voor leeftijd, geslacht, diabetesduur, SBP en HbA1c, en de gevoeligheidsmodellen omvatten daarnaast BMI, rookstatus, alcoholgebruik, gebruik van RASi, gebruik van SGLT2i, gebruik van GLP-1RA en statinegebruik.
5. Gegevensbeheer en voorverwerking
Na de-identificatie werden gegevens geëxporteerd als de analytische dataset, en een kort datawoordenboek met variabelennamen, definities, eenheden en coderingsregels voor de analytische variabelen werd opgenomen in Aanvullende Tabel 1. Dubbele records werden samengevoegd op indexdatum en sleutelvariabelen werden gecontroleerd op logische consistentie.
Wanneer zowel EMR- als LIS-tijdstempels beschikbaar waren, werd de tijd voor het verzamelen van LIS-monsters gebruikt als primaire tijdstempel voor temporele uitlijning; indien niet beschikbaar, werd de LIS-rapport/verificatietijd gebruikt, terwijl de bezoekdatum van het EPD alleen werd gebruikt om de indexdatum te definiëren. Wanneer er meerdere in aanmerking komende resultaten beschikbaar waren binnen het vooraf bepaalde venster, bleef het resultaat dat het dichtst bij de indexdatum lag; Als twee resultaten even dicht bij elkaar lagen, kreeg het resultaat op dezelfde dag prioriteit, waarbij urinemonsters verder werden geselecteerd volgens de vooraf gespecificeerde monstertyperegel.
Ontbrekende waarden werden behandeld volgens de vooraf gespecificeerde hiërarchische strategie: als het ontbrekende aandeel van een enkele covariaat ≤ 10% was, gebruikten de primaire analyses een complete-case benadering; Als deze meer dan 10% bedroeg, werd meervoudige imputatie met kettingvergelijkingen met 20 imputaties uitgevoerd, waarbij het imputatiemodel de blootstellingen, uitkomst en alle covariaten omvatte, en werden de geïmputeerde resultaten vergeleken met de volledige casusresultaten.
Deelnemers met hs-CRP > 10 mg/L of bewijs van acute ontsteking werden uitgesloten van de primaire analyses. Continue variabelen werden beoordeeld op verdeling met behulp van de Shapiro-Wilk-test en Q-Q-grafieken, en UACR was natuurlijk log-getransformeerd. Als hs-CRP of HOMA-IR duidelijke scheefheid vertoonde, werden rang- of logaritme-transformaties uitgevoerd in gevoeligheidsanalyses. Categorische variabelen werden gecodeerd als binaire of ordinaalvariabelen volgens vooraf gespecificeerde regels. Invloedrijke waarnemingen werden geïdentificeerd door absolute gestudentiseerde residuen > 3 of Cook's afstand > 4/n, en de primaire modellen werden herhaald nadat deze waarnemingen waren uitgesloten in gevoeligheidsanalyses.
Vooraf gespecificeerde subgroepanalyses omvatten seksesubgroepen en subgroepen gebaseerd op HbA1c < 7% en ≥ 7%. Vooraf gespecificeerde sensitiviteitsanalyses omvatten beperking tot ochtendurinemonsters, toevoeging van medicatievariabelen aan de primaire modellen, uitsluiting van invloedrijke waarnemingen, gebruik van robuuste standaardfouten in plaats van conventionele standaardfouten, en beperking tot deelnemers met eGFR 60–89 mL/min/1,73m2.
6. Statistische analyse
Alle statistische analyses werden uitgevoerd in R versie 4.3.2. Continue variabelen werden beoordeeld op verdeling met behulp van de Shapiro-Wilk-test en Q-Q-plots. Normaal verdeelde gegevens werden uitgedrukt als gemiddelde ± standaarddeviatie, niet-normaal verdeelde gegevens als mediaan (interkwartielbereik), en categorische variabelen als frequentie en percentage.
Basislijnkenmerken werden beschreven op basis van of de microalbuminurie-drempel was bereikt, en vergelijkingen tussen groepen werden uitgevoerd met behulp van de onafhankelijke steekproeven t-test, Mann-Whitney U-test, chi-kwadraattest of Fisher's exacte test, afhankelijk van variabeletype en verdeling. Spearman-rangcorrelatieanalyse werd uitgevoerd tussen ln-UACR en hs-CRP en tussen ln-UACR en HOMA-IR, en correlatiecoëfficiënten en hun 95% betrouwbaarheidsintervallen werden berekend, waarbij betrouwbaarheidsintervallen werden geschat door Fishers z-transformatie. Correlatiegrafieken werden gepresenteerd als spreidingsdiagrammen overheen met lokaal gewogen regressielijnen om de trend te tonen.
Meervariabele lineaire regressiemodellen werden geconstrueerd met ln-UACR als afhankelijke variabele. In de subgroep die geen exogeen insuline gebruikte, werden hs-CRP en HOMA-IR samen in het model opgenomen, met aanpassing voor leeftijd, geslacht, diabetesduur, SBP en HbA1c; In de totale steekproef werd een apart model met alleen HS-CRP geïnstalleerd om de algehele correlatie te onderzoeken. Gestandaardiseerde regressiecoëfficiënten, 95% betrouwbaarheidsintervallen en p-waarden werden gerapporteerd, en de verandering in de bepalingscoëfficiënt vóór en na opname van de blootstellingsvariabelen werd gerapporteerd. Collineariteit werd beoordeeld met behulp van variantie-inflatiefactoren, met een drempel van 5.
Nadat hs-CRP en HOMA-IR werden gecategoriseerd in steekproeftertielen, gedefinieerd door de empirische 33,3e en 66,7e percentielen van de bijbehorende analytische steekproef (hs-CRP: de totale steekproef voor de totale analyses en de niet-insuline subgroep voor subgroep/gezamenlijke analyses; HOMA-IR: alleen de niet-insuline subgroep), met T1 ≤ de onderste grens, T2 > de onderste grens tot ≤ de bovenste grens, en T3 > de bovenste grens, werd de mediane waarde van elke tertiel ingevoerd als een continue variabele om de lineaire trend te testen, en werden de aangepaste marginale gemiddelden van ln-UACR over de tertielen en de p-waarde voor trend gerapporteerd.
In de exploratieve analyse werden hoge hs-CRP en hoge HOMA-IR gedefinieerd door de bovenste tertiel cutoffs, afgeleid van de niet-insuline subgroep voor de 2 × 2 gewrichtsanalyse, en werd een 2 × 2 groepering geconstrueerd om de aangepaste marginale gemiddelden van ln-UACR tussen groepen te vergelijken. Er werd een interactieterm toegevoegd om statistische interactie te testen, en het aangepaste gemiddelde verschil tussen de groepen "hoog × hoog" en "laag × laag" werd gerapporteerd; Deze analyse werd niet causaal geïnterpreteerd.
Voor surrogaatuitkomstanalyse werd multivariabele logistische regressie geconstrueerd met het bereiken van de microalbuminurie-drempel als afhankelijke variabele, en werden de odds ratio en het 95% betrouwbaarheidsinterval gerapporteerd voor elke 1 standaarddeviatie-toename in hs-CRP of HOMA-IR.
Multivariabele lineaire regressie werd geconstrueerd met eGFR als afhankelijke variabele om de richting van de associatie tussen de blootstellingen en de glomerulaire filtratiesnelheid te onderzoeken.
Gevoeligheidsanalyses omvatten beperking tot ochtendurinemonsters, aanvullende aanpassing voor BMI, rookstatus, alcoholgebruik en de vier medicatiecategorieën in de primaire modellen, uitsluiting van invloedrijke waarnemingen, gebruik van Huber-White robuuste standaardfouten en herhaling van de primaire modellen in de beperkte populatie met UACR < 300 mg/g en eGFR ≥ 60 mL/min/1,73m2. Modeldiagnostiek gebruikte residual-versus-fitted plots, Q-Q plots en de Shapiro-Wilk test om residuele normaliteit en lineariteit te beoordelen, en de Breusch-Pagan test om homoscedasticiteit te beoordelen. Als heteroscedasticiteit aanwezig was, werden robuuste standaardfouten gerapporteerd; Als duidelijke niet-lineariteit werd gevonden, werden tertiele indicatorvariabelen gebruikt in plaats van continue variabelen in gevoeligheidsanalyses.
Correlatie- en regressieanalyses werden uitgevoerd met het stats-pakket, collineariteit werd beoordeeld met het autopakket, robuuste standaardfouten werden verkregen met de sandwich- en lmtest-pakketten, meervoudige imputatie werd uitgevoerd met het muizenpakket en cijfers werden gegenereerd met het ggplot2-pakket. Alle tests waren tweezijdig, de significantiedrempel werd vastgesteld op α = 0,05, en er werden 95% betrouwbaarheidsintervallen gerapporteerd.