key: cord-0018686-hmtuisdu
authors: Wang, Yanfeng; Ma, Qian; Huo, Zhenghao
title: Identification of hub genes, pathways, and related transcription factors in systemic lupus erythematosus: A preliminary bioinformatics analysis
date: 2021-06-25
journal: Medicine (Baltimore)
DOI: 10.1097/md.0000000000026499
sha: 6fe165650a98f8f67bddfca925cfd355ae4c8833
doc_id: 18686
cord_uid: hmtuisdu

BACKGROUND: Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by multiple organ damage and the production of a variety of autoantibodies. The pathogenesis of SLE has not been fully defined, and it is difficult to treat. Our study aimed to identify candidate genes that may be used as biomarkers for the screening, diagnosis, and treatment of SLE. METHODS: We used the GEO2R tool to identify the differentially expressed genes (DEGs) in SLE-related datasets retrieved from the Gene Expression Omnibus (GEO). In addition, we also identified the biological functions of the DEGs by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. Additionally, we constructed protein–protein interaction (PPI) networks to identify hub genes, as well as the regulatory network of transcription factors related to DEGs. RESULTS: Two datasets were identified for use from the GEO (GSE50772, GSE4588), and 34 up-regulated genes and 4 down-regulated genes were identified by GEO2R. Pathway analysis of the DEGs revealed enrichment of the interferon alpha/beta signaling pathway; GO analysis was mainly enriched in response to interferon alpha, regulation of ribonuclease activity. PPIs were constructed through the STRING database and 14 hub genes were selected and 1 significant module (score = 12.923) was obtained from the PPI network. Additionally, 11 key transcription factors that interacted closely with the 14 hub DEGs were identified from the gene transcription factor network. CONCLUSIONS: Bioinformatic analysis is an effective tool for screening the original genomic data in the GEO database, and a large number of SLE-related DEGs were identified. The identified hub DEGs may be potential biomarkers of SLE.

Systemic lupus erythematosus (SLE), a complex autoimmune disease that often involves multiple organ systems, is characterized by the production of a large number of autoantibodies, extensive deposition of immune complexes, and abnormalities of the innate and adaptive immune responses. [1] At present, the pathogenesis and detailed etiology of SLE remains unclear; as such, uncovering the potential molecular mechanisms and identifying reliable biomarkers is critical for early diagnosis and effective treatment, genetic deficiency in C1q is a strong susceptibility factor for SLE. [2] Current evidence suggests that mutations in single genes, such as TREX1, RNASEH2B, ADAR, IFIH1, and SAMHD1 contribute to the development of SLE, can induce chronic type I interferon responses, and are associated with SLE-like lesions. Nucleic acid degradation is associated with monogenic forms of lupus-like disease. [3, 4] Additionally, thrombosis is often induced in the process of SLE. Targeted platelet therapy can help manage the disease group most affected by platelet dysfunction and improve the quality of life of these patients. [5] It has been reported that radiomonoclonal antibodies have been developed for tumor diagnosis and targeted radiotherapy. [6] Human studies provide evidence that SLE is a polygenic disease in both patients with SLE and SLE mouse models, and involves defects in many cellular signaling pathways that lead to increased survival of self-reactive cells. [7] [8] [9] [10] High-throughput sequencing and microarray technology have been widely used to screen genes critical to the initiation and progression of various diseases, including SLE. Despite the identification of hundreds of differentially expressed genes (DEGs) through various methods, no suitable biomarkers have been reported, which may be a result of varying experimental conditions and study populations. In addition, early studies on genetic factors in SLE focused mainly on the contribution of single genes. Transcription factors (TFs) also play an important role in the development of SLE; however, early studies on TFs were done by microarray have not been reported for SLE. Thus, the application of bioinformatics methods to re-analyze currently available transcriptome data to gain new insights is warranted to better understand the molecular mechanisms of SLE and to identify reliable biomarkers for disease diagnosis. At present, an increasing number of studies have reported that the interaction of multiple genes in the development of SLE is closely related to the incidence of disease and the involvement of multiple organ systems. Therefore, mining the available databases for potential biomolecular mechanism of SLE is of great value for disease screening and diagnosis, and may provide a new way to explore drug treatment targets. With the rapid development of gene chip and RNA sequencing technology, bioinformatic analysis can play an important role in the screening of candidate biomarkers in the research of many diseases. [11] In this study, we aimed to identify the DEGs in SLE and the TF regulatory network associated with the DEGs from publicly available SLE-related datasets. Specifically, we downloaded 2 available databases, GSE50772 and GSE4588, and screened the DEGs by GEO2R. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was also performed for the functional analysis of the DEGs. A protein-protein interaction (PPI) network of the DEGs was constructed by STRING and was visualized by Cytoscape software. Furthermore, a gene-TF regulatory network of hub DEGs was constructed to assess the interactions between the TFs and hub DEGs. The workflow used in the present study is presented in Figure 1 . In our study, 38 SLE DEGs and TFs of key genes were identified. We found that TRIM28 and EED functioned as transcription factors of the IFI gene family and regulated interferon inducible 44 like (IFI44L). Additionally, IFITM3 was highly expressed in SLE. Our study provides a new methodology to study the molecular mechanisms of SLE from the genomic level and explore potential early diagnostic markers and drug treatment targets. Bioinformatic analysis can provide new clues and core data for the research of DEGs and microRNAs (miRNAs) of diseases.

A microarray search of GSE50772 and GSE4588 in SLE was conducted in the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/) with the keywords "SLE" or "Systemic lupus erythematosus." The entry type was restricted to "series" and the organism was filtered by "Homo sapiens." The criteria for inclusion were as follows:

(1) patients diagnosed with SLE, (2) at least 20 healthy and SLE samples from peripheral blood mononuclear cells (PBMCs), and (3) the research subjects can only be humans, not animal models.

The GSE50772 database contains 81 samples, including 20 control samples and 61 SLE patient samples.

The GSE4588 database contains 49 samples, including 23 control samples and 18 SLE patients. The remaining 8 rheumatoid arthritis samples were not selected in this study. No ethical approval was needed as this study is solely based on bioinformatic analysis.

Using the GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r) online data analysis tool to compare the gene expression of the healthy control and SLE group, the original data was downloaded, the DEGs were screened, and the screening standard was set to P < .05 and jlog 2 FCj > 1.0. After screening, the DEGs were regarded as the genes with significant differences in expression, and Sanger Box software (version 1.0.9) was used to draw the volcano plot of the DEGs. Heat maps were drawn using the web-based Morpheus tool (https://software.broad institute.org/morpheus/). GO analysis is a common method for large-scale functional enrichment research. GO database can describe the standardized gene products from biological process, molecular functions, and cellular components. KEGG analysis is a widely used database which stores a large number of data about the genome, biological pathways, diseases, chemicals, and drugs. We used Metascape (http://metascape.org/) as a database annotation and visualization tool.

We constructed PPI networks to analyze protein interaction information and further identify hub genes. The results showed that the interaction between protein and protein was significant. Then, the PPI network retrieved from the STRING database was visualized using the Cytoscape software (version 3.7.2, http://www. cytoscape.org/). Then, according to the interaction pair information, the PPI network was constructed and visualized by Cytoscape software. The Cytoscape plug-in Network Analyzer was used for further analysis, and the topological properties of the PPI network and node degree were calculated to search for hub genes from the PPI network. Subsequently, Molecular Complex Detection (MCODE) analysis in Cytoscape was performed to screen the significant modules of PPI network. A degree cutoff = 2, node score cutoff = 0.2, and K-Core = 2 were set in the advanced options.

Network Analyst, an online visual tool based on network gene expression profile analysis, can comprehensively analyze the functional pathway of hub genes and the effect of TFs on gene expression. In this study, we used the network analysis database to predict the TFs of hub genes, and constructed the TF regulatory network of genes, used Cytoscape software and Network Analyst to visualize the results (http://www.networka nalyst.ca/faces/home.xhtml).

Using the GEO2R tool to analyze the sample data from GSE50772 and GSE4588 in the GEO public database, we constructed a volcanic plot (Fig. 2 ) of the expression of DEGs identified in GSE50772 and GSE4588 using the Sanger Box software (version 1.0.9). Each color point represents an upregulated or down-regulated gene, wherein green represents the down-regulated genes, red represents the up-regulated genes, and black indicates the standard based on no differential gene expression. The selection criteria was set as a P value < .05 and jlog 2 FCj > 1.0. In total, 307 up-regulated and 281 down-regulated genes were identified from GSE50772, and103 up-regulated and 24 downregulated genes were identified from GSE4588. The top 100 important genes from the heat maps of GSE50772 and GSE3588 of DEGs are shown in Figure 3 . Cluster analysis showed that there were differences in the DEGs between the control group and the SLE group. Venn diagram analysis showed that 34 upregulated genes and 4 down-regulated genes overlapped in GSE50772 and GSE12452 (Fig. 4) . These results are also listed by significance in the gene table (Table 1 ). The 38 overlapping genes were identified as candidate genes for further analysis.

GO functional analysis and KEGG pathway enrichment analysis of the candidate DEGs were performed based on the Metascape database. Terms or pathways with a P-value < .01, min overlap www.md-journal.com genes = 3, and min enrichment factor > 1.5 were set as the cutoff criteria. There were 4 terms and 2 pathways involved in the enrichment analysis of the DEGs, indicating that these DEGs were mainly enriched in response to virus infection, regulation of ribonuclease activity, response to interferon alpha, and mito-chondrion organization. A previous study by Monti and Montecucco found that response to virus and response to interferon alpha is strongly correlated with SLE progression. [12] In addition, these enrichment terms are closely connected with each other and are integrated into a complete network (Fig. 5b) . The results of the KEGG enrichment analysis are also shown in Fig. 5a . Two significantly enriched pathways, interferon alpha/ beta signaling, and ISG15 antiviral mechanism were identified as being associated with the DEGs, suggesting that these signaling pathways have important significance for the pathogenesis of SLE. [24] 3.3. PPI network construction and identification of hub genes

The candidate DEGs were uploaded to the STRING online tool to obtain PPI information, and the results of the PPI analysis were visualized (Fig. 6a ) and analyzed with Cytoscape (Fig. 6b) . A total of 33 of the 38 shared DEGs were filtered into the DEG PPI network complex, which contained 34 nodes and 97edges (criteria of filtering degree ≥10).

Based on the MCDODE analysis of cell morphology, 4 important modules were obtained from the PPI network. We selected the most important module (MCODE score = 12.923) for further analysis (Fig. 6b) . This module includes 14 nodes (OAS1, IFI27, RSAD2, IFIT2, ISG15, OAS3, IFITM3, IFIT1,  IFIT3 , OASL, Myxovirus resistance protein 1, IFI44L, IFI6, and IFI44) and 84 edges, and these genes were the hub DEGs we identified. Functional enrichment analysis showed that this PPI module was significantly enriched in interferon alpha/beta signaling, response to virus, and mitochondrion organization.

For the 14 central genes identified, we constructed a TF regulatory network of 57 interacting pairs, 14 genes, and 9 TFs (Fig. 7) . For example, IFIT1 was found to be regulated by 11 TFs, IFI6 was regulated by 6TFs, IFI44L was regulated by 2TFs, and RSAD2 was regulated by 3 TFs. In addition, we also found that a number of TFs can regulate more than 1 hub gene; 9 TFs were identified as having a regulatory gene connectivity of ≥1 network, suggesting that these TFs have close interactions with the hub DEGs (Table 2) . For example, IRF1 was predicted to regulate IFIT2, IFI6, ISG15, IFIT3, IFIT1, OAS3, IFITM3, and RSAD2; hepatoma derived growth factor (HDGF) was predicted to regulate both IFI6, ISG15, OAS3, and IFITM3; and EED was predicted to regulate IFI44L, IFIT3, and IFIT2. Leung et al found that prolactin can activate IRF1, which is a TF involved in SLE and interferon response. [13] 4. Discussion SLE is a common autoimmune disease, and common infections have become the leading cause of death in patients with SLE.

Recently, an increasing number of studies have used microarray data profiling to elucidate the pathogenesis of SLE. [14, 15] The GEO database is a large database established and maintained by the National Biotechnology Information Center of the United States. It contains high-throughput genomic data information uploaded by many research institutions around the world, but a large number of data has not been fully mined. [16] Such a situation also exists in the study of SLE-related genomic data. Due to the limitations of research purposes and sample inclusion independence, a large amount of genomics data has not been vertically or horizontally compared. [16, 17] Effective biomarkers can provide accurate treatment options and prognosis for diseases. In the diagnosis of multiple sclerosis, surfaced-enhanced Raman spectroscopy characterizes tissues and cells by recognizing molecular chemical composition, which has significant differences or similarities in potential tumor markers or diagnosis. [18, 19] Based on this, using bioinformatics methods to analyze the gene chip and RNA sequencing data of clinical studies, we can screen the core markers of diseases and identify the potential biomarkers critical for the process of disease development. Gene therapy and molecular targeted therapy have become increasingly important in the clinical treatment of diseases. In this study, the gene chip expression profiles from 2 SLE studies (GSE50772, GSE4588) were analyzed using bioinformatic analysis, and a total of 38 DEGs were obtained, including 34 up-regulated genes and 4 down-regulated genes. According to the network topology parameters, 34 hub DEGs were obtained from the PPI network. Furthermore, 1 significant module with an MCODE score = 12.923 was screened from the PPI network, consisting of 14 hub DEGs. A number of early studies also applied bioinformatic analysis to SLE data to identify candidate genes as potential biomarkers for predicting the prognosis of patients with SLE. Yang et al analyzed the expression profiles of 3 SLE gene chips (GSE72509, GSE20864, and GSE39088) in the GEO database, and identified 4 hub genes, including RPL26L1, FBXW11, FOXO1, and SMAD7, [20] that play key roles in the pathogenesis and development of SLE. Additionally, the RIG-I-like receptor signaling pathway, antigen processing and presentation pathway, and p53 signaling pathway may be closely implicated in SLE pathogenesis. [20] In our study, we obtained 14 hub DEGs, all of which were all up-regulated in SLE. Functional and pathway enrichment analysis showed that the DEGs were mainly involved in functions related to responses to virus, regulation of ribonuclease activity, response to interferon alpha, and mitochondrion organization, and in the interferon alpha/beta signaling and ISG15 antiviral mechanism pathways. This has been proved to be the key to the treatment of SLE. In addition, through PPI analysis, these 14 hub DEGs were found to have a close relationship in SLE, and they are involved in key pathways in the biological process of SLE. A large number of studies have shown that the expression level of interferon-induced genes is closely related to the pathogenesis of SLE. IFI6, IFI27, IFI44, IFI44L, IFIT1, IFIT2, IFIT3, and IFITM3 are interferon (IFN)-a-related genes. For example, IFI44L (interferon inducible 44-like) is an interferon-induced gene that is located on the same chromosome as the known antiviral gene. [23] The expression of the IFI44L gene can be upregulated 33.33 times by IFN-a. Both IFI44L and IFIT3 are involved in the type I IFN signaling pathway, indicating that the type I IFN signaling pathway is closely related to the pathogenesis of SLE. Zhao et al showed that the level of DNA methylation of IFI44L in peripheral blood can be used to evaluate and diagnose SLE. [21] The specificity and sensitivity of IFI44L promoter methylation in differentiating SLE from healthy controls and other autoimmune diseases, such as rheumatoid arthritis and systemic sclerosis, are significantly better than the biomarkers currently in use. Additionally, many studies have shown that the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon gene (STING) signaling pathway plays an important role in the pathogenesis of SLE. [21] Wang et al found that increased levels of IFIT3 in the monocytes of patients with SLE are correlated with the overactivation of the CGAs/sting signaling pathway. As a novel therapeutic target for SLE, reducing the level of IFIT3 may inhibit the activation of the CGAs/sting signaling pathway, thereby reducing the production of type I IFN and other proinflammatory cytokines. [22] Fan et al reported that IFI44L was significantly up-regulated in patients with SLE, and it was found that differences in the incidence of SLE based on sex were due to estrogen-induction on IFI44L. [24] RSAD2 is an important core gene in the type I interferon signaling pathway, and is involved in viral defense response and negative regulation of viral genome replication in the process of viral infection. However, the mechanism of RSAD2 in the pathogenesis of SLE remains Table 2 Transcription factor analysis of key differentially expressed genes.

Genes Count   IRF1  IFIT2, IFI6, ISG15, IFIT3, IFIT1, OAS3, IFITM3, RSAD2  8  HDGF  IFI6, ISG15, OAS3, IFITM3  4  MTA2  OASL, IFIT2, IFI6  3  NR2C2  ISG15, IFIT3, IFI27  3  ZNF644  ISG15, IFITM3, IFIT2  3  WRNIP1  IFI6, FITM3, IFIT3  3  SMARCE1  IFI6, FITM3, ISG15  3  GTF2F1  IFI6, OASL, ISG15  3  ZNF12  OASL, OAS3, RSAD2  3  TRIM28  IFITM3, RSAD2, IFI44L  3  EED  IFI44L, IFIT3, unclear. Therefore, the study of RSAD2 expression in SLE is required to further elucidate the pathogenesis of SLE. [25, 26] OAS1, OAS2, and OASL are 3 subtypes of 2 0 5 0 -oligoadenylate synthetase (OAS). OAS is a type I IFN inducible enzyme with known antiviral functions. [30] The increase of OAS levels in peripheral blood lymphocytes of patients with SLE was observed in the early 1980s. [27] In previous array studies, we identified 3 subtypes of OAS, namely OAS1, OAS2, and OASL, that were upregulated in the newly diagnosed active lupus cohort. In addition, there is sufficient evidence to suggest that when different types of interferons or different doses of the same type of interferon the stage of effector SLE. [28, 29] The pattern of OAS subtype can provide useful information to distinguish the outbreak and infection of SLE and other studies have shown that OAS is associated with SLE, and several microarray studies of patients with SLE have been rediscovered to participate in the type I interferon pathway. Ye et al investigated the expression of OAS subtypes (OAS1, OAS2, and OASL) in patients with SLE to differentiate disease outbreak and infection. [30] It was found that the expression level of the OAS family, especially OAS2, was significantly higher in patients with SLE than that of healthy people and patients with non-SLE autoimmune disease. [31] ISG15, a member of the Ub1 family, is a 17 kDa secreted protein that has significant sequence homology to ubiquitin. [32] ISG15 may induce the synthesis and secretion of IFN-g from B-cell-depleted lymphocytes. It has been previously reported that ISG15 acts as a cytokine modulator in the immune response. [33, 34] It has also been previously reported that IFNa/IFN-g may induce the activation of the JAK-STAT1 pathway, regulate proliferation and activation of immunocytes, and lead to the abnormal activation of the immune system, which may affect the occurrence and development of SLE. [35] [36] [37] Therefore, ISG15 may be involved in the abnormal immune response of SLE. ISG15 was also found to be highly expressed in patients with SLE and was related to the degree of disease activity before treatment. ROC curve analysis showed that ISG15 may be a new marker for the effective diagnosis of SLE. [38] Myxovirus resistance protein 1 (Mx1) is one of the downstream targets of the type I IFN pathway and an important component of the early innate immune response. [39] Mx1 is a GTPase that inhibits the multiplication replication of several RNA viruses, including influenza virus, by association to viral nucleoprotein into ribonucleoprotein particles. [39, 40] Overexpression of Mx1 promotes cell death triggered by endoplasmic reticulum stress in mouse-derived fibroblasts in the course of viral infection. [41] Mx1 was found in exosomes from bovine uterine glandular epithelial cells and regulates exosome secretion. [42] Moreover, Mx1 has been used as a marker of type I IFN signature in basic/translational studies. [43] [44] [45] However, there have been few studies that analyzed the protein levels of Mx1 in association with more critical and clinically important organ damage. Mx1 levels were up-regulated in the peripheral blood of patients with lupus even when their disease activities were stable. On the other hand, Mx1 was highly expressed in kidneys from patients with Lupus nephritis (LN) before treatment, which decreased after immunosuppressive treatment. These results suggest that Mx1 is a potential marker for the diagnosis of SLE in the peripheral blood and also for the activity of lupus nephritis in the kidney.

The interferon alpha/beta signaling pathway can produce interferon-induced IFI expression, whose expression can indirectly reflect the amount of interferon in serum. The gene expression products induced by IFN-a include OAS1, OAS2, OASL, IFIT1 (interferon-induced protein 1), IFIT4, LY6E, MX, and ISG15. Interferon was the first discovered cytokine and can be divided into 3 types, namely a, b, and g, according to its source, physical and chemical properties, and biological activity. IFN-a is mainly produced by mononuclear phagocytes, IFN-b by fibroblasts, and IFN-g by activated T cells and NK cells. IFN-a and IFN-b are also called type I IFNs, and IFN-g is a type II IFN. Interferon I has antiviral, antitumor, and immunomodulatory effects, but its immunomodulatory effect is weaker than IFN-g. The main function of type II interferons is to regulate the immune response and its antiviral and antitumor effects are weaker than that of type I interferons. Interferons mainly act on lymphocytes and monocytes, and play an increasingly prominent role in the pathogenesis of SLE. Some studies have shown that a large number of IFN-a cells are found in lymph node biopsy samples of patients with active SLE and it is believed that IFN-a cells play an important role in the pathogenesis of SLE. [46] [47] [48] Many authors believe that the imbalance of interferon, especially IFN-a, is a central cause of the immune abnormalities in SLE. [49] [50] [51] [52] [53] [54] [55] Microarray analysis of mononuclear cells in the peripheral blood of patients with SLE has revealed up-regulated expression of many IFN stimulating genes, which constituted a whole IFN signal. [56, 57] TFs are key gene expression regulators and are related to the development and progression of human diseases. In the current study, we also identified a number of TFs that have close interactions with the hub DEGs. Interferon regulatory factor 1 (IRF1) is a member of a family of transcription factors that regulates the expression of interferons. It is believed that the IFN system plays a key role in the pathogenesis of SLE. It is important to note that the overexpression of inflammatory bodies in lupus monocytes is significantly reduced after IRF1 knockout or inhibition. [58] In patients with SLE, long-term exposure to type I interferons causes monocytes to activate inflammation in an IRF-1-dependent manner. Inhibition of IRF-1 may be a new target for the treatment of SLE-related inflammation and organ damage. [59] HDGF is an acidic heparin-binding protein, which can promote the growth of the human hepatoma cell line, HuH7 Syntharin. [60] It is known that HDGF not only promotes the proliferation of endothelial cells, smooth muscle cells, and a variety of malignant tumor cells, but also promotes cell proliferation and migration, angiogenesis, antiapoptosis, and neurotrophic isobiological activity. [61] [62] [63] 

In conclusion, 38 DEGs were identified from GSE50772 and GSE4588, including 34 up-regulated genes and 4 down-regulated genes. The analysis of the biological information in this study focuses on 2 data sets, comprising 43 healthy samples and 79 SLE samples. Therefore, the sample size is limited, and further investigation on a large sample size is needed to support the findings in this study. The 38s DEGs are mainly enriched in viral response, ribonuclease activity regulation, response to interferon alpha, and mitochondrion organization. KEGG enrichment analysis also showed that interferon alpha/beta signaling and ISG15 antiviral mechanism were 2 significantly enriched pathways, which are closely related to the occurrence of SLE. Key TFs, such as IRF1, are also closely related to the 14 hub DEGs in the TF gene network, suggesting that these hub DEGs may be potential biomarkers of SLE. However, further investigation on the biological functions and molecular mechanisms of these genes 

IFN-I mediates dysfunction of endothelial progenitor cells in atherosclerosis of systemic lupus erythematosus

SLE: novel postulates for therapeutic options

Mutations in the gene encoding the 3'-5' DNA exonuclease TREX1 are associated with systemic lupus erythematosus

Impaired receptor editing and heterozygous RAG2 mutation in a patient with systemic lupus erythematosus and erosive arthritis

Beyond thrombosis: the role of platelets in pulmonary hypertension

Anti-angiogenesis therapy of cancer cells using 153Sm-bevasesomab

Molecular aberrations in human systemic lupus erythematosus

Systemic lupus erythematosus

New insights into disease pathogenesis from mouse lupus genetics

Genomic view of systemic autoimmunity in MRLlpr mice

Microarray analysis in drug discovery and clinical applications

Can hydroxychloroquine protect patients with rheumatic diseases from COVID-19? Response to: 'does hydroxychloroquine prevent the transmission of COVID-19?' by Heldwein and Calado and 'SLE, hydroxychloroquine and no SLE patients with COVID-19: a comment' by Joob and Wiwanitkit

Prolactin activates IRF1 and leads to altered balance of histone acetylation: implications for systemic lupus erythematosus

Bioinformatics analysis of differentially expressed gene profiles associated with systemic lupus erythematosus

Identification of alterations in macrophage activation associated with disease activity in systemic lupus erythematosus

Predicting biomedical metadata in CEDAR: a study of Gene Expression Omnibus (GEO)

GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor

Raman-enhanced spectroscopy (RESpect) probe for childhood non-Hodgkin lymphoma

Evaluation of mir-377-3p expression in patients with multiple sclerosis

Bioinformatics identification of key candidate genes and pathways associated with systemic lupus erythematosus

IFI44L promoter methylation as a blood biomarker for systemic lupus erythematosus

Association of abnormal elevations in IFIT3 with overactive cyclic GMP-AMP synthase/stimulator of interferon genes signaling in human systemic lupus erythematosus monocytes

Identification of distinct gene expression profiles in the synovium of patients with systemic lupus erythematosus

Gender differences of B cell signature related to estrogen-induced IFI44L/BAFF in systemic lupus erythematosus

Gene expression analysis reveals novel shared gene signatures and candidate molecular mechanisms between pemphigus and systemic lupus erythematosus in CD4+ T cells

A key anti-viral protein, RSAD2/ VIPERIN, restricts the release of measles virus from infected cells

Interferon-induced 2 0 -5 0 adenylate synthetase in vivo and interferon production in vitro by lymphocytes from systemic lupus erythematosus patients with and without circulating interferon

Protein synthesis-dependent and -independent induction of p69 2 0 -5 0 -oligoadenylate synthetase by interferon-alpha

Transcriptional induction of the p69 isoform of 2',5'-oligoadenylate synthetase by interferon-beta and interferon-gamma involves three regulatory elements and interferonstimulated gene factor 3

Could 2 0 5 0 -oligoadenylate synthetase isoforms be biomarkers to differentiate between disease flare and infection in lupus patients? A pilot study

A T cell gene expression panel for the diagnosis and monitoring of disease activity in patients with systemic lupus erythematosus

ISG15 in the tumorigenesis and treatment of cancer: an emerging role in malignancies of the digestive system

A human 15-kDa IFN-induced protein induces the secretion of IFN-gamma

Network analysis identifies proinflammatory plasma cell polarization for secretion of ISG15 in human autoimmunity

STS-1 promotes IFN-a induced autophagy by activating the JAK1-STAT1 signaling pathway in B cells

Dipyrithione inhibits IFN-gamma-induced JAK/STAT1 signaling pathway activation and IP-10/CXCL10 expression in RAW264.7 cells

Activation of JAK-STAT1 signal transduction pathway in lesional skin and monocytes from patients with systemic Wang et al. Medicine (2021) 100:25 www.md-journal.com lupus erythematosus

Interferon-stimulated gene 15 expression in systemic lupus erythematosus: diagnostic value and association with lymphocytopenia

Human MxA protein: an interferon-induced dynamin-like GTPase with broad antiviral activity

GTP-bound human MxA protein interacts with the nucleocapsids of Thogoto virus (Orthomyxoviridae)

Interferoninducible antiviral protein MxA enhances cell death triggered by endoplasmic reticulum stress

The myxovirus-resistance protein, MX1, is a component of exosomes secreted by uterine epithelial cells

Association of increased interferoninducible gene expression with disease activity and lupus nephritis in patients with systemic lupus erythematosus

Activation of the interferon-alpha pathway identifies a subgroup of systemic lupus erythematosus patients with distinct serologic features and active disease

Local synthesis of interferonalpha in lupus nephritis is associated with type I interferons signature and LMP7 induction in renal tubular epithelial cells

The natural interferon-alpha producing cells in systemic lupus erythematosus

IFN-alpha skews monocyte differentiation into Toll-like receptor 7-expressing dendritic cells with potent functional activities

Functional assay of type I interferon in systemic lupus erythematosus plasma and association with anti-RNA binding protein autoantibodies

Type I interferons (alpha/beta) in immunity and autoimmunity

The interferon-alpha signature of systemic lupus erythematosus

Cytokines as therapeutic targets in SLE

Interferon alpha as a primary pathogenic factor in human lupus

The type I interferon system in the development of lupus

Type I interferon blockade in systemic lupus erythematosus: where do we stand?

Anti-interferon alpha treatment in SLE

Interferon and granulopoiesis signatures in systemic lupus erythematosus blood

Interferon-inducible gene expression signature in peripheral blood cells of patients with severe lupus

Enhanced inflammasome activity in systemic lupus erythematosus is mediated via type I interferon upregulation of interferon regulatory factor 1

Interferon regulatory factor-1 gene deletion decreases glomerulonephritis in MRL/lpr mice

Molecular cloning of complementary DNA for a novel human hepatoma-derived growth factor. Its homology with high mobility group-1 protein

The expression and clinical significance of HDGF in osteosarcoma

Hepatoma-derived growth factor is expressed after vascular injury in the rat and stimulates smooth muscle cell migration

Differential proteomic analysis reveals that EGCG inhibits HDGF and activates apoptosis to increase the sensitivity of non-small cells lung cancer to chemotherapy

We thank Prof. Jing Zhang (Medical College of Yan'an University, China) for the kind assistance. The authors thank AiMi Academic Services (www.aimieditor.com) for English language editing and review services.