تحديد وتوصيف وظيفي للجينات ذات التعبير المشترك في الموت الخلوي النحاسي والموت الخلوي الحديدي
تم دمج أربع مجموعات بيانات تدريبية تضم 421 عينة من التهاب القولون التقرحي (UC) و97 عينة ضابطة سليمة. وقد تم تحديد إجمالي 551 جينًا متباين التعبير، بما في ذلك 362 جينًا مرتفع التعبير و189 جينًا منخفض التعبير. وأسفر تحليل الارتباط عن 444 جينًا مرتبطًا بـ ferroptosis وcuproptosis، وأدى التقاطع مع الجينات متباينة التعبير إلى إنتاج 32 جينًا متباين التعبير مشترك التعبير بين cuproptosis وferroptosis (CF-DEGs).
أظهرت الجينات 32 CF-DEGs إثراءً في الاستجابة للجروح، وأيونات النحاس، ونقل الأحماض الدهنية، والليبوبوليساكاريد، ومرض التهاب الأمعاء، ومسارات إشارات NF-κB، ومسارات إشارات TNF، والموت الخلوي الحديدي (الجدول التكميلي 4).
تحديد أولويات المؤشرات الحيوية الأساسية من خلال التعلم الآلي المتكامل وتحليل شبكة共التعبير الجيني المترابط (WGCNA)
حدّد تحليل شبكة التعبير الجيني المشترك الموزون (WGCNA) الوحدات MEpurple وMEbrown وMEblack كأكثر الوحدات ارتباطاً بقوة بـ UC. واحتوت الوحدات المختارة مجتمعة على 1,426 جيناً (الجدول التكميلي 5).
اختار كل من عامل التقلص والاختيار الأقل مطلقاً (LASSO)، وآلة متجه الدعم مع استبعاد الميزات المتكرر (SVM-RFE)، والغابة العشوائية 21 و32 و19 ميزة على التوالي. وهناك ثلاثة عشر جيناً مشتركة بين النماذج الثلاثة جميعها. وقد أدى دمج توافق التعلم الآلي، ونواة التفاعل بين البروتينات المصنفة بواسطة CytoHubba، ووحدات WGCNA المختارة إلى الحصول على خمسة مؤشرات حيوية مرشحة مرتبطة بـ UC وهي: LCN2 وIDO1 وCXCL2 وNOS2 وCD274.
التقييم التشخيصي وإثراء المسارات لبصمة المؤشرات الحيوية
ارتبطت جميع المؤشرات الحيوية المرشحة الخمسة ارتباطاً إيجابياً وظهر ارتفاع في تنظيمها في حالات التهاب القولون التقرحي (UC) ضمن مجموعة التدريب (الشكل 2A، B). حقق كل مؤشر مساحة تحت منحنى خصائص التشغيل للمستقبل > 0.80 في مجموعة التدريب و > 0.75 في مجموعة التحقق المستقلة GSE47908 (الشكل 2C–E).
أظهر المخطط البياني (nomogram) المكون من خمس جينات معايرةً جيدة، وذلك وفقاً لاختبار Hosmer–Lemeshow P > 0.05 ومتوسط الخطأ المطلق < 0.1، والقدرة التمييزية في المجموعة التي خضعت للتحليل، حيث بلغت المساحة تحت المنحنى ومؤشر التطابق 0.945، مع فاصل ثقة 95% لمؤشر التطابق يتراوح بين 0.922 و0.968 (الشكل 3).
ربط تحليل إثراء مجموعة الجينات المرشحين الخمسة بمسارات تأشير JAK–STAT، وinterferon–RIPK1/3، وToll-like receptor–NF-κB (الشكل 4).

الشكل 2: تعبير المؤشرات الحيوية المرشحة وأداؤها التشخيصي. (A) خريطة حرارية للارتباط. (B, C) تعبير المؤشرات الحيوية ومنحنيات خصائص تشغيل المستقبل في مجموعة التدريب. (D, E) تعبير المؤشرات الحيوية ومنحنيات خصائص تشغيل المستقبل في مجموعة التحقق GSE47908. ROC: خصائص تشغيل المستقبل؛ AUC: المساحة تحت المنحنى؛ UC: التهاب القولون التقرحي. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 3: تطوير وتقييم نوموغرام UC. (أ) نوموغرام لخمسة جينات. (ب) مخطط المعايرة. (ج) منحنى خصائص تشغيل المستقبل. ROC: خصائص تشغيل المستقبل؛ UC: التهاب القولون التقرحي؛ AUC: المساحة تحت المنحنى؛ C-index: مؤشر التوافق. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 4تحليل إثراء مجموعة الجينات للمؤشرات الحيوية الخمسة المرشحة. (أ-هـنتائج تحليل إثراء مجموعة الجينات لكل من LCN2 وIDO1 وCXCL2 وNOS2 وCD274، على التوالي. GSEA: تحليل إثراء المجموعات الجينية. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
تحليل مشهد البيئة الدقيقة المناعية وشبكة التنظيم
أظهرت عينات UC زيادة في نسب الخلايا المتعادلة، وخلايا CD4⁺ T الذاكرة المنشطة، والبلعميات M1، والخلايا الصارية المنشطة، مع انخفاضات متبادلة في البلعميات M2، والخلايا الصارية الساكنة، والخلايا التغصنية الساكنة. وقد تكررت هذه الأنماط في مجموعة التحقق (الشكل 5A–E).
ارتبط المرشحون الخمسة إيجابياً بالخلايا المتعادلة، وخلايا CD4⁺ T الذاكرة المنشطة، وبالخلايا البلعمية من النوع M1، وارتبطوا عكسياً بالخلايا الصارية الساكنة والخلايا البلعمية من النوع M2 (الشكل 5F–J). احتوت شبكة الجينات وmiRNA على 98 عقدة و146 حافة. أظهر كل من hsa-miR-34a-5p وhsa-miR-16-5p أعلى اتصال للمؤشرات الحيوية، بينما كان AR وRELA أكثر عوامل النسخ اتصالاً (الشكل 6).

الشكل 5: تحليل الارتشاح المناعي. (A) تكوين الخلايا المناعية. (B) الفروق بين المجموعات في نسب الخلايا المناعية. (C) الخريطة الحرارية لارتباط الخلايا المناعية. (D, E) ارتشاح الخلايا المناعية في مجموعتي التدريب والتحقق. (F–J) الارتباطات بين المؤشرات الحيوية المرشحة ومجموعات الخلايا المناعية. UC: التهاب القولون التقرحي. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 6: الشبكات المتوقعة لـ miRNA وعوامل النسخ. (A) شبكة الجين–miRNA. تشير الدوائر إلى المؤشرات الحيوية المرشحة، وتشير المربعات إلى miRNAs. (B) شبكة عامل النسخ–الجين. تشير المعينات إلى المؤشرات الحيوية المرشحة، وتشير المثلثات المقلوبة إلى عوامل النسخ. miRNA: microRNA; TF: transcription factor. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
ديناميكيات التعبير الزمكاني بدقة الخلية الواحدة
بعد مراقبة الجودة، شكل 31,712 خلية 22 عنقوداً تم تصنيفها كـ تسع مجموعات خلوية رئيسية (الشكل 7A–C). وكانت نسب الخلايا الظهارية وخلايا البلازما أعلى في حالات UC. كما أظهرت الخلايا الظهارية إثراءً في LCN2 و NOS2، بينما أظهرت الخلايا النخاعية إثراءً في IDO1 و CXCL2 و CD274 (الشكل 7D–G).
حدد التصنيف الفرعي للخلايا الظهارية 11 مجموعة فرعية، مع توسع في الخلايا القولونية الالتهابية وإثراء في LCN2 وNOS2 في هذه المجموعة الفرعية (الشكل 8). وحدد التصنيف الفرعي للخلايا النخاعية سبع مجموعات فرعية، مع زيادة في الخلايا الوحيدة، وانخفاض في البلاعم، وإثراء في IDO1 وCXCL2 وCD274 في الخلايا الوحيدة (الشكل 9).
تراكمت الخلايا القولونية الالتهابية في وقت متأخر من المسار الظهاري، مع زيادة في تعبير LCN2 وNOS2 (الشكل 10). وأظهرت الخلايا الوحيدة مساراً متميزاً مرتبطاً بالتهاب القولون التقرحي (UC)، مع تعبير ديناميكي لكل من IDO1 وCXCL2 وCD274 (الشكل 11).
أظهرت الخلايا القولونية الالتهابية أقوى إشارات صادرة واتصالاً بارزاً مع الخلايا الوحيدة. وكان الزوج APP–CD74 من أهم أزواج الربيطة-المستقبل بين هذه المجموعات الفرعية (الشكل 12).

الشكل 7: التعبير عن المؤشرات الحيوية المرشحة عبر المجموعات الخلوية. (A) التجمعات الخلوية. (B) علامات التعليق التوضيحي. (C) تسع مجموعات خلوية مصنفة. (D, E) توزيعات الخلايا ونسبها في عينات الأصحاء وعينات UC. (F, G) التعبير عن المؤشرات الحيوية المرشحة موضحاً بواسطة UMAP ومخطط الفقاعات. UC: التهاب القولون التقرحي؛ UMAP: التقريب والإسقاط المتشعب الموحد. يرجى الضغط هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 8: التجميع الفرعي للخلايا الظهارية والتعبير المرشح. (A) التجمعات الأولية للخلايا الظهارية. (B) علامات التعليق التوضيحي. (C) المجموعات الفرعية المحددة من الخلايا الظهارية. (D) نسب المجموعات الفرعية في عينات الأصحاء وعينات UC. (E) تعبير LCN2 وNOS2 عبر المجموعات الفرعية للخلايا الظهارية. t-SNE: التضمين العشوائي للجيران الموزع t؛ UC: التهاب القولون التقرحي. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 9: التجميع الفرعي للخلايا النخاعية والتعبير المرشح. (A) عناقيد الخلايا النخاعية الأولية. (B) علامات التعليق التوضيحي. (C) المجموعات الفرعية للخلايا النخاعية الموضحة. (D) نسب المجموعات الفرعية في عينات الأصحاء وعينات التهاب القولون التقرحي UC. (E) تعبير IDO1 و CXCL2 و CD274 عبر المجموعات الفرعية للخلايا النخاعية. t-SNE: تضمين الجوار العشوائي الموزع t؛ UC: التهاب القولون التقرحي. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 10: تحليل الزمن الكاذب للخلايا الظهارية. (أ) مسار الزمن الكاذب وتعيينات الحالة. (ب) توزيع الخلايا الظهارية السليمة وتلك المصابة بـ UC على طول المسار. (ج) ديناميكيات تعبير LCN2 و NOS2 على طول الزمن الكاذب. UC: التهاب القولون التقرحي. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 11: تحليل الزمن الكاذب للخلايا النخاعية. (A) مسار الزمن الكاذب وتعيينات الحالة. (B) توزيع الخلايا النخاعية السليمة وخلايا خلايا التهاب القولون التقرحي على طول المسار. (C) ديناميكيات تعبير IDO1 و CXCL2 و CD274 على طول الزمن الكاذب. UC: التهاب القولون التقرحي. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 12: التواصل الظهاري-النخاعي. (A) عدد التفاعلات وقوتها. (B) قوة الإشارات الصادرة والواردة. (C) التفاعلات التي تشمل الخلايا القولونية الالتهابية. (D) خريطة حرارية لقوة التواصل. (E) أزواج الربيطة-المستقبل التي تشمل الخلايا القولونية الالتهابية. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
In vitro التحقق التجريبي من تدخلات الموت الخلوي الحديدي (ferroptosis) والموت الخلوي النحاسي (cuproptosis)
أدى التعرض لعديد السكاريد الشحمي (LPS) إلى تقليل حيوية خلايا Caco-2 وزيادة تعبير IL-6 وIL-1β مقارنة بمجموعة التحكم (الشكل 13A–C و الشكل التكميلي 5).
أدى LPS إلى زيادة Fe2⁺، والمالونديالدهيد، ووفرة الحمض الريبوزي المرسال (mRNA) لكل من LCN2 و IDO1 و CXCL2 و NOS2، بينما عكس ferrostatin-1 كل تغيير مرتبط بـ LPS (الشكل 13D–I). وبالنسبة للمقارنات بين LPS مقابل LPS + ferrostatin-1، تراوحت قيم P الثنائية الدقيقة من 9.45 × 10⁻5 إلى 0.0027 بعد حساب متوسط المكررات التقنية ضمن كل مكرر بيولوجي.
أدى LPS إلى تقليل الحيوية وزيادة كل من IL-6 وIL-1β وFDX1/DLAT وCD274 والفلورة الحساسة للنحاس، بينما عكس tetrathiomolybdate هذه التغيرات (الشكل 14 أ-هـ). بالنسبة للمقارنات بين LPS و LPS + tetrathiomolybdate، تم استخدام اختبار ثنائي الجانب دقيق P تراوحت القيم بين < 1 × 10⁻15 إلى 0.0008.
حددت التحليلات مجتمعةً خمسة مرشحات تشخيصية مرتبطة بالتهاب القولون التقرحي (UC)، وحصرت تعبيرها في المجموعات الظهارية والنخاعية، وأظهرت أن LCN2 وIDO1 وCXCL2 وNOS2 قد استجابت لتثبيط الموت الخلوي الحديدي (ferroptosis)، بينما استجاب CD274 لعملية مخلبية النحاس في خلايا Caco-2.

الشكل 13: يعمل Ferrostatin-1 على تخفيف التغيرات الالتهابية والتغيرات المرتبطة بـ ferroptosis والناتجة عن LPS في خلايا Caco-2. (A) حيوية الخلايا. (B, C) تعبير RNA المرسال لـ IL-6 وIL-1β. (D, E) مستويات malondialdehyde وFe2⁺ داخل الخلايا. (F–I) تعبير RNA المرسال لـ LCN2، وIDO1، وCXCL2، وNOS2. تم حساب متوسط المكررات التقنية ضمن كل مكرر بيولوجي مستقل؛ وتشير أشرطة الخطأ إلى الانحراف المعياري. كانت أعداد المكررات البيولوجية المتاحة للمجموعة الضابطة n = 4 في A وC، والمجموعة الضابطة n = 2 في F، وn = 6 لجميع تركيبات اللوحات/المجموعات الأخرى. تم تسجيل قيم P الدقيقة ثنائية الجانب وفحوصات الافتراضات في الجدول التكميلي 5. Fer-1: ferrostatin-1؛ LPS: lipopolysaccharide؛ MDA: malondialdehyde. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.

الشكل 14: مادة Tetrathiomolybdate تخفف من التغيرات المرتبطة بالالتهاب والموت النحاسي (cuproptosis) الناتجة عن LPS في خلايا Caco-2. (A) حيوية الخلايا. (B, C) تعبير الحمض النووي الريبوزي المرسال (mRNA) لـ IL-6 و IL-1β. (D) تعبير mRNA لـ FDX1 و DLAT. (E) تعبير mRNA لـ CD274. (F) الفلورة لـ Cu2⁺ داخل الخلايا؛ شريط المقياس = 50 µm. تم حساب متوسط ثلاث مكررات تقنية ضمن كل من ست مكررات بيولوجية مستقلة لكل مجموعة؛ تشير أشرطة الخطأ إلى الانحراف المعياري. قيم P الدقيقة ثنائية الجانب وفحوصات الافتراضات موضحة في الجدول التكميلي 5. TTM: tetrathiomolybdate؛ LPS: lipopolysaccharide. يرجى النقر هنا لعرض نسخة أكبر من هذا الشكل.
البيانات الخام والمعالجة، ونصوص التحليل البرمجية، وجداول بيانات مصدر تجارب الخلايا متاحة للعموم على الرابط https://doi.org/10.5281/zenodo.21202720.
الشكل التكميلي 1: المعالجة المسبقة للبيانات وتشخيصات التعبير التفاضلي. (A–F) توزيعات العينات قبل (A–C) وبعد (D–F) تصحيح الدفعة. (G) مخطط البركان للجينات ذات التعبير التفاضلي. (H) خريطة حرارية لأنماط التعبير المعيرة. DEGs: الجينات ذات التعبير التفاضلي؛ GEO: مستودع تعبير الجينات (Gene Expression Omnibus).يرجى النقر هنا لتحميل هذا الملف.
الشكل التكميلي 2: تحليل الإثراء للجينات ذات التعبير التفاضلي المشتركة في الموت النحاسي والموت الحديدي (CF-DEGs). (أ) إثراء الوجود الجيني (Gene Ontology). (ب) إثراء موسوعة كيوتو للجينات والجينومات (Kyoto Encyclopedia of Genes and Genomes). CF-DEGs: الجينات ذات التعبير التفاضلي المشتركة في الموت النحاسي والموت الحديدي؛ GO: الوجود الجيني؛ KEGG: موسوعة كيوتو للجينات والجينومات.يرجى النقر هنا لتنزيل هذا الملف.
الشكل التكميلي 3: تشخيصات فحص WGCNA. (A) اختيار العتبة اللينة. (B) تجميع الجينات الذاتية للوحدات. (C) المخطط الشجري للجينات والوحدات. (D) الارتباطات بين الوحدات والسمات. (E–G) علاقات الأهمية الجينية/عضوية الوحدة للوحدات الأرجوانية والبنية والسوداء. WGCNA: تحليل شبكة التعبير الجيني المشترك الموزون؛ UC: التهاب القولون التقرحي.يرجى النقر هنا لتحميل هذا الملف.
الشكل التكميلي 4: فحص المؤشرات الحيوية متعدد المراحل. (A, B) اختيار الميزات باستخدام LASSO. (C, D) اختيار الميزات باستخدام SVM-RFE. (E, F) ترتيب الميزات باستخدام الغابة العشوائية (Random-forest). (G) الإجماع بين النماذج الثلاثة. (H, I) شبكة PPI والعقد المركزية في CytoHubba. (J) دمج نتائج التعلم الآلي، وWGCNA، وPPI. (K) المواقع الكروموسومية للمرشحين الخمسة. LASSO: عامل تقليص واختيار أدنى مطلق؛ SVM-RFE: آلة متجه الدعم - استبعاد الميزات العودي؛ RF: الغابة العشوائية؛ PPI: التفاعل بين البروتينات؛ WGCNA: تحليل شبكة التعبير الجيني المشترك الموزون.يرجى النقر هنا لتنزيل هذا الملف.
الشكل التكميلي 5: التحقق من صحة نموذج LPS في خلايا Caco-2. (A) حيوية الخلايا. (B, C) تعبير الحمض الريبي النووي المرسال لـ IL-6 و IL-1β. تم حساب متوسط الثلاث مكررات التقنية ضمن ست مكررات بيولوجية مستقلة لكل مجموعة؛ تشير أشرطة الخطأ إلى الانحراف المعياري. LPS: عديد السكاريد الدهني.يرجى النقر هنا لتنزيل هذا الملف.
الجدول التكميلي 1: مجموعات بيانات GEO المستخدمة للاكتشاف والتحقق. أرقام الوصول إلى مجموعات البيانات، والمنصات، وتكوين العينات، وتعيينات المجموعات لمجموعات بيانات GEO المضمنة في تحليلات الاكتشاف والتحقق. GEO: Gene Expression Omnibus.يرجى النقر هنا لتنزيل هذا الملف.
الجدول التكميلي 2: تسلسلات بادئات تفاعل البوليميراز المتسلسل الكمي بنسخ عكسي. تسلسلات البادئات المستخدمة في تحليل تفاعل البوليميراز المتسلسل الكمي بنسخ عكسي للجينات التي تم تقييمها في هذه الدراسة.يرجى النقر هنا لتنزيل هذا الملف.
الجدول التكميلي 3: التفاصيل الإحصائية للأشكال 13 و14 والشكل التكميلي 1. أعداد المكررات البيولوجية، وفحوصات الافتراضات، وطرق التحليل الإحصائي، وقيم P الدقيقة للمقارنات التجريبية المحددة.يرجى النقر هنا لتحميل هذا الملف.
الجدول التكميلي 4: الجينات ذات التعبير التفاضلي المشتركة بين الموت النحاسي والموت الحديدي. قائمة بالجينات ذات التعبير التفاضلي المرتبطة بالموت الحديدي والنحاسي (CF-DEGs) التي تم تحديدها عن طريق تقاطع الجينات المرتبطة بالموت الحديدي والموت النحاسي مع الجينات ذات التعبير التفاضلي.يرجى النقر هنا لتحميل هذا الملف.
الجدول التكميلي 5: الجينات في وحدات WGCNA المختارة لفحص المؤشرات الحيوية. قائمة بالجينات الموجودة في وحدات WGCNA المختارة لفحص المؤشرات الحيوية اللاحق. WGCNA: تحليل شبكة التعبير الجيني المشترك الموزون.يرجى النقر هنا لتنزيل هذا الملف.