Exploring ubiquitination-related biomarkers in osteoarthritis based on single-cell and transcriptomic data
تدمج هذه الدراسة بيانات النسخ المتماثل أحادي الخلية والنسخ المتماثل الكلي لتحديد خمسة مؤشرات حيوية مرتبطة باليوبيكويتين (MYC، وTLR7، وCX3CR1، وPER1، وGADD45B) تشارك في إمراضية الفصال العظمي، مما يكشف عن إثراء هذه المؤشرات في مسار TNF-alpha/NF-kappaB، وارتباطها بمجموعات خلايا مناعية محددة، وارتفاع تعبيرها في الخلايا الليفية الزليلية، مما يقدم بذلك أهدافاً علاجية جديدة لهذا المرض.
المؤلفون الأصليون:Tian Lan, Huan Dai, DeGuang Li, ChongYu Zhao, XiaoYu Mu, Li Li, Chao Wang
تخيل جسدك كمدينة صاخبة وعالية التقنية. داخل هذه المدينة، كل خلية هي عامل، وللحفاظ على سير العمل بسلاسة، يعتمدون على نظام متطور لإدارة النفايات ومراقبة الجودة. أحد أهم الأدوات في هذا النظام هو علامة جزيئية صغيرة تسمى "يوبيكويتين" (ubiquitin). فكر في اليوبيكويتين كملصق ملاحظة صغير أو ملصق "لا تستخدم" تلصقه الخلية على البروتينات. أحياناً، تخبر هذه الملصقة الخلية بإعادة تدوير بروتين ما؛ وفي أحيان أخرى، تشير إلى أن البروتين معطل ويجب التخلص منه. عندما يعمل نظام الوسم هذا بشكل مثالي، تظل المدينة صحية. ولكن عندما يتم وضع هذه الملصقات بشكل غير صحيح — سواء كان ذلك بكثرة، أو بقلة، أو على عناصر خاطئة — تبدأ البنية التحتية للمدينة في الانهيار.
الآن، تخيل جزءاً محدداً من هذه المدينة: المفاصل. بمرور الوقت، قد تتآكل الوسائد الناعمة والمبطنة (الغضاريف) التي تسمح لركبتيك ووركيك بالانزلاق، مما يؤدي إلى حالة مؤلمة تسمى الفصال العظمي (OA). لفترة طويلة، عرف الأطباء أن الفصال العظمي مشكلة ناتجة عن التآكل والالتهاب، لكنهم لم يفهموا تماماً "الخلل" الجزيئي المسبب لهذا التدهور. تشير الأبحاث الحديثة إلى أن نظام وسم اليوبيكويتين قد يكون معطلاً في حالات الفصال العظمي، مما يترك خلايا المفاصل في حالة ارتباك وغير قادرة على إصلاح نفسها. ويتساءل العلماء الآن: إذا تمكنا من العثور على البروتينات المحددة التي تحصل على الوسوم الخاطئة، فهل يمكننا إصلاح النظام وإيقاف الألم؟
هذه هي المهمة بالضبط لدراسة جديدة أجراها باحثون من مستشفيات في كونمينغ، الصين. لقد قرروا استقصاء الفصال العظمي من خلال النظر في "علامات اليوبيكويتين" عبر مجهر عالي القدرة مصنوع من البيانات. وبدلاً من مجرد النظر إلى المفصل بأكمله، استخدموا أداتين قويتين: "علم النسخ المتعدد للكتلة" (bulk transcriptomics)، وهو يشبه التقاط صورة ضبابية للمدينة بأكملها لرؤية الاتجاهات العامة، و"تسلسل الحمض النووي الريبوزي وحيد الخلية" (single-cell RNA sequencing)، وهو يشبه التكبير لالتقاط صورة حادة لكل عامل بمفرده في المدينة لرؤية ما يفعلونه بالضبط.
بدأ الباحثون بجمع البيانات من ثلاث قواعد بيانات عامة مختلفة، والتي احتوت على معلومات جينية من مفاصل الركبة لأشخاص مصابين بالفصال العظمي وأشخاص يتمتعون بمفاصل سليمة. كما أعدوا قائمة تضم 79 جيناً معروفة بأنها تشارك في نظام وسم اليوبيكويتين. ومن خلال مطابقة هذه القوائم، وجدوا 7 جينات كانت جزءاً من نظام الوسم وتتصرف بغرابة لدى مرضى الفصال العظمي. لكنهم لم يتوقفوا عند هذا الحد؛ فقد استخدموا خدعة إحصائية ذكية للعثور على 321 جيناً آخر كانت "تتحدث" مع تلك الجينات السبعة الغريبة، مما خلق شبكة هائلة من التفاعلات.
ولإيجية الجناة الحقيقيين في هذه الشبكة، استخدم الفريق خوارزميات الكمبيوتر (تعلم الآلة) لتضييق القائمة. قاموا بتشغيل البيانات عبر نظامين مختلفين للتصفية، يشبهان استخدام نوعين مختلفين من أجهزة الكشف عن المعادن للعثور على أكثر القطع قيمة. أدت هذه العملية إلى تقليص القائمة إلى خمسة جينات رئيسية فقط: MYC، وTLR7، وCX3CR1، وPER1، وGADD45B. واختبر الباحثون هذه الجينات الخمسة في مجموعتين مختلفتين من المرضى للتأكد من أنها ليست مجرد صدفة. ووجدوا أن هذه الجينات كانت ممتازة في التمييز بين المفاصل السليمة والمفاصل المصابة بالفصال العظمي، حيث عملت مثل نظام إنذار موثوق.
ثم تعمقت الدراسة لمعرفة ما تفعله هذه الجينات بالفعل. اكتشفوا أن هذه الجينات الخمسة كانت تشارك بقوة في مسار التهابي رئيسي يسمى مسار TNF-α/NF-κB. يمكنك التفكير في هذا المسار كصافرة الطوارئ في المدينة؛ فعندما تنطلق، تخبر الجهاز المناعي بالاندفاع للقتال. وفي حالات الفصال العظمي، يبدو أن هذه الصافرة عالقة في وضع التشغيل "On"، مما يسبب التهاباً مزمناً يدمر المفصل. كما نظر الباحثون في الخلايا المناعية في المفصل ووجدوا أن أنواعاً معينة، مثل الخلايا البلازمية وخلايا بلعمية محددة، كانت أكثر شيوعاً في مجموعة الفصال العظمي، مما يغذي النيران على الأرج
ملخص تقني: استكشاف المؤشرات الحيوية المرتبطة باليوبيكويتين في التهاب المفاصل العظمي بناءً على بيانات الخلية الواحدة والنسخ المتماثل
بيان المشكلة يعد التهاب المفاصل العظمي (OA) مرضاً تنكسياً شائعاً في المفاصل، يتميز بتدهور الغضروف، والتهاب الغشاء الزلالي، وإعادة تشكيل العظام تحت الغضروفية. وبينما تُعرف التعديلات ما بعد الترجمة، وتحديداً عملية اليوبيكويتين (Ubiquitination)، بدورها في تنظيم مسارات رئيسية مرتبطة بالتهاب المفاصل العظمي مثل مسار NF-κB وتوازن المصفوفة خارج الخلوية، إلا أن الآليات الدقيقة التي تربط خلل تنظيم اليوبيكويتين بالإمراض السريري لالتهاب المفاصل العظمي لا تزال غامضة. علاوة على ذلك، هناك نقص في الاستكشاف المنهجي للمؤشرات الحيوية المرتبطة باليوبيكويتين عبر المجموعات الخلوية غير المتجانسة داخل أنسجة مفاصل المصابين بالتهاب المفاصل العظمي. تهدف هذه الدراسة إلى سد هذه الفجوة من خلال تحديد الجينات المتميزة في التعبير المرتبطة باليوبيكويتين (UR-DEGs)، وفك رموز أصولها الخلوية، والتحقق من الأهداف العلاجية المحتملة.
المنهجية استخدمت الدراسة نهجاً متعدد الأوميكس يدمج بين النسخ المتماثل الكلي (Bulk Transcriptomics) وتسلسل الحمض النووي الريبوزي أحادي الخلية (scRNA-seq):
الحصول على البيانات: تم استرداد مجموعات بيانات النسخ المتماثل الكلي (GSE55235 كمجموعة تدريب، وGSE55457 كمجموعة تحقق) ومجموعة بيانات الخلية الواحدة (GSE152805) من مستودع التعبير الجيني (GEO). تألف جميع العينات من أنسجة الركبة الزلالية. وتم تجميع قائمة تضم 79 جيناً مرتبطاً باليوبيكويتين (URGs) من الأدبيات الموجودة.
التعبير التفاضلي والفرز: تم تحديد الجينات المتميزة في التعبير (DEGs) في مجموعة التدريب. وأدى تقاطع الـ DEGs مع الـ URGs إلى الحصول على جينات اليوبيكويتين المتميزة في التعبير (DE-URGs). ولتوسيع النطاق الجيني، تم إجراء تحليل ارتباط سبيرمان بين الـ DE-URGs والـ DEGs لتحديد الـ UR-DEGs.
التحليل الشبكي والوظيفي: تم إنشاء شبكة تفاعل البروتين-بروتين (PPI) لتحديد الـ UR-DEGs المتفاعلة. كما أُجريت تحليلات إثراء "مصطلحات الأنطولوجيا الجينية" (GO) و"موسوعة كيوطو للجينات والمسارات" (KEGG) لتحديد المسارات الوظيفية.
تحديد المؤشرات الحيوية: تم تطبيق خوارزميات التعلم الآلي، وتحديداً خوارزمية Boruta وخوارزمية LASSO، لفرز الجينات الرئيسية المرشحة. وتم تعريف الجينات التي سجلت مساحة تحت المنحنى (AUC) أكبر من 0.7 في كل من مجموعتي التدريب والتحقق، والتي أظهرت اتجاهات تعبير متسقة، كـ "مؤشرات حيوية".
التحليل اللاحق:
الارتشاح المناعي: استُخدم CIBERSORT لتحليل البيئة المناعية الدقيقة وربط المؤشرات الحيوية بنسب الخلايا المناعية.
تحليل الخلية الواحدة: تمت معالجة بيانات scRNA-seq باستخدام Seurat للتجميع وتوصيف أنواع الخلايا. وتم تقييم تعبير المؤشرات الحيوية عبر المجموعات الفرعية للخلايا، مع التركيز على المجموعة الفرعية ذات التعبير الأعلى. كما أُجري تحليل الزمن الكاذب (Pseudotime) باستخدام Monocle وتحليل التواصل الخلي باستخدام CellPhoneDB.
الشبكات التنظيمية: تم إنشاء شبكة RNA تنافسية داخلية (ceRNA) تشمل lncRNAs وmiRNAs والمؤشرات الحيوية. وتم التنبؤ بالتفاعلات بين الدواء والجين باستخدام DrugBank وDGIdb.
التحقق التجريبي: تم التحقق من التعبير للمؤشرات الحيوية المحددة سريرياً باستخدام تفاعل البوليميراز المتسلسل للنسخ العكسي الكمي (RT-qPCR) على عينات الأنسجة الزلالية من مرضى التهاب المفاصل العظمي والمجموعات الضابطة.
النتائج الرئيسية
تحديد الجينات: من بين 2,822 جينًا متميزًا في التعبير و79 جينًا مرتبطًا باليوبيكويتين، تم تحديد 7 جينات DE-URGs في البداية. ووسع تحليل الارتباط هذا العدد إلى 321 جينًا (UR-DEGs)، منها 120 جينًا أظهرت تفاعلات PPI كبيرة.
المؤشرات الحيوية: تم تأكيد خمسة جينات كمؤشرات حيوية قوية: MYC، وTLR7، وCX3CR1، وPER1، وGADD45B. أظهرت هذه الجينات قيمة تشخيصية عالية (AUC > 0.7) وأنماط تعبير تفاضلي متسقة عبر مجموعات البيانات.
الإثراء الوظيفي: تم إثراء المؤشرات الحيوية بشكل مشترك في مسار إشارات TNF-alpha عبر NF-kappaB. كما ربط تحليل GSEA هذه الجينات بمسارات p38MAPK، ومسارات MAPK المنشطة بالإجهاد، ونقاط تفتيش الدورة الخلوية.
البيئة المناعية الدقيقة: أظهرت مجموعة التهاب المفاث العظمي ارتفاعاً ملحوظاً في الخلايا البلازمية، والخلايا التائية التنظيمية (Tregs)، والبلعميات من نوع M0، والخلايا الصارية (Mast cells) في حالة الراحة. وُجدت ارتباطات محددة بين المؤشرات الحيوية وهذه الجمهرات المناعية (على سبيل المثال، ارتبط TLR7 إيجابياً بالبلعميات M0 والخلايا البلازمية).
رؤى الخلية الواحدة: أظهرت المؤشرات الحيوية أعلى تعبير في الخلايا الليفية تحت الغشاء الزلالي (SSF). كشف تحليل الزمن الكالتد (Pseudotime) لـ SSF عن 9 حالات تمايز، حيث أظهرت الحالة C1 درجة أعلى من التمايز. أظهر تعبير GADD45B وMYC اتجاهاً تصاعدياً عبر مراحل التمايز، بينما زاد TLR7 خطياً في الحالة 7.
الشبكة التنظيمية: تم إنشاء شبكة ceRNA تشمل 8 من lncRNAs و6 من miRNAs تنظم المؤشرات الحيوية الخمسة. وتم تحديد أدوية علاجية محتملة تستهدف هذه المؤشرات، بما في ذلك حمض البنزويك، والكرياتين، والساكساجليبتين.
التحقق: أكدت نتائج RT-qPCR أن MYC وPER1 وGADD45B انخفض تعبيرها بشكل كبير في عينات التهاب المفاث العظمي، بينما ارتفع تعبير CX3CR1 بشكل ملحوظ. أظهر TLR7 اتجاهاً نحو الارتفاع، وإن لم يكن ذلك ذا دلالة إحصائية في المجموعة السريرية الصغيرة (p=0.1854).
الأهمية والادعاءات يزعم المؤلفون أن هذه الدراسة توفر مورداً شاملاً لفهم "اليوبيكويتينوم" (Ubiquitinome) في التهاب المفاث العظمي. ومن خلال دمج بيانات التعبير الكلي والخلية الواحدة، تحدد الدراسة MYC، وTLR7، وCX3CR1، وPER1، وGADD45B كمؤشرات حيوية جديدة مرتبطة باليوبيكويتين والخاصة بالتهاب المفاث العظمي. وتشير الأبحاث إلى أن هذه الجينات تلعب أدواراً حاسمة في الاستجابات الالتهابية، وتدهور الغضروف، والتنظيم المناعي، واضطراب الإيقاع اليوماوي في سياق التهاب المفاث العظمي.
تفترض الدراسة أن هذه المؤشرات الحيوية تقدم أهدافاً علاجية محتملة، خاصة نظراً لارتباطها بمسار TNF-α/NF-κB. ويذكر المؤلفون صراحة أنه بينما توفر النتائج أساساً نظرياً للتشخيص والعلاج، فإن الدراسة لها قيود، بما في ذلك نقص التحقق السريري الواسع خارج تجارب RT-qPCR الأولية، وتدعو الأبحاث المستقبلية للتركيز على تجارب ميكانيكية أعمق وتطبيقات سريرية أوسع.