استخدمت هذه الدراسة مجموعات بيانات متاحة للجمهور وسلالات خلايا بشرية متوفرة تجارياً. لم يتم إشراك مشاركين بشريين، أو أنسجة بشرية جُمعت حديثاً، أو حيوانات فقارية؛ ولذلك، لم يتطلب الأمر موافقة لجنة أخلاقيات المؤسسة. تدرج أدوات البحث والبرمجيات والموارد المتاحة عبر الإنترنت المستخدمة في البروتوكول في جدول المواد.
1. معالجة البيانات لتحليل التعبير التفاضلي
لتأسيس مجموعات الاكتشاف والتحقق من النسخ المتماثل، تم البحث في قاعدة بيانات Gene Expression Omnibus (GEO) باستخدام المصطلح "ulcerative colitis". استُخدمت أربع مجموعات من البيانات وهي GSE87466 وGSE92415 وGSE107499 وGSE75214 كمجموعة تدريب، بينما خُصصت GSE47908 للتحقق الخارجي. تمت معايرة مجموعات بيانات التدريب وتصحيح تأثير الدفعة (batch-corrected) باستخدام ComBat في برنامج Sangerbox 3.0، بينما تمت معالجة مجموعة بيانات التحقق بشكل مستقل.
تم الحصول على ما مجموعه 583 جيناً مرتبطاً بـ ferroptosis (FRGs) من قاعدة بيانات FerrDB، وتم استكمالها بجينات وردت في الدراسات السابقة لتقليل التحيز الناجم عن الاعتماد على مصدر واحد12. كما تم تجميع مجموعة مكونة من 96 جيناً مرتبطاً بـ cuproptosis (CRGs) من دراسات منشورة12,13,14. حُلل التعبير التفاضلي باستخدام حزمة limma بناءً على عتبات |log₂ fold change| > 1 وقيمة P معدلة بطريقة Benjamini–Hochberg < 0.05. تتوفر المخططات التشخيصية للتطبيع (normalization)، وتصحيح الدفعة (batch correction)، والتعبير التفاضلي، والتجميع (clustering) في الشكل التكميلي 1. كما تتوفر تفاصيل تكوين العينات، والمنصات المستخدمة، وتوزيع المجموعات في الجدول التكميلي 1.
2. تحديد الجينات المعبّر عنها تفاضلياً في حالة التليف الكيسي (CF-DEGs) في القولون التقرحي (UC)
أُجري تحليل ارتباط بيرسون باستخدام عتبات |r| > 0.5 و P < 0.05 لتحديد الجينات المرتبطة بالحديد (FRGs) والجينات المرتبطة بالنحاس (CRGs) ذات التعبير المشترك. ثم تم تقاطع مجموعة الجينات الناتجة مع الجينات ذات التعبير التفاضلي (DEGs) لتحديد الجينات ذات التعبير التفاضلي والمشترك بين الموت الخلوي بالنحاس والموت الخلوي بالحديد (CF-DEGs).
3. تحليل الإثراء
أُجريت تحليلات الإثراء للوجودية الجينية (GO) وموسوعة كيوتو للجينات والجينومات (KEGG) باستخدام clusterProfiler ومنصة Bioinformatics.com.cn. واعتُمدت قيمة P < 0.05 المعدلة وفق طريقة Benjamini–Hochberg كحد لعتبة الدلالة الإحصائية. وتتوفر مخططات الإثراء في الشكل التكميلي 2.
4. تحليل شبكة التعبير الجيني المشترك الموزون (WGCNA)
تم تحليل البنية التركيبية (modular structure) الكامنة وراء ملفات التعبير الجيني باستخدام تحليل شبكة التعبير الجيني المشترك الموزون (WGCNA). كما تم تحديد العينات غير الطبيعية واستبعادها من خلال التجميع الهرمي (hierarchical clustering) مقترناً بدالة goodSamplesGenes. وعقب ذلك، تم تحديد قوة عتبة لينة (β) مناسبة لتلبية معيار الطوبولوجيا الخالية من القياس (scale-free topology) اللازم لبناء الشبكة.
استُخدمت بيانات التعبير المعالجة لبناء مصفوفة التداخل الطوبولوجي (TOM) لتحليل التعبير الجيني المشترك. وباستخدام منهجية القطع الديناميكي للشجرة، نُظمت الجينات ذات أنماط التعبير المتسقة للغاية في وحدات جينية، بحيث تحتوي كل وحدة على 100 جين على الأقل. كما حُللت الارتباطات بين الجينات المميزة للوحدات (MEs) والأنماط الظاهرية لتحديد الوحدات المرتبطة بعملية المرض واستخلاص مجموعات الجينات المقابلة لها. تتوفر تشخيصات الفرز الخاصة بتحليل شبكة التعبير الجيني المشترك المرجح (WGCNA) في الشكل التكميلي 3.
5. اختيار المؤشرات الحيوية
بهدف تحديد أولويات المرشحين للمتابعة التجريبية، تم تطبيق ثلاث خوارزميات للتعلم الآلي بالتوازي15. استخدم انحدار مشغل الاختيار والتقليص المطلق الأصغر (LASSO) لاختيار الميزات، واستخدمت الغابة العشوائية (random forest) مع ntree = 500 لترتيب أهمية الجينات، بينما استخدم استبعاد الميزات العودي لآلة متجه الدعم (SVM-RFE) لتحديد المجموعة الفرعية من الميزات ذات أدنى خطأ في التحقق المتقاطع. تتوفر تشخيصات الفحص الكاملة في الشكل التكميلي 5.
تم تقييم أهمية الجينات بناءً على مساهمة الجينات الفردية في دقة التصنيف، وجُمعت الجينات ذات التصنيف العالي لإنشاء مجموعة جينات أولية16. وقامت خوارزمية SVM-RFE بفحص الميزات وتحسين النموذج من خلال الحذف التكراري للميزات، مما سمح بتحديد المجموعة الفرعية المثلى من الميزات للتصنيف17. واعتُبرت الجينات المشتركة بين الخوارزميات الثلاث مؤشرات حيوية مرشحة تم تحديدها من خلال التعلم الآلي.
تم استخدام الجينات المعبّر عنها تفاضلياً في التهاب القولون (CF-DEGs) لبناء شبكة تفاعلات البروتين-بروتين (PPI)، واستُخدم برنامج CytoHubba لتصنيف العقد المركزية. وقد حدد التقاطع بين وحدات WGCNA الرئيسية، ومخرجات نماذج تعلم الآلة الثلاثة، والمجموعة المركزية لشبكة PPI خمسة مؤشرات حيوية مرشحة متكاملة وهي: LCN2 وIDO1 وCXCL2 وNOS2 وCD274. صُنفت هذه الجينات كمرشحات تشخيصية مرتبطة بالتهاب القولون التقرحي (UC). وأدت تجارب المثبطات اللاحقة إلى تصنيف LCN2 وIDO1 وCXCL2 وNOS2 كعلامات مستجيبة لتثبيط الموت الخلوي الحديدي (ferroptosis)، بينما صُنّف CD274 كعلامة مستجيبة لتثبيط الموت الخلوي النحاسي (cuproptosis). ولم تُعتبر هذه الجينات منظمات عليا مثبتة لأي من المسارين.
6. تطوير النموذج التنبؤي والتقييم متعدد الطرق
لتقييم إمكانية الترجمة السريرية لبصمة الجينات المحددة، تم إنشاء مخطط ترقيمي (nomogram) باستخدام حزمة rms لتقييم الأداء التشخيصي. وقد تم تقييم أداء النموذج باستخدام منحنى المعايرة، واختبار Hosmer–Lemeshow (HL)، ومؤشر توافق Harrell (C-index).
تم تقييم معايرة النموذج باستخدام اختبار HL ومتوسط الخطأ المطلق (MAE). وقد تم تحديد قيمة P > 0.05 وقيمة MAE < 0.1 كحدود دنيا لملاءمة النموذج المقبولة ودقته، على التوالي. كما تم إنشاء منحنيات خصائص تشغيل المستقبل (ROC)، وحساب المساحة تحت المنحنى (AUC) لتقييم الأداء التنبؤي للمخطط البياني (nomogram) لمرض UC. واستُخدمت مجموعة بيانات مستقلة لتقييم مدى قابلية تعميم النموذج واستقراره.
7. تحليل إثراء مجموعة الجينات (GSEA)
في مجموعة التدريب، تم إجراء تحليل ارتباط سبيرمان بين كل مؤشر حيوي وجميع الجينات الأخرى. ورُتبت نتائج الارتباط لإنشاء مجموعة جينات مرتبة. كما أُجري تحليل إثراء مجموعة الجينات (GSEA) باستخدام مجموعة الجينات c2.cp.kegg_medicus.v2025.1.Hs.symbols من قاعدة بيانات التواقيع الجزيئية (MSigDB) كمرجع. واعتُبرت المسارات البيولوجية التي كانت قيمة P لها < 0.05 مثرية بشكل ملحوظ.
8. تحليل المشهد المناعي
تم استخدام أداة CIBERSORT ومصفوفة التوقيع LM22 لإجراء تحليل فك الالتفاف (deconvolution analysis) لبيانات التعبير الجيني وتقدير نسب 22 نوعاً فرعياً من الخلايا المناعية. بعد ذلك، تم إجراء تحليل ارتباط سبيرمان (Spearman correlation analysis) لتقييم الارتباطات بين المؤشرات الحيوية الخمسة والأنواع الفرعية للخلايا المناعية.
9. تحليل الشبكة التنظيمية
تم الحصول على معلومات حول تنظيم عوامل النسخ (TF) والحمض النووي الريبوزي الميكروي (miRNA) المرتبطة بالجينات من قاعدتي بيانات ChEA3 وTarBase 9.0 المتكاملتين عبر منصة NetworkAnalyst. وقد تم تضمين جزيئات miRNA المرتبطة بعقدتين على الأقل في الشبكة التنظيمية.
10. تسلسل الحمض النووي الريبوزي في خلية واحدة (Single-cell RNA-seq)
تم دمج مجموعات بيانات الخلايا المفردة لمرض التهاب القولون التقرحي UC، وهي GSE116222 التي تضم ثلاث عينات سليمة وثلاث عينات مصابة بـ UC، وGSE214695 التي تضم ست عينات سليمة وست عينات مصابة بـ UC. تم الاحتفاظ بالخلايا التي تحتوي على 200–5,000 جين مكتشف، و>1,000 معرف جزيئي فريد، و< 30% من النسخ الميتوكوندرية، و< 3% من نسخ الهيموجلوبين، و< 50% من النسخ الريبوسومية. استُخدمت دالة NormalizeData لإجراء عملية التسوية (normalization)، واستُخدم برنامج Harmony لتصحيح تأثير الدفعة (batch correction) قبل البدء بعمليتي التجميع (clustering) والتوصيف (annotation).
بعد عملية الترشيح، تم تطبيع البيانات باستخدام دالة NormalizeData. واستُخدم برنامج Harmony لتقليل تأثيرات الدفعة (batch effects) بين العينات. كما استُخدم كل من التقريب والمسقط المتعدد الأبعاد الموحد (UMAP) والتضمين العشوائي للجيران الموزع ت (t-SNE) لخفض الأبعاد وتجميع الخلايا لتوصيف الهياكل المحلية والعامة ضمن مجتمع الخلايا. وقد تم الاحتفاظ بإجمالي 31,712 خلية عالية الجودة لإجراء المزيد من التحليلات. كما تم تعيين أنواع الخلايا وفقاً لجينات واسمات الخلايا القياسية.
تم إجراء تحليل الوقت الزائف (Pseudotime analysis) باستخدام برنامج Monocle 2 لاستقصاء التغيرات في حالات الخلايا الظهارية والنخاعية من أنسجة التهاب القولون التقرحي (UC). كما تم تقييم مسارات تطور الخلايا والاختلافات بين المجموعات، وتصوير أنماط تعبير المؤشرات الحيوية على طول مسار الوقت الزائف.
استُخدم برنامج CellChat لتقييم التواصل بين الخلايا بين المجموعات الفرعية للخلايا الظهارية والخلايا النخاعية. كما دُمجت أنماط التعبير الخاصة بكل مجموعة فرعية لاستنتاج احتمالات التفاعل ومسارات الإشارات الرئيسية بين الخلايا الظهارية والنخاعية، ولتوصيف شبكات التنظيم البينية المحتملة بين هذه الخلايا.
11. التجارب الخلوية
تم استزراع خلايا سرطان القولون الغدي البشري Caco-2 في وسط Dulbecco’s modified Eagle المحتوي على 10% fetal bovine serum، و100 U/mL penicillin، و100 µg/mL streptomycin عند درجة حرارة 37 °C في جو يحتوي على 5% CO₂. تم اختيار خلايا Caco-2 كنموذج طبقة ظهارية أحادية قابلة للتكرار لاختبارات الحاجز المعوي والإصابة الالتهابية؛ ومع ذلك، فإن أصلها الخبيث يحد من إمكانية التعميم المباشر على ظهارة القولون غير الخبيثة. لم تكن معلومات رقم الممر (Passage) متاحة في السجلات التجريبية المؤرشفة.
تم إنشاء نموذج إصابة التهابية in vitro عن طريق تعريض الخلايا لـ 1 µg/mL من lipopolysaccharide (LPS) لمدة 24 h. تم تقسيم الخلايا إلى مجموعة ضابطة باستخدام dimethyl sulfoxide كمركبة، ومجموعة LPS فقط، ومجموعة LPS + 4 µM ferrostatin-1، ومجموعة LPS + 40 µM tetrathiomolybdate، مع مطابقة التركيزات النهائية للمذيب. عملت مجموعة LPS فقط كمجموعة ضابطة إيجابية للإصابة/النموذج لإجراء مقارنات الإنقاذ باستخدام المثبطات.
- حيوية الخلايا
تم تقييم حيوية الخلايا باستخدام مقايسة حيوية الخلايا وفقاً لتعليمات الشركة المصنعة. وبعد الحضن عند 37 °C تم قياس الامتصاصية عند طول موجي 450 نانومتر باستخدام قارئ الأطباق الدقيقة لمدة ساعتين تقريباً.
- قياس MDA و Fe²⁺ مستويات
المالونديالدهيد (MDA) والحديد (Fe) داخل الخلايا2⁺ قِيست المستويات باستخدام أطقم المقايسة المقابلة وفقاً لإرشادات الشركات المصنعة. كما تم تقدير كمية البروتين الكلي باستخدام مقايسة حمض البيسينكونينينيك، وعُيِّرت كل قراءة وفقاً لمحتوى البروتين في العينة المقابلة.
- مقايسة مسبار فلوري لأيونات النحاس
بعد إزالة وسط الزرع، يتم...2⁺تمت إضافة المسبار الفلوري النوعي بتركيز نهائي قدره 1 µM وحُضنت لمدة 30 دقيقة عند 37 °Cوتم الحصول على الصور على الفور دون غسيل باستخدام مجهر فلوري مقلوب عند أطوال موجية للإثارة والانبعاث بلغت 510 و578 نانومتر، على التوالي. كما تم تثبيت إعدادات الحصول على الصور عبر جميع المجموعات، وعُير شريط مقياس الرسم في الصورة عند 50 µm.
- تفاعل البوليميراز المتسلسل الكمي في الوقت الحقيقي مع النسخ العكسي (RT-qPCR)
تم عزل الحمض النووي الريبوزي الكلي (Total RNA) باستخدام كاشف استخلاص RNA يعتمد على الفينول والغوانيدينيوم. ثم تم تخليق الحمض النووي الريبوزي منقوص الأكسجين المتمم (Complementary DNA) من 1 µg للحمض النووي الريبوزي (RNA) باستخدام طقم تخليق الحمض النووي المتمم للشريط الأول، واستُخدمت كواشف قائمة على صبغة SYBR Green لتفاعل البوليميراز المتسلسل في الوقت الحقيقي (real-time PCR). وقد أُجري كل تفاعل في ثلاث مكررات تقنية. وتتوفر تسلسلات البادئات في الجدول التكميلي 2، وتم حساب الوفرة النسبية للحمض النووي الريبوزي المرسال باستخدام طريقة 2^-ΔΔCt⁻ΔΔعتبة الدورة (Ct) الطريقة
12. التحليل الإحصائي
تم تقييم التوزيع الطبيعي وتجانس التباين باستخدام اختباري Shapiro–Wilk وLevene المرتكز على الوسيط، على التوالي. وحُللت البيانات التي استوفت كلا الافتراضين باستخدام اختبار t ثنائي الطرف لعينات غير مرتبطة، أو تحليل التباين أحادي الاتجاه متبوعاً باختبار Tukey للفروق المعنوية الصادقة. وفي حال عدم تحقق أي من الافتراضين، استُخدم اختبار Welch t أو تحليل التباين لـ Welch متبوعاً بمقارنات Games–Howell.
تُعرض النتائج في صورة متوسط حسابي ± الانحراف المعياري. أعداد المكررات البيولوجية الفعلية، الاختبار الثنائي الدقيق P تُدرج القيم، وعمليات التحقق من الافتراضات في الجدول التكميلي 3المؤرشفة الشكل 13 تحتوي ورقة المصدر على عدد أقل من مكررات الضبط للأجزاء A (ن = 4)، ج (يرجى تزويد النص المراد ترجمته. = 4)، و F (ن = 2)، بينما تحتوي جميع تركيبات اللوحات/المجموعات الأخرى على نـ = 6.
أُجريت التحليلات الحسابية باستخدام لغة R وسير عمل كل من limma، وWGCNA، وclusterProfiler، وglmnet، وrandomForest، وe1071، وpROC، وrms، وSeurat v4، وHarmony، وMonocle 2، وCellChat، وCytoscape/CytoHubba. وشملت الموارد المتاحة عبر الإنترنت FerrDB، وGEO، وSTRING، وChEA3، وTarBase 9.0، وNetworkAnalyst، وCIBERSORT/LM22، وMSigDB، وBioinformatics.com.cn، وSangerbox 3.0.