מחקר רטרואקטיבי זה נבדק ואושר על ידי ועדת האתיקה הרפואית של בית החולים באוקואנלינג של קבוצת ביידהואנג ב-9 במאי 2025 (אישור מס' BQH-BDHG-EC-2025-056), והסכמה מראשית בוטלה. המחקר נערך בהתאם להצהרת הלסינקי ולמדיניות ניהול פרטיות המידע של בית החולים, וכל המידע האישי המזהה הוסר לפני הניתוח. כלי המחקר ששימשו במחקר זה מופיעים בטבלת החומרים.
1. עיצוב מחקר
זהו מחקר חתך רטרוספקטיבי במרכז יחיד שנערך במחלקה לאנדוקרינולוגיה ומטבוליזם בבית החולים באוקואנלינג בקבוצת ביידהואנג, וכלל 126 מטופלים זכאים עם סוכרת סוג 2 שהגיעו לבית החולים בין 1 בינואר 2022 ל-31 בדצמבר 2024. תאריך הביקור האמבולטורי או האשפוז הוגדר כתאריך האינדקס, ותוצאות בדיקות המעבדה שהתקבלו באותו תאריך שימשו כקו בסיס. מדידות חשיפה ותוצאה הושגו באותו זמן או בתוך 7 ימים לפני או אחרי תאריך האינדקס. כל הנתונים נגזרו מרשומות רפואיות היסטוריות ללא כל התערבות.
הדוח הוכן בהתאם להצהרת חיזוק הדיווח על מחקרים תצפיתיים באפידמיולוגיה (STROBE)11 למחקרים חתכיים, ושאלת המחקר, משתנים, שיטות סטטיסטיות ואסטרטגיות ניתוח רגישות הוגדרו מראש, כאשר התוצאות הוצגו בהתאם לתוכנית שנקבעה מראש.
2. אוכלוסיית המחקר
אוכלוסיית המחקר נבדקה ברצף ממערכת הרשומות הרפואיות האלקטרונית בהתאם לקריטריוני הכללה והחרגה שנקבעו מראש. קריטריוני ההכללה היו: גיל 18–80, אבחנה מתועדת של סוכרת סוג 2 בתיק הרפואי, תוצאות זמינות לגלוקוז בפלזמה בצום, אינסולין בצום, hs-CRP, אלבומין שתן, קריאטינין שתן וקריאטינין בסרום בתוך 7 ימים לפני או אחרי תאריך האינדקס, eGFR ≥ 60 מ"ל/דקה/1.73 מ' 2 (CKD-EPI 2021), דגימת שתן שנאספה טרייה ועובדה לפי הליך הבדיקה הסטנדרטי, ונתונים קליניים מלאים ומשתנים מרכזיים או חלק חסר שעומד בקריטריוני הטיפול שנקבעו מראש.
קריטריוני ההדרה כללו עדויות לדלקת חריפה או זיהום, hs-CRP > 10 מ"ג/ליטר, זיהום בדרכי השתן, המטוריה או פיוריה, הריון או הנקה, היסטוריה ברורה של מחלת כליות שאינה סוכרתית או ראיות דימות למחלת כליות מבנית, UACR ≥ 300 מ"ג/גרם, טיפול גלוקוקורטיקואידי מערכתית בשלושת החודשים הקודמים, אירועים אירועי לב וכלי דם חריפים או כלי דם מוחיים או ניתוח משמעותי בשלושת החודשים הקודמים. ומממאירות בקבלת כימותרפיה או אימונותרפיה.
לאחר שהוצאו רשומות שלא עמדו בקריטריונים, התקבלה המדגם הסופי. תחת α דו-צדדי = 0.05 ועוצמה של 80%, גודל המדגם הסופי של 126 במחקר זה התאים למקדם מתאם מוחלט שניתן לזיהוי של כ-0.25. הצהרה זו הייתה רק תיאור לאחר מכן של טווח ההשפעה הניתן לזיהוי ולא מהווה הערכת גודל מדגם מראש בלבד. ניתוחים שכללו חישובי HOMA-IR הוגבלו למשתתפים שלא קיבלו טיפול אינסולין חיצוני, וגודל המדגם של תת-קבוצה זו דווח כפי שנצפה.
3. מקורות נתונים ונהלי איסוף
לפי מילון נתונים מאוחד, החוקרים חילצו מידע דמוגרפי (גיל, מין), משך הסוכרת, גובה, משקל, מצב עישון, שימוש באלכוהול, לחץ דם סיסטולי (SBP) ולחץ דם דיאסטולי (DBP) (הקריאה השנייה לאחר שתי מדידות בתאריך האינדקס), רשומות מחלות נלוות ומידע על תרופות מהרשומות הרפואיות האלקטרוניות.
משתני התרופות כללו מעכבי מערכת רנין-אנגיוטנסין (RASi), מעכבי נתרן-גלוקוז קוטרנספורטר 2 (SGLT2i), אגוניסטים לקולטני פפטיד-1 דמויי גלוקגון (GLP-1RA), וסטטינים; שימוש רציף במשך ≥ 3 חודשים לפני תאריך האינדקס נרשם כ"שימוש". משך הסוכרת הוגדר כמספר השנים מהאבחון ועד לתאריך האינדקס. מדד מסת הגוף (BMI) חושב כמשקל (ק"ג)/גובה 2 (מטר2).
בדיקות מעבדה בוצעו על פלטפורמה מאוחדת והיו כפופות לבקרת איכות פנימית והערכת איכות חיצונית. hs-CRP נמדד באמצעות אימונוטורבידימטריה ברגישות גבוהה, עם מגבלת גילוי של ≤ 0.1 מ"ג/ליטר; FINS נמדדה באמצעות כימלומינסצנציה עם כיול בתוך אצווה באמצעות כיולים; FPG נמדד בשיטת הקסוקינאז; המוגלובין גליגליקטיבי (HbA1c) נמדד באמצעות כרומטוגרפיה נוזלית בעלת ביצועים גבוהים; אלבומין השתן נמדד באימונוטורבידימטריה, קריאטינין שתן בשיטה אנזימטית, והיחס הובע כ-mg/g; קריאטינין בסרום נמדד בשיטה אנזימטית עם כיול שניתן למעקב ב-IDMS.
כל הדגימות נבדקו תוך שעתיים לאחר האיסוף או לאחר אחסון קצר טווח בטמפרטורה של 4 מעלות צלזיוס. אם היו זמינות מספר תוצאות לאותו ביקור, התוצאות שהושגו באותו יום עם תאריך האינדקס קיבלו עדיפות.
4. הגדרות ומדידה של משתנים
משתני החשיפה כללו hs-CRP ו-HOMA-IR. hs-CRP הוצג ב-mg/L, הוכנס לניתוחים הראשיים כמשתנה רציף, וסווג לשלישיית דגימות לניתוח מגמות. FINS הובע במיקרוU/מיליטר, FPG ב-mmol/L, ונוסחת החישוב ל-HOMA-IR הייתה כדלקמן:6:
HOMA - IR = (FINS FPG)/22.5
ניתוחים שכללו HOMA-IR הוגבלו למשתתפים שלא השתמשו באינסולין חיצוני שקיבלו תוצאות FPG ו-FINS באותו יום. בתת-קבוצה זו, HOMA-IR ו-hs-CRP הוכנסו יחד למודל הרב-משתנים. UACR חושב באמצעות אלבומין שתן וקריאטינין שתן שנמדדו באותו מדגם, והובטא כ-mg/g 12. קריאטינין בשתן הותאם לגרמים לחישוב יחס בעת הצורך. במאגר נתונים זה, לא היה ערך UACR אפסי. עבור ערכי אלבומין בשתן מתחת לגבול הגילוי הנמוך (<2.0 מ"ג/ליטר), המעבדה דיווחה על התוצאה כ-<2.0 מ"ג/ליטר), וערכי הגילוי הנמוך הוחלפו במחצית מגבלת הגילוי הנמוכה (1.0 מ"ג/ליטר) לפני ניתוח הנתונים. כדי להפחית את השפעת ההטיה הימנית, הוחל טרנספורמציה טבעית של לוגריתם, והשתמשו ב-ln(UACR) כתוצאה העיקרית.
אם היו זמינות מספר בדיקות שתן לאותו משתתף בתוך 7 ימים לפני או אחרי תאריך האינדקס, נשמר רק המדגם הקרוב ביותר לתאריך האינדקס, המסומן כשתן בוקר או שתן ריק ראשון. אם הייתה רק דגימת שתן אקראית זמינה, סוג המדגם נרשם, וניתוחי רגישות הוגבלו לתת-הקבוצה של שתן הבוקר. מיקרואלבומינוריה הוגדרה כ-UACR ≥ 30 מ"ג/גרם ושימשה כתוצאה בינארית חלופית ברגרסיה לוגיסטית.
eGFR חושב באמצעות משוואת קריאטינין CKD-EPI2021 13. כאשר הקריאטינין בסרום דווח במיקרומול/ליטר, הוא הומר ל-mg/dL (μmol/L ÷ 88.4) לפני החישוב. eGFR הוצג כ-mL/min/1.73 m2 ושימש כתוצאה משנית רציפה בניתוחי רגישות. עקביות התוצאות נבדקה בתת-הקבוצה eGFR 60–89 mL/min/1.73m 2 . המודל הראשי הוגדר מראש להתאמה לגיל, מין, משך הסוכרת, SBP ו-HbA1c, ומודלים הרגישות כללו גם BMI, מצב עישון, שימוש באלכוהול, שימוש ב-RASi, שימוש ב-SGLT2i, שימוש ב-GLP-1RA ושימוש בסטטין.
5. ניהול נתונים ועיבוד מוקדם
לאחר הסרת הזיהוי, הנתונים יוצאו כמערך נתונים אנליטי, ומילון נתונים קצר המסכם שמות משתנים, הגדרות, יחידות וכללי קידוד למשתנים האנליטיים סופק בטבלה משלימה 1. רשומות כפולות אוחדו לפי תאריך אינדקס, ומשתנים מרכזיים נבדקו לעקביות לוגית.
כאשר היו זמינות חותמות זמן של EMR ו-LIS, זמן איסוף דגימות LIS שימש כחותמת הזמן העיקרית ליישור הזמן; אם לא היה זמין, נעשה שימוש בדוח/זמן האימות של LIS, בעוד שתאריך הביקור ב-EMR שימש רק להגדרת תאריך האינדקס. כאשר מספר תוצאות זכאיות היו זמינות בתוך חלון הזמן שנקבע מראש, התוצאה הקרובה ביותר לתאריך האינדקס נשמרה; אם שתי תוצאות היו קרובות באותה מידה, התוצאה באותו יום קיבלה עדיפות, כאשר דגימות השתן נבחרו בהתאם לכלל סוג הדגימה שנקבע מראש.
ערכים חסרים טופלו לפי האסטרטגיה ההיררכית שנקבעה מראש: אם החלק החסר של משתנה אחד היה ≤ 10%, הניתוחים הראשיים השתמשו בגישה של מקרה שלם; אם הוא עבר את ה-10%, בוצע הערכה מרובה באמצעות משוואות שרשרתיות עם 20 הערכות, כאשר מודל ההערכה כלל את החשיפה, התוצאה וכל המשתנים, והתוצאות שהוחשפו הושוו לתוצאות המקרה המלא.
משתתפים עם hs-CRP > 10 מ"ג/ליטר או עדויות לדלקת חריפה הוצאו מהניתוחים הראשוניים. משתנים רציפים הוערכו להתפלגות באמצעות מבחן שפירו-וילק וגרפים Q-Q, ו-UACR עבר טרנספורמציה לוגרימית טבעית. אם hs-CRP או HOMA-IR הראו סטייה משמעותית, בוצעו טרנספורמציות דרגה או לוגריתיות בניתוחי רגישות. משתנים קטגוריים קודדו כמשתנים בינאריים או סודרים לפי כללים מוגדרים מראש. תצפיות משפיעות זוהו על ידי שאריות סטודנטים מוחלטות > 3 או מרחק קוק > 4/n, והמודלים הראשוניים חזרו על עצמם לאחר שהוצאו תצפיות אלו בניתוחי רגישות.
ניתוחי תת-קבוצות מוגדרים מראש כללו תת-קבוצות מין ותת-קבוצות המבוססות על HbA1c < 7% ו-≥ 7%. ניתוחי רגישות מוגדרים מראש כללו הגבלה לדגימות שתן בבוקר, הוספת משתני תרופות למודלים הראשיים, הוצאת תצפיות משפיעות, שימוש בשגיאות תקן חזקות במקום שגיאות סטנדרטיות קונבנציונליות, והגבלה למשתתפים עם eGFR 60–89 m/min/1.73m 2.
6. ניתוח סטטיסטי
כל הניתוחים הסטטיסטיים בוצעו בגרסה R 4.3.2. משתנים רציפים הוערכו להתפלגות באמצעות מבחן שפירו-וילק וגרפים Q-Q. הנתונים המפוזרים בדרך כלל הוצגו כסטיית תקן ממוצעת ±, נתונים לא מבוזרים כחציון (טווח בין-רבעוני), ומשתנים קטגוריים כתדירות ואחוזים.
מאפייני הבסיס תוארו בהתאם לשאלה האם הושג סף המיקרואלבומינוריה, והשוואות בין קבוצות בוצעו באמצעות מבחן t דגימות עצמאיות, מבחן Mann-Whitney U, מבחן כי-ריבוע או מבחן מדויק של פישר, בהתאם לסוג המשתנה והתפלגות. ניתוח קורלציה בדירוג ספירמן בוצע בין ln-UACR ו-hs-CRP ובין ln-UACR ל-HOMA-IR, וחושבו מקדמי קורלציה ורווחי האמון של 95%, כאשר טווחי האמון הוערכו על ידי טרנספורמציית z של פישר. גרפי קורלציה הוצגו כקטעי פיזור עם קווים מותאמים עם רגרסיה משוקללת מקומית כדי להראות את המגמה.
מודלים של רגרסיה ליניארית מרובת משתנים נבנו כאשר ln-UACR הוא המשתנה התלוי במונחים. בתת-הקבוצה שלא השתמשה באינסולין אקסוגני, הוכנסו יחד למודל hs-CRP ו-HOMA-IR, עם התאמת גיל, מין, משך הסוכרת, SBP ו-HbA1c; במדגם הכולל, הותקן מודל נפרד שכלל רק HS-CRP כדי לבחון את הקורלציה הכוללת. דווחו מקדמי רגרסיה סטנדרטיים, רווחי ביטחון של 95% וערכי p, ודווח על השינוי במקדם הקביעה לפני ואחרי הכללת משתני החשיפה. הקולינאריות הוערכה באמצעות גורמי אינפלציה של שונות, עם סף של 5.
לאחר ש-hs-CRP ו-HOMA-IR סווגו לשלישיות מדגם, שהוגדרו על ידי האחוזון האמפירי 33.3 ו-66.7 של המדגם האנליטי המתאים (hs-CRP: המדגם הכולל לניתוחים כוללים ותת-הקבוצה הלא-אינסולינית לניתוחי תת-קבוצות/משותפים; HOMA-IR: תת-הקבוצה שאינה אינסולין בלבד), כאשר T1 ≤ החיתוך התחתון, T2 > החיתוך התחתון ≤ החיתוך העליון, ו-T3 > החיתוך העליון, הערך החציוני של כל שלישייה הוזן כמשתנה רציף לבדיקת מגמה ליניארית, ודווחו הממוצעים השוליים המתוקנים של ln-UACR בין שלישיות וערך p למגמה.
בניתוח החקר, hs-CRP גבוה ו-HOMA-IR גבוה הוגדרו על ידי חיתוך השלישייה העליונה, שמקורם בתת-הקבוצה הלא-אינסולינית לניתוח המפרק 2 × 2, ונבנה קבוצה של 2 × 2 להשוואת ממוצעים שוליים מותאמים של ln-UACR בין קבוצות. נוסף מונח אינטראקציה לבדיקת אינטראקציה סטטיסטית, ודווח על ההפרש הממוצע המותאם בין קבוצות "גבוה× גבוה" ו"נמוך × נמוך"; ניתוח זה לא פורש באופן סיבתי.
לניתוח תוצאות חלופיות, נבנתה רגרסיה לוגיסטית רב-משתנית עם הגעה לסף המיקרואלבומינוריה כמשתנה התלוי, ודווחו יחס הסיכויים ורווח האמון של 95% התואמים לכל עלייה של סטיית תקן אחת ב-hs-CRP או HOMA-IR.
רגרסיה ליניארית רב-משתנית נבנתה כאשר eGFR הוא המשתנה התלוי כדי לבחון את כיוון הקשר בין החשיפות לקצב הסינון הגלומרולרי.
ניתוחי הרגישות כללו הגבלה לדגימות שתן בבוקר, התאמות נוספות ל-BMI, מצב עישון, שימוש באלכוהול וארבע קטגוריות התרופות במודלים הראשוניים, הוצאת תצפיות משפיעות, שימוש בשגיאות תקן חזקות של Huber-White, וחזרה על מודלים ראשוניים באוכלוסייה מוגבלת עם UACR < 300 מ"ג/גרם ו-eGFR ≥ 60 מ"ל/דקה/1.73מ"ר. אבחון המודל השתמש בגרפים של שארית לעומת התאמה מותאמת, גרפים Q-Q, ומבחן שפירו-וילק להערכת נורמליות וליניאריות שארית, ובמבחן ברויש-פאגן להערכת הומוסקדסטיות. אם הייתה הטרוסדסטיות, דווחו שגיאות תקן חזקות; אם נמצאה אי-ליניאריות ברורה, השתמשו במשתני אינדיקטור שלישיים במקום משתנים רציפים בניתוחי רגישות.
בוצעו ניתוחי קורלציה ורגרסיה באמצעות חבילת הסטטיסטיקות, הוערכו קולינאריות באמצעות חבילת הרכב, שגיאות תקן חזקות התקבלו באמצעות חבילות הסנדוויץ' ו-lmtest, בוצעו אימפוטציות מרובות באמצעות חבילת העכברים, ונוצרו נתונים באמצעות חבילת ggplot2. כל המבחנים היו דו-צדדיים, סף המובהקות נקבע ב-α = 0.05, ודווחו טווחי ביטחון של 95%.