Identifikation von DEGs
Zwei Datensätze, GSE78721 und GSE12288, wurden für bioinformatische Analysen in dieser Studie verwendet. Im GSE78721 wurden insgesamt 2.266 DEGs zwischen T2DM-Patienten und normalen Proben identifiziert, darunter 1.209 hochregulierte und 1.057 herunterregulierte Gene (Abbildung 1A). Im Jahr GSE12288 wurden 924 DEGs zwischen ASCVD-Patienten und Kontrollgruppen verifiziert, darunter 403 überexprimierte Gene und 521 herunterregulierte Gene (Abbildung 1C). Die 50 wichtigsten DEGs, die in beiden Datensätzen identifiziert wurden, sind in Heatmaps dargestellt, die in Abbildung 1B,D dargestellt sind. Durch das Kreuzen dieser beiden Datensätze wurden 32 C-DEGs identifiziert, bestehend aus 20 hochregulierten und 12 herunterregulierten Genen (Abbildung 2A,B).
GO- und KEGG-Anreicherung der C-DEGs
Um die biologischen Funktionen und Schlüsselwege, die mit den 32 C-DEGs verbunden sind, zu erläutern, wurden GO- und KEGG-Analysen durchgeführt. Die zehn wichtigsten Begriffe im Zusammenhang mit biologischem Prozess (BP), zellulärer Komponente (CC) und molekularer Funktion (MF) unter den C-DEGs sind in Abbildung 3A dargestellt. Die BP-Funktionen umfassen überwiegend die Regulation der Plasminogenaktivierung, die Regulation der Fibrinolyse, die positive Regulation der Kardiozytendifferenzierung, die negative Regulation der Proteinbindung und die Kardioblastendifferenzierung. Die Funktionen des CC betreffen den Serin-Typ Endopeptidasekomplex, Proteinkomplexe, die an der Zell-Matrix-Adhäsion beteiligt sind, kondensiertes Chromosom, spezifische Granulatmembran und PML-Körper. Das MF umfasst strukturelle Bestandteile des Zytoskeletts, die Superoxid-erzeugende NADPH-Oxidase-Aktivatoraktivierung, die intramolekulare Transferase-Aktivität (Phosphotransferasen), die Serin-Typ-Exopeptidase-Aktivität und die IgG-Bindung (Abbildung 3A). Im Rahmen der KEGG-Analyse waren die C-DEGs hauptsächlich an kritischen Signalwegen beteiligt, darunter virale Myokarditis, arrhythmogene Rechtsventrikel-Kardiomyopathie, Komplement- und Koagulationskaskaden, Prostatakrebs, hypertrophe Kardiomyopathie, dilatierte Kardiomyopathie, Apoptose, Phagosom und Proteinverarbeitung im endoplasmatischen Retikulum (Abbildung 3B).
Aufbau von PPI-Netzwerken und Hub-Gen-Screening
Unter den 32 C-DEGs wiesen drei Gene keine bekannten Interaktionen mit anderen C-DEGs in der STRING-Datenbank auf und wurden daher vom PPI-Netzwerkaufbau ausgeschlossen. Die STRING-Plattform wurde für die PPI-Analyse der verbleibenden 29 C-DEGs eingesetzt und zeigte ein Netzwerk mit 29 Knoten und 39 Kanten, das mit Cytoscape visualisiert wurde (Abbildung 4A). Anschließend wurden acht Hub-C-DEGs (HSP90B1, PLAU, SLPI, TOP3A, NCF4, PRF1, TUBA1C und CS) basierend auf ihrem Konnektivitätsgrad identifiziert (Abbildung 4B).
Unterschiede in der Infiltration von Immunzellen und deren Korrelation zu Hub-C-DEGs
Diese Analyse konzentrierte sich ausschließlich auf den T2DM-Datensatz (GSE78721), um vorläufig Veränderungen im immunologischen Mikroumfeld im Kontext von T2DM und deren Beziehung zu Hub-Genen zu untersuchen. Violindiagramme wurden verwendet, um die Beziehung zwischen dem T2DM-Datensatz und den Infiltrationsniveaus von 22 Immunzelltypen zu visualisieren (Abbildung 5A; Details in ergänzender Tabelle 3). Die Analyse der Immunzellinfiltration zeigte einen höheren Anteil aktiverter dendritischer Zellen in der T2DM-Gruppe im Vergleich zur Kontrollgruppe (P = 0,049), während der Anteil der ruhenden dendritischen Zellen in der T2DM-Gruppe signifikant niedriger war (P = 0,023). Darüber hinaus wurde die Beziehung zwischen den Expressionsniveaus der acht Hub-C-DEGs und den Infiltrationsniveaus der 22 verschiedenen Immunzelltypen mithilfe einer Heatmap visualisiert (Abbildung 5B). Bemerkenswert ist, dass die TOP3A-Expression eine starke positive Korrelation mit dem Infiltrationsniveau regulatorischer T-Zellen (Tregs ) aufwies (P < 0,001), während die SLPI-Expression stark positiv mit der Gamma-Delta-T-Zellinfiltration korreliert war (P < 0,001). Umgekehrt zeigte die NCF4-Expression eine starke negative Korrelation mit dem Infiltrationsniveau aktiverter NK-Zellen (P < 0,001), die PRF1-Expression war stark negativ mit dem Infiltrationsniveau der M2-Makrophagen korreliert (P < 0,001) und die PLAU-Expression war stark negativ mit dem Infiltrationsniveau naiver B-Zellen verbunden (P < 0,001).
Kandidat-Arzneimittelziel-Netzwerk für Hub-C-DEGs
Unter den acht identifizierten Hub-C-DEGs wurden CS und TOP3A nicht als bekannte Arzneimittelziele in der DGIdb-Datenbank aufgeführt und wurden daher von der anschließenden Analyse des Kandidaten-Ziel-Netzwerks ausgeschlossen. Um Medikamente mit potenziellen therapeutischen Rollen bei T2DM und ASCVD zu identifizieren, wurden die verbleibenden sechs Hub-C-DEGs der DGIdb-Datenbank übermittelt (Details in ergänzender Tabelle 4). Anschließend wurde Cytoscape eingesetzt, um ein Kandidaten-Ziel-Netzwerk zu erstellen, das die Zusammenhänge zwischen Hub-C-DEGs und gezielten Kandidatenmedikamenten veranschaulicht (Abbildung 6). Dieses Netzwerk bestand aus sechs Hub-C-DEGs, 136 molekularen Medikamenten und 137 Kanten. Konkret wurden 82, 5, 5, 4, 1 und 42 Kandidaten-Zielmedikamente für TUBA1C, NCF4, PRF1, SLPI, HSP90B1 und PLAU identifiziert, wobei die Mehrheit als antineoplastische, entzündungshemmende und antithrombotische Medikamente klassifiziert wurde.
Validierung des Hub-C-DEG-Ausdrucks
Da CS und TOP3A als Kandidat-Wirkstoffziele in der DGIdb-Datenbank fehlen und ihre biologischen Funktionen weniger direkt mit therapeutischen Strategien im Kontext der T2DM- und ASCVD-Forschung zusammenhängen, konzentrierte sich die Studie auf die verbleibenden sechs Hub-C-DEGs (HSP90B1, PLAU, SLPI, NCF4, PRF1 und TUBA1C) für die experimentelle Validierung. Um die Rolle dieser Hub-C-DEGs weiter zu untersuchen, wurden qRT-PCR und Western-Blot-Assays durchgeführt, um die Expressionsniveaus von Kerngenen in verschiedenen Gruppen zu bewerten (Abbildung 7 und Abbildung 8). Die Expressionswerte von TUBA1C, NCF4, SLPI, HSP90B1 und PLAU wurden bei Patienten mit T2DM und ASCVD im Vergleich zu entsprechenden gesunden Kontrollgruppen signifikant hochreguliert, während die PRF1-Expression herunterreguliert war. Diese Ergebnisse deuten darauf hin, dass diese Kerngene eine bedeutende Rolle im Fortschreiten von Patienten mit T2DM und ASCVD spielen könnten.
DATENVERFÜGBARKEIT:Die Genexpressionsdatensätze, die GSE78721 und GSE12288 für die bioinformatische Analyse in dieser Studie verwendet werden, sind in der öffentlichen NCBI GEO-Datenbank verfügbar. Die Rohdaten aus dem experimentellen Validierungsteil, der in dieser Studie generiert wird, sind in Supplementary File 1 bereitgestellt.

Abbildung 1: Identifikation der DEGs in den Datensätzen GSE78721 und GSE12288. (A,C) Die Vulkandiagramme der DEGs in GSE78721 bzw. GSE12288. |logFC|≥ 0,5 und P < 0,05 wurden als Screening-Schwelle festgelegt; die hochregulierten Gene waren als rot markiert, die herunterregulierten Gene als grün. (B,D) Heatmaps der 50 besten DEGs in GSE78721 bzw. GSE12288. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 2: Venn-Diagramm der C-DEGs aus zwei Datensätzen. (A) Venn-Diagramm hochregulierter C-DEGs aus T2DM und ASCVD in GSE78721 und GSE12288 Datensätzen. (B) Venn-Diagramm der herunterregulierten C-DEGs aus T2DM und ASCVD in GSE78721 und GSE12288 Datensätzen. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 3: Funktionale und Weganreicherungsanalyse von C-DEGs. Balkendiagramme der GO (A) und KEGG (B) Anreicherungsanalysen. Vertikale Koordinaten geben die Beschreibung verschiedener GO-Elemente oder KEGG-Pfade an, und horizontale Koordinaten geben die Anzahl der angereicherten C-DEGs an. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 4: PPI-Netzwerk- und Hub-Gen-Screening von C-DEGs. (A) Das PPI-Netzwerk der C-DEGs. Rot und Blau stehen für die hochregulierten bzw. herunterregulierten C-DEGs. (B) Acht Hub-Gene wurden unter den C-DEGs identifiziert. Knoten stellen einen C-DEG dar, und Kanten stehen für die Interaktion zwischen C-DEGs. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 5: Analyse der Immunzellinfiltration. (A) Unterschiede in der Infiltration von 22 Immunzelltypen zwischen der T2DM-Gruppe und der Kontrollgruppe. Horizontale Achse: Subtypen von Immunzellen; Vertikale Achse: relative Infiltrationsstufe der Immunzellen. (B) Korrelations-Wärmekarte der Expression von 8 Hub-C-DEGs und der Infiltrationsniveaus von 22 Immunzelltypen. Rote Quadrate: positive Korrelation; Blaue Quadrate: negative Korrelation. Dunklere Farben deuten auf stärkere Korrelationen hin. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 6: Kandidat-Arzneimittelziel-Netzwerk für sechs Hub-C-DEGs. Rote Knoten stehen für Hub-C-DEGs, grüne Knoten für Kandidaten-Zielmedikamente und die Kanten zwischen ihnen für Interaktionspaare. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 7: Validierung der Kerngenexpression bei Patienten mit T2DM und ASCVD mittels qRT-PCR. Die mRNA-Expressionsniveaus von HSP90B1 (A), PLAU (B), SLPI (C), NCF4 (D), PRF1 (E) und TUBA1C (F) wurden durch qRT-PCR nachgewiesen. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 8: Validierung der Kerngenexpression bei Patienten mit T2DM und ASCVD durch Western Blot. Die Proteinexpressionswerte von HSP90B1 (A), PLAU (B), SLPI (C), NCF4 (D), PRF1 (E) und TUBA1C (F) wurden durch Western-Blot-Tests nachgewiesen. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Ergänzende Tabelle 1. Genexpressionsmatrix für den GSE78721-Datensatz, verwendet für die T2DM-Transkriptomanalyse in dieser Studie. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Tabelle 2. Genexpressionsmatrix für den GSE12288-Datensatz, verwendet für die ASCVD-Transkriptomanalyse in dieser Studie. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Tabelle 3: Analyse der Immunzellinfiltration im GSE78721-Datensatz. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Tabelle 4: Kandidaten-Arzneimittel-Ziel-Interaktionen für zentrale gemeinsame differenziell exprimierte Gene. Bitte klicken Sie hier, um diese Datei herunterzuladen.
Ergänzende Datei 1: Die Rohdaten aus dem experimentellen Validierungsteil, der in dieser Studie generiert wurde. Bitte klicken Sie hier, um diese Datei herunterzuladen.