Gastric cancer-specific neuroinflammation-related genes with prognostic significance: scRNA-seq with bulk RNA-seq data and experimental validation
تدمج هذه الدراسة بيانات تسلسل الحمض النووي الريبي (RNA) من النوع الضخم ومن الخلية الواحدة مع التحقق التجريبي لتحديد بصمة مرتبطة بالالتهاب العصبي مكونة من خمسة جينات (GRIA3، وNGF، وSLC17A7، وTGFB2، وTGFB3) تعمل كعلامة بيولوجية قوية للتنبؤ بالإنذار لسرطان المعدة، مما يربط هذه الجينات بإعادة تشكيل البيئة المجهرية للورم، والارتشاح المناعي، والنتائج السريرية.
يُعد سرطان المعدة مرضًا معقدًا حيث تتشابك أنظمة الدفاع في الجسم مع إشارات نمو الورم. وفي السنوات الأخيرة، أدرك العلماء أن البيئة المحيطة بالورم ليست مجرد خلفية سلبية، بل هي مشارك نشط في المرض. هذه البيئة، المعروفة باسم البيئة المجهرية للورم، هي حي مزدحم من الخلايا المناعية، وخلايا الدعم الهيكلي، والألياف العصبية. وبينما يُنظر إلى الالتهاب عادةً على أنه استجابة للإصابة أو العدوى، فإن نوعًا محددًا يسمى الالتهاب العصبي يتضمن تواصل النظام العصبي مع الجهاز المناعي، وأحيانًا بطرق تغذي السرطان بدلًا من محاربته. إن فهم كيفية حدوث هذا التداخل في الإشارات قد يكشف عن سبب بقاء بعض المرضى على قيد الحياة لفترات أطول، وقد يشير إلى طرق جديدة لعلاج المرض.
لقد شرع فريق من الباحثين في رسم خريطة لهذا الحوار الخفي في سرطان المعدة. بدأوا بجمع كمية هائلة من البيانات الجينية من آلاف المرضى، بحثًا عن جينات محددة تعمل كرسل بين الأعصاب والجهاز المناعي. ومن خلال مقارنة الملفات الجينية للأنسجة السرطانية بالأنسجة السليمة، واستخدام نماذج حاسوبية متقدمة لفرز الضجيج، حددوا مجموعة متميزة من الجينات المرتبطة بهذا النشاط الالتهابي العصبي. وجد الباحثون أن نشاط خمسة جينات فقط من هذه الجينات يمكن أن يروي قصة قوية عن مستقبل المريض. استُخدمت هذه الجينات، التي تشمل تعليمات لبروتينات تشارك في الإشارات العصبية والتنظيم المناعي، لبناء نموذج للمخاطر. وعندما طبق الفريق هذا النموذج على بيانات المرضى، نجح في تصنيف الأفراد إلى مجموعات عالية المخاطر ومنخفضة المخاطر. وكانت المجموعة ذات المخاطر العالية، والتي تميزت بأنماط محددة من هذه الجينات الخمسة، تعاني من فترات بقاء أقصر بكثير مقارنة بالمجموعة منخفضة المخاطر.
ولضمان أن هذه النتيجة لم تكن مجرد صدفة إحصائية، اختبر الفريق نموذجه على مجموعة منفصلة تمامًا من المرضى. وبينما أكدت النتائج قدرة النموذج العامة على التمييز بين مجموعات المخاطر، أشار الباحثون صراحةً إلى أن النموذج أظهر انخفاضًا في الفعالية التنبؤية في هذه المجموعة الخارجية للتحقق، مما يعني أن أداء تعميمه لم يكن مثاليًا. ورغم هذا القصور، ظل التوقيع الجيني مؤشرًا قيمًا للتنبؤ بسير المرض. ثم تعمق الباحثون لفهم ما يحدث داخل أورام المرضى ذوي المخاطر العالية. واكتشفوا أن هذه الأورام لم تكن ببساطة "باردة مناعيًا"، بل كانت تتميز باختراق واسع للخلايا المناعية، حيث أظهرت 16 نوعًا من الخلايا المناعية حضورًا أعلى بكثير في المجموعة عالية المخاطر. ومع ذلك، فإن هذا الاختراق الكثيف لم يترجم إلى قوة فعالة في محاربة السرطان؛ بل إن المشهد المناعي قد تغير بطريقة قد تعيق قدرة الجسم على محاربة السرطان، حيث ظهرت أنواع محددة مثل خلايا T المساعدة (CD4) النشطة والخلايا القاتلة الطبيعية أقل نشاطًا. كما أشارت الدراسة إلى أن هذه الأنماط الجينية يمكن أن تتنبأ بمدى استجابة المريض لأنواع معينة من العلاج الكيميائي، مما يلمح إلى مستقبل يمكن فيه تخصيص العلاج وفقًا للملف الجيني الخاص بكل مريض.
ولم يتوقف الاستقصاء عند المستوى الجيني. فمن أجل معرفة أي الخلايا تقوم فعليًا بتنفيذ هذه التعليمات، فحص الباحثون خلايا فردية من عينات الأورام باستخدام تقنية تقرأ الشفرة الجينية للخلايا المفردة. كشفت هذه الرؤية عالية الدقة أن نوعًا محددًا من الخلايا، وهو الخلايا الليفية (Fibroblast)، كان المحرك الرئيسي لهذه البيئة الالتهابية العصبية. الخلايا الليفية هي خلايا توفر الدعم الهيكلي في الحالات العادية، لكنها في هذه الأورام بدت وكأنها تعيد تشكيل الحي لصالح السرطان. وأظهرت الدراسة أن هذه الخلايا الليفية كانت تتواصل بكثافة مع الخلايا المناعية، وخاصة الخلايا التائية (T cells)، بطريقة قد تثبط الهجوم المناعي. علاوة على ذلك، تتبع الباحثون كيفية تغير هذه الخلايا الليفية بمرور الوقت، ووجدوا أن الجينات الخمسة الرئيسية تعمل وتتوقف في مراحل مختلفة من هذا التحول الخلوي، مما يشير إلى عملية ديناميكية لتطور الورم.
أخيرًا، أكد الفريق نتائجهم الحاسوبية في المختبر. فقد أخذوا عينات نسيجية فعلية من المرضى وقاسوا مستويات الجينات الرئيسية مباشرة. وطابقتت النتائج توقعاتهم: فالجينات المرتبطة بالمخاطر الأعلى كانت تُعبر عن نفسها بشكل مختلف في الأنسجة السرطانية مقارنة بالأنسجة السليمة. وبينما لم تحل الدراسة مشكلة سرطان المعدة، إلا أنها قدمت خريطة واضحة مدفوعة بالبيانات لكيفية تشكيل الالتهاب المرتبط بالأعصاب للمرض. ومن خلال تحديد هذه الجينات الخمسة المحددة والخلايا الليفية التي تحملها، قدم الباحثون إطارًا جديدًا لفهم سبب سلوك سرطان المعدة بهذا الشكل. ويشير هذا العمل إلى أن النظر في الحوار بين الأعصاب والجهاز المناعي قد يكون لا يقل أهمية عن النظر في الورم نفسه، مما يفتح بابًا جديدًا للتنبؤ بالنتائج وربما توجيه علاجات أكثر فعالية.
بيان المشكلة لا يزال سرطان المعدة (GC) سبباً رئيسياً للوفيات المرتبطة بالسرطان عالمياً، ويتميز بتباين شديد، وتشخيص في مراحل متأخرة، وانخفاض معدل البقاء على قيد الحياة لمدة 5 سنوات (OS)، لا سيما في المراحل المتقدمة. وبينما توجه المؤشرات الحيوية مثل PD-L1 وHER2 وCLDN18.2 العلاجات الحالية، لا تزال الحاجة قائمة لمؤشرات إنذارية موثوقة لتحسين تصنيف المخاطر والعلاج الفردي. تشير الأدلة الناشئة إلى أن الالتهاب العصبي — وهو التفاعل بين الإشارات العصبية والمناعية داخل البيئة المجهرية للورم (TME) — يلعب دوراً حاسماً في تطور سرطان المعدة، ومع ذلك، لا تزال القيمة الإنذارية المحددة للجينات المرتبطة بالالتهاب العصبي (NIRGs) في سرطان المعدة غير واضحة.
المنهجية استخدمت هذه الدراسة نهجاً تكاملياً متعدد الأوميات يجمع بين تسلسل الحمض النووي الريبي للكتلة (bulk RNA-seq) وتسلسل الحمض النووي الريبي أحادي الخلية (scRNA-seq) لتحديد والتحقق من السمات الإنذارية القائمة على الـ NIRGs.
مصادر البيانات: تم الحصول على بيانات bulk RNA-seq من مشروع (The Cancer Genome Atlas-Stomach Adenocarcinoma (TCGA-STAD) للتدريب، ومن (Gene Expression Omnibus (GSE62254 للتحقق الخارجي. واستُخدمت بيانات scRNA-seq (GSE183904 لتفكيك التباين الخلوي. كما تم استرجاع قائمة منسقة مكونة من 185 جيناً من الـ NIRGs من قاعدة بيانات التواقيع الجزيئية (MSigDB).
خط المعالجة المعلوماتية الحيوية:
التنميط الفرعي: تم تطبيق تجميع مصفوفة القيم غير السالبة (NMF) على ملفات التعبير الخاصة بالـ NIRGs لتحديد الأنماط الجزيئية الفرعية.
اختيار الجينات: تم تحديد الجينات المرشحة عن طريق تقاطع الجينات المعبر عنها تفاضلياً (DEGs) بين أنسجة سرطان المعدة والأنسجة الطبيعية، والـ DEGs بين الأنماط الفرعية لـ NMF، وقائمة الـ NIRG.
بناء النموذج: استُخدم تحليل انحدار كوكس أحادي المتغير، واختبار فرضية المخاطر النسبية، ونموذج غابة البقاء العشوائية (RSF) لاختيار أكثر الجينات الإنذارية أهمية وبناء نموذج درجة المخاطر.
التحقق والتحليل: تم التحقق من النموذج في مجموعة (GSE62254). وشملت التحليلات اللاحقة تقييم الارتشاح المناعي (ssGSEA)، والتنبؤ بحساسية الأدوية (IC50)، وإثراء المسارات (GSEA)، وبناء الشبكة التنظيمية (miRNA/lncRNA والتفاعلات بين المركب والهدف).
تحليل الخلية الواحدة: تمت معالجة بيانات scRNA-seq لتحديد أنواع الخلايا الرئيسية التي تعبر عن الجينات الإنذارية، متبوعاً بتحليل التواصل بين الخلايا (CellChat) وتحليل مسار الزمن الكاذب (Monocle).
التحقق التجريبي: تم إجراء تفاعل البوليميراز المتسلسل للنسخ العكسي الكمي (RT-qPCR) على عينات أنسجة سريرية (5 من حالات سرطان المعدة مقابل 5 طبيعية) للتحقق من تعبير الجينات المختارة.
المساهمات والنتائج الرئيسية
تحديد توقيع إنذاري مكون من خمسة جينات: حددت الدراسة خمسة جينات (NIRGs) مرتبطة بشكل كبير بإنذار المريض وهي: GRIA3 وNGF وSLC17A7 وTGFB2 وTGFB3. نجح نموذج المخاطر القائم على هذه الجينات في تصنيف المرضى إلى مجموعات عالية المخاطر (HRG) ومنخفضة المخاطر (LRG). في مجموعة TCGA-STAD، أظهرت جميع الجينات الخمسة نسب مخاطر (HR) أكبر من 1، مما يشير إلى أن التعبير العالي يرتبط بزيادة خطر الوفاة.
الأداء الإنذاري: أظهر نموذج المخاطر قدرة تنبؤية قوية، مع قيم المساحة تحت المنحنى (AUC) بلغت 0.911 و0.910 و0.851 لمدد بقاء 3 و5 و7 سنوات على التوالي في مجموعة TCGA. وحافظ النموذج على أداء ثابت في مجموعة التحقق المستقلة GSE62254.
الارتباطات السريرية والمناعية: كانت درجة المخاطر عاملاً إنذارياً مستقلاً، حيث ارتبطت بالسمات السريرية المرضية (العمر، المرحلة) والارتشاح المناعي. أظهرت المجموعة عالية المخاطر (HRG) ارتشاحاً أعلى لمختلف الخلايا المناعية (مثل الخلايا البائية المنشطة) ولكن مستويات أقل من خلايا T CD4+ المنشطة وخلايا NK CD56dim، مما يشير إلى بيئة مجهرية معقدة مع ضعف في المناعة الفعالة. كما تنبأ النموذج باختلاف حساسية الأدوية، حيث أظهر مرضى HRG قيم IC50 أعلى لعلاجات كيميائية قياسية (مثل 5-fluorouracil وpaclitaxel) وقيم IC50 أقل لبعض الوكلاء المستهدفين.
الرؤى الوظيفية: ربط تحليل إثراء مجموعة الجينات التوقيع بالحفاظ على القدرة الجذعية، والانتقال الظهاري الميزنشيمي (EMT)، وغزو الورم. ارتبط TGFB2 وTGFB3 بالمرض المتقدم وEMT، بينما ارتبط GRIA3 وSLC17A7 بمسارات التسرطن المبكرة.
الآلية الخلوية (scRNA-seq): حدد تحليل الخلية الواحدة الخلايا الليفية (fibroblasts) كنوع الخلايا الرئيسي الذي يقود إعادة تشكيل البيئة المجهرية للورم المرتبطة بالـ NIRG. أظهرت هذه الخلايا تعبيراً تفاضلياً كبيراً للجينات الإنذارية بين أنسجة سرطان المعدة والأنسجة الطبيعية. وكشف تحليل الزمن الكاذب عن أنماط تعبير ديناميكية لهذه الجينات أثناء تمايز الخلايا الليفية، حيث كان TGFB3 نشطاً في المراحل المبكرة وGRIA3/SLC17A7 مرتفعاً في المراحل اللاحقة. وسلط تحليل التواصل الخلوي الضوء على تعزيز إشارات MIF–(CD74+CXCR4) من الخلايا الليفية إلى الخلايا التائية في سرطان المعدة، مما يشير إلى آلية لإعادة التشكيل المناعي.
التحقق التجريبي: كشف تحليل RT-qPCR على مجموعة سريرية صغيرة (5 سرطان معدة مقابل 5 طبيعية) أن GRIA3 وNGF وSLC17A7 انخفض تعبيرها بشكل ملحوظ في أنسجة سرطان المعدة مقارنة بالأنسجة الطبيعية. وتتعارض هذه النتيجة مع الارتفاع الملاحظ في مجموعة بيانات TCGA-STAD، وهو تباين عزاه المؤلفون إلى صغر حجم عينة التحقق، والتباين الفردي، والاختلافات في خصائص العينات (المرحلة، علم الأمراض) بين المجموعات.
الأهمية والادعاءات يزعم المؤلفون أن هذه الدراسة توفر إطاراً إنذارياً جديداً لسرطان المعدة من خلال دمج تواقيع الالتهاب العصبي مع بيانات الأوميات المتعددة بشكل منهجي. يوفر نموذج المخاطر المقترح أداة محتملة لـ:
التصنيف: تحسين التنبؤ بالبقاء على قيد الحياة وتصنيف المخاطر بما يتجاوز العوامل السريرية المرضية التقليدية.
الفهم الآلي: توضيح دور "النيش العصبي-المناعي" في سرطان المعدة، وبالتحديد تسليط الضوء على دور الخلايا الليفية المرتبطة بالسرطان (CAFs) في توسط إشارات الالتهاب العصبي، وكبح المناعة، وتطور الورم.
التوجيه العلاجي: تقديم مراجع لاختيار العلاج الفردي بناءً على تنبؤات حساسية الأدوية وتحديد التفاعلات المحتملة بين المركب والهدف لتطوير الأدوية في المستقبل.
تخلص الدراسة إلى أنه بينما يبدو النموذج واعداً، فإن أداءه في التعميم يتطلب مزيداً من التحقق في مجموعات أكبر، كما تستوجب الآليات التنظيمية للجينات الموقعة مزيداً من الاستقصاء العميق. يهدف العمل إلى وضع أساس نظري للعلاجات التي تستهدف كلاً من الإشارات العصبية ونقاط التفتيش المناعية في سرطان المعدة.