Diese Studie verwendete öffentlich zugängliche Datensätze und kommerziell erhältliche humane Zelllinien. Es waren keine menschlichen Probanden, neu gewonnene menschliche Gewebe oder Wirbeltiere beteiligt; daher war keine Genehmigung durch einen institutionellen Ethikausschuss erforderlich. In der Tabelle der Materialien sind die im Protokoll verwendeten Forschungswerkzeuge, Software und Online-Ressourcen aufgelistet.
1. Datenverarbeitung für die differentielle Expressionsanalyse
Um die transkriptomischen Entdeckungs- und Validierungskohorten aufzustellen, wurde die Gene Expression Omnibus (GEO)-Datenbank mit dem Begriff „ulcerative colitis“ durchsucht. Vier Datensätze – GSE87466, GSE92415, GSE107499 und GSE75214 – wurden als Trainingsset verwendet, während GSE47908 für die externe Validierung reserviert wurde. Die Trainingsdatensätze wurden mit ComBat in Sangerbox 3.0 normalisiert und auf Batch-Effekte korrigiert, während der Validierungsdatensatz unabhängig verarbeitet wurde.
Insgesamt wurden 583 eisenabhängige Zelltod-assozierte Gene (FRGs) aus FerrDB gewonnen und durch Gene ergänzt, die in der Literatur beschrieben wurden, um eine Einzelquellenverzerrung zu reduzieren12. Ein Satz von 96 kupferabhängige Zelltod-assozierten Genen (CRGs) wurde aus veröffentlichten Studien zusammengestellt12,13,14. Die differentielle Expression wurde mit limma unter Verwendung der Schwellenwerte |log₂-Fold-Change| > 1 und benjamini-hochberg-korrigiertem P < 0,05 analysiert. Diagnostische Plots zur Normalisierung, Batch-Korrektur, differentiellen Expression und Clusterung sind in Abbildung 1, Zusatzmaterial enthalten. Die Probenzusammensetzung, Plattformen und Kohortenzuordnungen sind in Tabelle 1, Zusatzmaterial aufgeführt.
2. Identifizierung von CF-DEGs bei UC
Die Pearson-Korrelationsanalyse wurde unter Verwendung von Schwellenwerten von |r| > 0,5 und P < 0,05 durchgeführt, um koexprimierte FRGs und CRGs zu identifizieren. Der resultierende Gen-Satz wurde mit den differentiell exprimierten Genen (DEGs) geschnitten, um kuproptotisch-ferroptotisch koexprimierte differentiell exprimierte Gene (CF-DEGs) zu ermitteln.
3. Anreicherungsanalyse
Genontologie (GO)- und Kyoto Encyclopedia of Genes and Genomes (KEGG)-Anreicherungsanalysen wurden mit clusterProfiler und der Bioinformatics.com.cn-Plattform durchgeführt. Ein nach Benjamini–Hochberg korrigierter P < 0,05 diente als Signifikanzschwelle. Anreicherungsdiagramme sind in Abbildung 2, Zusatzmaterial enthalten.
4. Gewichtete Gen-koexpressions-Netzwerkanalyse (WGCNA)
Die der Genexpressionsprofile zugrunde liegende modulare Struktur wurde mittels gewichteter Gen-Koexpressionsnetzwerkanalyse (WGCNA) analysiert. Anormale Proben wurden durch hierarchisches Clustering in Kombination mit der goodSamplesGenes-Funktion identifiziert und ausgeschlossen. Anschließend wurde eine geeignete Soft-Thresholding-Potenz (β) bestimmt, um das Kriterium der skalenfreien Topologie für die Netzwerkkonstruktion zu erfüllen.
Die aufbereiteten Expressionsdaten wurden verwendet, um eine topologische Überlappungsmatrix (TOM) für die Gen-Koexpressionsanalyse zu erstellen. Mithilfe des dynamischen Baumzerschneidungsansatzes wurden Gene mit stark übereinstimmenden Expressionsmustern in Genmodule gruppiert, wobei jedes Modul mindestens 100 Gene enthielt. Die Assoziationen zwischen Moduleigenengenen (MEs) und Phänotypen wurden analysiert, um Module zu identifizieren, die mit dem Krankheitsprozess assoziiert sind, und die entsprechenden Gensätze zu extrahieren. Die WGCNA-Screening-Diagnosen sind in Abbildung 3 (Zusatzmaterial) dargestellt.
5. Auswahl von Biomarkern
Um Kandidaten für die experimentelle Weiterverfolgung zu priorisieren, wurden drei maschinelle Lernalgorithmen parallel angewandt15. Für die Merkmalsauswahl wurde eine Least-Absolute-Shrinkage-and-Selection-Operator-(LASSO-)Regression verwendet, ein Random-Forest-Modell mit ntree = 500 diente zur Rangfolgebestimmung der Genwichtigkeit und ein Support-Vektor-Maschine–rekursiver Merkmalsausschluss (SVM-RFE) wurde genutzt, um die Merkmalsuntergruppe mit dem geringsten Kreuzvalidierungsfehler zu identifizieren. Die vollständigen Screening-Diagnosen sind in der ergänzenden Abbildung 5 dargestellt.
Die Genwichtigkeit wurde anhand des Beitrags einzelner Gene zur Klassifizierungsgenauigkeit bewertet, und hochrangige Gene wurden gesammelt, um einen vorläufigen Gen-Satz zu generieren¹⁶. Der SVM-RFE-Algorithmus führte eine Merkmalsauswahl und Modelloptimierung durch iterative Merkmalseliminierung durch, wodurch die optimale Teilmenge an Merkmalen für die Klassifizierung identifiziert werden konnte¹⁷. Gene, die von allen drei Algorithmen gemeinsam identifiziert wurden, galten als Kandidaten-Biomarker, die mittels maschinellem Lernen ermittelt wurden.
Die CF-DEGs wurden verwendet, um ein Protein-Protein-Interaktionsnetzwerk (PPI) zu konstruieren, und CytoHubba wurde zur Rangfolgebestimmung der zentralen Knotenpunkte genutzt. Der Durchschnitt der Schlüsselmodule der WGCNA, der drei maschinellen Lernausgaben und des PPI-Kernsets ergab fünf integrierte Kandidaten-Biomarker: LCN2, IDO1, CXCL2, NOS2 und CD274. Diese Gene wurden als UC-assoziierte diagnostische Kandidaten klassifiziert. Nachfolgende Inhibitorexperimente klassifizierten LCN2, IDO1, CXCL2 und NOS2 zusätzlich als Marker, die auf Hemmung der Ferroptose ansprechen, und CD274 als Marker, der auf Hemmung der Cuproptose anspricht. Diese Gene wurden nicht als nachgewiesene upstream-Regulatoren beider Signalwege betrachtet.
6. Entwicklung eines Vorhersagemodells und multimethodische Bewertung
Um das klinisch übertragbare Potenzial der identifizierten Gen-Signatur zu bewerten, wurde mithilfe des rms-Pakets ein Nomogramm zur Beurteilung der diagnostischen Leistung erstellt. Die Modellleistung wurde anhand einer Kalibrierungskurve, des Hosmer–Lemeshow-(HL-)Tests und des Harrellschen Konkordanzindexes (C-Index) evaluiert.
Die Modellkalibrierung wurde anhand des HL-Tests und des mittleren absoluten Fehlers (MAE) bewertet. Ein P-Wert > 0,05 und ein MAE < 0,1 wurden als Schwellenwerte für eine akzeptable Modellanpassung bzw. Genauigkeit festgelegt. Es wurden Receiver-Operating-Characteristic-(ROC)-Kurven erstellt, und die Fläche unter der Kurve (AUC) wurde berechnet, um die Vorhersageleistung des Nomogramms für UC zu bewerten. Ein unabhängiger Datensatz wurde verwendet, um die Generalisierbarkeit und Stabilität des Modells zu überprüfen.
7. Analyse der Gen-Set-Verstärkung (GSEA)
Im Trainingsset wurde eine Spearman-Korrelationsanalyse zwischen jedem Biomarker und allen anderen Genen durchgeführt. Die Korrelationsergebnisse wurden rangiert, um ein geordnetes Gen-Set zu generieren. Eine Gen-Set-Verdünnungsanalyse (GSEA) wurde unter Verwendung des c2.cp.kegg_medicus.v2025.1.Hs.symbols Gen-Sets aus der Molecular Signatures Database (MSigDB) als Referenz durchgeführt. Biologische Signalwege mit P < 0,05 galten als signifikant angereichert.
8. Analyse des Immunlandschaftsprofils
Das CIBERSORT-Tool und die LM22-Signaturmatrix wurden verwendet, um eine Dekonvolutionanalyse der Expressionsdaten durchzuführen und die Anteile von 22 Immunzellsubtypen zu schätzen. Anschließend wurde eine Spearman-Korrelationsanalyse durchgeführt, um die Assoziationen zwischen den fünf Biomarkern und den Immunzellsubtypen zu bewerten.
9. Analyse des regulatorischen Netzwerks
Informationen zur genbezogenen Regulation durch Transkriptionsfaktoren (TF) und Mikro-RNAs (miRNA) wurden über die integrierten Datenbanken ChEA3 und TarBase 9.0 über die Plattform NetworkAnalyst gewonnen. In das regulatorische Netzwerk wurden miRNAs einbezogen, die mit mindestens zwei Knoten assoziiert waren.
10. Einzelzell-RNA-Seq
Die Einzelzell-Datensätze GSE116222, bestehend aus drei gesunden und drei UC-Proben, und GSE214695, bestehend aus sechs gesunden und sechs UC-Proben, wurden integriert. Zellen wurden beibehalten, wenn sie 200–5.000 nachgewiesene Gene, >1.000 eindeutige molekulare Identifikatoren, < 30 % mitochondriale Transkripte, < 3 % Hämoglobin-Transkripte und < 50 % ribosomale Transkripte enthielten. Die NormalizeData-Funktion wurde zur Normalisierung verwendet, und Harmony wurde vor der Clusterung und Annotation für die Batch-Korrektur eingesetzt.
Nach der Filterung wurden die Daten mithilfe der NormalizeData-Funktion normalisiert. Harmony wurde verwendet, um Batch-Effekte zwischen den Proben zu minimieren. Die uniforme Mannigfaltigkeitsapproximation und Projektion (UMAP) sowie die t-verteilte stochastische Nachbarn-Einbettung (t-SNE) dienten der Dimensionsreduktion und Zellclustering, um lokale und globale Strukturen innerhalb der Zellpopulation zu charakterisieren. Insgesamt wurden 31.712 hochwertige Zellen für weitere Analysen beibehalten. Die Zelltypen wurden basierend auf kanonischen Zellmarkergenen annotiert.
Die Pseudotime-Analyse wurde mit Monocle 2 durchgeführt, um Veränderungen in den Zuständen epithelialer und myeloischer Zellen aus UC-Geweben zu untersuchen. Zelluläre Entwicklungstrajektorien und Gruppenunterschiede wurden evaluiert, und die Expressionsmuster von Biomarkern entlang der Pseudotime-Trajektorie wurden visualisiert.
CellChat wurde verwendet, um die Zell-Zell-Kommunikation zwischen epithelialen und myeloiden Zellsubpopulationen zu bewerten. Subpopulationsspezifische Expressionsmuster wurden integriert, um Interaktionswahrscheinlichkeiten und zentrale Signalwege zwischen epithelialen und myeloiden Zellen abzuleiten und ihre potenziellen zellulären Regulierungsnetzwerke zu charakterisieren.
11. Zellversuche
Humanen Caco-2-Kolonadenokarzinomzellen wurden in Dulbeccos modifiziertem Eagle-Medium, das 10 % fetalen Kälberserum, 100 U/mL Penicillin und 100 µg/mL Streptomycin enthielt, bei 37 °C in 5 % CO₂ kultiviert. Caco-2-Zellen wurden als reproduzierbares epitheliales Einschichtmodell für Darmbarrieren- und Entzündungs-Schädigungs-Assays ausgewählt; jedoch beschränkt ihr maligner Ursprung eine direkte Übertragung auf nichtmalignes kolonales Epithel. Informationen zur Passage waren in den archivierten experimentellen Aufzeichnungen nicht verfügbar.
Ein in vitro-Entzündungs-Schädigungsmodell wurde etabliert, indem Zellen 24 h lang 1 µg/mL Lipopolysaccharid (LPS) ausgesetzt wurden. Die Zellen wurden einer Dimethylsulfoxid-Fahrzeugkontrollgruppe, einer LPS-allein-Gruppe, einer LPS + 4 µM Ferrostatin-1-Gruppe oder einer LPS + 40 µM Tetrathiomolybdat-Gruppe zugeordnet, jeweils mit angepassten endgültigen Lösungsmittelkonzentrationen. Die LPS-allein-Gruppe diente als positive Schädigungs-/Modellkontrolle für die Inhibitor-Rettungsvergleiche.
- Zellviabilität
Die Zellviabilität wurde mithilfe eines Zellviabilitätsassays gemäß den Anweisungen des Herstellers bewertet. Nach einer Inkubation bei 37 °C für etwa 2 h wurde die Absorption bei 450 nm mit einem Mikroplatten-Reader gemessen.
- Bestimmung der MDA- und Fe²⁺-Konzentrationen
Die intrazellulären Konzentrationen von Malondialdehyd (MDA) und Fe2⁺ wurden mit den entsprechenden Testkits gemäß den Anweisungen der Hersteller bestimmt. Die Gesamtproteinmenge wurde mittels eines Bicinchoninsäure-Assays quantifiziert, und jede Messgröße wurde auf den Proteingehalt der entsprechenden Probe normalisiert.
- Fluoreszenzsonden-Assay für Kupferionen
Nach dem Entfernen des Kulturmediums wurde eine Cu2⁺-spezifische fluoreszierende Sonde in einer Endkonzentration von 1 µM hinzugefügt und 30 min bei 37 °C inkubiert. Bilder wurden unmittelbar ohne Waschschritte mit einem invertierten Fluoreszenzmikroskop bei Anregungs- und Emissionswellenlängen von jeweils 510 und 578 nm aufgenommen. Die Aufnahmeparameter blieben zwischen den Gruppen konstant, und die Maßstabsleiste im Bild wurde auf 50 µm kalibriert.
- RT-qPCR
Gesamt-RNA wurde mit einem auf Phenol–Guanidinium basierenden RNA-Extraktionsreagenz isoliert. Komplementäre DNA wurde aus 1 µg RNA unter Verwendung eines Kits zur Synthese des ersten DNA-Strangs hergestellt, und SYBR-Green-basierte Reagenzien wurden für die Echtzeit-PCR verwendet. Jede Reaktion wurde technisch dreifach durchgeführt. Die Primersequenzen sind in Supplementary Table 2 angegeben, und die relative mRNA-Menge wurde nach der 2⁻ΔΔCt-Methode berechnet.
12. Statistische Analyse
Die Normalverteilung und die Varianzhomogenität wurden mithilfe des Shapiro-Wilk-Tests bzw. des medianzentrierten Levene-Tests überprüft. Daten, die beide Annahmen erfüllten, wurden mit einem ungepaarten zweiseitigen t-Test oder einer einfaktoriellen Varianzanalyse, gefolgt vom ehrlichen signifikanten Differenztest nach Tukey, analysiert. Wenn eine der Annahmen nicht erfüllt war, kamen der t-Test nach Welch oder die Varianzanalyse nach Welch, gefolgt von den Vergleichen nach Games-Howell, zum Einsatz.
Die Ergebnisse werden als Mittelwert ± Standardabweichung angegeben. Die tatsächlichen Anzahlen biologischer Wiederholungen, exakten zweiseitigen P-Werte und Annahmenprüfungen sind in ergänzender Tabelle 3 aufgeführt. Die archivierte Quelldatei zu Abbildung 13 enthält für die Bereiche A (n = 4), C (n = 4) und F (n = 2) weniger Kontrollwiederholungen, während alle anderen Kombinationen aus Bereich und Gruppe n = 6 aufweisen.
Die computergestützten Analysen wurden mit R sowie den limma-, WGCNA-, clusterProfiler-, glmnet-, randomForest-, e1071-, pROC-, rms-, Seurat-v4-, Harmony-, Monocle-2-, CellChat- und Cytoscape/CytoHubba-Workflows durchgeführt. Zu den verwendeten Online-Ressourcen gehörten FerrDB, GEO, STRING, ChEA3, TarBase 9.0, NetworkAnalyst, CIBERSORT/LM22, MSigDB, Bioinformatics.com.cn und Sangerbox 3.0.