Integrated analysis identifies prognostic biomarkers associated with mitochondrial dysfunction in head and neck squamous cell carcinoma
تحدد هذه الدراسة بصمة مرتبطة بخلل الوظائف الميتوكوندريا مكونة من 10 جينات في سرطان الخلايا الحرشفية في الرأس والعنق، والتي تعمل كعلامة بيولوجية مستقلة للتنبؤ بالإنذار، حيث تصنف المرضى بفعالية حسب مخاطر البقاء على قيد الحياة وتحدد السمات المناعية والجينية المتميزة المرتبطة بالاستجابة للعلاج المناعي.
يُعد سرطان الخلايا الحرشفية في الرأس والعنق مصطلحاً عاماً للسرطانات التي تنمو في الأغشية الرطبة للفم والحلق والحنجرة. وهو مرض يتسم بعدم القدرة على التنبؤ؛ فبعض المرضى يستجيبون جيداً للجراحة والإشعاع، بينما يواجه آخرون مستقبلاً قاتماً رغم العلاج المكثف. لعقود من الزمن، اعتمد الأطباء على الفحوصات البدنية والمسوحات لتخمين كيفية سلوك الورم، لكن هذه الأدوات غالباً ما تغفل القوى البيولوجية الخفية التي تقود المرض. إحدى تلك القوى الخفية هي الميتوكوندريا، وهي الهياكل الصغيرة داخل كل خلية والتي تعمل كمحطات طاقة، حيث تقوم بتحويل الغذاء إلى طاقة. وعندما تتعطل محطات الطاقة هذه، لا تفقد الخلايا الطاقة فحسب، بل يتغير سلوكها أيضاً، وغالباً ما تصبح أكثر عدوانية وتتعلم كيفية الاختباء من الجهاز المناعي للجسم. إن فهم كيف تؤثر محطات الطاقة المعطلة هذه على نمو السرطان قد يفتح طريقاً جديداً للتنبؤ بمن سيعيش ومن قد يستفيد من العلاجات المناعية الحديثة.
لقد شرع فريق من الباحثين في رسم هذه الخريطة للمناظر الطبيعية الخفية في سرطان الرأس والعنق من خلال النظر في التعليمات الجينية داخل الخلايا الورمية. بدأوا بمجموعة ضخمة من البيانات لآلاف المرضى، بحثاً عن الجينات المرتبطة بمشاكل الميتوكوندريا والمرتبطة أيضاً بمدة بقاء المريض على قيد الحياة. ومن خلال تمحيص هذه المعلومات، حددوا مجموعة محددة مكونة من ثمانية وتسعين جيناً ترتبط جميعها بخلل الميتوكوندريا وبالبقاء على قيد الحياة للمرضى. ومن هذه المجموعة الكبيرة، حصروا تركيزهم في قائمة دقيقة مكونة من عشرة جينات. تشكل هذه الجينات العشرة نوعاً من التوقيع البيولوجي الذي يمكنه تصنيف المرضى إلى مجموعتين متمايزتين: أولئ وجه منهم معرضون لخطر النتائج السيئة، والآخرون معرضون لخطر منخفض. ووجد الباحثون أن توقيع الجينات العشرة هذا كان متنبئاً موثوقاً بالبقاء على قيد الحياة، حيث أظهر أداءً جيداً باستمرار حتى عند اختباره على مجموعات مختلفة تماماً من المرضى من مستشفيات وقواعد بيانات أخرى.
كشفت الدراسة أن هاتين المجموعتين من المرضى لا يختلفان فقط في مدة بقائهما على قيد الحياة، بل في طبيعة أورامهما ذاتها. فقد أظهر المرضى في المجموعة منخفضة المخاطر، الذين كانت توقعاتهم أفضل، علامات على وجود جهاز مناعي قوي يحارب السرطان بنشاط. كانت أورامهم مليئة بالخلايا المناعية المعروفة بقدرتها على البحث عن السرطان وتدميره، كما أشارت بيئتهم الجينية إلى أنهم على الأرجح سيستجيبون جيداً للعلاج المناعي، وهو علاج يساعد الجهاز المناعي على التعرف على المرض. وفي المقابل، أظهرت المجموعة عالية المخاطر صورة مختلفة؛ حيث اتسمت أورامهم بوجود مناعي أضعف وأظهرت علامات على تغيرات جينية تؤدي غالباً إلى مقاومة العلاج. كما نظر الباحثون في الطفرات المحددة في الحمض النووي لهذه الأورام، ووجدوا أنه بينما كان الجين السرطاني الأكثر شيوعاً، TP53، متحوراً في كلتا المجموعتين، فإن النمط العام للأخطاء الجينية اختلف بشكل كبير بين المرضى عاليي ومخفضي المخاطر. وهذا يشير إلى أن توقيع الجينات العشرة يلتقط مزيجاً معقداً من العوامل الأيضية والمناعية التي تتجاوز مجرد الطفرات الجينية البسيطة.
ولضمان أن نتائجهم لم تكن مجرد تمرين حاسوبي، انتقل العلماء من تحليل البيانات إلى العمل المخبري. لقد اختاروا أحد الجينات العشرة، ويسمى PPP1R3A، لاختباره في خلايا سرطانية بشرية مزروعة في طبق. كان هذا الجين مثيراً للاهتمام لأنه، على عكس معظم الجينات الأخرى في قائمتهم التي كانت مفرطة النشاط، كان خامل النشاط في الأورام. نجح الباحثون في زيادة وتقليل مستويات هذا الجين في الخلايا لإثبات قدرتهم على التحكم فيه. وقد أكدت تجاربهم أنه يمكن التلاعب بالجين كما هو متوقع، مما أثبت الواقع البيولوجي وراء نموذجهم الحاسوبي. ورغم أن هذه التجربة المحددة لم تختبر ما إذا كان تغيير الجين سيوقف نمو السرطان، إلا أنها أثبتت أن الجين هو هدف حقيقي وملموس يمكن دراسته بشكل أكبر.
يشير العمل إلى أن صحة محطات الطاقة في الورم مرتبطة بعمق بكيفية رؤية الجهاز المناعي له ومهاجمته. فالمرضى الذين يعانون من خلل في الميتوكوندريا يؤدي إلى نمط جيني محدد، يبدو أن أورامهم تكون أكثر قدرة على الاختباء من الجهاز المناعي، مما يجعل علاجها أصعب. وعلى العكس من ذلك، يبدو أن أولئك الذين لديهم نمط مختلف يمتلكون أوراماً أكثر وضوحاً للدفاعات الجسدية. وبينما يحذر الباحثون من أن توقيع الجينات العشرة هذا ليس أداة معيارية للأطباء بعد ويتطلب المزيد من الاختبارات في التجارب السريرية، إلا أنه يقدم عدسة واعدة لفهم المرض. فمن خلال النظر في هذه الجينات العشرة، قد يتمكن الأطباء يوماً ما من إخبار المريض ليس فقط بمدى عدوانية سرطانه، بل أيضاً بما إذا كان ورم من المرجح أن يستجيب لأحدث وأقوى العلاجات المتاحة.
ملخص تقني: تحليل متكامل يحدد المؤشرات الحيوية التنبؤية المرتبطة بالخلل الوظيفي للميتوكوندريا في سرطان الخلايا الحرشفية في الرأس والعنق
بيان المشكلة يعد سرطان الخلايا الحرشفية في الرأس والعنق (HNSCC) من الأورام الخبيثة شديدة التغاير ذات الإنذار السيئ، لا سيالما لدى المرضى المصابين بأمراض مرتجعة أو نقيلية. وبينما يُعترف بالخلل الوظيفي للميتوكوندريا كمحرك لإعادة البرمجة الأيضية والتنظيم المناعي في السرطان، إلا أن قيمته التنبؤية المحددة وعلاقته بالبيئة المجهرية المناعية للورم (TME) والاستجابة للعلاج المناعي في سرطان (HNSCC) لا تزال غير واضحة. لقد ركزت الدراسات السابقة بشكل كبير على جينات ميتوكوندريا فردية أو مسارات محددة، مما أدى إلى نقص في التقييم المنهجي للاختلالات المنسقة للميتوكوندريا على المستوى النسخي.
المنهجية استخدمت الدراسة نهجاً يجمع بين المعلوماتية الحيوية والتحقق التجريبي:
الحصول على البيانات ومعالجتها: تم الحصول على البيانات النسخية والسريرية من مجموعة (TCGA-HSC) (522 عينة ورمية، 44 عينة طبيعية ضابطة) ومجموعة بيانات (GEO) المستقلة (GSE41613) التي تضم 97 عينة من سرطان الفم الحرشفي سلبي فيروس الورم الحليمي البشري (HPV-negative). تم استرداد الجينات المرتبطة بالخلل الوظيفي للميتوكوندريا (MDGs) من قاعدة بيانات (GeneCards) (8,957 جيناً).
فحص الجينات وبناء النموذج: تم تحديد الجينات المعبر عنها تفاضلياً (DEGs) بين الأنسجة الورمية والطبيعية. وتم تقاطع هذه الجينات مع الجينات التنبؤية (عبر تحليل رتبة لوك-رانج وانحدار كوكس) وجينات (MDGs) لتحديد 98 جيناً مرشحاً. تم تقسيم مجموعة (TCGA) إلى مجموعتي تدريب (70%) واختبار (30%). تم بناء نموذج تنبئي باستخدام انحدار (Cox) أحادي المتغير، وانحدار (LASSO)، والانحدار التدريجي، وانحدار (Cox) متعدد المتغيرات لاختيار البصمة الجينية المثلى.
حساب درجة المخاطر: تم حساب درجة المخاطر لكل مريض بناءً على مستويات التعبير للجينات المختارة ومعاملات الانحدار الخاصة بها. تم تصنيف المرضى إلى مجموعات عالية المخاطر ومنخفضة المخاطر بناءً على وسيط درجة المخاطر.
التحقق والتوصيف: تم التحقق من الأداء التنبئي للنموذج في مجموعة اختبار (TCGA)، ومجموعة (TCGA) الكاملة، ومجموعة (GSE41613) الخارجية باستخدام تحليل بقاء كابلان-ماير ومنحنيات (ROC) المعتمدة على الوقت.
التحليل البيولوجي والمناعي: استخدمت الدراسة خوارزميات (CIBRESORT)، و(ssGSEA)، و(ESTIMATE)، و(TIDE) لتحليل ارتشاح الخلايا المناعية، وخصائص البيئة المجهرية المناعية (TME)، والتعبير عن نقاط التفتيش المناعية، والاستجابة المحتملة للعلاج المناعي. كما تم تقييم المشاهد الطفرية الجسدية والعبء الطفري الورمي (TMB).
التحقق التجريبي: تم اختيار الجين PPP1R3A، وهو جين ممثل منخفض التعبير في البصمة الجينية، للتحقق منه مخبرياً. تم إنشاء نماذج الإفراط في التعبير والتعطيل في خطوط خلايا (CAL27) و(SCC-25)، مع تأكيد الفعالية عبر (qRT-PCR) واللطخة الغربية (Western blotting).
المساهمات والنتائج الرئيسية
تحديد بصمة مكونة من 10 جينات: حددت الدراسة بصمة تنبؤية مكونة من 10 جينات مرتبطة بالخلل الوظيفي للميتوكوندريا تشمل: DHRS2, PLAU, ULBP2, HOXB9, CELSR3, FST, TNFRSF4, RBP4, PPP1R3A, و IFI27. تم زيادة تعبير ثمانية جينات في الأورام، بينما انخفض تعبير RBP4 و PPP1R3A.
الأداء التنبئي: نجحت درجة المخاطر في تصنيف المرضى إلى مجموعات عالية ومنخفضة المخاطر مع اختلاف كبير في البقاء الإجمالي (OS) عبر جميع المجموعات (تدريب TCGA، اختبار TCGA، TCGA الكاملة، وGSE41613). أظهر النموذج دقة تنبؤية متوسطة إلى جيدة، حيث تراوحت قيم (AUC) للبقاء الإجمالي لـ 1 و3 و5 سنوات بين 0.651 و0.754. وأكد التحليل متعدد المتغيرات أن درجة المخاطر هي عامل تنبؤي مستقل.
الارتباط بالبيئة المجهرية المناعية: أظهر المرضى منخفضو المخاطر ملفاً مناعياً متميزاً يتميز بارتفاع ارتشاح الخلايا المناعية التكيفية (مثل خلايا CD8+ T، وخلايا CD4+ T الذاكرة النشطة، والخلايا البلازمية) وارتفاع درجات الوظائف المتعلقة بالمناعة (مثل النشاط السيتوليتي، وتحفيز الخلايا التائية المشترك، والاستجابة للإنترفيرون من النوع الثاني). وعلى العكس من ذلك، أظهرت المجموعات عالية المخاطر مستويات أعلى من خلايا CD4+ T الذاكرة المستقرة، والبلعميات من نوع M0/M1، والخلايا الصارية المنشطة.
التنبؤ بالاستجابة للعلاج المناعي: باستخدام خوارزمية (TIDE)، أظهرت المجموعة منخفضة المخاطر معدل استجابة متوقع أعلى لحصار نقاط التفتيش المناعية (36% مقابل 27% في المجموعة عالية المخاطر، P=0.028). وقد تعزز ذلك في مجموعة (IMvigor210) المستقلة للعلاج المناعي، حيث حقق المرضى منخفضو المخاطر نتائج بقاء أفضل وكانت درجات المخاطر المنخفضة مرتبطة باستجابات موضوعية (CR/PR).
التغيرات الجينومية: رغم أن الجينات العشرة المكونة للبصمة نفسها كانت ذات ترددات طفرية منخفضة، إلا أن مجموعات المخاطر أظهرت أنماطاً طفرية جسدية متميزة. فقد كانت المجموعة عالية المخاطر غنية بالطفرات في TP53, TTN, CDKN2A, FAT1, و NOTCH1، بينما أظهرت المجموعة منخفضة المخاطر توزيعات مختلفة تشمل PIK3CA و CSMD3.
التأكيد التجريبي: نجحت الدراسة في التحقق من تعديل تعبير الجين PPP1R3A في خطوط خلايا (HNSC)، مما أكد إمكانية التحكم في تعبير الجين بفعالية على مستويي الحمض النووي الريبي (mRNA) والبروتين، مما يدعم إمكانية استخدامه كهدف لمزيد من الاستقصاءات الوظيفية.
الأهمية والادعاءات يدعي المؤلفون أن هذه الدراسة تقدم تقييماً منهجياً للجينات المرتبطة بالخلل الوظيفي للميتوكوندريا في (HNSCC)، مما نتج عنه بصمة قوية مكونة من 10 جينات تدعم التقييم التنبؤي وتصنيف المخاطر. تُقدم هذه البصمة ليس فقط كمتنبئ للبقاء على قيد الحياة، بل أيضاً كأداة لتوصيف السمات المناعية والجينومية ذات الصلة بالاستجابة للعلاج المناعي. ومن خلال ربط الخلل الوظيفي للميتوكوندريا بالبيئة المجهرية المناعية للورم، تشير الدراسة إلى أن الحالات الأيضية للميتوكوندريا قد تؤثر على التهرب المناعي والحساسية للعلاج. إن التحقق التجريبي من PPP1R3A يوفر أساساً للدراسات الميكانيكية المستقبلية حول دور أيض الجليكوجين والوظيفة الميتوكوندريا في بيولوجيا (HNSCC).
يحافظ المؤلفون على نبرة متواضعة فيما يتعلق بنتائجهم، حيث يقرون بأن القدرة التنبؤية للنموذج متوسطة وأن الارتباطات المناعية تعتمد على التنبؤات الحسابية وتحليل المجموعات الاسترجاعية. ويؤكدون أن التحقق المستقبلي الاستباقي في مرضى (HNSCC) الذين يتلقون مثبطات نقاط التفتيش المناعية ومزيد من التجارب الميكانيكية ضروري لإثبات الفائدة السريرية والسببية البيولوجية للبصمة المحددة.