Voor deze studie zijn gebruikgemaakt van publiekelijk beschikbare datasets en commercieel verkrijgbare menselijke cellijnen. Er waren geen menselijke proefpersonen, recent verzameld menselijk weefsel of gewervelde dieren betrokken; daarom was goedkeuring van een institutionele ethische commissie niet vereist. De onderzoeksinstrumenten, software en online bronnen die in het protocol zijn gebruikt, staan vermeld in de Tabel met Materialen.
1. Gegevensverwerking voor differentieel expressieonderzoek
Om de transcriptomische discovery- en validatiecohorten vast te stellen, is de Gene Expression Omnibus (GEO) database doorzocht met de term “ulcerative colitis”. Vier datasets—GSE87466, GSE92415, GSE107499 en GSE75214—werden gebruikt als trainingsset, terwijl GSE47908 werd gereserveerd voor externe validatie. De trainingsdatasets werden genormaliseerd en batch-gecorrigeerd met ComBat in Sangerbox 3.0, terwijl de validatiedataset onafhankelijk werd verwerkt.
In totaal werden 583 ferroptose-gerelateerde genen (FRG's) verkregen uit FerrDB en aangevuld met genen die in de literatuur zijn beschreven om bias door een enkele bron te verminderen12Een set van 96 cuproptose-gerelateerde genen (CRG's) werd samengesteld op basis van gepubliceerde studies.12,13,14Differentiële expressie werd geanalyseerd met limma met drempelwaarden van |log₂ vouwverandering > 1 en Benjamini-Hochberg-gecorrigeerd P < 0.05. Diagnostische plots van normalisatie, batchcorrectie, differentiële expressie en clustering zijn opgenomen in Aanvullende figuur 1De monstersamenstelling, platforms en cohorttoewijzingen zijn vermeld in Aanvullende tabel 1.
2. Identificatie van CF-DEGs in UC
Pearson-correlatieanalyse werd uitgevoerd met drempelwaarden van |r| > 0,5 en P < 0,05 om co-expressie van FRGs en CRGs te identificeren. De resulterende genenset werd gekruist met de differentieel tot expressie gebrachte genen (DEGs) om cuproptose-ferroptose co-expressie differentieel tot expressie gebrachte genen (CF-DEGs) te identificeren.
3. Verrijkingsanalyse
Verrijkingsanalyses voor Gene Ontology (GO) en de Kyoto Encyclopedia of Genes and Genomes (KEGG) werden uitgevoerd met clusterProfiler en het Bioinformatics.com.cn platform. Een volgens Benjamini–Hochberg gecorrigeerde P < 0,05 werd gebruikt als significantiedrempel. De verrijkingsplots zijn opgenomen in Aanvullende Figuur 2.
4. Gewogen gen-co-expressienetwerkanalyse (WGCNA)
De modulaire structuur die ten grondslag ligt aan de genexpressieprofielen werd geanalyseerd met behulp van weighted gene co-expression network analysis (WGCNA). Abnormale monsters werden geïdentificeerd en uitgesloten via hiërarchische clustering in combinatie met de goodSamplesGenes-functie. Vervolgens werd een geschikte soft-thresholding power (β) bepaald om te voldoen aan het criterium van de schaalvrije topologie voor de netwerkconstructie.
De verwerkte expressiegegevens werden gebruikt om een topologische overlapmatrix (TOM) op te stellen voor gen-co-expressieanalyse. Met behulp van de dynamic tree-cutting benadering werden genen met zeer consistente expressiepatronen georganiseerd in genmodules, waarvan elke module ten minste 100 genen bevatte. Associaties tussen module eigengenes (MEs) en fenotypes werden geanalyseerd om modules te identificeren die geassocieerd zijn met het ziekteproces en om de overeenkomstige genensets te extraheren. De WGCNA-screeningsdiagnostiek zijn weergegeven in Aanvullende figuur 3.
5. Selectie van biomarkers
Om kandidaten voor experimentele follow-up te prioriteren, werden drie machine-learning-algoritmen parallel toegepast15. Least absolute shrinkage and selection operator (LASSO)-regressie werd gebruikt voor kenmerkselectie, random forest met ntree = 500 werd gebruikt om de belangrijkheid van genen te rangschikken, en support vector machine–recursive feature elimination (SVM-RFE) werd gebruikt om de subset van kenmerken met de laagste cross-validatiefout te identificeren. De volledige screeningsdiagnostiek is weergegeven in Aanvullende Figuur 5.
De belangrijkheid van genen werd geëvalueerd op basis van de bijdrage van individuele genen aan de classificatienauwkeurigheid, en hooggeplaatste genen werden verzameld om een voorlopige genenset te genereren16. Het SVM-RFE-algoritme voerde feature-screening en modeloptimalisatie uit via iteratieve feature-eliminatie, waardoor de optimale subset van features voor classificatie kon worden geïdentificeerd17. Genen die in alle drie de algoritmen voorkwamen, werden beschouwd als kandidaat-biomarkers geïdentificeerd via machine learning.
De CF-DEG's werden gebruikt voor de constructie van een eiwit-eiwitinteractienetwerk (PPI), waarbij CytoHubba werd aangewend om de kernknopen te rangschikken. De intersectie van de belangrijkste WGCNA-modules, de drie outputgegevens van machine learning en de PPI-kernset definieerde vijf geïntegreerde kandidaat-biomarkers: LCN2, IDO1, CXCL2, NOS2 en CD274. Deze genen werden geclassificeerd als UC-geassocieerde diagnostische kandidaten. Vervolgexperimenten met inhibitoren classificeerden LCN2, IDO1, CXCL2 en NOS2 verder als ferroptose-inhibitie-responsieve markers en CD274 als een cuproptose-inhibitie-responsieve marker. Deze genen werden niet beschouwd als bewezen upstream-regulatoren van een van beide pathways.
6. Ontwikkeling van het voorspellende model en evaluatie via meerdere methoden
Om het klinische translationele potentieel van de geïdentificeerde gen-signatuur te evalueren, is een nomogram geconstrueerd met behulp van het rms-pakket om de diagnostische prestaties te beoordelen. De prestaties van het model werden geëvalueerd met een kalibratiecurve, de Hosmer–Lemeshow (HL)-test en de concordantie-index (C-index) van Harrell.
De kalibratie van het model werd beoordeeld met behulp van de HL-test en de gemiddelde absolute fout (MAE). Een P-waarde > 0,05 en een MAE < 0,1 werden vastgesteld als drempelwaarden voor respectievelijk een aanvaardbare modelpassing en nauwkeurigheid. Er werden Receiver operating characteristic (ROC)-curven gegenereerd, en de oppervlakte onder de curve (AUC) werd berekend om de voorspellende prestaties van het nomogram voor UC te evalueren. Een onafhankelijke dataset werd gebruikt om de generaliseerbaarheid en stabiliteit van het model te beoordelen.
7. Gene set enrichment analyse (GSEA)
In de trainingset werd een Spearman-correlatieanalyse uitgevoerd tussen elke biomarker en alle andere genen. De correlatieresultaten werden gerangschikt om een geordende genenset te genereren. Gene set enrichment analysis (GSEA) werd uitgevoerd met de c2.cp.kegg_medicus.v2025.1.Hs.symbols genenset van de Molecular Signatures Database (MSigDB) als referentie. Biologische paden met P < 0,05 werden beschouwd als significant verrijkt.
8. Analyse van het immuunlandschap
De CIBERSORT-tool en de LM22-signatuurmatrix werden gebruikt om deconvolutieanalyse van de expressiegegevens uit te voeren en de proporties van 22 immuunceltypen te schatten. Vervolgens werd een Spearman-correlatieanalyse uitgevoerd om de associaties tussen de vijf biomarkers en de immuunceltypen te beoordelen.
9. Analyse van het regulatoire netwerk
Informatie over de regulatie door gen-gerelateerde transcriptiefactoren (TF) en microRNA (miRNA) werd verkregen uit de geïntegreerde ChEA3- en TarBase 9.0-databanken via het NetworkAnalyst-platform. MiRNA's die geassocieerd waren met ten minste twee knooppunten werden opgenomen in het regulatienetwerk.
10. RNA-seq van individuele cellen
De UC single-cell datasets GSE116222, bestaande uit drie gezonde en drie UC-monsters, en GSE214695, bestaande uit zes gezonde en zes UC-monsters, werden geïntegreerd. Cellen werden behouden wanneer zij 200–5.000 gedetecteerde genen, >1.000 unieke moleculaire identificatoren, < 30% mitochondriale transcripten, < 3% hemoglobinetranscripten en < 50% ribosomale transcripten bevatten. Voor de normalisatie werd de NormalizeData-functie gebruikt en Harmony werd ingezet voor batchcorrectie voorafgaand aan clustering en annotatie.
Na filtering werden de gegevens genormaliseerd met de functie NormalizeData. Harmony werd gebruikt om batch-effecten tussen monsters te minimaliseren. Uniform manifold approximation and projection (UMAP) en t-distributed stochastic neighbor embedding (t-SNE) werden toegepast voor dimensionaliteitsreductie en celclustering om lokale en globale structuren binnen de celpopulatie te karakteriseren. In totaal werden 31.712 hoogwaardige cellen behouden voor verdere analyse. Celtypen werden geannoteerd op basis van canonieke celmarkergenen.
Pseudotijd-analyse werd uitgevoerd met Monocle 2 om veranderingen in de toestanden van epitheliale en myeloïde cellen uit UC-weefsels te onderzoeken. Cellulaire ontwikkelingstrajecten en verschillen tussen groepen werden geëvalueerd, en expressiepatronen van biomarkers langs het pseudotijdstraject werden gevisualiseerd.
CellChat werd gebruikt om de cel-celcommunicatie tussen epitheliale en myeloïde cel-subsets te evalueren. Subpopulatie-specifieke expressiepatronen werden geïntegreerd om interactiekansen en belangrijke signaleringsroutes tussen epitheliale en myeloïde cellen af te leiden en om hun potentiële intercellulaire regulatoire netwerken te karakteriseren.
11. Celexperimenten
Humane Caco-2 colonadenocarcinoomcellen werden gekweekt in Dulbecco’s modified Eagle medium met 10% foetaal runderserum, 100 U/mL penicilline en 100 µg/mL streptomycine bij 37 °C in 5% CO₂. Caco-2-cellen werden gekozen als een reproduceerbaar epitheel monolayer-model voor intestinale barrière- en ontstekingsbeschadigingsassays; hun maligne oorsprong beperkt echter de directe generalisatie naar niet-maligne colonepitheel. Informatie over de passage was niet beschikbaar in de gearchiveerde experimentele verslagen.
Een in vitro ontsteking-schademodel werd vastgesteld door cellen 24 h lang bloot te stellen aan 1 µg/mL lipopolysacharide (LPS). Cellen werden toegewezen aan een dimethylsulfoxide-vehikelcontrolegroep, een LPS-only groep, een LPS + 4 µM ferrostatin-1 groep, of een LPS + 40 µM tetrathiomolybdaat groep, met overeenkomstige uiteindelijke oplosmiddelconcentraties. De LPS-only groep diende als positieve schade-/modelcontrole voor de inhibitor-rescue vergelijkingen.
- Celviabiliteit
De celviabiliteit werd beoordeeld met een celviabiliteitsassay volgens de instructies van de fabrikant. Na incubatie bij 37 °C gedurende ongeveer 2 u werd de absorbentie bij 450 nm gemeten met een microplate reader.
- Meting van MDA- en Fe²⁺-niveaus
Intracellulaire malondialdehyde (MDA)- en Fe2⁺-niveaus werden gemeten met de bijbehorende assaykits volgens de instructies van de fabrikanten. Totaal eiwit werd gekwantificeerd met een bicinchoninezuurassay, en elke meting werd genormaliseerd naar het eiwitgehalte van het betreffende monster.
- Fluorescente probe-assay voor koperionen
Na verwijdering van het cultuurmedium werd een Cu2⁺-specifieke fluorescente probe toegevoegd tot een eindconcentratie van 1 µM en 30 min geïncubeerd bij 37 °C. Beelden werden onmiddellijk zonder wassen vastgelegd met een omgekeerde fluorescentiemicroscoop bij excitatie- en emissiegolflengten van respectievelijk 510 en 578 nm. De acquisitie-instellingen werden constant gehouden voor alle groepen, en de schaalbalk van de afbeelding werd gekalibreerd op 50 µm.
- RT-qPCR
Totaal RNA werd geïsoleerd met een op fenol-guanidinium gebaseerd RNA-extractiereagens. Complementair DNA werd gesynthetiseerd uit 1 µg RNA met een first-strand complementair DNA-synthesekits, en op SYBR Green gebaseerde reagentia werden gebruikt voor real-time PCR. Elke reactie werd in technische triplicaat uitgevoerd. De primersequenties zijn opgenomen in Aanvullende Tabel 2, en de relatieve messenger RNA-overvloed werd berekend met de 2⁻ΔΔCt methode.
12. Statistische analyse
Normaliteit en homogeniteit van de variantie werden respectievelijk beoordeeld met de Shapiro–Wilk-toets en de mediaan-gecentreerde Levene-toets. Gegevens die aan beide aannames voldeden, werden geanalyseerd met een ongepairede tweezijdige t-toets of een eenweg-variantieanalyse gevolgd door de Tukey-toets voor eerlijk significante verschillen. Wanneer aan een van de aannames niet werd voldaan, werd de t-toets van Welch of de variantieanalyse van Welch gevolgd door Games–Howell-vergelijkingen gebruikt.
Resultaten worden weergegeven als gemiddelde ± standaarddeviatie. De werkelijke aantallen biologische replica's, exacte tweezijdige P-waarden en assumptiecontroles zijn vermeld in Aanvullende Tabel 3. Het gearchiveerde bronblad van Figuur 13 bevat minder controlereplica's voor panelen A (n = 4), C (n = 4) en F (n = 2), terwijl alle andere paneel-/groepcombinaties n = 6 bevatten.
Computationele analyses werden uitgevoerd met R en de workflows van limma, WGCNA, clusterProfiler, glmnet, randomForest, e1071, pROC, rms, Seurat v4, Harmony, Monocle 2, CellChat en Cytoscape/CytoHubba. Gebruikte online bronnen waren FerrDB, GEO, STRING, ChEA3, TarBase 9.0, NetworkAnalyst, CIBERSORT/LM22, MSigDB, Bioinformatics.com.cn en Sangerbox 3.0.