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 microARN provenant de TCGA-COAD et TCGA-READ, issues du GDC (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 plus élevée de miR-192 que les échantillons normaux (médiane de 16,05 contre 13,46 log2[RPM + 1], Mann-Whitney p = 9,96 x 10-8). La même tendance vers des valeurs plus élevées dans les tumeurs a été observée pour COAD seul (455 tumeurs et 8 échantillons normaux, p = 9,89 x 10-8) et pour READ seul (161 tumeurs et 3 échantillons normaux, p = 1,94 x 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 tumeur-normale = 3,13 log2[RPM + 1], Wilcoxon p = 0,0029). Ces résultats restaient discordants avec le résultat local obtenu par qRT-PCR appariée, qui indiquait 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 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 supplémentaire équilibré 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 comprenait seulement 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 la cohorte combinée COADREAD, la comparaison de Kaplan-Meier avec une découpe médiane n'était pas significative (test 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 (test log-rank p = 0,889 ; Cox p = 0,493) et dans READ (test log-rank p = 0,807 ; Cox p = 0,716). Ainsi, les données ne soutiennent pas une affirmation forte en faveur de 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 pour 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
Des 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 le 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 témoins (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 miARN et d'ARN-séquençage 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 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 représente le tissu normal et le rouge 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 CRC et les tissus appariés. (A) Analyse de séquençage d'ARNmi provenant 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 CRC par rapport aux échantillons normaux. (B) Analyse par qRT-PCR de prélèvements locaux appariés montrant une expression plus faible de la forme mature de miR-192-5p dans les tissus CRC par rapport aux 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-normale 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 le chevauchement entre les prédictions de gènes cibles des 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 géniques prédites ou documentées. Veuillez cliquer ici pour visualiser une version agrandie de cette figure.

Figure 5 : Analyse de criblage des gènes centraux à partir du réseau cible candidat. (A,B) Modules de réseau identifiés par criblage 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 du criblage basé sur les modules et sur le degré. 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 de fluorescence de cellules HT29 transfectées avec un mimétique de miR-192-5p. (C) Image de 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 l'incision à 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 visualiser une version agrandie de cette figure.

Figure 9 : Corrélation entre l'expression de miR-192 dans TCGA et les scores de signatures liées au système immunitaire dans les échantillons tumoraux de CCR. L'analyse a révélé des associations immunitaires hétérogènes, avec plusieurs corrélations négatives faibles 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 de 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 de 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 d'expression et de 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 initiales de Ct obtenues par qRT-PCR, les données de cytométrie en flux ou de transfection, les mesures de largeur des plaies, ainsi que 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.