An anesthesia-related gene signature is associated with prognosis and immune phenotypes in colorectal cancer: a retrospective multi-cohort bioinformatics study
حددت هذه الدراسة المعلوماتية الحيوية الاستعادية متعددة المجموعات بصمة مرتبطة بالتخدير مكونة من 15 جينًا ومرتبطة بإنذار سرطان القولون والمستقيم والأنماط الظاهرية المناعية، رغم أن أداءها التنبؤي كان غير متسق عبر مجموعات التحقق الخارجية ويتطلب مزيدًا من التحقق التجريبي قبل التطبيق السريري.
في كل عام، يخضع الملايين من الناس لعمليات جراحية لاستئصال سرطان القولون والمستقيم، وهو مرض يصيب القولون والمستقيم. وبالنسبة لجميع هؤلاء المرضى تقريبًا، تتطلب العملية التخدير العام، وهو الدواء الذي يجعلهم يغطون في النوم ويمنعهم من الشعور بالألم. لعقود من الزمن، تساءل العلماء عما إذا كانت العقاقير المحددة المستخدمة لإحداث هذا النوم قد تفعل أكثر من مجرد إبقاء المريض ساكنًا. تشير بعض الأبحاث إلى أن بعض عوامل التخدير قد تؤثر على كيفية سلوك الخلا السرطانية، فربما تساعدها على النمو أو الانتشار، بينما قد تساعد أخرى جسم الإنسان على استخدام جهازه المناعي للمقاومة. ومع ذلك، ظل الرابط البيولوجي بين الأدوية المستخدمة في غرفة العمليات والنتائج طويلة المدى للسرطان لغزًا. إنه سؤال يهم المرضى والأطباء على حد سواء: هل يغير اختيار التخدير من فرص النجاة؟
وللتحقيق في ذلك، لجأ فريق من الباحثين إلى الأرشيفات الضخمة للبيانات الجينية المتاحة من آلاف المرضى. ركزوا على مجموعة محددة من الجينات المعروفة بتفاعلها مع أدوية التخدير. هذه الجينات، التي تعمل مثل المفاتيح التي تفتح أو تغلق عند التعرض للتخدير، تم تحديدها من خلال مقارنة كيفية استجابة أنسجة الدماغ لأدوية مختلفة في دراسات سابقة ومن خلال مراجعة الأدبيات العلمية الموجودة. ثم نظر الباحثون في الملفات الجينية لما يقرب من 650 مريضًا مصابًا بسرطان القولون والمستقيم لمعرفة أي من هذه الجينات الحساسة للتخدير كانت تتصرف بشكل مختلف في الأورام مقارنة بالأنسجة السليمة. ووجدوا أن 78 من هذه الجينات كانت بالفعل متغيرة في عينات السرطان، حيث أصبح بعضها أكثر نشاطًا بكثير والبعض الآخر أصبح أقل نشاطًا.
ومن قائمة الجينات المتغيرة هذه، بنى الباحثون نموذجًا رياضيًا للتنبؤ بكيفية مسار حالة المريض. اختاروا 1NS جينًا محددًا بدت الأكثر أهمية للبقاء على قيد الحياة ودمجوا مستويات نشاطها في درجة مخاطر واحدة. وعندما اختبروا هذه الدرجة على المجموعة الأولية من المرضى، نجحت بشكل معقول في الفصل بين أولئك الذين سيعيشون لفترة أطول وأولئك الذين لن يفعلوا. المرضى الذين لديهم درجة مخاطر عالية، مما يعني أن أورامهم أظهرت نمطًا محددًا من هذه الجينات الـ 15، كانت لديهم فرصة أكبر بكثير في عودة السرطان أو التسبب في الوفاة مقارنة بأصحاب الدرجة المنخفضة. كان النموذج فعالًا للغاية لدرجة أنه استطاع التنبؤ بنتائج البقاء على قيد الحياة حتى عندما وضع الباحثون في الاعتبار عوامل أخرى مثل حجم الورم أو مدى انتشاره إلى العقد اللمفاوية.
ثم حاول الباحثون معرفة ما إذا كان هذا النمط الجيني يخبرهم بأي شيء عن البيئة داخل الورم. وجدوا أن الأورام عالية المخاطر كانت مليئة بأنواع محددة من الخلايا المناعية التي عادة ما تثبط دفاعات الجسم الطبيعية، بدلًا من مهاجمة السرطان. كما أظهرت هذه الأورام علامات على كونها أكثر نشاطًا من الناحية الأيضية، حيث تحرق الطاقة وتبني خلايا جديدة بمعدل أسرع. استخدمت الدراسة أيضًا محاكاة حاسوبية لتخمين أي الأدوية قد تعمل بشكل أفضل للمرضى ذوي المخاطر العالية مقابل ذوي المخاطر المنخفضة، مما يشير إلى أن التركيبة الجينية للورم يمكن أن تؤثر على كيفية استجابته للعلاج الكيميائي. ومع ذلك، كان الباحثون حذرين في الإشارة إلى أن تنبؤات الأدوية هذه استندت إلى نماذج حاسوبية باستخدام خطوط خلوية، وليس على مرضى حقيقيين يتلقون العلاج.
وعندما اختبر الفريق نموذجهم على مجموعتين أخريين كبيرتين من المرضى من مستشفيات مختلفة، كانت النتائج مختلطة. في إحدى المجموعتين، أدى النموذج أداءً مشابهًا لما حققه في المجموعة الأولى، حيث نجح في التمييز بين المرضى ذوي المخاطر العالية والمنخفضة. أما في المجموعة الثانية، فقد فشل النموذج في التفوق على الصدفة، حيث اقتربت مقاييس الأداء من مستوى لا يعد فيه التنبؤ أفضل من رمي عملة معدنية. يشير هذا التباين إلى أنه بينما يلتقط التوقيع الجيني نمطًا بيولوجيًا حقيقيًا لدى بعض المرضى، فإنه ليس موثوقًا بما يكفي بعد ليُستخدم للجميع. وخلص الباحثون إلى أن نتائجهم هي نقطة انطلاق واعدة تربط الجينات المرتبطة بالتخدير بسلوك السرطان، لكن الأداة تحتاج إلى الكثير من الاختبار والتحسين قبل أن يمكن استخدامها لتوجيه القرارات الطبية. وقد أكدوا أن عملهم لا يثبت أن التخدير يسبب نمو السرطان أو انكماشه، بل إن الجينات المتأثرة بالتخدير هي جزء من قصة بيولوجية أكبر تحدد مدى عدوانية الورم.
ملخص تقني: البصمة الجينية المرتبطة بالتخدير في سرطان القولون والمستقيم
بيان المشكلة تشير الأدلة السريرية وما قبل السريرية الناشئة إلى أن عوامل التخدير (مثل البروبوفول مقابل السيفوفلوران) قد تؤثر بشكل متفاوت على تكاثر الأورام، والانتقال الظهاري-المتوسطي، والمناعة الجهازية المضادة للأورام. ومع ذلك، لا تزال المسارات الجزيئية المحددة التي تربط التخدير أثناء الجراحة بتطور سرطان القولون والمستقيم (CRC) غير واضحة. وبينما تم تحديد جينات فردية تستجيب للتحفيز بالتخدير، إلا أن هناك نقصًا في البصمات الشاملة القائمة على النسخ المتماثلة (transcriptome-based) والمستمدة من الجينات المرتبطة بالتخدير (ARGs) لتصنيف مرضى سرطان القولون والمستقيم. علاوة على ذلك، يتطلب البحث في المنطق البيولوجي للجينات المستهدفة عصبياً (مثل مستقبلات الناقلات العصبية) في السرطانات الظهارية تحققاً منهجياً، حيث تُعرف هذه المستقبلات بتعبيرها غير الطبيعي في سرطان القولون والمستقيم وقد تعدل سلوك الورم.
المنهجية استخدمت هذه الدراسة نهجاً معلوماتياً حيوياً استرجاعياً متعدد المجموعات باستخدام مجموعات بيانات متاحة علنًا ومجهولة الهوية.
بناء كتالوج الجينات: تم تجميع كتالوج للجينات المرتبطة بالتخدير (ARG) من مصدرين:
مراجعة منهجية للأدبيات المتعلقة بالجينات المرتبطة بعوامل التخدير.
تحليل مصفوفة استكشافي (GSE179004) يقارن بين عينات الورم الدبقي المعرضة للبروبوفول-ريميفينتانيل مقابل السيفوفلوران-ريميفينتانيل (n=3 لكل مجموعة). ونظراً لصغر حجم العينة، تم اختيار الجينات بناءً على قيمة اسمية لـ P<0.05 و ∣log2FC∣>0.5، مع الإقرار بأن هذه خطوة لتوليد الفرضيات وليست اكتشافاً قوياً إحصائياً.
جمع البيانات: تم استرداد البيانات النسخية والسريرية من مشروع "The Cancer Genome Atlas" (TCGA-CRC؛ 649 ورماً، 51 نسيجاً طبيعياً) لتدريب النموذج. واستُخدمت مجموعتان خارجيتان من بيانات GEO (GSE17536 و GSE39582) كمجموعات تحقق.
تحليل التعبير التفاضلي: باستخدام DESeq2، تم تحديد الجينات ذات التعبير التفاضلي (DEGs) بين الأورام والأنسجة الطبيعية في TCGA (بناءً على ∣log2FC∣>1، و P<0.05 المعدل). وأسفر تقاطع الـ DEGs مع كتالوج الـ ARG عن 78 جيناً مرتبطاً بالتخدير وذا تعبير تفاضلي (DEARGs).
تطوير البصمة: تم بناء بصمة إنذارية مكونة من 15 جيناً باستخدام فحص Cox أحادي المتغير يليه انحدار LASSO (L1-penalized) لـ Cox مع التحقق المتقاطع ذي الـ 10 طيات على مجموعة TCGA. وتم تصنيف المرضى إلى مجموعات عالية المخاطر ومنخفضة المخاطر بناءً على وسيط درجة المخاطر.
التحقق والتوصيف:
الأداء الإنذاري: تم التقييم عبر تحليل بقاء Kaplan–Meier، واختبارات log-rank، ومعامل التوافق (C-index) في TCGA والمجموعات الخارجية لـ GEO. كما تم إجراء انحدار Cox متعدد المتغيرات مع ضبط المتغيرات السريرية (العمر، الجنس، مرحلة TNM).
البيئة المناعية الدقيقة: تم تقييم الارتشاح المناعي للورم باستخدام خوارزميات CIBERSORT (22 نوعاً من الخلايا المناعية) و ESTIMATE. وتم تقدير الاستجابة لمثبطات نقاط التفتيش المناعية (ICI) باستخدام أطر عمل TIDE و Immunophenoscore (IPS).
التحليل الوظيفي: حدد تحليل إثراء مجموعة الجينات (GSEA) الاختلافات في المسارات. وتم التنبؤ بالحساسية للأدوية باستخدام خوارزمية oncoPredict (المدربة على بيانات GDSC2). كما تمت مطابقة التعبير على مستوى البروتين مع أطلس البروتين البشري (HPA).
النتائج الرئيسية
تحديد الجينات: تم تحديد 78 جيناً من الـ DEARGs، بما في ذلك جينات مرتفعة التعبير مثل GABRP و GABRA3 و DRD4، وجينات منخفضة التعبير مثل CHRNA3 و EPHX2. وسلط إثراء مسار KEGG الضوء على المسارات ذات الصلة بآليات التخدير (مثل المشبك GABAergic، وإدمان النيكوتين، والمشبك الدوباميني).
البصمة الإنذارية: تضمنت البصمة النهائية المكونة من 15 جيناً مزيجاً من الجينات الواقية (معاملات سالبة: CACNA1D، CCNE1، EPHX2، CHRNA1، PHKG1) والجينات المرتبطة بالمخاطر (معاملات موجبة: DRD4، GABRD، VEGFA، إلخ).
أداء TCGA: أظهرت البصمة قدرة تمييز متوسطة (C-index = 0.682). وكان لدى المرضى في المجموعة عالية المخاطر بقاء إجمالي أسوأ بشكل ملحوظ (HR = 2.91، P<0.001). وظلت درجة المخاطرة عاملاً إنذارياً مستقلاً بعد ضبط متغيرات TNM والمتغيرات السريرية الأخرى (Multivariate HR = 2.27، 95% CI 1.51–3.43).
التحقق الخارجي: كان الأداء متوسطاً في GSE17536 (C-index 0.668–0.681) ولكنه كان ضعيفاً في GSE39582 (C-index 0.535)، مما يقترب من التمييز العشوائي.
الارتباطات المناعية والبيولوجية:
ارتبطت الأورام عالية المخاطر بارتفاع ارتشاح الخلايا المثبطة للمناعة (الخلايا البلعمية من النوع M2، الخلايا التائية التنظيمية) وارتفاع درجات السداة (stromal scores).
أشارت درجات TIDE و IPS إلى ملفات مناعية أقل ملاءمة في مجموعة المخاطر العالية، رغم أن هذه تعتبر مؤشرات حسابية دون بيانات استجابة مناعية سريرية متطابقة.
كشف تحليل GSEA أن الأورام عالية المخاطر أظهرت زيادة في تنظيم المسارات الأيضية والتكاثرية (مثل أهداف Myc، الأكسدة الفوسفورية، دورة TCA)، بينما أظهرت الأورام منخفضة المخاطر إثراءً في برامج السداة والأوعية الدموية.
الحساسية للأدوية: اقترح التحليل الحاسوبي (In silico) قيم IC50 متفاوتة المتوقعة لمختلف العوامل الكيميائية بين مجموعات المخاطر، وإن كانت هذه مجرد فرضيات تتطلب تحققاً تجريبياً.
تأكيد البروتين: أكدت بيانات HPA وجود تعبير بروتيني قابل للكشف لـ 10 من أصل 15 جيناً من نموذج الدراسة في أنسجة سرطان القولون والمستقیم، حيث تم تصنيف CDKN2A كعلامة إنذارية غير مواتية في قواعد بيانات مستقلة.
الأهمية والادعاءات يضع المؤلفون هذه الدراسة كتحليل معلوماتي حيوي لتوليد الفرضيات وليس كأداة سريرية نهائية.
القيمة الإنذارية: تظهر البصمة الجينية (1.5-gene ARG) ارتباطاً بالإنذار والأنماط المناعية المستنتجة في مجموعة التطوير وإحدى المجموعات الخارجية، وتحتفظ بالأهمية بعد الضبط السريري.
الرؤية البيولوجية: تدعم الدراسة فرضية أن المسارات الجزيئية المرتبطة بالتخدير (التي تشمل القنوات الأيونية، ومستقبلات الناقلات العصبية، ومنظمات دورة الخلية) ذات صلة ببيولوجيا سرطان القولون والمستقيم، مما قد يعكس العواقب النسخية اللاحقة لشراسة الورم بدلاً من التنظيم السببي المباشر بواسطة عوامل التخدير.
القيود والتحذير: يذكر المؤلفون صراحة أن أداء البصمة في GSE39582 اقترب من التمييز العشوائي، مما يشير إلى محدودية القدرة على النقل (transportability). كما يؤكدون أن المؤشرات المناعية (TIDE/IPS) وتوقعات الحساسية للأدوية هي تأكيدات حسابية وليست تحققاً تجريبياً.
الخلاصة: النموذج ليس جاهزاً للاستخدام السريري. ويخلص المؤلفون إلى أن التحقق المستقل في مجموعات بيانات RNA-seq، والمعايرة المحددة مسبقاً، والاختبارات المختبرية الوظيفية هي أمور ضرورية قبل إمكانية تفسير البصمة سريرياً أو استخدامها لتوجيه استراتيجيات التخدير في التجارب المستقبلية.