Bioinformatics analysis of vascular RNA-seq data revealed hub genes and pathways in a novel Tibetan minipig atherosclerosis model induced by a high fat/cholesterol diet

Background Atherosclerosis is a major contributor to cardiovascular events, however, its molecular mechanism remains poorly known. Animal models of atherosclerosis can be a valuable tool to provide insights into the etiology, pathophysiology, and complications of atherosclerosis. In particular, Tibetan minipigs are a feasible model for studying diet-related metabolic and atherosclerotic diseases. Methods We used vascular transcriptomics to identify differentially expressed genes (DEGs) in high fat/cholesterol (HFC) diet-fed Tibetan minipig atherosclerosis models, analyzed the DEGs gene ontology (GO) terms, pathways and protein-protein interactions (PPI) networks, and identified hub genes and key modules using molecular complex detection (MCODE), Centiscape and CytoHubba plugin. The identified genes were validated using the human carotid atherosclerosis database (GSEA 43292) and RT-PCR methods. Results Our results showed that minipigs displayed obvious dyslipidemia, oxidative stress, inflammatory response, atherosclerotic plaques, as well as increased low-density lipoprotein (LDL) and leukocyte recruitment after 24 weeks of HFC diet feeding compared to those under a regular diet. Our RNA-seq results revealed 1716 DEGs in the atherosclerotic/NC group, of which 1468 genes were up-regulated and 248 genes were down-regulated. Functional enrichment analysis of DEGs showed that the HFC diet-induced changes are related to vascular immune-inflammatory responses, lipid metabolism and muscle contraction, indicating that hypercholesterolemia caused by HFC diet can activate innate and adaptive immune responses to drive atherosclerosis development. Furthermore, we identified four modules from the major PPI network, which are implicated in cell chemotaxis, myeloid leukocyte activation, cytokine production, and lymphocyte activation. Fifteen hub genes were discovered, including TNF, PTPRC, ITGB2, ITGAM, VCAM1, CXCR4, TYROBP, TLR4, LCP2, C5AR1, CD86, MMP9, PTPN6, C3, and CXCL10, as well as two transcription factors (TF), i.e. NF-ĸB1 and SPI1. These results are consistent with the expression patterns in human carotid plaque and were validated by RT-PCR. Conclusions The identified DEGs and their enriched pathways provide references for the development and progression mechanism of Tibetan minipig atherosclerosis model induced by the HFC diet.