- 1Department of Neurology, Yonsei University College of Medicine, Seoul, Republic of Korea
- 2Brain Korea 21 PLUS Project for Medical Science, Yonsei University, Seoul, Republic of Korea
- 3Research Institute of Women's Disease, Sookmyumg Women's University, Seoul, Republic of Korea
- 4Department of Neurology, CHA Bundang Medical Center, School of Medicine, CHA University, Seongnam-si, Republic of Korea
- 5Department of Rehabilitation Medicine, Gangnam Severance Hospital, Yonsei University College of Medicine, Seoul, Republic of Korea
- 6Rehabilitation Institute of Neuromuscular Disease, Yonsei University College of Medicine, Seoul, Republic of Korea
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.
1 Introduction
Idiopathic inflammatory myopathies (IIMs) comprise a heterogeneous group of autoimmune conditions characterized by skeletal muscle weakness and inflammation. The main clinical subtypes of IIM are classified as follows: dermatomyositis (DM); polymyositis (PM); immune-mediated necrotizing myopathy (IMNM); overlap syndrome with myositis (overlap myositis [OM]), including anti-synthetase syndrome (ASS); and inclusion body myositis (IBM) (1). The pathogenic mechanisms underlying the different myositis types and subtypes remain unclear. In recent years, different approaches, including genomic, epigenetic, transcriptomic, proteomic, and autoantibody studies, have been used to understand IIM pathogenesis (2).
Studies involving expression profiling of patients with DM or PM showed that the genes related to major histocompatibility complex (MHC) class I (MHC-I), MHC class II (MHC-II), cytokines, adhesion molecules, the actin cytoskeleton, and immunoglobulins were differentially expressed in muscle tissue (3, 4). In particular, the type 1 interferon (IFN1) pathway is relevant to DM pathogenesis (5). Overexpression of IFN1-inducible genes has been detected in the muscle, peripheral blood, and skin of patients with DM (5–8). PM is generally considered a prototypic T cell–mediated autoimmune myopathy (9).
Dysferlinopathy, caused by a dysferlin gene mutation, usually presents in late adolescence with muscle weakness; the degenerative muscle changes are often accompanied by inflammatory infiltrates, resulting in misdiagnosis as PM (10). Dysferlin protein appears to play a role in calcium-dependent vesicle-mediated muscle membrane repair (11). The precise pathogenic mechanism underlying dysferlinopathy remains unclear; it is hypothesized that defective membrane repair in dysferlin-deficient muscle leads to muscular dystrophy associated with remarkable muscle inflammation (12). Previous research has shown that the inflammasome pathway components are upregulated and activated in the muscles of patients and mice with dysferlinopathy (13). Dysferlin-deficient muscles also exhibit ubiquitin–proteasomal pathway activation, complement system factor upregulation, and macrophage infiltration (14–16).
Ribonucleic acid sequencing (RNA-Seq) helps detect and quantify complementary deoxyribonucleic acid (cDNA) via next-generation sequencing. In recent years, bioinformatics analysis of gene expression profiles has played a critical role in studying human disease pathogenesis. Differences in the gene expression levels of disease and healthy groups can help detect cellular changes related to disease status, thereby suggesting pathomechanisms (17).
In this study, we aimed to understand the transcriptomic signature and pathomechanism of DM, PM, and dysferlinopathy. We generated transcriptome profiles of muscles from patients with DM, PM, or dysferlinopathy and analyzed their gene expression profiles by using bioinformatics. Additionally, we compared the transcriptome profiles of dysferlinopathy and PM to help identify potential biomarker genes for differentiating between these two clinically similar myopathies because dysferlinopathy is often misdiagnosed as PM.
2 Materials and methods
2.1 Participant enrollment
We reviewed the medical records from January 2002 to October 2016 in the myopathy database of Gangnam Severance Hospital. The study included 7, 8, 8, and 5 vastus lateralis muscle biopsy samples from 7 patients with DM, 8 patients with PM, 8 patients with dysferlinopathy, and 5 control participants, respectively. The 15 patients with DM or PM had been diagnosed on the basis of the 2017 European League Against Rheumatism/American College of Rheumatology classification criteria for adult and juvenile IIMs (18), by combining the data on clinical characteristics and serological, electrophysiological, and histological evaluations of skeletal muscle biopsies. The diagnosis of dysferlinopathy in the 8 patients had been genetically confirmed. We selected 5 control participants who fulfilled the following criteria: (i) normal muscle pathological features, (ii) normal serum creatine kinase (CK) level, (iii) absence of definite muscle weakness, and (iv) availability of vastus lateralis muscle biopsy sample.
2.2 Clinical data collection
Clinical records of all the participants were reviewed to obtain the following information: age at diagnosis, age at onset of symptoms, pattern of weakness on clinical examination, serum CK levels, electrophysiological results, and findings from muscle pathology analysis.
2.3 Ribonucleic acid sequencing
cDNA libraries were constructed using the TruSeq RNA library kit, with 1 μg total RNA. The protocol consisted of the following steps: polyA-selected RNA extraction, RNA fragmentation, random hexamer–primed reverse transcription, and 100-nt paired-end sequencing performed using the Illumina HiSeq2500 system. The libraries were quantified using quantitative polymerase chain reaction (qPCR) performed according to the qPCR Quantification Protocol Guide and qualified using an Agilent Technologies 2,100 Bioanalyzer. The sequencer reads were processed and then aligned to the hg19 reference genome by using the STAR v2.7.10a algorithm to perform the alignment (19). STAR is an algorithm for aligning high-throughput short RNA-Seq data to a reference genome; it was developed to overcome the speed and accuracy issue. The reference genome and annotation data were downloaded from the UCSC Table Browser.1
For gene annotation, the gencode version 19 annotation file was downloaded from GENCODE.2 While running the STAR alignment, “--quantMode” was enabled to create a gene count file. Appropriate read counts were selected on the basis of the sequencing protocol. These raw reads were deposited in the NCBI Sequence Read Archive (SRA) database (accession number: SRP149027).
2.4 Bioinformatics analysis
Adaptors were removed from the raw reads by using trim-galore (0.6.6). After the adaptors were trimmed, alignment was performed using an inhouse STAR alignment pipeline by following the recommended procedures (19). The GENCODE version 19 annotation database was used to annotate genes. Read counts, which were obtained from the STAR alignment results, were normalized to compare expression levels between samples by using the DESeq2 package. In addition, DESeq2 identified differentially expressed genes (DEGs). The criteria for identifying DEGs were false discovery rate (FDR) < 0.05 and log2 fold change ≥2 or ≤ −2. Principal component analysis (PCA) was performed to confirm how gene expression in each sample was characterized by pathogenic features. The overall expression pattern in each comparison group was observed using a gene expression heat map.
Gene enrichment analysis of DEGs was performed using Metascape (20).3 For coexpression analysis, the protein–protein interaction (PPI) network was constructed using the Metascape website; the results were further visualized using the Cytoscape software (version 3.9.1) (21). Next, molecular complex detection tool (MCODE; version 1.5.1) (22) was utilized to identify the most important modules in the PPI network. The cutoff values were set as follows: node score cutoff = 0.2, maximal depth = 100, K-Core = 2, and degree cutoff = 2. Significant genes were detected using MCODE analysis.
2.5 Ethics statement
The institutional review board of Gangnam Severance Hospital, South Korea, approved the research protocol (IRB No. 3–2021-0450). All participants provided written informed consent for clinical data collection and genetic analysis. Access to information that could identify individual participants was prevented during or after data collection.
3 Results
3.1 Clinical data
Clinical data review for a total of 15 patients with PM or DM confirmed that the diagnosis was acceptable according to the 2017 European League Against Rheumatism/American College of Rheumatology classification criteria (18). The clinical features of the patients with PM or DM have been provided in Table 1.
The clinical features of the eight patients with dysferlinopathy were as follows: miyoshi distal myopathy (three patients), limb–girdle muscular dystrophy (three patients), proximodistal type (one patient), and hyperCKemia (one patient). Immunohistochemical analysis of dysferlin showed total loss of expression in seven patients and decreased expression in one patient. Two patients with dysferlinopathy had a history of initial misdiagnosis as PM. The clinical features of the patients with dysferlinopathy have been provided in Table 2.
The control group consisted of two men and three women, with ages ranging from 28 to 47 years. Two patients had been diagnosed with cramps, and three patients had been diagnosed with psychogenic weakness (one patient), psychogenic dystonia (one patient), and myalgia (one patient).
3.2 Identification of DEGs
Bioinformatics analysis showed the presence of DEGs in all patient groups. PCA showed differences between the patient and control samples (Figure 1). Patients with DM could be clearly distinguished from those in the control group. Patients with dysferlinopathy were also observed as a distinct group, except for 1 patient (DYS8). In contrast, patients with PM could not be obviously distinguished from controls. The number of DEGs obtained in the analysis was 1,048, 179, and 3,807 for DM, PM, and dysferlinopathy, respectively (Figure 2).
Figure 1. PCA plot between each group of patients and the control. PCA visualized DM and dysferlinopathy as groups distinct from the control. The distinction between PM and the control was ambiguous. var., variance.
Figure 2. Heat maps of the DEGs of each patient group. (A) Heat map of the genes differentially expressed between DM and the control. (B) Heat map of the genes differentially expressed between PM and the control. (C) Heat map of the genes differentially expressed between dysferlinopathy and the control.
3.3 Enrichment analysis of DEGs
The DEGs were categorized by gene enrichment analysis, and the top 10 biological pathways of upregulated and downregulated DEGs in patients with DM, PM, or dysferlinopathy have been provided in Figure 3. The results revealed that 576 DEGs were upregulated in DM and were most abundant in the interferon alpha/beta signaling pathway (R-HSA-909733). Twenty-six DEGs mapped to that reactome pathway: IFI6, IFI27, IFI35, IFIT1, IFIT2, IFIT3, IFIT5, IFITM1, IFITM3, IRF1, MX1, MX2, OAS1, OAS2, OAS3, HLA-A, HLA-B, HLA-C, STAT2, ISG15, PSMB8, SOCS3, USP18, SAMHD1, XAF1, and RSAD2. Additional information regarding these genes is provided in Table 3. The findings showed that 472 DEGs were downregulated in DM; these genes were most abundant in the eukaryotic translation elongation pathway (R-HSA-156842). In addition, 110 DEGs mapped to that reactome pathway.
Figure 3. (A) Gene enrichment analysis of DEGs upregulated in DM. (B) Gene enrichment analysis of DEGs downregulated in DM. (C) Gene enrichment analysis of DEGs upregulated in PM. (D) Gene enrichment analysis of DEGs downregulated in PM. (E) Gene enrichment analysis of DEGs upregulated in dysferlinopathy. (F) Gene enrichment analysis of DEGs downregulated in dysferlinopathy.
Table 3. Upregulated DEGs involved in the interferon alpha/beta signaling pathway (R-HSA-909733) in DM.
In the case of PM, 112 DEGs were upregulated and were most abundant in the antigen processing and presentation of peptide antigen pathway (GO:0048002); 24 DEGs mapped to that GO pathway: B2M, CD74, PDIA3, HLA-C, HLA-DRA, TAP1, TAPBPL, CDKN1A, GADD45A, ITGAL, PSMD1, ICAM2, KIF22, SPTBN2, TNFRSF14, SOCS3, DTX3L, RHOC, ZNF627, CREB3L2, IFIT3, PLAU, GSDMD, and TBC1D10C. Additional information regarding these genes has been provided in Table 4. Sixty-seven DEGs were downregulated in PM; these genes were most abundant in the striated muscle contraction pathway (R-HSA-390522). Six DEGs (MYBPC2, TNNC2, TNNI2, TNNT3, SCN9A, and MYH1) mapped to that reactome pathway.
Table 4. Upregulated DEGs involved in the antigen processing and presentation of peptide antigen pathway (GO:0048002) in PM.
In the case of dysferlinopathy, 1,650 DEGs were upregulated; these genes were most abundant in the cellular response to stimuli (R-HSA-8953897). Furthermore, 261 DEGs were mapped to that reactome pathway. The other top-ranked biological pathways were as follows: VEGFA-VEGFR2 signaling pathway (WP3888), actin filament–based process (GO:0030029), muscle structure development (GO:0061061), pathway of neurodegeneration – multiple diseases (hsa05022), and vesicle-mediated transport (R-HSA-5653656). Furthermore, 2,157 DEGs were downregulated in dysferlinopathy; these genes were most abundant in the DNA metabolic process pathway (GO:0006259). Sixty-two DEGs were mapped to that GO pathway. The other top-ranked biological pathways were as follows: protein localization to organelle (GO:0033365), regulation of response to DNA damage stimulus (GO:2001020), cell cycle (R-HSA-1640170), organelle assembly (GO:0070925), and cellular amide metabolic process (GO:0043603).
3.4 Identification of key modules and genes, using co-expression analysis
In the case of the DEGs upregulated in DM, the PPI network comprised 370 DEGs and 1,311 interaction pairs. The MCODE analysis revealed 6 modules, including 71 nodes and 311 pairs. The key module (MCODE score = 11) contained 23 genes related to the interferon alpha/beta signaling pathway (R-HSA-909733), with 253 pairs (Figure 4A). The 23 genes in this module were closely linked to each other, and the degree score of all these genes was 22. These genes corresponded to DEGs mapped to the IFN1 signaling pathway in the gene enrichment analysis.
Figure 4. Co-expression analysis of DEGs in DM. (A) MCODE analysis of upregulated DEGs reveals the key module containing 23 genes related to the IFN1 signaling pathway (R-HSA-909733). (B) MCODE analysis of downregulated DEGs reveals the key module containing 41 genes related to the eukaryotic translation elongation pathway (R-HSA-156842). The 24 genes with the highest degree score of 40 are marked in red.
In the case of the DEGs downregulated in DM, the PPI network comprised 317 DEGs and 2,123 interaction pairs. The MCODE analysis revealed seven modules, including 78 nodes and 874 pairs. The key module (MCODE score = 19) contained 41 genes related to the eukaryotic translation elongation pathway (R-HSA-156842), with 794 pairs (Figure 4B). Of the 41 genes, the following 24 genes had the highest degree score of 40: RPL10A, RPL4, RPL7, RPL13, RPL21, RPL26, RPL27, RPL30, RPLP2, RPS3, RPS5, RPS11, RPS13, RPS14, RPS15, RPS15A, RPS16, RPS20, RPS21, RPS23, RPS25, RPS28, RPL35, and RPL13A.
In the case of DEGs upregulated in PM, the PPI network consisted of 59 DEGs and 63 interaction pairs. The MCODE analysis revealed 2 modules, including 7 nodes and 9 pairs. The key module (MCODE score = 1.5) contained 4 genes (PDIA3, HLA-C, B2M, and TAP1) related to the antigen processing and presentation pathway (GO:0048002; Figure 5A).
Figure 5. Co-expression analysis of DEGs in PM. (A) MCODE analysis of upregulated DEGs reveals the key module containing 4 genes related to the antigen processing and presentation of peptide antigen pathway (GO:0048002). (B) MCODE analysis of downregulated DEGs reveals the key module containing 5 genes related to the striated muscle contraction pathway (R-HSA-390522).
In the case of DEGs downregulated in PM, the PPI network consisted of 27 DEGs and 32 interaction pairs. The MCODE analysis results revealed 2 modules, including 8 nodes and 12 pairs. The key module (MCODE score = 1.7) contained 5 genes (MYBPC2, TNNI2, TNNT3, TNNC2, and MYH1) related to the striated muscle contraction pathway (R-HSA-390522; Figure 5B).
In the case of dysferlinopathy, no network was identified in the PPI analysis of up- or downregulated DEGs.
3.5 Identification of hub genes within important biological pathways in dysferlinopathy
We performed PPI analysis with 261 upregulated DEGs related to the cellular response to stimuli (R-HSA-8953897) to identify hub genes in this pathway. This PPI network comprised 261 DEGs and 5,420 interaction pairs. The MCODE analysis revealed 8 modules, including 185 nodes and 2,891 pairs. The key module (MCODE score = 4.12) contained 32 genes related to the cellular response to stress (R-HSA-2262752), with 134 pairs. Of the 32 genes, 4 (HSPA9, RPTOR, ATP6V0D1, and LAMTOR5) had the highest degree score of 13 and 3 (MTOR, LAMTOR1, and ATP6V0B) had a score of 12 (Figure 6A).
Figure 6. Co-expression analysis of DEGs mapped to important pathways in dysferlinopathy. (A) PPI with 261 upregulated DEGs related to the cellular response to stimuli (R-HSA-8953897). MCODE analysis reveals the key module containing 32 genes related to the cellular response to stress (R-HSA-2262752). The 4 genes with the highest degree score of 13 are marked in red, and the 3 genes with a score of 12 are marked in brown. (B) PPI with 62 downregulated DEGs related to the DNA metabolic process (GO:0006259). MCODE analysis reveals the key module containing 11 genes related to the cellular response to DNA damage stimulus (GO:0006974). The gene with the highest degree score of 8 is marked in red, the gene with a score of 7 in brown, and the 4 genes with a score of 6 in pink.
We performed PPI analysis with 62 downregulated DEGs related to the DNA metabolic process (GO:0006259) to identify hub genes in this pathway. The PPI network comprised 43 DEGs and 117 interaction pairs. The MCODE analysis revealed 3 modules, including 20 nodes and 40 pairs. The key module (MCODE score = 2.63) contained 11 genes (ATAD5, RAD9A, RAD17, TP53BP1, RFC5, RFC4, MCM9, WDR48, FAN1, EME1, and EME2) related to the cellular response to DNA damage stimulus (GO:0006974), with 29 pairs. Of the 11 genes, RFC4 was found to have the highest degree score of 8 (Figure 6B).
3.6 Comparative analysis between PM and dysferlinopathy
The results revealed that 2,869 DEGs were upregulated in dysferlinopathy and 1,444 DEGs in PM (Figure 7A). PPI analysis of the upregulated DEGs did not identify any networks for dysferlinopathy or PM. Gene enrichment analysis results for the DEGs upregulated in dysferlinopathy and PM have been presented in Figures 7B,C, respectively. The DEGs upregulated in dysferlinopathy were most abundant in the translation pathway (R-HSA-72766), and 775 DEGs were mapped to this reactome pathway. PPI analysis of the 775 DEGs mapped to the translation pathway (R-HSA-72766) did not identify any network. The DEGs upregulated in PM were most abundant in the protein localization to organelle pathway (GO:0033365), and 58 DEGs were mapped to this GO pathway. PPI analysis of 58 DEGs mapped to the protein localization to organelle pathway (GO:0033365) showed that the network was composed of 35 DEGs and 42 pairs. The MCODE analysis results revealed 1 module, including 4 nodes and 6 pairs (Figure 7D). These 4 genes (KDELR3, C0PB2, TMED7, and RAB1A) were related to endoplasmic reticulum (ER) to Golgi vesicle–mediated transport (GO:0006888).
Figure 7. Comparative analysis between PM and dysferlinopathy. (A) Heat map of the genes differentially expressed between PM and dysferlinopathy. (B) Gene enrichment analysis of DEGs upregulated in dysferlinopathy. (C) Gene enrichment analysis of DEGs upregulated in PM. (D) Co-expression analysis of DEGs mapped to the protein localization to organelle pathway (GO:0033365) in PM. MCODE analysis shows 1 module, including 4 genes related to endoplasmic reticulum to Golgi vesicle–mediated transport (GO:0006888).
4 Discussion
In this study, we performed bioinformatics analysis of the gene expression profiles obtained for the muscles of patients with DM, PM, or dysferlinopathy. Gene enrichment analysis revealed that upregulated DEGs were most abundant in the IFN1 signaling pathway in DM and the antigen processing and presentation of peptide antigen pathway in PM. The DEGs upregulated in dysferlinopathy were most abundant in the cellular response to stimuli.
4.1 IFN1 signaling pathway in DM
Analysis of muscle transcripts from patients with DM helped identify 26 upregulated DEGs involved in the interferon alpha/beta signaling pathway. Of these, 18 were IFN1-stimulated genes (ISGs), which included 2 ISG transcription factor–encoding genes (STAT2 and IRF1), 2 IFN1 negative regulator–encoding genes (SOCS3 and USP18), and 3 genes involved in MHC-I formation (HLA-A, HLA-B, and HLA-C). These upregulated DEGs have been schematized in Figure 8.
Figure 8. Upregulated DEGs associated IFN1 signaling pathway in DM. IFN1 signal via the IFNAR (interferon-alpha/beta receptor)-1 and IFNAR-2, which are expressed in muscle cell activate JAK/STAT pathway. Subsequent JAK1 and TYK2-dependent phosphorylation of IFNAR creates docking sites for STAT1 and STAT2. Phosphorylation of STAT1 and STAT2, which leads to heterodimerization and together with IRF9 to the formation of ISGF3. This heterotrimeric complex translocates to the nucleus and binds the ISRE (interferon-sensitive response element) sequence present in ISGs (interferon stimulated genes). Upregulated DEGs in this study is marked in red.
IFN1-related proteins are prominent in DM muscles and are considered a potential DM diagnostic biomarker (23, 24). Additionally, the IFN1 pathway biomarker in blood is highly correlated with clinical severity in DM (25). Analysis of different subgroups for DM muscle biopsies has shown that the IFN1 signature is a useful unifying pattern for all DM subgroups and can be used to differentiate DM from other diseases (26).
The ISGs such as MX, ISG15, OAS, IFITM, IFIT, and IFI observed in our results have also been reported as hub genes in a bioinformatics study performed using muscle from a patient with DM (27–29). Myxovirus resistance protein 1 (MxA) is a key mediator of the interferon-induced antiviral response against a variety of viruses. Sarcoplasmic MxA expression has been shown to be a DM hallmark, and the European NeuroMuscular Center (ENMC) 2018 DM classification criteria include perifascicular MxA overexpression as definite DM muscle biopsy findings (26). ISG15 is one of the most strongly induced ISGs; A previous study proposed that the ISG15 expression level in muscle alone could be used to reliably quantify IFN1 pathway activation in adult inflammatory myositis (30). In patients with juvenile DM, muscle ISG15 expression was found to be inversely correlated with the clinical and histological severity of the disease; these results suggest that ISG15 negatively regulates the IFN1 signaling pathway (31). Thus, ISG15 may be used as a biomarker for DM diagnosis and monitoring; further research on the use of ISG15 for this purpose is expected. The expression of USP18 and SOCS3, which are negative regulators of IFN1 pathway signaling (31–33), was also found to increase in the current study. The possibility of their involvement in regulating sustained IFN1 signaling pathway activation in DM can be considered. Among the other genes analyzed in this study, STAT2 and IRF1 which are ISG transcription factors, have been reported to be overexpressed in DM (34).
The mechanisms underlying IFN1 induction in DM remain unclear. IFN1 is secreted from dendritic cells because of Toll-like receptor (TLR) induction and from muscle cells because of retinoic acid–inducible gene 1 (RIG-1) activation (35). Direct response by a viral pathogen or secondary response associated with muscle tissue remodeling may cause TLR induction (36). An in vitro study showed that hypoxia triggers IFN1 production in muscles (37).
The IFN1 pathway plays an important role with respect to muscle fibers and vascular injury in DM pathophysiology (5, 38, 39). Autocrine IFN1 signaling was recently found to prevent muscle stem cell proliferation in DM, leading to muscle repair deficit (40). These in vitro pathogenic effects of IFN1 in muscle and endothelial cells can be prevented by pharmacologic blockade of IFN1 signaling by using JAK inhibitors or anti-IFN receptor antibodies (39, 40). A systematic literature review suggested that JAK inhibitors represent a viable treatment option for DM (41).
4.2 Antigen processing and presentation of peptide antigen pathway in PM
Analysis of muscle transcripts from patients with PM helped identify 4 upregulated DEGs (PDIA3, TAP1, B2M, and HLA-C) involved in antigen processing and peptide antigen presentation. All these genes encode proteins involved in MHC-I formation and quality control. Normal muscle fibers do not express MHC-I, but MHC-I overexpression has been histopathologically confirmed in the muscles of patients with PM or IBM, and the inflammatory response by the MHC-I–CD8 complex is a major pathological mechanism in PM (42, 43). Thus, our findings are consistent with those of previous studies.
Protein disulfide-isomerase A3 (PDIA3), transporter associated with antigen processing 1 (TAP1), and beta2-microglobulin (B2M) are involved in the MHC-I peptide loading complex (PLC) (44, 45). The PLC, a transient multisubunit membrane complex in the ER, is essential for establishing stable peptide–MHC-I complexes by stringent peptide proofreading and quality control processes (45). TAP, a central part of the PLC, ensures high local concentrations of translocated peptides in the ER lumen proximal to the PLC (45, 46). PDIA3 forms a disulfide-linked complex with tapasin as a component of the MHC-I loading complex, which is thought to either stabilize the complex or facilitate correct class I molecule assembly (47). These upregulated DEGs have been schematized in Figure 9.
Figure 9. Upregulated DEGs associated MHC-I peptide loading complex (PLC) in PM. Peptide-receptive MHC-I heterodimers are recruited by calreticulin to the PLC. The peptide transporter TAP1/2 shuttles peptides from the cytosol into a molecular basket formed by the editing modules. After tapasin-facilitated peptide loading and proofreading, MHC-I are released from the PLC and traffic via the Golgi compartment to the cell surface. Upregulated DEGs in this study is marked in red.
4.3 Cellular response to stress in dysferlinopathy
Analysis of muscle transcripts from patients with dysferlinopathy helped identify the following 7 upregulated DEGs involved in the cellular response to stress: HSPA9, RPTOR, MTOR, LAMTOR1, LAMTOR5, ATP6V0D1, and ATP6V0B.
4.3.1 HSPA9
HSPA9 encodes mortalin, a 74 kDa mitochondrial-resident protein also known as p66mot-1, mitochondrial stress-70 protein (mtHsp70), peptide-binding protein 74 (PBP74), and glucose-regulated protein 75 (GRP75) (48, 49). Mortalin is located in the mitochondria, ER, plasma membrane, cytoplasmic vesicles, and cytosol (50). Its activity and function are determined by its cellular localization and binding partners (48).
GRP75 (mortalin) levels in human skeletal muscles increase with oxidative stress (51, 52). GRP75 expression has been found to increase in the muscles of patients with inflammatory myositis. GRP75 levels were found to increase in nonregenerating myofibers of patients who had myositis and were positive for sarcolemmal MHC-I immunoreactivity (53). GRP75 is also involved in mediating endo/sarcoplasmic reticulum (ER/SR)–mitochondria Ca2+ transport by forming the inositol 1,4,5-triphosphate receptor (IP3R)–GRP75–voltage-dependent anion channel 1 (VDAC1) complex in skeletal muscle cells (54). Changes in ER/SR-mitochondria Ca2+ transport may be related to impaired skeletal muscle function, such as in Duchenne muscular dystrophy (DMD) (55). Dysferlin is involved in T-tubule biogenesis and Ca2+ homeostasis maintenance (56–58). HSPA9 overexpression may be related to changes in intramuscular Ca2+ homeostasis. HSPA9 has varying functions depending on the cell type and interacting proteins; therefore, further research is required to clarify the significance of intramuscular HSPA9 expression in dysferlinopathy. Hypotheses regarding upregulated HSPA9 in dysferlinopathy has been schematized in Figure 10.
Figure 10. Hypotheses regarding upregulated HSPA9 in dysferlinopathy. (A) Membrane damage to muscle fibers with membrane repair dysfunction causes a prolonged releases of damage associated molecular patterns. These molecules stimulate leukocytes to produce pro-inflammatory cytokines and activate the complement system. The sustained inflammatory signaling may cause increased GRP75 in muscle cell. (B) In proposed model of dysferlin as a Ca2+-sensitive signaling scaffold localized to the t-tubule membrane, that is designed to respond to changes in intracellular Ca2+ caused by t-tubule membrane stress and damage. The SR(sarcoplasmic reticulum)-mitochondria Ca2+ transport is mediated by the macromolecular complex IP3R-GRP75-VDAC1. HSPA9 overexpression may be associated with changes in Ca2+ homeostasis in dysferlinopathy.
4.3.2 mTOR1 pathway
Regulatory-associated protein of MTOR complex 1 (RPTOR; RPTOR), MTOR, and late endosomal/lysosomal adaptor, MAPK, and MTOR activator (LAMTOR) are associated with the mTOR pathway (59). mTORC1 acts as an upstream regulator of autophagy induction in skeletal muscle (60). Chronic mTORC1 activation in old muscle leads to muscle atrophy mainly because of inability to induce autophagy (61). Autophagy dysregulation is related to various hereditary myopathies, including autophagic vacuolar myopathies (60). Autophagy was found to be impaired in patients with DMD and in mdx mouse models. In mdx mice, mTOR was constitutively activated, leading to autophagy-inducing gene downregulation at the molecular level (62). In dysferlinopathy, the autophagy program has been found to be activated for degradation of mutant dysferlin aggregates (15). Inhibition of mutant dysferlin aggregate formation in the ER by rapamycin, an mTOR inhibitor, has been reported in a cellular model (63). A recently published bioinformatics analysis involving dysferlinopathy identified the genes involved in the ubiquitin–proteasome pathway as key genes (64). Thus, our results suggest that the mTORC1 and autophagy pathways play a role in dysferlinopathy pathogenesis. The hypothesis of upregulated mTOR1 signaling in dysferlinopathy has been schematized in Figure 11.
Figure 11. Hypothesis regarding upregulated mTOR signaling pathway in dysferlinopathy. In dysferlinopathy, the autophagy program is activated for the degradation of mutant dysferlin aggregates. Chronic mTOR1 signaling inhibits autophagic activity. Upregulated mTOR pathway may be associated with dysregulated autophagy in dysferlinopathy.
4.4 DEGs downregulated in DM, PM, and dysferlinopathy
Twenty-four ribosome-associated genes were found to be downregulated on analyzing muscle transcripts from patients with DM. Decreased ribosomal-associated gene expression has recently been detected in the muscle and skin in juvenile DM (65). Analysis of muscle transcripts from patients with PM helped identify 5 downregulated DEGs (MYBPC2, TNNI2, TNNT3, TNNC2, and MYH1) related to striated muscle contraction. These findings may reflect muscle degeneration in DM and PM.
Analysis of muscle transcripts from patients with dysferlinopathy identified 6 downregulated DEGs (RFC4, RFC5, ATAD5, TP53BP1, WDR48, and RAD17) related to the cellular response to DNA damage stimulus. Replication factor C (RFC) plays an important role in DNA replication and repair, cell proliferation, cell cycle checkpoint regulation, and cell growth under stress (66). However, the implications of reduced RFC expression in dysferlinopathy require further study because association between myopathy and RFC expression has not been reported.
4.5 Comparative analysis between PM and dysferlinopathy
Comparative analysis between PM and dysferlinopathy revealed that the DEGs upregulated in PM were most abundant in the protein localization to organelle pathway and helped identify four genes (KDELR3, C0PB2, TMED7, and RAB1A) that were related to ER to Golgi vesicle–mediated transport. However, no promising pathway or gene that could be used as a biomarker for differentiating PM from dysferlinopathy was identified.
4.6 Limitations of this study
This study has several limitations. First, myositis-specific antibodies were not tested in patients with DM or PM. Additionally, patients with PM were not obviously distinguished from control, which may be related to the debatable issue about the real existence of PM as a separate pathobiological entity. PM is reported to be an overdiagnosed entity as subgroups of inflammatory myopathy are subdivided. Therefore, confirmation with a complete myositis autoantibody panel is recommended for the diagnosis of PM (67). However, we could not perform autoantibody test because we did not have serum samples from the patients. Second, 1 patient with dysferlinopathy (DYS8) showed a different gene expression pattern among the dysferlinopathy groups. This result was consistent with the pathological findings for dysferlin expression. Dysferlin expression was completely lost in the other patients; however, only this patient (DYS8) showed reduced expression. For the interpretation of this results, it would be helpful to collect additional dysferlinopathy cases in which the dysferlin was expressed and analyze the intramuscular transcriptome profiles of these patients. Third, the sample size of the dataset used in this study was small; therefore, the results obtained from the bioinformatics analysis should be further verified. Fourth, several of the hub genes identified in this study have not been previously reported to be associated with dysferlinopathy; therefore, further in vitro and functional studies are required to verify the molecular biological mechanisms of the genes identified in this study.
5 Conclusion
We identified differential biological pathways on the basis of bioinformatics analysis of gene expression profiles obtained from the muscles of patients with DM, PM, or dysferlinopathy. Overexpression of genes related to the IFN1 signaling pathway was identified in DM, and overexpression of genes related to MHC-I formation and quality control was identified in PM. In dysferlinopathy, overexpression of HSPA9 and the mTORC1 signaling pathway genes was detected. Further studies are required to verify the role of HSPA9 and the mTORC1 pathway in the pathomechanism of dysferlinopathy.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, SRP149027.
Ethics statement
The studies involving humans were approved by the institutional review board of Gangnam Severance Hospital, South Korea, approved the research protocol (IRB No. 3–2021-0450). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
H-NJ: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing. TL: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing. HP: Conceptualization, Data curation, Methodology, Supervision, Writing – review & editing. YY: Supervision, Writing – review & editing. S-HO: Supervision, Writing – review & editing. S-WK: Supervision, Writing – review & editing. Y-CC: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Supervision, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
1. Schmidt, J. Current classification and Management of Inflammatory Myopathies. J Neuromuscul Dis. (2018) 5:109–29. doi: 10.3233/JND-180308
2. Gao, S, Luo, H, Zhang, H, Zuo, X, Wang, L, and Zhu, H. Using multi-omics methods to understand dermatomyositis/polymyositis. Autoimmun Rev. (2017) 16:1044–8. doi: 10.1016/j.autrev.2017.07.021
3. Zhou, X, Dimachkie, MM, Xiong, M, Tan, FK, and Arnett, FC. Cdna microarrays reveal distinct gene expression clusters in idiopathic inflammatory myopathies. Med Sci Monit. (2004) 10:Br191–7.
4. Suárez-Calvet, X, Gallardo, E, Nogales-Gadea, G, Querol, L, Navas, M, Díaz-Manera, J, et al. Altered rig-I/Ddx58-mediated innate immunity in dermatomyositis. J Pathol. (2014) 233:258–68. doi: 10.1002/path.4346
5. Greenberg, SA, Pinkus, JL, Pinkus, GS, Burleson, T, Sanoudou, D, Tawil, R, et al. Interferon-alpha/Beta-mediated innate immune mechanisms in dermatomyositis. Ann Neurol. (2005) 57:664–78. doi: 10.1002/ana.20464
6. Baechler, EC, Bauer, JW, Slattery, CA, Ortmann, WA, Espe, KJ, Novitzke, J, et al. An interferon signature in the peripheral blood of dermatomyositis patients is associated with disease activity. Mol Med. (2007) 13:59–68. doi: 10.2119/2006-00085.Baechler
7. Walsh, RJ, Kong, SW, Yao, Y, Jallal, B, Kiener, PA, Pinkus, JL, et al. Type I interferon-inducible gene expression in blood is present and reflects disease activity in dermatomyositis and polymyositis. Arthritis Rheum. (2007) 56:3784–92. doi: 10.1002/art.22928
8. Wong, D, Kea, B, Pesich, R, Higgs, BW, Zhu, W, Brown, P, et al. Interferon and biologic signatures in dermatomyositis skin: specificity and heterogeneity across diseases. PloS One. (2012) 7:e29161. doi: 10.1371/journal.pone.0029161
9. Venalis, P, and Lundberg, IE. Immune mechanisms in polymyositis and dermatomyositis and potential targets for therapy. Rheumatology (Oxford). (2014) 53:397–405. doi: 10.1093/rheumatology/ket279
10. Amato, AA, and Brown, RH Jr. Dysferlinopathies. Handb Clin Neurol. (2011) 101:111–8. doi: 10.1016/b978-0-08-045031-5.00007-4
11. Han, R, and Campbell, KP. Dysferlin and muscle membrane repair. Curr Opin Cell Biol. (2007) 19:409–16. doi: 10.1016/j.ceb.2007.07.001
12. Han, R. Muscle membrane repair and inflammatory attack in Dysferlinopathy. Skelet Muscle. (2011) 1:10. doi: 10.1186/2044-5040-1-10
13. Mariano, A, Henning, A, and Han, R. Dysferlin-deficient muscular dystrophy and innate immune activation. FEBS J. (2013) 280:4165–76. doi: 10.1111/febs.12261
14. Han, R, Frett, EM, Levy, JR, Rader, EP, Lueck, JD, Bansal, D, et al. Genetic ablation of complement C3 attenuates muscle pathology in Dysferlin-deficient mice. J Clin Invest. (2010) 120:4366–74. doi: 10.1172/jci42390
15. Fanin, M, Nascimbeni, AC, and Angelini, C. Muscle atrophy, ubiquitin-proteasome, and Autophagic pathways in Dysferlinopathy. Muscle Nerve. (2014) 50:340–7. doi: 10.1002/mus.24167
16. Roche, JA, Tulapurkar, ME, Mueller, AL, van Rooijen, N, Hasday, JD, Lovering, RM, et al. Myofiber damage precedes macrophage infiltration after in vivo injury in Dysferlin-deficient a/J mouse skeletal muscle. Am J Pathol. (2015) 185:1686–98. doi: 10.1016/j.ajpath.2015.02.020
17. Ji, F, and Sadreyev, RI. Rna-Seq: basic bioinformatics analysis. Curr Protoc Mol Biol. (2018) 124:e68. doi: 10.1002/cpmb.68
18. Lundberg, IE, Tjärnlund, A, Bottai, M, Werth, VP, Pilkington, C, de Visser, M, et al. 2017 European league against rheumatism/American College of Rheumatology Classification Criteria for adult and juvenile idiopathic inflammatory myopathies and their major subgroups. Arthritis Rheumatol. (2017) 69:2271–82. doi: 10.1002/art.40320
19. Dobin, A, Davis, CA, Schlesinger, F, Drenkow, J, Zaleski, C, Jha, S, et al. Star: ultrafast universal Rna-Seq aligner. Bioinformatics. (2013) 29:15–21. doi: 10.1093/bioinformatics/bts635
20. Zhou, Y, Zhou, B, Pache, L, Chang, M, Khodabakhshi, AH, Tanaseichuk, O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. (2019) 10:1523. doi: 10.1038/s41467-019-09234-6
21. Shannon, P, Markiel, A, Ozier, O, Baliga, NS, Wang, JT, Ramage, D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303
22. Bader, GD, and Hogue, CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. (2003) 4:2. doi: 10.1186/1471-2105-4-2
23. Suárez-Calvet, X, Gallardo, E, Pinal-Fernandez, I, De Luna, N, Lleixà, C, Díaz-Manera, J, et al. Rig-I expression in Perifascicular Myofibers is a reliable biomarker of dermatomyositis. Arthritis Res Ther. (2017) 19:174. doi: 10.1186/s13075-017-1383-0
24. Uruha, A, Nishikawa, A, Tsuburaya, RS, Hamanaka, K, Kuwana, M, Watanabe, Y, et al. Sarcoplasmic Mxa expression: a valuable marker of dermatomyositis. Neurology. (2017) 88:493–500. doi: 10.1212/wnl.0000000000003568
25. Huard, C, Gullà, SV, Bennett, DV, Coyle, AJ, Vleugels, RA, and Greenberg, SA. Correlation of cutaneous disease activity with type 1 interferon gene signature and interferon Β in dermatomyositis. Br J Dermatol. (2017) 176:1224–30. doi: 10.1111/bjd.15006
26. Mammen, AL, Allenbach, Y, Stenzel, W, and Benveniste, O. 239th Enmc international workshop: classification of dermatomyositis, Amsterdam, the Netherlands, 14-16 December 2018. Neuromuscul Disord. (2020) 30:70–92. doi: 10.1016/j.nmd.2019.10.005
27. Xiao, L, Xiao, W, and Lin, S. Ten genes are considered as potential biomarkers for the diagnosis of dermatomyositis. PloS One. (2021) 16:e0260511. doi: 10.1371/journal.pone.0260511
28. Chen, S, Li, H, Zhan, H, Zeng, X, Yuan, H, and Li, Y. Identification of hub biomarkers and immune cell infiltration in polymyositis and dermatomyositis. Aging (Albany NY). (2022) 14:4530–55. doi: 10.18632/aging.204098
29. Ouyang, X, Zeng, Y, Jiang, X, Xu, H, and Ning, Y. Identification of vital hub genes and potential molecular pathways of dermatomyositis by bioinformatics analysis. Biomed Res Int. (2021) 2021:9991726–13. doi: 10.1155/2021/9991726
30. Pinal-Fernandez, I, Casal-Dominguez, M, Derfoul, A, Pak, K, Plotz, P, Miller, FW, et al. Identification of distinctive interferon gene signatures in different types of myositis. Neurology. (2019) 93:e1193–204. doi: 10.1212/WNL.0000000000008128
31. Hou, C, Durrleman, C, Periou, B, Barnerias, C, Bodemer, C, Desguerre, I, et al. From diagnosis to prognosis: revisiting the meaning of muscle Isg15 overexpression in juvenile inflammatory myopathies. Arthritis Rheumatol. (2021) 73:1044–52. doi: 10.1002/art.41625
32. Carow, B, and Rottenberg, ME. Socs3, a major regulator of infection and inflammation. Front Immunol. (2014) 5:58. doi: 10.3389/fimmu.2014.00058
33. Zhao, L, Wang, Q, Zhou, B, Zhang, L, and Zhu, H. The role of immune cells in the pathogenesis of idiopathic inflammatory myopathies. Aging Dis. (2021) 12:247–60. doi: 10.14336/ad.2020.0410
34. Kuriyama, Y, Shimizu, A, Kanai, S, Oikawa, D, Motegi, SI, Tokunaga, F, et al. Coordination of retrotransposons and type I interferon with distinct interferon pathways in dermatomyositis, systemic lupus erythematosus and autoimmune blistering disease. Sci Rep. (2021) 11:23146. doi: 10.1038/s41598-021-02522-6
35. Bolko, L, Jiang, W, Tawara, N, Landon-Cardinal, O, Anquetil, C, Benveniste, O, et al. The role of interferons type I, ii and iii in myositis: a review. Brain Pathol. (2021) 31:e12955. doi: 10.1111/bpa.12955
36. Cappelletti, C, Baggi, F, Zolezzi, F, Biancolini, D, Beretta, O, Severa, M, et al. Type I interferon and toll-like receptor expression characterizes inflammatory myopathies. Neurology. (2011) 76:2079–88. doi: 10.1212/WNL.0b013e31821f440a
37. De Luna, N, Suárez-Calvet, X, Lleixà, C, Diaz-Manera, J, Olivé, M, Illa, I, et al. Hypoxia triggers Ifn-I production in muscle: implications in dermatomyositis. Sci Rep. (2017) 7:8595. doi: 10.1038/s41598-017-09309-8
38. Franzi, S, Salajegheh, M, Nazareno, R, and Greenberg, SA. Type 1 interferons inhibit Myotube formation independently of upregulation of interferon-stimulated gene 15. PloS One. (2013) 8:e65362. doi: 10.1371/journal.pone.0065362
39. Ladislau, L, Suárez-Calvet, X, Toquet, S, Landon-Cardinal, O, Amelin, D, Depp, M, et al. Jak inhibitor improves type I interferon induced damage: proof of concept in dermatomyositis. Brain. (2018) 141:1609–21. doi: 10.1093/brain/awy105
40. Gallay, L, Fermon, C, Lessard, L, Weiss-Gayet, M, Viel, S, Streichenberger, N, et al. Involvement of type I interferon signaling in muscle stem cell proliferation during dermatomyositis. Neurology. (2022) 98:e2108–19. doi: 10.1212/wnl.0000000000200271
41. Paik, JJ, Lubin, G, Gromatzky, A, Mudd, PN Jr, Ponda, MP, and Christopher-Stine, L. Use of Janus kinase inhibitors in dermatomyositis: a systematic literature review. Clin Exp Rheumatol. (2022) 41:348–58. doi: 10.55563/clinexprheumatol/hxin6o
42. Ikenaga, C, Kubota, A, Kadoya, M, Taira, K, Uchio, N, Hida, A, et al. Clinicopathologic features of myositis patients with Cd8-Mhc-1 complex pathology. Neurology. (2017) 89:1060–8. doi: 10.1212/wnl.0000000000004333
43. Dalakas, MC. Inflammatory muscle diseases. N Engl J Med. (2015) 372:1734–47. doi: 10.1056/NEJMra1402225
44. Margulies, DH, Jiang, J, and Natarajan, K. Structural and dynamic studies of Tapbpr and Tapasin reveal the mechanism of peptide loading of Mhc-I molecules. Curr Opin Immunol. (2020) 64:71–9. doi: 10.1016/j.coi.2020.04.004
45. Blees, A, Januliene, D, Hofmann, T, Koller, N, Schmidt, C, Trowitzsch, S, et al. Structure of the human Mhc-I peptide-loading complex. Nature. (2017) 551:525–8. doi: 10.1038/nature24627
46. Neefjes, J, Jongsma, ML, Paul, P, and Bakke, O. Towards a systems understanding of Mhc class I and Mhc class ii antigen presentation. Nat Rev Immunol. (2011) 11:823–36. doi: 10.1038/nri3084
47. Mahmood, F, Xu, R, Awan, MUN, Song, Y, Han, Q, Xia, X, et al. Pdia3: structure, functions and its potential role in viral infections. Biomed Pharmacother. (2021) 143:112110. doi: 10.1016/j.biopha.2021.112110
48. Londono, C, Osorio, C, Gama, V, and Alzate, O. Mortalin, apoptosis, and neurodegeneration. Biomol Ther. (2012) 2:143–64. doi: 10.3390/biom2010143
49. Kaul, SC, Taira, K, Pereira-Smith, OM, and Wadhwa, R. Mortalin: present and prospective. Exp Gerontol. (2002) 37:1157–64. doi: 10.1016/s0531-5565(02)00135-3
50. Wadhwa, R, Taira, K, and Kaul, SC. An Hsp70 family chaperone, Mortalin/Mthsp70/Pbp74/Grp75: what, when, and where? Cell Stress Chaperones. (2002) 7:309–16. doi: 10.1379/1466-1268(2002)007<0309:ahfcmm>2.0.co;2
51. Gjøvaag, TF, and Dahl, HA. Effect of training and detraining on the expression of heat shock proteins in M. Triceps Brachii of untrained males and females. Eur J Appl Physiol. (2006) 98:310–22. doi: 10.1007/s00421-006-0281-y
52. Dalla Libera, L, Ravara, B, Gobbo, V, Tarricone, E, Vitadello, M, Biolo, G, et al. A transient antioxidant stress response accompanies the onset of disuse atrophy in human skeletal muscle. J Appl Physiol (1985). (2009) 107:549–57. doi: 10.1152/japplphysiol.00280.2009
53. Vitadello, M, Doria, A, Tarricone, E, Ghirardello, A, and Gorza, L. Myofiber stress-response in myositis: parallel investigations on patients and experimental animal models of muscle regeneration and systemic inflammation. Arthritis Res Ther. (2010) 12:R52. doi: 10.1186/ar2963
54. Zhang, SS, Zhou, S, Crowley-McHattan, ZJ, Wang, RY, and Li, JP. A review of the role of Endo/sarcoplasmic reticulum-mitochondria ca(2+) transport in diseases and skeletal muscle function. Int J Environ Res Public Health. (2021) 18:3874. doi: 10.3390/ijerph18083874
55. Pauly, M, Angebault-Prouteau, C, Dridi, H, Notarnicola, C, Scheuermann, V, Lacampagne, A, et al. Er stress disturbs Sr/Er-mitochondria ca(2+) transfer: implications in Duchenne muscular dystrophy. Biochim Biophys Acta Mol basis Dis. (2017) 1863:2229–39. doi: 10.1016/j.bbadis.2017.06.009
56. Kerr, JP, Ward, CW, and Bloch, RJ. Dysferlin at transverse tubules regulates ca(2+) homeostasis in skeletal muscle. Front Physiol. (2014) 5:89. doi: 10.3389/fphys.2014.00089
57. Hofhuis, J, Bersch, K, Büssenschütt, R, Drzymalski, M, Liebetanz, D, Nikolaev, VO, et al. Dysferlin mediates membrane Tubulation and links T-tubule biogenesis to muscular dystrophy. J Cell Sci. (2017) 130:841–52. doi: 10.1242/jcs.198861
58. Hofhuis, J, Bersch, K, Wagner, S, Molina, C, Fakuade, FE, Iyer, LM, et al. Dysferlin links excitation-contraction coupling to structure and maintenance of the cardiac transverse-axial tubule system. Europace. (2020) 22:1119–31. doi: 10.1093/europace/euaa093
59. Wullschleger, S, Loewith, R, and Hall, MN. Tor signaling in growth and metabolism. Cells. (2006) 124:471–84. doi: 10.1016/j.cell.2006.01.016
60. Xia, Q, Huang, X, Huang, J, Zheng, Y, March, ME, Li, J, et al. The role of autophagy in skeletal muscle diseases. Front Physiol. (2021) 12:638983. doi: 10.3389/fphys.2021.638983
61. Castets, P, Lin, S, Rion, N, Di Fulvio, S, Romanino, K, Guridi, M, et al. Sustained activation of Mtorc1 in skeletal muscle inhibits constitutive and starvation-induced autophagy and causes a severe. Late-Onset Myopathy Cell Metab. (2013) 17:731–44. doi: 10.1016/j.cmet.2013.03.015
62. De Palma, C, Morisi, F, Cheli, S, Pambianco, S, Cappello, V, Vezzoli, M, et al. Autophagy as a new therapeutic target in Duchenne muscular dystrophy. Cell Death Dis. (2012) 3:e418. doi: 10.1038/cddis.2012.159
63. Fujita, E, Kouroku, Y, Isoai, A, Kumagai, H, Misutani, A, Matsuda, C, et al. Two endoplasmic reticulum-associated degradation (Erad) Systems for the Novel Variant of the mutant Dysferlin: ubiquitin/proteasome Erad(I) and autophagy/lysosome Erad(ii). Hum Mol Genet. (2007) 16:618–29. doi: 10.1093/hmg/ddm002
64. Xie, Y, Li, YH, Chen, K, Zhu, CY, Bai, JY, Xiao, F, et al. Key biomarkers and latent pathways of Dysferlinopathy: bioinformatics analysis and in vivo validation. Front Neurol. (2022) 13:998251. doi: 10.3389/fneur.2022.998251
65. Roberson, EDO, Mesa, RA, Morgan, GA, Cao, L, Marin, W, and Pachman, LM. Transcriptomes of peripheral blood mononuclear cells from juvenile dermatomyositis patients show elevated inflammation even when clinically inactive. Sci Rep. (2022) 12:275. doi: 10.1038/s41598-021-04302-8
66. Li, Y, Gan, S, Ren, L, Yuan, L, Liu, J, Wang, W, et al. Multifaceted regulation and functions of replication factor C family in human cancers. Am J Cancer Res. (2018) 8:1343–55.
Keywords: dermatomyositis, polymyositis, dysferlinopathy, transcriptome, bioinformatics
Citation: Jeong H-N, Lee TG, Park HJ, Yang Y, Oh S-H, Kang S-W and Choi Y-C (2023) Transcriptome analysis of skeletal muscle in dermatomyositis, polymyositis, and dysferlinopathy, using a bioinformatics approach. Front. Neurol. 14:1328547. doi: 10.3389/fneur.2023.1328547
Edited by:
Massimiliano Filosto, University of Brescia, ItalyReviewed by:
Massimiliano Mirabella, Catholic University of the Sacred Heart, ItalyWladimir Bocca Vieira De Rezende Pinto, Federal University of São Paulo, Brazil
Copyright © 2023 Jeong, Lee, Park, Yang, Oh, Kang and Choi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Young-Chul Choi, eWNjaG9pQHl1aHMuYWM=
†These authors have contributed equally to this work and share first authorship