Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

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

  • Arsen Arakelyan ,

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Software, Supervision, Writing – original draft, Writing – review & editing

    aarakelyan@sci.am

    Affiliations Bioinformatics Group, Institute of Molecular Biology, National Academy of Sciences RA, Yerevan, Armenia, Department of Bioinformatics and Bioengineering, Russian-Armenian University, Yerevan, Armenia

  • Lilit Nersisyan,

    Roles Formal analysis, Investigation, Software, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Bioinformatics Group, Institute of Molecular Biology, National Academy of Sciences RA, Yerevan, Armenia, Zaven and Sonia Akian College of Science and Engineering, American University of Armenia, Yerevan, Armenia

  • David Poghosyan,

    Roles Formal analysis, Writing – original draft, Writing – review & editing

    Affiliation Group of Immune Response Regulation, Institute of Molecular Biology, National Academy of Sciences RA, Yerevan, Armenia

  • Lusine Khondkaryan,

    Roles Formal analysis, Writing – original draft, Writing – review & editing

    Affiliation Group of Immune Response Regulation, Institute of Molecular Biology, National Academy of Sciences RA, Yerevan, Armenia

  • Anna Hakobyan,

    Roles Formal analysis, Software, Visualization, Writing – original draft, Writing – review & editing

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

  • Henry Löffler-Wirth,

    Roles Software, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Interdisciplinary Centre for Bioinformatics, University of Leipzig, Leipzig, Germany

  • Evie Melanitou,

    Roles Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Department of Parasitology and Insect Vectors, Institut Pasteur, Paris, France

  • Hans Binder

    Roles Conceptualization, Funding acquisition, Methodology, Software, Supervision, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Interdisciplinary Centre for Bioinformatics, University of Leipzig, Leipzig, Germany

Abstract

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 separation 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 studied diseases are affected by activation of common inflammatory processes, disease-specific 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 [13]. It is predicted that the incidence of chronic inflammatory disorders, particularly autoimmune diseases, 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 inflammation is being recognized as an important trigger and contributor to the development and progression 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) [610].Therefore, understanding the molecular mechanisms underlying 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 conditions affecting susceptible human subjects [13]. Autoimmune conditions may be systemic (targeting 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 simultaneous 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 episodes 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 autoreactive 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 system directly causes tissue inflammation, in AIs the innate immune system activates the adaptive 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 autoimmunity. There are several AIF associated with genetic mutations affecting the innate immune system, the primary defense system against foreign antigens [19]. Such genetic mutations 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 HIDS 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 [2527].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) [3133] 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.

thumbnail
Table 1. Description of the samples and datasets included in the study.

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

Our study offers a systems biology approach for comparative analysis of gene expression data from cohorts carrying apparent divergent diseases with common immune related 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 133A, 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 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 [3941]. 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 de-activation, 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].

thumbnail
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.

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].

thumbnail
Fig 2. Global landscape of pathway dysregulations in diseases.

The summary SOM portrait with spots transferred from individual disease portraits. The seven spots detected in the SOM portraits are assigned the letters A-G (in the center of figure). PSF-activity profile barplots (left and right panels) represent the mean PSF-signal for a given spot across the diseases. If that value passes the defined threshold the spot is marked with a “+” sign for up-regulation and “-” sign for downregulation. Bar coloring indicates the relation of diseases to similarity clusters (see Fig 3).

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

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 number of spots. Based on a walktrap community detection search on this graph structure, we have identified two clusters or communities of highly similar diseases (Fig 3).

thumbnail
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].

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

Cluster 1 contains four polygenic autoimmune diseases studied i.e. T1D, BD, SS and MS, as well as one autoinflammatory syndrome, the PAPA (PSTPIP1). This cluster is characterized by up-regulated spots D and F and collects mostly autoimmune diseases, which have been shown to largely share symptoms and clinical manifestations. Several published case reports documented that primary SS manifestations can mimic those of MS [4951]. Moreover, SS frequently accompanies diabetes in humans [52] and in the T1D animal model of Non-Obese Diabetic (NOD) mouse [53]. Studies in the NOD mice have identified that T1D resistance genes in specific chromosomal loci (Idd3 and Idd5) protected the organism from inflammation and also of dysfunction of the salivary glands [53]. This indicates that possibly a general mechanism exists for maintaining a self-tolerant anti-inflammatory state and it affects more than one AI disease [54]. MS shares similarities with neurological manifestations of BD [55]. BD in turn has been reported to accompany T1D [56,57].

Notably, despite the fact that PAPA is an autoinflammatory syndrome and mutations in PSTPIP1 are thought to be linked to inflammasome activation and IL-1β oveproduction, PAPA lacks spot C and clusters with autoimmune diseases. The observed difference between 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].

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.

thumbnail
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.

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

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.

thumbnail
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.

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

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 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 [7274]. 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 [8385]. 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 [9092]. 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.

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

(DOCX)

S2 File. Analysis of inflammasome activation.

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

(DOCX)

S5 File. Distribution of pathway sinks across spots.

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

(XLSX)

S6 File. Treatment effect on pathway activation profiles.

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

(DOCX)

S7 File. Analysis of pathway activation profile overlaps between tissue specific and systemic response.

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

(DOCX)

References

  1. 1. Bach J-F. The Effect of Infections on Susceptibility to Autoimmune and Allergic Diseases. N Engl J Med. 2002;347: 911–920. pmid:12239261
  2. 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. pmid:20066507
  3. 3. Mackay IR, Leskovsek N V, Rose NR. The odd couple: a fresh look at autoimmunity and immunodeficiency. J Autoimmun. 2010;35: 199–205. pmid:20817405
  4. 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. pmid:25676324
  5. 5. Lerner A, Matthias T. Rheumatoid arthritis–celiac disease relationship: Joints get that gut feeling. Autoimmun Rev. 2015;14: 1038–1047. pmid:26190704
  6. 6. Diakos CI, Charles KA, McMillan DC, Clarke SJ. Cancer-related inflammation and treatment effectiveness. Lancet Oncol. 2014;15: e493–503. pmid:25281468
  7. 7. Golia E, Limongelli G, Natale F, Fimiani F, Maddaloni V, Pariggiano I, et al. Inflammation and Cardiovascular Disease: From Pathogenesis to Therapeutic Target. Curr Atheroscler Rep. 2014;16: 435. pmid:25037581
  8. 8. Jin R, Yang G, Li G. Inflammatory mechanisms in ischemic stroke: role of inflammatory cells. J Leukoc Biol. 2010;87: 779–789. pmid:20130219
  9. 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. pmid:15714382
  10. 10. Zakharyan R, Boyajyan A. Inflammatory cytokine network in schizophrenia. World J Biol Psychiatry. 2014;15: 174–87. pmid:24041158
  11. 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. pmid:22041580
  12. 12. Zhernakova A, Withoff S, Wijmenga C. Clinical implications of shared genetics and pathogenesis in autoimmune diseases. Nat Rev Endocrinol. 2013;9: 646–59. pmid:23959365
  13. 13. Ramos PS, Shedlock AM, Langefeld CD. Genetics of autoimmune diseases: insights from population genetics. J Hum Genet. 2015;60: 657–664. pmid:26223182
  14. 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. pmid:22053914
  15. 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. pmid:22878274
  16. 16. Kastner DL, Aksentijevich I, Goldbach-Mansky R, Aksentijevich I, Torosyan Y, Samuels J, et al. Autoinflammatory disease reloaded: a clinical perspective. Cell. Elsevier; 2010;140: 784–90. pmid:20303869
  17. 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. pmid:20004775
  18. 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. pmid:20179416
  19. 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. pmid:15353551
  20. 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. pmid:10199409
  21. 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. pmid:10369262
  22. 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. pmid:10369261
  23. 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. pmid:14595024
  24. 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. pmid:11971877
  25. 25. Aganna E, Martinon F, Hawkins PN, Ross JB, Swan DC, Booth DR, et al. Association of mutations in theNALP3/CIAS1/PYPAF1 gene with a broad phenotype including recurrent fever, cold sensitivity, sensorineural deafness, and AA amyloidosis. Arthritis Rheum. 2002;46: 2445–2452. pmid:12355493
  26. 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. pmid:27548431
  27. 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. pmid:14630794
  28. 28. WTCC. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007;447: 661–678. pmid:17554300
  29. 29. Castiblanco J, Arcos-Burgos M, Anaya J-M. What is next after the genes for autoimmunity? BMC Med. 2013;11: 197. pmid:24107170
  30. 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. 31. Baughman RP, Lower EE. Medical therapy of sarcoidosis. Semin Respir Crit Care Med. 2014;35: 391–406. pmid:25007090
  32. 32. Ciccarelli F, De Martinis M, Ginaldi L. An update on autoinflammatory diseases. Curr Med Chem. Bentham Science Publishers; 2014;21: 261–9. pmid:24164192
  33. 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. pmid:24635166
  34. 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. pmid:27200087
  35. 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. 36. Gautier L, Cope L, Bolstad BM, Irizarry RA. affy—analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004;20: 307–15. pmid:14960456
  37. 37. Carvalho BS, Irizarry RA. A framework for oligonucleotide microarray preprocessing. Bioinformatics. 2010;26: 2363–2367. pmid:20688976
  38. 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. pmid:27200087
  39. 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. 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. pmid:26834984
  41. 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.
  42. 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. pmid:21794127
  43. 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. pmid:26063839
  44. 44. Csardi G, Nepusz T. The igraph software package for complex network research. InterJournal. 2006;Complex Sy: 1695. citeulike-article-id:3443126
  45. 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. pmid:15980575
  46. 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. pmid:23703215
  47. 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. pmid:21789182
  48. 48. Arakelyan A, Nersisyan L, Poghosyan D, Khondkaryan L, Löffler-Wirth H, Melanitou E, et al. PSF-SOM full results [Dataset]. Zenodo. 2017;
  49. 49. de Seze J, Devos D, Castelnovo G, Labauge P, Dubucquoi 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. pmid:11673571
  50. 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. pmid:10719822
  51. 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. pmid:23776474
  52. 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. pmid:2590806
  53. 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. 54. Todd JA, Wicker LS. Genetic Protection from the Inflammatory Disease Type 1 Diabetes in Humans and Animal Models. Immunity. 2001;15: 387–395. pmid:11567629
  55. 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. 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. pmid:24171669
  57. 57. Zapatero DC, Colin IM. Behcet disease and type 1 diabetes, a fortuitous association? Louv Med. 2008;127.
  58. 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. pmid:22161697
  59. 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. pmid:23293022
  60. 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. pmid:20506269
  61. 61. McGonagle D, McDermott MF. A proposed classification of the immunological diseases. PLoS Med. Public Library of Science; 2006;3: e297. pmid:16942393
  62. 62. Lisnevskaia L, Murphy G, Isenberg D. Systemic lupus erythematosus. Lancet (London, England). 2014;384: 1878–88.
  63. 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. pmid:27821552
  64. 64. Kahlenberg JM, Kaplan MJ. The inflammasome and lupus. Curr Opin Rheumatol. 2014;26: 475–481. pmid:24992143
  65. 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. pmid:23315075
  66. 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. pmid:27250627
  67. 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. pmid:24449582
  68. 68. Hawkins PT, Stephens LR. PI3K signalling in inflammation. Biochim Biophys Acta. 2015;1851: 882–97. pmid:25514767
  69. 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. pmid:24718556
  70. 70. Covarrubias AJ, Aksoylar HI, Horng T. Control of macrophage metabolism and activation by mTOR and Akt signaling. Semin Immunol. 2015;27: 286–296. pmid:26360589
  71. 71. Bashaw GJ, Klein R. Signaling from axon guidance receptors. Cold Spring Harb Perspect Biol. Cold Spring Harbor Laboratory Press; 2010;2: a001941. pmid:20452961
  72. 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. pmid:25946140
  73. 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 co-infection. Virology. 2014;450–451: 336–349. pmid:24503097
  74. 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. pmid:20096956
  75. 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. pmid:9927410
  76. 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. pmid:17938642
  77. 77. Fitzpatrick M, Young S. Metabolomics–a novel window into inflammatory disease. Swiss Med Wkly. 2013; pmid:23348753
  78. 78. Khor B, Gardet A, Xavier RJ. Genetics and pathogenesis of inflammatory bowel disease. Nature. NIH Public Access; 2011;474: 307–17. pmid:21677747
  79. 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. pmid:26379192
  80. 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. pmid:24858699
  81. 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. pmid:11971877
  82. 82. Mulders-Manders CM, Simon A. Hyper-IgD syndrome/mevalonate kinase deficiency: what is new? Semin Immunopathol. Springer; 2015;37: 371–6. pmid:25990874
  83. 83. Wilson SP, Cassel SL. Inflammasome-mediated autoinflammatory disorders. Postgrad Med. NIH Public Access; 2010;122: 125–33. pmid:20861596
  84. 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. pmid:24919149
  85. 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. pmid:27270401
  86. 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. pmid:10369261
  87. 87. Favier LA, Schulert GS. Mevalonate kinase deficiency: current perspectives. Appl Clin Genet. Dove Press; 2016;9: 101–10. pmid:27499643
  88. 88. Liu S, Trapnell C. Single-cell transcriptome sequencing: recent advances and remaining challenges. F1000Research. Faculty of 1000 Ltd; 2016;5. pmid:26949524
  89. 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. pmid:16503242
  90. 90. Tektonidou MG, Ward MM. Validation of new biomarkers in systemic autoimmune diseases. Nat Rev Rheumatol. 2011;7: 708–717. pmid:22045310
  91. 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. pmid:22647780
  92. 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. pmid:23031813
  93. 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. pmid:22427800
  94. 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. 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. pmid:17805542
  96. 96. Eltzschig HK, Carmeliet P. Hypoxia and inflammation. N Engl J Med. 2011;364: 656–65. pmid:21323543
  97. 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. pmid:18755889
  98. 98. Wessel A, Hsu A, Goldbach-Mansky R, Zilberman-Rudenko J, Siegel R, Hanson E. Inflammatory disease due to selectively impaired NF-kB 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. 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. pmid:23434408
  100. 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. pmid:23434408
  101. 101. Lawrence T. The nuclear factor NF-kappaB pathway in inflammation. Cold Spring Harb Perspect Biol. Cold Spring Harbor Laboratory Press; 2009;1: a001651. pmid:20457564
  102. 102. Singh D, Naik S. The role of toll-like receptors in autoimmune diseases. J Indian Rheuatol Assoc. 2005;13: 162–165.
  103. 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. pmid:23162288
  104. 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. pmid:11477402
  105. 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. pmid:18990074
  106. 106. Martinon F, Aksentijevich I. New players driving inflammation in monogenic autoinflammatory diseases. Nat Rev Rheumatol. 2015;11: 11–20. pmid:25247411
  107. 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. pmid:24391216
  108. 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. pmid:21533026
  109. 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. pmid:16704732
  110. 110. Nimmerjahn F, Ravetch J V. Fcγ receptors as regulators of immune responses. Nat Rev Immunol. Nature Publishing Group; 2008;8: 34–47. pmid:18064051
  111. 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. pmid:10022447
  112. 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