Identification des DEG
Deux ensembles de données, GSE78721 et GSE12288, ont été utilisés pour l’analyse bioinformatique dans cette étude. En GSE78721, un total de 2 266 DEG ont été identifiés entre des patients atteints de diabète de type 2 et des échantillons normaux, comprenant 1 209 gènes à la hausse et 1 057 à la baisse (Figure 1A). En GSE12288, 924 DEG ont été vérifiés entre les patients atteints de MCSCV et les témoins, incluant 403 gènes surexprimés et 521 sous-régulés (Figure 1C). Les 50 principaux DEG identifiés dans les deux ensembles de données sont représentés dans les cartes thermiques montrées à la Figure 1B,D. En croisant ces deux ensembles de données, 32 C-DEG ont été identifiés, comprenant 20 gènes à la hausse et 12 à la baisse (Figure 2A, B).
Enrichissement des C-DEG par GO et KEGG
Pour élucider les fonctions biologiques et les voies clés associées aux 32 C-DEG, des analyses GO et KEGG ont été réalisées. Les dix principaux termes liés au processus biologique (BP), au composant cellulaire (CC) et à la fonction moléculaire (MF) parmi les C-DEG sont présentés dans la Figure 3A. Les fonctions de la BP incluent principalement la régulation de l’activation du plasminogène, la régulation de la fibrinolyse, la régulation positive de la différenciation des cardiocytes, la régulation négative de la liaison des protéines et la différenciation des cardioblastes. Les fonctions CC concernent le complexe endopeptidase de type sérine, les complexes protéiques impliqués dans l’adhésion de la matrice cellulaire, le chromosome condensé, la membrane granulaire spécifique et le corps PML. La MF englobe les constituants structurels du cytosquelette, l’activité activatrice de la NADPH oxydase générant des superoxydes, l’activité intramoléculaire de la transferase (phosphotransférases), l’activité exopeptidase de type sérine et la liaison à l’IgG (Figure 3A). En ce qui concerne l’analyse KEGG, les C-DEG étaient principalement impliqués dans des voies critiques, notamment la myocardite virale, la cardiomyopathie du ventricule droit arrythmogène, les cascades de complément et de coagulation, le cancer de la prostate, la cardiomyopathie hypertrophique, la cardiomyopathie dilatée, l’apoptose, le phagosome et le traitement des protéines dans le réticulum endoplasmique (Figure 3B).
Construction de réseaux IPP et criblage génique central
Parmi les 32 C-DEG, trois gènes n’ont montré aucune interaction connue avec d’autres C-DEG dans la base de données STRING et ont donc été exclus de la construction du réseau PPI. La plateforme STRING a été utilisée pour l’analyse PPI des 29 C-DEG restants, révélant un réseau composé de 29 nœuds et 39 arêtes, visualisé à l’aide de Cytoscape (Figure 4A). Par la suite, huit C-DEG de hub (HSP90B1, PLAU, SLPI, TOP3A, NCF4, PRF1, TUBA1C et CS) ont été identifiés en fonction de leur degré de connectivité (Figure 4B).
Différences dans l’infiltration des cellules immunitaires et leur corrélation avec les C-DEG du hub
Cette analyse s’est concentrée exclusivement sur le jeu de données T2DM (GSE78721) afin d’étudier préliminairement les changements dans le microenvironnement immunitaire dans le contexte du DT2DM et leur relation avec les gènes centrals. Des graphiques de violon ont été utilisés pour visualiser la relation entre le jeu de données T2DM et les niveaux d’infiltration de 22 types de cellules immunitaires (Figure 5A ; détails dans le Tableau Supplémentaire 3). L’analyse d’infiltration des cellules immunitaires a indiqué une proportion plus élevée de cellules dendritiques activées dans le groupe T2DM par rapport au groupe témoin (P = 0,049), tandis que la proportion de cellules dendritiques au repos était significativement plus faible dans le groupe T2DM (P = 0,023). De plus, la relation entre les niveaux d’expression des huit C-DEG centraux et les niveaux d’infiltration des 22 différents types de cellules immunitaires a été visualisée à l’aide d’une carte thermique (Figure 5B). Notamment, l’expression de TOP3A a montré une forte corrélation positive avec le niveau d’infiltration des lymphocytes T régulateurs (Tregs) (P < 0,001), tandis que l’expression du SLPI était fortement corrélée positivement avec l’infiltration des lymphocytes T gamma delta (P < 0,001). Inversement, l’expression de NCF4 a montré une forte corrélation négative avec le niveau d’infiltration des cellules NK activées (P < 0,001), l’expression de PRF1 était fortement corrélée négativement avec le niveau d’infiltration des macrophages M2 (P < 0,001), et l’expression de PLAU était fortement négativement liée au niveau d’infiltration des cellules B naïves (P < 0,001).
Réseau candidat de médicaments cibles pour les C-DEG concentrées
Parmi les huit C-DEG du hub identifiés, CS et TOP3A n’étaient pas listés comme cibles médicamenteuses connues dans la base de données DGIdb et ont donc été exclus de l’analyse suivante du réseau candidat-cible. Pour identifier des médicaments ayant un rôle thérapeutique potentiel dans le DT2DM et l’ASCVD, les six C-DEG du hub restants ont été soumis à la base de données DGIdb (détails dans le Tableau Supplémentaire 4). Par la suite, Cytoscape a été utilisé pour construire un réseau candidat-médicament cible illustrant les associations entre les C-DEG centraux et les médicaments ciblés candidats (Figure 6). Ce réseau comprenait six C-DEG centrales, 136 médicaments moléculaires et 137 arêtes. Plus précisément, 82, 5, 3, 4, 1 et 42 médicaments cibles candidats ont été identifiés pour TUBA1C, NCF4, PRF1, SLPI, HSP90B1 et PLAU, respectivement, la majorité étant classés comme agents antinéoplasiques, anti-inflammatoires et antithrombotiques.
Validation de l’expression des C-DEG du hub
En raison de l’absence de CS et TOP3A comme cibles médicamenteuses candidates dans la base de données DGIdb et de leurs fonctions biologiques étant moins directement liées aux stratégies thérapeutiques dans le contexte de la recherche sur le DT2 et l’ASCVD, l’étude s’est concentrée sur les six C-DEG centraux restants (HSP90B1, PLAU, SLPI, NCF4, PRF1 et TUBA1C) pour une validation expérimentale. Pour approfondir les rôles de ces C-DEG hub, des tests qRT-PCR et western blot ont été réalisés afin d’évaluer les niveaux d’expression des gènes centraux entre différents groupes (Figure 7 et Figure 8). Les niveaux d’expression de TUBA1C, NCF4, SLPI, HSP90B1 et PLAU étaient significativement augmentés chez les patients atteints de DT2 et ASCVD par rapport aux groupes témoins sains correspondants, tandis que l’expression de PRF1 était à la baisse. Ces résultats suggèrent que ces gènes fondamentaux pourraient jouer un rôle significatif dans la progression des patients atteints de DT2 et d’ASCVD.
DISPONIBILITÉ DES DONNÉES :Les ensembles de données d’expression génique GSE78721 et GSE12288 utilisés pour l’analyse bioinformatique dans cette étude sont disponibles dans la base de données publique GEO du NCBI. Les données brutes issues de la partie validation expérimentale générées dans cette étude sont fournies dans le Fichier Supplémentaire 1.

Figure 1 : Identification des DEG dans les ensembles de données GSE78721 et GSE12288. (A,C) Les graphiques volcaniques des DEG dans GSE78721 et GSE12288, respectivement. |logFC|≥ 0,5 et P < 0,05 ont été établis comme seuil de dépistage ; Les gènes à la hausse étaient marqués comme rouges, et ceux à la baisse étaient marqués comme verts. (B,D) Cartes thermiques des 50 meilleurs DEG de GSE78721 et GSE12288, respectivement. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 2 : Diagramme de Venn des C-DEG à partir de deux jeux de données. (A) Diagramme de Venn des C-DEG à la hausse provenant de T2DM et ASCVD dans les ensembles de données GSE78721 et GSE12288. (B) Diagramme de Venn des C-DEG à la baisse à partir de T2DM et ASCVD dans les ensembles de données GSE78721 et GSE12288. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 3 : Analyse fonctionnelle et d’enrichissement des voies des C-DEG. Graphiques à barres des analyses d’enrichissement GO (A) et KEGG (B). Les coordonnées verticales indiquent la description de différents items GO ou des voies KEGG, et les coordonnées horizontales indiquent le nombre de C-DEG enrichis. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 4 : Dépistage génique du réseau et central des C-DEG. (A) Le réseau PPI des C-DEG. Le rouge et le bleu représentent respectivement les C-DEG à la hausse et à la baisse. (B) Huit gènes hubs ont été identifiés parmi les C-DEG. Les nœuds représentent un C-DEG, et les arêtes représentent l’interaction entre les C-DEG. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 5 : Analyse de l’infiltration des cellules immunitaires. (A) Différences dans l’infiltration de 22 types de cellules immunitaires entre le groupe T2DM et le groupe témoin. Axe horizontal : sous-types de cellules immunitaires ; Axe vertical : niveau relatif d’infiltration des cellules immunitaires. (B) Carte thermique de corrélation de l’expression de 8 C-DEG hubs et des niveaux d’infiltration de 22 types de cellules immunitaires. Carrés rouges : corrélation positive ; Carrés bleus : corrélation négative. Des couleurs plus foncées indiquent des corrélations plus fortes. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 6 : Réseau cible de médicaments candidats pour six C-DEG centrales. Les nœuds rouges représentent les C-DEG du hub, les nœuds verts représentent les médicaments cibles candidates, et les arêtes entre eux représentent les paires d’interactions. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 7 : Validation de l’expression génique centrale chez les patients atteints de DT2DM et ASCVD par qRT-PCR. Les niveaux d’expression de l’ARNm HSP90B1 (A), PLAU (B), SLPI (C), NCF4 (D), PRF1 (E) et TUBA1C (F) ont été détectés par qRT-PCR. Veuillez cliquer ici pour voir une version agrandie de cette figurine.

Figure 8 : Validation de l’expression génique centrale chez les patients atteints de DT2DM et d’ASCVD par Western blot. Les niveaux d’expression protéique HSP90B1 (A), PLAU (B), SLPI (C), NCF4 (D), PRF1 (E) et TUBA1C (F) ont été détectés par des tests Western blot. Veuillez cliquer ici pour voir une version agrandie de cette figurine.
Tableau supplémentaire 1. Matrice d’expression génique pour le jeu de données GSE78721, utilisé pour l’analyse transcriptomique T2DM dans cette étude. Veuillez cliquer ici pour télécharger ce fichier.
Tableau supplémentaire 2. Matrice d’expression génique pour le jeu de données GSE12288, utilisée pour l’analyse transcriptomique de l’ASCVD dans cette étude. Veuillez cliquer ici pour télécharger ce fichier.
Tableau complémentaire 3 : Analyse de l’infiltration des cellules immunitaires dans le jeu de données GSE78721. Veuillez cliquer ici pour télécharger ce fichier.
Tableau complémentaire 4 : Interactions candidates-cibles pour les gènes communs exprimés différemment. Veuillez cliquer ici pour télécharger ce fichier.
Fichier supplémentaire 1 : Les données brutes issues de la validation expérimentale générées dans cette étude. Veuillez cliquer ici pour télécharger ce fichier.