Identificação de DEGs
Dois conjuntos de dados, GSE78721 e GSE12288, foram utilizados para análise de bioinformática neste estudo. Em GSE78721, um total de 2.266 DEGs foram identificados entre pacientes com DM T2 e amostras normais, compreendendo 1.209 genes regulados para cima e 1.057 genes para baixo (Figura 1A). Em GSE12288, 924 DEGs foram verificados entre pacientes com DSCV e controles, incluindo 403 genes superexpressos e 521 que foram regulados negativamente (Figura 1C). Os 50 principais DEGs identificados em ambos os conjuntos de dados estão representados em mapas de calor mostrados na Figura 1B,D. Ao intersectar esses dois conjuntos de dados, foram identificados 32 C-DEGs, consistindo em 20 genes regulados para cima e 12 genes para baixo (Figuras 2A,B).
Enriquecimento dos C-DEGs em GO e KEGG
Para elucidar as funções biológicas e as principais vias associadas aos 32 C-DEGs, foram realizadas análises GO e KEGG. Os dez principais termos relacionados a processos biológicos (BP), componente celular (CC) e função molecular (MF) entre os C-DEGs são apresentados na Figura 3A. As funções do BP incluem predominantemente a regulação da ativação do plasminôgeno, regulação da fibrinólise, regulação positiva da diferenciação de cardiocitos, regulação negativa da ligação de proteínas e diferenciação de cardioblastos. As funções CC referem-se ao complexo endopeptidase do tipo serina, complexos proteicos envolvidos na adesão à matriz celular, cromossomo condensado, membrana granular específica e corpo PML. A MF abrange constituintes estruturais do citoesqueleto, atividade ativadora da oxidase de NADPH geradora de superóxidos, atividade intramolecular de transferase (fosfotransferases), atividade de exopeptidase do tipo serina e ligação à IgG (Figura 3A). Em termos de análise KEGG, os C-DEGs estiveram principalmente envolvidos em vias críticas, incluindo miocardite viral, cardiomiopatia arritmogênica do ventrículo direito, cascatas de complemento e coagulação, câncer de próstata, cardiomiopatia hipertrófica, cardiomiopatia dilatada, apoptose, fagosomo e processamento de proteínas no retículo endoplasmático (Figura 3B).
Construção de redes de IBP e triagem gênica central
Entre os 32 C-DEGs, três genes não apresentaram interações conhecidas com outros C-DEGs no banco de dados STRING e, portanto, foram excluídos da construção da rede PPI. A plataforma STRING foi empregada para análise de PPI dos 29 C-DEGs restantes, revelando uma rede composta por 29 nós e 39 arestas, visualizada usando o Cytoscape (Figura 4A). Posteriormente, oito C-DEGs do hub (HSP90B1, PLAU, SLPI, TOP3A, NCF4, PRF1, TUBA1C e CS) foram identificados com base em seu grau de conectividade (Figura 4B).
Diferenças na infiltração de células imunes e sua correlação com os C-DEGs hub
Esta análise focou exclusivamente no conjunto de dados T2DM (GSE78721) para investigar preliminarmente as mudanças no microambiente imunológico no contexto da T2DM e sua relação com genes hub. Gráficos de violino foram utilizados para visualizar a relação entre o conjunto de dados T2DM e os níveis de infiltração de 22 tipos de células imunes (Figura 5A; detalhes na Tabela Suplementar 3). A análise de infiltração de células imunes indicou uma proporção maior de células dendríticas ativadas no grupo T2DM em comparação com o grupo controle (P = 0,049), enquanto a proporção de células dendríticas em repouso foi significativamente menor no grupo T2DM (P = 0,023). Além disso, a relação entre os níveis de expressão dos oito C-DEGs hub e os níveis de infiltração dos 22 diferentes tipos de células imunes foi visualizada usando um mapa de calor (Figura 5B). Notavelmente, a expressão de TOP3A exibiu forte correlação positiva com o nível de infiltração das células T reguladoras (Tregs) (P < 0,001), enquanto a expressão do SLPI foi fortemente correlacionada positivamente com a infiltração das células T gama delta (P < 0,001). Por outro lado, a expressão de NCF4 mostrou forte correlação negativa com o nível de infiltração das células NK ativadas (P < 0,001), a expressão de PRF1 foi fortemente correlacionada negativamente com o nível de infiltração dos macrófagos M2 (P < 0,001), e a expressão de PLAU esteve fortemente relacionada negativamente ao nível de infiltração das células B naïves (P < 0,001).
Rede candidata de fármacos-alvo para C-DEGs hub
Entre os oito C-DEGs centrais identificados, CS e TOP3A não foram listados como alvos de medicamentos conhecidos no banco de dados DGIdb e, portanto, foram excluídos da análise subsequente de rede candidata-alvo de fármacos. Para identificar medicamentos com potenciais papéis terapêuticos no DT2DM e na DCASC, os seis C-DEGs hubs restantes foram submetidos ao banco de dados DGIdb (detalhes na Tabela Suplementar 4). Posteriormente, o Cytoscape foi empregado para construir uma rede candidata-fármaco-alvo ilustrando as associações entre C-DEGs hub e fármacos candidatos-alvo (Figura 6). Essa rede era composta por seis C-DEGs hub, 136 medicamentos moleculares e 137 arestas. Especificamente, 82, 5, 3, 4, 1 e 42 fármacos-alvo candidatos foram identificados para TUBA1C, NCF4, PRF1, SLPI, HSP90B1 e PLAU, respectivamente, com a maioria classificada como agentes antineoplásicos, anti-inflamatórios e antitrombóticos.
Validação da expressão dos C-DEGs hub
Devido à ausência de CS e TOP3A como alvos candidatos a fármacos no banco de dados DGIdb e à falta de relação direta com suas funções biológicas com estratégias terapêuticas no contexto da pesquisa sobre T2DM e DCASC, o estudo focou nos seis C-DEGs centrais restantes (HSP90B1, PLAU, SLPI, NCF4, PRF1 e TUBA1C) para validação experimental. Para explorar mais profundamente os papéis desses C-DEGs hub, foram realizados ensaios qRT-PCR e western blot para avaliar os níveis de expressão dos genes centrais em diferentes grupos (Figura 7 e Figura 8). Os níveis de expressão de TUBA1C, NCF4, SLPI, HSP90B1 e PLAU foram significativamente aumentados em pacientes com DM 2 e DCS em comparação com grupos controle saudáveis correspondentes, enquanto a expressão de PRF1 foi reduzida. Esses achados sugerem que esses genes centrais podem desempenhar um papel significativo na progressão de pacientes com DMO2 e DCASC.
DISPONIBILIDADE DE DADOS:Os conjuntos de dados de expressão gênica GSE78721 e GSE12288 usados para análise de bioinformática neste estudo estão disponíveis no banco de dados público NCBI GEO. Os dados brutos da parte de validação experimental gerados neste estudo são fornecidos no Arquivo Suplementar 1.

Figura 1: Identificação de DEGs nos conjuntos de dados GSE78721 e GSE12288. (A,C) Os gráficos vulcânicos dos DEGs em GSE78721 e GSE12288, respectivamente. |logFC|≥ 0,5 e P < 0,05 foram estabelecidos como limiar de triagem; Os genes regulados para cima foram marcados como vermelhos, e os genes com efeito negativo foram marcados como verdes. (B,D) Mapas de calor dos 50 principais DEGs em GSE78721 e GSE12288, respectivamente. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 2: Diagrama de Venn dos C-DEGs de dois conjuntos de dados. (A) Diagrama de Venn dos C-DEGs aumentados de T2DM e ASCVD em conjuntos de dados GSE78721 e GSE12288. (B) Diagrama de Venn dos C-DEGs sob regulação reduzida de T2DM e ASCVD em conjuntos de dados GSE78721 e GSE12288. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 3: Análise de enriquecimento funcional e de vias de C-DEGs. Gráficos de barras das análises de enriquecimento GO (A) e KEGG (B). Coordenadas verticais indicam a descrição de diferentes itens GO ou caminhos KEGG, e coordenadas horizontais indicam o número de C-DEGs enriquecidos. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 4: Rede de PPI e triagem gênica central de C-DEGs. (A) A rede de PPI dos C-DEGs. Vermelho e azul representam, respectivamente, os C-DEGs regulados para cima e para baixo. (B) Oito genes hub foram identificados entre os C-DEGs. Nós representam um C-DEG, e as arestas representam a interação entre C-DEGs. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 5: Análise de infiltração de células imunes. (A) Diferenças na infiltração de 22 tipos de células imunes entre o grupo T2DM e o grupo controle. Eixo horizontal: subtipos de células imunes; Eixo vertical: nível relativo de infiltração das células imunes. (B) Mapa de calor de correlação da expressão de 8 C-DEGs hub e os níveis de infiltração de 22 tipos de células imunes. Quadrados vermelhos: correlação positiva; Quadrados azuis: correlação negativa. Cores mais escuras indicam correlações mais fortes. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 6: Rede de fármacos candidatos para seis C-DEGs hub. Nós vermelhos representam os C-DEGs do hub, nós verdes representam drogas candidatas-alvo, e as bordas entre eles representam pares de interação. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 7: Validação da expressão gênica central em pacientes com DMT 2 e DCASCV por qRT-PCR. Os níveis de expressão de mRNA HSP90B1 (A), PLAU (B), SLPI (C), NCF4 (D), PRF1 (E) e TUBA1C (F) foram detectados por qRT-PCR. Por favor, clique aqui para ver uma versão ampliada desta figura.

Figura 8: Validação da expressão gênica central em pacientes com DM2 e DCASCV por Western blot. Os níveis de expressão proteica HSP90B1 (A), PLAU (B), SLPI (C), NCF4 (D), PRF1 (E) e TUBA1C (F) foram detectados por ensaios Western blot. Por favor, clique aqui para ver uma versão ampliada desta figura.
Tabela Suplementar 1. Matriz de expressão gênica para o conjunto de dados GSE78721, usada para a análise transcriptômica T2DM neste estudo. Por favor, clique aqui para baixar este arquivo.
Tabela Suplementar 2. Matriz de expressão gênica para o conjunto de dados GSE12288, utilizada para a análise transcriptômica de DCASCV neste estudo. Por favor, clique aqui para baixar este arquivo.
Tabela suplementar 3: Análise de infiltração de células imunes no conjunto de dados GSE78721. Por favor, clique aqui para baixar este arquivo.
Tabela Suplementar 4: Interações candidatas fármaco-alvo para genes comuns diferencialmente expressos. Por favor, clique aqui para baixar este arquivo.
Arquivo Suplementar 1: Os dados brutos da parte de validação experimental gerados neste estudo. Por favor, clique aqui para baixar este arquivo.