Kuproptozis–ferroptozis ko-ekspresyon gösteren genlerin tanımlanması ve fonksiyonel anotasyonu
421 UC örneği ve 97 sağlıklı kontrolden oluşan dört eğitim veri seti entegre edildi. 362'si artmış ve 189'u azalmış olmak üzere toplam 551 diferansiyel ekspresyon gösteren gen tanımlandı. Korelasyon analizi sonucunda 444 ferroptoz-kuproptoz ile ilişkili gen elde edildi ve diferansiyel ekspresyon gösteren genlerle yapılan kesişim analizi, 32 kuproptoz-ferroptoz ko-eksprese diferansiyel ekspresyon gösteren gen (CF-DEG) ortaya çıkardı.
32 CF-DEG; yaralanma yanıtları, bakır iyonları, yağ asidi taşınımı, lipopolisakkarit, inflamatuar bağırsak hastalığı, NF-κB sinyali, TNF sinyali ve ferroptozis süreçlerinde zenginleşmiştir (Ek Tablo 4).
Entegre makine öğrenimi ve WGCNA aracılığıyla temel biyobelirteçlerin önceliklendirilmesi
Ağırlıklı gen ko-ekspresyon ağ analizi (WGCNA), MEpurple, MEbrown ve MEblack modüllerini UC ile en güçlü şekilde ilişkili modüller olarak tanımlamıştır. Seçilen modüller toplamda 1.426 gen içermektedir (Ek Tablo 5).
En küçük mutlak büzülme ve seçim operatörü (LASSO), destek vektör makinesi-özyinelemeli özellik eleme (SVM-RFE) ve rastgele orman yöntemleri sırasıyla 21, 32 ve 19 özellik seçmiştir. On üç gen her üç model tarafından ortaklaşa seçilmiştir. Makine öğrenmesi konsensüsünün, CytoHubba ile sıralanmış protein-protein etkileşim çekirdeğinin ve seçilen WGCNA modüllerinin entegrasyonu; LCN2, IDO1, CXCL2, NOS2 ve CD274 olmak üzere UC ile ilişkili beş aday biyobelirteç sağlamıştır.
Biyobelirteç imzasının tanısal değerlendirmesi ve yol zenginleştirmesi
Eğitim kohortunda, beş aday biyobelirtecin tamamı pozitif korelasyon göstermiş ve ÜK'de yukarı regüle olmuştur (Şekil 2A, B). Her bir belirteç, alıcı işletim karakteristik eğrisi altında bir alan değerine ulaşmıştır > eğitim kohortunda 0,80 ve > bağımsız GSE47908 doğrulama kohortundaki 0,75 (Şekil 2C–E).
Beş genli nomogram, Hosmer–Lemeshow P > 0.05 ve ortalama mutlak hata < 0.1 ile uygun kalibrasyon; analiz edilen kohortta ise 0.945'lik bir eğri altı alan ve konkordans indeksi ile konkordans indeksi için 0.922–0.968'lik bir %95 güven aralığı ile ayırt edicilik göstermiştir (Şekil 3).
Gen seti zenginleştirme analizi, beş adayı JAK–STAT, interferon–RIPK1/3 ve Toll-benzeri reseptör–NF-κB sinyal yollarıyla ilişkilendirmiştir (Şekil 4).

Şekil 2: Aday biyobelirteç ekspresyonu ve tanısal performans. (A) Korelasyon ısı haritası. (B, C) Eğitim kohortundaki biyobelirteç ekspresyonu ve alıcı işletim karakteristiği eğrileri. (D, E) GSE47908 doğrulama kohortundaki biyobelirteç ekspresyonu ve alıcı işletim karakteristiği eğrileri. ROC: alıcı işletim karakteristiği; AUC: eğri altındaki alan; UC: ülseratif kolit. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 3: UC nomogramının geliştirilmesi ve değerlendirilmesi. (A) Beş genli nomogram. (B) Kalibrasyon grafiği. (C) Alıcı işletim karakteristik eğrisi. ROC: receiver operating characteristic; UC: ülseratif kolit; AUC: eğri altındaki alan; C-index: konkordans indeksi. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 4: Beş aday biyobelirtecin gen seti zenginleştirme analizi. (A–E) Sırasıyla LCN2, IDO1, CXCL2, NOS2 ve CD274 için gen seti zenginleştirme analizi sonuçları. GSEA: gen seti zenginleştirme analizi. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
İmmün mikroçevre manzarası ve düzenleyici ağ analizi
UC örnekleri; nötrofillerin, aktiveye memory CD4⁺ T hücrelerinin, M1 makrofajların ve aktive mast hücrelerinin oranlarında artış, buna karşılık M2 makrofajların, dinlenme halindeki mast hücrelerin ve dinlenme halindeki dendritik hücrelerin oranlarında azalma göstermiştir. Bu paternler doğrulama kohortunda da tekrarlanmıştır (Şekil 5A–E).
Beş aday; nötrofiller, aktive olmuş hafıza CD4⁺ T hücreleri ve M1 makrofajlar ile pozitif, dinlenme halindeki mast hücreleri ve M2 makrofajlar ile ise ters korelasyon gösterdi (Şekil 5F–J). Gen miRNA ağı 98 düğüm ve 146 kenar içeriyordu. hsa-miR-34a-5p ve hsa-miR-16-5p en yüksek biyobelirteç bağlantısallığını gösterirken, AR ve RELA en bağlantılı transkripsiyon faktörleriydi (Şekil 6).

Şekil 5: İmmün infiltrasyon analizi. (A) İmmün hücre kompozisyonu. (B) İmmün hücre oranlarındaki gruplar arası farklar. (C) İmmün hücre korelasyon ısı haritası. (D, E) Eğitim ve doğrulama kohortlarındaki immün hücre infiltrasyonu. (F–J) Aday biyobelirteçler ile immün hücre popülasyonları arasındaki korelasyonlar. UC: ülseratif kolit. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 6: Öngörülen miRNA ve transkripsiyon faktörü ağları. (A) Gen–miRNA ağı. Daireler aday biyobelirteçleri, kareler ise miRNA'ları göstermektedir. (B) Transkripsiyon faktörü–gen ağı. Baklava dilimleri aday biyobelirteçleri, ters üçgenler ise transkripsiyon faktörlerini göstermektedir. miRNA: mikroRNA; TF: transkripsiyon faktörü. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Tek hücre çözünürlüğünde spatiotemporal ekspresyon dinamikleri
Kalite kontrolünden sonra, 31.712 hücre, dokuz ana hücre popülasyonu olarak tanımlanan 22 küme oluşturmuştur (Şekil 7A–C). Epitel ve plazma hücresi oranları UC grubunda daha yüksekti. Epitel hücrelerinde LCN2 ve NOS2 zenginleşmişken, miyeloid hücrelerde IDO1, CXCL2 ve CD274 zenginleşmişti (Şekil 7D–G).
Epitel hücre alt kümeleme analizi 11 alt küme belirlemiş olup, bu alt kümede inflamatuar kolonositlerin genişlediği ve LCN2 ile NOS2'nin zenginleştiği görülmüştür (Şekil 8). Miyeloid hücre alt kümeleme analizi ise yedi alt küme belirlemiş olup, monositlerin arttığı, makrofajların azaldığı ve monositlerde IDO1, CXCL2 ve CD274'ün zenginleştiği tespit edilmiştir (Şekil 9).
İnflamatuar kolonositler, artan LCN2 ve NOS2 ekspresyonu ile epitelial yörüngenin geç evrelerinde birikmiştir (Şekil 10). Monositler; IDO1, CXCL2 ve CD274'ün dinamik ekspresyonu ile belirgin bir UC ile ilişkili yörünge sergilemiştir (Şekil 11).
İnflamatuar kolonositler en güçlü dışa dönük sinyalleşmeyi ve monositlerle belirgin bir iletişimi göstermiştir. APP–CD74, bu alt kümeler arasındaki baskın ligand-reseptör çiftiydi (Şekil 12).

Şekil 7: Hücre popülasyonları genelindeki aday biyobelirteç ekspresyonu. (A) Hücre kümeleri. (B) Anotasyon belirteçleri. (C) Anote edilmiş dokuz hücre popülasyonu. (D, E) Sağlıklı ve UC örneklerindeki hücre dağılımları ve oranları. (F, G) UMAP ve balon grafiği ile görselleştirilmiş aday biyobelirteç ekspresyonu. UC: ülseratif kolit; UMAP: uniform manifold approximation and projection. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 8: Epitel hücre alt kümelemesi ve aday ekspresyonu. (A) Başlangıç epitel hücre kümeleri. (B) Notasyon markörleri. (C) Notasyonlanmış epitel hücre alt kümeleri. (D) Sağlıklı ve UC örneklerindeki alt küme oranları. (E) Epitel hücre alt kümelerinde LCN2 ve NOS2 ekspresyonu. t-SNE: t-distributed stochastic neighbor embedding; UC: ülseratif kolit. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 9: Miyeloid hücre alt kümelemesi ve aday ekspresyonu. (A) Başlangıç miyeloid hücre kümeleri. (B) Annotasyon markerları. (C) Annote edilmiş miyeloid hücre alt kümeleri. (D) Sağlıklı ve UC örneklerindeki alt küme oranları. (E) Miyeloid hücre alt kümeleri genelinde IDO1, CXCL2 ve CD274 ekspresyonu. t-SNE: t-distributed stochastic neighbor embedding; UC: ülseratif kolit. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 10: Epitel hücre psödozaman analizi. (A) Psödozaman yörüngesi ve durum atamaları. (B) Sağlıklı ve UC epitel hücrelerinin yörünge boyunca dağılımı. (C) Psödozaman boyunca LCN2 ve NOS2 ekspresyon dinamikleri. UC: ülseratif kolit. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 11: Miyeloid hücre psödozaman analizi. (A) Psödozaman yörüngesi ve durum atamaları. (B) Sağlıklı ve UC miyeloid hücrelerinin yörünge boyunca dağılımı. (C) IDO1, CXCL2 ve CD274 ekspresyon dinamiklerinin psödozaman boyunca seyri. UC: ülseratif kolit. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 12: Epitel–miyeloid iletişimi. (A) Etkileşim sayısı ve gücü. (B) Giden ve gelen sinyalizasyon güçleri. (C) İnflamatuar kolonositleri içeren etkileşimler. (D) İletişim gücü ısı haritası. (E) İnflamatuar kolonositleri içeren ligand–reseptör çiftleri. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
In vitro ferroptozis ve kuproptozis müdahalelerinin deneysel doğrulaması
Lipopolisakkarit (LPS) maruziyeti, araç kontrolüne kıyasla Caco-2 hücre canlılığını azaltmış ve IL-6 ile IL-1β ekspresyonunu artırmıştır (Şekil 13A–C ve Ek Şekil 5).
LPS; Fe2⁺, malondialdehit ile LCN2, IDO1, CXCL2 ve NOS2 haberci RNA bolluğunu artırırken, ferrostatin-1 LPS ile ilişkili her bir değişikliği tersine çevirmiştir (Şekil 13D–I). LPS ile LPS + ferrostatin-1 karşılaştırmaları için, her biyolojik tekrar içindeki teknik tekrarların ortalaması alındıktan sonra kesin iki yönlü P değerleri 9,45 × 10⁻5 ile 0,0027 arasında değişmiştir.
LPS canlılığı azaltmış; IL-6, IL-1β, FDX1/DLAT, CD274 ve bakıra duyarlı floresansı artırmışken, tetratiyomolibdat bu değişiklikleri tersine çevirmiştir (Şekil 14 A–E). LPS ile LPS + tetratiyomolibdat karşılaştırmaları için kesin çift yönlü P değerleri < 1 × 10⁻15 ile 0,0008 arasında değişmiştir.
Analizler hep birlikte, UC ile ilişkili beş tanısal aday belirlemiş, bunların ekspresyonlarının epitelyal ve miyeloid popülasyonlarda yerelleştiğini saptamış ve Caco-2 hücrelerinde LCN2, IDO1, CXCL2 ve NOS2'nin ferroptoz inhibisyonuna, CD274'ün ise bakır şelasyonuna yanıt verdiğini göstermiştir.

Şekil 13: Ferrostatin-1, Caco-2 hücrelerinde LPS ile ilişkili inflamatuar ve ferroptozla ilgili değişiklikleri hafifletir. (A) Hücre canlılığı. (B, C) IL-6 ve IL-1β haberci RNA ekspresyonu. (D, E) Hücre içi malondialdehit ve Fe2⁺ seviyeleri. (F–I) LCN2, IDO1, CXCL2 ve NOS2 haberci RNA ekspresyonu. Teknik triplikatlar, her bağımsız biyolojik replika içinde ortalaması alınmıştır; hata çubukları standart sapmayı göstermektedir. Mevcut biyolojik replika sayıları A ve C'de Kontrol n = 4, F'de Kontrol n = 2 ve diğer tüm panel/grup kombinasyonları için n = 6'dır. Kesin iki yönlü P değerleri ve varsayım kontrolleri Ek Tablo 5'te bildirilmiştir. Fer-1: ferrostatin-1; LPS: lipopolisakkarit; MDA: malondialdehit. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 14: Tetratiyomolibdat, Caco-2 hücrelerinde LPS ile ilişkili inflamatuar ve kuproptozla ilgili değişiklikleri hafifletir. (A) Hücre canlılığı. (B, C) IL-6 ve IL-1β haberci RNA ekspresyonu. (D) FDX1 ve DLAT haberci RNA ekspresyonu. (E) CD274 haberci RNA ekspresyonu. (F) Hücre içi Cu2⁺ floresansı; ölçek çubuğu = 50 µm. Grup başına altı bağımsız biyolojik tekrardan her birinde teknik triplikaların ortalaması alınmıştır; hata çubukları standart sapmayı göstermektedir. Kesin iki taraflı P değerleri ve varsayım kontrolleri Ek Tablo 5'te bildirilmiştir. TTM: tetratiyomolibdat; LPS: lipopolisakkarit. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Ham ve işlenmiş veriler, analiz betikleri ve hücre deneyi kaynak tablolarına https://doi.org/10.5281/zenodo.21202720 adresinden açıkça erişilebilir.
Ek Şekil 1: Veri ön işleme ve diferansiyel ekspresyon tanıları. (A–F) Parti düzeltmesi öncesi (A–C) ve sonrası (D–F) örnek dağılımları. (G) Diferansiyel olarak eksprese edilen genlerin yanardağ grafiği. (H) Normalize edilmiş ekspresyon paternlerinin ısı haritası. DEGs: diferansiyel olarak eksprese edilen genler; GEO: Gene Expression Omnibus.Lütfen bu dosyayı indirmek için buraya tıklayın.
Ek Şekil 2: CF-DEG'lerin zenginleştirme analizi. (A) Gen Ontolojisi zenginleştirmesi. (B) Genler ve Genomlar Kyoto Ansiklopedisi zenginleştirmesi. CF-DEGs: kuproptozis-ferroptozis ortak ifade edilen farklı ifade edilen genler; GO: Gen Ontolojisi; KEGG: Genler ve Genomlar Kyoto Ansiklopedisi.Lütfen bu dosyayı indirmek için buraya tıklayın.
Ek Şekil 3: WGCNA tarama tanıları. (A) Yumuşak eşik seçimi. (B) Modül-eigengene kümelemesi. (C) Gen-modül dendrogramı. (D) Modül-özellik ilişkileri. (E–G) Mor, kahverengi ve siyah modüller için gen anlamlılığı/modül üyeliği ilişkileri. WGCNA: ağırlıklı gen ko-ekspresyon ağ analizi; UC: ülseratif kolit.Bu dosyayı indirmek için lütfen buraya tıklayın.
Ek Şekil 4: Çok aşamalı biyobelirteç taraması. (A, B) LASSO özellik seçimi. (C, D) SVM-RFE özellik seçimi. (E, F) Random-forest özellik sıralaması. (G) Üç model arasındaki konsensüs. (H, I) PPI ağı ve CytoHubba çekirdek düğümleri. (J) Makine öğrenimi, WGCNA ve PPI sonuçlarının entegrasyonu. (K) Beş adayın kromozomal konumları. LASSO: en küçük mutlak daralma ve seçim operatörü; SVM-RFE: destek vektör makinesi–özyinelemeli özellik eleme; RF: rastgele orman; PPI: protein–protein etkileşimi; WGCNA: ağırlıklı gen eş-ekspresyon ağ analizi.Lütfen bu dosyayı indirmek için buraya tıklayın.
Ek Şekil 5: Caco-2 hücrelerinde LPS model doğrulaması. (A) Hücre canlılığı. (B, C) IL-6 ve IL-1β haberci RNA ekspresyonu. Grup başına altı bağımsız biyolojik tekrarda teknik triplikelerin ortalaması alınmıştır; hata çubukları standart sapmayı göstermektedir. LPS: lipopolisakkarit.Bu dosyayı indirmek için lütfen buraya tıklayın.
Ek Tablo 1: Keşif ve doğrulama için kullanılan GEO veri setleri. Keşif ve doğrulama analizlerine dahil edilen GEO veri setleri için veri seti erişim numaraları, platformlar, örnek kompozisyonu ve kohort atamaları. GEO: Gene Expression Omnibus.Bu dosyayı indirmek için lütfen buraya tıklayınız.
Ek Tablo 2: Ters transkripsiyon kantitatif PCR primer dizileri. Bu çalışmada değerlendirilen genlerin ters transkripsiyon kantitatif PCR analizi için kullanılan primer dizileri.Bu dosyayı indirmek için lütfen buraya tıklayınız.
Ek Tablo 3: Şekil 13, 14 ve Ek Şekil 1 için istatistiksel ayrıntılar. Belirtilen deneysel karşılaştırmalar için biyolojik tekrarlama sayıları, varsayım kontrolleri, istatistiksel analiz yöntemleri ve kesin P değerleri.Bu dosyayı indirmek için lütfen buraya tıklayın.
Ek Tablo 4: Kuproptozis–ferroptozis ortak eksprese edilen diferansiyel eksprese olan genler. Ferroptozis–kuproptozis ile ilişkili genlerin, diferansiyel eksprese olan genlerle kesiştirilmesiyle belirlenen CF-DEG listesi.Lütfen bu dosyayı indirmek için buraya tıklayın.
Ek Tablo 5: Biyobelirteç taraması için seçilen WGCNA modüllerindeki genler. Sonraki biyobelirteç taraması için seçilen WGCNA modüllerinde yer alan genlerin listesi. WGCNA: ağırlıklı gen eş-ifade ağ analizi.Lütfen bu dosyayı indirmek için buraya tıklayın.