Identificazione dei DEG
Due dataset, GSE78721 e GSE12288, sono stati utilizzati per l'analisi bioinformatica in questo studio. Nel GSE78721, sono stati identificati in totale 2.266 DEG tra pazienti con DM 2 e campioni normali, comprendenti 1.209 geni con aumento e 1.057 geni sottoregolati (Figura 1A). Nel GSE12288, sono stati verificati 924 DEG tra pazienti con ASCVD e controlli, inclusi 403 geni sovraespressi e 521 sottoregolati (Figura 1C). I primi 50 DEG identificati in entrambi i dataset sono rappresentati nelle mappe di calore mostrate nella Figura 1B,D. Attraverso l'intersezione di questi due dataset, sono stati identificati 32 C-DEG, composti da 20 geni aumentati e 12 sottoregolati (Figura 2A,B).
Arricchimento GO e KEGG dei C-DEG
Per chiarire le funzioni biologiche e le vie chiave associate ai 32 C-DEG, sono state effettuate analisi GO e KEGG. I dieci termini principali relativi al processo biologico (BP), al componente cellulare (CC) e alla funzione molecolare (MF) tra i C-DEG sono presentati nella Figura 3A. Le funzioni BP includono prevalentemente la regolazione dell'attivazione del plasminogeno, la regolazione della fibrinolisi, la regolazione positiva della differenziazione cardiocita, la regolazione negativa del legame proteico e la differenziazione cardioblastica. Le funzioni CC riguardano il complesso endopeptidasi di tipo serina, i complessi proteici coinvolti nell'adesione cellulare-matriciale, il cromosoma condensato, la membrana granulosa specifica e il corpo PML. La MF comprende i costituenti strutturali del citoscheletro, l'attività attivatrice della NADPH ossidasi generatrice di superossidi, l'attività intramolecolare della transferasi (fosfotransferasi), l'attività di esopeptidasi di tipo serina e il legame all'IgG (Figura 3A). In termini di analisi KEGG, i C-DEG sono stati principalmente coinvolti in vie critiche, tra cui miocardite virale, cardiomiopatia ventricolare destra aritmogenica, cascate di complemento e coagulazione, cancro alla prostata, miocardiopatia ipertrofica, cardiomiopatia dilatata, apoptosi, fagosoma e elaborazione delle proteine nel reticolo endoplasmatico (Figura 3B).
Costruzione di reti PPI e screening genico hub
Tra i 32 C-DEG, tre geni non hanno mostrato interazioni note con altri C-DEG nel database STRING e sono stati quindi esclusi dalla costruzione della rete PPI. La piattaforma STRING è stata utilizzata per l'analisi PPI dei restanti 29 C-DEG, rivelando una rete composta da 29 nodi e 39 archi, visualizzata tramite Cytoscape (Figura 4A). Successivamente, otto C-DEG dell'hub (HSP90B1, PLAU, SLPI, TOP3A, NCF4, PRF1, TUBA1C e CS) sono stati identificati in base al loro grado di connettività (Figura 4B).
Differenze nell'infiltrazione delle cellule immunitarie e la loro correlazione con i C-DEG hub
Questa analisi si è concentrata esclusivamente sul dataset T2DM (GSE78721) per indagare preliminarmente i cambiamenti nel microambiente immunitario nel contesto del T2DM e la loro relazione con i geni hub. Sono stati utilizzati grafici di violino per visualizzare la relazione tra il dataset T2DM e i livelli di infiltrazione di 22 tipi di cellule immunitarie (Figura 5A; dettagli nella Tabella Supplementare 3). L'analisi di infiltrazione delle cellule immunitarie ha indicato una proporzione più alta di cellule dendritiche attivate nel gruppo T2DM rispetto al gruppo di controllo (P = 0,049), mentre la proporzione di cellule dendritiche a riposo era significativamente più bassa nel gruppo T2DM (P = 0,023). Inoltre, la relazione tra i livelli di espressione degli otto C-DEG hub e i livelli di infiltrazione dei 22 diversi tipi di cellule immunitarie è stata visualizzata utilizzando una mappa termica (Figura 5B). In particolare, l'espressione di TOP3A ha mostrato una forte correlazione positiva con il livello di infiltrazione delle cellule T regolatorie (Treg) (P < 0,001), mentre l'espressione della SLPI era fortemente correlata positivamente all'infiltrazione delle cellule T gamma delta (P < 0,001). Al contrario, l'espressione di NCF4 ha mostrato una forte correlazione negativa con il livello di infiltrazione delle cellule NK attivate (P < 0,001), l'espressione di PRF1 è stata fortemente correlata negativamente con il livello di infiltrazione dei macrofagi M2 (P < 0,001), e l'espressione di PLAU è stata fortemente correlata negativamente al livello di infiltrazione delle cellule B naïve (P < 0,001).
Rete candidata di farmaco target per hub C-DEG
Tra gli otto C-DEG hub identificati, CS e TOP3A non erano elencati come target di farmaci noti nel database DGIdb e quindi esclusi dall'analisi successiva di rete farmaco-target candidato. Per identificare farmaci con potenziali ruoli terapeutici nel T2DM e nell'ASCVD, i restanti sei C-DEG hub sono stati inviati al database DGIdb (dettagli nella Tabella Supplementare 4). Successivamente, Cytoscape è stato utilizzato per costruire una rete candidato farmaco-target che illustra le associazioni tra i C-DEG hub e i farmaci target candidati (Figura 6). Questa rete comprendeva sei C-DEG hub, 136 farmaci molecolari e 137 bordi. In particolare, sono stati identificati rispettivamente 82, 5, 3, 4, 1 e 42 farmaci target candidati per TUBA1C, NCF4, PRF1, SLPI, HSP90B1 e PLAU, con la maggioranza classificata come agenti antineoplastici, antinfiammatori e agenti antitrombotici.
Validazione dell'espressione dei C-DEG hub
A causa dell'assenza di CS e TOP3A come target candidati di farmaci nel database DGIdb e delle loro funzioni biologiche meno direttamente correlate alle strategie terapeutiche nel contesto della ricerca su T2DM e ASCVD, lo studio si è concentrato sui restanti sei centri C-DEG hub (HSP90B1, PLAU, SLPI, NCF4, PRF1 e TUBA1C) per la validazione sperimentale. Per approfondire i ruoli di questi C-DEG hub, sono stati condotti saghi qRT-PCR e western blot per valutare i livelli di espressione dei geni core tra diversi gruppi (Figura 7 e Figura 8). I livelli di espressione di TUBA1C, NCF4, SLPI, HSP90B1 e PLAU sono stati significativamente aumentati nei pazienti con T2DM e ASCVD rispetto ai gruppi di controllo sani corrispondenti, mentre l'espressione di PRF1 era sottoregolata. Questi risultati suggeriscono che questi geni fondamentali possano svolgere un ruolo significativo nella progressione dei pazienti con T2DM e ASCVD.
DISPONIBILITÀ DEI DATI:I dataset di espressione genica GSE78721 e GSE12288 utilizzati per l'analisi bioinformatica in questo studio sono disponibili nel database pubblico NCBI GEO. I dati grezzi della parte di validazione sperimentale generati in questo studio sono forniti nel File Supplementare 1.

Figura 1: Identificazione dei DEG nei dataset GSE78721 e GSE12288. (A,C) I grafici vulcanici dei DEG in GSE78721 e GSE12288, rispettivamente. |logFC|≥ 0,5 e P < 0,05 sono stati stabiliti come soglia di screening; i geni con regolazione al rialzo erano contrassegnati come rossi, mentre quelli con regolazione al ribasso erano contrassegnati come verdi. (B,D) Heatmap dei primi 50 DEG di GSE78721 e GSE12288, rispettivamente. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 2: Diagramma di Venn dei C-DEG da due dataset. (A) Diagramma di Venn dei C-DEG alzati da T2DM e ASCVD in dataset GSE78721 e GSE12288. (B) Diagramma di Venn dei C-DEG ridimensionati da T2DM e ASCVD in dataset GSE78721 e GSE12288. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 3: Analisi funzionale e di arricchimento delle vie dei C-DEG. Grafici a barre delle analisi di arricchimento GO (A) e KEGG (B). Le coordinate verticali indicano la descrizione di diversi elementi GO o percorsi KEGG, mentre le coordinate orizzontali indicano il numero di C-DEG arricchiti. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 4: Screening genico della rete PPI e hub dei C-DEG. (A) La rete PPI dei C-DEG. Il rosso e il blu rappresentano rispettivamente i C-DEG al rialzo e al rialzo. (B) Otto geni hub sono stati identificati tra i C-DEG. I nodi rappresentano un C-DEG, e gli archi rappresentano l'interazione tra i C-DEG. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 5: Analisi dell'infiltrazione delle cellule immunitari. (A) Differenze nell'infiltrazione di 22 tipi di cellule immunitarie tra il gruppo T2DM e il gruppo di controllo. Asse orizzontale: sottotipi di cellule immunitarie; Asse verticale: livello relativo di infiltrazione delle cellule immunitari. (B) Mappa termica di correlazione tra l'espressione di 8 C-DEG hub e i livelli di infiltrazione di 22 tipi di cellule immunitarie. Quadrati rossi: correlazione positiva; Quadrati blu: correlazione negativa. I colori più scuri indicano correlazioni più forti. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 6: Rete farmaco-target candidata per sei C-DEG hub. I nodi rossi rappresentano i C-DEG dell'hub, i nodi verdi rappresentano farmaci candidati mirati, e i bordi tra di essi rappresentano coppie di interazione. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 7: Validazione dell'espressione genica principale nei pazienti con T2DM e ASCVD tramite qRT-PCR. I livelli di espressione dell'mRNA di HSP90B1 (A), PLAU (B), SLPI (C), NCF4 (D), PRF1 (E) e TUBA1C (F) sono stati rilevati tramite qRT-PCR. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 8: Validazione dell'espressione genica centrale in pazienti con T2DM e ASCVD tramite Western blot. I livelli di espressione proteica di HSP90B1 (A), PLAU (B), SLPI (C), NCF4 (D), PRF1 (E) e TUBA1C (F) sono stati rilevati tramite saghi Western blot. Clicca qui per visualizzare una versione più grande di questa figura.
Tabella supplementare 1. Matrice di espressione genica per il dataset GSE78721, utilizzata per l'analisi trascritomica T2DM in questo studio. Clicca qui per scaricare questo file.
Tabella supplementare 2. Matrice di espressione genica per il dataset GSE12288, utilizzata per l'analisi trascritomica ASCVD in questo studio. Clicca qui per scaricare questo file.
Tabella supplementare 3: Analisi dell'infiltrazione delle cellule immunitarie nel dataset GSE78721. Clicca qui per scaricare questo file.
Tabella supplementare 4: Interazioni candidato farmaco-bersaglio per geni comuni espressi differenzialmente a centro. Clicca qui per scaricare questo file.
File supplementare 1: I dati grezzi della parte di validazione sperimentale generati in questo studio. Clicca qui per scaricare questo file.