פרוטוקול זה מתאר תהליך עבודה בר-שחזור של רנדומיזציה מנדלית (Mendelian randomization) בעזרת שני מדגמים, להערכת הקשר הסיבתי הפוטנציאלי בין שבץ איסכמי לדמנציה וסקולרית, תוך שימוש בסטטיסטיקות סיכום של מחקרים לאסוציאציה גנומית רחבה (GWAS) הזמינות לציבור.
מאמר שיטה
קשר סיבתי בין שבץ איסכמי לדמנציה וסקולרית: מחקר רנדומיזציה מנדלית
26 צפיות
⸱
11 בספטמבר 2026
במאמר זה
סיכום
תקציר
שבץ איסכמי (IS) הוא גורם מרכזי לנכות ותמותה ברחבי העולם, ודמנציה וסקולרית (VaD) היא תת-סוג נפוץ של דמנציה הקשור לפגיעה צרברו-וסקולרית. מחקרים תצפיתיים העלו קשר בין IS ל-VaD, אך מחקרים אלו חשופים להשפעות של משתנים מתערבים ולסיבתיות הפוכה. פרוטוקול זה מתאר זרימת עבודה הדירה של רנדומיזציה מנדלית (MR) לשני מדגמים להערכת הקשר הסיבתי הפוטנציאלי בין IS ל-VaD, תוך שימוש בסטטיסטיקות סיכום של מחקרי אסוציאציה גנומיים (GWAS) הזמינים לציבור. כלי גנטיים הקשורים ל-IS הופקו ממאגר נתוני GWAS ציבורי, ואסוציאציות של התוצאה עבור VaD הופקו ממאגר נתוני GWAS ציבורי של VaD. מזהי מערכי הנתונים המתאימים מופיעים בסעיף הפרוטוקול. לאחר התאמת התוצאות והרמוניזציה של האללים, 51 פולימורפיזמיים של נוקליאוטיד יחיד (SNPs) נשמרו לניתוח ה-MR הסופי. זרימת העבודה כוללת בחירת משתנה מכשיר, קיבץ (clumping) של חוסר שיווי משקל של קישור (linkage disequilibrium), הרמוניזציה של אללים, הערכת עוצמת המכשיר, ניתוח שקלול שונות הפוך (IVW), ניתוח חציון משוקלל, ניתוח MR-Egger, בדיקת הומוגניות, הערכת פליאוטרופיה אופקית וניתוח רגישות של "השמטת אחד" (leave-one-out). בניתוח המייצג, שיטת IVW הראתה אסוציאציה חיובית בין IS שנחזה גנטית לבין הסיכון ל-VaD, ושיטת החציון המשוקלל הניבה תוצאה תואמת בכיוונה. האומדן של MR-Egger היה עקבי בכיוונו אך לא הגיע למשמעות סטטיסטית. לפיכך, יש לפרש ממצאים אלו כראיות מרמזות על השפעה סיבתית אפשרית, ולא כהוכחה מוחלטת לסיבתיות. פרוטוקול זה עשוי לסייע לחוקרים ליישם זרימת עבודה שקופה והדירה של MR לבחינת תוצאות הקשורות למחלות צרברו-וסקולריות באמצעות נתוני GWAS ציבוריים.
מבוא
שבץ מוחי איסכמי (IS) הוא גורם מוביל לנכות חמורה ולמוות ברחבי העולם1. על פי הנתונים האחרונים ממחקר נטל המחלה העולמי (Global Burden of Disease Study), השכיחות העולמית של IS ירדה, אך היא נותרה נטל משמעותי, במיוחד במזרח אירופה, מזרח אסיה, מרכז אסיה ובאפריקה ממסביב לסהרה2,3. יתרה מכך, השכיחות של IS נמצאת במגמת עלייה במדינות בעלות מדד סוציו-דמוגרפי (SDI) נמוך, והיא צפויה להמשיך לעלות בין השנים 2020 ל-2030.
מחקרים אפידמיולוגיים רבים הצביעו על קשר פוטנציאלי בין מחלות צרברו-וסקולריות לבין דמנציה וסקולרית (VaD)4,5. בסין בפרט, VaD היא אחד מסוגי הדמנציה השכיחים ביותר בקרב קשישים, והשכיחות הגבוהה של IS קשורה קשר הדוק להופעת VaD. עם זאת, מרבית המחקרים הקיימים הם תצפיתיים, וייתכן שבלבול שיורי (residual confounding), סיבתיות הפוכה והבדלים בהגדרת המחלה מגבילים את הפרשנות הסיבתית. לכן, נדרשות גישות אנליטיות נוספות כדי להעריך האם IS מעורבת פוטנציאלית בהתפתחות של VaD.
שבץ מוחי איסכמי (IS) עשוי להשפיע על הופעת דמנציה וסקולרית (VaD) באמצעות מספר מנגנונים. נזק קרביו-וסקולרי הנגרם כתוצאה מ-IS עלול להוביל לפרפוזיה מוחית לא מספקת ולשינויים נוירודגנרטיביים עוקבים6. בנוסף, VaD עשויה להיות קשורה לפגיעה קרביו-וסקולרית כרונית, נזק מיקרו-וסקולרי ותגובות דלקתיות לאחר IS7,8. לכן, חקירת הקשר בין IS ל-VaD היא חשובה להבנת הקשרים הפוטנציאליים בין המחלות ולהכוונת אסטרטגיות מניעה עתידיות.
רנדומיזציה מנדלית (MR) היא שיטה אנליטית המשתמשת בווריאנטים גנטיים כמשתנים אינסטרומנטליים להערכת קשרי סיבתיות פוטנציאליים בין חשיפה לתוצאה9. בהשוואה למחקרים תצפיתיים קונבנציונליים, MR יכולה להפחית הטיות הנובעות ממשתנים מתערבים ומסיבתיות הפוכה, תחת הנחות תקפות למשתנים האינסטרומנטליים. בהשוואה לניסויים מבוקרים אקראיים או למודלים ניסויים, MR יכולה להשתמש בסטטיסטיקות סיכום קיימות ממחקרי אסוציאציה גנומיים (GWAS) כדי להעריך קשרים בעלי רלוונטיות קלינית כאשר מחקרי התערבות ישירים אינם ישמיים, אינם אתיים או קשים לביצוע. מחקרים גנטיים אחרונים שהשתמשו בגישות של רנדומיזציה מנדלית בחנו את הקשר בין IS ל-VaD10. עם זאת, תיאורים מפורטים ברמת הפרוטוקול של בחירת האינסטרומנטים לחשיפה, חילוץ נתוני התוצאה, הרמוניזציה של אללים, הערכת עוצמת האינסטרומנט, בדיקת ההטרוגניות, הערכת פליאוטרופיה ואנליזות רגישות, נותרים חשובים לצורך יישום הניתן לשחזור.
פרוטוקול זה מתאים כאשר קיימים נתוני סיכום (summary statistics) פומייים של GWAS הן עבור החשיפה והן עבור התוצאה הנבדקת, כאשר ניתן לבחור כלים גנטיים (genetic instruments) תקפים עבור החשיפה, וכאשר שאלת המחקר שואפת להעריך קשר סיבתי פוטנציאלי ולא לבסס ישירות מנגנונים ביולוגיים. לפיכך, מחקר זה מציג פרוטוקול MR של שני דגמים להערכת הקשר הסיבתי הפוטנציאלי בין IS ל-VaD באמצעות נתוני סיכום פומייים של GWAS. המטרה היא לספק זרימת עבודה אנליטית הניתנת לשחזור, ולא לבסס מנגנונים ביולוגיים סופיים. פרוטוקול זה מתאר בחירת כלים גנטיים, התאמת תוצאות, הרמוניזציה של אללים, הערכת MR ואנליזות רגישות, ומשתמש בקשר שבין IS ל-VaD כדוגמה מייצגת.
פרוטוקול
מחקר זה השתמש בנתונים סטטיסטיים מסוכמים של GWAS הזמינים לציבור. המחקרים המקוריים קיבלו אישור מוועדה לאתיקה מוסדית והסכמה מדעת מהמשתתפים. לא נדרש אישור אתי נוסף עבור ניתוח משני זה.
1. מערך המחקר והצהרת אתיקה
- השתמשו בתכנון MR בעל שני מדגמים כדי להעריך את הקשר הסיבתי הפוטנציאלי בין IS ל-VaD.
- השתמשו אך ורק בנתוני סיכום של GWAS הזמינים לציבור. מכיוון שכל המחקרים המקוריים שנכללו במאגרי הנתונים של ה-GWAS קיבלו אישור מוועדה לאתיקה מוסדית והסכמה מושכלת מהמשתתפים, אין צורך באישור אתי נוסף עבור ניתוח משני זה.
- הגדירו את IS כחשיפה (exposure) ואת VaD כתוצאה (outcome) לפני שליפת הנתונים.
2. הכנת התוכנה
- בצע את כל הניתוחים באמצעות תוכנת הסטטיסטיקה R.
- טען את חבילת TwoSampleMR לצורך חילוץ כלים (instruments), חילוץ נתוני תוצאה, הרמוניזציה של אללים, הערכת MR, בדיקת הטרוגניות וניתוחי רגישות. טען את חבילת MRPRESSO להערכת פליאוטרופיה אופקית גלובלית וזיהוי וריאנטים חריגים פוטנציאליים.
- השתמש בפונקציות הליבה הבאות במסגרת זרימת העבודה: extract_instruments, extract_outcome_data, harmonise_data, mr, generate_odds_ratios, mr_heterogeneity, mr_pleiotropy_test, mr_singlesnp, mr_leaveoneout, mr_scatter_plot, mr_forest_plot, mr_funnel_plot, ו-mr_leaveoneout_plot. השתמש ב-MRPRESSO::mr_presso לביצוע בדיקת MR-PRESSO גלובלית, ובמידת הצורך, לבדיקות של ערכים חריגים ועיוותים.
3. בחירת מערך נתוני החשיפה וסינון משתנים אינסטרומנטליים
- חלץ פולימורפיזמיים של נוקלאיד יחיד (SNPs) הקשורים ל-IS מתוך מערך הנתונים IEU OpenGWAS ebi-a-GCST90018864. מערך נתונים זה כולל 11,929 מקרי IS ו-472,192 ביקורות, עם גודל מדגם כולל של 484,121 משתתפים.
- בחר SNPs מועמדים הקשורים ל-IS באמצעות סף מובהקות של p < 5 × 10-6.
- סף ה- p < 5 × 10-6 משמש בפרוטוקול מייצג זה כדי לשמר מספר מספיק של מכשירים (instruments) בלתי תלויים עבור אנליזות רגישות. מכיוון שסף זה פחות מחמיר מסף המובהקות המקובל ברמת הגנום, יש לפרש ממצאים המבוססים על מכשירים אלו כחקרים ראשוניים ומצביעים ולא כסופיים.
- בצע קיבוב (clumping) של חוסר שיווי משקל במקומי (linkage disequilibrium) כדי לשמר SNPs בלתי תלויים. השתמש בסף r2 של 0.001 ובחלון קיבוב של 10,000 kb.
- השתמש בפקודת R הבאה כדי לחלץ מכשירים מועמדים:
extract_instruments(outcomes = "ebi-a-GCST90018864", p1 = 5e-06, clump = TRUE) - חשב את סטטיסטיקת F עבור כל SNP שנשמר באמצעות הנוסחה F = (β/SE)2, כאשר β ו-SE הם אומדן ההשפעה של ה-SNP על החשיפה והטעות הסטנדרטית, בהתאמה.
- החרג SNPs עם סטטיסטיקת F < 10 כדי להפחית את הסיכון להטיה של מכשיר חלש. סכם את התפלגות סטטיסטיקות F באמצעות המינימום, המקסימום, החציון והטווח הבין-רבעוני (IQR).
4. חילוץ סט נתוני התוצאות
- חלץ קשרי תוצאה עבור VaD מתוך מערך נתוני ה-GWAS של FinnGen המסומן כ-finn-b-F5_VASCDEM. מערך נתונים זה כולל 881 מקרי VaD ו-211,508 ביקרים.
- חלץ את קשרי התוצאה עבור כל ה-SNPs שנבחרו הקשורים ל-IS באמצעות הפונקציה extract_outcome_data.
- חפש proxy SNPs כאשר ה-SNPs המטרה אינם זמינים באופן ישיר במערך נתוני התוצאה. שמור וריאנטים התואמים ל-proxy כאשר מידע תקף על ה-proxy והתאמה של אללים זמינים.
- השתמש בפקודת R הבאה כדי לחלץ את מערך נתוני התוצאה:
extract_outcome_data(snps = exp_dat$SNP, outcomes = "finn-b-F5_VASCDEM", access_token = NULL)
5. הרמוניזציה של אללים ושימור SNP
- יש להרמוניזציה (Harmonize) את מערכי הנתונים של החשיפה והתוצאה לפני ניתוח MR.
- יישרו את אלל האפקט (effect allele) ואת האלל השני כך שהערכות האפקט של SNP–exposure ו-SNP–outcome יתאימו לאותו אלל אפקט.
- יש לבחון את ה-SNPs במהלך ההרמוניזציה כדי לזהות חוסר התאמה באללים, קשרי תוצאה (outcome associations) שאינם זמינים ויישור אללים מעורפל. יש להשאיר SNPs עם הערכות השפעה זמינות עבור החשיפה (exposure) והתוצאה לאחר ההרמוניזציה.
- השתמש בפקודת R הבאה לצורך הרמוניזציה של אללים: harmonise_data(exposure_dat = exp_dat, outcome_dat = out_dat, action = 1)
- יש לתעד את מספר ה-SNPs שחולצו בתחילה ממאגר הנתונים של החשיפה, את מספר ה-SNPs עם קשרי תוצאה זמינים, את מספר הווריאנטים שהותאמו כפרוקסי (proxy-matched variants), ואת מספר ה-SNPs שנשמרו לניתוח ה-MR הסופי.
6. ניתוח MR ראשוני וניתוחי רגישות
- בצע את ניתוח ה-MR הראשוני באמצעות שיטת השקלול של שונות הפוכה (IVW).
- בצע ניתוחי MR משלימים באמצעות שיטות החציוני המשוקלל (weighted median) ו-MR-Egger כדי להעריך את העקביות של האומדן הסיבתי תחת הנחות שונות.
- הפק יחסי סיכויים (ORs) ורווחי סמך (CIs) של 95% מהאומדני של ה-MR.
- השתמש בפקודות R הבאות כדי להריץ את ניתוח ה-MR ולחשב את ה-ORs:
res <- mr(dat)
res <- generate_odds_ratios(res) - מכיוון שמוערך קשר ראשוני אחד שנקבע מראש בין החשיפה לתוצאה, הגדר את סף תיקון בונפרוני (Bonferroni-corrected threshold) עבור הסקת ה-MR הראשונית כ-0.05/1 = 0.05. השתמש בערך p דו-צדדי < 0.05 כסף המובהקות הסטטיסטית עבור ניתוח ה-IVW הראשוני.
- כדי להעניק הקשר לדיוק הסטטיסטי של הניתוח הראשוני בהינתן המספר המוגבל של מקרי VaD, בצע חישוב עוצמה (power calculation) של מבחן Wald דו-צדדי בדיעבד (post hoc), המבוסס על אומדן ה-IVW הנצפה וטעות התקן שלו. חשב את העוצמה באמצעות פרמטר אי-מרכזיות βIVW/SEIVW ברמת α דו-צדדית של 0.05. חשב את האפקט המינימלי הניתן לזיהוי התואם לעוצמה של 80% כ-exp[(z0.975 + z0.80) × SE]. מכיוון ש-IS היא חשיפה בינארית וקשרי SNP–exposure מבוטאים בסולם log-odds, דווח על ניתוח זה כהערכה תיאורית של הדיוק המותנית באפקט, ולא כחישוב גודל מדגם שנקבע מראש11.
7. ניתוחי הטרוגניות, پلیאוטרופיה אופקית ואנליזות השמטת-דגימה (leave-one-out)
- העריכו את ההומוגניות בין ה-SNP באמצעות סטטיסטיקת Cochran’s Q.
- השתמשו בפקודת ה-R הבאה כדי לחשב סטטיסטיקות הומוגניות:
mr_heterogeneity(dat) - העריכו פליוטרופיה אופקית כוונית באמצעות מבחן ה-intercept של MR-Egger.
- השתמשו בפקודת ה-R הבאה כדי לחשב את ה-intercept של MR-Egger:
mr_pleiotropy_test(dat) - העריכו פליוטרופיה אופקית גלובליתת ווריאנטים חריגים פוטנציאליים באמצעות MR-PRESSO עם 10,000 סימולציות, סף מובהקות של 0.05, כאשר הן מבחני החריגים (outlier) והן מבחני העיוות (distortion) מופעלים12.
- השתמשו בפקודת ה-R הבאה לניתוח MR-PRESSO:
set.seed(20260729)
mr_presso(
BetaOutcome = "beta.outcome",
BetaExposure = "beta.exposure",
SdOutcome = "se.outcome",
SdExposure = "se.exposure",
OUTLIERtest = TRUE,
DISTORTIONtest = TRUE,
data = dat,
NbDistribution = 10000,
SignifThreshold = 0.05
) - בצעו ניתוח single-SNP כדי להעריך את ההשפעה של כל SNP בנפרד.
- השתמשו בפקודת ה-R הבאה לניתוח single-SNP:
res_single <- mr_singlesnp(dat) - בצעו ניתוח leave-one-out על ידי הסרה רצופית של SNP אחד בכל פעם וחזרה על ניתוח ה-IVW.
- השתמשו בפקודת ה-R הבאה לניתוח leave-one-out:
res_loo <- mr_leaveoneout(dat)
8. ויזואליזציה ופלט
- הפיקו תרשים פיזור (scatter plot) כדי להציג את הערכות ה-MR בשיטות שונות.
- הפיקו תרשים יער (forest plot) כדי להציג את ההערכות עבור SNP בודד ואת הערכת ה-MR הכוללת.
- הפיקו תרשים משפך (funnel plot) כדי להעריך באופן חזותי את הסימטריה של ההערכות הספציפיות ל-SNP.
- הפיקו תרשים השמטת-אחד (leave-one-out plot) כדי לבחון האם ההערכה הכוללת מושפעת מ-SNP בודד.
- השתמשו בפקודות R הבאות כדי להפיק את התרשימים:
mr_scatter_plot(res, dat)
mr_forest_plot(res_single)
mr_funnel_plot(res_single)
mr_leaveoneout_plot(res_loo) - ייצאו את תוצאות ה-MR, את טבלת המשתנים האינסטרומנטליים, את טבלת הקשר לתוצאה (outcome association) ואת פלטי ניתוח הרגישות לצורך דיווח ושחזור.
תוצאות
איור 1 מסכם את זרימת העבודה הכוללת של פרוטוקול הרנדומיזציה המנדלית לשני מדגמים, החל מבחירת מערכי נתונים וסינון כלים, דרך ניתוחי רגישות ועד להדמיית התוצאות. תוצרים מייצגים מכל שלב מוצגים להלן.
משתנים אינסטרומנטליים גנטיים ובקרת איכות של הרמוניזציה
בסך הכל 52 SNPs הקשורים ל-IS הופקו תחילה ממאגר נתוני ה-GWAS של החשיפה באמצעות סף מובהקות של p < 5 × 10-6. לאחר מכן, הופקו קשרי תוצאה ממאגר נתוני ה-GWAS של VaD. לאחר התאמת תוצאות, חיפוש פרוקסי (proxy) והרמוניזציה של אללים, נשמרו 51 SNPs לניתוח ה-MR הסופי. שני וריאנטים הותאמו באמצעות proxy SNPs. כלי המדידה הגנטיים שנשמרו מופיעים ב-טבלה 1. סטטיסטיקות ה-F של 51 הכלים שנשמרו נעו בין 20.91 ל-46.06, עם חציון של 25.08 ו-IQR של 22.95–30.71. לאף SNP שנשמר לא היה סטטיסטיקת F נמוכה מ-10, מה שמעיד על סבירות נמוכה להטיה של כלי מדידה חלשים.
הערכות MR עבור IS ו-VaD
תרשים הפיזור הציג את הכיוון ואת העוצמה של אומדני ה-MR הספציפיים ל-SNP בשיטות MR שונות (איור 2). שיטות ה-IVW והחציון המשוקלל הראו קשרים חיוביים, בעוד שאומדן ה-MR-Egger היה עקבי מבחינת הכיוון אך לא הגיע למשמעות סטטיסטית.
תרשים היער (forest plot) של אומדני ה-SNP הבודדים הראה את הקשרים הספציפיים ל-SNP עם הסיכון ל-VaD (איור 3). אומדני ה-MR הכלליים מסוכמים ב-איור 4. שיטת ה-IVW הראתה קשר חיובי מובהק סטטיסטית בין IS שנחזה גנטית לבין הסיכון ל-VaD (OR = 1.63, 95% CI: 1.21–2.21, p = 0.0013). שיטת החציון המשוקלל (weighted median) הניבה אומדן מובהק ותואם (OR = 1.60, 95% CI: 1.06–2.42, p = 0.0240). אומדן ה-MR-Egger היה עקבי בכיוונו אך לא הגיע למובהקות סטטיסטית (OR = 2.21, 95% CI: 0.97–5.01, p = 0.0645).
מחקר ה-GWAS לתוצאות VaD של FinnGen כלל 881 חולים ו-211,508 ביקורות. בחישוב בדיקת Wald דו-כיוונית לאחר המאבד (post hoc), המותנית בהערכת ה-IVW שנצפתה, העוצמה המוערכת עבור ניתוח ה-IVW הראשוני הייתה 89.5% ב-α = 0.05. האפקט המינימלי שניתן לזיהוי עבור עוצמה של 80% היה OR של 1.54. לעומת זאת, האפקט המינימלי שניתן לזיהוי עבור ניתוח ה-MR-Egger היה OR של 3.23, אשר עלה על הערכת ה-MR-Egger העקבית כיוונית שנצפתה (OR = 2.21). לפיכך, יש לפרש את תוצאת ה-MR-Egger שאינה מובהקת כמשקפת את הדיוק המוגבל של מעריך רגישות פחות יעיל זה, ולא כראיה כשלעצמה נגד כיוון הערכת ה-IVW הראשונית.
כיוון שבוחנה קשר ראשוני אחד שנקבע מראש בין חשיפה לתוצאה, סף תיקון בונפרוני (Bonferroni) עבור הסקת ה-MR הראשונית היה 0.05/1 = 0.05. לפיכך, תוצאת ה-IVW עמדה בסף המובהקות המתוקן. יחד, ממצאים אלו מספקים ראיות מרמזות לאפקט חיובי אפשרי של IS שנחזה גנטית על הסיכון ל-VaD, בעיקר בהתבסס על אומדן ה-IVW ובסיוע של ניתוח הרגישות של החציוני המשוקלל (weighted median sensitivity analysis). עם זאת, אומדן ה-MR-Egger היה עקבי מבחינת הכיוון אך לא הגיע למובהקות סטטיסטית; לפיכך, אין לפרש את התוצאות כראיה חלוטה לקשרי סיבתיות.
ניתוחי הטרוגניות ופלייוטרופיה אופקית
הטרוגניות בין-SNP הוערכה באמצעות סטטיסטיקת Q של Cochran. בדיקת ההטרוגניות של IVW הניבה Q = 57.46 עם 50 דרגות חופש (p = 0.218), ובדיקת ההטרוגניות של MR-Egger הניבה Q = 56.77 עם 49 דרגות חופש (p = 0.208). תוצאות אלו לא הצביעו על הטרוגניות משמעותית בין האומדנים הספציפיים ל-SNP. פלייוטרופיה אופקית כיוונית הוערכה באמצעות בדיקת ה-intercept של MR-Egger. ה-intercept של MR-Egger היה -0.0195 (SE = 0.0253, p = 0.444), דבר המעיד על היעדר ראיה סטטיסטית לפלייוטרופיה אופקית כיוונית. ניתוח MR-PRESSO בוצע באמצעות 51 ה-SNPs ההרמוניים שנשמרו עם 10,000 סימולציות. הבדיקה הגלובלית של MR-PRESSO לא הראתה ראיה לפלייוטרופיה אופקית גלובלית (RSSobs = 59.61; p אמפירי = 0.2388). מכיוון שהבדיקה הגלובלית לא הייתה מובהקת סטטיסטית, בדיקות החריגים והעיוותים האינדיבידואליות לא היו רלוונטיות, ולא הופק אומדן מתוקן לחריגים. תרשים המשפך (funnel plot) סיפק הערכה ויזואלית של הסימטריה של האומדנים הספציפיים ל-SNP (איור 5). הגדרות ותוצאות מפורטות של MR-PRESSO מדווחות ב-טבלה משלימה 1.
ניתוח רגישות מסוג השארת-איבר-אחד (Leave-one-out) ותיקוף פלט של הפרוטוקול
ניתוח השמטת-אחד (Leave-one-out analysis) בוצע כדי להעריך האם הערכת ה-MR הכוללת הושפעה מ-SNP בודד כלשהו. תרשים ההשמטת-אחד הראה כי הסרה רצופה של SNPs בודדים לא שינתה באופן מהותי את ההערכה הכוללת (איור 6), מה שמרמז כי אף כלי גנטי בודד לא שלט בקשר. יחד, תוצאות מייצגות אלו מדגימות את התוצר המעשי של זרימת העבודה של MR לשני מדגמים המתוארת בפרוטוקול זה. טבלה 1 מציג את המכשירים הגנטיים שנשמרו לאחר בחירת SNP והרמוניזציה. איור 2 מדגים את כיוון אומדני ה-MR על פני שיטות שונות, איור 3 מציג הערכות ספציפיות ל-SNP, איור 4 מסכם את הערכות ה-MR הכוללות, איור 5 תומך בהערכה חזותית של פליאוטרופיה או אסימטריה, ו- איור 6 מעריך את השפעתם של SNP בודדים. טבלה משלימה 2 מקשר בין כל שלב מרכזי בפרוטוקול לבין התוצר המתאיים המאמת אותו.
כל הסטטיסטיקות המסכמות הגולמיות של GWAS שנותחו במחקר זה זמינות לציבור. מערך הנתונים של החשיפה לשבץ מוחי איסכמי התקבל ממאגר הנתונים IEU OpenGWAS תחת מזהה מערך הנתונים ebi-a-GCST90018864. מערך הנתונים של תוצאת הדמנציה וסקולרית התקבל ממערך הנתונים FinnGen GWAS תחת מזהה מערך הנתונים finn-b-F5_VASCDEM. טבלת המשתנה הכלי השליפה, מערך הנתונים המהורמוני לניתוח, טבלאות תוצאות ה-MR וקבצי פלט רלוונטיים מסופקים כקבצים משלימים.
איור 1. זרימת העבודה של פרוטוקול MR בעל שני מדגמים. איור זה מסכם את השלבים המרכזיים של הפרוטוקול, הכוללים בחירת סט נתוני חשיפה, סינון משתנים אינסטרומנטליים, חילוץ נתוני תוצאה, הרמוניזציה של אללים, הערכת MR, בדיקת הומוגניות, הערכת פליאוטרופיה אופקית, ניתוח leave-one-out והדמיה. אנא לחצו כאן כדי לצפות בגרסה גדולה יותר של איור זה.
איור 2. תרשים פיזור של הערכות MR לקשר שבין IS ל-VaD. כל נקודה מייצגת הערכה ספציפית ל-SNP. הקווים המותאמים מייצגים את הקשר המוערך שהתקבל באמצעות שיטות MR שונות. אנא לחצו כאן כדי לצפות בגרסה מורחבת של איור זה.
איור 3. תרשים פורסט (Forest plot) של אומדני MR עבור SNP בודד. איור זה מציג את האומדנים האינדיבידואליים הספציפיים ל-SNP עבור הקשר בין וריאנטים גנטיים הקשורים ל-IS לבין סיכון ל-VaD. קווים אופקיים מייצגים רווחי סמך של 95%. אנא לחץ כאן כדי לצפות בגרסה גדולה יותר של איור זה.
איור 4. Forest plot של אומדני MR כוללים בשיטות שונות. איור זה מסכם את אומדני ה-MR הכוללים שהתקבלו באמצעות שיטות IVW, חציוני משוקלל (weighted median) ו-MR-Egger. קווים אופקיים מייצגים רווחי סמך של 95%. MR, רנדומיזציה מנדלית (Mendelian randomization); IVW, משוקלל שונות הפוכה (inverse variance weighted). אנא לחצו כאן כדי לצפות בגרסה גדולה יותר של איור זה.
איור 5. תרשים משפך של הערכות MR ספציפיות ל-SNP . איור זה מציג את התפלגות הערכות ה-MR הספציפיות ל-SNP ומאפשר הערכה חזותית של הסימטריה בין הכלים הגנטיים. MR, רנדומיזציה מנדלית; IVW, משקל שונות הפוכה. אנא לחצו כאן כדי לצפות בגרסה גדולה יותר של איור זה.
איור 6. ניתוח רגישות מסוג leave-one-out. איור זה מציג את אומדן ה-MR לאחר הסרה רצופית של כל SNP. התרשים שימש להערכה האם האומדן הכללי הושפע ממכשיר גנטי בודד כלשהו. אנא לחצו כאן כדי להציג גרסה מורחבת של איור זה.
| לא. | SNP (פולימורפיזם של נוקליאוטיד בודד) | גן | כרומ. | EA | דלקת מפרקים שחיקתית (OA) | EAF.IS | EAF.VD | IS β (SE) | VD β (טעות תקן - SE) | סטטיסטיקת F | R² |
| 1 | rs10886430 | GRK5 | 10 | ג' | A | 0.118244 | 0.09536 | 0.1147 (0.0220) | -0.0172 (0.0852) | 27.18 | 5.61E-05 |
| 2 | rs10936572 | LOC107986051 | 3 | ט | C | 0.179643 | 0.1235 | -0.0512 (0.0104) | 0.1444 (0.0743) | 24.24 | 5.01E-05 |
| 3 | rs11045239 | PDE3A | 12 | A | חלבון G | 0.492943 | 0.4629 | 0.0604 (0.0089) | 0.0303 (0.0494) | 46.06 | 9.51E-05 |
| 4 | rs11047532 | LOC105369698 | 12 | ג' | C | 0.305214 | 0.1513 | 0.0499 (0.0105) | 0.0159 (0.0689) | 22.59 | 4.66E-05 |
| 5 | rs11065836 | CUX2 | 12 | A | ג' | 0.159123 | 0.07442 | -0.0608 (0.0103) | 0.0245 (0.0929) | 34.84 | 7.20E-05 |
| 6 | rs11105378 | ATP2B1 | 12 | T | C | 0.203466 | 0.07522 | -0.0540 (0.0098) | -0.1174 (0.0916) | 30.36 | 6.27E-05 |
| 7 | rs117140252 | - | 14 | A | ג' | 0.041216 | 0.06988 | 0.1568 (0.0311) | 0.1498 (0.0967) | 25.42 | 5.25E-05 |
| 8 | rs117343276 | - | 10 | G | A | 0.030533 | 0.0431 | -0.0905 (0.0196) | -0.1449 (0.1233) | 21.32 | 4.40E-05 |
| 9 | rs11831940 | HDAC7 | 12 | A | G | 0.251914 | 0.3526 | 0.0631 (0.0130) | 0.1226 (0.0513) | 23.56 | 4.87E-05 |
| 10 | rs11880613 | DNM2 | 19 | A | ג' | 0.184857 | 0.1804 | -0.0613 (0.0110) | 0.1499 (0.0640) | 31.06 | 6.41E-05 |
| 11 | rs12445022 | LOC124903748 | 16 | A | G | 0.271667 | 0.3193 | 0.0605 (0.0112) | 0.0218 (0.0527) | 29.18 | 6.03E-05 |
| 12 | rs12509595 | - | 4 | C | T | 0.297929 | 0.3124 | 0.0577 (0.0091) | 0.1462 (0.0533) | 40.2 | 8.30E-05 |
| 13 | rs12633109 | - | 3 | ט | ג' | 0.360897 | 0.3263 | 0.0436 (0.0089) | 0.0395 (0.0526) | 24 | 4.96E-05 |
| 14 | rs1275980 | KCNK3 | 2 | T | C | 0.498545 | 0.5481 | -0.0582 (0.0094) | -0.0128 (0.0495) | 38.33 | 7.92E-05 |
| 15 | rs13123551 | - | 4 | A | T | 0.616442 | 0.5065 | 0.0552 (0.0097) | 0.0818 (0.0493) | 32.38 | 6.69E-05 |
| 16 | rs147871383 | MIR99AHG | 21 | A | G | 0.02368 | 0.01925 | 0.2259 (0.0494) | 0.2534 (0.1903) | 20.91 | 4.32E-05 |
| 17 | rs16918175 | - | 10 | C | T | 0.104037 | 0.0811 | -0.0649 (0.0131) | 0.0332 (0.0895) | 24.54 | 5.07E-05 |
| 18 | rs17182166 | ACVR1 | 2 | T | G | 0.12494 | 0.137 | 0.0751 (0.0156) | -0.0046 (0.0703) | 23.18 | 4.79E-05 |
| 19 | rs1906779 | - | 15 | A | ג' | 0.211386 | 0.1239 | 0.0438 (0.0095) | -0.0342 (0.0758) | 21.26 | 4.39E-05 |
| 20 | rs1948696 | ITGB5 | 3 | C | T | 0.650138 | 0.6041 | -0.0459 (0.0091) | 0.0219 (0.0501) | 25.44 | 5.25E-05 |
| 21 | rs1973765 | LSP1 | 11 | C | T | 0.486717 | 0.4124 | 0.0441 (0.0092) | -0.0229 (0.0505) | 22.98 | 4.75E-05 |
| 22 | rs2429123 | CACNA1C, DCP1B | 12 | C | T | 0.729926 | 0.725 | 0.0473 (0.0098) | 0.0190 (0.0549) | 23.3 | 4.81E-05 |
| 23 | rs2447561 | - | 8 | T | A | 0.873442 | 0.8382 | -0.0589 (0.0116) | -0.0570 (0.0674) | 25.78 | 5.33E-05 |
| 24 | rs245015 | MSH3 | 5 | A | ג' | 0.685061 | 0.6562 | -0.0456 (0.0089) | 0.0065 (0.0519) | 26.25 | 5.42E-05 |
| 25 | rs2501968 | CENPQ | 6 | ג׳ | A | 0.465302 | 0.4229 | -0.0490 (0.0086) | -0.0879 (0.0498) | 32.46 | 6.71E-05 |
| 26 | rs2526620 | - | 7 | חג | A | 0.252144 | 0.2459 | 0.0500 (0.0092) | 0.0505 (0.0576) | 29.54 | 6.10E-05 |
| 27 | rs284160 | TGFBR3 | 1 | A | ג' | 0.192644 | 0.1186 | 0.0572 (0.0096) | 0.1573 (0.0768) | 35.5 | 7.33E-05 |
| 28 | rs2842870 | PMF1, PMF1-BGLAP | 1 | C | T | 0.360398 | 0.3434 | -0.0446 (0.0087) | -0.0727 (0.0518) | 26.28 | 5.43E-05 |
| 29 | rs2880492 | NCOR2 | 12 | C | T | 0.035841 | 0.01984 | -0.1293 (0.0266) | -0.1917 (0.1781) | 23.63 | 4.88E-05 |
| 30 | rs35790371 | RBFOX1 | 16 | A | G | 0.005125 | 0.001903 | 0.5587 (0.1176) | 0.1447 (0.5420) | 22.57 | 4.66E-05 |
| 31 | rs5752720 | TTC28 | 22 | T | C | 0.28441 | 0.1844 | 0.0422 (0.0091) | -0.0311 (0.0635) | 21.51 | 4.44E-05 |
| 32 | rs57694670 | SH3PXD2A | 10 | ג' | A | 0.433638 | 0.3834 | -0.0507 (0.0086) | 0.0080 (0.0507) | 34.76 | 7.18E-05 |
| 33 | rs6462001 | - | 7 | T | ג' | 0.881494 | 0.8329 | -0.0765 (0.0151) | -0.0567 (0.0662) | 25.67 | 5.30E-05 |
| 34 | rs6843082 | - | 4 | A | ג' | 0.65294 | 0.6907 | -0.0438 (0.0092) | -0.0568 (0.0534) | 22.67 | 4.68E-05 |
| 35 | rs7091346 | SH3PXD2A | 10 | T | C | 0.45273 | 0.2664 | -0.0582 (0.0098) | -0.0143 (0.0552) | 35.27 | 7.28E-05 |
| 36 | rs7194129 | CFDP1 | 16 | T | C | 0.573691 | 0.5538 | 0.0399 (0.0084) | 0.0293 (0.0495) | 22.56 | 4.66E-05 |
| 37 | rs7341574 | ZFPM2 | 8 | ט | C | 0.313243 | 0.3895 | 0.0463 (0.0094) | -0.0774 (0.0505) | 24.26 | 5.01E-05 |
| 38 | rs7451833 | - | 6 | ג' | A | 0.103042 | 0.1166 | 0.1280 (0.0220) | 0.1253 (0.0777) | 33.85 | 6.99E-05 |
| 39 | rs74617384 | LPA | 6 | T | A | 0.070999 | 0.04573 | 0.1415 (0.0281) | -0.0532 (0.1162) | 25.36 | 5.24E-05 |
| 40 | rs74849463 | PIK3C2B | 1 | T | פחמן | 0.240949 | 0.2043 | 0.0445 (0.0092) | -0.0319 (0.0616) | 23.4 | 4.83E-05 |
| 41 | rs757241 | AFAP1-AS1 | 4 | סוף | G | 0.673377 | 0.6704 | -0.0681 (0.0145) | -0.0274 (0.0523) | 22.06 | 4.56E-05 |
| 42 | rs76099321 | CNNM2 | 10 | A | ג' | 0.052085 | 0.02979 | -0.0691 (0.0150) | -0.2215 (0.1464) | 21.22 | 4.38E-05 |
| 43 | rs7670136 | - | 4 | C | T | 0.571088 | 0.6269 | -0.0473 (0.0101) | 0.0391 (0.0509) | 21.93 | 4.53E-05 |
| 44 | rs77455924 | NTM | 11 | T | C | 0.056457 | 0.02996 | 0.1632 (0.0331) | 0.0444 (0.1477) | 24.31 | 5.02E-05 |
| 45 | rs7820334 | - | 8 | T | C | 0.246549 | 0.2949 | -0.0643 (0.0127) | -0.1327 (0.0648) | 25.63 | 5.29E-05 |
| 46 | rs7859727 | CDKN2B-AS1 | 9 | T | C | 0.518845 | 0.4151 | 0.0569 (0.0087) | 0.0281 (0.0498) | 42.77 | 8.83E-05 |
| 47 | rs7989823 | COL4A1, COL4A2 | 13 | C | A | 0.591615 | 0.6188 | 0.0527 (0.0089) | 0.0693 (0.0518) | 35.06 | 7.24E-05 |
| 48 | rs79960344 | - | 17 | ג' | T | 0.102808 | 0.1031 | 0.0651 (0.0130) | -0.1596 (0.0805) | 25.08 | 5.18E-05 |
| 49 | rs880315 | CASZ1 | 1 | C | T | 0.445525 | 0.4139 | 0.0416 (0.0089) | 0.0444 (0.0501) | 21.85 | 4.51E-05 |
| 50 | rs9112 | LOC100505841 | 5 | A | ג׳ | 0.436595 | 0.3573 | 0.0407 (0.0085) | -0.0274 (0.0523) | 22.93 | 4.74E-05 |
| 51 | rs979380 | - | 17 | A | ג' | 0.531844 | 0.6197 | -0.0417 (0.0084) | -0.0196 (0.0506) | 24.64 | 5.09E-05 |
| EAF.IS ו-EAF.VD מציינים את תדרי אלל-השפעה (effect-allele frequencies) במאגרי הנתונים של החשיפה לשבץ איסכמי ותוצאת הדמנציה הווסקולרית, בהתאמה. סטטיסטי F = (β/SE)². R² לכל SNP = F/(F + N − 2), כאשר N = 484,121. SNPs נבחרו באמצעות p < 5 × 10-6 וצמצום מקבצים (clumping) של חוסר שיווי משקל קישורי (linkage disequilibrium) באמצעות r² < 0.001 בתוך חלון של 10,000 kb. | |||||||||||
טבלה 1: כלים גנטיים לשבץ איסכמי והקשרים תואמים בין SNP לתוצאה עבור VaD. הטבלה מפרטת את 51 ה-SNPs הסופיים, כולל מזהה SNP, הגן הממופה, כרומוזום, אללים, EAF, הערכות SNP–trait, סטטיסטיקות F, ו-R2 לכל SNP. ה-SNPs נבחרו ב- p < 5 × 10-6 ובוצע עבורם clumping ב- r2 < 0.001 בטווח של 10,000 kb. EAF.IS ו-EAF.VD מציינים את השכיחויות של אללי האפקט במאגרי הנתונים של IS ו-VD, בהתאמה. Chr., כרומוזום; EA, אללי אפקט; OA, אליל אחר; EAF, שכיחות אללי אפקט; IS, שבץ איסכמי; VD, דמנציה וסקולרית. אנא לחצו כאן להורדת טבלה זו.
טבלה משלימה 1. הערכת MR-PRESSO של פלייוטרופיה אופקית עבור ניתוח השבץ האיסכמי – דמנציה וסקולרית.טבלה זו מסכמת את ניתוח ה-MR-PRESSO שבוצע כדי להעריך פלייוטרופיה אופקית גלובלית ווריאנטים פוטנציאליים שהם חריגים (outliers), תוך שימוש במערך הנתונים ההרמוני הסופי של 51 SNPs.אנא לחצו כאן להורדת קובץ זה.
טבלה נלווית 2. שלבי הפרוטוקול ותוצרי האימות התואמים. טבלה זו מקשרת בין כל שלב מרכזי בפרוטוקול לתוצר המייצג התואם לו ולמיקומו המדווח בכתב היד, ובכך מדגימה את היישום והאימות של זרימת העבודה האנליטית.אנא לחצו כאן להורדת קובץ זה.
דיון
מחקר זה מציג פרוטוקול MR של שני מדגמים להערכת הקשר הסיבתי הפוטנציאלי בין IS ל-VaD באמצעות סטטיסטיקות סיכום של GWAS הזמינות לציבור. בניתוח המייצג, שיטות ה-IVW והחציון המשוקלל תמכו בקשר חיובי בין IS שנחזה גנטית לבין סיכון ל-VaD, בעוד שהערכת ה-MR-Egger הייתה עקבית מבחינת הכיוון אך לא הגיעה למשמעות סטטיסטית. לפיכך, ממצאים אלו מספקים ראיות מרמזות לאפקט סיבתי אפשרי, ולא הוכחה חלוטה לסיבתיות.
הקשר שנצפה הוא סביר מבחינה ביולוגית בהקשר של פגיעה צרברו-וסקולרית. IS עלול להוביל לפגיעה נוירונלית אזורית, לליקוי בפרפוזיה המוחית ולשינויים נוירו-דגנרטיביים עוקבים. מחקרים קודמים הראו כי VaD מתרחש לעיתים קרובות במצב של זרימת דם מוחית מופחתת או נזק צרברו-וסקולרי13,14. שינויים צרברו-וסקולריים כרוניים בעקבות IS עשויים לתרום עוד יותר לירידה קוגניטיבית ולסיכון ל-VaD15. בנוסף, פגיעה מיקרווסקולרית לאחר שבץ, תגובות דלקתיות ודיספונקציה נוירו-וסקולרית כרונית עשויים לספק קשרים ביולוגיים פוטנציאליים בין IS ל-VaD. עם זאת, פרוטוקול MR הנוכחי אינו יכול לבסס מנגנונים אלו באופן ישיר, ויש לפרש את הממצאים יחד עם מחקרי תיקוף קליניים ומכניסטיים עתידיים.
הערכת MR-Egger לא הגיעה למשמעות סטטיסטית. ממצא זה צריך להיבחן בזהירות, ולא כסתירה ישירה לתוצאות של IVW והחציון המשוקלל (weighted median). שיטת MR-Egger יכולה לספק הערכות עמידות יותר לצורות מסוימות של פלייוטרופיה כיוונית (directional pleiotropy), אך בדרך כלל היא בעלת עוצמה סטטיסטית נמוכה יותר, במיוחד כאשר השפעות המכשירים (instrument effects) הן מתונות. בניתוח זה, נקודת החיתוך של MR-Egger לא הצביעה על פלייוטרופיה אופקית כיוונית בעלת משמעות סטטיסטית, ובדיקות ההומוגניות לא הראו הטרוגניות משמעותית בין ה-SNP. בנוסף, הבדיקה הגלובלית של MR-PRESSO לא הראתה עדות לפלייוטרופיה אופקית גלובלית. ניתוחי רגישות אלו מפחיתים את החשש מפני פלייוטרופיה אופקית כיוונית או גלובלית מדידה, אך הם אינם שוללים את כל מקורות ההטיה האפשריים או קובעים קשרי סיבתיות.
עם זאת, יש להביא בחשבון את מספרם הצנוע של מקרי ה-VaD בעת פירוש ממצאים אלו. למרות שניתוח ה-IVW הראשוני הראה עוצמה מוערכת בדיעבד של 89.5% בהתנייה לגודל האפקט שנצפה, תוצאה זו היא תיאורית ואין לפרש אותה כהצדקה מוגדרת מראש לגודל המדגם. מספר המקרים המוגבל הפחית במיוחד את הדיוק של אומדני רגישות פחות יעילים; האפקט המינימלי שניתן היה לזהות (minimum detectable effect) ב-80% עבור MR-Egger היה OR של 3.23, רמה שהייתה גבוהה מהאומדן שנצפה. בהתאם לכך, התוצאה הבלתי משמעותית של MR-Egger אינה מבטלת את חוסר הוודאות לגבי עוצמת הקשר. חישובי עוצמה אלו אינם מתייחסים גם להטיות פוטנציאליות הקשורות לתקינות המכשיר (instrument validity), חפיפה במדגמים, הטרוגניות של הפנוטיפ או פליאוטרופיה אופקית שיורית.
מספר שלבים אנליטיים הם קריטיים לאמינות של תהליך עבודה זה של MR. ראשית, יש לזהות בבירור את מערכי הנתונים של GWAS עבור החשיפה והתוצאה באמצעות מזהי מערכי נתונים או מידע על גישה (accession information). שנית, יש לבצע בחירת SNP וצבירה (clumping) של חוסר שיווי משקל בקישור (linkage disequilibrium) באמצעות ספים שנקבעו מראש כדי להשיג כלים גנטיים עצמאיים. שלישית, יש להעריך את חוזק הכלי באמצעות סטטיסטיקות F כדי להפחית את הסיכון להטיה של כלי חלש. רביעית, הרמוניזציה של אללים היא חיונית כדי להבטיח שהערכות של SNP–חשיפה ו-SNP–תוצאה תואמות לאותו אלל השפעה. לבסוף, יש להשתמש בניתוחים של הטרוגניות, פליאוטרופיה אופקית וניתוחי leave-one-out כדי להעריך האם ההערכה הראשית מושפעת מהשפעות SNP לא עקביות, פליאוטרופיה כיוונית או כלי גנטי דומיננטי יחיד.
מקורים נפוצים לכישלון או להטיה בזרימה זו כוללים מספר לא מספיק של SNPs הקשורים לחשיפה, כלים חלשים, חסרון של SNPs עבור התוצאה, יישור אללים מעורפל, וריאנטים פלינדרומיים, חוסר התאמה של SNP פרוקסי, הטרוגניות בין אומדנים ספציפיים ל-SNP, פליאוטרופיה אופקית, חוסר התאמה באוכלוסייה, הטרוגניות של הפנוטיפ וחפיפה אפשרית של דגימות. ניתן לטפל בבעיות אלו על ידי בדיקת מזהי מערכי הנתונים של ה-GWAS, החלה של ספי שיכרון וספי Clumping עקביים לבחירת SNP, תיעוד השימוש ב-SNP פרוקסי, בדיקת פלטי ההרמוניזציה, סיכום של F-statistics ופירוש אומדני MR יחד עם תוצאות של הטרוגניות, פליאוטרופיה ואנליזות leave-one-out. אם מזוהה הטרוגניות או פליאוטרופיה משמעותית, יש לפרש את אומדן ה-IVW הראשוני בזהירות ולשקול ניתוחי רגישות נוספים או מערכי נתונים חלופיים.
יש לשקול גם חפיפה פוטנציאלית של הדגימות. GWAS של החשיפה ו-GWAS של התוצאות מ-FinnGen נגזרו שניהם ממשאבים של מוצא אירופאי, ולא ניתן היה לכמת את מידת החפיפה בין המשתתפים מהסטטיסטיקות המסכמות הזמינות. חפיפה כזו עלולה להטות את הערכות ה-MR של שני דגימות לעבר קשרים תצפיתיים, במיוחד כאשר הכלים אינם קשורים חזק לחשיפה. מגבלה זו צריכה להילקח בחשבון בעת פרשנות התוצאה המייצגת.
בהשוואה למחקרי עקב (cohort) תצפיתיים קונבנציונליים או למחקרי מקרה-ביקורת, MR יכול להפחית הטיות הנובעות ממשתנים מתערבים ומסיבתיות הפוכה, תחת הנחות תקפות של משתנה מכשיר. עם זאת, מחקרים תצפיתיים נותרים בעלי ערך להערכת שכיחות מחלות, דפוסים זמניים ופרוגנוזה קלינית. בהשוואה לניסויים מבוקרים אקראים או למודלים ניסויים, MR יכול להעריך קשרי סיבתיות פוטנציאליים באמצעות נתונים גנטיים קיימים כאשר מחקרי התערבות ישירים אינם ישימים או אינם אתיים. עם זאת, MR אינו יכול לחשוף באופן ישיר מנגנונים תאיים או להחליף תיקוף מכניסטי. לפיכך, יש לראות פרוטוקול זה כמשלים למחקרי עקב, ניתוחים חתךים, מודלים של בעלי חיים וניסויי תאים, ולא כתחליף להם.
מחקרים קודמים בחנו גורמי סיכון קרייביליים ותוצאים הקשורים לדמנציה באמצעות גישות אפידמיולוגיות וגנטיות. לדוגמה, סוכרת ויתר לחץ דם נמצאו קשורים לסיכון לדמנציה במחקרי אוכלוסין ומחקרים קליניים קודמים16. מחקרים קיימים המבוססים על MR בחנו גם IS, VaD ומנגנוני תחלואה נלווים קשוריים. כתב יד זה נבדל ממחקרים אלו במטרתו העיקרית. במקום לשאוף לזיהוי מטרות מולקולריות חדשות או לאישור ניסויי של מנגנוני מחלה, מאמר זה מתמקד בהצגת פרוטוקול הניתן לשחזור לביצוע ניתוח MR של שני מדגמים באמצעות נתוני סיכום של GWAS הזמינים לציבור.
לפרוטוקול זה מספר מגבלות. ראשית, הניתוח התבסס על סטטיסטיקות סיכום של GWAS הזמינות לציבור, והתוצאות עשויות שלא להיות ניתנות להכללה מלאה לאוכלוסיות שאינן אירופאיות או למסדי נתונים עם הגדרות מקרה שונות. שנית, נעשה שימוש בסף בחירת SNP מקל (p < 5 × 10-6) כדי לשמר מספר מספיק של מכשירים עצמאיים עבור ניתוחי רגישות. למרות שכל המכשירים שנשמרו הראו F-statistics >10, סף זה פחות מחמיר מסף המובהקות המקובל ברמת הגנום, ועלול לכלול וריאנטים עם קשרי חשיפה חלשים או פחות אמינים. לפיכך, יש לפרש את הממצאים המייצגים כחקרים ראשוניים ומסמנים ולא כראיה חלוטה לקשר סיבתי. שלישית, למרות ש-MR יכול להפחית הטיות של משתנים מתערבים וסיבתיות הפוכה, הוא אינו יכול לשלול לחלוטין פליאוטרופיה אופקית, סיווג שגוי של פנוטיפים, חפיפת דגימות, או הטיות שנובעות ממכשירים לא תקפים. רביעית, הערכת MR-Egger שאינה מובהקת מעידה על כך שיש לנקוט משנה זהירות בפרשנות הסיבתית. חמישית, פרוטוקול זה אינו כולל ניסויי תאים, מודלים של בעלי חיים או תיקוף של קבוצות עצמאיות; גישות אלו יהיו בעלות ערך להבהרת מנגנונים ביולוגים וישימות קלינית.
בנוסף, IS מודל כפנוטיפ כולל. פרוטוקול מייצג זה לא תוכנן כדי לזהות אם קשר כלשהו מיוחס לשבץ מוחי התרתית בעורקים גדולים, לחסימה של כלי דם קטנים, לשבץ קרדיו-אמבולי או לתת-סוג אחר. כתוצאה מכך, אין לפרש את התוצאה כראיה לאפקט הספציפי לתת-סוג. עבודות עתידיות צריכות להחיל את אותו זרימת עבודה על מאגרי נתוני GWAS בעלי עוצמה סטטיסטית מספקת עבור תת-סוגים ספציפיים, תוך שימוש בתוכנית ניתוח מותאמת ומוגדרת מראש.
יישומים עתידיים של פרוטוקול זה עשויים להתרחב מעבר ל-IS ו-VaD. ניתן להתאים את אותו זרימת עבודה להערכה של קשרים אחרים בין חשיפה לתוצאה במחלות צרברו-ווסקולריות, ניוורו-דגנרטיביות, מטבוליות, דלקתיות או קרדיו-ווסקולריות, כאשר זמינים נתונים מסכמים מתאימים של GWAS. מחקרים עתידיים עשויים להרחיב זרימת עבודה זו על ידי שילוב של תיקוף GWAS עצמאי, MR דו-כיווני, MR רב-משתני, MR תיווך, או אינטגרציה עם ראיות ניסיוניות ומבוססות קוהורט.
לסיכום, פרוטוקול זה מספק זרימת עבודה ניתנת לשחזור לביצוע ניתוח MR של שני מדגמים באמצעות סטטיסטיקות סיכום של GWAS ממאגרים ציבוריים. התוצאות המייצגות מספקות ראיות עקיפות לקשר חיובי אפשרי בין IS שנחזה גנטית לבין הסיכון ל-VaD, כאשר התמיכה העיקרית הגיעה משיטות ה-IVW והחציון המשוקלל, בעוד שהערכת MR-Egger לא הגיעה למשמעות סטטיסטית. עם זאת, אין לפרש את הממצאים כהוכחה מכניסטית מוחלטת. נדרשים מחקרי תיקוף עצמאיים מבוססי אוכלוסייה וחקירות ניסיוניות כדי להבהיר יותר את המנגנונים הביולוגיים ואת הרלוונטיות הקלינית של קשר זה.
גילויים
המחברים מצהירים כי אין להם ניגודי אינטרסים כספיים. לא נעשה שימוש בכלי בינה מלאכותית (AI) או במודלי שפה גדולים (LLM) בכתיבה, בעריכה או בהכנה של כתב יד זה.
תודות
המחברים לא קיבלו מימון ספציפי למחקר זה.
חומרים
| שם | חברה | מספר קטלוג | הערות |
|---|---|---|---|
| סט נתוני FinnGen GWAS | FinnGen | finn-b-F5_VASCDEM | סט נתוני GWAS ציבורי ששימש כמקור לתוצאה עבור דמנציה וסקולרית |
| מאגר IEU OpenGWAS | MRC Integrative Epidemiology Unit, University of Bristol | ebi-a-GCST90018864 | סט נתוני GWAS ציבורי ששימש כמקור לחשיפה עבור שבץ איסכמי |
| תוכנת R | R Foundation for Statistical Computing | Not Applicable | תוכנה סטטיסטית ששימשה לניתוח MR |
| RStudio | Posit Software, PBC | Not Applicable | סביבת פיתוח משולבת ששימשה להרצת סקריפטים של R |
| חבילת TwoSampleMR | MRC Integrative Epidemiology Unit, University of Bristol | Not Applicable | חבילת R ששימשה לחילוץ מכשירים, הרמוניזציה, ניתוח MR וניתוחי רגישות |
| חבילת MRPRESSO | CRAN | Version 1.0 | חבילת R ששימשה להערכת פליאוטרופיה אופקית גלובלית ומכשירים חריגים פוטנציאליים |
| חבילת writexl | CRAN | Not Applicable | חבילת R ששימשה לייצוא טבלאות תוצאות כקבצי xlsx |
מקורות
- Sharma R, Lee K. Advances in treatments for acute ischemic stroke. BMJ. 2025;389:e076161.
- Zhu H, et al. Interleukins and ischemic stroke. Front Immunol. 2022;13:828447.
- Feske SK. Ischemic stroke. Am J Med. 2021;134(12):1457-1464.
- Lin HF, et al. Apolipoprotein E polymorphism in ischemic cerebrovascular diseases and vascular dementia patients in Taiwan. Neuroepidemiology. 2004;23(3):129-134.
- Semplicini A, et al. Hypertension and cerebrovascular diseases: a specific role of vascular protection for the prevention of dementia. J Cardiovasc Pharmacol. 2001;38(Suppl 2):S79-S82.
- Bulwa Z, et al. Management of blood pressure after acute ischemic stroke. Curr Neurol Neurosci Rep. 2019;19(6):29.
- Hosoki S, Tanaka T, Ihara M. Diagnostic and prognostic blood biomarkers in vascular dementia: from the viewpoint of ischemic stroke. Neurochem Int. 2021;146:105015.
- Yamagata K. Docosahexaenoic acid inhibits ischemic stroke to reduce vascular dementia and Alzheimer's disease. Prostaglandins Other Lipid Mediat. 2023;167:106733.
- Yeung S, Luo S, Iwagami M, Goto A. Introduction to Mendelian randomization. Ann Clin Epidemiol. 2025;7(1):27-37.
- Luo J, et al. Cardiovascular diseases and risk of dementia in the general population. Eur J Prev Cardiol. 2025:zwaf129. Available from: https://academic.oup.com/eurjpc
- Burgess S. Sample size and power calculations in Mendelian randomization with a single instrumental variable and a binary outcome. Int J Epidemiol. 2014;43(3):922-929. Available from: https://academic.oup.com/ije/article/43/3/922/757361
- Verbanck M, et al. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693-698. Available from: https://www.nature.com/articles/s41588-018-0099-7
- Ng S, et al. Updates on vascular dementia. Stroke Vasc Neurol. 2025;10(5):542-550.
- Sanders AE, Schoo C, Kalish VB. Vascular Dementia. StatPearls Publishing; Treasure Island (FL); 2025. Available from: https://www.ncbi.nlm.nih.gov/books/
- Ho JP, Powers WJ. Contemporary management of acute ischemic stroke. Annu Rev Med. 2025;76(1):417-429.
- Yen FS, et al. Diabetes, hypertension, and the risk of dementia. J Alzheimers Dis. 2022;89(1):323-333.