Identificatie van DEG's
Twee datasets, GSE78721 en GSE12288, werden in deze studie gebruikt voor bio-informatica-analyse. In GSE78721 werden in totaal 2.266 DEG's geïdentificeerd tussen T2DM-patiënten en normale monsters, bestaande uit 1.209 upreguliere en 1.057 downgereguleerde genen (Figuur 1A). In GSE12288 werden 924 DEG's geverifieerd tussen ASCVD-patiënten en controlegroepen, waaronder 403 genen die overexpressief waren en 521 genen die werden verlaagd (Figuur 1C). De top 50 DEG's die in beide datasets zijn geïdentificeerd, worden weergegeven in heatmaps in Figuur 1B,D. Door deze twee datasets te kruisen, werden 32 C-DEG's geïdentificeerd, bestaande uit 20 opgereguleerde en 12 neergereguleerde genen (Figuur 2A,B).
GO- en KEGG-verrijking van de C-DEG's
Om de biologische functies en belangrijke routes die geassocieerd zijn met de 32 C-DEGs te verduidelijken, werden GO- en KEGG-analyses uitgevoerd. De top tien termen gerelateerd aan biologisch proces (BP), cellulaire component (CC) en moleculaire functie (MF) onder de C-DEG's worden weergegeven in Figuur 3A. De BP-functies omvatten voornamelijk regulatie van plasminogeenactivatie, regulatie van fibrinolyse, positieve regulatie van cardiocytendifferentiatie, negatieve regulatie van eiwitbinding en cardioblastdifferentiatie. De functies van de CC hebben betrekking op het serine-type endopeptidasecomplex, eiwitcomplexen die betrokken zijn bij cel-matrixadhesie, gecondenseerd chromosoom, specifiek granulemembraan en PML-lichaam. De MF omvat structurele bestanddelen van het cytoskelet, superoxide-genererende NADPH oxidase-activatoractivatie, intramoleculaire transferaseactiviteit (fosfotransferases), serine-type exopeptidase-activiteit en IgG-binding (Figuur 3A). Wat betreft KEGG-analyse waren de C-DEG's voornamelijk betrokken bij kritieke routes, waaronder virale myocarditis, aritmogene rechterventrikel cardiomyopathie, complement- en stollingscascades, prostaatkanker, hypertrofische cardiomyopathie, gedilateerde cardiomyopathie, apoptose, fagosoom en eiwitverwerking in het endoplasmatisch reticulum (Figuur 3B).
Constructie van PPI-netwerken en hubgen-screening
Van de 32 C-DEG's vertoonden drie genen geen bekende interacties met andere C-DEG's in de STRING-database en werden daarom uitgesloten van de PPI-netwerkconstructie. Het STRING-platform werd gebruikt voor PPI-analyse van de resterende 29 C-DEG's, waarbij een netwerk van 29 knooppunten en 39 randen werd onthuld, dat werd gevisualiseerd met Cytoscape (Figuur 4A). Vervolgens werden acht hub C-DEG's (HSP90B1, PLAU, SLPI, TOP3A, NCF4, PRF1, TUBA1C en CS) geïdentificeerd op basis van hun connectiviteitsgraad (Figuur 4B).
Verschillen in immuuncelinfiltratie en hun correlatie met hub C-DEG's
Deze analyse richtte zich uitsluitend op de T2DM-dataset (GSE78721) om voorlopig veranderingen in de immuunmicro-omgeving te onderzoeken in de context van T2DM en hun relatie met hubgenen. Vioolgrafieken werden gebruikt om de relatie tussen de T2DM-dataset en de infiltratieniveaus van 22 immuunceltypen te visualiseren (Figuur 5A; details in Aanvullende Tabel 3). De infiltratieanalyse van immuuncellen toonde een hoger aandeel geactiveerde dendritische cellen in de T2DM-groep vergeleken met de controlegroep (P = 0,049), terwijl het aandeel rustende dendritische cellen significant lager was in de T2DM-groep (P = 0,023). Bovendien werd de relatie tussen de expressieniveaus van de acht hub C-DEG's en de infiltratieniveaus van de 22 verschillende immuunceltypen gevisualiseerd met behulp van een heatmap (Figuur 5B). Opmerkelijk is dat TOP3A-expressie een sterke positieve correlatie vertoonde met het infiltratieniveau van regulerende T-cellen (Tregs) (P < 0,001), terwijl SLPI-expressie sterk positief gecorreleerd was met gamma-delta T-celinfiltratie (P < 0,001). Omgekeerd vertoonde NCF4-expressie een sterke negatieve correlatie met het infiltratieniveau van geactiveerde NK-cellen (P < 0,001), was PRF1-expressie sterk negatief gecorreleerd met het infiltratieniveau van M2-macrofagen (P < 0,001), en was de expressie van PLAU sterk negatief gerelateerd aan het infiltratieniveau van naïeve B-cellen (P < 0,001).
Kandidaat geneesmiddelendoelnetwerk voor hub C-DEG's
Van de acht geïdentificeerde hub C-DEGs werden CS en TOP3A niet vermeld als bekende geneesmiddeldoelen in de DGIdb-database en werden ze daarom uitgesloten van de daaropvolgende kandidaat-geneesmiddel-doelgroep netwerkanalyse. Om geneesmiddelen met potentiële therapeutische rollen in T2DM en ASCVD te identificeren, werden de resterende zes hub C-DEG's ingediend in de DGIdb-database (details in aanvullende tabel 4). Vervolgens werd Cytoscape ingezet om een kandidaat-geneesmiddel-doelnetwerk op te bouwen dat de verbanden tussen hub C-DEG's en kandidaat-gerichte geneesmiddelen illustreert (Figuur 6). Dit netwerk bestond uit zes hub C-DEG's, 136 moleculaire geneesmiddelen en 137 randen. Specifiek werden 82, 5, 5, 4, 1 en 42 kandidaat-doelgeneesmiddelen geïdentificeerd voor respectievelijk TUBA1C, NCF4, PRF1, SLPI, HSP90B1 en PLAU, waarbij de meerderheid werd geclassificeerd als antineoplastische middelen, ontstekingsremmende middelen en antitrombotische middelen.
Validatie van hub C-DEGs-expressie
Vanwege het ontbreken van CS en TOP3A als kandidaat-geneesmiddeldoelwitten in de DGIdb-database en hun biologische functies minder direct gerelateerd aan therapeutische strategieën in de context van T2DM- en ASCVD-onderzoek, richtte de studie zich op de overige zes hub-C-DEG's (HSP90B1, PLAU, SLPI, NCF4, PRF1 en TUBA1C) voor experimentele validatie. Om de rollen van deze hub C-DEGs verder te onderzoeken, werden qRT-PCR en western blot-assays uitgevoerd om de expressieniveaus van kerngenen over verschillende groepen te beoordelen (Figuur 7 en Figuur 8). De expressieniveaus van TUBA1C, NCF4, SLPI, HSP90B1 en PLAU waren significant geupreguleerd bij patiënten met T2DM en ASCVD vergeleken met corresponderende gezonde controlegroepen, terwijl de PRF1-expressie werd gedownreguleerd. Deze bevindingen suggereren dat deze kerngenen een belangrijke rol kunnen spelen in de progressie van patiënten met T2DM en ASCVD.
BESCHIKBAARHEID GEGEVENS:De genexpressiedatasets GSE78721 en GSE12288 gebruikt voor bio-informatische analyse in deze studie zijn beschikbaar in de openbare NCBI GEO-database. De ruwe gegevens uit het experimentele validatiegedeelte dat in deze studie wordt gegenereerd, worden verstrekt in Aanvullend Bestand 1.

Figuur 1: Identificatie van DEG's in de datasets GSE78721 en GSE12288. (A,C) De vulkaangrafieken van DEG's in respectievelijk GSE78721 en GSE12288. |logFC|≥ 0,5 en P < 0,05 werden vastgesteld als screeningsdrempel; De upreguliere genen werden als rood gemarkeerd, en de downgereguleerde genen als groen. (B,D) Heatmaps van de top 50 DEG's in respectievelijk GSE78721 en GSE12288. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 2: Venn-diagram van C-DEG's uit twee datasets. (A) Venn-diagram van opgereguleerde C-DEG's van T2DM en ASCVD in GSE78721 en GSE12288 datasets. (B) Venn-diagram van gedownreguleerde C-DEG's van T2DM en ASCVD in GSE78721 en GSE12288 datasets. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 3: Analyse van functionele en padverrijking van C-DEG's. Staafdiagrammen van GO (A) en KEGG (B) verrijkingsanalyses. Verticale coördinaten geven de beschrijving aan van verschillende GO-items of KEGG-paden, en horizontale coördinaten geven het aantal verrijkte C-DEGs aan. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 4: PPI-netwerk- en hubgenscreening van C-DEG's. (A) Het PPI-netwerk van C-DEG's. Rood en blauw vertegenwoordigen respectievelijk de upreguleerde en downgereguleerde C-DEG's. (B) Acht hubgenen werden geïdentificeerd onder de C-DEG's. Knooppunten vertegenwoordigen een C-DEG, en randen vertegenwoordigen de interactie tussen C-DEG's. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 5: Analyse van immuuncelinfiltratie. (A) Verschillen in de infiltratie van 22 immuunceltypen tussen de T2DM-groep en de controlegroep. Horizontale as: subtypen immuuncellen; Verticale as: relatieve infiltratieniveau van immuuncellen. (B) Correlatie-warmtekaart van de expressie van 8 hub C-DEG's en de infiltratieniveaus van 22 immuunceltypen. Rode vierkanten: positieve correlatie; Blauwe vierkanten: negatieve correlatie. Donkerdere kleuren duiden op sterkere correlaties. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 6: Kandidaat geneesmiddel-doelnetwerk voor zes hub C-DEG's. Rode knooppunten vertegenwoordigen hub C-DEG's, groene knooppunten geven kandidaatgerichte geneesmiddelen voor en de randen ertussen vertegenwoordigen interactieparen. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 7: Validatie van kerngenexpressie bij patiënten met T2DM en ASCVD door qRT-PCR. De mRNA-expressieniveaus van HSP90B1 (A), PLAU (B), SLPI (C), NCF4 (D), PRF1 (E) en TUBA1C (F) werden gedetecteerd door qRT-PCR. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 8: Validatie van kerngenexpressie bij patiënten met T2DM en ASCVD door Western blot. De eiwitexpressieniveaus van HSP90B1 (A), PLAU (B), SLPI (C), NCF4 (D), PRF1 (E) en TUBA1C (F) werden gedetecteerd door Western blot-assays. Klik hier om een grotere versie van deze figuur te bekijken.
Aanvullende Tabel 1. Genexpressiematrix voor de GSE78721 dataset, gebruikt voor de T2DM-transcriptomische analyse in deze studie. Klik hier om dit bestand te downloaden.
Aanvullende Tabel 2. Genexpressiematrix voor de GSE12288 dataset, gebruikt voor de ASCVD-transcriptomische analyse in deze studie. Klik hier om dit bestand te downloaden.
Aanvullende tabel 3: Analyse van immuuncelinfiltratie in de GSE78721 dataset. Klik hier om dit bestand te downloaden.
Aanvullende tabel 4: Kandidaat-geneesmiddel-doelinteracties voor hub-gemeenschappelijke differentieel expressieve genen. Klik hier om dit bestand te downloaden.
Aanvullend Bestand 1: De ruwe gegevens van het experimentele validatiegedeelte dat in deze studie is gegenereerd. Klik hier om dit bestand te downloaden.