key: cord-0050554-4k7bh40t authors: Khayer, Nasibeh; Mirzaie, Mehdi; Marashi, Sayed-Amir; Jalessi, Maryam title: Rps27a might act as a controller of microglia activation in triggering neurodegenerative diseases date: 2020-09-17 journal: PLoS One DOI: 10.1371/journal.pone.0239219 sha: a1e87159c44bf00b33e393ffe11af38306a40697 doc_id: 50554 cord_uid: 4k7bh40t Neurodegenerative diseases (NDDs) are increasing serious menaces to human health in the recent years. Despite exhibiting different clinical phenotypes and selective neuronal loss, there are certain common features in these disorders, suggesting the presence of commonly dysregulated pathways. Identifying causal genes and dysregulated pathways can be helpful in providing effective treatment in these diseases. Interestingly, in spite of the considerable researches on NDDs, to the best of our knowledge, no dysregulated genes and/or pathways were reported in common across all the major NDDs so far. In this study, for the first time, we have applied the three-way interaction model, as an approach to unravel sophisticated gene interactions, to trace switch genes and significant pathways that are involved in six major NDDs. Subsequently, a gene regulatory network was constructed to investigate the regulatory communication of statistically significant triplets. Finally, KEGG pathway enrichment analysis was applied to find possible common pathways. Because of the central role of neuroinflammation and immune system responses in both pathogenic and protective mechanisms in the NDDs, we focused on immune genes in this study. Our results suggest that "cytokine-cytokine receptor interaction" pathway is enriched in all of the studied NDDs, while "osteoclast differentiation" and "natural killer cell mediated cytotoxicity" pathways are enriched in five of the NDDs each. The results of this study indicate that three pathways that include "osteoclast differentiation", "natural killer cell mediated cytotoxicity" and "cytokine-cytokine receptor interaction" are common in five, five and six NDDs, respectively. Additionally, our analysis showed that Rps27a as a switch gene, together with the gene pair {Il-18, Cx3cl1} form a statistically significant and biologically relevant triplet in the major NDDs. More specifically, we suggested that Cx3cl1 might act as a potential upstream regulator of Il-18 in microglia activation, and in turn, might be controlled with Rps27a in triggering NDDs. Neurodegenerative diseases (NDDs) are increasing serious menaces to human health in recent years [1] . In spite of the considerable body of research on NDDs such as Alzheimer's disease (AD), amyotrophic lateral sclerosis (ALS), Huntington's disease (HD), multiple sclerosis (MS), Parkinson's disease (PD) and schizophrenia (SCHIZ), effective treatment for preventing the progressive neuronal loss is yet to be available [2] . Despite exhibiting different clinical phenotypes and selective neuronal loss in neurodegenerative abnormalities, there are many common pathological features that contribute to these disorders, which suggest the existence of certain common pathways at the molecular level. A few of such common pathological features are oxidative stress [3] , mitochondrial dysfunction [4] , glutamate excitotoxicity [5] , abnormal protein aggregation [6] , abnormal cellular transport [7] , iron accumulation [8] and neuroinflammation [9, 10] . Therefore, it would be interesting to survey the NDDs at once to find such a potential similarity at the molecular level. The number of comparative studies in NDDs is still rather low. Therefore, our knowledge is limited to the underling common pathways and central genes in these disorders. In a few previous studies, the major NDDs were comparatively surveyed at the systemic level. Durrenberger and coworker generated comparative genome-wide gene expression data for six main NDDs, namely AD, ALS, HD, MS, PD, and SCHIZ. They performed whole-genome expression analysis to look for potential common molecular pathogenic mechanisms and genes, in addition to exclusive disease-specific changes. They reported that no dysregulated gene was found in common across all six NDDs. Furthermore, they reported that "neuronal homeostatic", "survival activity" and "synaptic plasticity" pathways are shared among five and four of the diseases, respectively [11] . On the other hand, Godini and coworkers performed a comparative system-level analysis on NDDs microarray data that was generated by Durrenberger and coworkers. They constructed a protein-protein interaction network based on the co-expression approach and performed gene set enrichment analysis to find out the central factors that are related to each disease and, in addition, those features that can be considered as the similarities among them. Finally, they suggested a set of key genes that might contribute to all NDDs. Nevertheless, in this study, similar to the above mentioned, no pathway or gene was found to be in common across all six major NDDs [12] . In the present work, we applied the three-way interaction model to find possibly shared genes and pathways among the major NDDs (Fig 1) . High-throughput gene expression microarray data is a form of genome-scale data that provide an opportunity to trace gene expression patterns. There are numerous gene expression patterns which contain diverse information about gene relations [13, 14] . One of them is cross-shaped co-expression pattern. Such pattern is observed when the expression levels of two genes are either directly or inversely correlated depending on the expression level (or genotype) of a third gene (i.e., the switch gene). Three-way interaction mechanism can be inferred from a cross-shaped coexpression pattern. In this model, the changing nature of co-expression relationship of two genes{X 1 , X 2 } is captured by introducing a third gene (X 3 ) which modulates the X 1 -X 2 expression dependence. Therefore, such model can potentially shed light on some sophisticated gene interactions that cannot be traced trace with the simplistic correlative model. In addition, three-way interaction model can provide valuable information about gene regulatory relations [13, 15] . A significant body of evidence shows that neuroinflammation and immune system responses have both pathogenic and protective roles in the neurological disorders [9, [16] [17] [18] . Microglia and other innate immune cells such as astrocytes with a morphologically activated phenotype are observed at relatively early stages in large numbers of NDDs. Such cells can be activated in response to misfolded proteins, aberrantly localized nucleic acids or reactive oxygen species (ROS) [19] . Activation of innate immune cells leads to sever inflammatory responses in the central nervous system (CNS), renewed promoting neural damage. More specifically, Toll-like receptors (TLRs) and other pattern recognition receptors, which are expressed on microglia, play central roles in triggering neuroinflammation. Neural damage and subsequent inflammatory responses result in a feed-forward loop of neurodegenaration [19] . In the present study, we focus on the immune genes to look for NDDs related mechanisms and central switch genes that are involved in the immune responses using the three-way interaction model. In the present work, we exploited the gene expression data reported by Durrenberger et al. [11] . The selected dataset includes comparative genome-wide gene expression data for AD, ALS, HD, MS, PD and SCHIZ from 118 cases of postmortem brain tissues, belonging to entorhinal cortex, cervical spinal cord, ventral head of the caudate nucleus, substantia nigra, grey matter in brodmann area, respectively. The data were generated using the Illumina human-Ref-8 v2.0 expression bead chip platform and are available in GEO dataset under accession number GSE26927 [11] . A brief description of the dataset was presented in S1 Table in S1 Material. The duplicated probes were removed using the genefilter package in Bioconductor [20] . For this purpose, the mean of expression levels of duplicated probes were retained for each gene. Finally, the probes with "immune" annotations were selected based on the Biological Processes of the Gene Ontology annotation for the rest analysis [21] . Additionally, in order to identify differentially expressed genes (DEGs), for each disease, theempirical Bayes t-test was performed using the LmFit function in the Limma package [22] for each disease separately. Three-way interactions among all passible triplet combinations of genes in the dataset were calculated using fastMLA function in fast Liquid Association package [23] . This package uses a modified liquid association algorithm for determining changes in co-expression relations of two genes, X 1 and X 2 , based on the expression level of a third gene (X 3 ). Before running liquid association algorithm, performing two preprocessing steps on the data is necessary: (i) normal quntile transformation [24]; (ii) standardization of mean and standard deviation to o and 1, respectively [25] for each variable. The first preprocess was performedby an in-house implementation and the second one by using CTT package [26] . False discovery rate (FDR) was estimated by using the Bonferroni correction method [27] . Subsequently, the three-way interactions with FDR < 0.01 were reported as the statistically significant triplets. A gene regulatory network (GRN) model provides a systematic knowledge about the complex molecular mechanisms that control gene expression under diverse biological conditions [28] . Here, we used geWorkbench_2.6.0 software to construct GRN among all genes in dataset as hub markers by considering p-value <0.05. Specifically, ARACNE (Algorithm for the Reconstruction of Accurate Cellular Networks) was applied [29] to construct the GRN. This algorithm uses mutual information as a measure to detect direct regulatory interactions between each transcriptional regulator and its potential target(s). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis is a popular statistical method to determine the validation of the shared association of a set of genes using predefined pathways [30] . Here, KEGG pathway enrichment analysis was performed separately for each of the six sets of DEGs using the KEGG database [31] . Above-mentioned analyses were performed using ClueGO plugin [32] within the Cytoscape v.3.3.0 environment [33] , and in addition, two-sided hypergeometric test with the Bonferroni step down correction method and a kappa score threshold of 0.4 were considered for statistical significance calculation. The presented study is a computational research using a public GEO dataset (accession number: GSE26927). Therefore, any ethics approval and consent statements were not required for this study. After preprocessing and extracting the immune-related genes from dataset, 822 genes remained for further analysis. In order to determine differentially-expressed genes, eBayes ttest was performed on the normalized data. By considering p-value <0.05, we selected 21, 150, 171, 111, 145 and 83 genes as significant differentially-expressed genes in AD, ALS, HD, MS, PD and SCHIZ, respectively. The detailed results of this analysis are available in S2 Table in S1 Material. All 822 immune-related genes were considered as potential switch genes in the fast liquid association analysis. The changes in FDR versus -log(p-value) for the first 300000 triplets were showed in S3 Fig in S1 Material. In order to survey the accuracy of fastLA analysis, observed event rate of X 3 position (switch) genes in the wide range of the significant fastLA p-values were compared with random event rate (see Fig 2) . On the other hand, the relevancy of X 3 position (switch) genes with DEGs was examined. As shown in Fig 3, the percentage of the DEGs in all of the immune-related genes and the percentage of triplets that X 3 position of them is DEG in all statistically significant triplets is 25% and 53%, respectively. With considering FDR <10 −2 , a set of significant cross-shaped triplets which consists of 1717 triplet combinations were selected for further analyses. In the wide range of the significant fastLA p-values, observed event rate of X 3 position (switch) genes were compared with random event rate. As shown, the observed event rate for switch genes is far greater than the random one, which confirms the validity of fastLA analysis. https://doi.org/10.1371/journal.pone.0239219.g002 We used GSEA in order to find common pathway(s) in NDDs in addition to detect the biologically-relevant triplets. By considering FDR <0.05 and a minimum of 3 genes in each groups, three pathways were enriched simultaneously at least in 5 of the NDDs. All of the significantly enriched terms based on "KEGG pathway" were presented in Fig 4A, and in addition, a detailed information about KEGG pathway enrichment analysis was reported in S4 Table in S1 Material. Based on the results, "cytokine-cytokine receptor interaction" was enriched in all 6 NDDs, while "osteoclast differentiation", "natural killer cell mediated cytotoxicity" were enriched in 5 of the NDDs each ( Fig 4B) . Based on the definition of three-way interactions in switch mechanism model, it is expected that in biologically-relevant triplets, X 1 and X 2 are in The proportion of triplets in which the X 3 position is a DEG. As shown, the percentage of DEGs in A and B is 25% and 53%, respectively. The prevalence of DEGs in the X 3 position of triplets suggest that the liquid association model is appropriate for understanding the role of DEGs in neuroinflammation. https://doi.org/10.1371/journal.pone.0239219.g003 the same pathway. By detecting the statistically significant triplets in the enriched terms, 89 triplets were found to be biologically relevant based on GSEA. These triplets are reported in S5 Table in S1 Material. As another attempt for analyzing the functional relevance of statistically significant triplets, we reconstructed a GRN by ARACNE. By considering p-value <0.05, the GRN included 3199 edges and 822 nodes. Consequently, the statistically significant triplets were mapped on to the GRN. The results showed that the regulatory relations of 242 triplets are detectable in this GRN. Overall, the biological relevance of the 14 statistically significant triplets was defined by both GRN and GSEA methods. The position of such triplets are shown as a subnetwork of GRN in Fig 5. In addition, the liquid association analysis information of these triplets were reported in Table 1 . One of these triplets is the 1133 nd triplet, which consists of genes Rps27a and {Il18, Cx3cl1}. Although some triplets have a higher liquid association score than 1133 nd triplet, there is a feature that distinguishes such triplets from the other ones. 1133 nd triplet is the only statistically significant triplet that is not only biologically relevant based on both GRN and GSEA analyses but also it is involved in 5 diseases (i.e., across all six NDDs except ALS). On the other hand, the biological relevance of 1328 nd triplet, which is involved in common across all 6 NDDs, was confirmed only by GSEA method. This triplets consists of genes Bcl2 and {Ltbr, Il18}. The scatter plots of 1133 nd and 1328 nd triplets in three different ranges of expression level of their associated X 3 are presented in Fig 6. As shown, a considerable change in the correlation of X 1 and X 2 occurs as a result of the changes in X 3 expression level. Although different NDDs are clinically characterized with certain distinctive features, these diseases exhibit share common neuropathological characteristics [34] . One such characteristic is neuroinflammation which is associated with the immune system processes such as inflammation and immune response. Previously, a number of studies investigated the molecular factor(s) of such common changes [11, 12] . However, to the best of our knowledge, no pathway or gene was reported to be in common across all major NDDs. In the present study, the three-way interaction model was used for to investigate the immune-related switch mechanisms in NDDs. Our results suggests that a number of immune-related switch genes act variety of different NDDs. The liquid association analysis information of the 14 triplets that are statistically, as well as biologically, significant. The Pearson's correlation coefficient of X 1 and X 2 when X 3 is at a high gene expression level is indicated by ρ high (i.e., rho{X 1 , X 2 }|X 3 ") and when X 3 is at a low gene expression level is indicated by ρ low (i.e., rho{X 1 , X 2 }| X 3 "). In addition, ρ diff is the absolute value of their difference (i.e., ρ diff = |ρ high-ρ low |). Arrows indicate direct ("") or inverse ("#) correlation of X 1 and X 2 in the specific state of X 3 . https://doi.org/10.1371/journal.pone.0239219.t001 As the mentioned above, the purpose of the present study was to perform a comparative analysis on the involved mechanisms in NDDs. Because distinct regions of brain are affected in different NDDs [35] , there was an inevitable diversity among the brain regions from which the samples are taken. On the other hand, several studies show that the properties of cell population in samples depended on the region of brain [36] . Therefore, in order to reduce the variability between samples, a preprocessing and normalization step as the first stage was performed. More specifically, the quantile normalization method [37] was utilized to minimize the batch effects. The results of eBayes t-test analysis suggest that the number of DEGs are strikingly different among these diseases. The lowest and the highest number of immune-related DEGs belong to AD and HD, respectively. This observation is consistent with the idea that beside a number of morphological similarities, fundamental differences exist between the natures of these diseases. As shown in Fig 2, the observed event rate of X 3 position (switch) genes in a wide range of significant fastLA p-values is far greater than what is expected by chance. In other words, most X 3 positions in the statistically significant triplets are occupied by certain genes. Therefore such result confirms the validity of fastLA analysis. In addition, in Fig 3, we show that such certain X 3 positions genes are frequently belong to DEGs. As shown in Fig 3, the percentage of the DEGs in the all of the immune genes is 25%, while the percentage of DEGs in the all of the X 3 position genes is 53%. This results shows that DEGs play a central role in the switch mechanisms. The results of the KEGG enrichment analysis (Fig 4) show that three pathways are commonly involved in at least five of the NDDs. These pathways include "osteoclast differentiation", "natural killer cell mediated cytotoxicity" and "cytokine-cytokine receptor interaction". In the following, the relationship between these pathways and NDDs is explained in detail. From its name, the osteoclast differentiation pathway seems to be an unrelated pathway to the NDDs at the first sight. It is interesting to note that the link between osteoclast differentiation pathway and the NDDs is suggested in numerous studies [38] [39] [40] [41] [42] [43] Another enriched biological pathway for NDDs is "Natural killer cell mediated cytotoxicity" (Fig 4) . Natural killer (NK) cells are large granular lymphocytes that have cytotoxic properties [56] . The cytotoxicity of NK cells is carried out a collectively by inhibitory and stimulatory receptors that are expressed on NK cells surface [57] . In addition, NK cells have a crucial role in innate and adaptive immune system [58] . The central role of the NK cell activity in AD, MS, PD and SCHIZ has been reported (see below). Some evidence suggests that the NK cell activity in AD patients is significantly lower than that in normal controls [59] . Moreover, the cytotoxic activity of NK cells is associated with neurological inflammation observed in neuronal degeneration [60] . In addition, the NK cells activity has been suggested as a biomarker for AD progression [61] . Likewise, the role of NK cell in MS disease has been well studied [56, 62] . Uchida and colleagues showed that NK cell activity in MS patients is significantly lower than that in normal ones. They also reported that the interferon treatment of lymphocytes leads to an increase in the activity of NK cells that is less than in MS patients compared to normal ones [63] . Chanvillard and colleagues showed that mitoxantrone, a cytotoxic drug that is approved for treatment of progressive MS, induces maturation of NK cells in patients with progressive MS [64] . Rodriguez-Martin and colleagues reported that all NK cell subsets significantly increased in MS and other inflammatory neurological diseases compared to in non-inflammatory neuropathies [65] . In addition, Segal reported that in MS, NK cells curb neuroinflammation by inhibiting the auto-reactive T cell responses. In other words, NK cells by regulating of some adaptive immune responses, restricts neuroinflammation in MS [66] . Furthermore, the role of NK cells in PD has been also studied. Mihara and colleagues reported that expression of inhibitory receptor CD94/NKG2A was remarkably lower in PD patients than those in non-PD individuals. Actually, inhibitory receptor CD94/NKG2A suppresses activation signaling processes, which trigger after binding NK cells with its ligand on the normal cells, to avoid destruction of normal bystander cells [67] . Therefore, because of limited number of such inhibitory receptor in PD patient and subsequently a deficiency in inhibition of NK cells activation in the normal cells, the PD patients are prone to disruption of normal cells. [68] . Jiang and colleagues conducted a systematic review to found relevance between lymphocyte subsets/NK cell and the risk of PD. Their results showed that the number of NK cells significantly increased in PD; and consequently, such cells may be associated with the risk of PD [69] . Additionally, there are numerous studies on the activation of NK cells in SCHIZ. Schindler and colleagues compared the activity of NK cells in the blood of SCHIZ patients and healthy controls, but they did not find significant differences between the two groups [70] . On the other hand, Abdeljaber and colleagues reported that the activity of NK cells in SCHIZ patients decline compared to the control group [71] , although similar results are not observed in other studies. Another biological pathway that is enriched in the six major NDDs is "cytokine-cytokine receptor interaction" (Fig 4) . Cytokines are produced by peripheral immunocompetent, glial cells and neurons [72] . During immune responses, cytokines mediate specific signaling communications between neurons and immune cells. Additionally, cytokines play pleiotropic roles in the central nervous system, which include roles in synaptic plasticity, neurotransmission, nuclear signal transduction, neurogenesis, and inflammatory responses [73, 74] . While the pathophysiological role of cytokines in NDDs is well known [75] , to our knowledge, the cytokine-cytokine receptor interaction (CCRI) pathway was reported only in MS, SCHIZ, and AD in previous studies. See below. Xu and colleagues have previously indicated that the change in of Il-18 pathway may contribute to psychopathology of SCHIZ. They also showed that the Il-18 pathway was strongly associated with the CCRI pathway [74] . The relevance of the pathway to MS has been also reported. Guo and colleagues, by comparison between blood transcriptome of MS and healthy individuals, showed that genes involved in the CCRI pathway were potentially associated with MS [76] . Likewise, about AD, Sattlecker and colleagues reported that the members of CCRI pathway change during AD in proteome-level quantification [77] . On the other hand, Khayer and colleagues reported that CCRI pathway is a critical switch pathway in triggering AD [15] . The 1133 nd triplet which contains of genes Rps27a and {Il-18, Cx3cl1} is one of the most statistically significant triplets in our analysis. Furthermore, this triplet is found to be biologically relevant, based on GSEA and GRN analyses ( Table 1 and Fig 5) . As shown in Fig 6, when normalized expression level of Rps27a gene is between 0.40 and 2.39, there is an inverse correlation between Il-18 and Cx3cl1 expression levels. In contrast, when normalized expression level of Rps27a gene is between -2.39 and -0.40, there is a direct correlation between Il-18 and Cx3cl1 expression levels. In other words, changes in expression level of Rps27a gene can act as the switch factor for altering the regulated genes. To the best of our knowledge, no previous study has investigated the role of Rps27a gene in NDDs. Nevertheless, there is considerable evidence to support that Rps27a gene may be associated with this category of diseases. Rps27a gene encodes a fusion protein consisting of ubiquitin (Ub) at the N terminus and ribosomal protein S27a at the C terminus [78] . Rps27a protein is cleaved to free ubiquitin monomer and ribosomal protein S27a by a cysteine protease, namely UCHL1 [79] . In other words, UCHL1 is a deubiquitinating enzymes (DUBs). Ub is a versatile regulatory protein that plays diverse roles in cells [80] . In addition, the cellular Ub pool, which is composed of free Ub and Ub conjugates, are in dynamic equilibrium inside cells. On the other hand, the maintenance of cellular levels of Ub (i.e., Ub homeostasis) is substantial for appropriate ubiquitination of substrate proteins [81] . Previous studies indicated that disruption of Ub homeostasis plays a critical role in the pathophysiology of various NDDs [82] [83] [84] [85] [86] . There are at least two ways for maintaining Ub homeostasis in cells. One of them is regulation of DUBs that convert Ub conjugates into their monomeric free forms, and consequently increase the levels of Ub pool [87] . UCHL1, a DUB enzyme for Rps27a (see previous paragraph), dysfunction is associated with a variety of NDDs, including AD [88, 89] , PD [90] [91] [92] , HD [93, 94] and SD [95] . The other way is regulation of gene expression of monomeric Ub-ribosomal fusion genes as well as polyubiquitin genes. As mentioned above, Rps27a is a monomeric Ub-ribosomal fusion gene; therefore regulation of such gene can be a way to maintain cellular Ub pool and in turn, might contribute to NDDs. As it was mentioned above, Il-18 and Cx3cl1 genes that are involved in triplet 1133 nd are directly correlated in some conditions. Surprisingly, the positive relationship between Il-18 and Cx3cl1 was reported in previous studies. Alboni and colleagues demonstrated that fluoxetine, as a chronic antidepressant drug, that is involved in anti-inflammatory, anti-apoptotic and antioxidant activity [96] , decreases hypothalamic expression of Il-8 and Cx3cl1 genes simultaneously [97] . On the other hand, both Il-18 and Cx3cl1 genes were suggested as interferon-related biomarkers for Listeria monocytogenes by Koopmans and colleagues [98] . In addition, Kasama and colleagues reported that gene expression levels of Cx3cl1 and Il-18 are positively correlated with each other in adult-onset Still's disease [99] . Furthermore, Bian and collegues indicated that the interaction between Cx3cl1 and Il-18 signaling plays a critical role in the development of allodynia [100] . Some evidence suggests that Cx3cl1 might be as an upstream regulator of Il-18 in microglia cells. See below. Cx3cl1 and Il-18 genes are involved in spinal long term potentiation [101] as well as microglia activation [102] . Bian and colleagues, by utilizing double immunostaining technique, showed that Cx3cl1 receptor (Cx3cr1) is co-localized with Il-18 in the spinal cord and predominately expressed in spinal microglia in rats. They also reported that blockade of Cx3cl1/ Cx3cr1 signaling and/or Il-18, obviously suppresses the spinal long term potentiation [101] . On the other hand, Miyoshi and colleagues determined that production of Il-18 in the spinal cord was regulated by p38MAPK [103] . Furthermore, Zhang and colleagues reported that the p38MAPK signaling was activated in spinal microglia by exposing to exogenous Cx3cl1 [104] . Therefore, it is reasonable to infer that Cx3cl1 might act as a potential upstream regulator of Il-18 in microglia cells. Besides, some previous studies indicate that microglia activation is involved in spinal long term potentiation [105, 106] as well as in NDDs [107] . It should be mentioned that the central role of Il-18 gene in AD [108] , ALS [109, 110] , MS [111] and SD [74, 112, 113] , and Cx3cl1 gene in all of the six surveyed diseases (i.e. AD [114, 115] , ALS [116] , HD [117] , MS [118] , PD [119] [120] [121] [122] and SD [123, 124] ) has been determined. However, as far as our knowledge goes, the relation between Il-18 and Cx3cl1 in triggering of NDDs is not reported in the literature yet. Microglia, the central participator in the immune system of CNS, are versatile cells. Such cells, depending on environment signals, can acquire different phenotypes [125] . Two novel known subtypes of microglia that occur in mouse models of NDDs, are disease-associated microglia (DAM) [126] and primed microglia [127] . In this section, we asked whether induction of a DAM-like and/or primed-like expression program is traced in microglia of human NDDs. To this end, we compared our results with recently published gene expression data of DAM and primed microglia. See below. Recently, Keren-Shaul and co-authors [126] reported a protective sub-type of microglia associated with NDDs, termed disease-associated microglia (DAM), in AD and ALS mouse model. They utilized single-cell RNA sequencing technique on microglia isolated from mouse model of AD (5XFAD) and ALS (mSOD1 [G93A] mice) as well as the corresponding wildtype mice. Subsequently, using graph-based clustering analysis, they identified a DAM gene signature for mouse model of AD and ALS. We compared our results with DAM characteristic genes to investigate the DAM-like expression program in human NDDs. Of the 22 DAM characteristic genes, reported by Keren-Shaul and co-authors, 7 genes, including Cx3cr1, Tyrobp, B2m, Fth1, Trem2, Axl and Csf1, were found in our results as involved genes in statistically significant triplets. Additionally, Trem2, a key gene that enable the microglial progression to a full activation DAM profile, was found in our results as a switch gene. On the other hand, they highlighted lysosomal, phagocytic, and lipid metabolism pathways for DAM characteristic genes. In a comparative aspect, none of DAM-related pathways were enriched in our study for any disease. Taken together, there are not a striking overlap between DAM characteristic genes and key genes reported in our study. Two possible reasons for such observation may be: (i) differences in applied transcriptomic analysis techniques, which can affect the accuracy of the study; (ii) differences in the nature of investigated samples. To the best of our knowledge, a slight genetic differences between mouse and human might be translated in to a salient distinction [128] . Another subtype of microglia that is reported for mouse model of NDDs is primed microglia. Such subtype of microglia acquires a hypersensitive phenotype, becoming hyper-reactive and secreting large amount of inflammatory mediators in response to signal from a systemic inflammatory event [129] . Holtman and co-authors [127] , to determine the gene expression signature of priming, analyzed the transcriptome of microglia in mouse model of AD and ALS using weighted gene co-expression network analysis. They identified four hub genes specific for the primed microglia expression network, including Lgals3, Igf1, Csf1, and Axl. Interestingly, all of hub genes except Igf1 were found in our results. These genes are involved in X 3 position and/or X 1 , X 2 of significant triplets. Hence, it can be expected that a primed-like expression program triggers in human neurodegenerative diseases as in mouse model occurs. However, our knowledge for such claim is not sufficient and more evidence is needed to support it. Recent study is the first time that was utilized three-way interaction approach to survey common pathways in the NDDs. Three-way interaction approach is a promising strategy that provides new type of information on the specific relationships among genes. In conclusion, three common switch pathways that include "osteoclast differentiation", "natural killer cell mediated cytotoxicity" and "cytokine-cytokine receptor interaction" were found in the major NDDs. More specifically, our analysis showed that Rps27a as the switch gene and the gene pair {Il-18, Cx3cl1} form a statistically significant and biologically relevant triplet in the major NDDs. Additionally, we suggested that Cx3cl1 might act as a potential upstream regulator of Il-18 in microglia activation, and in turn, might be controlled with Rps27a in triggering NDDs. The coming crisis: obtaining care for the growing burden of neurodegenerative conditions Drug development for CNS disorders: strategies for balancing risk and reducing attrition Oxidative stress and neurodegenerative diseases: a review of upstream and downstream antioxidant therapeutic options Mitochondrial dynamics-fusion, fission, movement, and mitophagy-in neurodegenerative diseases Excitotoxicity: bridge to various triggers in neurodegenerative disorders Protein aggregation and neurodegenerative disease Axonal transport defects in neurodegenerative diseases Iron, brain ageing and neurodegenerative disorders Inflammation in neurodegenerative diseases Inflammation in neurodegenerative disease-a double-edged sword Selection of novel reference genes for use in the human central nervous system: a BrainNet Europe Study A comparative system-level analysis of the neurodegenerative diseases Three-way interaction model with switching mechanism as an effective strategy for tracing functionally-related genes Detecting differentially coexpressed genes from labeled expression data: a brief review Three-way interaction model to trace the mechanisms involved in Alzheimer's disease transgenic mice Innate and adaptive immune responses in neurodegeneration and repair Neurodegeneration in autoimmune CNS inflammation Innate and Adaptive Immune-Mediated Neuroinflammation and Neurodegeneration in Parkinson's Disease. Neuroinflammation and Neurodegeneration Innate immune activation in neurodegenerative disease Genefilter: methods for filtering genes from highthroughput experiments Natural killer cells and their receptors in multiple sclerosis Natural killer cell activity in patients with multiple sclerosis: interferon and plasmapheresis Mitoxantrone induces natural killer cell maturation in patients with secondary progressive multiple sclerosis Natural killer cell subsets in cerebrospinal fluid of patients with multiple sclerosis The role of natural killer cells in curbing neuroinflammation The cell biology of the human natural killer cell CD94/NKG2A inhibitory receptor Natural killer cells of Parkinson's disease patients are set up for activation: a possible role for innate immunity in the pathogenesis of this disease The correlation of lymphocyte subsets, natural killer cell, and Parkinson's disease: a meta-analysis Studies of cellular immunity, serum interferon titers, and natural killer cell activity in schizophrenic patients Depressed natural killer cell activity in schizophrenic patients Cytokine production in the central nervous system Cytokine-to-brain communication: a review & analysis of alternative mechanisms Potential involvement of the interleukin-18 pathway in schizophrenia Cytokines and acute neurodegeneration Mining gene expression data of multiple sclerosis Longitudinal Protein Changes in Blood Plasma Associated with the Rate of Cognitive Decline in Alzheimer's Disease Database Resources of the National Center for Biotechnology Information The functions of UCH-L1 and its relation to neurodegenerative diseases Regulatory mechanisms involved in the control of ubiquitin homeostasis Cellular ubiquitin pool dynamics and homeostasis Walking the tightrope: proteostasis and neurodegenerative disease The ubiquitin-proteasome system in neurodegenerative diseases: precipitating factor, yet part of the solution SOD1(A4V) aggregation alters ubiquitin homeostasis in a cell model of ALS Ubiquitin homeostasis: from neural stem cell differentiation to neuronal development UBA1: At the Crossroads of Ubiquitin Homeostasis and Neurodegeneration Cellular ubiquitin pool dynamics and homeostasis It Is All about (U)biquitin: Role of Altered Ubiquitin-Proteasome System and UCHL1 in Alzheimer Disease MicroRNA-922 promotes tau phosphorylation by downregulating ubiquitin carboxy-terminal hydrolase L1 (UCHL1) expression in the pathogenesis of Alzheimer's disease Association between the ubiquitin carboxyl-terminal esterase L1 gene (UCHL1) S18Y variant and Parkinson's Disease: a HuGE review and meta-analysis Ubiquitination of alpha-synuclein and autophagy in Parkinson's disease The effect of Parkinson's-disease-associated mutations on the deubiquitinating enzyme UCH-L1 Polymorphism of HD and UCHL-1 genes in Huntington's disease The S18Y polymorphism in the UCHL1 gene is a genetic modifier in Huntington's disease Alterations of ubiquitin related proteins in the pathology and development of schizophrenia: Evidence from human and animal studies Anti-inflammatory, antiapoptotic, and antioxidant activity of fluoxetine Chronic antidepressant treatments resulted in altered expression of genes involved in inflammation in the rat hypothalamus Cerebrospinal fluid inflammatory markers in patients with Listeria monocytogenes meningitis Correlation of serum CX3CL1 level with disease activity in adult-onset Still's disease and significant involvement in hemophagocytic syndrome Up-regulation of interleukin-23 induces persistent allodynia via CX3CL1 and interleukin-18 signaling in the rat spinal cord after tetanic sciatic stimulation Involvement of CX3CL1/CX3CR1 signaling in spinal long term potentiation Fractalkine Attenuates Microglial Cell Activation Induced by Prenatal Stress Interleukin-18-mediated microglia/astrocyte interaction in the spinal cord enhances neuropathic pain processing after nerve injury Role of the CX3CR1/p38 MAPK pathway in spinal microglia for the development of neuropathic pain following nerve injury-induced cleavage of fractalkine Involvement of microglia and interleukin-18 in the induction of long-term potentiation of spinal nociceptive responses induced by tetanic sciatic stimulation Changing face of microglia Microglia in neurodegenerative disease Interleukin 18 gene polymorphisms predict risk and outcome of Alzheimer's disease Thymic involution, a comorbidity factor in amyotrophic lateral sclerosis NLRP3 inflammasome is expressed by astrocytes in the SOD1 mouse model of ALS and in human sporadic ALS patients Interleukin 18 gene polymorphism is a risk factor for multiple sclerosis Association between promoter variants of interleukin-18 and schizophrenia in a Han Chinese population Comprehensive evaluation of positional candidates in the IL-18 pathway reveals suggestive associations with schizophrenia and herpes virus seropositivity Changes in the expression of genes related to neuroinflammation over the course of sporadic Alzheimer's disease progression: CX3CL1, TREM2, and PPARgamma CX3CL1/CX3CR1 in Alzheimer's Disease: A Target for Neuroprotection Downregulated Glia Interplay and Increased miRNA-155 as Promising Markers to Track ALS at an Early Stage Network analysis of human post-mortem microarrays reveals novel genes, microRNAs, and mechanistic scenarios of potential importance in fighting huntington's disease The role of fractalkine (CX3CL1) in regulation of CD4(+) cell migration to the central nervous system in patients with relapsing-remitting multiple sclerosis Fractalkine signaling regulates the inflammatory response in an α-synuclein model of Parkinson disease New evidences for fractalkine/CX3CL1 involved in substantia nigral microglial activation and behavioral changes in a rat model of Parkinson's disease CX3CL1 reduces neurotoxicity and microglial activation in a rat model of Parkinson's disease The soluble isoform of CX3CL1 is necessary for neuroprotection in a mouse model of Parkinson's disease Chemokines and chemokine receptors in mood disorders, schizophrenia, and cognitive impairment: a systematic review of biomarker studies Differential expression of brain immune genes and schizophrenia-related behavior in C57BL/6N and DBA/2J female mice Alternative activation of macrophages A unique microglia type associated with restricting development of Alzheimer's disease Induction of a common microglia gene expression signature by aging and neurodegenerative conditions: a co-expression meta-analysis A comparison of human and mouse gene co-expression networks reveals conservation and divergence at the tissue, pathway and disease levels Microglial priming in neurodegenerative disease We would like to thank Institute for Research in Fundamental Sciences (IPM) for providing us with research facilities.