# Autoimmunity and autoinflammation: A systems view on signaling pathway dysregulation profiles

Arsen Arakelyan, Nersisyan Lilit, David Poghosyan, Lusine Khondkaryan, Anna Hakobyan, Henry Lö Ffler-Wirth, Evie Melanitou, Hans Binder

<figure><img src="image_1.png" /></figure>

## To cite this version:

Arsen Arakelyan, Nersisyan Lilit, David Poghosyan, Lusine Khondkaryan, Anna Hakobyan, et al. Autoimmunity and autoinflammation: A systems view on signaling pathway dysregulation profiles. PLoS ONE, 2017, 12 (11), pp.e0187572. 10.1371/journal.pone.0187572. pasteur-01643667

HAL Id: pasteur-01643667

https://pasteur.hal.science/pasteur-01643667

Submitted on 21 Nov 2017

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH ARTICLE

# Autoimmunity and autoinflammation: A systems view on signaling pathway dysregulation profiles

Arsen Arakelyan¹,²*, Lilit Nersisyan¹,³, David Poghosyan⁴, Lusine Khondkaryan⁴, Anna Hakobyan¹, Henry Löffler-Wirth⁵, Evie Melanitou⁶, Hans Binder⁵

1. Bioinformatics Group, Institute of Molecular Biology, National Academy of Sciences RA, Yerevan, Armenia,

2. Department of Bioinformatics and Bioengineering, Russian-Armenian University, Yerevan, Armenia,

3. Zaven and Sonia Akian College of Science and Engineering, American University of Armenia, Yerevan, Armenia, 4. Group of Immune Response Regulation, Institute of Molecular Biology, National Academy of Sciences RA, Yerevan, Armenia, 5. Interdisciplinary Centre for Bioinformatics, University of Leipzig, Leipzig, Germany, 6. Department of Parasitology and Insect Vectors, Institut Pasteur, Paris, France

*aarakelyan@sci.am

<figure><img src="image_3.png" /></figure>

## Abstract

<figure><img src="image_4.png" /></figure>

### OPEN ACCESS

**Citation:** Arakelyan A, Nersisyan L, Poghosyan D, Khondkaryan L, Hakobyan A, Löffler-Wirth H, et al. (2017) Autoimmunity and autoinflammation: A systems view on signaling pathway dysregulation profiles. PLoS ONE 12(11): e0187572. https://doi.org/10.1371/journal.pone.0187572

**Editor:** Xia Li, College of Bioinformatics Science and Technology, CHINA

**Received:** May 30, 2017

**Accepted:** October 23, 2017

**Published:** November 3, 2017

**Copyright:** © 2017 Arakelyan et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This study was supported by the State Committee of Science MES RA projects no. SCS 13YR-1F0022, 16GE-025, and 15T-1F150 (to AA) and German Federal Ministry of Education and Science (BMBF) WTZ ARM II-010 and 01ZX1304A (to HB) and FFE-0034 (to HL-W).

**Competing interests:** The authors have declared that no competing interests exist.

## Introduction

Autoinflammatory and autoimmune disorders are characterized by aberrant changes in innate and adaptive immunity that may lead from an initial inflammatory state to an organ specific damage. These disorders possess heterogeneity in terms of affected organs and clinical phenotypes. However, despite the differences in etiology and phenotypic variations, they share genetic associations, treatment responses and clinical manifestations. The mechanisms involved in their initiation and development remain poorly understood, however the existence of some clear similarities between autoimmune and autoinflammatory disorders indicates variable degrees of interaction between immune-related mechanisms.

## Methods

Our study aims at contributing to a holistic, pathway-centered view on the inflammatory condition of autoimmune and autoinflammatory diseases. We have evaluated similarities and specificities of pathway activity changes in twelve autoimmune and autoinflammatory disorders by performing meta-analysis of publicly available gene expression datasets generated from peripheral blood mononuclear cells, using a bioinformatics pipeline that integrates Self Organizing Maps and Pathway Signal Flow algorithms along with KEGG pathway topologies.

## Results and conclusions

The results reveal that clinically divergent disease groups share common pathway perturbation profiles. We identified pathways, similarly perturbed in all the studied diseases, such as PI3K-Akt, Toll-like receptor, and NF-kappa B signaling, that serve as integrators of signals guiding immune cell polarization, migration, growth, survival and differentiation. Further, two clusters of diseases were identified based on specifically dysregulated pathways: one
gathering mostly autoimmune and the other mainly autoinflammatory diseases. Cluster sep-
aration was driven not only by apparent involvement of pathways implicated in adaptive
immunity in one case, and inflammation in the other, but also by processes not explicitly
related to immune response, but rather representing various events related to the formation
of specific pathophysiological environment. Thus, our data suggest that while all of the stud-
ied diseases are affected by activation of common inflammatory processes, disease-spe-
cific variations in their relative balance are also identified.

**Introduction**

Epidemiological studies provide increasing evidence for the rise in prevalence of immune-
related disorders including autoimmune and allergic diseases in Western countries [1–3]. It is
predicted that the incidence of chronic inflammatory disorders, particularly autoimmune dis-
eases, such as type 1 diabetes, Crohn’s disease, rheumatoid arthritis and multiple sclerosis, will
grow even more rapidly during the next several decades [1,4,5]. Moreover, chronic inflamma-
tion is being recognized as an important trigger and contributor to the development and pro-
gression of various other human complex diseases, such as certain cancers, atherosclerosis,
strokes and ischemic heart diseases, and even psychiatric disorders (schizophrenia and post-
traumatic stress disorder) [6–10]. Therefore, understanding the molecular mechanisms under-
lying the development of inflammatory disorders will have significant impact on public health.

Autoimmune diseases(AI) are characterized by dysfunction of the immune system leading
to loss of immune tolerance against self-tissues, by the presence of autoreactive T and B cells,
and by a complex pathogenesis of multifactorial etiology, whereas genetics and environmental
factors together are responsible for disease onset [11,12]. There are more than 80 such condi-
tions affecting susceptible human subjects [13]. Autoimmune conditions may be systemic (tar-
geting multiple organs and tissues), as is the case of lupus, or tissue specific, as is the case of
multiple sclerosis (against myelin) or type 1 diabetes (against pancreatic beta cells). The simul-
taneous presence of several autoimmune diseases is observed in some cases, pointing at the
possibility of a shared origin and/or mechanisms [14].

Autoinflammatory diseases (AIF) are a relatively new and expanding group of self-directed
inflammatory disorders, clinically described as periodic fever syndromes but also with epi-
sodes of acute inexplicable inflammation involving the innate immune system [15,16]. They
are characterized by inflammatory episodes at disease-prone sites, in the absence of autoreac-
tive T cells and high autoantibody titers [17,18]. Despite the differences in primary players,
they share common characteristics with AI diseases, such as self-tissue directed inflammation
in the absence of an obvious infectious trigger or injury. While in AIFs the innate immune sys-
tem directly causes tissue inflammation, in AIs the innate immune system activates the adap-
tive system and this later activates the inflammatory process [15]. The autoinflammatory
syndromes include a subset of hereditary conditions characterized by recurrent episodes of
fever and self-resolving attacks of systemic inflammation without microbial infection or auto-
immunity. There are several AIF associated with genetic mutations affecting the innate
immune system, the primary defense system against foreign antigens [19]. Such genetic muta-
tions have been identified affecting genes encoding for the tumor necrosis factor receptor
(TNFR1) as is the case for the autosomal dominant TNF receptor associated syndrome
(TRAPS) [20] or the gene encoding for mevalonate kinase (MVK) responsible for the AIDS
syndrome. Mutations of the *MVK* gene are responsible for an increase of mevalonic acid
concentrations, elevated levels of IgD in the serum and also an increased secretion of IL-1β [21,22]. Similarly several autoinflammatory syndromes are related with the gene encoding for pyrin [22]. The PAPA (pyogenic arthritis, pyoderma gangrenosum and acne) syndrome, a dominantly inherited autoinflammatory condition is associated with mutations in the gene encoding for proline serine threonine phosphatase-interacting protein (PSTPIP1) that interacts with pyrin [23]. Pyrin and PSTPIP1 proteins are associated with the cytoskeleton in myeloid/monocytes and their interaction contributes to increased IL-1β production, NFkB activation and apoptosis [24]. Finally the connection between pyrin and autoinflammation is observed in the Cryopyrin-Associated Periodic Syndromes (CAPS) in which mutations of genes encoding for the components of the proteins involved in the inflammasome (NLRP3) are implicated [25-27]. The phenotypic heterogeneity characteristic of AI and AIF diseases, does not necessarily reflect fundamental genetic or mechanistic differences between these groups. Indeed, while some gene variants and SNPs are specific to a particular disease [11,28], others predispose an individual to the development of multiple disorders, which indicates that shared mutations may affect genes or pathways implicated in the pathogeneses of several diseases [29,30]. Moreover, while the mechanisms and causes of development may differ in various diseases, the downstream effects occurring after the disease onset may be similar. This similarity can potentially underlie responsiveness of different types of diseases to the same type of treatment (e.g. glucocorticoids) [31-33] indicating the presence of shared drug targets. However, the global picture of molecular mechanisms underlying the similarities and specificities of chronic inflammatory disorders is not completely understood and is addressed in this study, using gene expression data to compare activation patterns of selected pathways in a subset of AI and AIF diseases.

Luckily, the advances in high-throughput biological data measurements, explosive growth of data in various public repositories and development of new bioinformatics algorithms for data analysis, already provide an opportunity to address the above mentioned issues from a systems biology viewpoint. Previously, we had developed a bioinformatics pipeline for pathway perturbation based analysis and similarity/specificity assessment of disease groups, and have applied it to a number of lung diseases [34]. The obtained results had confirmed the validity of our methodology, as well as had led us to draw new conclusions on pathological characteristics of lung diseases. In the present work, we have used our approach for global assessment of similarities and specificities of downstream molecular events of chronic inflammation through evaluation of pathway activity changes in autoimmune and autoinflammatory disorders.

We have combined several publicly available datasets for autoimmune and autoinflammatory diseases in an attempt to uncover their common and specific pathobiological features. Our goal was to include as many conditions as possible, while in the meantime minimizing heterogeneity of microarray platforms and source tissue used for transcriptome measurement. These criteria left us with four monogenic autoinflammatory, four polygenic autoinflammatory disorders sharing characteristics with autoinflammation and autoimmunity, and four polygenic autoimmune conditions (Table 1). Monogenic autoimmune disorders were not included in our study, as to our knowledge, respective transcriptome data meeting our selection criteria were not available. Even though the apparent heterogeneity of the disorders selected in our study presented a challenge for identification of common mechanisms, the existing similarities and, in particular, the common inflammatory component of these disorders, remained intriguing and encouraged us to undertake a comparative gene expression analysis of a large set of existing transcriptome data.

Our study offers a systems biology approach for comparative analysis of gene expression data from cohorts carrying apparent divergent diseases with common immune related
<table><caption>Table 1. Description of the samples and datasets included in the study.</caption><thead><tr><th>Disease</th><th>Type*</th><th>Abbrev.</th><th>GEO Acc.</th><th>Number of Cases</th><th>Number of Controls</th><th>Array Platform</th></tr></thead><tbody><tr><td>Type 1 Diabetes</td><td>Al</td><td>T1D</td><td>GSE55100</td><td>12</td><td>10</td><td>Affymetrix Human Genome U133 Plus 2.0 Array</td></tr><tr><td>Multiple sclerosis</td><td>Al</td><td>MS</td><td>GSE21942</td><td>12</td><td>15</td><td>Affymetrix Human Genome U133 Plus 2.0 Array</td></tr><tr><td>Systemic lupus erythematosus</td><td>Al</td><td>SLE</td><td>GSE50772</td><td>61</td><td>20</td><td>Affymetrix Human Genome U133 Plus 2.0 Array</td></tr><tr><td>Sjögren's syndrome</td><td>Al</td><td>SS</td><td>GSE48378</td><td>11</td><td>16</td><td>Affymetrix Human Exon 1.0 ST Array</td></tr><tr><td>Cryopyrin associated periodic syndrome</td><td>AIF</td><td>CAPS</td><td>GSE43553</td><td>23</td><td>20</td><td>Affymetrix Human Genome U133A 2.0 Array</td></tr><tr><td>Hyper IgD Syndrome (HIDS) (mutations in MVK gene)</td><td>AIF</td><td>MVK</td><td>GSE43553</td><td>8</td><td>20</td><td>Affymetrix Human Genome U133A 2.0 Array</td></tr><tr><td>PAPA Syndrome (mutations in PSTPIP1 gene)</td><td>AIF</td><td>PSTPIP1</td><td>GSE43553</td><td>6</td><td>20</td><td>Affymetrix Human Genome U133A 2.0 Array</td></tr><tr><td>TNF receptor associated periodic syndrome</td><td>AIF</td><td>TRAPS</td><td>GSE43553</td><td>29</td><td>20</td><td>Affymetrix Human Genome U133A 2.0 Array</td></tr><tr><td>Behcet's disease</td><td>PAIF</td><td>BD</td><td>GSE17114</td><td>15</td><td>14</td><td>Affymetrix Human Genome U133 Plus 2.0 Array</td></tr><tr><td>Crohn's disease</td><td>PAIF</td><td>CD</td><td>GSE3365</td><td>59</td><td>42</td><td>Affymetrix Human Genome U133A Array</td></tr><tr><td>Ulcerative colitis</td><td>PAIF</td><td>UC</td><td>GSE3365</td><td>26</td><td>42</td><td>Affymetrix Human Genome U133A Array</td></tr><tr><td>Juvenile Idiopathic Arthritis</td><td>PAIF</td><td>JIA</td><td>GSE67596</td><td>22</td><td>15</td><td>Affymetrix Human Genome U133 Plus 2.0 Array</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0187572.t001

manifestations. Despite the limitations discussed below, such studies may contribute to identification of common pathways in autoimmunity and autoinflammation, and have the potential to facilitate development of new experimental strategies for understanding the molecular mechanisms underlying the interplay between those processes.

## Material and methods

### Data sources

We have used data deposited at Gene Expression Omnibus (GEO), which is the largest repository for microarray gene expression studies [35]. The search was performed with keywords “autoinflammation”, “autoimmunity” and limited to human samples. We were aimed at selection of the most similar study designs (sample source tissue type, availability of controls) and microarray platforms. Thus, only datasets derived from peripheral blood mononuclear cells (PBMCs), as the most frequent sample source; and Affymetrix platforms, as the most prevalent ones, were chosen.

Our selection resulted in eight microarray datasets containing gene expression profiles of PBMC samples from patients with autoimmune and (or) autoinflammatory diseases, as well as healthy subjects. All of the platforms were Affymetrix microarrays, of which five were Human Genome U133 Plus 2.0, four were Human Genome U133A 2.0, two were Human Genome U133A, while only one was from Human Exome 1.0 ST Array platforms (Table 1). Data on Sjögren's Syndrome was the only dataset in our analysis that was not derived from Human Genome U133 series platforms, since there were no alternative datasets for this condition in

<footnote>* Al—autoimmunity, AIF—monogenic autoinflammation, PAIF—polygenic autoinflammation (20, 21)</footnote>
GEO at the time of conducting this study. The effect of possible systematic biases introduced by biological and technical variation in our dataset collection is addressed below.

Only untreated samples were included in the analyses. Dataset accessions, descriptions and sample counts are presented in Table 1. For detailed information on subject characteristics, disease stage, sample preparation and scan protocol and low level data processing, the reader is referred to the GEO database pages under sample accessions provided in Table 1.

## Dataset preprocessing

Raw Affymetrix CEL files have been downloaded for all the datasets, except for GSE3365 (gene expression in Crohn's disease and Ulcerative colitis), for which no raw data were available. Probe signal intensity conversions, RMA normalization and chip annotation were performed using the "affy" package for R [36]. The GSE48378 exon array dataset (gene expression in Sjögren's syndrome) was preprocessed using the "oligo" package [37], since "affy" is not intended for analysis of Affymetrix Human Exon arrays. The GSE3365 dataset was annotated based on platform information provided by GEO (GPL96: Affymetrix Human Genome U133A Array), and preprocessed using Affymetrix MAS5 algorithm.

All the genes with known Entrez IDs were annotated, and expression values for multiple probes of the same gene were averaged. The mean gene expression values of the controls within each dataset were used as reference to calculate the log fold changes (logFC) for each gene. The logFC values were then anti-logged to linear scale FC values.

## Pathway signal flow calculation

Pathway signal flow analysis was performed as described previously [38]. Briefly, 168 signaling and metabolic pathway maps were obtained from the KEGG pathways database, with the exclusion of disease and drug response pathways. The KEGG pathways are represented as graphs, where nodes are gene groups with similar functions and edges are interactions between them. The interactions are classified into "activation" and "inhibition" types. Each pathway has more than one input node, and more than one output (sink) node. The network of interactions that end with a single sink node is called a pathway branch. Note that each pathway may have many branches ending in specific sink nodes, where each may be associated with a different biological event or outcome. Thus, it is reasonable to evaluate activities of pathway branches individually [38].

The Pathway Signal Flow (PSF) algorithm for calculation of the activity values of pathway branches is described in detail in [39–41]. In short, PSF first assigns gene expression fold change (FC) values of the member genes to pathway nodes, and then sequentially updates the target node values for each source-target interaction: starting from the input nodes and finishing at the sinks. The target node values are updated depending on the type of interaction: 'activation' means that the greater the FC value at the source is, the greater the FC value at the target will be; whilst 'inhibition' has the opposite effect. In the end, the PSF algorithm returns the PSF values of the sink nodes for each pathway branch. The PSF values estimate the pathway activity fold change compared to the controls. Thus, values less than 1, are indicative for deactivation, while values greater than one for activation of the respective pathway branch.

We have evaluated in total 1825 sinks for the 168 pathways analyzed. Thus, each dataset is described with a vector of 1825 PSF values. Each sink is described by a profile of PSF values across all samples.

To estimate the possible systematic effect of biological and technical variation, we have substantially extended our analyses based on additional datasets (S1 File). We have demonstrated that there is a high correlation between PSF values obtained from datasets under abovementioned
condition variations. PSF analysis minimizes the systematic bias introduced by variation in experimental procedures, since PSF considers relative fold change values of gene expression in disease versus control.

## Pathway signal flow-self organizing maps analysis (PSF-SOM)

The PSF profiles described above were clustered and visualized with the Self Organizing Maps (SOM) algorithm, implemented in the R package "oposSOM" [42,43]. This algorithm was initially intended for application to gene expression data, but has also been applied for analysis of PSF values (PSF-SOM) [38]. The PSF-SOM algorithm takes as input the sample-wise PSF profiles of the 1825 sink nodes, and aggregates the similar profiles into meta-profiles, or meta-sinks. It arranges the meta-sinks on a two-dimensional grid in a way that correlated meta-sinks are located close to each other. In our analysis, we have obtained a 35x35 grid of meta-sinks. Using the meta-sink locations on the grid, we can then obtain a SOM portrait for each sample, where the mean PSF values of the meta-sinks are visualized with blue-green-to-red color scale (Fig 1). In other words, all the SOM portraits have the same arrangement of meta-sinks, but the PSF values differ based on the sample. Clusters of co-regulated meta-sinks are then detected as so-called spots in the SOM portraits using a detection threshold described in [42].

<figure><img src="image_5.png" /><figcaption>**Fig 1.** Disease specific SOM portraits. PSF profiles for each disease were mapped on a 35x35 SOM grid, and visualized as 2D maps, where colors indicate the actual state of pathway sink dysregulations. Red to green color gradient indicates up-regulation, while blue to green one indicates down-regulation.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0187572.g001
## Spot based disease similarity assessment

The similarity between diseases was assessed based on the number of shared spots between disease pairs. For this, we have constructed a graph object, with nodes representing the individual SOM-portraits of each disease, and edges connecting the disease pairs that share the highest number of spots: in case of ties, more than one edge was introduced. Community search within this graph was performed with the walktrap algorithm implemented in the R "igraph" package [44].

## Functional annotation of sinks

Functional annotation of sink dysregulation was assessed using overrepresentation approach implemented in the WebGestalt web-tool [45,46]. The pathway sink genes were tested against Gene ontology (GO) database. The minimum number of genes in a category was set to 5, and the significance threshold to p<0.05 after Benjamini & Hochberg multiple test correction. In order to achieve higher level of outlook on the biological processes conveyed by GO terms we performed semantic-similarity based summarization and removal of redundant terms using the REVIGO tool [47].

## Results and discussion

### PSF-SOM disease portraits

The SOM method transforms the multidimensional PSF data into a series of two-dimensional images, called "portraits", which visualize the activities of the pathway sink nodes in each studied disease (Fig 1). These PSF-activity portraits show blue and red spot like regions corresponding to down- and up-regulated sink nodes, respectively. It should be noted that the portraits depict changes in pathway activities compared to the healthy state.

At-a-glance examination of the spot distribution across disease portraits revealed considerable similarities, though some differences were observed as well: almost all the diseases were characterized by the presence of upregulated spot areas at the lower left corner, except for MVK. In addition, MS had a characteristic red spot in upper right corner, while MVK, CAPS, MS and to lesser extent PAPA (PSTPIP1) were characterized by the presence of an additional spot near the bottom right corner of corresponding SOM portraits. Finally, it should be noted that three monogenic autoinflammatory diseases (CAPS, TRAPS, MVK), as well as JIA and MS, and to lesser extent UC and CD, contained an additional upregulated spot near the lower left corner of their SOM portraits. In order to proceed with comparative analysis and functional annotation, all the up- and down-regulated spots from the averaged portraits of each disease were transferred into a single summary map for an overview (Fig 2). We have identified seven spot clusters assigned with capital letters A-G, each showing a specific profile of PSF-values in the studied diseases (Fig 2). The full report of PSF-SOM is available at [48].

It should be noted that oposSOM scales data individually, before drawing the expression portraits, but applies spot detection thresholds on the global summary map, using all of the datasets at once (Fig 2). Hence, there may be inconsistency between visual perception of spot presence on the individual SOM portraits (Fig 1), and their threshold-based detection results (Fig 2) [42,43].

### Disease similarity assessment

To assess similarities between the PSF portraits we have created a graph object where the nodes represent the diseases, and the edges connect the disease portraits sharing the highest
PAPA and other autoinflammatory disorders has also been reported in the literature. So, patients with PAPA syndrome are less responsive to anti-inflammatory treatment targeted against IL-1β and TNF alpha signaling [58]. Moreover, it has been shown that in a mouse model of PAPA, PSTPIP1 did not regulate inflammasome activation, suggesting alternative effects of PSTPIP1 mutations in PAPA [59]. Furthermore, spot C collects sinks related to “actin filament-based movement”, while mutations in PSTPIP1 lead to impaired podosome (actin-containing adhesive and invasive structure) formation [60]. All in all, this result indicates that PAPA, being an autoinflammatory syndrome, has specificities in molecular events that distinguish it from other autoinflammatory disorders.

Based on unexpected clustering of PAPA with autoimmune disorders, and the observed difference between autoinflammatory diseases described by inflammasome activation, we have looked into PSF signals of NOD-like receptor signaling pathway in all the studied conditions, since it collects inflammasome related signaling events. The results show expected clustering of MVK, PAPA, TRAPS and CAPS that were separate from other diseases (S2 File). This in turn indicates that the specificities observed for PAPA syndrome are not related to inflammasome formation, but rather to dysregulation of other pathways.

Cluster 2 collects CD, CAPS, UC, JIA, TRAPS and MVK and is specifically characterized by up-regulated spots C and D (Fig 3). This cluster is more homogeneous, and consists of monogenic (MVK, CAPS, TRAPS) and polygenic (JIA, CD, UC) autoinflammatory disorders [32,61].

<figure><img src="image_6.png" /><figcaption>**Fig 3.** The spot-similarity graph of the diseases. The nodes represent the diseases; the edges connect the disease portraits sharing the highest number of spots. The portraits are provided as thumbnails to each node. The graph is divided into the pink and the blue clusters using the walktrap community search algorithm, implemented in R package "igraph" [44].</figcaption></figure>

https://doi.org/10.1371/journal.pone.0187572.g003
Surprisingly, SLE, being an autoimmune disorder, also fell into this cluster. The systemic inflammatory nature of SLE whereas the immune system attacks healthy tissues in several parts of the body may explain the shared spots with the other autoinflammatory syndromes in the cluster [62]. Recent findings indicate that SLE shares molecular signatures with type I interferon-mediated monogenic autoinflammatory disorders [63]. Moreover, there is an indication on the emerging role of the inflammasome in SLE pathogenesis [64]. Recent studies by Shin *et al* [65] and Zhang *et al* [66] demonstrated that anti-self dsDNA antibodies in monocytes can induce the activation of NLRP3 inflammasome, leading to production of IL-1β. The latter in turn induces Th17 cell response in SLE. Another study by Kahlenberg *et al* [67] has found that caspase 1, the central enzyme involved in inflammasome activation, is essential for development of SLE in a mouse model of inducible lupus.

The obtained results suggest that regardless of the disease initiation event, the downstream pathway activity profiles share considerable similarities that can be attributed to chronic inflammation, which is an essential pathological event implicated in all the studied diseases. Although this is not surprising, a formal demonstration of shared inflammation-related pathways between these diseases suggests that the observed pathway dysregulations may create conditions for the development of comorbid syndromes.

## Functional context of dysregulated sinks and pathways

For functional annotation of pathway sinks associated with the spots of the summary map, we have performed GO term enrichment analysis using the WebGestalt and REVIGO programs [46,47]. According to the results, all spots were enriched with GO terms related to various aspects of immune response (Fig 4). The obtained results confirm the redundancy of main pathological events leading to development of chronic inflammation, which is indicated by the presence of spot D in all the diseases (except for MVK). Meanwhile, we have also observed some degree of specificity driven by the presence of specific spots, each accumulating distinct sets of processes that were unique for each disease cluster. Spot F relevant to the disease cluster 1 (autoimmunity) shows enrichment of GO terms associated with immune response, while spot C, which is unique to the disease cluster 2 (autoinflammation) is enriched with terms associated with inflammation.

Two other spots, B and E, seem to serve as supporting nodes that are either amplifiers (as, for example, in the case of CAPS or MS) or provide alternative mechanisms for immune response activation (as in the case of MVK). The complete list of GO terms associated with the pathway sinks and the enrichment p values can be found in the S3 File.

In total, 131 from 168 signaling and metabolic pathways with at least one significantly dysregulated sink in at least one disease were detected in the spots indicating that almost all pathways were dysregulated in the studied diseases (the complete list of pathway sinks is available in the S4 File). In other words, characterization of these diseases in terms of dysregulation of specific pathways seems inappropriate. Moreover, we have observed a large overlap of pathways in different spots (on average, 3 spots share sinks for a given pathway), suggesting an important role of pathway branching and the diversity of functional events associated with a single pathway (see S5 File). We have also noted that the proportion of affected branches in different pathways is markedly variable. The top pathways characterized with the highest proportion of dysregulated sinks are presented in the Table 2. Interestingly, pathway branches of top dysregulated pathways are not spread across all the spots, but are mainly concentrated in the spots C, D, E and F (Fig 5). All these pathways are directly linked to chronic inflammation, indicating their central role in the etiology of the studied diseases, as expected.

In contrast to the above mentioned pathways, PI3K signaling and Axon guidance pathways had dysregulated branches in all the spots. The PI3K-Akt signaling pathway essentially serves
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Functional annotation of biological processes associated with spots and their contribution to diseases. The pink (left) and blue (right) background indicate the two clusters of disease obtained by spot-similarity analysis. The sinks in each spot were annotated with Webgestalt and REVIGO programs to reveal associated GO term and summarize them with non-redundant descriptors.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0187572.g004

<table><caption>**Table 2.** Top pathways with dysregulated sinks.</caption><thead><tr><td>Pathway</td><td>Dysregulated/ Total sinks</td><td>Affected spots (% of dysregulated sinks in node)</td><td>Functional impact</td><td>References</td></tr></thead><tbody><tr><td>HIF-1 signaling pathway</td><td>31/33</td><td>A(3), E(3.5), F(90), G(3.5)</td><td>links hypoxia with inflammation; switch that modulates the immune system response by regulating metabolism and apoptosis depending on the microenvironment.</td><td>[96]</td></tr><tr><td>NF-kappa B signaling pathway</td><td>26/28</td><td>A(4), B(15), D(73), G(8)</td><td>a major inflammatory signal transduction route; participates in immunodeficiency, autoimmunity and autoinflammation, characterized by inflammation of several organs; non-canonical branch is associated with multi-organ autoinflammation</td><td>[97–101]</td></tr><tr><td>Toll-like receptor signaling pathway</td><td>18/19</td><td>D(94), E(6)</td><td>part of innate immune response toward various pathogens; recognizes endogenous agents which are released upon cell damage and necrosis; a link between innate and acquired immunity and is implicated in the development of both autoimmunity and autoinflammation</td><td>[[102–106]</td></tr><tr><td>ErbB signaling pathway</td><td>16/16</td><td>C(25), D(44), G(31)</td><td>implicated in cancer development; expression of receptors for this pathway has been reported on lymphocytes; activation leads to downstream activation of calcium, PI3K-Akt, mTOR signaling</td><td>[107–109]</td></tr><tr><td>Fc gamma R-mediated phagocytosis</td><td>13/14</td><td>A(8), D(15), E(62), F(15)</td><td>able to triggering both activating and inhibiting signals, leads to initiation or inhibition of a range of inflammation related events, including release of cytokines</td><td>[110]</td></tr><tr><td>GnRH signaling pathway</td><td>9/9</td><td>B(22), D(22), E(56)</td><td>Expression of GnRH and its receptors is detected in lymphocytes; leads to increased IL-2R; affects c-jun and ELK dependent gene expression and activation of arachidonic acid metabolism</td><td>[111,112]</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0187572.t002
<figure><img src="image_8.png" /><figcaption>**Fig 5. Sink distribution of top dysregulated pathways in spots.** The pathways with the highest proportions of dysregulated sinks are represented in the graph. The proportion of dysregulated sinks in each spot is depicted in the heights of respective bars. Each pathway is indicated with corresponding color.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0187572.g005

as an integrator of signals from different pathways and controls the number of key aspects of inflammatory response, such as cell movement, growth, survival and differentiation [68,69]. Furthermore, it is known that the PI3K pathway mediates interactions between polarizing signals and glucose metabolism, as well as respiratory burst in immune cells [70].

The axon guidance pathway is implicated in formation of neuronal networks [71]. However, this pathway frequently appears in lists of pathways dysregulated in immune system related conditions [72–74]. It includes mitogen activating kinases and genes encoding proteins involved in cytoskeleton rearrangement (KEGG Pathway ID hsa04360); both essential for PBMCs to react to stimuli [75]. Moreover, it has also been demonstrated that the axon guidance genes are involved in T-cell dependent B-cell maturation [76].

Along with the expected dysregulations of immune system signal transduction pathways, our analysis points out to massive dysregulations in metabolic pathways, including energy metabolism, metabolism of biomolecules, vitamins and xenobiotics. Our results are in line with previous observations on inflammation-related metabolic changes analyzed with GWAS and metabolomics approaches [77,78]. These data demonstrate that the mechanisms of the studied diseases expand far beyond the immune system dysregulations and affect pathways that regulate the basic cell functions. It should be noted, that in contrast to immune system related pathways, dysregulations in metabolic pathways were concentrated in specific spots: 28 out of 43 metabolic pathway dysregulations were spot specific (see S5 File), compared to 4 out of 14 immune system pathways.

In order to evaluate the relevance of the observed pathway deregulations in disease context, we have also evaluated the impact of treatment and disease outcome on pathway activation profiles in the studied diseases (S6 File). Strong association was observed between the known disease-drug gene sets and the spot-related pathway sink nodes identified in our analyses. The results suggest that the output nodes associated with identified pathway perturbations in this study may eventually serve as drug targets. Finally, we have also shown that the pathways commonly perturbed in PBMCs and in target tissues of organ-specific autoimmune disorders are partly resolved after disease treatment or during disease remission stages.

## Dysregulated pathway-infrastructure in autoimmune and autoinflammatory diseases

Considering that in both AI and AIF diseases the common characteristics lie in self-directed inflammation and disturbed homeostasis of canonical cytokine cascades, a comparative
analysis of the pathways involved could lead to transition to a systems level understanding of the observed phenomena. The results of previous studies based on analysis of genetic variations using genome-wide scans, gene expression analysis data and response to similar drugs, have already suggested that chronic inflammatory diseases share common mechanisms or pathways [32,79,80]. However, such studies are usually limited in the examined diseases, as well as in their attempt to obtain a general image of the mechanisms involved.

Here we have employed a systematic approach that demonstrates the universality of pathway dysregulations in autoimmune and autoinflammatory disorders, as well as uncovers novel aspects of similarity in the pathomechanisms of these diseases that were not reported previously. The results of this study clearly suggest that there is no single mechanism governing the development of the studied conditions. On one hand, the picture is rather complex and conditioned by a big number of variables. On the other hand, our pathway activity based disease clustering showed that diseases characterized by highly divergent clinical manifestations largely share the underlying pathomechanisms, at least at pathway activation levels. In some cases, disease linkage in the obtained clusters was supported by previous observations, and in others, the obtained results were novel. For example, we have addressed in our study the specifics of pathway dysregulation in PAPA [81] and HIDS (MVK) syndromes [82]. Previous studies have clearly implicated NLRP3 and pyrin inflammasomes in activating caspase-1 and IL-1β production in both syndromes [83–85]. Interestingly, while HIDS, in contrast to other autoinflammatory diseases, is caused by mutations in a metabolizing enzyme (MVK) [21,85,86], it was later demonstrated that mutations in the *MVK* gene disrupt the action of mevalonate kinase (MVK) resulting in a cascade of processes that involve pyrin inflammasome activation and enhanced IL-1β secretion [83,87]. This may partially explain the difference in spot enrichment observed in this disease. However, less is known about dysregulation states of other pathways involved in pathogenesis of these disorders. Our data suggest that PAPA is characterized by pathway activation patterns close to that of autoimmune diseases (Fig 3, cluster 1), while HIDS clusters with the classical autoinflammatory syndromes (Fig 3, cluster 2). Our analysis allowed us to move beyond the traditional gene-centered approach of description of molecular pathomechanisms, and to summarize the diseases at a systems biology level. Compared to the popular Gene Set Enrichment Analysis, PSF-SOM was able to mine far more functional information [38], providing a global picture of self-directed inflammatory syndromes and generate novel hypothesis for data validation to be undertaken in additional studies.

## Advantages and limitations of PSF-SOM analysis of PBMC

Our aim in this study was to compare diseases that share self-inflammatory phenotypes without prior restriction for taking in consideration genetic determinants specific to each disease. Therefore the most appropriate approach was the comparative analysis of existing microarray transcriptome data by our PSF-SOM method, allowing to tackle quantitative variations in pathways rather than in single genes. We appreciate that utilization of RNA sequencing data would further improve the analysis results, as it has a number of advantages over microarrays and even more valuable results could be obtained with the use of single cell transcriptome data. Unfortunately, current availability of RNA sequencing data for cell subpopulations is limited, and it is impossible to collect homogeneous data in terms of sequencing platform and sample types for all the studied diseases. Single cell data are even scarcer and still possess technological and algorithmic challenges [88]. However, our approach is technology independent and makes it possible to replicate the study once such data is available.

In our analysis we have considered only PBMC. While this choice offers data homogeneity, together with the facility of a non-invasive and easily obtainable tissue, it may reveal inflammatory-
related and/or other unrelated processes, reflecting a steady state of the corresponding organism including influences from a disease-unrelated environmental status. However, taking into consideration the concept of an immunological continuum, whereas tissue perturbations at the target sites of inflammation, rather than the immune system *per se*, are the key to disease expression, such perturbations may be reflected in the peripheral blood. Previous reports of transcriptome analysis in the peripheral blood have demonstrated that it is a well suited surrogate tissue, reflecting distal organismal perturbations and providing a large sensitive pool of gene transcripts with quantitative fluctuations detectable as gene expression modifications [89].

In order to show if the processes identified in PBMC at least to some extent reflect the processes activated in the target tissues, we have performed additional analysis to explore the overlap of pathway activities between PBMCs and target tissues for organ-specific autoimmune disorders, namely, diabetes type 1, multiple sclerosis, Crohn's disease and ulcerative colitis, where immune responses are directed against antigens present in a particular organ (S7 File). We have found that the activated pathways in PBMCs and in target tissues considerably overlap (34-59%, depending on the disease). This suggests that the processes identified in PBMC show common features with the disease characteristics of the tissue lesions.

It is worth pointing out that the results of this study should be interpreted with some degree of caution. The difference in the number of samples for each condition may lead to uneven variation of pathway activity values across datasets, which may bias the interpretation of the results. Another issue concerns the fact that according to the dataset descriptions available in the GEO, the samples were collected at different stages of disease progression, and not at the stage of initiation, thus making it perhaps more difficult to identify disease-specific determinants. Indeed, the initiating factors of each disease can be easier identified at the early stages, while those can be masked at later stages of the disease. This could lead to clustering of diseases otherwise being of different molecular origin. Nevertheless, the difficulty of collecting human tissues at early disease time points, due to the absence of early phenotypic landmarks and/or known early markers, renders such studies not feasible for the moment. A few studies have addressed this issue in both AI and AIF disorders [90-92]. One promising approach for the identification of non-invasive biomarkers is the exploration of microRNAs easily accessible in many body fluids [93].

## Summary and conclusions

The data presented in this report support the notion that diseases could be distinguished by pathways of the adaptive or innate immune responses, with the majority of conditions connecting by variable degrees of interaction between these two systems [94,95]. Our previous study on lung diseases using a similar systems biology approach has resulted in clear separation between cancer and other chronic lung diseases, each group characterized by a specific set of dysregulated pathways [38]. The separation we observed in that study was attributed to a distinct set of pathomechanisms implicated in cancers (cell proliferation, metabolism) and other lung diseases (immune/inflammatory response and fibrotic tissue remodelling) [38]. Here we have limited our scope to study only immune-system related diseases, and expected to get higher degree of interrelatedness at the level of pathway activity perturbations, compared to the previous study.

Overall, our study offers a systems biology approach addressing the complexity and redundancy of the immune system mechanisms implicated in the inflammatory origin of the pathogenesis of the studied diseases. Identification of pathways shared between different diseases, has the potential to generate novel hypotheses and elaborate methods to subvert the immune dysfunction by modulating networks through regulation of one gene or miRNAs. Moreover, it
may contribute to our understanding of similarities and differences in responses to similar treatment depending on specific patterns of pathway alterations. The selection of pathway level resolution provides a more aggregated level of information, as compared to single gene level analysis. We believe that both gene- and pathway- level analyses are important, one complementing the other.

**Supporting information**

**S1 File. Analysis for batch effects and platform independence in PSF calculations.**
(DOCX)

**S2 File. Analysis of inflammasome activation.**
(DOCX)

**S3 File. GO term enrichment.**
(XLSX)

**S4 File. Sink-spot association.**
(XLSX)

**S5 File. Distribution of pathway sinks across spots.**
(XLSX)

**S6 File. Treatment effect on pathway activation profiles.**
(DOCX)

**S7 File. Analysis of pathway activation profile overlaps between tissue specific and systemic response.**
(DOCX)

**Author Contributions**

**Conceptualization:** Arsen Arakelyan, Hans Binder.

**Formal analysis:** Arsen Arakelyan, Lilit Nersisyan, David Poghosyan, Lusine Khondkaryan, Anna Hakobyan, Evie Melanitou.

**Funding acquisition:** Arsen Arakelyan, Hans Binder.

**Investigation:** Lilit Nersisyan, Evie Melanitou.

**Methodology:** Arsen Arakelyan, Evie Melanitou, Hans Binder.

**Project administration:** Arsen Arakelyan.

**Software:** Arsen Arakelyan, Lilit Nersisyan, Anna Hakobyan, Henry Löffler-Wirth, Hans Binder.

**Supervision:** Arsen Arakelyan, Hans Binder.

**Visualization:** Lilit Nersisyan, Anna Hakobyan, Henry Löffler-Wirth, Hans Binder.

**Writing – original draft:** Arsen Arakelyan, Lilit Nersisyan, David Poghosyan, Lusine Khondkaryan, Anna Hakobyan, Henry Löffler-Wirth, Evie Melanitou, Hans Binder.

**Writing – review & editing:** Arsen Arakelyan, Lilit Nersisyan, David Poghosyan, Lusine Khondkaryan, Anna Hakobyan, Henry Löffler-Wirth, Evie Melanitou, Hans Binder.
References

1. Bach J-F. The Effect of Infections on Susceptibility to Autoimmune and Allergic Diseases. N Engl J Med. 2002; 347: 911-920. https://doi.org/10.1056/NEJMra020100 PMID: 12239261

2. Eaton WW, Pedersen MG, Atladóttir HÓ, Gregory PE, Rose NR, Mortensen PB. The prevalence of 30 ICD-10 autoimmune diseases in Denmark. Immunol Res. 2010; 47: 228-231. https://doi.org/10.1007/s12026-009-8153-2 PMID: 20066507

3. Mackay IR, Leskovsek N V, Rose NR. The odd couple: a fresh look at autoimmunity and immunodeficiency. J Autoimmun. 2010; 35: 199-205. https://doi.org/10.1016/j.jaut.2010.06.008 PMID: 20817405

4. Lerner A, Matthias T. Changes in intestinal tight junction permeability associated with industrial food additives explain the rising incidence of autoimmune disease. Autoimmun Rev. 2015; 14: 479-489. https://doi.org/10.1016/j.autrev.2015.01.009 PMID: 25676324

5. Lerner A, Matthias T. Rheumatoid arthritis–celiac disease relationship: Joints get that gut feeling. Autoimmun Rev. 2015; 14: 1038-1047. https://doi.org/10.1016/j.autrev.2015.07.007 PMID: 26190704

6. Diakos CI, Charles KA, McMillan DC, Clarke SJ. Cancer-related inflammation and treatment effectiveness. Lancet Oncol. 2014; 15: e493–503. https://doi.org/10.1016/S1470-2045(14)70263-3 PMID: 25281468

7. Golia E, Limongelli G, Natale F, Fimiani F, Maddaloni V, Pariggiano I, et al. Inflammation and Cardio-vascular Disease: From Pathogenesis to Therapeutic Target. Curr Atheroscler Rep. 2014; 16: 435. https://doi.org/10.1007/s11883-014-0435-z PMID: 25037581

8. Jin R, Yang G, Li G. Inflammatory mechanisms in ischemic stroke: role of inflammatory cells. J Leukoc Biol. 2010; 87: 779-789. https://doi.org/10.1189/jlb.1109766 PMID: 20130219

9. Mullenix PS, Andersen CA, Starnes BW, al. et, Vao J, Vague P. Atherosclerosis as inflammation. Ann Vasc Surg. Lippincott-Raven, Philadelphia; 2005; 19: 130-8. https://doi.org/10.1007/s10016-004-0153-z PMID: 15714382

10. Zakharyan R, Boyajyan A. Inflammatory cytokine network in schizophrenia. World J Biol Psychiatry. 2014; 15: 174-87. https://doi.org/10.3109/15622975.2013.830774 PMID: 24041158

11. Costenbader KH, Gay S, Alarcón-Riquelme ME, Iaccarino L, Doria A. Genes, epigenetic regulation and environmental factors: Which is the most relevant in developing autoimmune diseases? Autoimmun Rev. 2012; 11: 604-609. https://doi.org/10.1016/j.autrev.2011.10.022 PMID: 22041580

12. Zhernakova A, Withoff S, Wijmenga C. Clinical implications of shared genetics and pathogenesis in autoimmune diseases. Nat Rev Endocrinol. 2013; 9: 646-59. https://doi.org/10.1038/nrendo.2013.161 PMID: 23959365

13. Ramos PS, Shedlock AM, Langefeld CD. Genetics of autoimmune diseases: insights from population genetics. J Hum Genet. 2015; 60: 657-664. https://doi.org/10.1038/jhg.2015.94 PMID: 26223182

14. Zakka LR, Reche PA, Ahmed AR. The molecular basis for the presence of two autoimmune diseases occurring simultaneously—preliminary observations based on computer analysis. Autoimmunity. 2012; 45: 253–63. https://doi.org/10.3109/08916934.2011.632454 PMID: 22053914

15. Doria A, Zen M, Bettio S, Gatto M, Bassi N, Nalotto L, et al. Autoinflammation and autoimmunity: Bridging the divide. Autoimmun Rev. 2012; 12: 22-30. https://doi.org/10.1016/j.autrev.2012.07.018 PMID: 22878274

16. Kastner DL, Aksentijevich I, Goldbach-Mansky R, Aksentijevich I, Torosyan Y, Samuels J, et al. Auto-inflammatory disease reloaded: a clinical perspective. Cell. Elsevier; 2010; 140: 784-90. https://doi.org/10.1016/j.cell.2010.03.002 PMID: 20303869

17. Goldbach-Mansky R, Kastner DL. Autoinflammation: the prominent role of IL-1 in monogenic autoinflammatory diseases and implications for common illnesses. J Allergy Clin Immunol. 2009; 124: 1141-9-1. https://doi.org/10.1016/j.jaci.2009.11.016 PMID: 20004775

18. Kambe N, Nakamura Y, Saito M, Nishikomori R. The inflammasome, an innate immunity guardian, participates in skin urticarial reactions and contact hypersensitivity. Allergol Int. 2010; 59: 105-13. https://doi.org/10.2332/allergolint.09-RAI-0160 PMID: 20179416

19. Stehlik C, Reed JC. The PYRIN connection: novel players in innate immunity and inflammation. J Exp Med. The Rockefeller University Press; 2004; 200: 551-8. https://doi.org/10.1084/jem.20032234 PMID: 15353551

20. McDermott MF, Aksentijevich I, Galon J, McDermott EM, Ogunkolade BW, Centola M, et al. Germline Mutations in the Extracellular Domains of the 55 kDa TNF Receptor, TNFR1, Define a Family of Dominantly Inherited Autoinflammatory Syndromes. Cell. 1999; 97: 133-144. https://doi.org/10.1016/S0092-8674(00)80721-7 PMID: 10199409

21. Drenth JP, Cuisset L, Grateau G, Vasseur C, van de Velde-Visser SD, de Jong JG, et al. Mutations in the gene encoding mevalonate kinase cause hyper-IgD and periodic fever syndrome. International
Hyper-IgD Study Group. Nat Genet. 1999; 22: 178–181. https://doi.org/10.1038/9696 PMID: 10369262

22. Houten SM, Kuis W, Duran M, de Koning TJ, van Royen-Kerkhof A, Romeijn GJ, et al. Mutations in MVK, encoding mevalonate kinase, cause hyperimmunoglobulinaemia D and periodic fever syndrome. Nat Genet. 1999; 22: 175–177. https://doi.org/10.1038/9691 PMID: 10369261

23. Shoham NG, Centola M, Mansfield E, Hull KM, Wood G, Wise CA, et al. Pyrin binds the PSTPIP1/CD2BP1 protein, defining familial Mediterranean fever and PAPA syndrome as disorders in the same pathway. Proc Natl Acad Sci U S A. National Academy of Sciences; 2003; 100: 13501–6. https://doi.org/10.1073/pnas.2135380100 PMID: 14595024

24. Wise CA, Gillum JD, Seidman CE, Lindor NM, Veile R, Bashiardes S, et al. Mutations in CD2BP1 disrupt binding to PTP PEST and are responsible for PAPA syndrome, an autoinflammatory disorder. Hum Mol Genet. Oxford University Press; 2002; 11: 961–969. https://doi.org/10.1093/hmg/11.8.961 PMID: 11971877

25. Aganna E, Martinon F, Hawkins PN, Ross JB, Swan DC, Booth DR, et al. Association of mutations in the NALP3/CIAS1/PYPAF1 gene with a broad phenotype including recurrent fever, cold sensitivity, sensorineural deafness, and AA amyloidosis. Arthritis Rheum. 2002; 46: 2445–2452. https://doi.org/10.1002/art.10509 PMID: 12355493

26. Mortimer L, Moreau F, MacDonald JA, Chadee K. NLRP3 inflammasome inhibition is disrupted in a group of auto-inflammatory disease CAPS mutations. Nat Immunol. 2016; 17: 1176–86. https://doi.org/10.1038/ni.3538 PMID: 27548431

27. Neven B, Callebaut I, Prieur A-M, Feldmann J, Bodemer C, Lepore L, et al. Molecular basis of the spectral expression of CIAS1 mutations associated with phagocytic cell-mediated autoinflammatory disorders CINCA/NOMID, MWS, and FCU. Blood. 2004; 103: 2809–15. https://doi.org/10.1182/blood-2003-07-2531 PMID: 14630794

28. WTCC. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007; 447: 661–678. https://doi.org/10.1038/nature05911 PMID: 17554300

29. Castiblanco J, Arcos-Burgos M, Anaya J-M. What is next after the genes for autoimmunity? BMC Med. 2013; 11: 197. https://doi.org/10.1186/1741-7015-11-197 PMID: 24107170

30. Melanitou E, Fain P, Eisenbarth GS. Genetics of type 1A (immune mediated) diabetes. J Autoimmun. 2003; 21: 93–8. Available: http://www.ncbi.nlm.nih.gov/pubmed/12935776 PMID: 12935776

31. Baughman RP, Lower EE. Medical therapy of sarcoidosis. Semin Respir Crit Care Med. 2014; 35: 391–406. https://doi.org/10.1055/s-0034-1376401 PMID: 25007090

32. Ciccarelli F, De Martinis M, Ginaldi L. An update on autoinflammatory diseases. Curr Med Chem. Ben-tham Science Publishers; 2014; 21: 261–9. https://doi.org/10.2174/09298673113206660303 PMID: 24164192

33. Zhao L, Ma H, Jiang Z, Jiang Y, Ma N. Immunoregulation therapy changes the frequency of interleukin (IL)-22 + CD4 + T cells in systemic lupus erythematosus patients. Clin Exp Immunol. 2014; 177: 212–218. https://doi.org/10.1111/cei.12330 PMID: 24635166

34. Arakelyan A, Nersisyan L, Petrek M, Löffler-Wirth H, Binder H. Cartography of Pathway Signal Perturbations Identifies Distinct Molecular Pathomechanisms in Malignant and Chronic Lung Diseases. Front Genet. 2016; 7. https://doi.org/10.3389/fgene.2016.00079 PMID: 27200087

35. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002; 30: 207–210. Available: http://www.ncbi.nlm.nih.gov/pubmed/11752295 PMID: 11752295

36. Gautier L, Cope L, Bolstad BM, Irizarry RA. affy—analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004; 20: 307–15. https://doi.org/10.1093/bioinformatics/btg405 PMID: 14960456

37. Carvalho BS, Irizarry RA. A framework for oligonucleotide microarray preprocessing. Bioinformatics. 2010; 26: 2363–2367. https://doi.org/10.1093/bioinformatics/btq431 PMID: 20688976

38. Arakelyan A, Nersisyan L, Petrek M, Löffler-Wirth H, Binder H. Cartography of Pathway Signal Perturbations Identifies Distinct Molecular Pathomechanisms in Malignant and Chronic Lung Diseases. Front Genet. 2016; 7: 79. https://doi.org/10.3389/fgene.2016.00079 PMID: 27200087

39. Arakelyan A, Aslanyan L, Boyajyan A. High-throughput Gene Expression Analysis Concepts and Applications. In: iConcept Press Ltd., editor. Genomics II—Bacteria, Viruses and Metabolic Pathways. 2013. pp. 71–95. Available: https://www iconceptpress com/books/genomics-ii-bacteria-viruses-and-metabolic-pathways/14410061/

40. Nersisyan L, Johnson G, Riel-Mehan M, Pico A, Arakelyan A. PSFC: a Pathway Signal Flow Calculator App for Cytoscape. F1000Research. 2015; 4: 480. https://doi.org/10.12688/f1000research.6706.2 PMID: 26834984
41. Nersisyan L, Löffler-Wirth H, Arakelyan A, Binder H. Gene Set- and Pathway- Centered Knowledge Discovery Assigns Transcriptional Activation Patterns in Brain, Blood, and Colon Cancer: Int J Knowl Discov Bioinforma. IGI Global; 2014; 4: 46–69. https://doi.org/10.4018/IJKDB.2014070104

42. Wirth H, Löffler M, von Bergen M, Binder H. Expression cartography of human tissues using self organizing maps. BMC Bioinformatics. 2011; 12: 306. https://doi.org/10.1186/1471-2105-12-306 PMID: 21794127

43. Löffler-Wirth H, Kalcher M, Binder H. oposSOM: R-package for high-dimensional portraying of genome-wide expression landscapes on bioconductor: Fig 1. Bioinformatics. 2015; 31: 3225–3227. https://doi.org/10.1093/bioinformatics/btv342 PMID: 26063839

44. Csardi G, Nepusz T. The igraph software package for complex network research. InterJournal. 2006; Complex Sy: 1695. citeulike-article-id:3443126

45. Zhang B, Kirov S, Snoddy J. WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. 2005; 33: W741–8. https://doi.org/10.1093/nar/gki475 PMID: 15980575

46. Wang J, Duncan D, Shi Z, Zhang B. WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013. Nucleic Acids Res. 2013; 41: W77–W83. https://doi.org/10.1093/nar/gkt439 PMID: 23703215

47. Supek F, Bošnjak M, Škunca N, Šmuc T, Rivals I, Personnaz L, et al. REVIGO Summarizes and Visualizes Long Lists of Gene Ontology Terms. Gibas C, editor. PLoS One. Public Library of Science; 2011; 6: e21800. https://doi.org/10.1371/journal.pone.0021800 PMID: 21789182

48. Arakelyan A, Nersisyan L, Poghosyan D, Khondkaryan L, Löffler-Wirth H, Melanitou E, et al. PSF-SOM full results [Dataset]. Zenodo. 2017; https://doi.org/10.5281/ZENODO.291404

49. de Seze J, Devos D, Castelnovo G, Labauge P, Dubucquois S, Stojkovic T, et al. The prevalence of Sjögren syndrome in patients with primary progressive multiple sclerosis. Neurology. Lippincott Williams & Wilkins; 2001; 57: 1359–63. https://doi.org/10.1212/WNL.57.8.1359 PMID: 11673571

50. Jung SM, Lee BG, Joh GY, Cha JK, Chung WT, Kim KH. Primary Sjögren's syndrome manifested as multiple sclerosis and cutaneous erythematous lesions: a case report. J Korean Med Sci. Korean Academy of Medical Sciences; 2000; 15: 115–8. https://doi.org/10.3346/jkms.2000.15.1.115 PMID: 10719822

51. Solomon AJ, Hills W, Chen Z, Rosenbaum J, Bourdette D, Whitham R. Autoantibodies and Sjogren's Syndrome in multiple sclerosis, a reappraisal. Reindl M, editor. PLoS One. 2013; 8: e65385. https://doi.org/10.1371/journal.pone.0065385 PMID: 23776474

52. Binder A, Maddison PJ, Skinner P, Kurtz A, Isenberg DA. Sjögren's syndrome: association with type-1 diabetes mellitus. Br J Rheumatol. Oxford University Press; 1989; 28: 518–20. https://doi.org/10.1093/RHEUMATOLOGY/28.6.518 PMID: 2590806

53. Brayer J, Lowry J, Cha S, Robinson CP, Yamachika S, Peck AB, et al. Alleles from chromosomes 1 and 3 of NOD mice combine to influence Sjögren's syndrome-like autoimmune exocrinopathy. J Rheumatol. 2000; 27: 1896–904. Available: http://www.ncbi.nlm.nih.gov/pubmed/10955330 PMID: 10955330

54. Todd JA, Wicker LS. Genetic Protection from the Inflammatory Disease Type 1 Diabetes in Humans and Animal Models. Immunity. 2001; 15: 387–395. https://doi.org/10.1016/S1074-7613(01)00202-3 PMID: 11567629

55. Ashjazadeh N, Borhani Haghighi A, Samangooie S, Moosavi H. Neuro-Behcet's disease: a masquerader of multiple sclerosis. A prospective study of neurologic manifestations of Behcet's disease in 96 Iranian patients. Exp Mol Pathol. 2003; 74: 17–22. PMID: 12645628

56. Song GG, Kim J-H, Lee YH. The chemokine receptor 5 delta32 polymorphism and type 1 diabetes, Behcet's disease, and asthma: a meta-analysis. Immunol Invest. 2014; 43: 123–36. https://doi.org/10.3109/08820139.2013.847457 PMID: 24171669

57. Zapatero DC, Colin IM. Behcet disease and type 1 diabetes, a fortuitous association? Louv Med. 2008; 127.

58. Demidowich AP, Freeman AF, Kuhns DB, Aksentijevich I, Gallin JI, Turner ML, et al. Brief report: genotype, phenotype, and clinical course in five patients with PAPA syndrome (pyogenic sterile arthritis, pyoderma gangrenosum, and acne). Arthritis Rheum. NIH Public Access; 2012; 64: 2022–7. https://doi.org/10.1002/art.34332 PMID: 22161697

59. Wang D, Hoing S, Patterson HC, Ahmad UM, Rathinam VAK, Rajewsky K, et al. Inflammation in Mice Ectopically Expressing Human Pyogenic Arthritis, Pyoderma Gangrenosum, and Acne (PAPA) Syndrome-associated PSTPIP1 A230T Mutant Proteins. J Biol Chem. 2013; 288: 4594–4601. https://doi.org/10.1074/jbc.M112.443077 PMID: 23293022

60. Cortesio CL, Wernimont SA, Kastner DL, Cooper KM, Huttenlocher A. Impaired podosome formation and invasive migration of macrophages from patients with a PSTPIP1 mutation and PAPA syndrome.
Arthritis Rheum. NIH Public Access; 2010; 62: 2556-8. https://doi.org/10.1002/art.27521 PMID: 20506269

61. McGonagle D, McDermott MF. A proposed classification of the immunological diseases. PLoS Med. Public Library of Science; 2006; 3: e297. https://doi.org/10.1371/journal.pmed.0030297 PMID: 16942393

62. Lisnevskaia L, Murphy G, Isenberg D. Systemic lupus erythematosus. Lancet (London, England). 2014; 384: 1878-88. https://doi.org/10.1016/S0140-6736(14)60128-8

63. Rodero MP, Crow YJ. Type I interferon-mediated monogenic autoinflammation: The type I interferonopathies, a conceptual overview. J Exp Med. 2016; 213: 2527-2538. https://doi.org/10.1084/jem.20161596 PMID: 27821552

64. Kahlenberg JM, Kaplan MJ. The inflammasome and lupus. Curr Opin Rheumatol. 2014; 26: 475-481. https://doi.org/10.1097/BOR.0000000000000088 PMID: 24992143

65. Shin MS, Kang Y, Lee N, Wahl ER, Kim SH, Kang KS, et al. Self double-stranded (ds)DNA induces IL-1β production from human monocytes by activating NLRP3 inflammasome in the presence of anti-dsDNA antibodies. J Immunol. NIH Public Access; 2013; 190: 1407-15. https://doi.org/10.4049/jimmunol.1201195 PMID: 23315075

66. Zhang H, Fu R, Guo C, Huang Y, Wang H, Wang S, et al. Anti-dsDNA antibodies bind to TLR4 and activate NLRP3 inflammasome in lupus monocytes/macrophages. J Transl Med. 2016; 14: 156. https://doi.org/10.1186/s12967-016-0911-z PMID: 27250627

67. Kahlenberg JM, Yalavarthi S, Zhao W, Hodgin JB, Reed TJ, Tsuji NM, et al. An essential role of caspase 1 in the induction of murine lupus and its associated vascular damage. Arthritis Rheumatol (Hoboken, NJ). 2014; 66: 152-62. https://doi.org/10.1002/art.38225 PMID: 24449582

68. Hawkins PT, Stephens LR. PI3K signalling in inflammation. Biochim Biophys Acta. 2015; 1851: 882-97. https://doi.org/10.1016/j.bbalip.2014.12.006 PMID: 25514767

69. Xie S, Chen M, Yan B, He X, Chen X, Li D, et al. Identification of a Role for the PI3K/AKT/mTOR Signaling Pathway in Innate Immune Cells. Ryffel B, editor. PLoS One. Public Library of Science; 2014; 9: e94496. https://doi.org/10.1371/journal.pone.0094496 PMID: 24718556

70. Covarrubias AJ, Aksoylar HI, Horng T. Control of macrophage metabolism and activation by mTOR and Akt signaling. Semin Immunol. 2015; 27: 286-296. https://doi.org/10.1016/j.smim.2015.08.001 PMID: 26360589

71. Bashaw GJ, Klein R. Signaling from axon guidance receptors. Cold Spring Harb Perspect Biol. Cold Spring Harbor Laboratory Press; 2010; 2: a001941. https://doi.org/10.1101/cshperspect.a001941 PMID: 20452961

72. Dorr C, Wu B, Guan W, Muthusamy A, Sanghavi K, Schladt DP, et al. Differentially expressed gene transcripts using RNA sequencing from the blood of immunosuppressed kidney allograft recipients. PLoS One. Public Library of Science; 2015; 10: e0125045. https://doi.org/10.1371/journal.pone.0125045 PMID: 25946140

73. Gupta P, Liu B, Wu JQ, Soriano V, Vispo E, Carroll AP, et al. Genome-wide mRNA and miRNA analysis of peripheral blood mononuclear cells (PBMC) reveals different miRNAs regulating HIV/HCV coinfection. Virology. 2014; 450-451: 336-349. https://doi.org/10.1016/j.virol.2013.12.026 PMID: 24503097

74. Mutez E, Larvor L, Leprêtre F, Mouroux V, Hamalek D, Kerckaert J-P, et al. Transcriptional profile of Parkinson blood mononuclear cells with LRRK2 mutation. Neurobiol Aging. 2011; 32: 1839-1848. https://doi.org/10.1016/j.neurobiolaging.2009.10.016 PMID: 20096956

75. Sánchez-Madrid F, del Pozo MA. Leukocyte polarization in cell migration and immune interactions. EMBO J. European Molecular Biology Organization; 1999; 18: 501-11. https://doi.org/10.1093/emboj/18.3.501 PMID: 9927410

76. Yu D, Cook MC, Shin D-M, Silva DG, Marshall J, Toellner K-M, et al. Axon growth and guidance genes identify T-dependent germinal centre B cells. Immunol Cell Biol. Nature Publishing Group; 2008; 86: 3-14. https://doi.org/10.1038/sj.icb.7100123 PMID: 17938642

77. Fitzpatrick M, Young S. Metabolomics—a novel window into inflammatory disease. Swiss Med Wkly. 2013; https://doi.org/10.4414/smw.2013.13743 PMID: 23348753

78. Khor B, Gardet A, Xavier RJ. Genetics and pathogenesis of inflammatory bowel disease. Nature. NIH Public Access; 2011; 474: 307-17. https://doi.org/10.1038/nature10209 PMID: 21677747

79. Chimenti MS, Triggianese P, Conigliaro P, Candi E, Melino G, Perricone R. The interplay between inflammation and metabolism in rheumatoid arthritis. Cell Death Dis. 2015; 6: e1887. https://doi.org/10.1038/cddis.2015.246 PMID: 26379192
80. Lee RW, Nicholson LB, Sen HN, Chan C-C, Wei L, Nussenblatt RB, et al. Autoimmune and autoinflammatory mechanisms in uveitis. Semin Immunopathol. 2014; 36: 581-594. https://doi.org/10.1007/s00281-014-0433-9 PMID: 24858699

81. Wise CA, Gillum JD, Seidman CE, Lindor NM, Veile R, Bashiardes S, et al. Mutations in CD2BP1 disrupt binding to PTP PEST and are responsible for PAPA syndrome, an autoinflammatory disorder. Hum Mol Genet. Oxford University Press; 2002; 11: 961-969. https://doi.org/10.1093/hmg/11.8.961 PMID: 11971877

82. Mulders-Manders CM, Simon A. Hyper-IgD syndrome/mevalonate kinase deficiency: what is new? Semin Immunopathol. Springer; 2015; 37: 371-6. https://doi.org/10.1007/s00281-015-0492-6 PMID: 25990874

83. Wilson SP, Cassel SL. Inflammasome-mediated autoinflammatory disorders. Postgrad Med. NIH Public Access; 2010; 122: 125-33. https://doi.org/10.3810/pgm.2010.09.2209 PMID: 20861596

84. Xu H, Yang J, Gao W, Li L, Li P, Zhang L, et al. Innate immune sensing of bacterial modifications of Rho GTPases by the Pyrin inflammasome. Nature. 2014; 513: 237-241. https://doi.org/10.1038/nature13449 PMID: 24919149

85. Park YH, Wood G, Kastner DL, Chae JJ. Pyrin inflammasome activation and RhoA signaling in the autoinflammatory diseases FMF and HIDS. Nat Immunol. 2016; 17: 914-921. https://doi.org/10.1038/ni.3457 PMID: 27270401

86. Waterham HR, Poll-The BT, Houten SM, Kuis W, Duran M, de Koning TJ, et al. Mutations in MVK, encoding mevalonate kinase, cause hyperimmunoglobulinaemia D and periodic fever syndrome. Nat Genet. 1999; 22: 175-177. https://doi.org/10.1038/9691 PMID: 10369261

87. Favier LA, Schulert GS. Mevalonate kinase deficiency: current perspectives. Appl Clin Genet. Dove Press; 2016; 9: 101-10. https://doi.org/10.2147/TACG.S93933 PMID: 27499643

88. Liu S, Trapnell C. Single-cell transcriptome sequencing: recent advances and remaining challenges. F1000Research. Faculty of 1000 Ltd; 2016; 5. https://doi.org/10.12688/f1000research.7223.1 PMID: 26949524

89. Liew C-C, Ma J, Tang H-C, Zheng R, Dempsey AA. The peripheral blood transcriptome dynamically reflects system wide biology: a potential diagnostic tool. J Lab Clin Med. 2006; 147: 126-132. https://doi.org/10.1016/j.lab.2005.10.005 PMID: 16503242

90. Tektonidou MG, Ward MM. Validation of new biomarkers in systemic autoimmune diseases. Nat Rev Rheumatol. 2011; 7: 708-717. https://doi.org/10.1038/nrrheum.2011.157 PMID: 22045310

91. Maecker HT, Lindstrom TM, Robinson WH, Utz PJ, Hale M, Boyd SD, et al. New tools for classification and monitoring of autoimmune diseases. Nat Rev Rheumatol. 2012; 8: 317-328. https://doi.org/10.1038/nrrheum.2012.66 PMID: 22647780

92. Melanitou E, Tekaia F, Yeramian E. Investigation of secreted protein transcripts as early biomarkers for type 1 diabetes in the mouse model [Internet]. Gene. 2013. https://doi.org/10.1016/j.gene.2012.09.055 PMID: 23031813

93. Gallo A, Tandon M, Alevizos I, Illei GG, Pritchard C. The Majority of MicroRNAs Detectable in Serum and Saliva Is Concentrated in Exosomes. Afarinkia K, editor. PLoS One. Public Library of Science; 2012; 7: e30679. https://doi.org/10.1371/journal.pone.0030679 PMID: 22427800

94. Aziz A, Wong CH, Emery P, McDermott MF. Comparison of autoinflammation and autoimmunity using TRAPS and rheumatoid arthritis as prototypes of these conditions. Int J Clin Rheumtol. 2009; 4: 681-695.

95. McGonagle D, Savic S, McDermott MF. The NLR network and the immunological disease continuum of adaptive and innate immune-mediated inflammation against self. Semin Immunopathol. 2007; 29: 303-313. https://doi.org/10.1007/s00281-007-0084-1 PMID: 17805542

96. Eltzschig HK, Carmeliet P. Hypoxia and inflammation. N Engl J Med. 2011; 364: 656-65. https://doi.org/10.1056/NEJMra0910283 PMID: 21323543

97. Gavanescu I, Benoist C, Mathis D. B cells are required for Aire-deficient mice to develop multi-organ autoinflammation: A therapeutic approach for APECED patients. Proc Natl Acad Sci U S A. 2008; 105: 13009-14. https://doi.org/10.1073/pnas.0806874105 PMID: 18755889

98. Wessel A, Hsu A, Goldbach-Mansky R, Zilberman-Rudenko J, Siegel R, Hanson E. Inflammatory disease due to selectively impaired NF-κB activation and type I interferon response resulting from a de novo human NEMO hypomorphic mutation (HUM2P.333). J Immunol. American Association of Immunologists; 2014; 192: 53.6-53.6.

99. Sun S-C, Chang J-H, Jin J. Regulation of nuclear factor-κB in autoimmunity. Trends Immunol. NIH Public Access; 2013; 34: 282-9. https://doi.org/10.1016/j.it.2013.01.004 PMID: 23434408
100. Sun S-C, Chang J-H, Jin J, Hayden MS, Ghosh S, Sun SC, et al. Regulation of nuclear factor-κB in autoimmunity. Trends Immunol. Elsevier; 2013; 34: 282–289. https://doi.org/10.1016/j.it.2013.01.004 PMID: 23434408

101. Lawrence T. The nuclear factor NF-kappaB pathway in inflammation. Cold Spring Harb Perspect Biol. Cold Spring Harbor Laboratory Press; 2009; 1: a001651. https://doi.org/10.1101/cshperspect. a001651 PMID: 20457564

102. Singh D, Naik S. The role of toll-like receptors in autoimmune diseases. J Indian Rheumatol Assoc. 2005; 13: 162–165.

103. Pradhan VD, Das S, Surve P, Ghosh K. Toll-like receptors in autoimmunity with special reference to systemic lupus erythematosus. Indian J Hum Genet. Medknow Publications; 2012; 18: 155–60. https://doi.org/10.4103/0971-6866.100750 PMID: 23162288

104. Akira S, Takeda K, Kaisho T. Toll-like receptors: critical proteins linking innate and acquired immunity. Nat Immunol. Nature Publishing Group; 2001; 2: 675–680. https://doi.org/10.1038/90609 PMID: 11477402

105. Netea MG, van de Veerdonk FL, Kullberg BJ, Van der Meer JWM, Joosten LAB. The role of NLRs and TLRs in the activation of the inflammasome. Expert Opin Biol Ther. Taylor & Francis; 2008; 8: 1867–1872. https://doi.org/10.1517/14712590802494212 PMID: 18990074

106. Martinon F, Aksentijevich I. New players driving inflammation in monogenic autoinflammatory diseases. Nat Rev Rheumatol. 2015; 11: 11–20. https://doi.org/10.1038/nrrheum.2014.158 PMID: 25247411

107. Lu N, Wang L, Cao H, Liu L, Van Kaer L, Washington MK, et al. Activation of the epidermal growth factor receptor in macrophages regulates cytokine production and experimental colitis. J Immunol. NIH Public Access; 2014; 192: 1013–23. https://doi.org/10.4049/jimmunol.1300133 PMID: 24391216

108. Menon R, Farina C. Shared molecular and functional frameworks among five complex human disorders: a comparative study on interactomes linked to susceptibility genes. PLoS One. Public Library of Science; 2011; 6: e18660. https://doi.org/10.1371/journal.pone.0018660 PMID: 21533026

109. Palmer C, Diehn M, Alizadeh AA, Brown PO, Parks D, Herzenberg L, et al. Cell-type specific gene expression profiles of leukocytes in human peripheral blood. BMC Genomics. BioMed Central; 2006; 7: 115. https://doi.org/10.1186/1471-2164-7-115 PMID: 16704732

110. Nimmerjahn F, Ravetch J V. Fcγ receptors as regulators of immune responses. Nat Rev Immunol. Nature Publishing Group; 2008; 8: 34–47. https://doi.org/10.1038/nri2206 PMID: 18064051

111. Chen HF, Jeung EB, Stephenson M, Leung PC. Human peripheral blood mononuclear cells express gonadotropin-releasing hormone (GnRH), GnRH receptor, and interleukin-2 receptor gamma-chain messenger ribonucleic acids that are regulated by GnRH in vitro. J Clin Endocrinol Metab. 1999; 84: 743–50. https://doi.org/10.1210/jcem.84.2.5440 PMID: 10022447

112. Kraus S, Naor Z, Seger R. Intracellular signaling pathways mediated by the gonadotropin-releasing hormone (GnRH) receptor. Arch Med Res. 32: 499–509. PMID: 11750725