From Transcriptomics to Game-Theoretic Models of Treatment Resistance in Glioblastoma Organoids
تقدم هذه الدراسة إطاراً لبيولوجيا الأنظمة القائم على البيانات، والذي يدمج تحليل النسخ المتماثل، والتعلم الآلي، والنمذجة القائمة على نظرية الألعاب للتنبؤ بالاستجابات التطورية للإشعاع في العضوانيات الورمية للأورام الدبقية، كاشفةً عن نقص الأكسجة كمحرك رئيسي للمقاومة قصيرة المدى، وموضحةً كيف يمكن لهذه النهج أن تُحسّن العلاجات السرطانية الشخصية.
المؤلفون الأصليون:Louise Spekking, Christer Lohk, Weronika Jung, Maikel Verduin, Sepinoud Azimi, Christopher Hubert, Marc Vooijs, Rachel Cavill, Kateřina Staňková
السرطان ليس عدوًا واحدًا ثابتًا؛ بل هو مجتمع متغير من الخلايا التي تتطور تحت الضغط. عندما يعالج الأطباء ورمًا بالإشعاع أو العلاج الكيميائي، فإنهم يسلطون قوة هائلة تقتل الخلايا الأكثر عرضة للتأثر. ومع ذلك، غالبًا ما يعمل هذا الضغط كمرشح (فلتر)، مما يسمح للخلايا القليلة التي تصادف كونها مقاومة بالبقاء والتكاثر. بمرور الوقت، ينمو الورم مجددًا، لكن هذه النسخة الجديدة تكون أصعب في القتل لأنها تتكون من هؤلاء الناجين المقاومين. وللتغلب على هذه العملية، يتجه العلماء إلى مجال يسمى "نظرية الألعاب التطورية". يعامل هذا النهج الورم ليس ككتلة من الأنسجة الميتة، بل كمجتمع من أنواع الخلايا المختلفة التي تتفاعل مع بعضها البعض. تمامًا كما تتنافس الحيوانات في الطبيعة على الموارد أو تتعاون من أجل البقاء، تتنافس خلايا السرطان وتتعاون داخل الورم. ومن خلال فهم هذه التفاعلات، يأمل الباحثون في تصميم علاجات تتوقع كيف سيتغير السرطان، بدلًا من مجرد الاستجابة له بعد أن يكون قد تكيف بالفعل.
لقد كان التحدي دائمًا هو معرفة ما يوجد بالضبط داخل الورم. للتنبؤ بكيفية تطور السرطان، يحتاج الأطباء إلى معرفة نسب أنواع الخلايا المختلفة الموجودة وكيف تتغير هذه النسب بمرور الوقت. تقليديًا، تطلب الحصول على هذا المستوى من التفاصيل إجراء تسلسل جيني أحادي الخلية مكلف ومعقد، وهو ما يحلل كل خلية فردية في العينة. وغالبًا ما يكون هذا مكلفًا للغاية ويستغرق وقتًا طويلاً للاستخدام الروتيني في العيادات. وبدلاً من ذلك، قرر فريق من الباحثين بقيادة لويز سبيكنيج وكريستر لوك معرفة ما إذا كان بإمكانهم الحصول على نفس المعلومات من البيانات "الكلية" (bulk data)، حيث يتم خلط المادة الوراثية لآلاف الخلايا وتحليلها دفعة واحدة. لقد ركزوا على الورم الأرومي الدبقي (glioblastoma)، وهو نوع عدواني من سرطان الدماغ، باستخدام نسخ صغيرة مزروعة مخبريًا من أورام المرضى تسمى "العضويات" (organoids). تعمل هذه العضويات كساحة اختبار واقعية، حيث تنمو في طبق لمحاكاة سلوك الورم الفعلي دون الحاجة إلى عمليات خزعة غازية من المريض.
بدأ الباحثون بتغذية البيانات الجينية لهذه العضويات في نوعين مختلفين من البرامج الحاسوبية لمعرفة أنواع الخلايا الموجودة. كان النوع الأول من البرامج غير موجه (unsupervised)، مما يعني أنه طُلب منه البحث عن الأنماط دون تعليمات مسبقة عما يجب العثหา. أما النوع الثاني فكان موجهًا (supervised)، مما يعني أنه تم تدريبه على خرائط موجودة لأنواع الخلايا تم إنشاؤها من دراسات أخرى. كان النهج غير الموجه يشبه محاولة فرز قطع أحجية مختلطة دون وجود صورة على الصندوق؛ فقد وجد ثلاث مجموعات متميزة من الخلايا. ومع ذلك، عندما فحص الباحثون الجينات التي تقود هذه المجموعات، أدركوا أن النتائج كانت مربكة. بدت إحدى المجموعات مثل الخلايا المناعية، وأخرى مثل خلايا الدماغ السليمة، والثالثة مثل الخلايا السرطانية. كانت هذه مشكلة لأن العضويات نمت في طبق بدون خلايا مناعية أو أنسجة دماغية سليمة، لذا فإن العثور عليها في البيانات يشير إلى أن الطريقة كانت تخلط بين الإشارات. أما النهج الموجه، الذي استخدم خريطة معروفة لأنواع الخلايا، فقد قدم قصة أكثر وضوحًا. فقد كشف أن العضويات كانت تتكون بالكامل تقريبًا من نوعين محددين من الخلايا السرطانية: نوع يزدهر في البيئات منخفضة الأكسجين، ونوع لا يزدهر فيها.
مع هذه الصورة الأكثر وضوحًا لمجتمعات الخلايا، انتقل الفريق إلى الخطوة التالية: نمذجة كيفية سلوك هذه الخلايا تحت ضغط الإشعاع. قاموا بمعالجة عضويات من ثلاثة مرضى بجرعتين مختلفتين من الإشعاع، 4 جيراي و10 جيراي، وراقبوا كيف يتغير التوازن بين الخلايا المحبة للأكسجين والخلايا الكارهة للأكسجين بمرور الوقت. ثم استخدموا إطارًا رياضيًا لوصف هذه التحولات، معتبرين التفاعل بين نوعي الخلايا بمثابة لعبة يحاول فيها كل نوع تعظيم بقائه الخاص. نجح النموذج في التنبؤ بالتركيبة المستقبلية للورم بناءً على التغييرات المبكرة الملحوظة بعد العلاج. وأظهرت النتائج أنه في معظم الحالات، كانت الخلايا التي تزدهر في الأكسجين المنخفض هي التي نجت وأصبحت في النهاية النوع السائد في الورم. وهذا أمر منطقي بيولوجيًا، حيث تحمي البيئات منخفضة الأكسجين الخلايا من الضرر الناتج عن الإشعاع. ومع ذلك، وجد الباحثون أيضًا أن النتيجة تعتمد بشكل كبير على المريض المحدد وجرعة الإشعاع. بالنسبة لأحد المرضى، أدت جرعة أقل من الإشعاع في الواقع إلى زيادة نسبة الخلايا الحساسة، مما يشير إلى أن تكلفة كون الخلية مقاومة قد تكون مرتفقة جدًا بالنسبة لها عندما يكون الضغط خفيفًا.
تسلط الدراسة الضوء على تمييز حاسم في كيفية فهمنا لمقاومة السرطان. وجد الباحثون أن التغيرات في مجتمعات الخلايا كانت مدفوعة بالبيئة المباشرة، وتحديدًا وجود الأكسجين، بدلاً من تحول الخلايا ببطء إلى أشكال جديدة عبر الطفرات. وهذا يشير إلى أن المقاومة يمكن أن تكون حالة مؤقتة، وهي مفتاح يُفتح ويُغلق اعتمادًا على الظروف، وليست تغييرًا جينيًا دائمًا. ومن خلال الجمع بين البيانات الجينية ونماذج نظرية الألعاب، أثبت الفريق طريقة لتتبع هذه التحولات دون الحاجة إلى تسلسل كل خلية على حدة. يقدم هذا النهج مسارًا عمليًا نحو الطب الشخصي. فإذا استطاع الأطباء استخدام هذه النماذج للتنبؤ بأنواع الخلايا التي ستسيطر على الورم بعد علاج معين، فيمكنهم تعديل العلاج في الوقت الفعلي. وبدلاً من إعطاء الجرعة القصوى التي تقتل كل شيء ولكنها تترك أقوى الناجين، يمكنهم استخدام استراتيجية تبقي الورم في حالة يسهل التحكم فيها، مما يمنعه من التطور إلى شكل لا يمكن علاجه. يظل هذا العمل بمثابة إثبات للمفهوم، يوضح أن هذه النماذج المعقدة يمكن بناؤها من البيانات القياسية، ولكنه يضع الأساس لمستقبل يكون فيه علاج السرطان حوارًا ديناميكيًا مع تطور الورم نفسه.
ملخص تقني: من النسخ المتجانس للترانسكريبتوم إلى النماذج القائمة على نظرية الألعاب لمقاومة العلاج في عضوانيات الورم الأرومي الدبقي
بيان المشكلة يمثل ظهور المقاومة المستحثة بالعلاج تحديًا رئيسيًا في علاج السرطانات المتقدمة، لا سيما الورم الأرومي الدبقي (Glioblastoma). وبينما أظهرت النماذج القائمة على نظرية الألعاب (العلاج التطوري أو التكيفي) نتائج واعدة في التجارب السريرية من خلال استباق التطور الورمي، إلا أن تطبيقها غالبًا ما يعوقه نقص البيانات الدقيقة والآنية حول تكوين الورم. وتعتمد الطرق الحالية التي تعتمد على تسلسل الحمض النووي الريبوزي أحادي الخلية (single-cell sequencing) على عمليات معقدة ومكلفة للمراقبة السريرية الروتينية. وبناءً على ذلك، هناك حاجة لاستخدام بيانات "الترانسكريبتوم" الكلية (bulk transcriptomic data) الأكثر سهولة في الوصول إليها لتقدير نسب أنواع خلايا السرطان المختلفة ومساراتها التطورية استجابةً للعلاج. تعالج هذه الدراسة الفجوة بين بيانات الترانسكريبتوم الكلية والنماذج الرياضية الميكانيكية لتطور الورم، مع التركيز بشكل خاص على العضوانيات (organoids) المشتقة من المرضى والمصابة بالورم الأرومي الدبقي.
المنهجية قام المؤلفون بتطوير إطار عمل لبيولوجيا الأنظمة القائم على البيانات، والذي يدمج تحليل الترانسكريبتوم، وتعلم الآلة، ونمذجة نظرية الألعاب. سارت المنهجية عبر ثلاث مراحل رئيسية:
تفكيك بيانات الترانسكريبتوم الكلية:
النهج غير الخاضع للإشراف (Unsupervised Approach): طبق المؤلفون خوارزميات تحليل المصفوفات غير السالبة (NMF) على بيانات RNA-seq الكلية من مجموعة بيانات TCGA-GBM. حددوا الرتبة المثلى للتحليل (والتي تم تحديدها كـ 3) واستخلصوا بصمات أنواع الخلايا دون افتراضات مسبقة حول ملفات تعريف الخلايا. وتم تحليل هذه البصمات عبر تحليل الإفراط في تمثيل مسارات KEGG لاستنتاج الهوية البيولوجية.
النهج الخاضع للإشراف (Supervised Approach): استخدم المؤلفون إطار InstaPrism، وهو إطار عمل لتفكيك البيانات يعتمد على المراجع ومحسن من حيث الكفاءة الحسابية. استخدمت هذه الطريقة بيانات RNA-seq أحادية الخلية من دراسة Neftel et al. كمرجع لتعيين أنواع الخلايا بناءً على الوحدات الميتالية (meta-modules) وهي: (MES-like 1, MES-like 2, AC, OPC, NPC1, NPC2). طُبق هذا النهج على بيانات المصفوفات الدقيقة (microarray) لعضوانيات الورم الأرومي الدبقي المشتقة من المرضى والتي عولجت بالإشعاع بجرعة 4 جراي و10 جراي.
التتبع الزمني:
تم إشعاع عضوانيات من ثلاثة مرضى، وجُمعت العينات قبل العلاج، ومباشرة بعد العلاج، وعند 24 و72 ساعة بعد العلاج.
تتبعت الدراسة التغيرات الزمنية في نسب مجموعات خلايا محددة (تحديدًا الخلايا الشبيهة بالميزنشيم المستقلة عن نقص الأكسجة مقابل الخلايا الشبيهة بالميزنشيم المعتمدة على نقص الأكسجة) والتي تم استنتاجها عبر طرق التفكيك.
نمذجة نظرية الألعاب:
بناءً على النتائج الخاضعة للإشراف، والتي أشارت إلى أن العضوانيات تتكون أساسًا من نوعين من الخلايا الميزنشيمية (MES-like 1: مستقلة عن نقص الأكسجة؛ MES-like 2: معتمدة على نقص الأكسجة)، قام المؤلفون بملاءمة نموذج ديناميكيات المكرر (replicator dynamics).
تم إنشاء مصفوفة لياقة (A) بنظام 2×2 حيث تمثل العناصر aij احتمالية التكاثر لنوع الخلية i عند تفاعلها مع النوع j. تم ضبط العناصر القطرية لتكون صفرًا لتبسيط الديناميكيات.
تمت ملاءمة معاملات النموذج مع بيانات السلاسل الزمنية المرصودة لاستنتاج المسارات التطورية وحساب الاستراتيجيات المستقرة تطوريًا (ESS).
النتائج الرئيسية
أداء النهج غير الخاضع للإشراف مقابل الخاضع للإشراف: حدد نهج NMF غير الخاضع للإشراف ثلاث مجموعات تقابل الخلايا المناعية، والخلايا السرطانية التكاثرية، والخلايا العصبية. ومع ذلك، أشار المؤلفون إلى أن هذه الطريقة واجهت صعوبة في عكس بيئة العضوانيات في المختبر بدقة (التي تفتقر إلى الخلايا المناعية والأوعية الدموية) وتطلبت تحديدًا تجريبيًا صعبًا لعدد أنواع الخلايا.
التفكيك الخاضع للإشراف: أثبت منهج InstaPrism الخاضع للإشراف أنه أكثر قابلية للتفسير والمتانة لهذا التطبيق المحدد. وكشف أن العضوانيات تتكون أساسًا (>85%) من خلايا شبيهة بالميزنشيم.
ديناميكيات الاستجابة للإشعاع:
المريض 1 و3: أظهرا اتجاهًا واضحًا حيث زادت نسبة الخلايا المعتمدة على نقص الأكسجة (MES-like 2) بعد الإشعاع، متجاوزة الخلايا المستقلة عن نقص الأكسجة (MES-like 1). وهذا يتوافق مع الفهم البيولوجي بأن الخلايا التي تعاني من نقص الأكسجة تكون أكثر مقاومة للإشعاع.
المريض 2: أظهر ديناميكيات تعتمد على الجرعة، حيث فضلت الجرعة المنخفضة (4 جراي) الخلايا المستقلة عن نقص الأكسجة، بينما فضلت الجرعة الأعلى (10 جراي) الخلايا المعتمدة على نقص الأكسجة، مما يشير إلى آليات مقاومة معقدة خاصة بكل مريض.
ملاءمة النموذج: نجح نموذج ديناميكيات المكرر في ملاءمة البيانات المرصودة مع متوسط مربع خطأ (MSE) ≤0.001. وتنبأت النماذج الملائمة بأنه في المدى الطويل (ESS)، ستشكل الخلايا المعتمدة على نقص الأكسجة الأغلبية عمومًا، رغم وجود كلا النوعين في حالة توازن مشترك.
تباين الاستراتيجية المستقرة تطوريًا (ESS): تباينت قيم ESS المحسوبة حسب المريض وجرعة الإشعاع، مما يشير إلى أن النتيجة التطورية ليست موحدة عبر جميع حالات الورم الأرومي الدبقي.
الأهمية والادعاءات يدعي البحث توفير منهجية لاستخدام بيانات الترانسكريبتوم الكلية لتقدير الديناميكيات التطورية للسرطان، مما يدعم تطوير علاجات تطورية مشخصنة.
دمج البيانات والنظرية: تثبت الدراسة أن بيانات الترانسكريبتوم من عضوانيات المرضى يمكن أن توفر بارامترات فعالة لنماذج نظرية الألعاب (تحديدًا ديناميكيات المكرر) لاستنتاج المسارات التطورية العابرة والتوازنات المستقرة.
المقارنة المنهجية: خلص المؤلفون إلى أنه بينما تعتبر الطرق غير الخاضعة للإشراف جذابة لعدم وجود افتراضات مسبقة، فإن النهج الخاضعة للإشراف (مثل InstaPran) توفر تقديرات أكثر دقة وقابلية للتفسير لديناميكيات التطور عندما تُعامل المقاومة العلاجية كسمة نوعية، لا سيما في سياق العضوانيات حيث تغيب أنواع خلايا معينة (مثل الخلايا المناعية).
دور نقص الأكسجة: تعزز النتائج أن نقص الأكسجة هو محرك رئيسي للاستجابات قصيرة المدى للإشعاع (من خلال اختيار الخلايا المعتمدة على نقص الأكسجة الأكثر مقاومة للإشعاع)، رغم أن تأثيره على النتوات التطورية طويلة المدى قد يتم تعديله بواسطة الاختلافات الجينية الخاصة بكل مريض وآليات المقاومة.
الفائدة المستقبلية: يذكر المؤلفون أن هذا الإطار يعمل كإثبات مفهوم لسلسلة عمليات (pipeline) حيث توفر بيانات الترانسكريبتوم معلومات للنماذج الميكانيكية. تشير الدراسة صراحة إلى أنه بينما تمت ملاءمة النماذج لبيانات العضوانيات، فإن تنبؤات الـ ESS لم يتم التحقق من صحتها بقياسات إضافية، حيث لم يكن هذا هو الهدف الأساسي للعمل الحالي. وبمجرد التحقق من صحتها باستخدام المزيد من البيانات الزمنية، يمكن استخدام هذه النماذج للتنبؤ بتأثير استراتيجيات العلاج البديلة (التوقيت والجرعات) لتأخير أو منع مقاومة العلاج. ويُبرز استخدام العضوانيات كجسر حاسم، مما يسمح بالتحقق من صحة تنبؤات هذه النماذج تجريبيًا قبل تطبيقها على المرضى الحقيقيين.
تظل الورقة متواضعة فيما يتعلق بنطاقها الحالي، حيث تشير إلى أنها لا تتحقق بعد من صحة تنبؤات الـ ESS باستخدام قياسات إضافية، كما أن الجوانب المكانية للتغاير الورمي (مثل نقص الأكسجة في المركز مقابل الحافة) لم تُدرج بعد في نموذج ديناميكيات المكرر الحالي الذي يفترض اختلاط العناصر (well-mixed).