Cette étude a utilisé des ensembles de données publiques et dé-identifiées (TCGA, GTEx et CGGA) et a établi des lignées cellulaires commerciales. Aucun participant humain nouvellement recruté, aucune information identifiable sur les patients ou des spécimens dérivés des patients n’a été impliqué. Toutes les analyses ont été réalisées conformément aux directives institutionnelles et internationales pertinentes pour l’utilisation des données publiques. Par conséquent, aucune approbation éthique institutionnelle supplémentaire ni consentement éclairé n’a été requis pour cette étude.
1. Flux de travail bioinformatique et prétraitement des données
Des profils d’expression génique et des annotations cliniques correspondantes ont été obtenus à partir de ressources publiques, notamment The Cancer Genome Atlas (TCGA), le projet Genotype-Tissue Expression (GTEx) et le Chinese Glioma Genome Atlas (CGGA). Les matrices d’expression et les métadonnées cliniques ont été importées dans R et harmonisées selon les identifiants de l’échantillon. Les entrées dupliquées ou non appariées ont été retirées, et des échantillons dépourvus d’informations essentielles sur la survie ont été exclus des analyses de survie. Les valeurs d’expression de MRC2 et les variables clinicopathologiques correspondantes ont été extraites pour les analyses en aval. Pour les comparaisons tumorales-normales, l’expression de MRC2 a été comparée entre types de cancer à partir de données TCGA uniquement lorsque des tissus normaux appariés étaient disponibles, et en utilisant des données intégrées TCGA–GTEx lorsque des contrôles normaux supplémentaires étaient nécessaires. Pour les analyses axées sur le gliome, les patients des cohortes TCGA-LGG et CGGA ont été stratifiés en groupes à haute et faible expression selon le niveau médian d’expression de la cohorte MRC2. L’analyse de survie de Kaplan-Meier, l’analyse des caractéristiques opérationnelles du récepteur (ROC) dépendantes du temps et la régression des risques proportionnels de Cox ont ensuite été réalisées pour évaluer la pertinence pronostique de MRC2. Les modèles de Cox incluaient les covariables disponibles telles que l’âge, le sexe, le grade de l’OMS et d’autres variables clinicopathologiques, selon la complétude des données au sein de chaque cohorte. Des analyses de sous-groupes ont été réalisées selon des variables telles que l’âge, le sous-type histologique, le grade OMS, le statut de mutation IDH, le statut de codélétion 1p/19q et les informations liées au traitement, lorsque disponibles. Les figures ont été générées en utilisant des flux de travail de tracé statistique standard dans R et GraphPad Prism.
2. Analyses pan-cancéreuses et analyses de survie
L’expression de MRC2 a été comparée entre tumeurs et tissus normaux entre types de cancer à l’aide des ensembles de données TCGA et GTEx. Les différences ont été visualisées à l’aide de boîtes et résumées à l’aide de diagrammes radar. Des analyses de régression des risques proportionnels de Cox ont été réalisées pour évaluer les associations entre l’expression de MRC2 et la survie globale (OS), la survie spécifique à la maladie (DSS) et l’intervalle sans progression (PFI) entre les types de cancer. Les rapports de risque et les niveaux de signification correspondants ont été visualisés à l’aide de parcelles forestières et de cartes thermiques.
3. Validation pronostique et analyses de sous-groupes dans les cohortes de gliomes
Dans la cohorte de gliome de grade inférieur (TCGA-LGG), les patients ont été stratifiés en groupes à haute et faible expression MRC2. Des analyses de survie de Kaplan–Meier ont été réalisées pour comparer OS, DSS et PFI entre groupes. Des courbes de caractéristiques de fonctionnement du récepteur (ROC) dépendantes du temps ont été générées pour évaluer la performance prédictive. Des analyses de régression de Cox univariées et multivariées ont été réalisées pour évaluer la valeur pronostique indépendante de MRC2. La validation externe a été réalisée à l’aide du jeu de données de gliomes CGGA, incluant l’analyse de survie de Kaplan–Meier, l’analyse ROC et la régression de Cox. Des analyses de sous-groupes ont été réalisées pour évaluer l’expression de MRC2 à travers des variables clinicopathologiques, notamment l’âge, le grade tumoral, le sous-type histologique, le statut de mutation IDH, le statut de codélétion 1p/19q et les variables liées au traitement.
4. Expression de MRC1 dans les gliomes mutants IDH versus les gliomes sauvages de type IDH
L’association entre l’expression de MRC1 et le statut de mutation IDH a été évaluée à l’aide de la plateforme web GlioVis et du jeu de données TCGA GBMLGG. Les cas ont été regroupés selon le statut annoté de l’IDH en gliomes mutants IDH et de type sauvage IDH. L’expression de l’ARNm MRC1 a été interrogée dans le jeu de données sélectionné, et les niveaux d’expression entre les deux groupes ont été comparés à l’aide de la fonction d’analyse statistique intégrée à la plateforme. Les distributions d’expressions par groupes étaient affichées telles qu’elles étaient générées par la plateforme, et la signification statistique était définie comme P < 0,05.
5. Culture cellulaire
Les lignées cellulaires de gliome humain et les astrocytes humains normaux ont été prélevés à partir de dépôts cellulaires établis. Les cellules ont été cultivées dans le milieu d’Eagle modifié (DMEM) de Dulbecco, complétés avec un sérum fœtal bovin (FBS) à 10 % et maintenus à 37 °C dans un incubateur humidifié contenant 5 % de CO₂ et environ 95 % d’humidité relative. Les cellules U87 et U251 ont été sélectionnées pour des expériences de perte de fonction car le dépistage de base sur le panel cellulaire de gliome testé a montré que ces deux lignées cellulaires exprimaient des niveaux relativement plus élevés de MRC2 et étaient donc adaptées aux analyses fonctionnelles basées sur le knockdown.
6. Réduction lentivirale de MRC2
Des ARN en épingle courte (shRNA) ciblant MRC2 et un témoin non ciblant ont été clonés en vecteurs lentiviraux. Les particules lentivirales étaient produites par co-transfection de cellules HEK293T avec des plasmides de transfert et en emballant les plasmides à l’aide d’un réactif de transfection standard. Les surnatants viraux étaient collectés à 48 h et 72 h après la transfection, centrifugés à 2 100 × g pendant 5 minutes pour éliminer les cellules d’emballage et les débris cellulaires, puis passés à travers un filtre de 0,45 μm avant utilisation ou stockage immédiat à −80 °C 19. Les surnadants contenant des virus étaient utilisés immédiatement ou aliquotés et stockés à −80 °C. Les cycles répétés de gel-dégel ont été évités. Les cellules de gliome ont été infectées par des particules lentivirales en présence de polybrène. Après l’infection, des cellules ont été sélectionnées avec de la puromycine afin d’établir des lignées cellulaires de knockdown stables. Les séquences utilisées étaient les suivantes : sh-MRC2-1 : 5′-GAAATGAATGAGCAGCAAGAA-3′ ; sh-MRC2-2 : 5′-CCGGTATTGCTATAAGGTGTT-3′ ; sh-NC : 5′-TTCTCCGAACGTGTCACGT-3′.
7. Validation du knockdown de MRC2 par RT-qPCR et immunoblotting
L’ARN total a été extrait à l’aide d’un réactif TRIzol, et l’ADNc a été synthétisé à l’aide d’un kit de transcription inverse selon les instructions du fabricant. La PCR quantitative a été réalisée en utilisant la chimie verte SYBR dans les conditions de cycle suivantes : dénaturation initiale à 95 °C pendant 30 s, suivie de 40 cycles à 95 °C pendant 5 s et 60 °C pendant 30 s, avec une analyse de la courbe de dissociation effectuée à la fin de l’amplification. Les séquences d’amorces étaient les suivantes : MRC2 avant, 5′-GGCAAGGAAGAAGTGCGTGT-3′ ; MRC2 en marche arrière, 5′-CTTTGGTGTTGCTGCTGCTT-3′ ; GAPDH avant, 5′-TCCACCCATGGCAAATTCC-3′ ; et GAPDH inverse, 5′-TCGCCCCACTTGATTTTGG-3′. L’expression relative de l’ARNm a été calculée à l’aide de la méthode comparative Ct (2^-ΔΔCt) avec le GAPDH comme contrôleinterne 20. Un immunoblotting a été réalisé en parallèle pour confirmer la réduction de MRC2 au niveau protéique. Pour la RT-qPCR, trois réplications biologiques indépendantes ont été analysées, et chaque échantillon a été analysé en triplique technique technique.
8. Essai de prolifération cellulaire
La prolifération cellulaire a été évaluée à l’aide du test CCK-8. Les cellules ont été ensemencées en plaques de 96 puits à une densité appropriée et cultivées pour les périodes indiquées. À chaque moment, un réactif CCK-8 était ajouté à chaque puits et incubé selon les instructions du fabricant, après quoi l’absorption à 450 nm était mesurée à l’aide d’un lecteur de microplaques. Au moins trois réplications biologiques indépendantes ont été réalisées.
9. Essai de formation de la colonie
Les cellules ont été ensemées à faible densité en plaques à six puits et cultivées pendant 10 à 14 jours. Les colonies ont été fixées avec du paraformoldéhyde à 4 % pendant 15 minutes à température ambiante et colorées à 0,1 % de violet cristallin pendant 15 minutes. Les colonies ont été lavées avec du sérum physiologique tamponné phosphate (PBS), séchées à l’air et comptées manuellement ou à l’aide du logiciel d’analysed’image 21. Au moins trois réplications biologiques indépendantes ont été réalisées.
10. Tests de migration et d’invasion cellulaires
La migration et l’invasion cellulaire ont été évaluées à l’aide de tests à insertion membranaireporeuse 22. Pour les tests d’invasion, les inserts supérieurs étaient pré-recouverts de gel de matrice extracellulaire ; Pour les tests de migration, des inserts non revêtus étaient utilisés. Des cellules suspendues dans un milieu sans sérum étaient ajoutées dans la chambre supérieure, et un milieu contenant du sérum était placé dans la chambre inférieure comme chimiottractant. Après l’incubation, les cellules non migratrices de la surface supérieure de la membrane ont été délicatement retirées à l’aide d’un coton-tige. Les cellules ayant migré ou envahi la surface inférieure ont été fixées avec 4 % de paraformoldéhyde pendant 15 minutes, colorées au violet cristallin à 0,1 % pendant 15 minutes, rincées avec du sérum salubre tamponné de phosphate froid, séchées à l’air et imagées au microscope inversé à un × 200 de grossissement. Les cellules ont été comptées dans cinq champs microscopiques sélectionnés au hasard par insert. Au moins trois réplications biologiques indépendantes ont été réalisées.
11. Analyse du Western blot
Les cellules ont été lysées dans un tampon RIPA complété avec des inhibiteurs de protéase, et la concentration en protéines a été déterminée à l’aide d’un test BCA. Des quantités égales de protéines totales ont été séparées par SDS-PAGE et transférées dans les membranes de PVDF. Les membranes ont été bloquées dans 5 % de lait écrémé pendant 1 heure à température ambiante et incubées avec des anticorps primaires toute la nuit à 4 °C. Les anticorps primaires et dilutions de travail suivants ont été utilisés pour le western blotting : anti-MRC2 (1:5 000), anti-CDK4 (1:1 000), anti-CDK6 (1:2 000), anti-β-caténine (1:5 000), anti-E-cadherine (1:10 000), anti-N-cadherine (1:3 000) et anti-β-actine (1:5 000). Après lavage, les membranes ont été incubées avec des anticorps secondaires anti-IgG de chèvre conjugués à la peroxydase de raifort, ou avec des anticorps secondaires anti-lapins de chèvre (1:5 000) pendant 1 h à température ambiante. Les signaux étaient visualisés à l’aide de la chimiluminescence, et les intensités des bandes étaient quantifiées et normalisées en β-actine. Au moins trois réplications biologiques indépendantes ont été réalisées.
12. Analyse statistique
Des analyses statistiques ont été réalisées à l’aide de R et de GraphPad Prism. Les courbes de survie ont été comparées à l’aide du test log-rank. Des modèles de régression des allocations proportionnelles de Cox ont été utilisés pour les analyses univariées et multivariées. Pour les expériences cellulaires, les données sont présentées en moyenne ± SD. Avant le test paramétrique, la normalité était évaluée à l’aide du test de Shapiro–Wilk. Pour les données distribuées normalement, les différences entre deux groupes ont été évaluées à l’aide du test t de Student, et des comparaisons entre plusieurs groupes ont été réalisées à l’aide d’ANOVA unidirectionnelle. Si les hypothèses de normalité n’étaient pas satisfaites, les tests non paramétriques correspondants étaient utilisés. Les valeurs p ont été ajustées pour plusieurs tests à l’aide de la procédure de taux de fausse découverte de Benjamini–Hochberg 23,24. Un P à deux faces < 0,05 était considéré comme statistiquement significatif.