$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
מקורות נתונים של FAERS
הנתונים האמיתיים למחקר זה התקבלו ממסד הנתונים של FAERS (https://fis.fda.gov/extensions/FPD-QDE-FAERS/FPD-QDE-FAERS.html). זהו מאגר נתונים נגיש לציבור ואנונימי, ולכן לא נדרש אישור אתי למחקר זה. מידע על חמשת התרופות הכלולות מוצג בטבלה 1. החיפוש בוצע תחילה על ידי הגבלת תופעות הלוואי לשבר בירך, ומשך החיפוש היה מרבעון הראשון של 2014 ועד הרבעון הרביעי של 2023. כדי להבטיח נתונים אמינים ויציבים, המחקר תקן את המונחים של תופעות לוואי מדווחות באמצעות מילון MedDRA גרסה 26.142. קבצי ASCII רבעוניים של FAERS מה-Q1 2014 עד Q4 2023 הורדו ויובאו לניתוח. טבלאות FAERS שהופקו כללו DEMO, DRUG, REAC, THER, RPSR ו-OUTC. טבלאות אלו אוחדו בכל הרבעים לפני ההקרנה. הדוחות קושרו באמצעות CASEID ו-PRIMARYID כדי להבטיח עקביות בין מידע דמוגרפי, תרופות, תגובה, טיפול, מדווח ותוצאות. התרופות המטרה כללו אתנרספט, אדלימומאב, אינפליקסימאב, צ'רטוליזומאב פגול וגולימומאב. שמות התרופות בטבלת DRUG היו סטנדרטיים על ידי המרת טקסט לאותיות גדולות, הסרת רווחים נוספים ובדיקת וריאציות איות במידת הצורך. תרופות מטרה זוהו באמצעות שמות גנריים סטנדרטיים בטבלת DRUGS, והגבלת תפקיד התרופה בוצעה באמצעות ROLE_COD = "PS", המצביע על התרופה החשודה העיקרית43. מפת הדרכים של טכניקת סינון המחקרים הספציפית מוצגת באיור 2. לאחר ניקוי הנתונים וסינון, התקבל מאגר נתונים ייחודי וניתוח של דוחות שבר ירך זכאים לניתוחים הבאים.
דוחות FAERS כפולים הוסרו לפני זיהוי אותות. כפילויות זוהו לפי CASEID ו-PRIMARYID. כאשר מספר דוחות שיתפו את אותו CASEID, הדוח האחרון נשמר לפי FDA_DT. אם מספר דיווחים כללו את אותו CASEID ו-FDA_DT, הדוח עם ה-PRIMARYID הגבוה ביותר נשמר. לאחר הסרת כפילות, כל CASEID תרם רק רשומה אחת למאגר הנתונים האנליטי הסופי. דוחות נכללו אם עמדו בכל הקריטריונים הבאים: תאריך דיווח בין רבעון 1 2014 לרבעון 4 2023; האירוע הלוואי סומן כ"שבר בירך"; לפחות אחד מחמשת מעכבי TNF נרשם בטבלת התרופות; ותפקיד הסמים סומן כחשוד העיקרי. הדיווחים הוצאו מכלל אם היו רשומות כפולות, חסרו מידע תקף על CASEID או PRIMARYID, לא היה רישום תואם של DRUG או REAC, לא כללו את אירוע הלוואי היעד, או שציינו את מעכב TNF היעד רק כתרופה נלווה או חשודה משנית.
מקורות נתוני GWAS לאקראיות מנדלית
נתוני החשיפה ל-TNF-α באקראיות המנדלית של מחקר זה התקבלו ממסד הנתונים של IEU OpenGWAS (https://gwas.mrcieu.ac.uk/), עם מזהה GWAS prot-c-3722_49_2 ממחקר של סוהר K ואחרים. אוכלוסיית המחקר הייתה ממוצא אירופי, ומספר ה-SNPs היה 501,428,44.
נתוני החשיפה ל-sTNFR1 באקראיות המנדלית של מחקר זה התקבלו ממאגר IEU OpenGWAS (https://gwas.mrcieu.ac.uk/), עם מזהה GWAS prot-c-2654_19_1 ממחקר של סוהר K ואחרים. אוכלוסיית המחקר הייתה ממוצא אירופי, ומספר ה-SNPs היה 501,428,44.
נתוני החשיפה ל-sTNFR2 באקראיות המנדלית למחקר זה התקבלו ממאגר הנתונים של IEU OpenGWAS (https://gwas.mrcieu.ac.uk/), עם מזהה GWAS prot-c-3152_57_1 ממחקר של Suhre K ואחרים. אוכלוסיית המחקר הייתה ממוצא אירופי, ומספר ה-SNPs היה 501,428,44.
נתוני התוצאה של שבר הירך, GWAS ID GCST90161240, שהופקדו בקטלוג GWAS (https://www.ebi.ac.uk/gwas/studies/GCST90161240), הם נתונים ממטא-אנליזה של GWAS רחב היקף שכלל 11,516 מקרים של שבר בירך ו-723,838 מקרים בביקורת, 45. סיווג המחלות תואם את הסיווג הבינלאומי למחלות (ICD; קודי ICD-10 S72.0–S72.2 וקוד ICD-9 820).
אישור אתי והסכמה מדעת התקבלו במחקרי GWAS המקוריים. מכיוון שהמחקר הנוכחי השתמש בנתוני FAERS אנונימיים זמינים לציבור ובסטטיסטיקות סיכום GWAS זמינות לציבור, לא נדרשה אישור אתי נוסף.
סביבת תוכנה ומימוש זרימת עבודה
כל הניתוחים בוצעו באמצעות גרסה 4.3.2 של R. ייבוא נתוני FAERS, ניקוי, מיזוג וטבלאות בוצעו באמצעות תהליכי עבודה מבוססי R. טבלאות נתונים יובאו באמצעות פונקציות כמו data.table::fread() או readr::read_delim(), אוחדו באמצעות CASEID ו-PRIMARYID, ועובדו באמצעות פונקציות dplyr. סטטיסטיקות תיאוריות ו-2 טבלאות × 2 אפשרויות נוצרו באמצעות סקריפטים מותאמים אישית של R.
ניתוחי אקראיות מנדליאניים נערכו באמצעות גרסה 0.5.6 של TwoSampleMR. מכשירי החשיפה הופקו באמצעות סף מובהקות של P < 1 × 10⁻5 או פורסמו מתוך סטטיסטיקות סיכום GWAS באמצעות מבני קלט תואמים ל-TwoSampleMR. איסוף מכשירים בוצע באמצעות clump_data() כאשר clump_r2 = 0.001 ו-clump_kb = 10,000. נתוני התוצאה הופקו או פורסמו באמצעות extract_outcome_data() או read_outcome_data(), בהתאם לפורמט המקור. מערכי הנתונים של חשיפה ותוצאה תואמו באמצעות harmonise_data(). הערכות סיבתיות נוצרו באמצעות mr() בשיטות האקראי המנדליאניות הבאות: MR-Egger, חציון משוקלל, משוקלל שונות הפוכה, מצב פשוט, ומצב משוקלל. ההטרוגניות הוערכה באמצעות mr_heterogeneity(), והפליאוטרופיה האופקית הוערכה באמצעות mr_pleiotropy_test(). כל מערכי הנתונים יובאו, נוקו, התאמרו ונותחו בסביבת תוכנה זו כדי להבטיח תהליך אנליטי עקבי וניתן לשחזור.
ניתוח פרמקוויג'ילנס
ניתוחים תיאוריים שימשו לסיכום תופעות לוואי הקשורות לשבר הירך הקשורים לחמש התרופות. לאחר מכן בוצעו ניתוחי זיהוי אותות באמצעות ארבעה אלגוריתמים של אי-פרופורציונליות, כולל יחס סיכויי דיווח (ROR), יחס דיווח פרופורציונלי (PRR), מתכווץ פואסון גמא רב-פריטי (MGPS) ורשת עצבית של הפצת ביטחון בייסיאנית (BCPNN). הקריטריונים לארבעת האלגוריתמים המרכזיים מוצגים בטבלה 246.
ניתוח רנדומיזציה מנדלית
סטטיסטיקות סיכום ל-TNF-α, sTNFR1 ו-sTNFR2 נאספו כמאגרי חשיפה, וסטטיסטיקות סיכום של שבר הירך הופקו כמאגר הנתונים לתוצאות. הניתוחים הוגבלו למאגרי נתוני מוצא אירופי כאשר היו זמינים, כדי להפחית את הטיית הסטרטיפיקציה של האוכלוסייה.
כדי למזער הטיה הנגרמת על ידי חוסר שוויון קישור ומכשירים חלשים, הוחלו הקריטריונים הבאים: סף מובהקות כלל-גנומי P < 1 × 10⁻5, סף אי-שוויון קישור r2 < 0.001, חלון גושים של 10,000 קילובייט, וסטטיסטיקת F > 20. הסטטיסטיקה F חושבה עבור כל משתנה אינסטרומנטלי שנשמר כבטא2/se2 להערכת חוזק המכשיר. SNPs עם סטטיסטיקה של F ≤ 20 הוצאו מהניתוחים במורד הזרם.
לאחר בחירת SNP, מערכי הנתונים של חשיפה ותוצאה אוחדו כדי ליישר את האללים של האפקט. במהלך ההרמוניזציה, אללי ההשפעה ואללים אחרים היו מיושרים בין מערכי הנתונים של החשיפה לתוצאה. SNPs עם אללים לא תואמים הוסרו, ו-SNPs פלינדרומיים עם תדירויות אללים לא ברורות הוצאו כאשר לא ניתן היה לקבוע את כיוון הגדילים. לאחר ההרמוניזציה, נבדקו SNPs שנשמרו כדי לאשר שמקדמי בטא תואמים לאותו אלל אפקט בשני מערכי הנתונים. מספר ה-SNPs שנשמרו לאחר התלכדות והרמוניזציה נרשם בכל חשיפה כנקודת ביקורת ביניים לשחזוריות.
חמש שיטות רנדומיזציה מנדליאניות יושמו, כולל MR-Egger, חציון משוקלל, משוקלל שונות הפוכה, מצב פשוט ומצב משוקלל. ההטרוגניות הפוטנציאלית של משתנים אינסטרומנטליים הוערכה באמצעות מבחן Q של קוקרן, ו-P < 0.05 נחשב לאינדיקטור להטרוגניות משמעותית. הפליאוטרופיה האופקית הפוטנציאלית הוערכה באמצעות יירוט MR-Egger, ו-P < 0.05 הצביע על פליאוטרופיה, מה שמרמז על ירידה באמינות ההערכה הסיבתית47. ניתוחים אלו יצרו הערכות השפעות סיבתיות יחד עם סטטיסטיקות הטרוגניות ופליאוטרופיות לכל חשיפה.
נקודות ביקורת ביניים לשכפול
נקודות ביקורת ביניים נרשמו לאחר כל שלב עיבוד עיקרי כדי להבטיח שחזור זרימת עבודה. בתהליך העבודה של FAERS, נקודות ביקורת כללו את מספר רשומות ה-DEMO שיובאות, מספר הרשומות הייחודיות לאחר הסרת כפילות, מספר הדיווחים הכוללים שבר בירך כאירוע לוואי היעד, מספר הדיווחים על חמשת מעכבי TNF, והמספר הסופי של דוחות זכאים שבהם נרשמו מעכבי TNF כתרופות חשודות עיקריות. בתהליך העבודה של רנדומיזציה מנדליאנית, נקודות בדיקה כללו את מספר ה-SNPs שהופקו בכל חשיפה, מספר SNPs שנשמרו לאחר אי-שיווי משקל קישורי, מספר SNPs הזמינים במאגר הנתונים של התוצאות, מספר SNPs שנשמרו לאחר ההרמוניזציה, ומספר המשתנים האינסטרומנטליים הסופיים ששימשו בכל ניתוח אקראי מנדלי.
דיווח סטטיסטי
דווחו תוצאות רציפות עם הערכות אפקט מתאימות, רווחי ביטחון של 95% (רווח בר-סמך 95%) וערכי P . אלא אם צוין אחרת, המובהקות הסטטיסטית הוגדרה כ-P דו-צדדי < 0.05. לניתוח פרמקוביז'ילנס, דווחו ספירות תיאוריות והערכות אי-פרופורציונליות עבור כל מעכב TNF בודד ועבור קבוצת מעכבי TNF המאוחדת. לניתוח אקראיות מנדליאנית, דווחו הערכות סיבתיות, שגיאות סטנדרטיות, רווחי ביטחון של 95%, ערכי P , סטטיסטיקות הטרוגניות, תוצאות בדיקות פליאוטרופיה, ומספר SNPs שנשמרו בכל חשיפה.