Analyse d'expression dans des jeux de données publics et des tissus cliniques appariés
Le profil d'expression de miR-192 a été réévalué à l'aide des données de séquençage des ARNmi provenant de GDC, issues des cohortes TCGA-COAD et TCGA-READ (Figure 1). Dans ces fichiers, l'entrée au niveau de la base de données était annotée sous la forme hsa-mir-192 plutôt que directement comme l'espèce mature hsa-miR-192-5p ; par conséquent, ce résultat provenant de la base de données publique a été interprété comme un signal TCGA de miR-192. Dans la cohorte COADREAD combinée, 616 tumeurs primitives et 11 échantillons normaux de tissu solide ont été analysés. Les échantillons tumoraux présentaient une expression de miR-192 plus élevée que les échantillons normaux (médiane de 16,05 contre 13,46 log2[RPM + 1], test de Mann-Whitney p = 9,96 × 10-8). La même tendance vers une expression plus élevée dans les tumeurs a été observée séparément dans COAD (455 tumeurs et 8 échantillons normaux, p = 9,89 × 10-8) et dans READ (161 tumeurs et 3 échantillons normaux, p = 1,94 × 10-5). Une analyse appariée TCGA a également montré une expression tumorale plus élevée dans 11 cas appariés (différence médiane entre tumeur et tissu normal = 3,13 log2[RPM + 1], test de Wilcoxon p = 0,0029). Ces résultats sont en désaccord avec le résultat local obtenu par qRT-PCR sur des échantillons appariés, qui indique une réduction de l'expression de la forme mature miR-192-5p dans les tissus de CCR (Figure 2).
Performance diagnostique et associations clinicopathologiques
L'analyse ROC basée sur TCGA a été répétée en utilisant la définition de discrimination entre échantillons tumoraux et normaux décrite précédemment (Figure 3). La valeur d'AUC était de 0,968 pour le groupe combiné COADREAD, de 0,947 pour COAD seul et de 0,994 pour READ seul. Un échantillonnage rééquilibré supplémentaire a donné des AUC médianes de 0,983 pour COADREAD, de 0,984 pour COAD et de 1,000 pour READ ; toutefois, les intervalles de confiance à 95 % étaient larges pour COADREAD (0,868–1,000) et pour COAD (0,781–1,000), et READ ne comprenait que trois échantillons normaux. Une analyse laissant de côté un échantillon normal à la fois a également montré que l'estimation ROC était sensible à la petite taille du groupe normal. Par conséquent, le résultat ROC a été conservé uniquement comme preuve exploratoire de séparation entre tumeur et tissu normal, et non comme preuve suffisante pour un biomarqueur diagnostique clinique. Une analyse d'association clinicopathologique basée sur l'expression médiane de miR-192 dans TCGA a révélé des associations significatives avec le stade N (p = 0,019) et l'âge (p = 0,004), tandis que le stade T (p = 0,110) et le stade M (p = 0,932) n'étaient pas significatifs (Table 1).
Analyse de survie
L'analyse de la survie globale n'a pas révélé d'association pronostique significative pour l'expression de miR-192 du TCGA dans le CCR (Figure 3). Dans le groupe COADREAD combiné, la comparaison de Kaplan-Meier avec une découpe médiane n'était pas significative (log-rank p = 0,857), et la régression de Cox univariée n'a également montré aucune association significative (HR = 0,934 par unité d'expression log2, p = 0,374). La même conclusion a été observée dans COAD (log-rank p = 0,889 ; Cox p = 0,493) et dans READ (log-rank p = 0,807 ; Cox p = 0,716). Ainsi, les données ne soutiennent pas une affirmation forte concernant miR-192 comme biomarqueur pronostique dans le CCR.
Gènes cibles candidats et enrichissement fonctionnel
Une analyse des gènes cibles candidats a été réalisée en intégrant les prédictions de miRWalk avec les gènes différentiellement exprimés provenant des jeux de données GEO GSE89076 et GSE156355. Ce criblage a identifié 93 gènes candidats destinés à des analyses ultérieures de réseau et d'enrichissement (Figure 4). L'analyse des interactions protéiques a retenu KIF20A, TPX2, CDCA5, CCNB1, CDK1, PLP1, NRXN1, GRIK3 et KIF5C comme gènes centraux exploratoires (Figure 5). Une analyse de sensibilité basée sur les données d'ARN-séquençage du projet TCGA a montré que KIF20A, TPX2, CDCA5, CCNB1 et CDK1 étaient surexprimés dans le tissu tumoral par rapport au tissu normal, tandis que PLP1, NRXN1, GRIK3 et KIF5C étaient sous-exprimés dans le tissu tumoral. L'analyse de corrélation avec l'expression de miR-192 dans TCGA a donné des résultats mitigés : CDCA5, CCNB1 et CDK1 présentaient des corrélations positives, tandis que PLP1, NRXN1 et KIF5C présentaient des corrélations négatives. L'analyse d'enrichissement a mis en évidence des termes représentatifs liés au métabolisme des médicaments par le cytochrome P450, au métabolisme des xénobiotiques par le cytochrome P450, au métabolisme du rétinol, aux termes associés à la synapse glutamatergique et à l'organisation de la présynapse (Figure 6). Ces résultats soutiennent la pertinence biologique du réseau de gènes centraux, mais n'établissent pas de ciblage direct par miR-192-5p.
Dosage de la migration cellulaire
Les cellules HT29 ont été transfectées avec des mimiques de miR-192 afin d'évaluer les effets fonctionnels de la surexpression de miR-192. La microscopie à fluorescence a confirmé une transfection réussie dans les groupes transfectés avec le mimique de miR-192-5p et dans le groupe contrôle négatif (Figure 7). Dans le test de cicatrisation de la plaie, le taux de migration des cellules HT29 transfectées avec miR-192 était inférieur à celui des cellules HT29 non transfectées et des cellules transfectées avec le contrôle négatif, à la fois à 24 h et à 48 h (Figure 8). À 24 h, le taux de migration était de 9,75 +/- 2,43 % dans le groupe transfecté avec miR-192, contre 15,69 +/- 3,47 % dans le groupe non transfecté et 15,45 +/- 3,92 % dans le groupe contrôle négatif (P = 0,010). À 48 h, le taux de migration était de 20,04 +/- 2,54 % dans le groupe transfecté avec miR-192, contre 26,84 +/- 7,65 % et 27,98 +/- 6,78 % respectivement dans les deux groupes contrôles (P = 0,021). Ces résultats suggèrent qu'une surexpression de miR-192 supprime la migration des cellules HT29 in vitro.
Analyse de corrélation des signatures immunitaires
L'association entre l'expression de miR-192 dans TCGA et les signatures liées à l'immunité a été réévaluée à l'aide de données d'ARN-séquençage provenant de tumeurs primitives (Figure 9). Un total de 425 échantillons tumoraux dotés de données appariées de séquençage d'ARN et de miARN a été inclus dans l'analyse de corrélation des scores de signatures. L'expression de miR-192 présentait de faibles corrélations négatives avec plusieurs signatures de cellules immunitaires, notamment les cellules Tfh (rho = -0,205, p = 2,15e-05), les cellules dendritiques activées (rho = -0,186, p = 0,00011), les cellules Treg (rho = -0,185, p = 0,00013), les cellules dendritiques plasmacytoïdes (rho = -0,169, p = 0,00046), les macrophages (rho = -0,149, p = 0,0020), les cellules Tcm (rho = -0,136, p = 0,0050), les cellules NK CD56bright (rho = -0,132, p = 0,0065) et les cellules Tem (rho = -0,129, p = 0,0078). Des corrélations positives faibles ont été observées pour les signatures des éosinophiles (rho = 0,139, p = 0,0040) et des cellules Th17 (rho = 0,121, p = 0,0129). Ainsi, l'association immunitaire de miR-192 dans le CCR semble hétérogène et ne devrait pas être décrite comme uniformément positive parmi les infiltrats immunitaires.
Disponibilité des données
Les jeux de données publics analysés dans cette étude sont disponibles sur le portail GDC/TCGA et dans la base de données Gene Expression Omnibus (GEO) sous les numéros d'accès GSE89076 et GSE156355. Les analyses in silico et les réanalyses générées au cours de cette étude sont fournies dans le Tableau supplémentaire 1. Les données expérimentales brutes à l'origine des analyses de qRT-PCR et de cicatrisation, incluant les valeurs brutes de Ct en qRT-PCR, les mesures de largeur de blessure, les images microscopiques originales et les données expérimentales associées, sont fournies dans le Dossier supplémentaire 1.

Figure 1 : Aperçu de l'expression pancancéreuse du miR-192 dans le TCGA. L'axe des abscisses indique les abréviations des tumeurs, et l'axe des ordonnées indique l'expression relative du miR-192. Le bleu indique le tissu normal, et le rouge indique le tissu tumoral. La signification statistique est indiquée comme suit : * = p < 0,05, ** = p < 0,01, *** = p < 0,001 ; NS = non significatif. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 2 : expression de miR-192/miR-192-5p dans les jeux de données de CCR et les tissus appariés. (A) Analyse de séquençage d'ARNmi à partir de TCGA/GDC montrant un signal plus élevé d'hsa-mir-192 au niveau de la base de données dans les échantillons tumoraux de CCR par rapport aux échantillons normaux. (B) Analyse par qRT-PCR d'échantillons appariés de tissus locaux montrant une expression plus faible de la forme mature de miR-192-5p dans les tissus de CCR comparée à celle des tissus non cancéreux adjacents. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 3 : Analyses diagnostiques exploratoires et analyses de survie basées sur l'expression de hsa-mir-192 dans TCGA. (A–C) Courbes de survie de Kaplan-Meier. L'analyse de survie n'a pas révélé d'association pronostique significative. (D) Courbe ROC pour la discrimination tumeur-normal dans la cohorte TCGA. La valeur AUC doit être interprétée avec prudence en raison du faible nombre d'échantillons normaux. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 4 : Analyse de criblage des gènes cibles candidats et analyse des interactions protéiques. (A) Diagramme de Venn montrant les recouvrements entre les prédictions de gènes cibles de miARN et les gènes différentiellement exprimés provenant des séries GSE89076 et GSE156355. (B) Réseau d'interactions protéine-protéine des 93 gènes candidats. Les nœuds représentent les gènes, et les arêtes représentent des associations de gènes prédites ou validées. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 5 : Analyse des gènes centraux à partir du réseau cible candidat. (A,B) Modules du réseau identifiés par l'analyse des modules. (C) Classement des gènes centraux selon la connectivité du réseau. (D) Les cibles centrales candidates sont conservées après intégration de l'analyse basée sur les modules et sur le degré de connexion. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 6 : Termes d'enrichissement GO et KEGG sélectionnés pour les gènes candidats et les modules significatifs du réseau. Le graphique en bulles résume les fonctions biologiques et les voies métaboliques enrichies représentatives. Les résultats d'enrichissement sont exploratoires et nécessitent une validation expérimentale. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 7 : Images de transfection des cellules HT29. (A) Cellules HT29 représentatives. (B) Image en fluorescence de cellules HT29 transfectées avec un mimétique de miR-192-5p. (C) Image en fluorescence de cellules HT29 transfectées avec un contrôle négatif con238. Échelle = 100 µm. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 8 : Essai de cicatrisation des plaies dans les cellules HT29. (A) Images représentatives de rayures à 0 h, 24 h et 48 h chez les cellules HT29 non transfectées, les cellules HT29 transfectées avec un mimétique de miR-192-5p et les cellules HT29 transfectées avec con238. (B) Quantification des taux de cicatrisation à 24 h et 48 h, exprimés en pourcentage. Échelle = 100 µm. Veuillez cliquer ici pour consulter une version agrandie de cette figure.

Figure 9 : Corrélation entre l'expression de miR-192 selon TCGA et les scores de signatures liées à l'immunité dans les échantillons tumoraux de CCR. L'analyse a révélé des associations immunitaires hétérogènes, avec plusieurs faibles corrélations négatives et, principalement pour les signatures des éosinophiles et des cellules Th17, des corrélations positives. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.
| Caractéristiques | Expression faible de TCGA miR-192 | Expression élevée de TCGA miR-192 | Valeur p | Méthode |
| n | 308 | 308 | | |
| Stade T, n (%) | 0,11 | Chisq.test |
| T1 | 9 (1,5 %) | 11 (1,8 %) | | |
| T2 | 44 (7,2 %) | 59 (9,6 %) | | |
| T3 | 209 (34,1 %) | 211 (34,4 %) | | |
| T4 | 43 (7 %) | 27 (4,4 %) | | |
| Stade N, n (%) | 0,019 | Chisq.test |
| N0 | 156 (25,5 %) | 190 (31 %) | | |
| N1 | 88 (14,4 %) | 63 (10,3 %) | | |
| N2 | 61 (10 %) | 54 (8,8 %) | | |
| Stade M, n (%) | 0,932 | Chisq.test |
| M0 | 214 (39,9 %) | 234 (43,7 %) | | |
| M1 | 41 (7,6 %) | 47 (8,8 %) | | |
| Âge, médiane (EIQ) | 66 (56, 75) | 69 (60, 77) | 0,004 | Wilcoxon |
Tableau 1 : Association entre les groupes d'expression de miR-192 du TCGA et les caractéristiques clinicopathologiques dans le CRC. Les groupes d'expression élevée et faible ont été définis d'après le niveau médian d'expression de miR-192 du TCGA. Dans le tableau récapitulatif, le stade N et l'âge étaient significatifs, tandis que le stade T et le stade M n'étaient pas significatifs.
Tableau supplémentaire 1 : Résultats de la réanalyse in silico générés pour cette étude. Le classeur comprend des comparaisons d'expression TCGA/GDC, des analyses de sensibilité ROC équilibrées, des analyses de type « laisser-un-normal-de-côté », des analyses d'association clinicopathologiques, des analyses de survie, des résumés de l'expression et de la corrélation des gènes centraux, ainsi que des résultats de corrélation des signatures immunitaires.Veuillez cliquer ici pour télécharger ce fichier.
Dossier supplémentaire 1 : Données expérimentales brutes soutenant les expériences de qRT-PCR, de transfection cellulaire et de cicatrisation des plaies. Ces fichiers comprennent les valeurs de Ct initiales de la qRT-PCR, les données de cytométrie en flux/transfection, les mesures de largeur des plaies et des images microscopiques représentatives utilisées pour les analyses présentées dans le manuscrit.Veuillez cliquer ici pour télécharger ce fichier.