Transcriptome analysis of skeletal muscle in dermatomyositis, polymyositis, and dysferlinopathy, using a bioinformatics approach

Abstract

Background: Polymyositis (PM) and dermatomyositis (DM) are two distinct subgroups of idiopathic inflammatory myopathies. Dysferlinopathy, caused by a dysferlin gene mutation, usually presents in late adolescence with muscle weakness, degenerative muscle changes are often accompanied by inflammatory infiltrates, often resulting in a misdiagnosis as polymyositis.

Objective: To identify differential biological pathways and hub genes related to polymyositis, dermatomyositis and dysferlinopathy using bioinformatics analysis for understanding the pathomechanisms and providing guidance for therapy development.

Methods: We analyzed intramuscular ribonucleic acid (RNA) sequencing data from seven dermatomyositis, eight polymyositis, eight dysferlinopathy and five control subjects. Differentially expressed genes (DEGs) were identified by using DESeq2. Enrichment analyses were performed to understand the functions and enriched pathways of DEGs. A protein-protein interaction (PPI) network was constructed, and clarified the gene cluster using the molecular complex detection tool (MCODE) analysis to identify hub genes.

Results: A total of 1,048, 179 and 3,807 DEGs were detected in DM, PM and dysferlinopathy, respectively. Enrichment analyses revealed that upregulated DEGs were involved in type 1 interferon (IFN1) signaling pathway in DM, antigen processing and presentation of peptide antigen in PM, and cellular response to stimuli in dysferlinopathy. The PPI network and MCODE cluster identified 23 genes related to type 1 interferon signaling pathway in DM, 4 genes (PDIA3, HLA-C, B2M, and TAP1) related to MHC class 1 formation and quality control in PM, and 7 genes (HSPA9, RPTOR, MTOR, LAMTOR1, LAMTOR5, ATP6V0D1, and ATP6V0B) related to cellular response to stress in dysferliniopathy.

Conclusion: Overexpression of genes related to the IFN1 signaling pathway and major histocompatibility complex (MHC) class I formation was identified in DM and PM, respectively. In dysferlinopathy, overexpression of HSPA9 and the mTORC1 signaling pathway genes was detected.