Omitting tumor proliferation from prognostic gene expression models: a preregistered cross-cancer test, a transcriptome-wide measurement, and a survey of how often the omission is made
تُظهر هذه الدراسة المسجلة مسبقاً أن إغفال تكاثر الورم من نماذج التعبير الجيني الإنذارية في سرطان الخلايا الكبدية وغيره من السرطانات يؤدي بشكل منهجي إلى تحيز جميع تقديرات الجينات تقريباً، ومع ذلك فإن توقعات المؤلف المحددة فيما يتعلق باتساق وحجم هذا التحيز عبر المجموعات المختلفة قد فُندت إلى حد كبير بواسطة البيانات.
في عالم أبحاث السرطان، غالبًا ما يبحث العلماء عن أنماط في التعليمات الجينية داخل الخلايا الورمية للتنبؤ بكيفية سلوك مرض المريض. إنهم يبحثون عن جينات محددة تعمل كعلامات تحذيرية، آملين في العثور على قائمة من المؤشرات التي يمكن أن تخبر الطبيب ما إذا كان الورم من المرجح أن ينمو بسرعة أو ينتشر. وللقيام بذلك، يبني الباحثون نماذج إحصائية تزن نشاط آلاف الجينات مقابل فترات بقاء المرضى على قيد الحياة. ومع ذلك، يجب أن تكون هذه النماذج حذرة في مراعاة عوامل أخرى تؤثر على البقاء، مثل عمر المريض أو حجم الورم. فإذا نسي الباحث تضمين عامل رئيسي يؤثر على كل من الجينات ووقت البقاء، فقد ينتج عن النموذج نتيجة مضللة؛ إذ قد يبدو جين معين هو السبب في النتيجة السيئة، بينما هو في الواقع مجرد راكب في حافلة يقودها قوة أخرى غير مقاسة. ومن أقوى هذه القوى هو "تضاعف الورم" (tumor proliferation)، وهو ببساطة مقياس لمدى سرعة انقسام وتكاثر الخلايا السرطانية. ولأن الخلايا سريعة الانقسام غالبًا ما تمتلك نشاطًا جينيًا مختلفًا عن تلك بطيئة النمو، فإن تجاهل هذه السرعة يمكن أن يشوه الصورة الكاملة لما تفعله الجينات حقًا.
لقد شرع باحث يدعى سون وو نام في استقصاء مدى تكرار هذا الخطأ المحدد وما هي التكلفة التي يترتب عليها عندما يحدث. بدأ بالنظر في خمسين ورقة علمية حديثة ركزت على سرطان الكبد. قرأها بعمق ليرى ما إذا كان المؤلف قد ضمن مقياسًا لسرعة انقسام الخلايا في نماذجه. كانت النتيجة صارخة: لم تكن هناك ورقة واحدة من بين واحد وأربعين ورقة استطاع قراءتها بالكامل قد عدلت لسرعة انقسام الخلايا. وفي الواقع، ذكرت أكثر من نصف تلك الدراسات أن مجموعاتها عالية الخطورة أظهرت علامات على انقسام سريع للخلايا، ومع ذلك تعاملت مع تلك السرعة كفضول بيولوجي بدلاً من كونها متغيرًا يجب التحكم فيه في حساباتهم الرياضية. وقد أشار هذا إلى أن عددًا هائلًا من الاستنتاقات المنشورة حول أي الجينات تتنبأ بالبقاء على قيد الحياة قد بُنيت على أساس هش، مفتقدة قطعة حاسمة من اللغز.
لفهم ما فُقد بالضبط بسبب ترك هذا العامل، لجأ نام إلى قاعدة بيانات عامة ضخمة تحتوي على بيانات جينية وبيانات بقاء لمئات المرضى المصابين بسرطان الكبد. أجرى اختبارًا شاملًا على كل جين تقريبًا كان نشطًا في تلك الأورام. ولكل جين، بنى نموذجين مختلفين للتنبؤ بالبقاء على قيد الحياة. النموذج الأول كان قياسيًا، يأخذ في الاعتبار العمر، والجنس، ومرحلة السرطان. أما النموذج الثاني فكان متطابقًا، باستثناء أنه تضمن أيضًا درجة تمثل مدى سرعة انقسام الأورام. ثم قارن النتائج. وعند إضافة سرعة انقسم الخلايا إلى المعادلة، تغيرت الأهمية المقدرة لكل جين تقريبًا. ولم يكن حجم هذا التغيير عشوائيًا؛ بل اعتمد مباشرة على مدى ارتباط نشاط ذلك الجين المحدد بسرعة الانقسام. فإذا كان الجين مرتبطًا بقوة بالنمو السريع، فإن تأثيره المتوقع على البقاء على قيد الحياة يتغير بشكل كبير بمجرد مراعاة سرعة النمو. وفي بيانات سرطان الكب، كان هذا التعديل قويًا لدرجة أنه تسبب في انهيار عدد الجينات التي بدت ذات دقة إحصائية من مئات الجينات إلى جين واحد فقط.
ثم طرح نام سؤالاً أعمق: هل يظل هذا النمط صحيحًا لجميع أنواع السرطان، أم أنه فريد لسرطان الكبد؟ وقبل النظر في البيانات، سجل تنبؤًا رسميًا لدى سجل علمي عام، موضحًا أنه يتوقع حدوث هذا الانخفاض في القيم إلى ما لا يقل عن تسعين بالمائة من أنواع السرطان المختلفة. اختبر هذا التنبؤ عبر عشر مجموعات بيانات مختلفة لأنواع السرطان، تتراوح من سرطان الرئة والكلى إلى سرطان القولون والمبيض. كانت النتيجة تفنيدًا لفرضيتة. لم يثبت النمط في كل مكان؛ ففي نصف أنواع السرطان التي اختبرها، كان اتجاه التغيير عكس ما توقعه. في بعض أنواع السرطان، جعلت إضافة درجة انقسام الخلايا تقديرات الجينات ترتفع، بينما في أنواع أخرى، جعلتها تنخفض. لم يكن اتجاه التغيير قاعدة بيولوجية ثابتة، بل كان محددًا بمدى قوة عمل انقسام الخلايا نفسه كمؤشر للبقاء على قيد الحياة في ذلك النوع المحدد من السرطان. ففي أنواع السرطان التي كان فيها الانقسام السريع علامة قوية على سوء الحالة، بدا الجينات المرتبطة بالانقسام أقل خطورة بمجرد مراعاة سرعة الانقسام. أما في أنواع السرطان التي لم يتنبأ فيها الانقسام بقوة بالنتيجة، فكان التغيير طفيفًا أو معكوسًا.
يكشف هذا الاكتشاف أن التحيز الناتج عن تجاهل تضاعف الورم ليس خطأً بسيطًا يدفع النتائج دائمًا في اتجاه واحد. بدلاً من ذلك، هو تشويه مرن يعتمد تمامًا على البيولوجيا المحددة للسرطان قيد الدراسة. كما سلطت الدراسة الضوء على مشكلة عملية في كيفية تقديم الأبحاث العلمية؛ فبسبب عدم ذكر معظم الأوراق العلمية العوامل التي عدلتها، لا يمكن للقارئ أن يعرف ما إذا كان الجين المذكور هو محرك حقيقي للمرض أم مجرد انعكاس لانقسام خلايا غير مقاس. ويقترح الباحث حلاً بسيطًا: يجب على العلماء الإبلاغ عن قوة إشارة انقسام الخلايا في مجموعتهم الخاصة من المرضى، وإظهار كيف تتغير تقديرات الجينات لديهم عند تضمين هذه الإشارة. لن يؤدي هذا بالضرورة إلى حل مشكلة معرفة أي الجينات هي المسببة فعليًا، ولكنه سيجعل التبعية مرئية، مما يسمح للأطباء والباحثين الآخرين برؤية السياق الكامل للأرقام التي يقرؤونها. ويعد هذا العمل تذكيرًا بأنه في المشهد المعقد لعلم وراثة السرطان، لا يمكن فصل الإجابة على سؤال حول جين واحد عن سلوك الورم ككل.
ملخص تقني: إغفال تكاثر الورم من نماذج التعبير الجيني الإنذارية
بيان المشكلة عادةً ما تقوم دراسات التعبير الجيني الإنذارية بفحص الجينات المرشحة مقابل نتائج البقاء على قيد الحياة، ومن ثم ملاءمة نماذج "كوكس" متعددة المتغيرات مع ضبط المتغيرات السريرية (العمر، الجنس، المرحلة). وتتمثل الفجوة الحرجة في هذه الممارسة التحليلية في الإغفال المتكرر لتكاثر الورم، وهو متغير يرتبط بقوة بكل من نتائج البقاء على قيد الحياة والتعبير عن العديد من الجينات. وعندما يكون المتغير مرتبطًا بكل من النتائج والمرشحات (الجينات المرشحة)، فإن إغفاله يؤدي إلى تحيز في نسب المخاطر (HRs) المقدرة لتلك المرشحات. وبينما ثبت أن التكاثر هو عامل إنذاري ويتم التقاطه عبر مجموعات جينية عشوائية، إلا أن حجم واتجاه التحيز الناتج عن إغفال التكاثر عبر مسح شامل للترانسكريبتوم، وما إذا كان هذا النمط يظل ثابتًا عبر أنواع السرطان المختلفة، لم يتم قياسهما بعد. علاوة على ذلك، ليس من الواضح مدى تكرار تعديل الباحثين لهذا المتغير في الدراسات المنشورة.
المنهجية تستخدم الدراسة نهجًا ثلاثي المحاور يجمع بين المسح الأدبي، والقياس الشامل للترانسكريبتوم في سرطان الخلايا الكبدية (HCC)، واختبار عبر أنواع السرطان تم تسجيله مسبقًا.
المسح الأدبي: تم سحب عينة من 50 ورقة بحثية حديثة حول التعبير الجيني الإنذاري في سرطان الخلايا الكبدية (HCC) من قاعدة بيانات PubMed. ومن بين 41 ورقة أمكن الوصول إلى نصوصها الكاملة، قام المؤلف بتصنيف ما إذا كان نموذج البقاء قد عدّل لمتغير التكاثر (باستخدام مصطلحات مثل Ki-67، أو MKI67، أو العد الميتوزي، أو درجات التكاثر) وما إذا كانت المتغيرات قد ذُكرت صراحة في النص.
القياس الشامل للترانسكريبتوم (الاكتشاف): في مجموعة (TCGA-HCC) (الاكتشاف) ومجموعة (GSE14520) المستقلة (التكرار)، قام المؤلف بملاءمة نموذجين من نماذج "كوكس" لكل جين مُعبر عنه:
النموذج (أ): معدل لـ (العمر، الجنس، والمرحلة).
النموذج (ب): النموذج (أ) بالإضافة إلى درجة تكاثر معيارية (متوسط ثمانية جينات مرتبطة بالتكاثر: MKI67, TOP2A, CCNB1, PCNA, BUB1, CCNA_2, AURKA, RRM2).
المقياس: عُرِّف "الانزياح" (shift) بأنه نسبة نسب المخاطر (HR) من النموذج (ب) إلى النموذج (أ). أما "المشهد" (landscape) فكان عبارة عن انحدار المربعات الصغرى العادية لـ (لوغاريتم الانزياح) على ارتباط "سبيرمان" الخاص بالجين مع درجة التكاثر.
الاختبار عبر أنواع السرطان المسجل مسبقًا: قبل تحليل البيانات، قام المؤلف بتسجيل فرضية على منصة (Open Science Framework - OSF). كان التنبؤ هو أنه في 90% على الأقل من أنواع السرطان المدرجة، سيكون ميل المشهد سالبًا، مع معامل تحديد (R2) وسيط ≥0.50. تم تطبيق الاختبار على مجموعات (TCGA) التي تستوفي معايير اشتمال محددة (عدد العينات ≥200، عدد الوفيات ≥60، بيانات كاملة). ميز التحليل بين المجموعات الأولية (المنفصلة) والمجموعات المجمعة لضمان الاستقلال الإحصائي.
التحليلات اللاحقة (Post Hoc):
عدم القابلية للانهيار (Non-collapsibility): تم تقييمها عبر ريجريشن (انحدار) لوغاريتم الانزياحات على لوغاريتم نسب المخاطر (HRs) للجينات ذات الارتباط القريب من الصفر بالتكاثر ("النطاق الصفري") للفصل بين عدم القابلية للانهيار الحسابي وبين الارتباك (confounding).
التعبير المفقود (Omitted Variable Expression): تم اختبار ما إذا كان تقريب خطي من الدرجة الأولى (حاصل ضرب معامل التكاثر والارتباط الجزئي للجين مع الدرجة) يمكنه تفسير الانزياحات الملحوظة.
مثال عملي: تم تحليل جين محدد (CYB5R3) لتوضيح الآلية واختبار المتانة ضد درجات تكوين الأنسجة وانتهاكات التناسب في المخاطر.
النتائج الرئيسية
المسح الأدبي: لم يقم أي من الأوراق الـ 41 التي تم مسحها بتعديل التكاثر (0%؛ فاصل ثقة 95% 0.0–8.6). ومع ذلك، ذكرت 56.1% منها نشاط مسار التكاثر كخاصية لمجموعات المخاطر العالية، ولم تذكر 68.3% قائمة المتغيرات صراحة في النص.
مشهد سرطان الخلايا الكبدية (HCC): في مجموعة الاكتشاف، كان ميل المشهد -0.278 (R2=0.785)، مما يشير إلى أن الجينات المرتبطة إيجابيًا بالتكاثر شهدت انخفاضًا في نسب المخاطر (HRs) عند إضافة التكاثر، بينما شهدت الجينات المرتبطة سلبيًا زيادة في نسب المخاطر. وقد تكرر هذا النمط في مجموعة (GSE14520) (الميل -0.180).
التأثير على المعنوية: في سرطان الخلايا الكبدية، انهار عدد الجينات التي نجت من التحكم في الخطأ التجريبي (q < 0.10) من 338 (بدون التكاثر) إلى جين واحد فقط (مع التكاثر).
دحض التنبؤ المسجل مسبقًا: تم دحض التنبؤ المسجل مسبقًا بأن الميل سيكون سالبًا في ≥90% من أنواع السرطان. فمن بين ثماني مجموعات أولية (منفصلة)، كان الميل سالبًا في 50% فقط (4/8). وكان الوسيط للاتفاق الزوجي للانزياحات عبر المجموعات سالبًا (-0.268)، مما يناقض التنبؤ باتساق الاتجاهية.
آلية الدحض: وُجد أن إشارة ميل المشهد تتحدد تمامًا بإشارة نسبة المخاطر الخاصة بدرجة التكاثر نفسها في تلك المجموعة المحددة. فإذا كان التكاثر ضارًا (HR > 1)، كان الميل سالبًا؛ وإذا كان التكاثر واقيًا أو معدوم الأثر (HR ≤ 1)، كان الميل موجبًا أو قريبًا من الصفر. وقد صمد هذا عبر 13 مجموعة (بما في ذلك المجموعات الاستكشافية).
التفكيك الكمي:
عدم القابلية للانهيار: ساهمت بنسبة صغيرة في الانزياح (زيادة تقارب 6.6% في HR لكل وحدة لوغاريتم HR) ولكنها لم تكن كافية لتفسير المشهد.
التعبير من الدرجة الأولى: فسر حوالي 84% من التباين في الانزياحات في سرطان الخلايا الكبدية، وهو ما جاء دون عتبة الـ 95% المسجلة مسبقًا. لقد بالغ التعبير بشكل منهجي في تقدير حجم الانزياح بنحو 10%.
البواقي (Residuals): لم يكن التباين غير المفسر ناتجًا فقط عن الجينات ذات الارتباط العالي مع الدرجة؛ إذ لم يؤدِ حذف هذه الجينات إلى استعادة القدرة التنبؤية للتعيير الخطي.
الأهمية والادعاءات يدعي البحث أن إغفال تكاثر الورم من النماذج الإنذارية يؤدي إلى تحيز منهجي وقابل للتنبؤ يحرك تقديرات كل جين تقريبًا. إن اتجاه هذا التحرك ليس خاصية ثابتة للترانسكريبتوم، بل يعتمد على القوة الإنذارية واتجاه التكاثر داخل المجموعة المحددة التي يتم تحليلها.
تتمثل المساهمة الأساسية للدراسة في إثبات أن افتراض "المقاس الواحد الذي يناسب الجميع" بشأن اتجاه تحيز الارتباك هو أمر خاطئ. لقد فشل التنبؤ المسجل مسبقًا لأن المؤلف عامل خاصية لوحظت في مجموعة الكبد (حيث التكاثر ضار بقوة) كقانون عالمي. وتخلص الورقة إلى ما يلي:
معايير الإبلاغ: يجب أن توضح التقارير الإنذارية ما إذا كان قد تم تعديل التكاثر أم لا.
الشفافية: يجب على المؤلفين الإبلاح عن نسبة المخاطر (HR) لدرجة التكاثر في مجموعتهم الخاصة، وتقديم تقديرات الجينات مع وبدون حد التكاثر، جنباً إلى جنب مع ارتباط الجين بالدرجة.
الغموض السببي: بينما تُظهر البيانات وجود انزياح منهجي، إلا أنها لا تستطيع التمييز بشكل قاطع ما إذا كان التكاثر يعمل كعامل مربك (confounder)، أو وسيط (mediator)، أو مُصادم (collider). ومع ذلك، فإن الادعاء الوصفي — وهو أن التقدير يعتمد بشكل منهجي على إدراج الحد — يظل قائمًا بغض النظر عن البنية السببية.
تؤكد الورقة أن هذه النتائج هي ارتباطية وخاصة بالممارسة التحليلية للمسح عالي الأبعاد. وهي لا تدعي تحديد علامات بيولوجية محددة، بل تهدف إلى قياس عواقب إغفال منهجي شائع. ويعمل هذا العمل كنقد لـ "الميتا-بحث" (meta-research)، داعيًا الباحثين إلى قياس السياق المحدد لبياناتهم بدلاً من الاعتماد على افتراضات مستمدة من أنواع أخرى من الأورام.