Identification of diagnostic biomarkers and immune infiltration analysis in prostate cancer bone metastasis based on EMT-NM- related genes
تحدد هذه الدراسة S100A8 وPLOD2 كعلامات بيولوجية تشخيصية قوية لسرطان البروستاتا المنتشر في العظام من خلال دمج الجينات المتعلقة بالتحول الظهاري المتوسطي (EMT) واستقلاب النيوكليوتيدات، وبناء مخطط اسمي عالي الأداء، والكشف عن ارتباط معنوي بين تعبير S100A8 وتسلل الخلايا المثبطة المشتقة من النخاع في البيئة المناعية للورم.
المؤلفون الأصليون:Ke Chen, Jiayang Yue, Zhaogui Ding, Min Cheng, Ning Yao, Tiansheng Ye, Jia Sha, Zhanqiang Hua
يُعد سرطان البروستاتا مرضاً شائعاً بين الرجال، ولكن أخطر أشكاله يحدث عندما ينتشر السرطان من غدة البروستاتا إلى العظام. هذا الانتشار، المعروف باسم "النقيلة" (metastasis)، يصعب علاجه ويقلل بشكل كبير من عمر المريض. إن رحلة الخلية السرطانية من جزء إلى آخر في الجسم ليست عشوائية؛ بل تتضمن تغير شكل الخلية وسلوكها لتنفصل، وتنتقل عبر الدم، وتستقر في موقع جديد. لقد عرف العلماء منذ فترة طويلة أن هذه العملية تعتمد على عاملين رئيسيين: قدرة الخلية على تغيير شكلها الفيزيائي، وحاجتها إلى إعادة صياغة كيميائها الداخلية لتوفير الوقود للنمو السريع. وتحديداً، تسمح عملية تسمى "الانتقال الظهاري-المتوسطي" (epithelial-mesenchymal transition) للخلايا بأن تصبح متحركة، بينما يوفر التحول في كيفية معالجة الخلايا لبنات البناء الخاصة بالحمض النووي الريبوزي منقوص الأكسجين (DNA) والحمض النووي الريبوزي (RNA) الطاقة اللازمة للبقاء على قيد الحياة في بيئة جديدة. إن إيجاد طريقة للكشف عن هذه التغييرات مبكراً أمر بالغ الأهمية، حيث يمكن أن يسمح للأطباء بالتدخل قبل أن يتمكن السرطان من الاستقرار في الهيكل العظمي، إلا أن الأساليب الحالية غالباً ما تغفل عن علامات الإنذار المبكر هذه.
لقد شرع فريق من الباحثين في حل هذه المشكلة من خلال البحث عن بصمات جينية محددة تظهر عندما ينتقل سرطان البر prostatate إلى العظام. بدأوا بجمع البيانات من قواعد بيانات طبية عامة، والتي احتوت على معلومات جينية من مرضى يعانون من سرطان البروستاتا الموضعي ومن أولئك الذين انتشر سرطانهم بالفعل إلى العظام. كان هدفهم هو العثور على الجينات التي تنشط في كل من التحول الفيزيائي للخلايا وعملية التمثيل الغذائي المتغيرة لديها. ومن خلال مقارنة الملفات الجينية للمجموعتين، حددوا مجموعة صغيرة من الجينات التي برزت بشكل لافت. ومن بين مجموعة كبيرة من المرشحين، قلصوا القائمة إلى جينين فقط كانا أكثر نشاطاً باستمرار لدى المرضى المصابين بنقائل عظمية مقارنة بأولئك الذين يعانون من المرض المحصور في البروستاتا.
وُجد أن هذين الجينين، واللذين سُميا S100A8 وPLOD2، يعبران عن نفسهما بشكل مفرط بشكل ملحوظ في عينات النقائل العظمية. ولضمان موثوقية هذه النتائج، اختبرها الباحثون مقابل مجموعة منفصلة من بيانات المرضى، مؤكدين أن هذا النمط ظل ثابتاً. ثم قاموا ببناء نموذج رياضي، وهو في الأساس أداة تسجيل نقاط، تستخدم مستويات هذين الجينين للتنبؤ باحتمالية حدوث النقائل العظمية. وعندما اختبروا هذه الأداة، أثبتت دقة عالية في التمييز بين المرضى الذين لديهم انتشار والذين ليس لديهم، حيث تفوقت في أدائها على العديد من طرق التشخيص الموجودة. واقترح النموذج أن قياس هذين الجينين يمكن أن يعمل كنظام إنذار مبكر قوي للأطباء.
وبعيداً عن مجرد التشخيص، بحث الباحثون فيما تفعله هذه الجينات داخل الجسم. واكتشفوا أن وجود جين S100A8 مرتبط بقوة بتراكم نوع معين من الخلايا المناعية يسمى "الخلايا النخاعية المثبطة" (myeloid-derived suppressor cells). وتشتهر هذه الخلايا بقدرتها على تثبيط الدفاعات الطبيعية للجسم، مما يخلق أساساً لدرع يحمي السرطان من التعرض للهجوم من قبل الجهاز المناعي. ووجدت الدراسة ارتباطاً قوياً جداً بين المستويات العالية من S100A8 والعدد المرتفع لهذه الخلايا المناعية الواقية في بيئة العظام. وهذا يشير إلى أن الجين لا يكتفي فقط بتمييز وجود السرطان، بل قد يساعد بنشاط في إخفاء الورم عن الاستجابة المناعية للجسم، مما يسهل على السرطان ترسيخ وجوده في العظام.
أما الجين الآخر، PLOD2، فقد وُجد أنه يشارك في مسارات تساعد الخلايا على الالتصاق بمحيطها والتحرك، وهي خطوات أساسية لتمكين السرطان من غزو أنسجة جديدة. ومن خلال الجمع بين المعلومات من كلا الجينين، رسم الباحثون صورة شاملة لكيفية تفاعل السرطان مع بيئته. وأشاروا إلى أنه بينما تبدو نتائجهم واعدة، إلا أن الدراسة استندت إلى بيانات موجودة بالفعل، وسوف تحتاج إلى اختبار في بيئات سريرية أكبر وواقعية لتأكيد فائدتها لكل مريض. ومع ذلك، فإن تحديد هذين الجينين يقدم هدفاً ملموساً وجديداً للكشف المبكر، ويشير أيضاً إلى أن العلاجات المصممة لتعطيل نشاطهما أو تعطيل الخلايا المناعية التي يجذبانها يمكن أن توقف السرطان من الانتشار إلى العظام في المقام الأول.
ملخص تقني: تحديد المؤشرات الحيوية التشخيصية وتحليل الارتشاح المناعي في سرطان البروستاتا المنتشر في العظام بناءً على الجينات المرتبطة بالانتقال الظهاري الميزنشيمي (EMT) واستقلاب النوكليوتيدات (NM)
بيان المشكلة يُعد سرطان البروستاتا (PCa) من الأورام الخبيثة الرائدة لدى الرجال، ويُعتبر العظم الموقع الأكثر شيوعاً للانتشار البعيد. وبينما تتميز حالات سرطان البروستاتا في مراحله المبكرة بمعدل بقاء مرتفع، فإن تطور النقائل العظمية (BMPCa) يخفض معدل البقاء لمدة 5 سنوات بشكل حاد إلى حوالي 25%. تفتقر الممارسة السريرية الحالية إلى مؤشرات حيوية جزيئية عالية الحساسية والنوعية للتشخيص المبكر لسرطان البروستاتا المنتشر في العظام، مما يؤدي غالباً إلى تفويت نوافذ التدخل العلاجي. علاوة على ذلك، لا يزال التفاعل بين الانتقال الظهاري الميزنشيمي (EMT) — وهي عملية تدفع نحو غزو الورم — واستقلاب النوكليوتيدات (NM) — الذي يدعم التكاثر السريع — غير مستكشف بشكل كافٍ في سياق سرطان البروسات المنتشر في العظام. بالإضافة إلى ذلك، تتطلب ميزات البيئة الدقيقة المناعية المحددة المرتبطة بهذه العمليات الأيضية والانتقالية مزيدًا من التوضيح في حالات النقائل العظمية.
المنهجية استخدمت هذه الدراسة نهجاً قائماً على المعلوماتية الحيوية يدمج تحليل البيانات النسخية مع خوارزميات تعلم آلي متعددة لتحديد المؤشرات الحيوية التشخيصية.
مصادر البيانات: تم الحصول على الملفات النسخية من قاعدة بيانات GEO، باستخدام المجموعة التدريبية GSE32269 (22 عينة سرطان بروستاتا و29 عينة سرطان بروستاتا منتشر في العظام) ومجموعة التحقق GSE77930 (22 عينة سرطان بروستاتا و20 عينة سرطان بروستاتا منتشر في العظام).
اختيار الجينات: استرجع الباحثون 5,596 جيناً مرتبطاً بـ EMT و605 جيناً مرتبطاً بـ NM من قاعدة بيانات GeneCards.
تحليل التعبير التفاضلي: باستخدام حزمة "limma"، تم تحديد الجينات ذات التعبير التفاضلي (DEGs) بين مجموعتي سرطان البروستاتا وسرطان البروسات المنتشر في العظام (القيمة الاحتمالية المعدلة P < 0.05، و|log2FC| > 1.5).
فحص المرشحين: أدى تقاطع الجينات ذات التعبير التفاضلي (DEGs)، والجينات المرتبطة بـ EMT، والجينات المرتبطة بـ NM إلى الحصول على سبعة جينات مرشحة.
تكامل تعلم الآلة: تم تطبيق ثلاث خوارمازميات — وهي انحدار LASSO، واختيار الميزات Boruta، والغابة العشوائية (RF) — على المجموعة التدريبية لفرز ميزات تشخيصية قوية. تم تحديد جينات الميزات النهائية من خلال تقاطع نتائج هذه النماذج الثلاثة.
التحقق والنمذجة: تم التحقق من أنماط التعبير للجينات المختارة في كل من مجموعتي التدريب والتحقق باستخدام اختبار رتبة ويلكوكسون. كما تم إنشاء منحنيات خصائص تشغيل المستقبِل (ROC) لتقييم القيمة التشخيصية. تم بناء نموذج "نوموجرام" (Nomogram) بناءً على المؤشرات الحيوية النهائية وتم تقييمه باستخدام تحليل ROC، ومنحنيات المعايرة، وتحليل منحنى القرار (DCA).
التحليل الوظيفي والمناعي: تم إجراء تحليل إثراء مجموعة الجينات (GSEA) لاستك�شاف المسارات الإشارية المرتبطة. كما استُخدم تحليل إثراء مجموعة الجينات أحادي العينة (ssGSEA) لقياس ارتشاح الخلايا المناعية، يليه تحليل الارتباط بين المؤشرات الحيوية وتجمعات الخلايا المناعية.
النتائج الرئيسية
تحديد المؤشرات الحيوية: من أصل مجموعة أولية مكونة من سبعة جينات مرشحة، حدد تقاطع خوارزميات تعلم الآلة خمسة جينات ميزة هي: S100A8، وCENPF، وPLOD2، وTYMS، وKPNA2. وأكد تحليل التعبير الإضافي أن جينَي S100A8 وPLOD2 فقط كانا يعبران بشكل مفرط بشكل ملحوظ في عينات سرطان البروستات المنتشر في العظام عبر مجموعتي التدريب والتحقق.
الأداء التشخيصي: أظهر كل من S100A8 وPLOD2 إمكانات تشخيصية قوية. في المجموعة التدريبية، حقق S100A8 وPLOD2 قيم المساحة تحت المنحنى (AUC) بلغت 0.82 و0.83 على التوالي. وفي مجموعة التحقق، بلغت قيم AUC لـ S100A8 وPLOD2 نحو 0.76 و0.72.
بناء النوموجرام: تم بناء نوموجرام يدمج مستويات تعبير S100A8 وPLOD2. وقد أظهر أداءً تشخيصياً ممتازاً بمساحة تحت المنحنى (AUC) بلغت 0.89 في مجموعة التدريب. وأظهرت منحنيات المعايرة توافقاً جيداً بين النتائج المتوقعة والملاحظة (Hosmer–Lemeshes P > 0.05، ومتوسط الخطأ المطلق = 0.064)، كما أكد تحليل منحنى القرار (DCA) الفائدة السريرية للنموذج.
إثراء المسارات: كشف تحليل GSEA أن جين S100A8 كان غنياً بالمسارات المتعلقة بتكوين الريبوسومات، والأكسدة الفسفورية، والسمية الخلوية بوساطة الخلايا القاتلة الطبيعية. بينما كان جين PLOD2 غنياً بالمسارات التي تتضمن التفاعل بين مستقبلات المصفوفة خارج الخلية (ECM) والخلية، ودورة الخلية، والالتصاق البؤري.
الارتشاح المناعي: حدد تحليل الارتشاح المناعي 14 نوعاً من الخلايا المناعية المختلفة بشكل كبير بين مجموعتي سرطان البروستاتا وسرطان البروستات المنتشر في العظام. ولوحظ وجود ارتباط إيجابي قوي بين تعبير S100A8 وارتشاح الخلايا المثبطة المشتقة من النخاع (MDSCs) (معامل الارتباط r = 0.803).
الأهمية والادعاءات يزعم المؤلفون أن هذه الدراسة توفر أساساً نظرياً جديداً للتشخيص المبكر والعلاج الدقيق الموجه للبيئة الدقيقة لسرطان البروستات المنتشر في العظام. ومن خلال دمج مجموعات الجينات المرتبطة بـ EMT وNM، نجحت الدراسة في تحديد S100A8 وPLOD2 كمؤشرات حيوية تشخيصية قوية. ويُقدم النوموجرام الذي تم بناؤه كأداة قيمة للتنبؤ السريري، مما قد يكمل التشخيص المرضي التقليدي لتسهيل الفحص المبكر للمرضى المعرضين لمخاطر عالية.
علاوة على ذلك، فإن اكتشاف الارتباط القوي بين S100A8 وارتشاح الخلايا المثبطة المشتقة من النخاع (MDSCs) يقدم رؤى جديدة حول البيئة المناعية الدقيقة لسرطان البروستات المنتشر في العظام. ويشير المؤلفون إلى أن S100A8 قد يعمل كإشارة كيميائية تجتذب الخلايا المثبطة المشتقة من النخاع، مما يؤدي بالتالي إلى تثبيط وظيفة الخلايا التائية وتسهيل التهرب المناعي. يسلط هذا الاكتشاف الضوء على S100A8 كهدف محتمل للكشف المبكر والتدخل العلاجي الدقيق.
المحددات يقر المؤلفون بأن الدراسة محدودة بسبب صغر حجم العينة والتباين المحتمل في السكان، مما يستلزم تحققاً خارجياً صارماً من خلال دراسات استباقية متعددة المراكز وكبيرة النطاق. بالإضافة إلى ذلك، ونظراً لعدم وجود بيانات متابعة طويلة الأمد، لم يكن من الممكن تقييم القيمة التنبؤية المستقلة لهذه المؤشرات الحيوية للنتائج السريرية طويلة الأمد (مثل البقاء على قيد الحياة دون تقدم المرض).