Abstract
We aim to identify the common genes, biological pathways, and treatment targets for primary Sjögren’s syndrome patients with varying degrees of fatigue features. We select datasets about transcriptomic analyses of primary Sjögren’s syndrome (pSS) patients with different degrees of fatigue features and normal controls in peripheral blood. We identify common differentially expressed genes (DEGs) to find shared pathways and treatment targets for pSS patients with fatigue and design a protein-protein interaction (PPI) network by some practical bioinformatic tools. And hub genes are detected based on the PPI network. We perform biological pathway analysis of common genes by Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. Lastly, potential treatment targets for pSS patients with fatigue are found by the Enrichr platform. We discovered that 27 DEGs are identified in pSS patients with fatigue features and the severe fatigued pSS-specific gene is RTP4. DEGs are mainly localized in the mitochondria, endosomes, endoplasmic reticulum, and cytoplasm and are involved in the biological process by which interferon acts on cells and cells defend themselves against viruses. Molecular functions mainly involve the process of RNA synthesis. The DEGs of pSS are involved in the signaling pathways of viruses such as hepatitis C, influenza A, measles, and EBV. Acetohexamide PC3 UP, suloctidil HL60 UP, prenylamine HL60 UP, and chlorophyllin CTD 00000324 are the four most polygenic drug molecules. PSS patients with fatigue features have specific gene regulation, and chlorophyllin may alleviate fatigue symptoms in pSS patients.
1. Introduction
Primary Sjögren’s syndrome (pSS) is an all-body autoimmune disease that mainly affects middle-aged women [1]. The main clinical feature of the disease is dryness of the mouth and eyes, and the pathophysiology is characterized by focal lymphocyte infiltration in exocrine glands [2, 3]. Fatigue is commonly seen in pSS patients as an extraglandular manifestation and closely links with poor life quality [4–6]. Fatigue affects approximately 70% of pSS patients [7, 8]. Normally, fatigue and depression are considered manifestations of psychological disorders and interact with physical pain and discomfort, which creates a vicious cycle. Fatigue in pSS is induced and regulated by genetic and molecular mechanisms, with the innate immune system playing an important role in the production of fatigue [9–11]. Although pSS always comes with fatigue, not all patients exhibit fatigue, which provides a good model for exploring the underlying biological mechanisms.
High-throughput methods play an increasingly essential role in biology spheres, and microarray data analysis highlights its advantage in large-scale analysis of gene expression among high-throughput applications [12, 13]. Former studies [14, 15] have shown the high-throughput sequencing analysis result for pSS patients with fatigue features but do not offer further analysis based on varying degrees of fatigue. This study tries to present characteristic genes and biological pathways in pSS patients with manifestations of fatigue, as well as drugs of potential benefit.
The GSE66795 dataset from the GPL10558 platform on the GEO database is selected for gene expression of pSS with fatigue. The GSE66795 dataset was first identified for differentially expressed genes (DEGs) in pSS patients with different levels of fatigue, and based on the coexpressed genes, further analyses including Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway are performed to understand the biological process. The top ten target genes from the protein-protein interaction (PPI) network will be obtained to identify potential drugs that may alleviate fatigue in pSS patients.
2. Materials and Methods
2.1. Dataset Collection
We search “Primary Sjögren’s syndrome” and “fatigue” in the GEO database [16] and select the dataset (GSE66795) demonstrating gene expression in pSS patients with varying degrees of fatigue characteristics and normal controls. The GSE66795 dataset is extracted from the GPL10558 platform (Illumina HumanHT-12 V4.0 expression microbead chip) for RNA sequence analysis. The data of GSE66795 is obtained from the UK registry of primary Sjögren’s syndrome. It includes whole genome microarray profiles of pSS patients with varying degrees of fatigue characteristics and normal controls in peripheral blood. One hundred and thirty-one patients with pSS are involved, including 21 patients with mild fatigue, 74 patients with moderate fatigue, 36 patients with severe fatigue, and 29 normal controls.
2.2. Differential Expression Analysis
Differential expression analysis is performed using the online analysis tool GEO2R; gene expression profiles of pSS patients with mild, moderate, and severe fatigue were compared with normal controls separately to identify DEGs. values and adjusted values are calculated using -tests. Genes with the following criteria were retained for each sample: (1) log2-fold change (log2FC) absolute value greater than 1 and (2) adjusted value less than 0.05. After identifying DEGs in pSS patients with varying degrees of fatigue, the online website (https://www.xiantao.love/gds) is used to plot a Venn diagram.
2.3. Gene Ontology and Pathway Discovery in Gene Set Enrichment Analysis
Gene set enrichment analysis is used to understand the general biological function and the chromosomal location of a gene [17]. For gene product annotation, the terms of Gene Ontology (GO) are used, including biological process (BP), molecular function (MF), and cellular component (CC) [18]. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways are commonly used to describe metabolic pathways [19]. GO terms and KEGG pathways were gotten through the platform Enrichr (https://amp.pharm.mssm.edu/Enrichr/) based on the DEGs [20].
2.4. Protein-Protein Interaction (PPI) Network
The information generated from the PPI network improves the understanding of protein function [21]. PPI networks are made by STRING (https://string-db.org/) after inputting the common DEGs. We analyze PPIs through Cytoscape (https://cytoscape.org/) to further present the network and identify target genes.
2.5. Transcription Factor- (TF-) Gene Interactions
We use NetworkAnalyst (https://www.networkanalyst.ca/) to identify interactions of TF-genes with DEGs [22]. NetworkAnalyst plays a comprehensive network platform for gene expression across a wide range of species and enables them to be subjected to a meta-analysis [23].
2.6. Identification of Potential Treatment Targets
Identification of drug molecules is a vital component of genomics research. We input the DEGs in the Drug Signature Database (DSigDB). Then, we get the designed drug molecules, which may have promising clinical application. DSigDB is obtained through the Enrichr (https://amp.pharm.mssm.edu/Enrichr/) platform. Enrichr is primarily used as an enrichment analysis platform, providing extensive visual details of the common functions of inputted genes [24].
3. Results
3.1. DEG Identifications
We use the GSE66795 dataset to identify the DEGs of pSS with fatigue. 37, 29, and 33 DEGs are obtained for pSS with mild, moderate, and severe fatigue, respectively. The collected DEGs are further compared by using the online website (https://www.xiantao.love/gds) for gathering common genes in pSS with varying degrees of fatigue. And 27 (OAS1, OAS2, GBP1, IRF7, EIF2AK2, IFIT2, USP18, SAMD9L, HES4, IFI44L, SERPING1, IFIT3, IFITM3, IFI6, XAF1, MX1, OASL, OTOF, HERC5, LY6E, EPSTI1, OAS3, ISG15, IFIT1, RSAD2, IFI44, and IFI27) common DEGs are identified. The specific genes to pSS with mild fatigue are DDX60, IFIH1, GBP5, LAP3, and TIMM10. The specific genes to pSS with moderate fatigue are HLA-DRB4 and HLA-DRB6, and that to pSS with severe fatigue is RTP4. The Venn diagram (Figure 1) shows that common DEGs accounted for 67.5% out of a total of 40 DEGs.

3.2. GO Terms and KEGG Pathways
We analyzed 27 common DEGs for both GO and KEGG pathways. Both of the results are taken from the top 10 GO entries. GO terms in Table 1 suggest that DEGs are mainly localized in the mitochondria, endosomes, endoplasmic reticulum, and cytoplasm. They are involved in the biological processes of interferon action on cells and cellular defense against viruses. And the molecular functions are mainly engaged in the process of RNA synthesis. KEGG pathways in Table 2 suggest that the DEGs of pSS with fatigue are involved in the signaling pathways of viruses such as hepatitis C, influenza A, measles, and EBV. Both are seen in Figures 2(a) and 2(b).

(a)

(b)
3.3. Identification of Hub Genes by PPI Networks
We put common DEGs into the STRING website, and the files generated after analysis are further entered into Cytoscape software for visual analysis. PPI networks are designed to detect hub genes for identifying drug molecules for pSS with fatigue. PPI networks involve 24 nodes and 552 edges, which are shown in Figure 3(a). We present the top 20 genes in Figure 3(b) and Table 3.

(a)

(b)
3.4. TF-Gene Interactions
The interactions of TF and genes are shown in Figure 4. The network has 60 nodes and 108 edges. Sixteen TF-genes regulate IFIT1, and IFIT3 is handled by 14 TF-genes. The network involves 60 TF-genes. Figure 4 shows the network of TF-gene interactions.

3.5. Identification of Drug Candidates
We identify drug molecules for the top 10 hub genes on the Enrichr platform. We collect drug candidates judged on adjusted values. The analysis reveals that acetohexamide PC3 UP, suloctidil HL60 UP, prenylamine HL 60 UP, and chlorophyllin CTD 00000324 are the four most polygenic drug molecules that interact with genes. Figure 5 and Table 4 present the drug candidates in DSigDB.

4. Discussion
Fatigue is an annoying experience that means physical and mental tiredness [25]. Mengshoel et al. [26] reveal that most pSS patients literally suffer from fluctuating fatigue out of control regardless of their health condition. Fatigue has a significant influence on patients’ daily life, and patients must adapt to their behavior and lives. Although the underlying mechanisms are still unclear, former studies take depression and pain as the prominent factors associated with fatigue [5, 27]. Currently, growing evidence suggests that fatigue has a molecular and genetic basis on its production and regulation. Therefore, most scholars view fatigue as a biological and brain phenomenon [9–11].
IL-1β tends to increase rapidly secreted from macrophages to activate the immune system when meeting tissue injury or infection. IL-1β plays its role by binding with the IL-1 receptor coming with the downstream of IL-1 response [28]. Then, immune and inflammation systems are activated, which induce the behavior of disease, with fatigue being involved as an important component [29]. All these inflammatory signaling pathways go on working and turn fatigue into a chronic state. In the brain, IL-1 β signaling pathways may explain the ultimate pathway of fatigue [30, 31], and IL-1 blocker treatment may effectively release fatigue [32, 33]. Thus, fatigue and other unpleasant mood in those patients with autoimmune disease not only should be understood by the unfortunate development of chronic illness but also may be related to some signaling pathways and activation of genes that regulate the mood in the cerebral system.
Genome-wide association analysis of pSS patients has been conducted, and a gene (RTP4) is identified as highly relevant. Similarly, we confirm that RTP4 is highly expressed in pSS patients with severe fatigue through bioinformatic analysis, suggesting that this gene is critical in the mechanisms of fatigue. RTP4 encodes a protein associated with the expression of opioid receptors on the cell surface. These receptors are also expressed in the lymph system and pain-regulated pathways in the brain [34]. However, the former study did not stratify pSS based on the degree of fatigue, and it is unclear which degree of fatigue expresses the RTP4 gene. Our study finds that pSS patients with severe fatigue specifically express the RTP4 gene, providing clues for further studies on the genomics of fatigue features in pSS patients.
OAS1, a coexpressed gene for pSS in our study, has been established in previous studies as a risk locus of pSS and impacts the flaw of virus clearance because of the altering response of IFN [35]. Our gene pathway analysis points out that DEGs for pSS with fatigue are mainly localized intracellularly and involved in signaling pathways of common viruses in the respiratory and digestive tracts, suggesting that pSS is a systemic disease with an uncertain etiology and that viral infection may be a predisposing factor.
Fatigue always accompanies pSS patients, but it is hard work to manage these bad feelings [36]. The clinical practice guidelines (CPG) committee emphasizes the many causes of fatigue in pSS; therefore, the comprehensive evaluation for diagnosis is essential. So far, the treatment for fatigue in pSS with solid recommendation is mere taking exercise, which is also practical in other autoimmune diseases [37]. In America, hydroxychloroquine (HCQ) is the most widely used drug therapy for pSS with fatigue, but the recommendation strength is not strong enough [34]. It is not recommended to release fatigue in pSS using dehydroepiandrosterone (DHEA) [34]. Both the tumor necrosis factor inhibitor is discouraged for the treatment of fatigue in pSS [38, 39]. Our bioinformatic study reveals that besides chloroquine and testosterone drugs that help improve fatigue, chlorophyllin, the sulphonylurea hypoglycaemic drug acetylhexane, and the antiallergic drug terfenadine may have improved fatigue in pSS. However, chloroquine and testosterone are not strongly recommended as we mentioned before. Acetohexamide has been discontinued in the American market due to its significant hypoglycaemic risk. Terfenadine is not suitable for long-term use since its central depression as an antiallergy drug. And chlorophyllin appears to hold some promise for reducing fatigue in pSS.
Chlorophyll is an ingredient of the derifil drug which is available as an over-the-counter medicine [40]. And chlorophyllin, obtained by hydrolyzing chlorophyll to remove phytyl alcohol, is a water-soluble derivative. Chlorophyll has been shown to exert its anticancer properties by playing a role as an antioxidant [41], a CYP inhibitor [42], an apoptosis inducer [43], a phase II enzyme stimulator [44], and a carcinogen transport modulator [45]. Currently, COVID-19 has swept the world and may last for a long time because of its rapid mutation. Almost 5,000,000 people have died in this epidemic [46], and the reduction of lymphocytes in COVID-19 patients is considered an important risk factor for poor prognosis [47–49]. Recent studies suggest that the chlorophyll derivative sodium copper chlorophyllin (SCC) may improve survival in critically ill COVID-19 patients by increasing the total number of lymphocytes [50]. Increasing consumers choose dietary chlorophyll which is derived from SCC for diet supplements for the sake of keeping healthy [51, 52]. Dietary chlorophyll is safe and has been shown to have a higher absorption rate in the human body, which may trigger ionic compound chelation [53, 54]. Zeng et al. [55] cognize one functional food called barley grass powder which is rich in chlorophyll, and other nutrients can effectively alleviate fatigue in chronic patients. The mechanism of chlorophyll’s role in relieving fatigue in pSS patients is unclear. It may be related to the nature of the hepatic enzyme inhibitors that increase the concentrations of immunosuppressant like hydroxychloroquine, which has better control of fatigue. And the capacity of scavenging the oxygen radical as an antioxidant may somewhat improve the fatigue of body.
We have identified gene expression profiles in peripheral blood specific to pSS with fatigue characteristics. The analysis of identified DEGs and pathways in this study will deepen our understanding of the essence of fatigue in pSS. The discovery that chlorophyllin may improve fatigue symptoms provides a theoretical basis for better improving the quality of life in pSS patients. And a preprint has previously been published [56].
Data Availability
The dataset supporting the conclusions of this article is available in the UK registry of primary Sjögren’s syndrome repository and in the hyperlink (https://www.ncbi.nlm.nih.gov/geo/geo2r/?acc=GSE66795).
Ethical Approval
GEO belongs to public databases. The patients we choose involved in the database have obtained ethical approval. It is available for all users to download relevant data for free. Our study is based on open-source data, so there is no need to offer ethics approval.
Consent
There is no need for consent to participate.
Conflicts of Interest
The authors declare that they have no competing interests.
Authors’ Contributions
Guangshu Chen and Li Che contributed equally to this work.
Acknowledgments
We sincerely acknowledge the GEO database for providing their platforms and contributors for uploading their meaningful datasets. This study was financially supported by the Guangdong Science and Technology Project Fund for Key Scientific Research Base under Grant no. 2019B020230001.