Bu çalışmada halka açık veri setleri ve ticari olarak temin edilebilen insan hücre hatları kullanılmıştır. Herhangi bir insan katılımcı, yeni toplanmış insan dokusu veya omurgalı hayvan kullanılmadığı için kurumsal etik kurul onayı gerekmemiştir. Protokolde kullanılan araştırma araçları, yazılımlar ve çevrimiçi kaynaklar Materyal Tablosu'nda listelenmiştir.
1. Diferansiyel ekspresyon analizi için veri işleme
Transkriptomik keşif ve doğrulama kohortlarını oluşturmak için Gene Expression Omnibus (GEO) veri tabanında “ulcerative colitis” terimi kullanılarak arama yapılmıştır. GSE87466, GSE92415, GSE107499 ve GSE75214 olmak üzere dört veri seti eğitim seti olarak kullanılırken, GSE47908 harici doğrulama için ayrılmıştır. Eğitim veri setleri Sangerbox 3.0 içindeki ComBat kullanılarak normalize edilmiş ve grup etkisinden (batch effect) arındırılmış, doğrulama veri seti ise bağımsız olarak işlenmiştir.
Tek kaynak yanlılığını azaltmak amacıyla FerrDB'den toplam 583 ferroptozis ilişkili gen (FRG) elde edilmiş ve literatürde bildirilen genlerle desteklenmiştir12. Yayınlanmış çalışmalardan 96 kuproptozis ilişkili gen (CRG) içeren bir set derlenmiştir12,13,14. Diferansiyel ekspresyon analizi, |log₂ fold change| > 1 ve Benjamini-Hochberg düzeltmeli P < 0.05 eşik değerleri kullanılarak limma ile gerçekleştirilmiştir. Normalizasyon, grup düzeltmesi (batch correction), diferansiyel ekspresyon ve kümeleme analizlerine ait tanısal grafikler Ek Şekil 1'de sunulmuştur. Örnek kompozisyonu, platformlar ve kohort atamaları Ek Tablo 1'de verilmiştir.
2. UC'de CF-DEGs'lerin Belirlenmesi
Ortak eksprese edilen FRG'ler ve CRG'leri belirlemek amacıyla |r| > 0.5 ve P < 0.05 eşik değerleri kullanılarak Pearson korelasyon analizi gerçekleştirilmiştir. Elde edilen gen seti, kuproptoz-ferroptoz ortak eksprese edilen diferansiyel eksprese genleri (CF-DEG'ler) belirlemek için diferansiyel eksprese genler (DEG'ler) ile kesiştirilmiştir.
3. Zenginleştirme analizi
Gen Ontolojisi (GO) ve Kyoto Gen ve Genom Ansiklopedisi (KEGG) zenginleştirme analizleri, clusterProfiler ve Bioinformatics.com.cn platformu kullanılarak gerçekleştirilmiştir. Anlamlılık eşiği olarak Benjamini–Hochberg düzeltmeli P < 0.05 değeri kullanılmıştır. Zenginleştirme grafikleri Ek Şekil 2'de sunulmuştur.
4. Ağırlıklı gen ko-ekspresyon ağ analizi (WGCNA)
Gen ekspresyon profillerinin temelini oluşturan modüler yapı, ağırlıklı gen eş-ekspresyon ağ analizi (WGCNA) kullanılarak incelenmiştir. Anormal örnekler, hiyerarşik kümeleme ve goodSamplesGenes fonksiyonu kombinasyonuyla belirlenmiş ve analiz dışı bırakılmıştır. Ardından, ağ kurulumu için ölçeksiz topoloji kriterini karşılayacak uygun bir yumuşak eşikleme gücü (β) belirlenmiştir.
İşlenmiş ekspresyon verileri, gen ko-ekspresyon analizi için bir topolojik örtüşme matrisi (TOM) oluşturmak amacıyla kullanıldı. Dinamik ağaç kesme yaklaşımı kullanılarak, oldukça tutarlı ekspresyon modellerine sahip genler, her biri en az 100 gen içeren gen modülleri halinde organize edildi. Hastalık süreciyle ilişkili modülleri belirlemek ve ilgili gen setlerini çıkarmak için modül özgenleri (ME'ler) ve fenotipler arasındaki ilişkiler analiz edildi. WGCNA tarama tanıları Ek Şekil 3'te sunulmuştur.
5. Biyobelirteçlerin seçimi
Deneysel takip için adayları önceliklendirmek amacıyla, üç makine öğrenimi algoritması paralel olarak uygulanmıştır15. Özellik seçimi için en küçük mutlak küçültme ve seçme operatörü (LASSO) regresyonu, gen önem sırasını belirlemek için ntree = 500 ile rastgele orman (random forest) ve en düşük çapraz doğrulama hatasına sahip özellik alt kümesini tanımlamak için destek vektör makinesi-özyinelemeli özellik eleme (SVM-RFE) kullanılmıştır. Tam tarama tanıları Ek Şekil 5'te sunulmuştur.
Gen önemi, tekil genlerin sınıflandırma doğruluğuna olan katkısına göre değerlendirilmiş ve öncül bir gen seti oluşturmak için yüksek sıralamadaki genler toplanmıştır16. SVM-RFE algoritması, iteratif özellik elemesi yoluyla özellik taraması ve model optimizasyonu gerçekleştirerek sınıflandırma için optimal özellik alt kümesinin belirlenmesini sağlamıştır17. Üç algoritma arasında ortak olan genler, makine öğrenimi yoluyla belirlenen aday biyobelirteçler olarak kabul edilmiştir.
CF-DEG'ler bir protein-protein etkileşim (PPI) ağı oluşturmak için kullanılmış ve çekirdek düğümleri sıralamak için CytoHubba kullanılmıştır. Temel WGCNA modüllerinin, üç makine öğrenimi çıktısının ve PPI çekirdek setinin kesişimi; LCN2, IDO1, CXCL2, NOS2 ve CD274 olmak üzere beş entegre aday biyobelirteci tanımlamıştır. Bu genler, UC ile ilişkili tanı adayları olarak sınıflandırılmıştır. Müteakip inhibitör deneyleri, LCN2, IDO1, CXCL2 ve NOS2'yi ferroptoz inhibisyonuna yanıt veren belirteçler, CD274'ü ise kuproptoz inhibisyonuna yanıt veren bir belirteç olarak daha ileri düzeyde sınıflandırmıştır. Bu genler, her iki yolun da kanıtlanmış üst regülatörleri olarak değerlendirilmemiştir.
6. Öngörücü model geliştirme ve çok yöntemli değerlendirme
Belirlenen gen imzasının klinik translasyonel potansiyelini değerlendirmek amacıyla, tanısal performansı ölçmek için rms paketi kullanılarak bir nomogram oluşturulmuştur. Model performansı; bir kalibrasyon eğrisi, Hosmer–Lemeshow (HL) testi ve Harrell uyum indeksi (C-indeksi) kullanılarak değerlendirilmiştir.
Model kalibrasyonu, HL testi ve ortalama mutlak hata (MAE) kullanılarak değerlendirildi. Kabul edilebilir model uyumu ve doğruluğu için eşik değerleri sırasıyla P değeri > 0,05 ve MAE < 0,1 olarak belirlendi. Nomogramın UC için öngörücü performansını değerlendirmek amacıyla alıcı işletim karakteristik (ROC) eğrileri oluşturuldu ve eğri altında kalan alan (AUC) hesaplandı. Modelin genellenebilirliğini ve stabilitesini değerlendirmek için bağımsız bir veri seti kullanıldı.
7. Gen seti zenginleştirme analizi (GSEA)
Eğitim setinde, her bir biyobelirteç ile diğer tüm genler arasında Spearman korelasyon analizi gerçekleştirilmiştir. Korelasyon sonuçları, sıralı bir gen seti oluşturmak amacıyla derecelendirilmiştir. Referans olarak Molecular Signatures Database (MSigDB) bünyesindeki c2.cp.kegg_medicus.v2025.1.Hs.symbols gen seti kullanılarak Gen Seti Zenginleştirme Analizi (GSEA) yapılmıştır. P < 0.05 olan biyolojik yolaklar anlamlı düzeyde zenginleşmiş olarak kabul edilmiştir.
8. İmmün manzara analizi
Ekspresyon verilerinin dekonvolüsyon analizini gerçekleştirmek ve 22 immün hücre alt tipinin oranlarını tahmin etmek için CIBERSORT aracı ve LM22 imza matrisi kullanılmıştır. Beş biyobelirteç ile immün hücre alt tipleri arasındaki ilişkileri değerlendirmek için ardından Spearman korelasyon analizi yapılmıştır.
9. Düzenleyici ağ analizi
Genle ilişkili transkripsiyon faktörü (TF) ve mikroRNA (miRNA) regülasyonuna ilişkin bilgiler, NetworkAnalyst platformu aracılığıyla entegre ChEA3 ve TarBase 9.0 veri tabanlarından elde edilmiştir. En az iki düğümle ilişkili olan miRNA'lar regülasyon ağına dahil edilmiştir.
10. Tek hücreli RNA-seq
Üç sağlıklı ve üç UC örneğinden oluşan GSE116222 ile altı sağlıklı ve altı UC örneğinden oluşan GSE214695 UC tek hücre veri setleri entegre edildi. 200–5.000 arası saptanan gen, >1.000 benzersiz moleküler tanımlayıcı, < %30 mitokondriyal transkript, < %3 hemoglobin transkripti ve < %50 ribozomal transkript içeren hücreler tutuldu. Normalizasyon için NormalizeData fonksiyonu, kümeleme ve anotasyon öncesindeki grup düzeltmesi (batch correction) için ise Harmony kullanıldı.
Filtreleme işleminden sonra, veriler NormalizeData fonksiyonu kullanılarak normalize edildi. Örnekler arasındaki batch etkilerini minimize etmek için Harmony kullanıldı. Hücre popülasyonu içindeki lokal ve global yapıları karakterize etmek amacıyla boyut indirgeme ve hücre kümeleme için Uniform manifold approximation and projection (UMAP) ve t-distributed stochastic neighbor embedding (t-SNE) kullanıldı. Daha ileri analizler için toplam 31.712 yüksek kaliteli hücre tutuldu. Hücre tipleri, kanonik hücre belirteç genlerine göre tanımlandı.
UC dokularındaki epitelyal ve miyeloid hücrelerin durumlarındaki değişiklikleri incelemek için Monocle 2 kullanılarak psödozaman analizi gerçekleştirildi. Hücre gelişim yörüngeleri ve gruplar arası farklar değerlendirildi ve psödozaman yörüngesi boyunca biyobelirteç ekspresyon modelleri görselleştirildi.
Epitel ve miyeloid hücre alt kümeleri arasındaki hücre-hücre iletişimini değerlendirmek için CellChat kullanılmıştır. Epitel ve miyeloid hücreler arasındaki etkileşim olasılıklarını ve temel sinyal yolaklarını belirlemek ve potansiyel hücreler arası düzenleyici ağlarını karakterize etmek için alt popülasyona özgü ekspresyon paternleri entegre edilmiştir.
11. Hücre deneyleri
İnsan Caco-2 kolon adenokarsinom hücreleri; %10 fetal bovine serum, 100 U/mL penisilin ve 100 µg/mL streptomisin içeren Dulbecco’s modified Eagle besiyerinde, 37 °C sıcaklıkta ve %5 CO₂ ortamında kültüre edilmiştir. Caco-2 hücreleri, intestinal bariyer ve inflamatuar hasar analizleri için tekrarlanabilir bir epitelyal monolayer modeli olarak seçilmiştir; ancak, malign kökenleri bu sonuçların malign olmayan kolon epitelyumuna doğrudan genellenmesini sınırlar. Arşivlenmiş deneysel kayıtlarda pasaj bilgisine ulaşılamamıştır.
Hücrelerin 24 saat boyunca 1 µg/mL lipopolisakkarit (LPS) maruziyetine bırakılmasıyla bir in vitro inflamatuar hasar modeli oluşturuldu. Hücreler; eşleşmiş nihai çözücü konsantrasyonları ile bir dimetil sülfoksit taşıyıcı kontrol grubuna, yalnızca LPS grubuna, LPS + 4 µM ferrostatin-1 grubuna veya LPS + 40 µM tetratiyomolibdat grubuna ayrıldı. Yalnızca LPS grubuna, inhibitör kurtarma karşılaştırmaları için pozitif hasar/model kontrolü olarak hizmet edildi.
- Hücre canlılığı
Hücre canlılığı, üretici talimatlarına uygun olarak bir hücre canlılığı testi kullanılarak değerlendirilmiştir. Yaklaşık 2 saat boyunca 37 °C'de inkübasyon sonrası, 450 nm'deki absorbans bir mikroplaka okuyucu kullanılarak ölçülmüştür.
- MDA ve Fe²⁺ seviyelerinin ölçümü
Hücre içi malondialdehit (MDA) ve Fe2⁺ seviyeleri, üretici talimatlarına uygun olarak ilgili analiz kitleri kullanılarak ölçülmüştür. Toplam protein miktarı bir bisinkoninik asit testi kullanılarak belirlenmiş ve her okuma, ilgili örneğin protein içeriğine göre normalize edilmiştir.
- Bakır iyonu floresan prob analizi
Kültür ortamı uzaklaştırıldıktan sonra, 1 µM final konsantrasyonunda Cu2⁺-spesifik bir floresan prob eklenmiş ve 37 °C'de 30 dakika inkübe edilmiştir. Görüntüler, yıkama işlemi yapılmadan, sırasıyla 510 ve 578 nm eksitasyon ve emisyon dalga boylarında bir ters floresans mikroskobu kullanılarak anlık olarak elde edilmiştir. Görüntüleme ayarları tüm gruplar için sabit tutulmuş ve görüntü ölçek çubuğu 50 µm olarak kalibre edilmiştir.
- RT-qPCR
Toplam RNA, fenol-guanidinyum tabanlı bir RNA ekstraksiyon reaktifi kullanılarak izole edilmiştir. Komplementer DNA, 1 µg RNA'dan bir birinci zincir komplementer DNA sentez kiti kullanılarak sentezlenmiş ve gerçek zamanlı PCR için SYBR Green tabanlı reaktifler kullanılmıştır. Her reaksiyon teknik olarak üç tekrarlı (triplicate) gerçekleştirilmiştir. Primer dizileri Ek Tablo 2'de verilmiş olup, göreceli haberci RNA bolluğu 2⁻ΔΔCt yöntemi kullanılarak hesaplanmıştır.
12. İstatistiksel analiz
Normallik ve varyans homojenliği, sırasıyla Shapiro-Wilk ve medyan merkezli Levene testleri kullanılarak değerlendirildi. Her iki varsayımı da karşılayan veriler, eşleşmemiş iki yönlü t-testi veya ardından Tukey'nin dürüst anlamlı fark testi ile tek yönlü varyans analizi kullanılarak analiz edildi. Varsayımlardan herhangi birinin karşılanmadığı durumlarda, Welch t-testi veya ardından Games-Howell karşılaştırmaları ile Welch varyans analizi kullanıldı.
Sonuçlar ortalama ± standart sapma olarak sunulmuştur. Gerçek biyolojik replika sayıları, kesin iki yönlü P değerleri ve varsayım kontrolleri Ek Tablo 3'te raporlanmıştır. Arşivlenen Şekil 13 kaynak sayfası, A (n = 4), C (n = 4) ve F (n = 2) panelleri için daha az kontrol replikası içerirken, diğer tüm panel/grup kombinasyonları n = 6 içermektedir.
Hesaplamalı analizler R ve limma, WGCNA, clusterProfiler, glmnet, randomForest, e1071, pROC, rms, Seurat v4, Harmony, Monocle 2, CellChat ve Cytoscape/CytoHubba iş akışları kullanılarak gerçekleştirilmiştir. Çevrimiçi kaynaklar arasında FerrDB, GEO, STRING, ChEA3, TarBase 9.0, NetworkAnalyst, CIBERSORT/LM22, MSigDB, Bioinformatics.com.cn ve Sangerbox 3.0 yer almıştır.