Abstract

Objective. Systemic lupus erythematosus (SLE) is a chronic autoimmune disease that can affect all organs in the body. It is characterized by overexpression of antibodies against autoantigen. Although previous bioinformatics analyses have identified several genetic factors underlying SLE, they did not discriminate between naive and individuals exposed to anti-SLE drugs. Here, we evaluated specific genes and pathways in active and recently diagnosed SLE population. Methods. GSE46907 matrix downloaded from Gene Expression Omnibus (GEO) was analyzed using R, Metascape, STRING, and Cytoscape to identify differentially expressed genes (DEGs), enrichment pathways, protein-protein interaction (PPI), and hub genes between naive SLE individuals and healthy controls. Results. A total of 134 DEGs were identified, in which 29 were downregulated, whereas 105 were upregulated in active and newly diagnosed SLE cases. GO term analysis revealed that transcriptional induction of the DEGs was particularly enhanced in response to secretion of interferon-γ and interferon-α and regulation of cytokine production innate immune responses among others. KEGG pathway analysis showed that the expression of DEGs was particularly enhanced in interferon signaling, IFN antiviral responses by activated genes, class I major histocompatibility complex (MHC-I) mediated antigen processing and presentation, and amyloid fiber formation. STAT1, IRF7, MX1, OASL, ISG15, IFIT3, IFIH1, IFIT1, OAS2, and GBP1 were the top 10 DEGs. Conclusions. Our findings suggest that interferon-related gene expression and pathways are common features for SLE pathogenesis, and IFN-γ and IFN-γ-inducible GBP1 gene in naive SLE were emphasized. Together, the identified genes and cellular pathways have expanded our understanding on the mechanism underlying development of SLE. They have also opened a new frontier on potential biomarkers for diagnosis, biotherapy, and prognosis for SLE.

1. Introduction

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease that can affect all organs in the body. It is characterized by hyperactivation of the immune system, resulting in abnormal production of autoantibodies, particularly antinuclear antibodies [1]. The clinical manifestations and severity of the disease are diverse, but are usually life-threatening. Currently, the global incidence of SLE is about 40/100,000 [2]. Women in their reproductive years are the most susceptible group for the development of SLE. The pathogenesis of SLE is complex, and the whole pathogenesis cascade is mediated by immune disorders. A lot of effort has been made to study the pathogenesis of SLE. In effect, hereditary, environmental, and hormonal dynamics have been identified as three major factors in SLE etiopathogenesis [3]. However, genetic factors are the major variables that mediate the pathogenesis of SLE [4]. Nonetheless, specific genes and pathways mediating the pathogenesis of the disease are not well understood. This in effect has undermined progress in treatment and prognosis prediction of SLE.

High-throughput sequencing is a new omics technology that has opened up a new frontier in disease genomics [5]. In the recent past, the combination of microarray technology and bioinformatics has facilitated the understanding of molecular mechanisms and reliable diagnostic and therapeutic targets for diseases. Through these techniques, many new genes involved in complex human diseases such as cancer and autoimmune complications have also been discovered [6, 7]. Several studies have identified disease-related genes and functional pathways by analyzing the expression profile between SLE and healthy controls [812]. However, in the majority of the studies, SLE patients have a history of exposure to anti-SLE drugs, a variable that can cause considerable transcriptional changes [13]. Consequently, findings from such studies may mask some features that would have otherwise been apparent. Therefore, in order to eliminate the interference of drugs on transcriptomes, the GSE46907 dataset from NCBI-GEO encompassed 5 newly diagnosed SLE patients without history of exposure to anti-SLE drugs and 5 controls was selected. Peripheral monocyte transcript samples were then obtained from the identified individuals for evaluation. Differentially expressed genes between the two groups were identified and analyzed based on the principles of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. Protein-protein interaction (PPI) network of the DEGs and hub genes was constructed using Cytoscape. Our study is aimed at identifying dominant genes associated with the pathogenesis of naive SLE in view of providing new targets and knowledge for the diagnosis and treatment of the disease.

2. Materials and Methods

2.1. Affymetrix Microarray Data

GSE46907 gene expression profile for peripheral blood monocytes both from SLE patients and controls was downloaded from the Gene Expression Omnibus database (available at https://www.ncbi.nlm.nih.gov/geo/) [14]. This dataset included gene expression profiles for 5 pediatric SLE patients naive for anti-SLE drugs. The matrix was based on GPL96 platforms (HG_U133A; Affymetrix Human Genome U133A Array).

2.2. Identification of Differentially Expressed Genes (DEGs)

DEGs between SLE patients and controls were identified using the LIMMA package downloaded from Bioconductor (http://www.bioconductor.org/packages/release/bioc/html/limma.html) in R software (version 3.4.0; https://www.r-project.org/). Student’s -test in the LIMMA package was used to analyze significant difference between the groups. DEGs with adjusted and were considered to be statistically significant. The genetic profiles for the peripheral blood monocytes were visualized using a heat map of hierarchical cluster analysis (HCA) and volcano plot.

2.3. Gene Ontology and Pathway Enrichment Analyses of DEGs

GO (Gene Ontology) terms [15] and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analyses are two widely used approaches in understanding high-level utilities of biological systems in bioinformatics research [16]. Metascape (http://metascape.org) was used to perform functional enrichment analysis. This is an online analytical tool integrated with several ontology sources including the KEGG pathway, GO biological processes, Reactome gene sets, canonical pathways, and CORUM [17]. Pathway and process enrichment for each DEGs gene were analyzed for statistical significance in each biological process. Only terms with , a minimum count of 3 and an enrichment were considered significant. Genes were clustered according to their pathways. Bar charts were used to visualize the results for the GO terms and KEGG pathway enrichment analyses of the DEGs.

2.4. PPI Network Analysis and Hub Gene

The STRING online database (http://string-db.org) predicted the protein-protein interaction (PPI) [18] between DEGs. A combined was considered significant. The PPI network was visualized using Cytoscape software (version 3.4.0; http://www.Cytoscape.org). It simplified the complex network of PPI. Nodes with a high degree from the Network Analyzer plugin (cytoHubba) [19] were classified as hub genes. These are genes potentially mediating key physiological regulatory functions in SLE.

3. Results

3.1. DEGs

This study is aimed at identifying DEGs between normal and individuals with SLE. Given that exposure to anti-SLE drugs can interfere normal transcription, the selected SLE patients in GSE46907 were healthy, newly diagnosed with the disease, and naive for autoimmune medications. The expressed transcripts in SLE peripheral blood monocytes were obtained and analyzed. Compared with healthy control, a total of 134 DEGs were identified in individuals with SLE, in which 29 were downregulated, whereas 105 were upregulated. The heat map and volcano plot for the DEGs are shown in Figures 1 and 2, respectively. Further details on the DEGs between the two groups are shown in Supplementary Materials in Table 1.

3.2. Gene Ontology and Pathway Enrichment Analyses of DEGs

Signaling pathway enrichment analysis of DEGs was performed using Metascape. The top 20 clusters for enriched sets are shown in Figure 3. For GO term enrichment analysis, expression of DEGs was mainly enhanced in response to secretion of interferon-γ and interferon-α, regulation of cytokine production, regulation of innate immune response, cellular viral responses, and regulation of I-kappaB kinase/NF-kappaB signaling. KEGG analysis revealed that the expression of DEGs was mostly enhanced in interferon signaling, antiviral responses by IFN-stimulated genes, and class I major histocompatibility complex (MHC-I) mediated antigen processing and presentation and formation of amyloid fibers. Specific genes associated with each pathway are shown in Table 1.

3.3. PPI Network Analysis and Hub Genes

A PPI network for the DEGs was constructed on the STRING website. The results were then visualized using Cytoscape (Figure 4). Generally, there were 80 nodes and 947 edges in the PPI network. Most DEGs in the PPI network were upregulated, with only six found to be downregulated. The top 10 genes with the highest degrees in the PPI network were identified using a plugin tool (cytoHubba) in the Cytoscape. They were as follows: signal transducer and activator of transcription 1 (STAT1, degree: 54), interferon regulatory factor 7 (IRF7, degree: 50), MX dynamin-like GTPase 1 (MX1, degree: 49), 2-5-oligoadenylate synthetase-like (OASL, degree: 49), ISG15 ubiquitin-like modifier (ISG15, degree: 48), interferon-induced protein with tetratricopeptide repeats 3 (IFIT3, degree: 47), interferon induced with helicase C domain 1 (IFIH1, degree: 47), interferon-induced protein with tetratricopeptide repeats 1 (IFIT1, degree: 46), 2-5-oligoadenylate synthetase 2 (OAS2, degree: 46), and guanylate binding protein 1 (GBP1, degree: 46). The interactions of the top 10 hub genes are shown in Figure 5. Detailed information on the listed genes and their functions is shown in Table 2.

4. Discussion

SLE is a heterogeneous autoimmunity disease characterized by overexpression of antibodies against autoantigens [20]. Genetic factors play critical roles in the development of SLE [21]. The role of single and cluster genes in the development and pathogenesis of SLE is described [22]. Recent technological advancements such as GWAS and microarray have expanded our understanding on gene expression, including those associated with the pathogenesis of SLE [23]. The previous studies however focused on the expression signature of blood IFN-α [2426]. Even so, the majority of individuals sampled for most of these studies did not discriminate naive and anti-SLE-exposed individuals, a factor that may influence gene transcription profile. This study regularized this limitation. GSE46907 matrix was downloaded from the GEO databases analyzed with integrated bioinformatics methods with a view of identifying reliable genes and pathway associated with SLE. Notably, all patients involved were naive for anti-SLE drug therapy. In the end, 134 DEGs were identified; most of which were upregulated. The upregulated category consisted of genes involved in the activation of the immune system rather than immunosuppression. Among them, 90 DEGs (10 downregulated, 80 upregulated) constituted a PPI network. Further analysis using Cytoscape revealed top 10 hub genes: STAT1, IRF7, MX1, OASL, ISG15, IFIT3, IFIH1, IFIT1, OAS2, and GBP1. Hub genes may mediate the development of SLE; thus, are potential targets and biomarkers for disease treatment and prognosis, respectively.

Previous bioinformatics analysis showed that IFN signature and IFN signaling pathways against viral infections were strongly associated with the pathogenesis and development of SLE. High levels of serum IFN together with overexpression of IFN-inducible genes have been found in individuals with SLE. The level of IFN correlated with the severity of the disease [12]. The role of type I IFN, paticularly IFN-α in the development and pathogenasis of SLE has been described [26]. To a greater extent, findings of this study were consistent with reports from previous studies. First, our analysis revealed that innate IFNs against viral infection were critical in the development and or pathogenesis of SLE. The overall interferon pathway was the most enhanced pathway in the pathogenesis of SLE. The top 10 hub genes identified in our studies are all associated with excessive production of IFN in viral infection and innate immune response (detailed in Table 2). Interestingly, 9 (STAT1, IRF7, MX1, OASL, ISG15, IFIT3, IFIH1, IFIT1, and OAS2) out of the 10 hub genes have been reported to be associated with SLE (except GBP1) previously [4, 11, 27, 28]. Second, class I MHC-mediated antigen processing and presentation and I-kappaB kinase/NF-kappaB signaling pathway were also critical in our study. One study found that class I MHC was strongly expressed among SLE individuals [29]. In 16/6ld-immunized mouse models, MHC molecules mediated induction of SLE, whereas repressing transcription of class I MHC prevents clinical manifestations of SLE [30]. Interestingly, the expression of MHC I molecules is enhanced by type I IFN [31]. Here, the interplay between type I IFN and MHC I molecules enhances peptide recognition by cytotoxic T cells. This intern promotes the induction of cell-mediated SLE. Similarly, the NF-κB signaling pathway also participates in regulating the immune and inflammatory responses in individuals with SLE [32]. Therefore, inhibiting the NF-κB signaling pathway could ameliorate symptoms of lupus nephritis or lupus-prone MRL/lpr mice [33, 34]. Because our study supports previous findings, it further strengthens our understanding on the pathogenesis of SLE.

Remarkably, our findings underscored the importance of IFN-γ in the pathogenesis of innate SLE, as opposed to IFN-α which has been emphasized in the previous bioinformatics analysis. IFN-γ is the only member of type II interferons, mainly produced by Th1-type T cells and NK cells [35]. One study found that in SLE patients, increase in the secretion of IFN-γ and IFN-α was positively correlated, and that types I and II IFN partially share signal pathways and target genes [36]. This crossinteraction facilitates coordination of these IFNs in executing specific functions in the immunopathogenesis of the SLE. Besides this, Groettrup et al. found that IFN-γ crossinteracts with MHC molecules, enhancing the quantity, quality, and repertoire of peptides bound by both class I and class II MHC [37]. IFN-γ can activate the transcription of both class I and class II MHC molecules, which later contribute to the development and severity of SLE. In one study, it was found that secretion of IFN-γ was elevated before the production of IFN-α and lupus-associated autoantibodies [35]. Because IFN-γ links innate and adaptive immunity, it can drive secretion of both type I IFN and BLyS (B lymphocyte stimulator) [38]. Chodisetti et al. found that IFN-γ is indispensable for TLR7-promoted development of autoreactive B cells and systemic autoimmunity [39]. IFN-γ can also regulate the production of Ig immunoglobulins and B cell class switching [40]. One study found that IFN-γ can promote aberrant activation and imbalanced polarization of macrophages [41], which also play an important role in SLE [42]. From a bioinformatics perspective, our study solidified the role of IFN-γ in the development of SLE. Based on previous bioinformatics analyses, we hypothesized that use of drugs such glucocorticoids may mainly reduce the initial IFN-γ signature. According to our analyses, one new hub gene (GBP1) that participates in the innate SLE was discovered. GBP1 is a critical interferon-stimulated gene (ISG). Interferon stimulation, particularly IFN-γ, activates overexpression of this gene [43]. Some studies suggest that GBP1 can be used as a marker for IFN-γ response [44]. GBP1 gene product can inhibit cellular proliferation and apoptosis in the early cellular inflammatory response [45]. GBP1 is also essential in intracellular defense in varied settings such as in oxidative or other cellular damages [46]. Furthermore, special attention has been put on the role of GBP1 in cancer pathology [47]. GBP1 is a potential target for cancer inflammatory responses [46]. In our case, there are limited researches on the relationship between GBP1 and SLE. However, our study underscores the significant role played by GBP1 and IFN-γ in the early onset of SLE; thus, the role of GBP1 on SLE should be further investigated.

5. Conclusion

In conclusion, this study strengthens the critical role played by IFN signatures and IFN signaling pathway in SLE. Different from previous bioinformatics analyses, we analyzed the transcription profile of treatment of the naive SLE patients. We found that IFN-γ and its response GBP1 genes play critical roles at the onset of SLE. These findings gave a new insight in the pathogenesis of the disease. It provided valuable markers or therapeutic targets for GBPI. However, this study had several limitations. First, the sample size used in our analysis was small, reducing the reliability of our findings. Second, given that some studies found that the expression of IFN-γ in peripheral blood mononuclear cells (PBMCs) in individuals with SLE was high and the dataset selected for this research measured the transcription profiles of peripheral blood monocyte, our findings may be biased. Therefore, our results require further validation.

Data Availability

The data used for analysis in this study are available from the Gene Expression Omnibus database freely.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this article.

Acknowledgments

We gratefully thank many researchers for providing technical assistance and all the patients who were enrolled in the study.

Supplementary Materials

Supplementary Materials of Table 1: the results of differentially expressed genes (DEGs) of healthy and newly diagnosed SLE population. (Supplementary Materials)