دراسة العلاقة ثنائية الاتجاه بين آلام أسفل الظهر، وألم الورك، وألم الركبة، واضطرابات المشي: دراسة عشوائية مندلية

ملخص

تشير تحليلات IVW إلى أن آلام أسفل الظهر والورك قد تكون مرتبطة باضطرابات في المشي، في حين أن ألم الركبة والارتباطات العكسية لا تزال غير حاسمة.

الملخص

ترتبط اضطرابات المشي بتغيرات في ميكانيكا الأطراف السفلية وضعف وظيفي في اضطرابات الجهاز العضلي الهيكلي. على الرغم من أن آلام أسفل الظهر وألم الورك وألم الركبة قد تم ربطها بضعف في المشي في الدراسات الرصدية، إلا أن علاقاتها السببية لا تزال غير واضحة. أجرينا دراسة عشوائية مندلية ثنائية الاتجاه (MR) لتقييم الارتباطات السببية المحتملة بين الألم الخاص بالموقع واضطرابات المشية. تم استخدام تعدد أشكال النيوكليوتيد المستقلة من دراسات الارتباط الجينومي واسعة النطاق كأدوات جينية لآلام أسفل الظهر، وألم الورك، وآلام الركبة. تم استخدام تشوهات المشي، التي عرفت في GWAS المصدر بأنها صعوبة مشي تم الإبلاغ عنها ذاتيا، كنتيجة في تحليلات الرنين المغناطيسي الأمامي وكتعرض في تحليلات الرنين المغناطيسي العكسي. كانت طريقة التحليل الوزني بالتباين العكسي (IVW) هي التحليل الأساسي، مدعومة بتحليلات الحساسية للتنوع العددي الأفقي والثبات. في التحليلات المستقبلية، ارتبطت ألم أسفل الظهر المتوقع وراثيا (OR = 1.53، فترة الثقة 95٪ = 1.058–2.219، P = 0.024) وألم الورك (OR = 1.55، CI 95٪ = 1.067–2.252، P = 0.021) بزيادة خطر حدوث تشوهات في المشي في تحليلات IVW. ومع ذلك، لم تكن الطرق الأربع الإضافية للرنين المغناطيسي ذات دلالة إحصائية، رغم أن التقديرات كانت متسقة عموما من حيث الاتجاه. أظهر ألم الركبة المتوقع وراثيا عدم وجود ارتباط ذي دلالة (OR = 2.15، فترة الثقة 95٪ = 0.234–19.789، P = 0.498)، مع انعكاس عدم يقين كبير في فترة الثقة الواسعة. لم تجد تحليلات الرنين المغناطيسي العكسي أي دليل على أن المسؤولية الوراثية عن تشوهات المشي تزيد من خطر آلام أسفل الظهر أو الورك أو الركبة، رغم توفر أربعة تشوهات فقط في الشذوذ في المشية. لم يتم اكتشاف أي تباين جوهري أو تعدد الأبعاد الأفقي. بشكل عام، تقدم هذه الدراسة أدلة جينية محدودة قائمة على IVW على أن آلام أسفل الظهر وألم الورك قد يرتبطان باضطرابات في المشي. يجب تفسير النتائج بحذر والتحقق منها باستخدام مجموعات بيانات GWAS أكبر مع أنماط ظاهرية محسنة.

المقدمة

يمكن أن تنبع اضطرابات المشي من إعاقات حركية أو حسية، وخصائصها السريرية تعتمد على موقع وطبيعة المرض الأساسي. قد توفر أنماط المشي غير الطبيعية المحددة أدلة تشخيصية مهمة لأمراض معينة. ومن الجدير بالذكر أنه في اضطرابات الجهاز العضلي الهيكلي، تظهر اضطرابات المشي غالبا بخصائص بيوميكانيكية مميزة وترتبط ارتباطا وثيقا بالألم، وضعف المفاصل، وانخفاض الحركة ^1,2,3. أظهرت الدراسات أن اضطرابات المشي طويلة الأمد يمكن أن تزيد من سوء حالة مرضى التهاب المفاصل العظمي⁴. علاوة على ذلك، قد تزيد المشية غير الطبيعية المستمرة من عدم استقرار المفاصل الحركية وتغير تحميل المفاصل؛ قد تؤدي الأنشطة المتكررة للمفاصل تحت الميكانيكا غير الطبيعية إلى تفاقم إصابة المفاصل، مما يسبب الألم، ويعزز ضمور العضلات تدريجيا، وينتج عنه فقدان أقوى عضلي بشكل واضح. قد تتفاعل هذه العوامل لتشكيل حلقة مفرغة، مما يساهم في تفاقم الألم وضعف الوظائف لدى المرضى الذين يعانون من اضطرابات عضلية هيكلية⁵. لذلك، قد يساعد توضيح اتجاه وسببية الروابط بين الألم العضلي الهيكلي الخاص بالموقع واضطرابات المشي في تفسير العلاقات الرصدية بين الألم والمشي وإعلام الأبحاث المستقبلية حول خلل وإعادة تأهيل المشية.

ومع ذلك، وعلى الرغم من تراكم الأدلة الرصدية التي تربط الألم المخصص للموقع—مثل آلام أسفل الظهر، وألم الورك، وألم الركبة—بضعف أداء المشي، إلا أن الدراسات الرصدية التقليدية معرضة للارتباك، وخطأ القياس، والسببية العكسية، مما يصعب تحديد اتجاه وسببية هذه الارتباطات. بالإضافة إلى ذلك، قد تشمل العلاقة بين الألم والمشي تكيف الجهاز العصبي المركزي وتعديل المعلومات الحسية للألم المستلمة من المحيط، مما يعقد تفسير الروابط الملحوظة بين الألم والمشي. العشوائية المندلية (MR) هي نهج وبائي مصمم لتجاوز القيود الرئيسية للدراسات الرصدية وقد تم تطبيقه على نطاق واسع في أبحاث الاستدلال السببي. يستخدم MR تعدد أشكال النيوكليوتيدات الأحادية المستقلة (SNPs) كمتغيرات أدوات لاستنتاج علاقات سببية محتملة بين التعرضات والنتائج⁶. على سبيل المثال، بدلا من تعيين الأفراد لآلام أسفل الظهر، يمكن للطب المغناطيسي فحص ما إذا كانت المسؤولية الجينية لآلام أسفل الظهر مرتبطة بنتائج متعلقة بالمشي في مجموعات بيانات دراسات الارتباط الجينومي الواسعة (GWAS). نظرا لأن المتغيرات الجينية تخصص عشوائيا عند الحمل، أي عند تكوين الجنين، يمكن للرنين المغناطيسي تقليل التحيز الناتج عن عوامل مربكة مثل العمر، ومؤشر كتلة الجسم، والنشاط البدني، والحالات العضلية الهيكلية المصاحبة. هذا النهج مناسب لاستكشاف العلاقات السببية على مستوى السكان عندما تتوفر بيانات GWAS قوية وأدوات جينية صالحة. ومع ذلك، كانت الأنماط الظاهرية المستخدمة في هذه الدراسة سمات GWAS واسعة تم الإبلاغ عنها ذاتيا وقد لا تعكس بشكل كامل شدة الألم، مدة الأعراض، التشخيصات المؤكدة من الأطباء أو معايير المشي الموضوعية. من خلال الاستفادة من التوزيع العشوائي للمتغيرات الجينية عند الحمل، يمكن للرنين المغناطيسي أن يقلل بفعالية من التحيز الناتج عن المراوغات ويعكس السببية، وبالتالي يقدم أدلة سببية أكثر منطقية.

كان هدف هذه الدراسة إجراء تحليل بالرنين المغناطيسي ثنائي الاتجاه باستخدام بيانات ملخصة لآلام أسفل الظهر، وألم الورك، وآلام الركبة التي تم الإبلاغ عنها ذاتيا، إلى جانب بيانات GWAS عن اضطرابات المشي من قاعدة بيانات IEU Open GWAS. في هذه الدراسة، مثلت الأنماط الظاهرية للألم سمات ألم واسعة تم الإبلاغ عنها ذاتيا حسب الموقع، وكانت الاضطرابات في المشي تمثل صعوبة في المشي المبلغ عنها ذاتيا بدلا من تشخيصات مؤكدة من الأطباء أو معايير المشي المقاسة في المختبرات أو درجات وظيفية موحدة. هدف هذا التصميم ثنائي الاتجاه إلى تقييم اتجاه الارتباطات السببية المحتملة بين الألم العضلي الهيكلي الخاص بالموقع واضطرابات المشي باستخدام إحصائيات ملخص GWAS⁷.

البروتوكول

استخدمت هذه الدراسة بيانات ملخصة على مستوى دراسة الارتباط الجينومي (GWAS) المتاحة للجمهور من قاعدة بيانات IEU OpenGWAS وفهرس GWAS. لم يتم الوصول إلى بيانات المشاركين على مستوى الأفراد، ولم يتم تجنيد أي مشاركين بشريين جدد. تم الحصول على الموافقة الأخلاقية والموافقة المكتوبة المستنيرة في الدراسات الأصلية المساهمة في GWAS. لذلك، لم تكن هناك حاجة لموافقة إضافية من مجلس مراجعة المؤسسات للتحليل الثانوي الحالي للإحصاءات الملخصة المتاحة للجمهور.

مبدأ الدراسة والتصميم

تشمل الخطوات الأساسية لدراسة MR الحصول على بيانات ملخصة GWAS، فحص وتقييم SNPs، إجراء تحليل إحصائي، وتنفيذ تدابير مراقبة الجودة. كما هو موضح في الشكل 1، تعتمد دقة تحليل الرنين المغناطيسي على تحقيق ثلاثة افتراضات حاسمة: (1) افتراض الصلة: يجب أن تكون الوريدات ذات صلة بأنماط التعرض، وهي آلام أسفل الظهر، ألم الورك، وألم الركبة⁸؛ (2) افتراض الاستقلالية: IVs غير ذات صلة بالعوامل المربكة التي تؤثر على "نتيجة التعرض"؛ (3) افتراض الحصرية: تؤثر الوريدات على النتيجة فقط من خلال التعرض وليس من خلال مسارات أخرى⁹. سيتم تبادل عوامل التعرض ومتغيرات النتائج في كل تحليل لتحديد ما إذا كان هناك سبب عكسي بين الاثنين.

لتحسين قابلية التكرار، يتم توفير سير العمل الحسابي الكامل المستخدم لاستيراد إحصائيات GWAS المحلية الملخصة، واستخراج SNP، وتكتل LD، واستخراج SNP النتائج، وتناغم الأليلات، وتحليل الرنين المغناطيسي، وتحليل الحساسية، ورسم التشخيص، وتوفير المخرجات كرمز تكميلي 1. يصف النص الرئيسي خطوات البروتوكول الرئيسية، بينما يوفر الملف التكميلي 1 تنفيذ R على مستوى الأوامر والوظيفة المقابل.

مصادر البيانات

يتم الحصول على البيانات المتعلقة بآلام أسفل الظهر، وألم الورك، وآلام الركبة من قاعدة بيانات IEU openGWAS، والموقع الإلكتروني هو: https://gwas.mrcieu.ac.uk/. يمكن العثور على أسئلة محددة تستخدم لتعريف هذه الحالات في الملف التكميلي 2. تم تعريف آلام أسفل الظهر، وألم الورك، وألم الركبة بناء على عناصر الألم المحددة التي تم الإبلاغ عنها في مجموعات بيانات GWAS الأصلية. رقم معرف مجموعة بيانات آلام أسفل الظهر GWAS هو: ebi-a-GCST90018797، بحجم عينة يبلغ 468269، تشمل 22,413 مريضا يعانون من آلام أسفل الظهر، و445,856 ضابطا في المجموعة، مع إجمالي عدد SNP 24174741. رقم معرف مجموعة بيانات آلام الورك في GWAS هو: ebi-a-GCST90013968، بحجم عينة 407746 وإجمالي عدد SNP يبلغ 11039206. رقم معرف مجموعة بيانات آلام الركبة GWAS هو: ukb-b-16254، بحجم عينة يبلغ 461857، وتشمل 98704 مريضا يعانون من ألم الركبة، بالإضافة إلى مجموعة ضابطة تبلغ 363153، مع إجمالي عدد SNPs يبلغ 9851867. يتم أيضا الحصول على بيانات GWAS لاضطرابات المشي من قاعدة بيانات openGWAS الخاصة ب IEU، ومجموعة بيانات شذوذات المشي رقم معرف GWAS هو: finn-b-R18 _ABNORMALITI_GAIT_MOBIL، بحجم عينة يبلغ 210717، بما في ذلك 1348 حالة من تشوهات المشي ومجموعات 209369 من التحكم، وتم تحديد 1638044 مواقع SNPs. التفاصيل موضحة في الجدول 1.

للحصول على بيانات قابلة لإعادة الإنتاج، يجب على الباحثين البحث في كل معرف GWAS مدرج أعلاه في قاعدة بيانات IEU OpenGWAS، وتحميل ملف الملخص-الإحصائي المقابل لGWAS حيثما توفر، وحفظ الملف كنص محلي أو ملف مفصول بفاصل قبل استيراده إلى R. في سير العمل الحسابي المستخدم في هذه الدراسة، تم استيراد ملفات ملخص-إحصائية محلية لنظام GWAS إلى R باستخدام read.table() أو read.csv()، حسب تنسيق ملف المصدر. تم اعتماد سير العمل القائم على الملفات المحلية لضمان إمكانية معالجة نفس ملفات الملخص-الإحصائية GWAS التي تم تنزيلها بشكل متكرر باستخدام تعيينات أعمدة ثابتة وأوامر تحليل متطابقة. الأوامر الدقيقة R لاستيراد ومعالجة الملفات المحلية متوفرة في الملف التكميلي 1.

الفحص وإدخال المرشحين الوطنيين

تم استخراج SNPs المرتبطة بشكل ملحوظ بكل نمط تعرض من مجموعات بيانات ملخص GWAS باستخدام عتبة دلالة جينومية واسعة النطاق تبلغ P < 5 × 10⁻⁸. إذا كان أقل من ثلاثة SNPs مستقلين متاحين عند هذه العتبة، يتم تخفيف العتبة إلى P < 5 × 10⁻⁶ للحصول على أدوات وراثية كافية لتحليل الرنين المغناطيسي. تم تقييم عدم توازن الربط باستخدام r² < 0.001 ضمن نافذة زمنية قدرها 10,000 كيلوبايت لضمان استقلالية المتغيرات الأداتية.

على مستوى السكربت، تم استيراد إحصائيات ملخص GWAS للتعرض أولا إلى R وتصفيتها وفقا لحد قيمة P المحددة مسبقا. تم بعد ذلك حفظ SNPs المختارة كملف تعرض وتنسيقها وفقا لمتطلبات حزمة TwoSampleMR. شملت خرائط الأعمدة المطلوبة معرف SNP، تقدير التأثير، الخطأ المعياري، أليل التأثير، أليل آخر، وقيمة P. في سير عمل R القابل للتكرار، تم استيراد بيانات التعريض باستخدام read_exposure_data()، وتم إجراء تجميع LD عن طريق تعيين clump = TRUE. هيكل أوامر R المقابل كان: read_exposure_data(اسم الملف = "exposure.csv"، sep = ","، snp_col = "rsids"، beta_col = "beta"، se_col = "sebeta"، effect_allele_col = "alt"، other_allele_col = "ref"، pval_col = "pval"، clump = TRUE). تم تعديل أسماء الأعمدة وفقا لملفات GWAS التي تم تنزيلها فعليا.

تم تقييم قوة الجهاز باستخدام إحصائية F:

F = [(N − k − 1)/k] × [R²/(1 − R²)],

حيث يمثل N حجم العينة، وk يمثل عدد المتغيرات الأداة، وR² يمثل نسبة تباين التعرض المفسرة بواسطة SNPs. تم حساب R² كالتالي

Σ[2 × MAF × (1 − MAF) × β²/(SE² × N)],

حيث أن MAF هو تردد الأليل الصغير، β هو تقدير تأثير الأليل، وSE هو الخطأ المعياري.

تم الاحتفاظ بوحدات SNP ذات إحصائيات F > 10. تم تحديد SNPs المرتبطة بعوامل محتملة مربكة، بما في ذلك التشوهات التشريحية الخلقية ومؤشر كتلة الجسم، باستخدام PhenoScanner V2 وإزالتها.

تم استخراج SNPs للنتائج عن طريق دمج SNPs المتجمعة مع ملف ملخص-إحصائي GWAS المحلي المقابل وفقا لمعرف SNP. إذا أبلغ الناتج عن GWAS −log10(P)، يتم تحويل قيمة P باستخدام P = 10^(−LP). ثم تم استيراد SNPs المستخرجة باستخدام read_outcome_data(). تم تنسيق مجموعات بيانات التعرض والنتائج باستخدام harmonise_data() لمحاذاة أليلات التأثير وإزالة SNPs غير المناسبة للتحليل بالرنين المغناطيسي. تم استبعاد SNPs باليندرومية ذات توجه الأليل الغامض، وتم الاحتفاظ بSNPs ذات mr_keep = TRUE لتحليلات MR لاحقة¹¹. النص التفصيلي لاستخراج SNP الناتج وتحويل قيمة P والتوحيد متوفر في الملف التكميلي 1.

التحليل الإحصائي

تم تقدير التأثيرات السببية بين التعرض ومتغيرات النتائج في هذه الدراسة باستخدام خمس طرق: طريقة IVW، طريقة الانحدار MR-Egger، طريقة الوسيط الموزون، طريقة النمط الموزون، وطريقة النمط البسيط¹². تعتبر طريقة IVW الطريقة القياسية لتحليل الرنين المغناطيسي، والتي تفترض أن جميع المتغيرات الآلية صالحة، وتعتمد على مبدأ دمج تقديرات نسبة والد، ويستخدم نموذج التأثيرات العشوائية إذا كان هناك تغاير، ويستخدم نموذج التأثيرات الثابتة إذا لم يكن هناك تغاير. تكتشف طريقة الانحدار MR-Egger إمكانية التعدد التعاوني وتأخذ في الاعتبار وجود مصطلح اعتراض في الانحدار¹³. تتطلب طريقة الوسيط المرجح أن تكون أكثر من 50٪ من المتغيرات الآلية SNPs صالحة. تتطلب طريقة الوضع المرجح أحجام عينات أصغر من الطرق الأخرى وتضمن انحياز أقل ومعدلات خطأ من النوع الأول أقل. تسمح طريقة الوضع البسيط بتجميع SNPs ذات التأثيرات المتشابهة بناء على ما إذا كانت التأثيرات السببية المقدرة متشابهة أم لا¹⁴.

على مستوى الوظيفة، أجريت تحليلات MR باستخدام دالة mr() في حزمة TwoSampleMR. تم تحديد قائمة الطرق المستخدمة في التحليل على النحو التالي: c("mr_ivw"، "mr_egger_regression"، "mr_weighted_median"، "mr_simple_mode"، "mr_weighted_mode"). تم توليد نسب الاحتمالات وفترات ثقة 95٪ باستخدام generate_odds_ratios(). يتم توفير تنفيذ R الكامل، بما في ذلك قائمة الطرق الدقيقة وأوامر حفظ الإخراج، في الملف التكميلي 1.

مراقبة الجودة وتحليل الحساسية

لاختبار استقرار وموثوقية نتائج الرنين المغناطيسي، تم أولا استخدام اختبار كوكران Q لتقييم تباين SNPs، وإذا كان اختبار كوكران Q ذا دلالة إحصائية فقد أشار إلى وجود تباين ذي دلالة في النتائج المحللة¹⁵. ثانيا، استخدم اختبار اعتراض MR-Egger لتقييم التنوع الأفقي المحتمل. يشير اعتراض MR-Egger ذو دلالة إحصائية إلى وجود تعدد الذرات الأفقي الاتجاهي في تحليل الرنين المغناطيسي. ثالثا، تم تطبيق العشوائية المندلية لجمع المجموع والعدد الشاذ المتبقي متعدد الأبعاد (MR-PRESSO) للبحث عن وجود SNPs الشاذة في النتائج، وإذا كانت موجودة، تم استبعادها وإعادة تحليلها¹⁶. رابعا، تم تطبيق طريقة "ترك واحد خارج" لاختبار متانة النتائج، وتم حساب التأثير المشترك لبقية SNPs باستبعاد SNPs واحدة تلو الأخرى لتقييم تأثير SNPs الفردية على العلاقة بين التعرض ومتغيرات النتائج¹⁷.

على مستوى الوظيفة، تم تقييم التباين باستخدام mr_heterogeneity()، وتم تقييم التنوع الاتجاهي باستخدام mr_pleiotropy_test()، وتقييم تأثيرات SNP المفرد باستخدام mr_singlesnp()، وأجري تحليل ترك واحد خارج باستخدام mr_leaveoneout(). تم إنشاء مخططات تشخيصية باستخدام mr_scatter_plot()، mr_funnel_plot()، و mr_leaveoneout_plot(). تم إجراء MR-PRESSO عندما كان عدد المتغيرات الآلية المتاحة كافيا؛ إذا لم يكن بالإمكان إجراء MR-PRESSO بسبب نقص عدد SNPs، كان يتم تسجيل ذلك والإبلاغ عنه. تم إجراء إجراء تحليل MR وإجراءات مراقبة الجودة باستخدام الإصدار R 4.3.2 وحزمة TwoSampleMR الإصدار 0.5.8، مع مستوى دلالة α = 0.05. الأوامر الكاملة مقدمة في الملف الإضافي 1.

حفظ المخرجات، التحقق من الجودة، والتقارير

تم حفظ جميع المخرجات الوسيطة والنهائية لإعادة إنتاجها. شملت هذه العوامل التعريض المختارة، وأدوات التعرض المتجمعة، ونتائج SNP المستخرجة، ومجموعات البيانات المتوافقة، وقوائم SNP المحتفظة، وتقديرات MR، ونسب الأرجاء، ونتائج التباين، ونتائج اعتراض MR-Egger، ونتائج MR-PRESSO، ونتائج SNP الفردي، ونتائج ترك واحد خارج، ومخططات تشخيصية، ومعلومات جلسات R. قبل التقرير، تم فحص SNPs المحتفظ بها من حيث الدلالة على مستوى الجينوم، واستقلالية LD، وحالة توافق الأليليل، والمتغيرات الباليندرومية الغامضة، وقوة الأجهزة، والتباين، وتعدد الروبوتية الأفقية الاتجاهية، وSNPs الشاذة. تم الإبلاغ عن نتائج التصوير بالرنين المغناطيسي الأولية كنسب احتمالات مع فترات ثقة 95٪ وقيم P. تم الإبلاغ عن نتائج تحليل الحساسية والمخططات التشخيصية في الجداول والأشكال التكميلية. يتم توفير أوامر حفظ الإخراج والتحقق من الملفات في الملف الإضافي 1.

النتائج

تم تحليل بيانات GWAS لثلاثة أنماط ظاهرية للألم العضلي الهيكلي الخاصة بموقع، بما في ذلك آلام أسفل الظهر، وألم الورك، وألم الركبة، ونمط ظاهري واحد من اضطرابات المشي في هذه الدراسة. يمكن العثور على معلومات مفصلة عن الوريد لكل عامل تعرض في الملف التكميلي 3.

تحليل بالرنين المغناطيسي لألم الجهاز العضلي الهيكلي الخاص بالموقع على اضطرابات المشي

في هذه الدراسة، تم استبعاد SNPs ذات عدم توازن الربط والبنية الب...

المناقشة

يمكن أن تغير اضطرابات المشي استراتيجيات تحميل وحركة مفاصل الأطراف السفلية، مما قد يسرع عمليات 821+_osteoarthritic ويساهم في قيود وظيفية تؤثر بشكل كبير على جودة الحياة. وقد جذب هذا الارتباط السريري اهتماما متزايدا بالتفاعل بين خلل المشي وألم الجهاز العضلي الهيكلي. تشير الدراسات الرصدية إلى أن الألم العضلي الهيكلي المرتبط بالموقع قد يرتبط بضعف في المشي؛ على سبيل المثال، أظهرت دراسة مقطعية شملت كبار السن الذين يعيشون في المجتمع وهم تبلغ أعمارهم 65 عاما فأكثر أن الألم متعدد المواقع مرتبط ببط...

المواد

قائمة المواد المستخدمة في هذه المقالة

الاسم	الشركة	رقم فهرسي	التعليقات
قاعدة بيانات IEU OpenGWAS	وحدة علم الأوبئة التكاملي التابعة لمجلس البحوث الطبية (MRC)، جامعة بريستول	لا يوجد	مصدر إحصائيات ملخص GWAS للتعرض والأنماط الظاهرية للنتائج.
حزمة MR-PRESSO R	ماري فيربانك / حزمة R	لا يوجد	يستخدم للكشف عن القيم الشافية وتقييم التنوع الأفقي.
فينو سكانر V2	جامعة كامبريدج	V2	يستخدم المورد على الويب لتحديد SNPs المرتبطة بالمراوغين المحتملين.
R	مؤسسة R للحوسبة الإحصائية	الإصدار 4.3.2	بيئة الحوسبة الإحصائية المستخدمة لتحليل التصوير المغناطيسي ومراقبة الجودة.
حزمة TwoSampleMR R	وحدة علم الأوبئة التكاملي التابعة لمجلس البحوث الطبية (MRC)، جامعة بريستول	الإصدار 0.5.8	يستخدم في تحليلات IVW، MR-Egger، الوسيط الموزون، النمط المرزون، والنمط البسيط.