$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
حللت هذه الدراسة إحصائيات دراسة الارتباط الجينومي الواسع (GWAS) غير المحددة على المستوى الملخص والمتاحة للجمهور. وفقا لسياسات المستودع والموافقات التي حصل عليها الباحثون الأصليون، لم تكن هناك حاجة لموافقة لجنة مراجعة مؤسسية جديدة أو موافقة فردية إضافية مستنيرة لهذا التحليل الثانوي. جميع المشاركين في GWAS أبلغوا عن إجراءات الموافقة والأخلاقيات والموافقة في منشوراتهم المصدرية. أجريت جميع التحليلات وفقا للإرشادات المؤسسية وإعلان هلسنكي.
نظرة عامة والمنطق
نفذت الدراسة إطار عمل عشوائي مندلي ثنائي الاتجاه (MR) مكون من عينتين مقتصرا على إحصاءات ملخص النسب الأوروبية لتقييم العلاقات السببية المحتملة بين التصلب المتعدد (MS) والأورام الدموية الخبيثة (HM). يلتزم التصميم بثلاثة افتراضات أساسية للبحث المغناطيسي: ملاءمة الأداة، الاستقلالية عن المربكات، وتقييد الاستبعاد. لذلك، يشمل سير العمل (1) الوصول إلى مجموعات البيانات وتنظيمها، (2) اختيار الأدوات عند دلالة جينومية شاملة مع تجمع عدم توازن الربط (LD)، (3) فحص الارتباك باستخدام PhenoScanner، (4) تنسيق الأليل مع التعامل الصريح مع متغيرات الباليندروم، (5) تقييم الاتجاه باستخدام اختبار Steiger¹²، (6) تقدير MR الأولي مع طرق تكميلية، (7) مجموعة كاملة من تشخيصات الحساسية، و(8ii) توليد الأشكال والجداول الموحدة تحت اختبار متعدد للتحكيم. يتم وصف كل خطوة من هذه الخطوات بالتفصيل في الأقسام الفرعية التالية للبروتوكول، ويتم عرض نظرة عامة على خط الأنابيب في الشكل 1.
المواد والبرمجيات وأجهزة RRID
أجريت التحليلات في النسخة R 4.3.1 (RRID:SCR_001905) باستخدام RStudio/Posit 2023.12+ (RRID:SCR_000432). عند تنفيذ التجميع محليا، استخدم PLINK v1.9 (النسخة 2.3؛ RRID:SCR_001757)13. تم استخدام تقدير MR واستخراج البيانات حزمة R TwoSampleMR v0.5.7 10؛ البحث عن أجهزة للمربكين المحتملين باستخدام Phenoscanner v1.0؛ استخدمت عملية اكتشاف وتصحيح الشاذة باستخدام MRPRESSO v1.0. يتم الإبلاغ عن النسخ الدقيقة للحزم التي لا تحتوي على RRIDs.
مصادر البيانات والوصول إليها
تم الحصول على إحصائيات ملخصة للتصلب المتعدد من تحليل التلوي لاتحاد الوراثة المتعدد الدولي الذي يضم 47,429 حالة تصلب متعدد و68,374 حالة ضابطة مع ضبط جودة موحد عبر 15 مجموعة. تم الحصول على إحصائيات ملخصة HM من FinnGen (إجمالي n = 218,792؛ >16 مليون متغير) وشملت لمفوما هودجكين (HL)، لمفوما الخلايا B الكبيرة المنتشرة (DLBCL)، لمفوما الجريبات (FL)، لمفوما الخلايا التائية/NK الناضجة (MTNKL)، لمفوما غير هودجكين أخرى أو غير محددة (NHL)، لوكيميا لمفاوي، لوكيميا نخاعي، لوكيميا من نوع خلايا غير محدد، وأورام متعددة من الخلايا البلازمية/14. تم الوصول إلى مجموعات البيانات من خلال بوابة IEU OpenGWAS باستخدام معرفات الوصول الموثقة15. لذلك، استندت جميع التحليلات في هذه الدراسة حصريا إلى مجموعات بيانات GWAS على مستوى ملخص متاحة للجمهور؛ لم تستخدم أو يتم إنشاء بيانات من المجموعات المؤسسية الداخلية أو بيانات المرضى على مستوى الفرد. وبما أننا لم نحدد GWAS إضافية مع تعريفات فرعية متناسقة للتصلب المتعدد والأورام الدموية تسمح بتكرار كامل لخط الأنابيب، لم يتم إجراء التحقق الخارجي المستقل باستخدام مجموعة بيانات منفصلة وتم الاعتراف به كقيد. البروتوكول مكتوب بحيث يمكن إعادة تطبيقه مباشرة على مجموعات بيانات GWAS المستقبلية للتحقق المستقل.
اختيار الآلات وتجمع LD
لكل تعرض، تم اختيار تعدد أشكال النيوكليوتيد الواحد (SNPs) عند دلالة جينومية واسعة (P < 5 × 10-8) باستخدام دالة extract_instruments في TwoSampleMR المطبقة على مجموعات بيانات OpenGWAS. لضمان استقلالية الأجهزة، تم إجراء تجميع LD ضد لوحة مرجعية للأصول الأوروبية باستخدام أدوات التجمع الداخلية ل TwoSampleMR أو محليا باستخدام PLINK، مع عتبة r² 0.001 ونافذة فيزيائية تبلغ 10,000 كيلوبايس. عند استخدام PLINK، تم تعيين معلمات سطر الأوامر على عتبة دلالة أولية 5 × 10-8، r² = 0.001، ونافذة 10 ميجابايت بحيث تتطابق الأجهزة المجمعة تماما مع هذه المعايير. تم تقييم قوة الجهاز باستخدام إحصائية F المستمدة من تقدير تأثير التعرض وخطأها القياسي (F ≈ β²/SE²); تم استبعاد المتغيرات التي تحتوي على F < 10 من مجموعات الأجهزة النهائية، وتم نقل باقي SNPs إلى فحص PhenoScanner.
فحص التشويش باستخدام PhenoScanner
لتقليل التنوع المتعدد الأفقي عبر عوامل الخطر المعروفة، تم استعلام كل أداة مرشحة في PhenoScanner V2 عبر كتالوج GWAS باستخدام حزمة Phenoscanner R (v1.0)16,17. لكل SNP، طلبنا جميع الارتباطات المبلغ عنها عند P < 1 × 10⁻5 وفحصنا الصفات العائدة يدويا. الارتباطات التي تشير إلى وجود روابط مع عوامل خطر الأورام الدموية الخبيثة المعروفة - مثل التعرض المرتبط بالتدخين أو صفات السمنة/الأنثروبومترية (مثل مؤشر كتلة الجسم، محيط الخصر، ومقاييس الدهون في الجسم) - أو الارتباطات المباشرة مع الأنماط الظاهرية للأورم الدموي الخبيثة دفعت إلى استبعاد SNP المقابل من مجموعة الأدوات18. كانت فئات السمات التي اعتبرت استبعادية مبنية على أدلة سابقة تربط السمنة والتدخين بخطر اللوكيميا أو اللمفوما أو الورم النخاعي 18,19,20. استخدمت الاستعلامات جذور كلمات مفتاحية واسعة (مثل الدخان، السجائر، مؤشر كتلة الجسم، السمنة، الخصر، السمنة، السمنة، الأورام الدموية الخبيثة، اللمفوما، اللوكيميا، الورم النخاعي). تم توثيق جميع عمليات الإزالة في جدول تتبع مع صفة PhenoScanner التي أدت إلى الاستبعاد، ثم تم تمرير قوائم الأدوات التي تم تنظيفها إلى خطوة التناغم.
التوافقية والتعامل مع الباليندرومي
تم تنسيق أليلات التأثير لكل SNP بين بيانات التعرض والنتائج باستخدام دالة harmonise_data في حزمة TwoSampleMR (v0.5.7، R). قمنا بمحاذاة جميع أليلات النتائج مع أليل تأثير التعرض بحيث كانت معاملات بيتا الإيجابية دائما تتوافق مع نفس الأليل في كلا مجموعتي البيانات. تم التعامل مع المتغيرات الباليندرومية (A/T أو C/G) ذات ترددات تأثير أليل متوسطة (0.42-0.58) في لوحة مرجعية OpenGWAS على أنها غامضة عن السلسلة وتمت إزالتها تلقائيا عن طريق ضبط إجراء التوافق لإسقاط SNPs غامضة. تم الاحتفاظ ومحاذاة SNPs باليندرومية ذات ترددات الأليل التأثير خارج هذا النطاق باستخدام ترددات الأليل المبلغ عنها. نظرا لاختلاف توفر الأليلات والحالة الباليندرومية قليلا بين نتائج فينجين، تم إجراء التوحيد بشكل منفصل لكل نمط ظاهري من HM، وتم استخراج العدد النهائي من الأدوات التي تدخل كل تحليل خاص بالنتيجة من كائنات R المنسقة والإبلاغ عنها في الجداول.
تقييم الاتجاهية (ترشيح ستايجر)
تم تقييم الاتجاهية باستخدام نهج ستايجر كما هو مطبق في وظيفة steiger_filtering ل TwoSampleMR. لكل SNP، كانت الدالة تحسب أولا التباين المفسر (R²) في التعرض والنتيجة من معامل بيتا GWAS، والخطأ المعياري، وحجم العينة. ثم أزالت الدراسة الأدوات التي كان R² أكبر في النتيجة مقارنة بالتعرض، مما يشير إلى احتمال عكس اتجاه التأثير. تم تطبيق ترشيح ستايغر بشكل منفصل لكل مجموعة بيانات نتائج، وتم حفظ الأدوات المتبقية (صفوف تحتوي على steiger_dir == صحيح) واستخدامها في تحليلات MR اللاحقة. تم تسجيل أعداد الأدوات بعد اختبار ستايجر لكل نتيجة ويتم الإبلاغ عنها جنبا إلى جنب مع تقديرات التصوير بالرنين المغناطيسي.
تقدير الرنين المغناطيسي الأولي والتحكم في الاختبارات المتعددة
تم الحصول على تقديرات سببية أولية باستخدام MR المرجح بالتباين العكسي (IVW) ضمن نموذج التأثيرات الثابتة باستخدام دالة MR في TwoSampleMR، مع تحديد الطرق ك "mr_ivw"، "mr_egger_regression"، و"mr_weighted_median". لكل نتيجة من نتائج HM، تم تمرير الأدوات المتناغمة والمفلترة بستايجر إلى MR، وتم استخراج نسب الاحتمالات اللوغاريتمية والأخطاء القياسية وزيادتها للحصول على نسب احتمالات (ORs) مع فترات ثقة 95٪ (CIs) للصفاتالثنائية 21. لفحص المتانة للانتهاكات الطفيفة لافتراض عدم تعدد التروبي، طبقنا أيضا مقدرات الانحدار الوسيط المرجح وMR-Egger 22,23، المطبقة في نفس الحزمة. عندما أشار اختبار Q الخاص بكوكران (من mr_heterogeneity) إلى وجود تباين كبير (P < 0.05)، قامت الدراسة أيضا بمطابقة نماذج IVW العشوائية المتعددة وقدمت نتائج ثابتة وعشوائية. تم التحكم في الخطأ حسب العائلة عبر تسع نتائج HM باستخدام تصحيح بونفيروني بمعدل α = 0.05/9 = 5.56 × 10-3؛ تم اعتبار الارتباطات التي تحتوي على قيم P أقل من هذا الحد ذات دلالة إحصائية، بينما تم تفسير الارتباطات التي لديها 0.0056 ≤ P < 0.05 على أنها موحية وتم وصفها بحذر.
تشخيصات الحساسية: التباين، التعدد العددي، والقيم الشاذة
تم استخدام إحصائية Q الخاصة بكوكران لتقييم التباين بين الأدوات لكل من نماذج IVW وMR-Egger، وتم تنفيذها عبر دالة mr_heterogeneity في TwoSampleMR. تم تقييم التسلسل الأفقي الاتجاهي باستخدام اختبار اعتراض MR-Egger (mr_pleiotropy_test) والاختبار العالمي في حزمة MR-PRESSOرقم 24. تم تشغيل MR-PRESSO24 مع الإعدادات الموصى بها في R (NbDistribution ≥ 5,000، SignifThreshold = 0.05) لاكتشاف المؤشرات الشاذة المؤثرة وقياس التشوهات المحتملة من خلال مقارنة تقديرات IVW قبل وبعد إزالة القيمالشاذة 25. تم إجراء تحليلات "ترك واحد خارج" (mr_leaveoneout) لكل زوج تعرض-نتيجة لتحديد ما إذا كان هناك SNP واحد يؤثر بشكل غير متناسب على التقدير العام. لأجل الشفافية وقابلية التكرار، تم تصدير جميع المخرجات التشخيصية من R وتم الإبلاغ عنها مع أعداد الأجهزة المقابلة بعد التوحيد والترشيح Steiger وإزالة الشاذ MR-PRESSO.
قوة الأجهزة وتقييم NOME
تم قياس قوة الجهاز ل MR-Egger باستخدام إحصائية I2GX، التي حسبت ك 1 ناقص متوسط مربعات الأخطاء القياسية لارتباطات التعرض مع SNP مقسومة على تباينها عبر الأدوات26. تشير القيم الأقرب إلى 1 إلى امتثال أفضل لافتراض عدم وجود خطأ في القياس (NOME)؛ تشير القيم المنخفضة إلى احتمال تخفيف الانحدار وتفسير حذر سريع لنتائج MR-Egger. تم حساب وإبلاغ I2GX لكل تحليل خاص بالنتائج.
العشوائية المندلية العكسية
تم تكرار خط الأنابيب الكامل في الاتجاه العكسي من خلال التعامل مع كل نوع فرعي من HM كتعرض للتعرض وMS كنتيجة. عندما كانت الأدوات ذات الدلالة الجينومية على مستوى الجينوم غير كافية لتعرض معين لمرض HM، تم السماح بعتبة اختيار مخففة من P < 5 × 10-6 مع الحفاظ على نفس معايير تجمع LD، وفحص PhenoScanner، وإجراءات التوافق، وترشيح Steiger، وتشخيص الحساسية. التحليلات التي استخدمت عتبات مرنة كانت موسمت بوضوح في الجداول والأساطير المشابهة للصور.
التصوير البصري وتصدير الأشكال
تم إنشاء مخططات التشتت، والغابة، والقمع، وقطع الرسم مع وضع الأساطير أسفل اللوحات وتعديل أحجام الخطوط لضمان عدم إخفاء البيانات المرسومة على الملصقات. تم توحيد حدود المحور عبر نتائج مماثلة لتسهيل المقارنة البصرية. تم تصدير المجسمات بدقة لا تقل عن 300 نقطة في البوصة في البوصة بصيغ بدون فقدان مثل TIFF أو PNG. تم التحقق من جميع القيم الرقمية المرسومة مع التقديرات المبلغ عنها لضمان الاتساق بين النص والجداول والأشكال.
قابلية التكرار ومشاركة البيانات
تم تثبيت البذور العشوائية حيثما كان ذلك مناسبا، وتسجيل نسخ البرمجيات، وأرشفت سكريبتات التحليل مع الكائنات الوسيطة للسماح بإعادة تشغيل جميع الخطوات. تم توثيق معرفات الوصول إلى مجموعات البيانات وتعريفات الأنماط الظاهرية، وتم إعداد قوائم الأدوات في كل مرحلة تصفية - بعد التجمع، بعد التوحيد (ما بعد التجميع)، التصفية بعد ستايجر، و"ما بعد MR-PRESSO" للرفع كملفات بيانات وفقا لإرشادات المجلة.