Identification of Feature Genes and Key Biological Pathways in Immune-Mediated Necrotizing Myopathy: High-Throughput Sequencing and Bioinformatics Analysis

Kai Chen,Feng Xiao,Song Tan,Qiao Zhou,Li Zeng,Chun-yan Zhu,Jia-ying Bai

DOI: https://doi.org/10.1016/j.csbj.2023.03.019

IF: 6.155

2023-03-17

Computational and Structural Biotechnology Journal

Abstract:Background Immune-mediated necrotizing myopathy (IMNM), a subgroup of idiopathic inflammatory myopathies (IIMs), is characterized by severe proximal muscle weakness and prominent necrotic fibers but no infiltration of inflammatory cells. IMNM pathogenesis is unclear. This study investigated key biomarkers and potential pathways for IMNM using high-throughput sequencing and bioinformatics technology. Methods RNA sequencing was conducted in 18 IMNM patients and 10 controls. A combination of weighted gene coexpression network analysis (WGCNA) and differentially expressed gene (DEG) analysis was conducted to identify IMNM-related DEGs. Feature genes were screened out by employing the proteinprotein interaction (PPI) network, support vector machine-recursive feature elimination (SVM-RFE), and least absolute shrinkage selection operator (LASSO). Quantitative real-time polymerase chain reaction (qRTPCR) was performed to verify their differential expression, and the receiver operating characteristic curve (ROC) was used to evaluate their diagnostic efficiency. Functional enrichment analysis was applied to reveal the hidden functions of feature genes. Furthermore, 28 immune cell abundance patterns in IMNM samples were measured. Results We identified 193 IMNM-related DEGs that were aberrantly upregulated in the IMNM population and were closely associated with immune-inflammatory responses, regulation of skeletal and cardiac muscle contraction, and lipoprotein metabolism. With the help of the PPI network and the LASSO and SVM-RFE algorithms, three feature genes, LTK , MYBPH , and MYL4, were identified and further confirmed by qRTPCR. ROC curves among IMNM, dermatomyositis (DM), inclusion body myositis (IBM), and polymyositis (PM) samples validated the LTK and MYL4 genes as IMNM-specific feature markers. In addition, all three genes had a notable association with the autophagylysosome pathway and immune-inflammatory responses. Ultimately, IMNM displayed a marked immune-cell infiltrative microenvironment. The most significant correlation was found between CD4 T cells, CD8 T cells, macrophages, natural killer (NK) cells, and dendritic cells (DCs). Conclusions LTK , MYBPH , and MYL4 were identified as potential key molecules for IMNM and are believed to play a role in the autophagylysosome pathway and muscle inflammation. Graphical abstract Download : Download high-res image (142KB) Download : Download full-size image

biochemistry & molecular biology