A subscription to JoVE is required to view this content. Sign in or start your free trial.
This study investigated the relationships between non-alcoholic fatty liver disease (NAFLD) and myocardial infarction (MI) through co-expressed genes, identifying Thrombospondin 1 (THBS1) as a biomarker. Immuno-infiltration analysis revealed CD8+ T cells and neutrophils as key factors, with THBS1 showing potential as a diagnostic tool for NAFLD and MI.
Non-alcoholic fatty liver disease (NAFLD) and myocardial infarction (MI) are two major health burdens with significant prevalence and mortality. This study aimed to explore the co-expressed genes to understand the relationship between NAFLD and MI and identify potential crucial biomarkers of NAFLD-related MI using bioinformatics and machine learning. Functional enrichment analysis was conducted, a co-protein-protein interaction (PPI) network diagram was constructed, and support vector machine-recursive feature elimination (SVM-RFE) and least absolute shrinkage and selection operator (LASSO) techniques were employed to identify one differentially expressed gene (DEG), Thrombospondin 1 (THBS1). THBS1 demonstrated strong performance in distinguishing NAFLD patients (AUC = 0.981) and MI patients (AUC = 0.900). Immuno-infiltration analysis revealed significantly lower CD8+ T cells and higher neutrophil levels in patients with NAFLD and MI. CD8+ T cells and neutrophils were effective in distinguishing NAFLD/MI from healthy controls. Correlation analysis showed that THBS1 was positively correlated with CCR (chemokine receptor), MHC class (major histocompatibility complex class), neutrophils, parainflammation, and Tfh (follicular helper T cells), and negatively correlated with CD8+ T cells, cytolytic activity, and TIL (tumor-infiltrating lymphocytes) in NAFLD and MI patients. THBS1 emerged as a novel biomarker for diagnosing NAFLD/MI in comparison to healthy controls. The results indicate that CD8+ T cells and neutrophils could serve as inflammatory immune features for differentiating patients with NAFLD/MI from healthy individuals.
Non-alcoholic fatty liver disease (NAFLD) is a major public health issue with a prevalence of 25%-30%1. It has been reported that the prevalence of NAFLD is high among patients with diabetes2. However, the significance of NAFLD in non-diabetic patients is not yet clear. Studies have suggested that NAFLD plays an independent role in the pathogenesis of atherosclerosis3,4. Additionally, a meta-analysis has shown that NAFLD is closely associated with coronary artery calcification, endothelial dysfunction, and atherosclerosis, and has emerged as an independent cardiovascular risk factor5. The connection between NAFLD and cardiovascular disease still requires further investigation.
Myocardial infarction (MI) is a catastrophic disease that threatens health and imposes a significant economic burden on patients and their families worldwide6. MI is also a major cause of death in patients with NAFLD. A clinical study published in the British Medical Journal has shown that the risk of myocardial infarction in NAFLD patients is 1.17 times greater than in non-NAFLD patients7,8. Some studies have identified molecular pathways, including inflammation, oxidative stress, and lipid metabolism, that contribute to NAFLD-related MI9,10,11. However, the underlying mechanisms linking NAFLD and MI remain unclear. It is crucial to identify novel biomarkers related to the prognosis of NAFLD and MI.
The increasing prevalence of NAFLD, which affects a broad segment of the population, underscores a significant public health issue, particularly given its association with diabetes. However, the impact of NAFLD on non-diabetic patients remains poorly understood. NAFLD has been implicated in the pathogenesis of atherosclerosis and is recognized as an independent cardiovascular risk factor, being closely linked to coronary artery calcification, endothelial dysfunction, and atherosclerosis. Despite these associations, the precise mechanisms bridging NAFLD and cardiovascular diseases such as myocardial infarction (MI) require further elucidation. MI is one of the leading causes of mortality worldwide and imposes a significant economic burden. The risk of MI in NAFLD patients is notably higher than in those without NAFLD, highlighting the need for a deeper understanding of the molecular pathways connecting these conditions. While inflammation, oxidative stress, and lipid metabolism have been suggested as contributing factors, the exact mechanisms remain unclear. There is an urgent need to identify novel biomarkers that can provide insights into the prognosis and management of NAFLD-related MI.
Consequently, in this study, RNA microarray datasets for NAFLD and MI were downloaded from the National Center for Biotechnology Information-Gene Expression Omnibus (NCBI-GEO, see Table of Materials) to identify and analyze the interaction of differentially expressed genes (DEGs) between NAFLD and MI. Enrichment analysis, protein-protein interaction (PPI) network diagram construction, support vector machine-recursive feature elimination (SVM-RFE), and least absolute shrinkage and selection operator (LASSO) algorithms were used to identify hub genes12,13,14,15,16,17,18,19. Immuno-infiltration analysis was performed to examine immune cells in NAFLD and MI patients. Ultimately, these methods were integrated to elucidate the relationship between NAFLD and MI. Figure 1 illustrates the design sequence followed in this study. By combining bioinformatics, machine learning, and immuno-infiltration analysis, this study aims to contribute to the development of a novel medical decision support platform.
The key contributions of this article are: (1) Identification of co-expressed genes: The study highlights the relationship between NAFLD and MI by identifying co-expressed genes, offering a deeper understanding of the molecular links between these two conditions. (2) Application of bioinformatics and machine learning: Using bioinformatics and machine learning techniques, including support vector machine-recursive feature elimination (SVM-RFE)17 and least absolute shrinkage and selection operator (LASSO)19, the study identifies THBS1 as a differentially expressed gene. THBS1 demonstrates high performance in distinguishing NAFLD and MI patients from healthy controls. (3) Immuno-infiltration analysis: The study conducts an immuno-infiltration analysis, revealing significantly lower levels of CD8+ T cells and higher levels of neutrophils in patients with NAFLD and MI. (4) Correlation analysis: The research demonstrates that THBS1 is positively correlated with several immune factors, including CCR (chemokine receptor), MHC class I (major histocompatibility complex class I), neutrophils, parainflammation, and Tfh (follicular helper T) cells. It is negatively correlated with CD8+ T cells, cytolytic activity, and tumor-infiltrating lymphocytes (TIL).
The details of the databases, weblinks, and software/packages used are listed in the Table of Materials. The simulation parameters used are provided in Table 1.
1. Obtaining RNA microarray datasets
2. Identification of DEGs
3. Enrichment analysis
4. PPI analysis by constructing PPI networks
5. Candidate hub DEGs screening by applying machine learning algorithms
6. ROC curve construction for assessing diagnostic performance
7. Immuno-infiltration analysis to explore immune infiltration
8. Statistical analysis
The key findings of the proposed study are presented here, encompassing various analyses conducted to elucidate the molecular mechanisms underlying NAFLD and MI.
Identification of DEGs
In the GSE89632 dataset, 76 up-regulated and 20 down-regulated genes were identified as NAFLD-DEGs (Figure 2B,D), while the GSE66360 dataset revealed 118 up-regulated and 8 down-regulated genes as MI-DEGs (Figure 2C
The method described in this study has significant implications for research into the molecular mechanisms underlying NAFLD and MI. By identifying key biomarkers such as THBS1, the proposed protocol offers potential targets for both diagnostic and therapeutic interventions. This approach can be extended to other complex diseases involving multiple pathways and immune responses, facilitating the discovery of novel biomarkers and therapeutic targets. Moreover, the integration of bioinformatics and machine learning techniqu...
None.
This study was supported by the National Natural Science Foundation of China (No. 62271511, U21A200949), Yucai Foundation of General Hospital of the Southern Theatre Command (2022NZC011), the Guangzhou Science and Technology Program Project (2023A03J0170), the National Clinical Research Center for Geriatrics (NCRCG-PLAGH-2023006) and Guangdong Basic and Applied Basic Research Foundation (No.2020A1515010288, No.2021A1515220101).
Name | Company | Catalog Number | Comments |
Cytoscape | Cytoscape Consortium | Version 3.6.1 | Used for visualizing protein-protein interaction (PPI) networks |
MI dataset GSE66360 (https://www.ncbi.nlm.nih.gov/geo/). | NCBI-GEO database | - | To collect RNA microarray datasets for analysis |
R package clusterProfiler | Bioconductor | - | Used for GO, KEGG, and DO enrichment analyses |
R package ggplot2 | CRAN | - | Used for creating Venn diagrams and other visualizations |
R package GSEABase | Bioconductor | - | Used in conjunction with GSVA for gene set enrichment analysis |
R package GSVA | Bioconductor | - | Used for single-sample gene set enrichment analysis (ssGSEA) |
R package limma | Bioconductor | - | Used for identifying differentially expressed genes (DEGs) |
R package pheatmap | CRAN | - | Used for generating heatmaps |
R package venn | CRAN | - | Used for creating Venn diagrams |
RNA microarray datasets (GSE66360, GSE89632) | NCBI-GEO | - | Publicly available RNA microarray datasets used for analysis |
RStudio | RStudio, PBC | Version 1.4.1717 | Integrated development environment for R |
String database | STRING (www.string-db.org/) | - | Online tool for constructing PPI networks |
Request permission to reuse the text or figures of this JoVE article
Request PermissionExplore More Articles
This article has been published
Video Coming Soon
Copyright © 2025 MyJoVE Corporation. All rights reserved