Cette étude a utilisé des ensembles de données accessibles au public et des lignées cellulaires humaines disponibles commercialement. Aucun participant humain, tissu humain nouvellement collecté ou animal vertébré n'a été impliqué ; par conséquent, aucune approbation par un comité d'éthique institutionnel n'était requise. Les outils de recherche, logiciels et ressources en ligne utilisés dans le protocole sont répertoriés dans le Tableau des matériaux.
1. Traitement des données pour l'analyse d'expression différentielle
Pour établir les cohortes de découverte et de validation transcriptomiques, la base de données Gene Expression Omnibus (GEO) a été interrogée à l’aide du terme « ulcerative colitis ». Quatre jeux de données — GSE87466, GSE92415, GSE107499 et GSE75214 — ont été utilisés comme ensemble d’apprentissage, tandis que GSE47908 a été réservé à la validation externe. Les jeux de données d’apprentissage ont été normalisés et corrigés par lot à l’aide de ComBat dans Sangerbox 3.0, tandis que le jeu de données de validation a été traité indépendamment.
Un total de 583 gènes liés à la ferroptose (FRG) a été obtenu à partir de FerrDB et complété par des gènes rapportés dans la littérature afin de réduire le biais d'une source unique12. Un ensemble de 96 gènes liés à la cuproptose (CRG) a été compilé à partir d'études publiées12,13,14. L'expression différentielle a été analysée à l'aide de limma avec des seuils de |variation du log₂| > 1 et une valeur P ajustée selon Benjamini–Hochberg < 0,05. Les graphiques diagnostiques de la normalisation, de la correction par lot, de l'expression différentielle et du regroupement sont fournis dans la figure supplémentaire 1. La composition des échantillons, les plateformes et les attributions des cohortes sont indiquées dans le tableau supplémentaire 1.
2. Identification des gènes différentiellement exprimés dans la FC au cours de la RC
Une analyse de corrélation de Pearson a été réalisée en utilisant des seuils de |r| > 0,5 et P < 0,05 afin d'identifier les gènes liés à la ferroptose (FRGs) et les gènes liés au cuproptose (CRGs) coexprimés. L'ensemble génique obtenu a été croisé avec les gènes différentiellement exprimés (DEGs) afin d'identifier les gènes différentiellement exprimés coexprimés impliqués dans le cuproptose–ferroptose (CF-DEGs).
3. Analyse d'enrichissement
Des analyses d'enrichissement de l'ontologie génique (GO) et de l'Encyclopédie de Kyoto des gènes et des génomes (KEGG) ont été réalisées à l'aide de clusterProfiler et de la plateforme Bioinformatics.com.cn. Un seuil de significativité défini par une valeur P ajustée selon Benjamini–Hochberg < 0,05 a été utilisé. Les graphiques d'enrichissement sont fournis dans la Figure supplémentaire 2.
4. Analyse de réseau de co-expression génique pondéré (WGCNA)
La structure modulaire sous-jacente aux profils d'expression génique a été analysée à l'aide de l'analyse de réseau de co-expression génique pondérée (WGCNA). Des échantillons anormaux ont été identifiés et exclus par regroupement hiérarchique combiné à la fonction goodSamplesGenes. Une puissance de seuil doux (β) appropriée a ensuite été déterminée afin de satisfaire le critère de topologie sans échelle pour la construction du réseau.
Les données d'expression traitées ont été utilisées pour construire une matrice de superposition topologique (TOM) destinée à l'analyse de la co-expression génique. À l'aide de l'approche de découpage dynamique de l'arbre, les gènes présentant des profils d'expression fortement cohérents ont été regroupés en modules géniques, chacun contenant au moins 100 gènes. Les associations entre les eigengènes des modules (ME) et les phénotypes ont été analysées afin d'identifier les modules liés au processus de la maladie et d'extraire les jeux de gènes correspondants. Les diagnostics de dépistage WGCNA sont fournis dans la figure supplémentaire 3.
5. Sélection des biomarqueurs
Pour hiérarchiser les candidats en vue d'un suivi expérimental, trois algorithmes d'apprentissage automatique ont été appliqués en parallèle¹⁵. Une régression par l'opérateur de sélection et de réduction absolue (LASSO) a été utilisée pour la sélection des caractéristiques, une forêt aléatoire avec ntree = 500 a permis de classer l'importance des gènes, et une machine à vecteurs de support avec élimination récursive des caractéristiques (SVM-RFE) a été employée pour identifier le sous-ensemble de caractéristiques présentant l'erreur de validation croisée la plus faible. Le diagnostic complet du criblage est fourni dans la figure supplémentaire 5.
L'importance des gènes a été évaluée en fonction de la contribution de chaque gène à la précision de la classification, et les gènes les mieux classés ont été regroupés afin de constituer un ensemble préliminaire de gènes¹⁶. L'algorithme SVM-RFE a effectué un criblage des caractéristiques et une optimisation du modèle par élimination itérative des caractéristiques, permettant ainsi d'identifier le sous-ensemble optimal de caractéristiques pour la classification¹⁷. Les gènes communs aux trois algorithmes ont été considérés comme des biomarqueurs candidats identifiés par apprentissage automatique.
Les gènes différentiellement exprimés associés à la CF ont été utilisés pour construire un réseau d'interactions protéine–protéine (PPI), et CytoHubba a permis de classer les nœuds centraux. L'intersection des modules clés de WGCNA, des trois sorties d'apprentissage automatique et de l'ensemble central du réseau PPI a permis d'identifier cinq biomarqueurs candidats intégrés : LCN2, IDO1, CXCL2, NOS2 et CD274. Ces gènes ont été classés comme des candidats diagnostiques associés à la colite ulcéreuse. Des expériences ultérieures avec des inhibiteurs ont ensuite classé LCN2, IDO1, CXCL2 et NOS2 comme des marqueurs sensibles à l'inhibition de la ferroptose, et CD274 comme un marqueur sensible à l'inhibition de la cuproptose. Ces gènes n'ont pas été considérés comme des régulateurs en amont prouvés de l'une ou l'autre voie.
6. Développement du modèle prédictif et évaluation multiméthode
Afin d'évaluer le potentiel clinique de translation de la signature génique identifiée, un nomogramme a été construit à l'aide du package rms pour évaluer la performance diagnostique. La performance du modèle a été évaluée à l'aide d'une courbe de calibration, du test de Hosmer–Lemeshow (HL) et de l'indice de concordance de Harrell (C-index).
L'étalonnage du modèle a été évalué à l'aide du test de Hosmer-Lemeshow (HL) et de l'erreur absolue moyenne (MAE). Une valeur P > 0,05 et une MAE < 0,1 ont été établies comme seuils pour un ajustement et une précision du modèle acceptables, respectivement. Des courbes caractéristiques de fonctionnement du récepteur (ROC) ont été générées, et l'aire sous la courbe (AUC) a été calculée afin d'évaluer les performances prédictives du nomogramme pour la CU. Un jeu de données indépendant a été utilisé pour évaluer la généralisabilité et la stabilité du modèle.
7. Analyse d'enrichissement de l'ensemble de gènes (GSEA)
Dans le jeu d'apprentissage, une analyse de corrélation de Spearman a été réalisée entre chaque biomarqueur et tous les autres gènes. Les résultats de corrélation ont été classés afin de générer un ensemble de gènes ordonné. Une analyse d'enrichissement des ensembles de gènes (GSEA) a été effectuée en utilisant l'ensemble de gènes c2.cp.kegg_medicus.v2025.1.Hs.symbols de la base de données Molecular Signatures Database (MSigDB) comme référence. Les voies biologiques ayant une valeur P < 0,05 ont été considérées comme significativement enrichies.
8. Analyse du paysage immunitaire
L'outil CIBERSORT et la matrice de signature LM22 ont été utilisés pour effectuer une analyse de déconvolution des données d'expression et estimer les proportions des 22 sous-types de cellules immunitaires. Une analyse de corrélation de Spearman a ensuite été réalisée pour évaluer les associations entre les cinq biomarqueurs et les sous-types de cellules immunitaires.
9. Analyse du réseau régulatoire
Les informations relatives à la régulation des gènes par les facteurs de transcription (TF) et les microARN (miARN) ont été obtenues à partir des bases de données intégrées ChEA3 et TarBase 9.0 via la plateforme NetworkAnalyst. Les miARN associés à au moins deux nœuds ont été inclus dans le réseau régulateur.
10. Séquençage ARN unicellulaire
Les jeux de données unicellulaires de la maladie de Crohn GSE116222, comprenant trois échantillons sains et trois échantillons atteints de la maladie de Crohn, et GSE214695, comprenant six échantillons sains et six échantillons atteints de la maladie de Crohn, ont été intégrés. Les cellules ont été conservées si elles contenaient entre 200 et 5 000 gènes détectés, >1 000 identifiants moléculaires uniques, < 30 % des transcrits mitochondriaux, < 3 % de transcrits d'hémoglobine, et < 50 % des transcrits ribosomaux. La fonction NormalizeData a été utilisée pour la normalisation, et Harmony a été utilisée pour la correction par lot avant le regroupement et l'annotation.
Après filtration, les données ont été normalisées à l'aide de la fonction NormalizeData. Harmony a été utilisée pour minimiser les effets de lot entre les échantillons. L'approximation uniforme des variétés et la projection (UMAP) ainsi que l'embedding stochastique t-distribué des voisins les plus proches (t-SNE) ont été utilisées pour la réduction de dimensionnalité et le regroupement cellulaire afin de caractériser les structures locales et globales au sein de la population cellulaire. Un total de 31 712 cellules de haute qualité a été conservé pour des analyses supplémentaires. Les types cellulaires ont été annotés selon les gènes marqueurs canoniques des cellules.
Une analyse de pseudotemporel a été réalisée à l'aide de Monocle 2 afin d'étudier les changements dans les états des cellules épithéliales et myéloïdes provenant de tissus atteints de rectocolite ulcéreuse. Les trajectoires développementales cellulaires et les différences entre groupes ont été évaluées, et les profils d'expression des biomarqueurs le long de la trajectoire pseudotemporelle ont été visualisés.
CellChat a été utilisé pour évaluer la communication cellule-cellule entre les sous-populations épithéliales et myéloïdes. Les profils d'expression spécifiques aux sous-populations ont été intégrés afin d'inférer les probabilités d'interaction et les principales voies de signalisation entre les cellules épithéliales et myéloïdes, et de caractériser leurs réseaux régulateurs intercellulaires potentiels.
11. Expériences sur les cellules
Les cellules humaines de carcinome colorectal Caco-2 ont été cultivées dans un milieu de Dulbecco modifié (Eagle) contenant 10 % de sérum bovin fœtal, 100 U/mL de pénicilline et 100 µg/mL de streptomycine, à 37 °C en atmosphère à 5 % de CO₂. Les cellules Caco-2 ont été choisies comme modèle reproductible de monolame épithéliale pour les tests de barrière intestinale et de lésion inflammatoire ; toutefois, leur origine cancéreuse limite l’extrapolation directe à l’épithélium colique non cancéreux. Les informations concernant les passages cellulaires n’étaient pas disponibles dans les dossiers expérimentaux archivés.
Un modèle in vitro d'inflammation-lésion a été établi en exposant les cellules à 1 µg/mL de lipopolysaccharide (LPS) pendant 24 h. Les cellules ont été réparties entre un groupe témoin recevant le véhicule diméthylsulfoxyde, un groupe LPS seul, un groupe LPS + 4 µM de ferrostatine-1 ou un groupe LPS + 40 µM de tétrathiomolybdate, avec des concentrations finales de solvant appariées. Le groupe LPS seul a servi de témoin positif de lésion/modèle pour les comparaisons de récupération par inhibition.
- Viabilité cellulaire
La viabilité cellulaire a été évaluée à l’aide d’un test de viabilité cellulaire conformément aux instructions du fabricant. Après une incubation à 37 °C pendant environ 2 h, l’absorbance à 450 nm a été mesurée à l’aide d’un lecteur de microplaques.
- Mesure des niveaux de MDA et de Fe²⁺
Les niveaux intracellulaires de malondialdéhyde (MDA) et de Fe2⁺ ont été mesurés à l’aide des kits d’analyse correspondants, conformément aux instructions des fabricants. La quantité totale de protéines a été déterminée par un dosage à l’acide bichinconinique, et chaque résultat a été normalisé par rapport à la teneur en protéines de l’échantillon correspondant.
- Dosage au moyen d’une sonde fluorescente spécifique des ions cuivre
Après élimination du milieu de culture, une sonde fluorescente spécifique du Cu2⁺ a été ajoutée à une concentration finale de 1 µM et incubée pendant 30 min à 37 °C. Les images ont été acquises immédiatement sans lavage à l’aide d’un microscope à fluorescence inversé, avec des longueurs d’onde d’excitation et d’émission de 510 et 578 nm, respectivement. Les paramètres d’acquisition ont été maintenus constants entre les groupes, et l’échelle de la barre de mesure des images a été calibrée à 50 µm.
- RT-qPCR
L’ARN total a été isolé à l’aide d’un réactif d’extraction d’ARN à base de phénol–guanidinium. L’ADN complémentaire a été synthétisé à partir de 1 µg d’ARN à l’aide d’un kit de synthèse du brin complémentaire, et des réactifs basés sur le vert SYBR ont été utilisés pour la PCR en temps réel. Chaque réaction a été réalisée en triple exemplaire technique. Les séquences des amorces sont indiquées dans le Tableau 2 supplémentaire, et l’abondance relative des ARN messagers a été calculée selon la méthode 2⁻ΔΔCt.
12. Analyse statistique
La normalité et l'homogénéité des variances ont été évaluées à l'aide des tests de Shapiro–Wilk et de Levene centré sur la médiane, respectivement. Les données satisfaisant aux deux hypothèses ont été analysées à l'aide d'un test t non apparié bilatéral ou d'une analyse de la variance à un facteur suivie du test de différence honnêtement significative de Tukey. Lorsqu'une des hypothèses n'était pas respectée, le test t de Welch ou l'analyse de la variance de Welch suivie de comparaisons de Games–Howell a été utilisé.
Les résultats sont présentés sous forme de moyenne ± écart type. Les nombres réels de réplicats biologiques, les valeurs P exactes bilatérales et les vérifications des hypothèses sont indiqués dans le Tableau supplémentaire 3. La feuille source archivée de la Figure 13 contient moins de réplicats témoins pour les panneaux A (n = 4), C (n = 4) et F (n = 2), tandis que toutes les autres combinaisons de panneaux/groupes contiennent n = 6.
Les analyses computationnelles ont été réalisées à l'aide de R et des pipelines limma, WGCNA, clusterProfiler, glmnet, randomForest, e1071, pROC, rms, Seurat v4, Harmony, Monocle 2, CellChat, et Cytoscape/CytoHubba. Les ressources en ligne comprenaient FerrDB, GEO, STRING, ChEA3, TarBase 9.0, NetworkAnalyst, CIBERSORT/LM22, MSigDB, Bioinformatics.com.cn et Sangerbox 3.0.