Study on the differentially expressed genes and signaling pathways in dermatomyositis using integrated bioinformatics method.

Dermatomyositis is a common connective tissue disease. The occurrence and development of dermatomyositis is a result of multiple factors, but its exact pathogenesis has not been fully elucidated. Here, we used biological information method to explore and predict the major disease related genes of dermatomyositis and to find the underlying pathogenic molecular mechanism.The gene expression data of GDS1956, GDS2153, GDS2855, and GDS3417 including 94 specimens, 66 cases of dermatomyositis specimens and 28 cases of normal specimens, were obtained from the Gene Expression Omnibus database. The 4 microarray gene data groups were combined to get differentially expressed genes (DEGs). The gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichments of DEGs were operated by the database for annotation, visualization and integrated discovery and KEGG orthology based annotation system databases, separately. The protein-protein interaction networks of the DEGs were built from the STRING website. A total of 4097 DEGs were extracted from the 4 Gene Expression Omnibus datasets, of which 2213 genes were upregulated, and 1884 genes were downregulated. Gene ontology analysis indicated that the biological functions of DEGs focused primarily on response to virus, type I interferon signaling pathway and negative regulation of viral genome replication. The main cellular components include extracellular space, cytoplasm, and blood microparticle. The molecular functions include protein binding, double-stranded RNA binding and MHC class I protein binding. KEGG pathway analysis showed that these DEGs were mainly involved in the toll-like receptor signaling pathway, cytosolic DNA-sensing pathway, RIG-I-like receptor signaling pathway, complement and coagulation cascades, arginine and proline metabolism, phagosome signaling pathway. The following 13 closely related genes, XAF1, NT5E, UGCG, GBP2, TLR3, DDX58, STAT1, GBP1, PLSCR1, OAS3, SP100, IGK, and RSAD2, were key nodes from the protein-protein interaction network.This research suggests that exploring for DEGs and pathways in dermatomyositis using integrated bioinformatics methods could help us realize the molecular mechanism underlying the development of dermatomyositis, be of actual implication for the early detection and prophylaxis of dermatomyositis and afford reliable goals for the curing of dermatomyositis.

Introduction

Dermatomyositis is a non-suppurative inflammatory disease mainly involving rhabdomyositis and lymphocyte infiltration.[ It can be accompanied or not accompanied by multiple skin lesions.[ This disease is characterized by symmetrical weakness of limb band, neck and pharynx muscles, often involving multiple organs, and may be associated with tumors and other connective tissue diseases.[ The exact cause of the disease is not yet clear, and it is generally believed to be related to heredity and viral infection.[ The incidence of polymyositis and dermatomyositis has significant racial differences.[ African Americans had the highest incidence, with the proportion of blacks to whites ranging from 3 to 4:1.[ The incidence of dermatomyositis in children is higher in Asia and Africa than in Europe and America.[ At present, the diagnosis of dermatomyositis is mainly based on symmetrical proximal muscle weakness, pain and tenderness, accompanied by characteristic skin lesions such as periorbital purple-red edema, Gottron sign and nail fold rigid dilated capillary erythema.[ Increased muscle enzyme activity is one of the important serological indicators for the diagnosis of this disease.[ The serum myoglobin level in most patients with myositis increased and was parallel to the condition of the disease.[ Patients with dermatomyositis may develop autoantibodies such as anti-nuclear antibodies and anti-Jo-1 antibodies.[ However, the positive rate of diagnosis is low and early diagnosis is difficult. Glucocorticoid, methotrexate, azathioprine, and other immunosuppressive agents are commonly used medicines for this disease.[ For patients with dermatomyositis complicated by heart or lung disease, the condition is often serious, and the treatment effect is poor.[ In recent years, although there are some new immunotherapy programs for dermatomyositis, its efficacy still needs to be improved.[ Therefore, it is very important for the prevention and treatment of dermatomyositis to find new early diagnostic indicators and explore possible therapeutic targets.

The dynamic expression of tens of thousands of genes or gene transcripts can be detected on a large scale by DNA microarray technology.[ In recent years, DNA microarray technology has developed rapidly. By comparing gene expression differences between normal tissue cells and pathological tissue cells, a subset of characteristic genes closely related to disease types can be obtained.[ By analyzing the differential expression of genes in different pathological tissues and observing the correlation between genes, we can find the mechanism of pathological changes from the differential expression of genes.[ However, due to the high dimension, small sample size and non-linear distribution of gene expression profile data, direct data analysis often results in “dimension disaster” and “over-fitting” phenomena, which make the research progress slow, and the results are not ideal or even cannot be used.[ Therefore, it is necessary to develop an efficient feature gene selection method to process the data and select a gene subset that is of guiding significance for pathological changes, to facilitate assistant diagnosis.

Robust rank aggregation (RRA) method is an algorithm using probability model aggregation.[ It was proposed by KOLDE in 2012.[ For ordered gene list, this algorithm can make good use of the influence of parameters, outlier reduction, and noise of Geran model optimization algorithm on data results, leaving only the genes related to statistics in the final gene list.[ In research practice, many bioinformatics analysis tools can output gene lists sorted by some sort for us. When using different analysis methods or different data sets for analysis, many different gene lists can be output.[ RRA method can be used to correlate these different gene lists and obtain a new gene list that eliminates interference.[ In this study, RRA was used to analyze the differentially expressed genes (DEGs) of dermatomyositis, hoping to find out the DEGs of dermatomyositis, and provide new ideas for the early diagnosis and gene therapy of dermatomyositis.

First of all, we downloaded 4 original microarray datasets, GDS1956,[ GDS2153,[ GDS2855,[ and GDS3417,[ from the NCBI-Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/).[ The DEGs in Dermatomyositis and normal samples were screened by R software,[ and the gene ontology (GO) pathway enrichment analysis of DEGs was carried out by the database for annotation, visualization and integrated discovery (DAVID) (https://david.ncifcrf.gov/)[ and KEGG orthology based annotation system (KOBAS)-Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (http://kobas.cbi.pku.edu.cn/).[ Then, the STRING online database[ protein–protein interaction (PPI) network was used to analyze the correlation of DEGs, and the intermolecular interactions related to tumorigenesis were found. In summary, the gene fragments related to dermatomyositis occurrence and development were screened by dermatomyositis geographic database and analyzed comprehensively. The biological functions and key signaling pathways of these gene fragments were discussed, and the interaction network of coding proteins was analyzed. Our study produced results of causative factors for early prediction and judgment of dermatomyositis and provides effective drug targets for the remedy of dermatomyositis.

Methods

GEO database preserves many gene expression data generated by the scientific community. Researchers can observe a great deal of evidence to show the behavior of a gene or a group of genes.[ Users can mine different genes to obtain evidence to confirm laboratory findings and can also search for candidate genes in GEO database that deserve further research in the laboratory. Combining sequence information with expression information is helpful for functional annotation and characterization of unknown genes, or for discovering new roles of characteristic genes.[ These data are also valuable for genome-wide research, enabling biologists to examine global gene expression in different cell types and states, compare homologous genes with other species, and find repetitive patterns that help to form hypothetical co-regulatory transcripts in functional networks and pathways.[ Ethical approval was not necessary in our study because the expression profiles were downloaded from the public database and no new experiments in patients or animals were performed.

The gene expression profiles of GDS1956, GDS2153, GDS2855, and GDS3417 were downloaded using the keyword “dermatomyositis” searched in the GEO database. GDS1956 is a GPL96 platform, [HG-U133A] Affymetrix Human Genome U133A Array, including 10 samples of normal human skeletal muscle subsets and 21 samples of skeletal muscle subsets of dermatomyositis patients. The platform of GDS2153 is GPL96, [HG-U133A] Affymetrix Human Genome U133A Array, which contains 4 samples of normal human skeletal muscle subsets and 5 samples of skeletal muscle subsets of dermatomyositis patients. The platform of GDS2855 is GPL97, [HG-U133B] Affymetrix Human Genome U133B Array, which consists of 10 samples of normal human skeletal muscle subsets and 21 samples of skeletal muscle subsets of dermatomyositis patients. The platform of GDS3417 is GPL96, [HG-U133A] Affymetrix Human Genome U133A Array, which contains 4 samples of normal human skeletal muscle subsets and 19 samples of skeletal muscle subsets of dermatomyositis patients. The platform and series matrix files are downloaded as TXT files. Inter chip standardization can determine a reference chip, assume some invariants, and calibrate the data of other chips as a whole. The microarray data were processed by the expresso function in affy software package in 3 steps: background correction, standardization and summary, and the standardized gene expression matrix was obtained.

Screening for DEGs

The purpose of screening DEGs is to find genes with significant changes in different experimental samples. It is of great significance to find molecular markers, therapeutic targets, and molecular biological mechanisms that are helpful for disease diagnosis. The obtained DEG platform and matrix file series are converted by using R language software and annotation package. The ID corresponding to probe name is transformed to the international canonical name of gene (gene symbol) and stored in TXT file. Gene diverse expression analysis was made using Limma software package in Bioconductor Package (available online: http://www.bioconductor.org/).[

The correlative manipulating instruction codes were entered into R, and the DEGs in dermatomyositis and normal specimens of the 4 microarray datasets were calculated by the limma software package. Specimens with an adjusted P-value of <.05 and log fold change (FC) >2 were treated as DEGs. The TXT outcomes were saved for later calculation.

Integration of microarray data

The list of DEGs from the 4 microarray datasets calculated by limma packet analysis was stored as a TXT file. The RRA software package was downloaded, and R was applied to run the command code. A portion of genes that were up- or downregulated in the 4 chips were saved for later research. The RRA approach is openly workable in the Comprehensive R Network (http://cran.r-project.org/).[

GO and KEGG pathway enrichment analyses of DEGs

GO ontology analysis is to annotate and classify genes by biological process (BP), molecular function (MF), and cellular component. By GO analysis, not only the functions of differential genes can be identified, but also the functions of DEGs can be studied. Through the statistical analysis of the function enrichment degree of DEGs, we can calculate the p-value of the corresponding function of DEGs and the FDR value of P-value, that is, Q-value. Finally, combined with the research purpose, we can locate the most likely function related to DEGs. The analysis of KEGG pathway enrichment is based on KEGG pathway. Using fish statistical test, we can find out the signal pathway of significant enrichment of DEGs, which helps us to further clarify the biological function of genes, understand which metabolic pathway genes are enriched in, and screen out the key signal pathway.

The DAVID database is a bioinformatics database, which integrates biological data and analysis tools, provides systematic and comprehensive bioinformatics annotations for large-scale gene or protein lists, and identifies the most significant bioinformatics annotations by statistical methods. GO enrichment analysis were carried out by a DAVID online tool for tracing DEGs. KEGG pathway annotation of DEGs was carried out by the KOBAS online computing database. We logged in the DAVID and KOBAS online computing database, inputted the gene name of the DEGs, screened out the path with high enrichment degree according to the threshold value of “P < .05”, and then used R software for visualization.

PPI network integration

The STRING database (http://string-db.org/) is a software system for searching for protein-protein interactions, which helps to mine the core regulatory genes. It includes not only the direct physical interaction between proteins, but also the indirect functional correlation between proteins. Each node in the network graph represents a protein, and the lines between nodes represent the interaction between 2 proteins. Different colors correspond to different interaction types.[ All protein-protein interaction data are weighted, integrated, and there will be a calculated reliable value, and all calculated predictions can be viewed through a specific view. The corresponding proteins in the central node may be core proteins or key candidate genes with important physiological regulatory functions.[ The DEGs were input into the string database to construct the interaction network of differential gene coding proteins. By analyzing the network structure, the local subnetworks formed in the whole network and the key node proteins were obtained.

Results

Microarray data information and identification of DEGs in dermatomyositis

The dermatomyositis expression microarray datasets GDS1956, GDS2153, GDS2855, and GDS3417 were standardized, and the results are shown in Figure 1. When the GDS1956 dataset was screened by the limma package (corrected P-value < .05, logFC > 2), 681 DEGs were obtained. Among them, 468 upregulated genes and 195 downregulated gene were identified. 1362 DEGs were screened from the GDS2153 dataset, including 733 upregulated genes and 629 downregulated genes. 900 DEGs were screened from the GDS2855 dataset, including 600 upregulated genes and 300 downregulated genes. Additionally, 2675 DEGs were screened from the GDS3417 dataset, including 1376 upregulated genes and 1299 downregulated genes. The differential expression of multiple genes from the 2 sets of sample data included in each of the 4 microarrays is shown in Figure 2. The cluster heatmaps of the top 20 DEGs are shown in Figure 3.

Figure 1

Standardization of gene expression. (A) The standardization of GDS1956 data, (B) the standardization of GDS2153 data, (C) the standardization of GDS2855 data, (D) the standardization of GDS3417 data. The green bar represents the data before normalization, and the red bar represents the normalized data.

Figure 2

Differential expression of data between 2 sets of samples. (A) GDS1956 data, (B) GDS2153 data, (C) GDS2855 data, (D) GDS3417 data. The red points represent upregulated genes screened based on |fold change| >2.0 and a corrected P-value of <.05. The green points represent downregulation of the expression of genes screened based on |fold change| >2.0 and a corrected P-value of <.05. The black points represent genes with no significant difference. FC is the fold change.

Figure 3

Hierarchical clustering heatmap of DEGs screened on the basis of |fold change| >2.0 and a corrected P-value < .05. (A) GDS1956 data, (B) GDS2153 data, (C) GDS2855 data, (D) GDS3417 data. Red indicates that the expression of genes is relatively upregulated, green indicates that the expression of genes is relatively downregulated, and black indicates no significant changes in gene expression; gray indicates that the signal strength of genes was not high enough to be detected. DEGs = differentially expressed genes.

DEGS in dermatomyositis was identified by comprehensive bioinformatics. The microarray data of 4 kinds of dermatomyositis gene expression were analyzed by limma software package. The data were classified according to logarithmic variation, and then by RRA analysis was performed (corrected P-value < .05). The RRA method assumes that each gene in each experiment is randomly arranged. If a gene ranks higher in all experiments, the smaller its P-value, the greater the possibility of differential gene expression. Using R-heatmap software, the Heatmap of Top 20 up- and downregulated genes of each dataset was drawn, as shown in Figure 3.

GO term enrichment analysis of DEGs

Biological annotation of the DEGs in dermatomyositis identified from an integrated analysis of microarray data was performed using the DAVID online analysis tool and GO functional enrichments of up- and downregulated genes with a P-value of <.05 were obtained. GO analysis of DEGs was divided into 3 functional groups, including MF, biological processes, and cell composition.[ The results are shown in Figures 4 and 5. Significant results of the GO enrichment analysis of DEGs in dermatomyositis are shown in Table 1. The biological functions of DEGs focused primarily on response to virus, type I interferon signaling pathway and negative regulation of viral genome replication. The main cellular components include extracellular space, cytoplasm and blood microparticle. The MFs include protein binding, double-stranded RNA binding and MHC class I protein binding. These results indicate that most DEGs are significantly enriched in cell proliferation, binding, cell cycle regulation, and transcriptional activity.

Figure 4

GO enrichment analysis of DEGs in dermatomyositis. (A) GO analysis divided DEGs into 3 functional groups: molecular function, biological processes, and cell composition. (B) GO enrichment significance items of DEGs in different functional groups. DEGs = differentially expressed genes, GO = gene ontology.

Figure 5

Distribution of DEGs in dermatomyositis for different GO-enriched functions. The color of the stripe that each gene points to indicates the enriched gene pathway. The color depth of each gene block represents the amount of logFC. DEGs = differentially expressed genes, GO = gene ontology.

Table 1

GO analysis of DEGs associated with dermatomyositis.

KEGG pathway analysis of DEGs

Using the KOBAS online analysis database (http://kobas.cbi.pku.edu.cn/) to analyze the DEGs identified from dermatomyositis integrated gene microarrays, the most significantly enriched pathways of the DEGs were submitted to KEGG analysis. The results are shown in Table 2. The signaling pathways of DEGs were mainly enriched in the toll-like receptor signaling pathway, cytosolic DNA-sensing pathway, RIG-I-like receptor signaling pathway, complement and coagulation cascades, arginine and proline metabolism, phagosome signaling pathway in dermatomyositis. The data were imported into Cytoscape to calculate the topological characteristics of the network and determine each node.

Table 2

KEGG pathway analysis of DEGs associated with dermatomyositis.

Analysis of common genes in 4 microarray datasets

The number of unique and common genes in 4 kinds of microarray datasets was analyzed by Venn diagram,[ and the similarity and overlap of gene composition among different microarray samples were further compared. The results showed that the number of identical genes in the 4 microarray datasets was 13, as shown in Figure 6. The common genes are XAF1, NT5E, UGCG, GBP2, TLR3, DDX58, STAT1, GBP1, PLSCR1, OAS3, SP100, IGK, and RSAD2.

Figure 6

Venn diagram of 4 microarray datasets. Each color represents 1 expression microarray dataset. There are 13 common genes in the intersection of the 4 expression microarray datasets.

Analyzing DEGs in dermatomyositis using a PPI network

13 genes are common to 4 microarray datasets. The common expression products in dermatomyositis were constructed using the STRING database to construct PPI networks. After removing the isolated and partially connected nodes, a complex network of DEGs was constructed, as shown in Figure 7.

Figure 7

PPI network. Circles represent genes, lines represent the interaction of proteins between genes, and the results within the circle represent the structure of proteins. Line color represents evidence of the interaction between the proteins. PPI = protein–protein interaction.

Discussion

We built a PPI network of protein encoded by DEGs and screened out the following 13 genes that have distinct relation with dermatomyositis: XAF1, NT5E, UGCG, GBP2, TLR3, DDX58, STAT1, GBP1, PLSCR1, OAS3, SP100, IGK, and RSAD2. The proteins expressed by these genes are critical nodes in the PPI network. Pathway enrichment analysis pointed out that the genes mostly participated in the toll-like receptor signaling pathway, cytosolic DNA-sensing pathway, RIG-I-like receptor signaling pathway, complement and coagulation cascades, arginine and proline metabolism, phagosome signaling pathway.

X chromosome linked inhibitor of apoptosis protein related factor 1 (XAF1) is a newly discovered apoptotic promoter, which is located on chromosome 17p13.2.[ X chromosome linked inhibitor of apoptosis protein inhibits apoptosis by binding to caspases and inhibiting the function of caspases.[ XAF1 is widely expressed in all embryonic and adult normal tissues, but it is low in cancer cell lines and tumor tissues. It was found that apoptotic nuclei could be seen in the muscle tissue of dermatomyositis, and the expression of apoptotic factors was upregulated, suggesting that apoptotic cells might be involved in the pathological process of dermatomyositis.[ XAF1-related excessive apoptosis may exist in dermatomyositis, which may be one of the causes of dermatomyositis.[

NT5E gene is located on chromosome 6 q14-q21, and its expression product NT5E is a membrane nuclease anchored by GPI with a size of about 69 kDa.[ NT5E not only hydrolyzes adenine as a nucleotidase, but also has many physiological functions. Because the activity of NT5E in peripheral blood T cells is 10 times higher than that in thymus, and the activity of NT5E in peripheral B cells is 5 to 6 times higher than that in fetal spleen or umbilical cord blood, it is a marker of lymphocyte maturation.[ In addition, NT5E has cell-specific function during neuronal development and plasticity. Ring et al. found that CD4+ CD25+ Treg cells can downregulate the expression of E-selectin and P-selectin in vascular endothelial cells through NT5E, prevent the adhesion of effector T cells to vascular endothelial cells, and alleviate local skin swelling reaction.[ Another mouse model study showed that CD4+ CD25+Treg cells could inhibit ATP-induced CD62L exfoliation on the surface of CD8+ T cells through NT5E, thereby interfering with the entry and exit of CD8+ T cells into regional lymph nodes and alleviating inflammation.[

UDP-glu-cose ceramide glucosyltransferase (UGCG) can reduce ceramide levels in vivo and induce cells to escape ceramide-induced apoptosis.[ UGCG gene can regulate apoptosis-related pathways. Increased expression of UGCG can induce increased expression of anti-apoptotic gene Bcl-2, decreased expression of Bax and Caspase-3.[ Bcl-2 protein blocks cell necrosis by maintaining mitochondrial membrane potential, which can delay or prevent Fas induced myocyte necrosis in patients with dermatomyositis.[ Bax and caspase-3 can cause the cleavage of PARP and apoptosis of muscle cells, leading to the occurrence of dermatomyositis.[

Guanylate-binding protein 2 (GBP2) is a kind of Lage GTPase protein with molecular weight of 65–67 kDa induced by interferon (IFN), which can mediate cellular active immunity and participate in many biological reactions such as anti-pathogenic microorganisms.[ GBP2 can increase the secretion of interleukin (IL)-6, IL-12, tumor necrosis factor-α (TNF-α), and other inflammatory factors, which has a strong immune promoting effect in the pathogenesis of dermatomyositis.[

Toll like receptors (TLR) are important proteins involved in non-specific immunity (natural immunity).[ TLR3 specifically recognizes ds-RNA, the intermediate product of viral replication, to activate NF-кB and IFN-β precursors.[ Doyle et al confirmed that anti TLR3 monoclonal antibody can inhibit the production of IFN-β in fibroblasts.[ TLR3 pathway can enhance the expression of inflammatory cytokines such as IL-1β and IL-6, amplify the local inflammatory response of patients with dermatomyositis, and aggravate the disease.[ TLR2, TLR4, and TLR9 are elevated in muscle tissue of patients with dermatomyositis, and can activate TNFα, promote the transcription of inflammatory factors.[

The DDX58 gene, also known as RIG-I, recognizes viral RNA and activates innate antiviral responses by inducing type I interferon and inflammatory cytokines.[DDX58 gene plays an important role in activating NF-κ B pathway and plays an important role in the pathogenesis of dermatomyositis.[ It has been found that DDX58 gene is significantly up regulated in psoriasis patients and psoriasis like mice.[

Signal transducer and activator of transcription 1 (STAT1) is involved in many basic biological processes, such as cell proliferation, apoptosis, angiogenesis, and immune response.[ Type I interferon can participate in the myopathy of patients with dermatomyositis through STAT1 signaling pathway.[ IFNα and IFNβ can promote the expression of type I interferon induced gene by binding with IFNAR1 and IFNAR2 to activate STAT1 signaling pathway, leading to muscle tissue damage.[ Recently, it has been confirmed that STAT1 expression in the muscle tissue of patients with dermatomyositis is significantly higher than that of the normal control group.[

Guanylate binding protein 1 (GBP1) is a kind of monomer engine protein composed of 593 amino acids.[ There is no expression of GBP1 in normal human skin, but when the skin is invaded by inflammation, the expression of GBP1 is significantly increased and the aggravation of inflammatory diseases is inhibited, GBP1 may be a potential early-warning signal of inflammation.[ In endothelial cells, GBP1 can be strongly induced by IFN-γ, IFN-α, IFN-β, TNF-α, ILs, and other inflammatory factors, and inhibit the proliferation and invasion of endothelial cells in inflammatory diseases.[

Phospholipid scramblase 1 (PLSCR1) has many functions including signal transduction, pathogen invasion/absorption, intracellular transport, secretion, and endocytosis. PLSCR1 can promote the activation of Caspase-3 and protein kinase C and increase the apoptosis induced by etoposide.[ The main pathological reaction in the course of dermatomyositis is not only degeneration and necrosis of muscle tissue, but also closely related to apoptosis.[ PLSCR1 induced apoptosis plays an important role in dermatomyositis.[

Oligoadenylate synthetase 3 (OAS3) is an important antiviral protein, which is induced by IFN, and plays a key role in regulating protein synthesis and immune response.[ It has been confirmed that OAS3is related to Sjogren syndrome complicated with hepatitis B virus infection.[ At present, there are few studies on the relationship between OAS3 and dermatomyositis, which need further experimental verification.

Anti-solid acid nuclear protein of 100,000 antibodies (SP100) is one of the antinuclear antibodies.[ SP100 has been confirmed to be involved in viral infection, self-ubiquitination regulation, and play a role in interferon, p53, and other signaling pathways.[ SP100 is clinically associated with autoimmune hepatitis.[ Clinically, it often occurs before liver injury in patients with dermatomyositis.[

Immunoglobulin κ light chain (IGK) gene is a kind of highly conserved housekeeping gene, which is produced and secreted only by plasma cells under normal conditions.[ It is highly expressed in some immune related diseases such as multiple myeloma, leukemia, lymphoma, autoimmune diseases, and infectious diseases.[ At present, there are few studies on the role of IGK in the pathogenesis of dermatomyositis. It is suggested that the overexpression of IGK may be involved in the pathogenesis of dermatomyositis by activating NF-κB pathway.

Radical S-adenosyl methionine domain-containing protein 2 (RSAD2) is an endoplasmic reticulum associated virus suppressor gene mediated by interferon I.[ RSAD2 promotes Th0 cells to differentiate into Th2 cells, and then induces and regulates B lymphocyte activation and autoantibody secretion.[ RSAD2 has been found to be related to a variety of autoimmune diseases, such as rheumatoid arthritis, systemic lupus erythematosus, ankylosing spondylitis, and so on.[

Conclusions

We retrieved several microarray datasets on the GEO website and got 4 microarray datasets. Then, we applied R software and statistics to merge these datasets. We have recognized 4097 candidate DEGs, which may participate in the pathogenic process of dermatomyositis. By enriching the GO and KEGG pathways, we got that DEGs were primarily enriched in the toll-like receptor signaling pathway, cytosolic DNA-sensing pathway, RIG-I-like receptor signaling pathway, complement and coagulation cascades, arginine and proline metabolism, phagosome signaling pathway in dermatomyositis, which offer a foundation for exploring the pathogenic processes of dermatomyositis. We ultimately formed a PPI network of DEGs in dermatomyositis and searched many key genes encoding proteins in the network that participate in the pathological mechanism of dermatomyositis in the form of molecular populations. Following research of this network would be instructive for cognition the relationship between DEGs. These researches ameliorate our cognition of the molecular mechanism of dermatomyositis and the appearance and development of the underlying pathogenic factors. Our study has considerable clinical meaning for the early diacrisis and remedy, in addition to the warding off, of dermatomyositis and obtains available targets for the medication of dermatomyositis. Nevertheless, further in vivo and in vitro verification experiments are demanded to demonstrate the effect of the recognized genes bound up with dermatomyositis.

Author contributions

Wei Liu developed the concept and algorithm. Yuan-Hao Wu and Wen-Jia Zhao integrated the databases, implemented the codes, and established the webservice. Wei Liu, Yuan-Hao Wu, and Wen-Jia Zhao analyzed and evaluated the results. Yuan-Hao Wu and Wen-Jia Zhao drafted the manuscript. All authors have read and approved the final version of the manuscript for publication.