key: cord-0817728-rqvdqi0f authors: Chen, Guangshu; Che, Li; Cai, Xingdong; Zhu, Ping; Ran, Jianmin title: Bioinformatic Analysis Identifies Biomarkers and Treatment Targets in Primary Sjögren's Syndrome Patients with Fatigue date: 2022-01-15 journal: Biomed Res Int DOI: 10.1155/2022/7697558 sha: 3f84ef110f3b3dbeb4036c648d227157cc8d2fe4 doc_id: 817728 cord_uid: rqvdqi0f 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. 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] [5] [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 produc-tion of fatigue [9] [10] [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.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. 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. P values and adjusted P values are calculated using t-tests. Genes with the following criteria were retained for each sample: (1) log2-fold change (log2FC) absolute value greater than 1 and (2) adjusted P 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. 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. 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] . 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. 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). 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 2 BioMed Research International 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 . 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 P 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. 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] [10] [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 painregulated 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 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] [48] [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 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 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] . 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). 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. There is no need for consent to participate. The authors declare that they have no competing interests. Primary Sjogren syndrome: clinical and immunopathologic features Epidemiological studies in incidence, prevalence, mortality, and comorbidity of the rheumatic diseases A new model for classification of disease manifestations in primary Sjögren's syndrome: evaluation in a retrospective long-term study A five-year prospective study of fatigue in primary Sjögren's syndrome Fatigue in primary Sjögren's syndrome: clinical, laboratory, psychometric, and biologic associations The prevalence of severe fatigue in rheumatic diseases: an international study Fatigue in the general population: a translation and test of the psychometric properties of the Norwegian version of the fatigue severity scale Fatigue in primary Sjögren's syndrome is associated with lower levels of proinflammatory cytokines Epigenome-wide DNA methylation patterns associated with fatigue in primary Sjögren's syndrome Heat shock proteins and chronic fatigue in primary Sjögren's syndrome The neuroimmune basis of fatigue Genesis: cluster analysis of microarray data Importance of replication in microarray gene expression studies: statistical methods and evidence from repetitive cDNA hybridizations Identification of differentially expressed genes in primary Sjögren's syndrome Meta-analysis of differentially expressed genes in primary Sjogren's syndrome by using microarray The Gene Expression Omnibus database Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles GoPubMed: exploring PubMed with the Gene Ontology KEGG: Kyoto Encyclopedia of Genes and Genomes Enrichr: a comprehensive gene set enrichment analysis web server BioMed Research International update Kernel methods for predicting protein-protein interactions Bioinformatic identification of candidate biomarkers and related transcription factors in nasopharyngeal carcinoma NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool Mechanisms and management of fatigue in progressive neurological disorders Primary Sjögren's syndrome: fatigue is an ever-present, fluctuating, and uncontrollable lack of energy Fatigue in rheumatoid arthritis reflects pain, not disease activity Treating inflammation by blocking interleukin-1 in a broad spectrum of diseases Biological basis of the behavior of sick animals Neuroinflammation drives anxiety and depression in relapsing-remitting multiple sclerosis Associations of tumor necrosis factor-α and interleukin-1β levels and polymorphisms with post-stroke depression The effect of interleukin-1 blockade on fatigue in rheumatoid arthritis-a pilot study Inhibition of IL-1beta improves fatigue in type 2 diabetes Treatment guidelines for rheumatologic manifestations of Sjögren's syndrome: use of biologic agents, management of fatigue, and inflammatory musculoskeletal pain Identification of a Sjögren's syndrome susceptibility locus at OAS1 that influences isoform switching, protein expression, and responsiveness to type I interferons Fatigue in primary Sjogren's syndrome Effects of exercise on aerobic capacity and fatigue in women with primary Sjogren's syndrome Etanercept in Sjögren's syndrome: a twelve-week randomized, double-blind, placebo-controlled pilot clinical trial Inefficacy of infliximab in primary Sjögren's syndrome: results of the randomized, controlled Trial of Remicade in Primary Sjögren's Syndrome (TRIPSS) Amelioration of radiation-induced hematopoietic syndrome by an antioxidant chlorophyllin through increased stem cell activity and modulation of hematopoiesis Chlorophyllin as a protector of mitochondrial membranes against γ-radiation and photosensitization Non-specific inhibition of cytochrome P450 activities by chlorophyllin in human and rat liver microsomes The chlorophyllininduced cell cycle arrest and apoptosis in human breast cancer MCF-7 cells is associated with ERK deactivation and cyclin D1 depletion Chlorophyll, chlorophyllin and related tetrapyrroles are significant inducers of mammalian phase 2 cytoprotective genes Effects of chlorophyllin on transport of dibenzo(a, l)pyrene, 2-amino-1-methyl-6-phenylimidazo-[4,5-b ]pyridine, and aflatoxin B 1 across Caco-2 cell monolayers COVID-19 Dashboard Hematological findings and complications of COVID-19 Lymphopenia is associated with severe coronavirus disease 2019 (COVID-19) infections: a systemic review and meta-analysis Functional exhaustion of antiviral lymphocytes in COVID-19 patients COVID-19 therapy: could a copper derivative of chlorophyll a be used to treat lymphopenia associated with severe symptoms of SARS-CoV-2 infection? An evidencebased systematic review of chlorophyll by the Natural 9 Phycobilins and phycobiliproteins used in food industry and medicine Sodium copper chlorophyllin: in vitro digestive stability and accumulation by Caco-2 human intestinal cells Galacto-oligosaccharides may directly enhance intestinal barrier function through the modulation of goblet cells Preventive and therapeutic role of functional ingredients of barley grass for chronic diseases in human beings Bioinformatic analysis identifies biomarkers and treatment targets in primary Sjögren's syndrome patients with fatigue 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. Guangshu Chen and Li Che contributed equally to this work.