Questo studio ha utilizzato set di dati pubblicamente disponibili e linee cellulari umane reperibili in commercio. Non sono stati coinvolti partecipanti umani, tessuti umani appena raccolti o animali vertebrati; pertanto, non è stata richiesta l'approvazione da parte del comitato etico istituzionale. Gli strumenti di ricerca, il software e le risorse online utilizzati nel protocollo sono elencati nella Tabella dei Materiali.
1. Elaborazione dei dati per l'analisi dell'espressione differenziale
Per definire i cohorti di scoperta e di validazione transcriptomici, è stata effettuata una ricerca nel database Gene Expression Omnibus (GEO) utilizzando il termine „colite ulcerosa“. Quattro dataset — GSE87466, GSE92415, GSE107499 e GSE75214 — sono stati utilizzati come set di addestramento, mentre il dataset GSE47908 è stato riservato per la validazione esterna. I dataset di addestramento sono stati normalizzati e corretti per effetti di batch mediante ComBat in Sangerbox 3.0, mentre il dataset di validazione è stato processato in modo indipendente.
Un totale di 583 geni associati alla ferroptosi (FRGs) è stato ottenuto da FerrDB e integrato con geni riportati in letteratura per ridurre il bias derivante da un'unica fonte12. Un insieme di 96 geni associati alla cuproptosi (CRGs) è stato raccolto da studi pubblicati12,13,14. L'espressione differenziale è stata analizzata utilizzando limma con soglie di |variazione del log₂| > 1 e P corretto con il metodo di Benjamini–Hochberg < 0,05. I grafici diagnostici relativi alla normalizzazione, alla correzione del batch, all'espressione differenziale e al clustering sono forniti nella Figura Supplementare 1. La composizione dei campioni, le piattaforme e le assegnazioni dei cohort sono riportate nella Tabella Supplementare 1.
2. Identificazione di CF-DEG nella UC
È stata eseguita un'analisi della correlazione di Pearson utilizzando soglie di |r| > 0,5 e P < 0,05 per identificare i FRG e i CRG co-espressi. L'insieme genico risultante è stato incrociato con i geni differenzialmente espressi (DEG) per identificare i geni differenzialmente espressi co-espressi in cuproptosi e ferroptosi (CF-DEG).
3. Analisi di arricchimento
Le analisi di arricchimento di Gene Ontology (GO) e Kyoto Encyclopedia of Genes and Genomes (KEGG) sono state eseguite utilizzando clusterProfiler e la piattaforma Bioinformatics.com.cn. È stato utilizzato un valore di P corretto con il metodo di Benjamini–Hochberg < 0,05 come soglia di significatività. I grafici di arricchimento sono forniti nella Figura supplementare 2.
Analisi della rete di co-espressione genica pesata (WGCNA)
La struttura modulare sottostante i profili di espressione genica è stata analizzata mediante l'analisi di rete di co-espressione genica pesata (WGCNA). I campioni anomali sono stati identificati ed esclusi attraverso il clustering gerarchico combinato con la funzione goodSamplesGenes. Successivamente è stata determinata una potenza di soglia morbida (β) adeguata per soddisfare il criterio di topologia libera da scala nella costruzione della rete.
I dati di espressione elaborati sono stati utilizzati per costruire una matrice di sovrapposizione topologica (TOM) per l'analisi della co-espressione genica. Mediante l'approccio di taglio dinamico dell'albero, i geni con pattern di espressione altamente coerenti sono stati raggruppati in moduli genici, ciascuno contenente almeno 100 geni. Le associazioni tra eigengeni dei moduli (ME) e fenotipi sono state analizzate per identificare i moduli associati al processo della malattia ed estrarre i corrispondenti insiemi genici. I risultati diagnostici della selezione WGCNA sono riportati nella Figura supplementare 3.
5. Selezione dei biomarcatori
Per stabilire una gerarchia di priorità tra i candidati da sottoporre a ulteriori analisi sperimentali, sono stati applicati in parallelo tre algoritmi di apprendimento automatico¹⁵. È stata utilizzata la regressione con operatore di riduzione e selezione della contrazione assoluta (LASSO) per la selezione delle caratteristiche, la foresta casuale con ntree = 500 per determinare l'ordine di importanza dei geni e la macchina a vettori di supporto con eliminazione ricorsiva delle caratteristiche (SVM-RFE) per identificare il sottoinsieme di caratteristiche con l'errore minore nella validazione incrociata. I risultati completi delle analisi di screening sono riportati nella Figura 5 supplementare.
L'importanza dei geni è stata valutata in base al contributo di ciascun gene all'accuratezza della classificazione, e i geni con rango più elevato sono stati raccolti per generare un insieme preliminare di geni¹⁶. L'algoritmo SVM-RFE ha eseguito la selezione delle caratteristiche e l'ottimizzazione del modello mediante eliminazione iterativa delle caratteristiche, consentendo l'identificazione del sottoinsieme ottimale di caratteristiche per la classificazione¹⁷. I geni comuni ai tre algoritmi sono stati considerati biomarcatori candidati identificati attraverso l'apprendimento automatico.
I CF-DEGs sono stati utilizzati per costruire una rete di interazione proteina-proteina (PPI), e CytoHubba è stato impiegato per classificare i nodi centrali. L'intersezione tra i moduli chiave di WGCNA, i tre risultati dell'apprendimento automatico e il set centrale PPI ha identificato cinque biomarcatori candidati integrati: LCN2, IDO1, CXCL2, NOS2 e CD274. Questi geni sono stati classificati come candidati diagnostici associati alla colite ulcerosa. Successivi esperimenti con inibitori hanno ulteriormente classificato LCN2, IDO1, CXCL2 e NOS2 come marcatori responsivi all'inibizione della ferroptosi e CD274 come marcatore responsivo all'inibizione della cuproptosi. Tali geni non sono stati considerati regolatori primari dimostrati di nessuno dei due percorsi.
6. Sviluppo del modello predittivo e valutazione multimetodo
Per valutare il potenziale clinico traslazionale della firma genica identificata, è stato costruito un nomogramma utilizzando il pacchetto rms per analizzare le prestazioni diagnostiche. Le prestazioni del modello sono state valutate mediante una curva di calibrazione, il test di Hosmer–Lemeshow (HL) e l'indice di concordanza di Harrell (C-index).
La calibrazione del modello è stata valutata utilizzando il test di Hosmer-Lemeshow (HL) e l'errore assoluto medio (MAE). Sono stati stabiliti come soglie per un adeguato adattamento del modello e per l'accuratezza rispettivamente un valore di P > 0,05 e un MAE < 0,1. Sono state generate curve caratteristiche operative del ricevitore (ROC) e calcolata l'area sotto la curva (AUC) per valutare le prestazioni predittive del nomogramma per la colite ulcerosa (UC). Un dataset indipendente è stato utilizzato per valutare la generalizzabilità e la stabilità del modello.
7. Analisi dell'arricchimento di geni (GSEA)
Nel set di addestramento, è stata eseguita un'analisi della correlazione di Spearman tra ciascun biomarcatore e tutti gli altri geni. I risultati della correlazione sono stati ordinati per generare un insieme di geni ordinato. È stata effettuata un'analisi dell'arricchimento degli insiemi genici (GSEA) utilizzando l'insieme genico c2.cp.kegg_medicus.v2025.1.Hs.symbols del Molecular Signatures Database (MSigDB) come riferimento. I percorsi biologici con P < 0,05 sono stati considerati significativamente arricchiti.
8. Analisi del paesaggio immunitario
Lo strumento CIBERSORT e la matrice di firma LM22 sono stati utilizzati per eseguire l'analisi di deconvoluzione dei dati di espressione e stimare le proporzioni di 22 sottotipi di cellule immunitarie. Successivamente è stata effettuata un'analisi di correlazione di Spearman per valutare le associazioni tra i cinque biomarcatori e i sottotipi di cellule immunitarie.
9. Analisi della rete regolatoria
Le informazioni sulla regolazione genica da parte di fattori di trascrizione (TF) e microRNA (miRNA) sono state ottenute dalle basi di dati integrate ChEA3 e TarBase 9.0 attraverso la piattaforma NetworkAnalyst. I miRNA associati ad almeno due nodi sono stati inclusi nella rete di regolazione.
10. RNA-seq a singola cellula
I dataset di singole cellule UC GSE116222, costituiti da tre campioni sani e tre campioni UC, e GSE214695, costituiti da sei campioni sani e sei campioni UC, sono stati integrati. Le cellule sono state mantenute se contenevano da 200 a 5.000 geni rilevati, >1.000 identificatori molecolari unici, < 30% di trascritti mitocondriali, < 3% di trascritti di emoglobina e < 50% di trascritti ribosomiali. La funzione NormalizeData è stata utilizzata per la normalizzazione e Harmony è stato impiegato per la correzione del batch prima del clustering e dell'annotazione.
Dopo il filtraggio, i dati sono stati normalizzati utilizzando la funzione NormalizeData. Harmony è stato utilizzato per minimizzare gli effetti di batch tra i campioni. L'approximazione uniforme della varietà (UMAP) e l'incorporamento stocastico di vicini distribuiti in t (t-SNE) sono stati impiegati per la riduzione della dimensionalità e il raggruppamento cellulare al fine di caratterizzare le strutture locali e globali all'interno della popolazione cellulare. Un totale di 31.712 cellule di alta qualità è stato mantenuto per ulteriori analisi. I tipi cellulari sono stati annotati in base ai geni marcatore canonici.
L'analisi del pseudotempo è stata eseguita utilizzando Monocle 2 per studiare i cambiamenti negli stati delle cellule epiteliali e mieloidi provenienti da tessuti di colite ulcerosa. Sono state valutate le traiettorie di sviluppo cellulare e le differenze tra gruppi, e i modelli di espressione dei biomarcatori lungo la traiettoria del pseudotempo sono stati visualizzati.
CellChat è stato utilizzato per valutare la comunicazione tra cellule epiteliali e sottopopolazioni di cellule mieloidi. I profili di espressione specifici delle sottopopolazioni sono stati integrati per inferire le probabilità di interazione e le principali vie di segnalazione tra cellule epiteliali e mieloidi e per caratterizzare le potenziali reti di regolazione intercellulare.
11. Esperimenti sulle cellule
Le cellule umane di adenocarcinoma colon-rettale Caco-2 sono state coltivate in terreno di Dulbecco modificato di Eagle contenente il 10% di siero fetale bovino, 100 U/mL di penicillina e 100 µg/mL di streptomicina a 37 °C in atmosfera al 5% di CO₂. Le cellule Caco-2 sono state scelte come modello riproducibile di monocapa epiteliale per saggi sulla barriera intestinale e su quelli infiammatori e lesivi; tuttavia, la loro origine maligna limita l'estrapolazione diretta all'epitelio colonnare non maligno. Le informazioni sui passaggi cellulari non erano disponibili nei registri sperimentali archiviati.
Un modello in vitro di infiammazione-danno è stato stabilito esponendo le cellule a 1 µg/mL di lipopolisaccaride (LPS) per 24 h. Le cellule sono state assegnate a un gruppo controllo con dimetilsulfossido come veicolo, a un gruppo con solo LPS, a un gruppo con LPS + 4 µM di ferrostatina-1 o a un gruppo con LPS + 40 µM di tetratiomolibdato, con concentrazioni finali di solvente equivalenti. Il gruppo con solo LPS è stato utilizzato come controllo positivo del danno/modello per i confronti di recupero con inibitori.
- Viabilità cellulare
La viabilità cellulare è stata valutata mediante un saggio di viabilità cellulare secondo le istruzioni del produttore. Dopo un'incubazione a 37 °C per circa 2 h, l'assorbanza a 450 nm è stata misurata utilizzando un lettore di micropiastre.
- Misurazione dei livelli di MDA e Fe²⁺
I livelli intracellulari di malondialdeide (MDA) e Fe2⁺ sono stati misurati utilizzando i rispettivi kit di saggio secondo le istruzioni dei produttori. Le proteine totali sono state quantificate mediante un saggio con acido bicinconinico e ogni valore ottenuto è stato normalizzato rispetto al contenuto proteico del campione corrispondente.
- Saggio con sonda fluorescente per ioni rame
Dopo la rimozione del mezzo di coltura, è stata aggiunta una sonda fluorescente specifica per Cu2⁺ a una concentrazione finale di 1 µM e incubata per 30 min a 37 °C. Le immagini sono state acquisite immediatamente senza lavaggi mediante un microscopio a fluorescenza invertito con lunghezze d'onda di eccitazione e di emissione di 510 e 578 nm, rispettivamente. Le impostazioni di acquisizione sono state mantenute costanti tra i diversi gruppi e la barra di scala delle immagini è stata calibrata a 50 µm.
- RT-qPCR
L'RNA totale è stato isolato utilizzando un reagente di estrazione dell'RNA basato su fenolo–guanidinio. Il DNA complementare è stato sintetizzato da 1 µg di RNA utilizzando un kit per la sintesi del DNA complementare del primo filamento e reagenti a base di SYBR Green sono stati impiegati per la PCR in tempo reale. Ogni reazione è stata eseguita in triplicato tecnico. Le sequenze dei primer sono riportate nella Tabella Supplementare 2 e l'abbondanza relativa dell'RNA messaggero è stata calcolata con il metodo 2⁻ΔΔCt.
12. Analisi statistica
La normalità e l'omogeneità della varianza sono state valutate utilizzando rispettivamente i test di Shapiro–Wilk e di Levene centrato sulla mediana. I dati che soddisfacevano entrambe le assunzioni sono stati analizzati mediante un test t a due code non appaiato o un'analisi della varianza ad un fattore seguita dal test di Tukey per la differenza significativa onesta. Quando una delle due assunzioni non era soddisfatta, è stato utilizzato il test t di Welch o l'analisi della varianza di Welch seguita dai confronti di Games–Howell.
I risultati sono presentati come media ± deviazione standard. I conteggi effettivi dei replicati biologici, i valori P esatti a due code e i controlli delle assunzioni sono riportati nella Tabella Supplementare 3. Il foglio sorgente dell'Figura 13 archiviato contiene un numero minore di replicati di controllo per i pannelli A (n = 4), C (n = 4) e F (n = 2), mentre tutte le altre combinazioni pannello/gruppo contengono n = 6.
Le analisi computazionali sono state eseguite utilizzando R e i flussi di lavoro limma, WGCNA, clusterProfiler, glmnet, randomForest, e1071, pROC, rms, Seurat v4, Harmony, Monocle 2, CellChat e Cytoscape/CytoHubba. Le risorse online incluse sono FerrDB, GEO, STRING, ChEA3, TarBase 9.0, NetworkAnalyst, CIBERSORT/LM22, MSigDB, Bioinformatics.com.cn e Sangerbox 3.0.