Identification of feature genes and key biological pathways in immune-mediated necrotizing myopathy: High-throughput sequencing and bioinformatics analysis

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 proteinprotein interaction (PPI) network, support vector machine-recursive feature elimination (SVM-RFE), and least absolute shrinkage selection operator (LASSO). Quantitative real-time polymerase chain reaction (qRTPCR) 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 qRTPCR. 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 autophagylysosome 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 autophagylysosome pathway and muscle inflammation.

Keywords

• IMNM
• RNA sequencing
• Biomarkers
• Immune cell infiltration
• Autophay-lysosome pathway