دراسة الارتباط الشامل بالنسخ التكاملي والعشوائية المندلية تحدد LYNX1 و MS4A14 كأهداف علاجية لالتهاب العظم

Summary

تدمج هذه الدراسة تحليل الارتباط على مستوى النسخ النصي والعشوائية المندلية لتحديد أهداف الأدوية لالتهاب العظم النقي بشكل منهجي. وباستخدام قواعد البيانات الجينية، حدد المؤلفون جينات LYNX1 وMS4A14 كجينات مرشحة رئيسية، مما أبرز المسارات الجزيئية الحيوية وقدم استراتيجيات علاجية دقيقة لإدارة التهابات العظام.

Abstract

التهاب العظام العظمي هو مرض معدي شديد يتميز بالتهاب عميق في العظم ونخاع العظم، ويسببه بشكل رئيسي مسببات أمراض بكتيرية مثل المكورات العنقودية الذهبية (S. aureus). ومع ذلك، فإن أنظمة العلاج الحالية تعاني بشكل كبير من ظهور مقاومة المضادات الحيوية والعدوى المتكررة، مدفوعة بتكوين الأغشية الحيوية القوي واستمرار البكتيريا داخل الخلية. لذا، هناك حاجة ملحة لتحديد أهداف علاجية جديدة تتجاوز المضادات الميكروبية التقليدية. تم تطبيق نهج حسابي تكاملي يجمع بين دراسة الارتباط على مستوى النسخ (TWAS) وتحليل العشوائية المندلية القائمة على الملخص (SMR) لتحديد جينات الحساسية السببية الكامنة وراء التهاب العظم النقي بشكل منهجي. تم استخدام المتغيرات الجينية واسعة النطاق كمتغيرات أداة للتنبؤ بتعبير الجينات عبر الأنسجة، مما أتاح استكشاف الروابط السببية بين هذه الأهداف ومخاطر المرض. حدد الإطار التكاملي LYNX1 و MS4A14 كجينات مرشحة رئيسية وأهداف علاجية محتملة. على وجه التحديد، ارتبط LYNX1 بزيادة القابلية لالتهاب العظم النقي، بينما أظهر MS4A14 خصائص وقائية محتملة. تسلط هذه النتائج الضوء على الأدوار التنظيمية لهذه الجينات في الاستجابة المناعية للمضيف والتعديل الالتهابي أثناء عدوى العظام. تربط هذه الدراسة الفجوة بين نتائج الارتباط الجينومي العام والتفسير البيولوجي، مما يعزز فهم الأساس الجيني لالتهاب العظم والنقي، ويدعم تطوير استراتيجيات طب دقيقة مستهدفة لتعزيز دفاع المضيف والتغلب على المقاومة العلاجية.

Introduction

التهاب العظم هو مرض معدي شديد يتميز بالتهاب عميق في العظم ونخاع العظم. غالبا ما يبدأ بسبب غزو البكتيريا، حيث تعتبر المكورات العنقودية الذهبية (S. aureus) مرضا شائعا عبر العروض الدموية الحادة والمزمنة ^1,2. بينما يكون التهاب العظام الذهبي بارزا، فإن التهاب العظم النخاعي غالبا ما يكون متعدد الميكروبات، والقضاء على S. aureus وحده لا يحل دائما العدوى بسبب استمرار الكائنات الدقيقة غير القابلة للزراعة والتفاعلات الديناميكية مع الغشية الحيوية^3،4. في المضيفات الحساسة، تستخدم مسببات الأمراض استراتيجيات تفادي معقدة، بما في ذلك الاستمرارية داخل الخلايا داخل البلعم وتكوين أغشية حيوية كثيفة. تتكون هذه الأغشية الحيوية من الحمض النووي خارج الخلية، والبوليساكاريدات، والبروتينات المتخصصة، وتخلق حواجز فيزيائية قوية تقاوم دفاعات المضيف المناعية وتحد من اختراق العوامل المضادة للميكروبات التقليدية ^5,6. وبالتالي، على الرغم من أن الإدارة السريرية عادة ما تتضمن إزالة الألياف الجراحية العدوانية والعلاج بالمضادات الحيوية لفترة طويلة، فإن وجود الأغشية الحيوية المرنة والمجتمعات متعددة الميكروبات يقلل بشكل كبير من فعالية العلاج، مما يؤدي إلى استمرار المراضة وتكرار العدوى^. يبقى توضيح البنية الجزيئية الكامنة وراء تفاعلات المضيف والممرض أمرا بالغ الأهمية لتحديد أهداف علاجية جديدة⁹.

لقد طورت دراسات الارتباط الجينومي (GWAS) فهم التهاب العظم النخاعي من خلال تحديد مواقع وراثية متعددة مرتبطة بحساسية المرض¹⁰. ومع ذلك، وبسبب عتبات الاختبار المتعددة الصارمة والقيود في القوة الإحصائية وحجم العينة، من المرجح أن تظل العديد من المواقع الحرجة المؤثرة على حساسية التهاب العظم النخاعي^{غير مكتشفة}. علاوة على ذلك، لا يزال تفسير نتائج GWAS تحديا، حيث توجد معظم المتغيرات المعروفة في مناطق غير مشفرة أو بين الجينات، مما يشير إلى تأثيرات تنظيمية على التعبير الجيني بدلا من تغييرات مباشرة في بنية البروتين¹¹. بالإضافة إلى ذلك، تفتقر بيانات GWAS الأولية إلى خصوصية الأنسجة المباشرة، مما يستلزم إجراء تحليلات تكميلية لاحقة لتحديد السياقات الخلوية التي تؤثر فيها هذه المتحورات. تعالج دراسات الارتباط الواسع النطاق بالنسخات (TWAS) هذه القيود من خلال دمج التعبير الجيني المتوقع وراثيا من لوحات المراجع في مجموعات بيانات GWAS الحالية، مما يمكن من التعرف المنهجي على الجينات المرتبطة بالمرض ذات الصلة وظيفيا عبر عدة أنسجة^.

لتعزيز الاستدلال السببي، يستخدم العشوائية المندلية (MR) المتغيرات الجينية كمتغيرات أداتية، مستفيدة من التخصيص العشوائي للأليلات عند الحمل لتقليل الخلط وعكس السببية^13،14. تربط تحليلات مواقع الصفات الكمية للتعبير (eQTL) متغيرات GWAS بنسخ^{الجينات 15}، بينما تدمج العشوائية المندلية القائمة على الملخص (SMR) مجموعات بيانات GWAS و QTL لإعطاء الأولوية للجينات السببية وتمييزها عن الارتباطات الناتجة عن اختلال التوازن في الربط من خلال اختبار التغاير¹⁶.

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

Protocol

جميع الإحصائيات الملخصة المستخدمة في تحليلات العشوائية المندلية (MR) ودراسة الارتباط الواسع للنسخ (TWAS) تم اشتقاقها فقط من مجموعات بيانات منشورة سابقا وغير محددة. يتم توثيق الموافقة الأخلاقية والموافقة الفردية على الدراسات الأصلية في منشوراتهم الخاصة. وبالتالي، تم التنازل عن الموافقة الأخلاقية الإضافية لهذه الدراسة التنقيب عن البيانات من قبل مجلس المراجعة المؤسسية لمستشفى تونغده في مقاطعة تشجيانغ (Zhe Tongde Lunshen 2024 [Yan] رقم 028-JY). الأدوات المستخدمة في هذا البحث مدرجة في جدول المواد.

1. جمع ومعالجة بيانات RNA-seq

تم الحصول على بيانات النسخ من قاعدة بيانات التعبير الجيني (GEO) (GSE272198) لتقييم حفظ المسارات المناعية الفطرية عبر أنواع الثدييات للتحقق الأولي¹⁷. أصيبت البلعم المشتقة من نخاع العظم (BMDMs) ب S. aureus (تعدد العدوى، MOI = 10) لمدة ساعة واحدة، تلاها العلاج بالليزوستافين (20 ميكروغرام/مل) وجنتاميسين (50 ميكروغرام/مل) لإزالة البكتيريا خارج الخلية. بعد ثلاث غسلات باستخدام محلول ملحي مخزن بالفوسفات (PBS)، تم زراعة BMDMs لمدة 24 ساعة، وتحللها في كاشف استخلاص RNA كلي، وتم تسلسلها.

تم تقييم جودة الحمض النووي الريبي باستخدام نظام الرحلان الكهربائي الآلي لضمان السلامة والنزاهة. تم إعداد المكتبات من ثلاث تجارب مستقلة وتم تسلسلها على منصة تسلسل عالية الإنتاجية. تم محاذاة القراءات الخام مع جينوم الفأر (GRCm38, mm10) باستخدام STAR (v2.7.10a). تم تحديد الجينات المعبر عنها بشكل تفاضلي (DEGs) باستخدام DESeq2 (v1.38.0). للتخفيف من الإيجابيات الكاذبة، تم تعريف الدلالة الإحصائية بأنها قيمة p معدلة (FDR) < 0.05 و |log₂ تغير الطية| > 1. تم إجراء تحليل الأنطولوجيا الجينية (GO) باستخدام clusterProfiler (v4.6.0)، وأجري تحليل إثراء مجموعات الجينات (GSEA) باستخدام GseaVis (v0.0.5). تم توليد خرائط الحرارة باستخدام حزمة pheatmap (v1.0.12) في R (v4.2.0).

تحليل TWAS

تم الحصول على بيانات تسلسل الحمض النووي الريبي الكامل للدم وتسلسل الجينوم الكامل (WGS) من مشروع تعبير النمط الجيني للأنسجة (GTEx) (V8)¹⁸. تم استخدام نماذج تعبير جينات مدربة مسبقا من مستودع عام (https://doi.org/10.5281/zenodo.3842289). تم استرجاع إحصائيات ملخص التهاب العظم لمرض TWAS من اتحاد فينجين، الذي يضم 2,336 حالة و473,264 حالة ضابطة¹².

تم إجراء TWAS باستخدام ثلاث خوارزميات: حساب أنسجة المفاصل (JTI)، PrediXcan¹⁹، وUTMOST ^12,20. تقدر JTI تشابه التعبير الجيني وسهولة الكروماتين فوق الجيني لتحسين دقة التنبؤ. يطبق PrediXcan الانحدار الصافي المرن مع التحقق المتقاطع الخماسي، بينما يعزز UTMOST الدقة من خلال الاستفادة من بيانات التعبير النسيجية المتعددة باستخدام LASSO المجموعات المتناثرة. الإطار المعدل ل UTMOST الذي وصفه تشو ^وآخرون يحدد معايير فائقة للتقدير غير المتحيز. تم الاحتفاظ بالجينات ذات درجات التحقق المتقاطع المستقرة—المعرفة مسبقا كمعامل ارتباط r > 0.1 ودلالة تنبؤية p < 0.05²¹—كجينات قابلة للنسب. تم إنشاء نماذج النسخ الكامل للدم باستخدام مصفوفات التغاير SNP من مجموعة بيانات مرجعية ل 1000 جينوم.

تم تحليل الروابط بين التعبير الجيني المتوقع وخطر التهاب العظم النقيعي لاحقا. لأخذ الاختبارات المتعددة في الاعتبار، تم تحديد الدلالة الإحصائية ل TWAS بشكل أساسي باستخدام عتبة معدل الاكتشاف الكاذب (FDR) < 0.05. نظرا لطبيعة توليد الفرضيات لهذه الدراسة متعددة المراحل، تم أيضا أولوية المواقع التي تحقق عتبة موحية (اسمية) ل p < 0.05 لتحليلات العشوائية المندلية (SMR) وتحليل التنقل المشترك. تهدف هذه الاستراتيجية التكاملية إلى تعظيم استغلال المحركات التنظيمية المحتملة مع الاعتماد على التحقق المتقاطع متعدد الأوميكس (TWAS + SMR) لضمان متانة المرشحين المفضلين.

تحليل SMR

التزمت هذه الدراسة بإرشادات تعزيز الإبلاغ عن الدراسات الرصدية في علم الأوبئة (STROBE)²². لتعريف النمط الظاهري الذي يمثل الاستعداد الجيني لخلل الميتوكوندريا (ويطلق عليه لاحقا "الميتوديات" لأغراض التحليل)، تم استخراج نسخ تمثل جميع الجينات المعروفة المرتبطة بالميتوكوندريا من قاعدة بيانات MitoCarta3.0²³. وقد خدمت هذه المجموعة الجينية كأساس محدد مسبقا ومستنعا بعلم الأحياء للتنبؤ بالمخاطر متعددة الجينات لاحقا. جميع التفسيرات الوظيفية اللاحقة المتعلقة ب "الميتوديات" مشتقة من هذا الاستنتاج الحسابي ويجب اعتبارها تنبؤية ومولدة فرضيات.

تم توليد أجهزة مواقع السمات الكمية للتعبير (eQTL) باستخدام متغيرات ضمن 1000 كيلوبايت من تسلسلات الترميز (cis-eQTLs). تم الحصول على الإحصائيات الملخصة من اتحاد eQTLGen وGTEx V8²⁴. تم اختيار ما مجموعه 8,932,843 SNP مرتبطة ب 1,013 نسخا مرتبطة بالميتوديهات بناء على عتبة دلالة على مستوى الجينوم P < 5E-8. تم الحصول على إحصائيات GWAS الأساسية لنتائج التهاب العظم النخاعي من FinnGen²⁰.

تم إجراء تحليل العشوائية المندلية (SMR) المستندة إلى بيانات ملخصة باستخدام SMR (الإصدار 1.0.3) مع معلمات افتراضية لتقدير الارتباطات متعددة الأبعاد بين صفات التعبير الجيني ونتائج التهاب العظم النقي. يمثل التأثير السببي beta_mitodys–التهاب العظم والنخاع حجم التأثير التقديري لاضطراب الميتوكوندريا على التهاب العظم والنقي، ويتم حسابه كالتالي:

تمثل نسب الاحتمالات (ORs) التغير لكل وحدة زيادة لوغاريتمية طبيعية واحدة في مستويات التعبير الجيني الموحدة. تم تقييم التواجد المشترك بشكل إضافي باستخدام اختبار التغاير في الأجهزة التابعة (HEIDI).

Results

توصيف تسلسل الحمض النووي الريبي الضخم والتخصيبات الوظيفية

وفقا للمنهجية المستخدمة لعزل BMDMs غير المصابة (الضابطة) مقابل مصابة ب S. aureus، تم تحديد 3,775 DEGs في المجموعة التجريبية، تضم 1,686 جينا مرتفعا و2,089 جينا منخفضا. كشف تحليل GO عن شبكة غنية بقوة لتنشيط المناعة الفطرية، بما في ذلك التنظيم الإيجابي (GO:0045089) والتنشيط (GO:0002218) للاستجابة المناعية الفطرية، مما يعكس آليات الدفاع المعتمدة عن مسببات الأمراض المدفوعة بالسيتوكي...

Discussion

من خلال دمج تحليلات TWAS وSMR عبر مجموعات البيانات الترجمية، حددت هذه الدراسة شبكة أساسية من المواقع الجينية المرتبطة بأنماط التهاب العظم والنخاع الظاهري. يمثل تحديد LYNX1 وMS4A14 كأهداف جينية مرشحة مساهمة رئيسية، حيث يوفر رؤى ميكانيكية تتجاوز نتائج GWAS التقليدية من خلال دمج السياق النصي. يبرز هذا النهج التكاملي إمكانات استراتيجيات الطب الدقيق لتعزيز الإدارة السريرية للالتهابات العظمية الشديدة.

منهجيا، يمثل استخدام العشوائية المندلية لاستنتاج العلاقات السببية...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
clusterProfiler	فريق Guangchuang Yu (جامعة صن يات سين)	NA	يُستخدم لتحليل الإثراء الوظيفي لـ Gene Ontology (GO) للجينات التي تختلف في التعبير.
DESeq2	مشروع Bioconductor (مصدر مفتوح)	NA	يُستخدم لتحديد الجينات التي تختلف في التعبير (DEGs)
GSEA / GseaVis	معهد Broad / مجتمع المصدر المفتوح	NA	يُستخدم لتحليل إثراء مجموعة الجينات (GSEA) لتحديد مجموعات الجينات التي يتم التعبير عنها بشكل منسق. من المحتمل أن يتم استخدام GseaVis للتصور.
اختبار HEIDI (الاختلاف في الأدوات التابعة)	وظيفة متكاملة داخل برنامج SMR	NA	يُستخدم لاختبار ما إذا كانت الارتباطات المرصودة في SMR مدفوعة باختلال التوافق أكثر من المتغير السببي المشترك. P_HEIDI > 0.05 يشير إلى عدم وجود اختلاف كبير
JTI (الاستدلال المشترك للأنسجة)	Zhou et al. (نموذج مصدر مفتوح)	NA	خوارزمية TWAS تقدر تشابه التعبير الجيني وقابلية الوصول إلى الكروماتين الوراثي لتحسين دقة التنبؤ.
pheatmap	Raivo Kolde (حزمة R)	NA	يُستخدم لإنشاء خرائط حرارية للجينات التي تختلف في التعبير لتصور البيانات.
PrediXcan	Gamazon et al. (نموذج مصدر مفتوح)	NA	خوارزمية TWAS تطبق انحدار الشبكة المرنة مع التحقق المتبادل الخماسي لتقدير التعبير الجيني واختبار الارتباط بمخاطر المرض.
SMR (التحليل العشوائي المبني على الملخص)	تم تطويره من قبل فريق في جامعة Fudan (برنامج مصدر مفتوح)	NA	يُستخدم لإجراء تحليل التمايز العشوائي المبني على الملخص لاستكشاف الارتباطات السببية بين التعبير الجيني ومخاطر التهاب العظم والنقي. تم استخدام الإصدار 1.0.3.
STAR	Alexander Dobin (مختبر Cold Spring Harbor)	NA	يُستخدم لمحاذاة قراءات RNA-seq مع الجينوم المرجعي (جينوم الفأرة GRCm38/mm10).
UTMOST	Hu et al. (نموذج مصدر مفتوح)	NA	خوارزمية TWAS تعزز دقة التنبؤ من خلال الاستفادة من بيانات التعبير متعددة الأنسجة باستخدام إطار عمل LASSO الجماعي المتناثر.