$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
קיבלנו אישור אתי והסכמה מדעת מוועדת האתיקה למחקר ביו-רפואי של בית החולים הראשון המסוים באוניברסיטת נאנצ'אנג. מספר אתיקה: (2025)CDYFYYLK(08-007).
ניתוח MR
שליפת נתונים
נתוני pQTL הפלזמה התקבלו ממחקר של Zheng ואח' 14, ששילב חמישה מערכי GWAS 15,16,17,18,19, וממחקר של פרקינגסטאד ואחרים. קריטריוני ההכללה לנתונים היו: (i) אסוציאציות מובהקות כלל-גנומיות (p < 5 × 10⁻⁸); ו-(ii) חלבוני פלזמה כמטרות טיפוליות פוטנציאליות ל-OA. עיצוב המחקר מסוכם באיור 1. ראשית, זיהינו יעדים טיפוליים מועמדים באמצעות נתוני GWAS מ-IEU OpenGWAS ונתוני pQTL פלזמה ממחקרים של Zheng14 ו-Ferkingstad20 (טבלה משלימה S1 וטבלה משלימה S2). לאחר מכן בוצעו סינון סטייגר וסריקת פנוטיפ כדי לאמת את עמידות התוצאות. ה-IEU OpenGWAS (https://gwas.mrcieu.ac.uk/) שימש לקבלת סטטיסטיקות סיכום עבור OA בירך או ברך (n = 417,596), OA בברך (n = 403,124), ו-OA ירך (n = 393,873)21.
פקודות סינון SNP
SNPs עם מובהקות כלל-גנומית (p < 5 × 10⁻⁸) עברו תהליך התגבשות (r² < 0.001, F-statistics > 10, גודל חלון = 10,000 קילובייט) לפני ניתוח MR.
ניתוח MR
כדי לחקור מטרות תרופות פוטנציאליות, בוצע ניתוח MR באמצעות חלבוני פלזמה כחשיפות ו-OA כתוצאה, שיושם באמצעות חבילת "TwoSampleMR" ב-R (v4.3.1). כאשר היה רק pQTL יחיד זמין לחלבון, נעשה שימוש ביחס וולד; אחרת, הוחל MR משוקלל השונות ההפוכה (MR-IVW), ואחריו הערכות הטרוגניות ופליאוטרופיה. תיקון בונפרוני שימש לביצוע בדיקות מרובות, עם סף של p < 5.63 × 10⁻⁵ לתעדוף חלבונים.
סינון סטייגר וסריקת פנוטיפ
כדי להעריך סיבתיות הפוכה, ביצענו סינון שטייגר. תוצאה של "TRUE" עם p < 0.05 הראתה שאין סיבתיות הפוכה. סריקת פנוטיפ בוצעה באמצעות LDtrait (https://ldlink.nih.gov/?tab=ldtrait#home-tab)22 כדי לבחון את הקשרים בין pQTLs לתכונות אחרות. הספים היו R² = 0.1 וחלון זוגות בסיס של ±500,000. השפעות פליאוטרופיות הוקצו ל-pQTLs שעמדו בשני התנאים הבאים: (i) אסוציאציה מובהקת כלל-גנומית (p < 5 × 10⁻⁸), ו-(ii) קשר עם גורמי סיכון ידועים ל-OA.
מחקר אסוציאציה רחב-תופעה
כדי לקחת בחשבון פליאוטרופיה גנטית והשפעות מחוץ למטרה, ערכנו מחקר אסוציאציה רחב-פנומן (PheWAS) באמצעות פורטל AstraZeneca PheWAS (https://azphewas.com/), שמכיל 15,500 פנוטיפים בינאריים ו-1,500 פנוטיפים רציפים מ-~450,000 משתתפי Biobank בבריטניה23. הספים הוגדרו לערכי ברירת מחדל כדי למזער חיוביים שגויים.
רשת אינטראקציה בין חלבון לחלבון (PPI)
כדי להמחיש אינטראקציות בין מטרות חלבון פוטנציאליות שזוהו על ידי MR, השתמשנו ב-GeneMANIA (https://genemania.org/) לניתוח אינטראקציה בין חלבון לחלבון ולהצגת תוצאות24.
ניתוח העשרה
כדי לחקור רלוונטיות ביולוגית, ביצענו ניתוח העשרה באמצעות כלי ביו-אינפורמטיקה מ-https://www.bioinformatics.com.cn לניתוח נתונים והדמיה.
תהליך עבודה טרנסקריפטומי
סך ה-RNA הופק באמצעות ערכת התגובה להפקת RNA בהתאם להנחיות היצרן. איכות ה-RNA הוערכה באמצעות מערכת הערכת איכות RNA אוטומטית; נעשה שימוש רק בדגימות עם RIN ≥7.0. האיכות אושרה על ידי אלקטרופורזה של ג'ל אגרוז ללא RNase (1.5% ג'ל). mRNA אאוקריוטי הועשר באמצעות חרוזי Oligo(dT); mRNA פרוקריוטי הועשר באמצעות ערכת מגנטית להסרת RNA. mRNA פורק (200-700 נטרון) והוסב ל-cDNA באמצעות ערכת ההכנה של ספריית ה-RNA. ספריית ה-cDNA תוקנה בקצה, עם זנב A, קשורה למתאמים, טוהרה באמצעות חרוזים מגנטיים מטהרים DNA (1.0×), והוגברה ב-PCR. הרצף בוצע על פלטפורמת רצף מהדור הבא עם תפוקה גבוהה. גנים שבוטאו באופן שונה הוגדרו על ידי log₂FC > 1 והתאימו p < 0.05.
פרמקולוגיה רשת
כדי לזהות תרופות פוטנציאליות לחלבוני מטרה, השתמשנו ב-BATMAN-TCM (http://bionet.ncpsb.org.cn/batman-tcm/index.php)25. סף ציון של 0.74 (LR = 32.5) שימש לבחירת תרכובות ידועות וחזויות. רכיבים צמחיים נאספו מ-TCMSP (https://old.tcmsp-e.com/index.php) וסוננו עם OB > 30% ו-DL > 0.1826.
עגינה מולקולרית
עגינה מולקולרית שימשה להערכת אינטראקציות קשירה. מבני חלבונים נאספו מה-PDB (https://www.rcsb.org/). UCSF Chimera שימשה לעיבוד מוקדם של מבנים על ידי הסרת ליגנדים וממסים. כלי AutoDock שימשו לחישוב מטעני גסטייגר ולהגדרת מרכזי קופסאות וגדלים. מבני התרופות התקבלו מ-PubChem (https://pubchem.ncbi.nlm.nih.gov/) ועובדו מראש באופן דומה. העגינה בוצעה באמצעות AutoDock Vina. מידות הקופסה משתנות לפי המטרה. חושבו זיקות קשירה, והתוצאות הוצגו בכימרה של UCSF.