What is shared across 33 cancer types is the tissue, not the tumour: a non-malignant axis of the microenvironment
يكشف تحليل أكثر من 10,000 عينة ورمية من أرشيف (TCGA) أن هياكل الشبكات السببية القابلة للتكرار عبر 33 نوعاً من السرطان مدفوعة بمكونات البيئة الدقيقة للأنسجة غير الخبيثة المشتركة بدلاً من المحركات الخاصة بالأورام، مما يشير إلى أن الاستدلال السببي الشامل لأنواع السرطان من البيانات الضخمة يعكس في المقام الأول تكوين العينات بدلاً من التنظيم الجوهري للورم.
إن السرطان ليس مرضاً واحداً، بل هو مجموعة هائلة من مئات الحالات المختلفة، ينشأ كل منها من نوع محدد من الأنسجة في الجسم. وبينما تمثل الخلايا التي تشكل الورم الجزء الأكثر وضوحاً من المشكلة، إلا أن الورم لا يوجد في فراغ؛ فهو محاط ببيئة معقدة من الخلايا السليمة، بما في ذلك المدافعون المناعيون، وخلايا الدعم الهيكلي، والأوعية الدموية. لقد حاول العلماء طويلاً رسم خرائط للأسلاك الداخلية لهذه الأورام، بحثاً عن التعليمات المحددة التي تخبر الخلايا السرطانية بكيفية النمو والانتشار. وكان الأمل يكمن في أنه من خلال دراسة أنواع عديدة ومختلفة من السرطان في آن واحد، يمكن للباحثين العثوف على مجموعة عالمية من القواعد — برنامج تنظيمي مشترك — يحرك جميع هذه الأنواع. وإذا وجد مثل هذا البرنامج، فسيقدم هدفاً قوياً لعلاجات جديدة يمكن أن تعمل بشكل شامل، بدلاً من تطلب حل فريد لكل نوع من أنواع السرطان على حدة.
تتحدى دراسة جديدة هذا الأمل الذي طال أمدُه عبر النظر إلى البيانات بطريقة لم تُستخدم من قبل. فقد حلل الباحثون نشاط الجينات من أكثر من 10,000 عينة ورم تمثل 33 نوعاً مختلفاً من السرطان. واستخدموا طريقة مصممة لإيجاد علاقات السبب والنتيجة بين الجينات، متسائلين أي الجينات تقوم بتشغيل أو إيقاف جينات أخرى. ومع ذلك، اكتشفوا خللاً جوهرياً في كيفية إجراء هذه المقارنات عادةً. فبسبب استخدام أنواع السرطان المختلفة لمجموعات مختلفة من الجينات طبيعياً لأداء وظائفها، فإن مقارنتها مباشرة تشبه محاولة مقارنة خريطتين لمدينتين مختلفتين من خلال النظر إلى نفس الإحداثيات على الصفحة؛ فقد يكون الموقع هو نفسه، لكن إحدى الخريطتين تظهر حديقة بينما تظهر الأخرى طريقاً سريعاً. ولإصلاح ذلك، أجبر الباحثون جميع أنواع السرطان الـ 33 على أن تُحلل باستخدام نفس قائمة الجينات المكونة من 100 جين. وفقط عندما قاموا بمحاذاة البيانات بهذه الطريقة، استطاعوا رؤية ما هو مشترك حقاً عبر الأمراض المختلفة.
ما وجدوه كان مفاجئاً. فبدلاً من الكشف عن مجموعة مشتركة من التعليمات التي تحرك الخلايا السرطانية نفسها، أشارت البنية المشتركة بالكامل إلى الأنسجة السليمة المحيطة بالورم. وكان النمط الأكثر اتساقاً الذي وجدوه هو مجموعة من ثلاثة جينات تعمل كإشارات لاستدعاء الخلايا المناعية، وتحديداً نوع من خلايا الدم البيضاء التي تحارب العدوى. ظهرت هذه الجينات معاً في كل أنواع السرطان تقريباً، لكنها لم تكن جزءاً من الآلية الداخلية للسرطان. بدلاً من ذلك، كانت علامات على بيئة الورم، تشير إلى مقدار النشاط المناعي والأنسجة الهيكلية الموجودة في العينة. وأظهرت الدراسة أن الجينات الأكثر استخداماً لبناء خرائط السرطان هذه هي الجينات التي تتفاوت أكثر في نشاطها. وهذه الجينات شديدة التباين هي بالصدفة ذاتها التي تصف مزيج الخلايا السليمة في العينة، وليس الجينات التي تحرك السرطان. أما الجينات الفعلية المعروفة بأنها تسبب السرطان، وهي "المحركات"، فقد كانت مستقرة جداً في نشاطها لدرجة أنها لم تظهر حتى ضمن أفضل 100 جين تم اختيارها للتحليل.
كما كشف الباحثون عن قاعدة عملية حول مقدار البيانات المطلوبة للوثوق بهذه الخرائط. فقد وجدوا أنه عندما تحتوي الدراسة على عدد قليل جداً من العينات، يبدأ البرنامج الحاسوبي المستخدم لإيجاد الروابط في ابتكار أنماط من الضوضاء العشوائية. ومن خلال خلط البيانات لمعرفة ما سيحدث عندما لا يتبقى أي علم حقيقي، حددوا عتبة واضحة: تحتاج الدراسة إلى 172 عينة على الأقل لإيجاد روابط حقيقية بشكل موثوق عند البحث في 100 جين. وتحت هذا العدد، من المرجح أن تكون النتائج مجرد حوادث إحصائية. ويشير هذا الاكتشاف إلى أن العديد من الدراسات السابقة ذات أحجام العينات الصغيرة قد تكون أبلغت عن روابط زائفة. علاوة على ذلك، أكدت الدراسة أن الروابط القليلة التي ظهرت باستمرار عبر جميع أنواع السرطان الـ 33 لم تكن تتعلق بتواصل الخلايا السرطانية مع بعضها البعض، بل كانت تتعلق بالأنسجة السليمة التي تستجيب لوجود الورم.
إن تداعيات هذا العمل كبيرة فيما يتعلق بكيفية تفسير العلماء لبيانات السرطان. فهي تشير إلى أنه عندما يرى الباحثون نفس الروابط الجينية تظهر في أنواع عديدة من السرطان، فلا ينبغي لهم افتراض أنهم وجدوا آلية عالمية للسرطان. بل إنهم على الأرجح يشاهدون انعكاساً للأنسجة السليمة المشتركة التي تحيط بكل ورم. لا تقول الدراسة إن الخلايا السرطانية تفتقر إلى أسلاكها الفريدة الخاصة، بل تشير إلى أن الأدوات القياسية المستخدمة لإيجاد تلك الأسلاك مضبوطة حالياً لسماع ضجيج البيئة المحيطة بدلاً من ذلك. ولإيجاد المحركات الحقيقية للسرطان، سيتعين على العلماء النظر في نطاق أوسع بكثير من الجينات والتأكد من امتلاكهم لعدد كافٍ من العينات للفصل بين الإشارات الحقيقية والصدفة العشوائية. وحتى ذلك الحين، فإن الميزة الأكثر قابلية للتكرار والمشتركة عبر جميع أنواع السرطان ليست الورم نفسه، بل النسيج الذي يعيش فيه.
ملخص تقني: "ما يتشارك فيه 33 نوعًا من السرطان هو النسيج، وليس الورم: محور غير خبيث للبيئة الدقيقة"
بيان المشكلة تفترض الدراسات السببية الشاملة للسرطان (pan-cancer) روتينياً أن الحواف التنظيمية المستنتجة في نوع واحد من الأورام تنتقل إلى أنواع أخرى، مفسرةً الحواف المشتركة كدليل على برنامج تنظيمي مشترك. ومع ذلك، هناك عقبة منهجية أساسية تمنع الاختبار الصارم لهذا الافتراض عبر مجموعات بيانات "أطلس جينوم الورم" (TCGA): إذ تختار التحليلات القياسية لوحات جينية خاصة بكل نسيج (على سبيل المثال، أعلى 100 جين تباينًا لكل مجموعة). وبما أن هذه اللوحات تختلف باختلاف أنواع السرطان، فإن مصفوفات التجاور (adjacency matrices) تشير إلى أزواج جينية مختلفة عند نفس المواضع، مما يجعل المقارنة العنصرية (element-wise comparison) غير معرفة رياضياً. إن النهج القائم حالياً إما يحلل المجموعات بشكل منعزل (مما يحول دون المقارنة)، أو يدمج العينات (بافتراض التجانس الذي يناقض الواقع البيولوجي)، أو يقتصر على مقارنة المنظمات المشتركة دون فحص الحواف نفسها. تعالج هذه الورقة السؤال التالي: إذا تم تثبيت لوحة جينية مشتركة عبر 33 مجموعة من بيانات TCGA، فما هو الهيكل السببي العابر للمجموعات الذي يمكن استرداده، وهل يصف هذا الهيكل تنظيم الخلايا الورمية أم تكوين العينات؟
المنهجية تستخدم الدراسة إطار عمل لتفكيك سببي متعدد الدفعات ومكون من مرحلتين باستخدام خوارزمية NOTEARS (Zheng et al., 2018) للتحسين المستمر للرسوم البيانية الموجهة غير الحلقية (DAGs).
البيانات والمحاذاة: استخدم المؤلف بيانات تعبير جيني (R-SEM normalized) من نوع HiSeqV2 لـ 33 نوعاً من السرطان (N=10,459 عينة). ولتمكين المقارنة، قام ببناء "لوحة محاذية" ثابتة مكونة من 100 جين تم اختيارها بناءً على الانحراف المطلق للوسيط (MAD) عبر جميع المجموعات، بدلاً من استخدام قوائم الـ 100 جين الأعلى لكل مجموعة على حدة.
التقدير السببي: تم تطبيق رسم بياني NOTEARS لكل من الـ 33 نوعاً من السرطان مقابل هذه اللوحة المشتركة. يعمل النموذج على تقليل خطأ إعادة البناء باستخدام عقوبة L1 وقيد عدم حلقي قابل للتفاضل.
التفكيك: صاغ المؤلف مصفوفة التجاور السببية لكل سرطان k كـ Wk=W0+Δk، حيث تمثل W0 عموداً فقرياً سببياً مشتركاً، وتمثل Δk الانحرافات الخاصة بكل سرطان. تم تقدير W0 عبر الوسيط العنصري لمصفوفات التجاور لكل مجموعة، مع إسقاطها على فضاء الرسم البياني الموجه غير الحلقي (DAG manifold).
التشخيص والتحقق:
ضوابط العدم (Null Controls): تم خلط تسميات العينات داخل المجموعات لاختبار ما إذا كان التكرار ناتجاً عن حجم العينة أو خصائص المقدر.
ضبط التكوين (Composition Control): تم إجراء انحدار لدرجة مركبة من علامات غير خبيثة (سدوية، مناعية، بطانية) من جينات اللوحة لاختبار ما إذا كان الهيكل المتكرر يتم بوساطة تكوين الأنسجة الكلية.
التحقق الخارجي: تمت مقارنة الحواف المتكررة بقواعد بيانات STRING (التفاعلات الفيزيائية) وTRRUST (العوامل النسخية وأهدافها).
تحليل المحركات (Driver Analysis): تمت مقارنة مئوية تشتت الجينات المحركة الكلاسيكية مع عتبة الاختيار للوحة الـ 100 جين.
النتائج الرئيسية
المحاذاة شرط مسبق للمقارنة: في التصميم الأصلي غير المحاذي (أعلى 100 جين لكل مجموعة)، كانت مصفوفة الوسيط التجاورية هي صفرية تماماً (أكبر معامل وسيط بلغ 0.048). وعند المحاذاة للوحة مشتركة، ظهر هيكل غير صفري (أكبر معامل وسيط بلغ 0.731)، مما يثبت أن المقارنة عبر المجموعات لا تكون معرفة إلا بعد المحاذاة.
الهيكل المتكرر هو بيئة دقيقة، وليس ورماً خبيثاً: تكررت 14 زوجاً جينياً موجهاً فقط في 10 أو أكثر من أصل 33 مجموعة. لم تتضمن هذه الأزواج الجينات المحركة الكلاسيكية، بل شكلت وحدة متماسكة متمحورة حول ثلاثية (CXCL9, CXCL10, CXCL11) المحفزة بـ إنترفيرون غاما (γ-inducible) والخاصة بـ مستقبلات CXCL3. وخرطت الحواف المتكررة المتبقية إلى ستة محاور غير خبيثة: الكيموكينات المتعلقة بالخلاضات المتعادلة/CXCR2، والمحتوى الدموي (الهيموجلوبين)، والخلايا البلازمية/الغلوبولين المناعي، والسدوية/المتممة، وإشارات الخلايا النخاعية/السيتوكين.
المحور المهيمن يعكس تكوين النسيج: ارتبط المكون الرئيسي الأول (PC1) للمصفوفة المجمعة والمحاذية ارتباطاً قوياً (r=0.803) مع العلامات السدوية والمناعية الخارجية عبر جميع المجموعات الـ 33. كما أن الجينات المحركة الكلاسيكية، التي تقع في المئوية 50.8 من تشتت النسخ، كانت مستبعدة هيكلياً من لوحة الـ 100 جين (التي حُددت عند المئوية 0.487). ولم تدخل المحركات في اللوحة بشكل منهجي إلا عندما توسعت عرض اللوحة لتصل إلى k≳1,000.
الاعتماد على حجم العينة وعتبة التحديد: ارتبط عدد الحواف المستنتجة لكل مجموعة ارتباطاً عكسياً قوياً مع حجم العينة (r=−0.652) عبر المجموعة الكاملة، لكن هذا الارتباط تلاشى في المجموعات التي بلغ عدد عيناتها n≥200. وقد أدى خلط تسميات العينات إلى إلغاء المجموعة المتكررة (من 14 إلى 0 حافة) مع الحفاظ على الاعتماد على حجم العينة. حدد هذا "أرضية تحديد" تجريبية عند n∗≈1.7d (حيث d=100)؛ فالمجموعات التي تقل عن هذه العتبة تعيد حواف مدفوعة بالضجيج بدلاً من البيولوجيا.
الاتجاه والتكرار: بينما أظهرت الحواف المتكررة دلالة إحصائية مقابل نماذج العدم لتبديل الجينات، إلا أن اتجاهاتها لم تُحدد بواسطة البيانات الملاحظة الغاوسية (على سبيل المثال، وحدة CXCL9-CXCL10-CXCL11 هي مجموعة محفزة معاً، وليست تسلسلاً هرمياً). وقد استُعيدت ثلاثة من أصل 14 زوجاً في كلا الاتجاهين.
الأهمية والادعاءات تدعي الورقة مساهمات رئيسية ثلاث تعيد تعريف تفسير الشبكات السببية الشاملة للسرطان:
متطلب المحاذاة: المقارنة السببية عبر المجموعات غير معرفة رياضياً بدون لوحة سمات ثابتة. بدون المحاذاة، يكون "الهيكل المشترك" صفرياً بحكم التصميم.
التكوين فوق التنظيم: في اللوحات المختارة بناءً على التباين (أعلى 100 جين)، يصف الهيكل السببي المتكرر عبر أنواع السرطان تكوين النسيج (الارتشاح المناعي، المحتوى السدوي، المحتوى الدموي) بدلاً من البنية التنظيمية للخلايا الورمية. إن غياب الجينات المحركة هو نتيجة هيكلية لمعيار اختيار الجينات، وليس مصادفة بيولوجية.
أرضية حجم العينة التشغيلية: إن عدد الحواف المستنتجة بواسطة المقدرات السببية المتفرقة (sparse estimators) هو خاصية مرتبطة بحجم العينة للمقدر، وليس إشارة بيولوجية. تضع الدراسة قاعدة تجريبية (n≳lập1.7d) لتحديد ما إذا كانت المجموعة كبيرة بما يكفي للتمييز بين الإشارة والضجيج في تعلم الهيكل.
خلص المؤلف إلى أن الدراسات الشاملة للسرطان التي تبلغ عن حواف محفوظة يجب أن تفسرها كاعتمادات مشروطة بالبيئة الدقيقة. ولاستعادة الهياكل السببية على مستوى المحركات، يجب على الدراسات إما توسيع اللوحة الجينية بشكل كبير (مع مراعاة قيود حجم العينة) أو الشرط صراحةً على نقاء الورم. الهيكل المتكرر هو النسيج، وليس الورم.