Toutes les procédures animales ont été menées conformément aux directives institutionnelles et approuvées avant le début de l’étude par le Comité éthique de l’Université de médecine traditionnelle chinoise de Chengdu (Approbation n° 2022-126). Seuls des ensembles de données GEO et GWAS disponibles publiquement ont été utilisés pour les analyses de bioinformatique ; par conséquent, aucun consentement éclairé supplémentaire ni approbation éthique n’a été requis pour la composante computationnelle de cette étude. Toutes les procédures impliquant des produits chimiques dangereux, des échantillons biologiques et des déchets animaux ont été réalisées conformément aux réglementations institutionnelles sur la sécurité des laboratoires et aux procédures institutionnelles pour l’élimination des déchets dangereux.
Acquisition et prétraitement des ensembles de données GEO
Les ensembles de données transcriptomiques humaines liés à la colite ulcéreuse (CU) ont été extraits de la base de données GEO en utilisant le mot-clé « colite ulcéreuse ». Après avoir sélectionné des ensembles de données de profilage d’expression avec des groupes UC et témoins sains clairement définis et une taille d’échantillon suffisante, GSE8746620 et GSE9241521, tous deux générés sur la plateforme GPL13158, ont été sélectionnés comme cohorte de formation. Ces ensembles de données comprenaient respectivement 87 échantillons de CU et 21 témoins sains, ainsi que 53 échantillons d’UC et 21 témoins sains. Pour la validation externe,GSE87473 22 (106 échantillons de coline ulcéreuse et 21 témoins sains) et GSE1687923 (48 échantillons de coline ultérieuse et 12 témoins sains ; GPL570) étaient incluses. Des fichiers bruts ou en matrice série étaient téléchargés et importés dans R pour le prétraitement. Les identifiants de sonde étaient convertis en symboles génétiques selon les fichiers d’annotation correspondants de la plateforme. Lorsque plusieurs sondes étaient mappées sur le même symbole génique, la valeur moyenne d’expression était conservée. Les échantillons dépourvus d’informations claires sur les groupes ont été exclus avant l’analyse en aval. Les deux ensembles de données d’entraînement ont été fusionnés par des symboles génétiques communs. Les valeurs d’expression étaient transformées log2 lorsque nécessaire, normalisées à l’aide d’une procédure standard de normalisation entre tableaux, et corrigées pour les effets batch à l’aide d’une méthode d’ajustement batch après avoir défini l’origine du jeu de données comme variable batch. L’analyse des composants principaux et l’inspection par boxplot étaient utilisées avant et après la normalisation/correction par lots comme points de contrôle qualité. Le prétraitement réussi a été défini par un meilleur chevauchement des distributions d’échantillons entre les ensembles de données et l’atténuation du clustering piloté par les ensembles de données. Un total de 902 gènes liés à la PANoptose ont été compilés à partir de la littérature publiée et utilisés comme ensemble de gènes de référence pour une analyse intégrative.
Identification des gènes exprimés différemment et construction du réseau WGCNA
La matrice d’entraînement fusionnée a été analysée dans R à l’aide d’un flux de travail d’expression différentielle. Gènes avec|log 2 multiplication| ≥ 0,585 et P ajusté < 0,05 ont été définis comme des gènes différenciellement exprimés (DEG). Des graphiques volcaniques et des cartes thermiques ont été générés comme résultats intermédiaires pour confirmer que les critères de filtrage produisaient des différences d’expression biologiquement interprétables entre les groupes UC et témoins. Pour l’analyse pondérée du réseau de co-expression génique (WGCNA), les gènes ont été classés par variance entre les échantillons, et les 25 % les gènes les plus variables ont été conservés comme entrées. Un arbre de regroupement d’échantillons a d’abord été examiné pour identifier des échantillons potentiels d’anomalies ; aucun cas aberrant évident n’a été conservé pour la construction du réseau, sauf si la suppression était justifiée par des critères de contrôle qualité. Une puissance de seuil souple a été sélectionnée en fonction de l’indice d’ajustement topologique sans échelle, et la plus petite puissance obtenant un réseau approximativement sans échelle a été utilisée pour construire la matrice d’adjacence. La matrice d’adjacence a ensuite été transformée en une matrice de chevauchement topologique, et un regroupement hiérarchique a été réalisé en utilisant la dissimilarité basée sur le TOM. Les modules ont été identifiés grâce à la coupe dynamique de l’arbre avec une taille minimale de module de 100 gènes. Des modules étroitement apparentés ont été fusionnés lorsque leur corrélation d’eigengènes dépassait le critère de fusion prédéfini. Les propres gènes des modules, la signification des modules et les valeurs de signification génique ont été calculés pour identifier les modules les plus fortement associés au phénotype de la CU. Une étape WGCNA réussie a été définie par une séparation stable des modules, des tailles biologiquement plausibles, et un ou plusieurs modules montrant une corrélation claire avec le statut de la maladie.
Identification des cibles candidates dans GQD
Des composés candidats et cibles correspondant aux quatre composants herbacés de la Décoction Gegen Qinlian (GQD) — Pueraria lobata, Scutellaria baicalensis, Coptis chinensis et Glycyrrhiza uralenses — ont été collectés à partir de TCMSP24 et BATMAN-TCM25. Les seuils de dépistage ont été fixés comme suit : biodisponibilité orale ≥ 30 % et ressemblance au médicament ≥ 0,18 pour les dossiers dérivés d’une seule base de données ; score de confiance ≥ 0,84 et ajusté P < 0,05 pour les enregistrements dérivés d’une autre base de données. Des composés supplémentaires ont été complétés à partir d’une base de données de médecine à base de plantes lorsque les dossiers pertinents n’ont pas été capturés lors de la recherche initiale. Les chaînes canoniques SMILES de composés retenus ont été soumises à une plateforme d’évaluation ADME en ligne. Les composés ont été conservés lorsqu’ils ont montré une forte absorption gastro-intestinale et ont satisfait au moins deux règles de ressemblance de médicaments chez Lipinski, Ghose, Veber, Egan et Muegge. Les cibles présumées ont ensuite été prédites à l’aide d’une plateforme de prédiction cible avec une probabilité > 0,1 et standardisées aux symboles officiels des gènes à l’aide d’une base de données d’annotation protéique. Les composés et cibles dupliqués ont été supprimés après la fusion de la base de données. Le résultat attendu de cette étape était un ensemble de données GQD composé-cible non redondant, adapté à l’analyse de chevauchement.
Identification de l’action des gènes qui se chevauchent et performance de l’analyse fonctionnelle
Les cibles thérapeutiques potentielles étaient définies comme l’intersection entre cibles liées à la GQD, les gènes liés à la CU identifiés à partir de l’expression différentielle et du WGCNA, ainsi que l’ensemble de gènes liés à la PANoptose. Les relations de chevauchement ont été visualisées à l’aide d’un diagramme de Venn. Les gènes intersectants ont ensuite été soumis à une base de données d’interactions protéine-protéine (IPP), et les paires d’interactions avec un score de confiance ≥ 0,4 ont été conservées. Le réseau PPI résultant a été exporté et visualisé dans une plateforme d’analyse de réseau. Une étape PPI réussie était définie par la présence d’un réseau connecté contenant la majorité des gènes intersectant, plutôt que par des nœuds isolés seulement. L’analyse d’enrichissement fonctionnel, incluant l’ontologie génique (GO) et l’analyse de l’Encyclopédie des Gènes et Génomes de Kyoto (KEGG), a été réalisée en R avec q < 0,05 comme seuil de signification. Les conditions ou voies ne respectant pas ce seuil ont été exclues. Les processus biologiques enrichis par le haut et les voies de signalisation ont été résumés afin d’inférer le mécanisme potentiel par lequel la GQD pourrait réguler l’UC. Des graphiques en barres et des diagrammes à bulles ont été générés comme résultats intermédiaires pour vérifier si les résultats d’enrichissement étaient dominés par des voies liées à l’inflammation, à l’immunité ou à la lésion épithéliale, comme on pouvait s’y attendre dans le contexte biologique.
Analyse des schémas d’expression, de la localisation chromosomique et des corrélations
Les niveaux d’expression des gènes partagés ont été comparés entre les groupes CU et témoins à l’aide d’un flux de travail standard de visualisation d’expression différentielle dans R. Des boxplots et des cartes thermiques ont été générés, et P < 0,05 a été jugé statistiquement significatif. Seuls les gènes présentant une direction d’expression interprétable et un comportement cohérent au sein de la cohorte d’entraînement ont été conservés pour une évaluation mécaniste ultérieure. La localisation chromosomique de certains gènes centraux a été cartographiée après formatage des fichiers d’annotation génétique dans la structure d’entrée requise, et les positions génomiques ont été visualisées à l’aide d’un flux de travail de tracé chromosomique circulaire. Les matrices de corrélation entre les gènes de noyau ont ensuite été calculées par analyse de corrélation par paires. Le résultat attendu de cette étape était l’identification de paires de gènes présentant des schémas d’expression coordonnés pouvant soutenir une régulation partagée ou une interaction fonctionnelle.
Construction et évaluation de modèles d’apprentissage automatique
Le développement du modèle diagnostique a été réalisé à l’aide de la cohorte d’entraînement fusionnée (GSE87466 + GSE92415), tandis que GSE16879 et GSE87473 ont été utilisés comme ensembles de données externes indépendants de validation. Douze algorithmes d’apprentissage automatique — Lasso, Ridge, Stepglm, XGBoost, Random Forest, Enet, plsRglm, GBM, NaiveBayes, LDA, glmBoost et SVM — ont été utilisés pour générer 113 combinaisons demodèles, dont 26. Une validation croisée dix fois plus grande a été appliquée au sein de la cohorte de formation. Dans chaque pliage, l’entraînement du modèle, la sélection des caractéristiques et l’estimation des performances ont été effectués en utilisant la même logique de partitionnement d’entraînement/test afin d’éviter la fuite d’informations. La performance du modèle a été principalement évaluée par la zone sous la courbe caractéristique de fonctionnement du récepteur (AUC). Le modèle optimal a été défini comme le modèle avec le taux d’AUC moyen le plus élevé dans les ensembles de données d’entraînement et de validation, plutôt que le seul AUC le plus élevé d’entraînement. Les modèles montrant un AUC d’entraînement très élevé mais de faibles performances de validation externe n’ont pas été conservés. Des courbes ROC ont été générées pour chaque modèle conservé, et un nomogramme a été construit à partir des biomarqueurs finaux sélectionnés. Une étape de modélisation réussie a été définie par une performance stable en validation croisée et la préservation de la capacité discriminative dans les ensembles de données externes.
Performance de l’expression génique GSVA et du profil central
L’analyse de variation des ensembles de gènes (GSVA) a été réalisée pour comparer l’enrichissement des voies entre les groupes à haute expression et à faible expression définis par la valeur médiane d’expression de chaque gène central. Un P ajusté < 0,05 a été utilisé comme seuil pour les différences de voie significatives. Les voies ne respectant pas ce seuil ont été exclues de toute interprétation. Des cartes thermiques et des graphiques de corrélation ont été générés pour illustrer l’activité des voies et les schémas de co-expression. Le succès intermédiaire de cette étape a été indiqué par des profils d’enrichissement compatibles avec une dysrégulation inflammatoire, immunitaire ou métabolique pertinente pour la CU.
Analyse du paysage immunitaire et construction du réseau ceRNA
La composition des cellules immunitaires a été estimée à l’aide d’un algorithme de déconvolution avec 1 000 permutations. Les échantillons ne respectant pas le critère de signification interne de la déconvolution ont été exclus de l’analyse immunitaire ultérieure. Pour compléter cette analyse, l’analyse d’enrichissement de l’ensemble de gènes à échantillon unique (ssGSEA) a également été utilisée pour évaluer les schémas d’infiltration immunitaire. Les corrélations entre les fractions des cellules immunitaires et l’expression des gènes de base ont ensuite été calculées et visualisées. Pour la construction d’un réseau de ceRNA, des miARN ciblant les ARNm de base ont été identifiés à partir de bases de données d’interactions publiques, et les RNncRNA candidats interagissant avec ces miARN ont ensuite été dépistés. Seules les interactions supportées par des enregistrements de prédiction ou d’annotation de la base de données étaient conservées. Le réseau régulateur mRNA-miARN-lncRNA a été visualisé sur une plateforme d’analyse de réseau. Le résultat attendu de cette étape était un réseau structuré de ceRNA avec des régulateurs en amont biologiquement plausibles liés aux gènes centraux.
Performance de l’analyse à cellule unique
Le jeu de données à une seule cellule GSE214695 a été importé dans R et converti en objets à une seule cellule en utilisant un flux de travail standard d’analyse à une seule cellule. Le contrôle de qualité a été effectué en filtrant les cellules selon des métriques standard, incluant les cellules à faible caractéristique, les gouttelettes vides potentielles et les cellules de mauvaise qualité des transcrits. Les cellules ayant réussi le contrôle qualité ont été normalisées, et les 2 000 gènes les plus variables ont été sélectionnés pour l’analyse des composantes principales. Les principaux composants contribuant à la variance biologique significative ont été conservés pour le regroupement non supervisé, et l’UMAP a été utilisé pour la visualisation bidimensionnelle. Les gènes marqueurs de cluster ont été identifiés par analyse d’expression différentielle entre groupes. L’annotation par type de cellule était ensuite réalisée par classification basée sur la référence, et les annotations étaient vérifiées avec l’expression canonique des marqueurs lorsque cela était possible. Une étape unicellulaire réussie a été définie par la séparation claire des principales populations cellulaires sur le graphique UMAP, des gènes marqueurs interprétables pour chaquegroupe 27, et une annotation cohérente des principaux types cellulaires immunitaires et stromales pertinents pour la CU.
Performance de l’analyse de randomisation mendélienne
Une randomisation mendélienne (RM) sur deux échantillons a été réalisée afin d’évaluer les possibles relations causales entre l’expression génique centrale et la susceptibilité à la CU. Des variables instrumentales ont été extraites des ensembles de données de locus quantitatifs de traits d’expression, et les statistiques résumées d’études d’association à l’échelle du génome UC ont été utilisées comme ensemble de données de résultats. Les SNP ont été dépistés selon des critères standards de contrôle qualité de la RM, incluant leur pertinence pour le trait d’exposition et la suppression des variantes ambiguës ou dupliquées lorsque cela était exigé par l’harmonisation. La méthode pondérée par l’inverse de la variance a été utilisée comme estimateur causal principal. Le test Q de Cochran a été appliqué pour évaluer l’hétérogénéité. Les analyses MR-Egger et MR-PRESSO ont été utilisées comme analyses de sensibilité pour évaluer et, si nécessaire, corriger la pléiotropie horizontale. Une étape d’IRM réussie a été définie par une harmonisation valide de l’instrument, l’absence d’hétérogénéité majeure ou depléiotropie 28, et des estimations directionnellement cohérentes entre les méthodes IRM complémentaires.
Performance du docking moléculaire de composés candidats
Les structures tridimensionnelles des protéines cibles ont été obtenues à partir de la Protein Data Bank, et les structures de ligands ont été obtenues à partir d’une base de données publique de petites molécules. Avant l’amarrage, les protéines et les ligands étaient prétraités en retirant les molécules d’eau lorsque cela était approprié, en ajoutant des hydrogènes, en définissant les types d’atomes et en convertissant les formats de fichiers nécessaires à l’amarrage. L’amarrage a ensuite été effectué selon un flux de travail d’amarrage moléculaire, et les énergies de liaison ont été calculées pour chaque paire ligand-cible. Les affinités de liaison inférieures à -5 kcal/mol ont été interprétées comme une liaison favorable, tandis que des valeurs inférieures à -7 kcal/mol ont été interprétées comme une liaison relativement forte. Des cartes thermiques d’amarrage ont été tracées pour comparer les schémas de liaison globales entre cibles. Les conformations les mieux classées ont été visualisées en trois dimensions pour inspecter les liaisons hydrogène, les interactions hydrophobes et la proximité des résidus. La réussite de cette étape était définie par des fichiers de sortie d’amarrage stables, des poses plausibles de ligands dans la région de liaison cible, et des énergies de liaison compatibles avec la priorisation des candidats.
Performance de la validation expérimentale in vivo
Préparation au GQD
Le GQD a été préparé à partir de quatre plantes médicinales traditionnelles chinoises dans les proportions brutes suivantes : 24 g de Puerariae Lobatae Radix (Gegen), 9 g de Scutellariae Radix (Huangqin), 9 g de Coptidis Rhizoma (Huanglian), et 6 g de Glycyrrhizae Radix et Rhizoma (Gancao). Les détails des matériaux botaniques et de leurs sources sont fournis dans le Tableau des matériaux. Les herbes étaient broyées en poudre fine, bien mélangées, puis trempées dans 390 mL d’eau distillée pendant 1 heure. Le mélange a été décocté deux fois. Pour la première décoction, les herbes ont été bouillies à feu vif puis mijotées à feu doux pendant 30 minutes. Après filtration à travers la gaze, les résidus étaient rébouillonnés avec un volume dix fois plus grand que l’eau. Les deux filtrats ont été combinés et concentrés dans un bain-marie jusqu’à une concentration finale de drogue brute de 5 g/mL. L’extrait était ensuite filtré, refroidi et conservé à 4 °C jusqu’à utilisation.
Établissement du modèle souris UC induit par le DSS
Dix-huit souris mâles BALB/c (âgés de 8 semaines) ont été utilisés dans cette étude, et toutes les procédures animales ont été approuvées par le Comité d’éthique de l’Université de médecine traditionnelle chinoise de Chengdu (Approbation n° 2022-126). Les détails sur la source animale sont fournis dans le tableau des matériaux. Les souris étaient hébergées dans des conditions spécifiques sans agents pathogènes à 20–22 °C et 55 % d’humidité lors d’un cycle clair/obscur de 12 heures. Après une période d’acclimatation d’une semaine, les animaux ont été divisés aléatoirement en trois groupes (n = 6 par groupe) : NC, UC et GQD. Le groupe NC a reçu de l’eau potable stérile tout au long de l’expérience, tandis que les groupes UC et GQD ont reçu 1,5 % (p/v) de dextran sulfate de sodium (MW 36-50 kDa) dans de l’eau potable pendant 5 jours consécutifs, suivis de 2 jours d’eau ordinaire. Ce cycle a été répété trois fois pour établir une colite chronique. Pendant la période d’intervention, les groupes NC et CU ont été traités deux fois par jour avec 200 μL de carboxyméthyl cellulose à 0,5 %, tandis que le groupe GQD a reçu la GQD par gavage oral deux fois par jour.
La prescription clinique de GQD comprenait 24 g de Gegen, 9 g de Huangqin, 9 g de Huanglien et 6 g de Gancao, correspondant à une quantité totale de drogue brute de 48 g/jour pour les adultes. En supposant un poids corporel standard adulte de 60 kg, la dose équivalente adulte était de 0,8 g/kg/jour. Sur la base du coefficient de conversion humain-souris conventionnel de 1:12, la dose équivalente pour la souris était de 9,6 g/kg/jour. Selon le ratio prédéfini de doses faible, moyenne et élevée de 1:2:4, les doses calculées étaient respectivement de 9,6, 19,2 et 38,4 g/kg/jour. Compte tenu de la gravité et de la persistance de l’inflammation du côlon dans le modèle répété de colite induite par le DSS, le schéma à dose élevée a été choisi afin d’assurer une intervention pharmacologique adéquate lors de la validation in vivo. Pour soutenir une exposition pharmacologique suffisante pendant la phase inflammatoire active, la dose quotidienne totale de 38,4 g/kg/jour a été divisée en deux administrations égales. Tout au long de l’expérience, le poids corporel, la consistance des selles et les saignements rectaux ont été surveillés quotidiennement. À la fin du protocole, les souris ont été euthanasiées et des tissus du côlon ont été prélevés pour des analyses ultérieures.
Évaluation du poids corporel et de l’indice d’activité de la maladie
Le poids corporel était enregistré quotidiennement et exprimé en pourcentage du poids corporel initial au jour 0. La gravité de la maladie a été évaluée à l’aide de l’indice d’activité de la maladie (DAI), qui comprenait trois paramètres : perte de poids corporel, consistance des selles et saignements rectaux. Chaque paramètre a été noté indépendamment sur une échelle de 0 à 4 selon la gravité des manifestations cliniques, les scores plus élevés indiquant une activité de maladie plus sévère. Le DAI final a été calculé comme la moyenne des trois scores individuels, donnant un score total allant de 0 à 4. En résumé, un score de 0 n’indiquait aucune anomalie, tandis que des scores de 1 à 4 reflétaient une perte de poids progressivement plus importante, une consistance de selles plus molles et des saignements rectaux plus sévères. Les critères de notation détaillés sont présentés dans le tableau 1.
Performance de la coloration H&E et calcul de l’indice d’activité histologique
Pour l’évaluation histologique, les tissus du côlon ont été fixés dans 4 % de paraformaldéhyde, déshydratés par de l’éthanol gradué, intégrés dans la paraffine et coupés à une épaisseur de 4 μm. Les coupes ont été déparaffinées, réhydratées, colorées à l’hématoxyline, contre-colorées à l’éosine, déshydratées, nettoyées et montées. Des images histologiques ont été observées à l’aide d’un microscope en champ clair et scannées à l’aide d’un système de balayage par lames. Les informations spécifiques sur l’équipement sont listées dans le tableau des matériaux. La lésion histopathologique a été évaluée selon trois paramètres, notamment la déplétion des cellules du calice, les dommages architecturaux de la crypte et l’infiltration des cellules inflammatoires. Chaque paramètre a été noté selon le degré de lésion tissulaire, et la somme des trois scores de paramètres a été définie comme l’indice d’activité histologique (HAI), allant de 0 à 10, comme résumé dans le tableau 2. Pour la quantification des cellules gobelet, cinq champs non chevauchants sélectionnés au hasard par section ont été analysés à un grossissement de 400×. Le nombre de cellules goblet par champ était compté manuellement de manière aveugle, et le comptage moyen par champ était calculé pour chaque section colonique représentative.
Performance de la coloration immunohistochimique
Après la déparaffinisation et la réhydratation, les sections tissulaires ont subi une prélèvement antigénique avec 0,1 % de trypsine à 37 °C pendant 30 minutes, suivies d’une incubation avec du peroxyde d’hydrogène à 3 % pour bloquer l’activité endogène de la peroxydase. Les sections ont ensuite été incubées toute la nuit à 4 °C avec des anticorps primaires de lapin contre STAT3 (1:200, HUABIO, ET1605-45), TIMP1 (1:1000, Proteintech, 26847-1-AP), SPHK2 (1:500, Proteintech, 17096-1-AP) et HSPA5 (1:1000, Cell Signaling Technology, #3177), dilués à 1 % de BSA. Les détails des anticorps et de leurs sources sont fournis dans le Tableau des matériaux. Après lavage avec PBS, les sections ont été incubées avec l’anticorps secondaire correspondant (1:2000) à température ambiante pendant 1,5 h. L’immunoréactivité a été visualisée à l’aide d’un réactif de détection chromogénique, et les noyaux ont été contre-colorés à l’hématoxyline. L’analyse quantitative des images a été réalisée à l’aide de logiciels d’analyse d’images listés dans le tableau des matériaux.
Performance de l’analyse qRT-PCR
L’ARN total a été extrait à l’aide d’un réactif d’extraction d’ARN à base de phénol et transcrit inversement à l’aide d’un kit de synthèse d’ADNc de première brin. La transcription inversée a été réalisée à 25 °C pendant 10 minutes, 55 °C pendant 15 minutes et 85 °C pendant 5 minutes. La qRT-PCR a ensuite été réalisée à l’aide d’un système de PCR en temps réel et d’un mélange SYBR Green qPCR. Les conditions d’amplification étaient les suivantes : dénaturation initiale à 95 °C pendant 60 s, suivie de 40 cycles à 95 °C pendant 15 s, 60 °C pendant 15 s, et 72 °C pendant 45 s. Les séquences d’amorces ont été conçues à l’aide d’un logiciel de conception d’amorce et sont listées dans le Tableau 3. Le GAPDH a été utilisé comme référence interne, et les niveaux relatifs d’expression des ARNm ont été calculés à l’aide de la méthode 2^-ΔΔCt . Chaque échantillon a été analysé en triple exemplaire, et chaque expérience a été répétée au moins trois fois. Les détails des réactifs, instruments et logiciels sont fournis dans le tableau des matériaux.
Performance de l’analyse du western blot
Les tissus du côlon ont été homogénéisés dans un tampon de lyse contenant un inhibiteur de protéase sur glace et lisés soigneusement pendant 2 heures. Les concentrations en protéines ont été déterminées à l’aide d’un kit de test protéique. Des quantités égales de protéines ont été mélangées avec un tampon de chargement, dénaturées avant l’électrophorèse, séparées par 10 % de SDS-PAGE, puis transférées sur les membranes PVDF à 100 V pendant 1 heure. Après un blocage avec 5 % de lait écrémé pendant 1 h à température ambiante, les membranes ont été incubées toute la nuit à 4 °C avec des anticorps primaires de lapin contre STAT3 (1:200, HUABIO, ET1605-45), TIMP1 (1:1000, Proteintech, 26847-1-AP), SPHK2 (1:500, Proteintech, 17096-1-AP), HSPA5 (1:1000, Cell Signaling Technology, #3177) et GAPDH à la dilution de travail correspondante. Après lavage avec TBST (20 mM Tris, 150 mM NaCl, 0,1 % Tween-20 ; pH 7,4), les membranes ont été incubées avec l’anticorps secondaire correspondant (1:5000) à température ambiante pendant 1 h. Les bandes protéiques ont été visualisées à l’aide d’un réactif de détection chimiluminescent. Les intensités des bandes ont été analysées à l’aide d’un logiciel d’analyse d’images, et les niveaux d’expression des protéines cibles ont été normalisés en GAPDH. Les poids moléculaires attendus des protéines détectées étaient d’environ 88 kDa pour STAT3, 21 kDa pour TIMP1, 28 kDa pour SPHK2, 70 kDa pour HSPA5 et 36 kDa pour GAPDH. Les détails des anticorps, membranes, réactifs et logiciels sont fournis dans le Tableau des matériaux. Les noms spécifiques des produits, les fabricants, les numéros de catalogue, les modèles d’instruments et les détails logiciels sont fournis dans le tableau des matériaux, tandis que des descriptions génériques sont utilisées dans le texte principal pour améliorer la neutralité et l’applicabilité générale.
Performance de l’analyse statistique
Les données sont présentées en moyenne ± écart-type. Les tests statistiques étaient réalisés à l’aide de logiciels statistiques. L’ANOVA unidirectionnelle suivie du test post hoc de Tukey a été utilisée pour des comparaisons multi-groupes. Une valeur P ≤ 0,05 était considérée comme statistiquement significative.