- Open Access
Transcriptome analysis of normal-appearing white matter reveals cortisol- and disease-associated gene expression profiles in multiple sclerosis
Acta Neuropathologica Communications volume 7, Article number: 60 (2019)
Inter-individual differences in cortisol production by the hypothalamus–pituitary–adrenal (HPA) axis are thought to contribute to clinical and pathological heterogeneity of multiple sclerosis (MS). At the same time, accumulating evidence indicates that MS pathogenesis may originate in the normal-appearing white matter (NAWM). Therefore, we performed a genome-wide transcriptional analysis, by Agilent microarray, of post-mortem NAWM of 9 control subjects and 18 MS patients to investigate to what extent gene expression reflects disease heterogeneity and HPA-axis activity. Activity of the HPA axis was determined by cortisol levels in cerebrospinal fluid and by numbers of corticotropin-releasing neurons in the hypothalamus, while duration of MS and time to EDSS6 served as indicator of disease severity. Applying weighted gene co-expression network analysis led to the identification of a range of gene modules with highly similar co-expression patterns that strongly correlated with various indicators of HPA-axis activity and/or severity of MS. Interestingly, molecular profiles associated with relatively mild MS and high HPA-axis activity were characterized by increased expression of genes that actively regulate inflammation and by molecules involved in myelination, anti-oxidative mechanism, and neuroprotection. Additionally, group-wise comparisons of gene expression in white matter from control subjects and NAWM from (subpopulations of) MS patients uncovered disease-associated gene expression as well as strongly up- or downregulated genes in patients with relatively benign MS and/or high HPA-axis activity, with many differentially expressed genes being previously undescribed in the context of MS. Overall, the data suggest that HPA-axis activity strongly impacts on molecular mechanisms in NAWM of MS patients, but partly also independently of disease severity.
Multiple sclerosis (MS) is a chronic neurodegenerative disease that is hallmarked by the presence of inflammatory demyelinating lesions in the central nervous system (CNS). Therefore, major efforts have been dedicated to identify the molecular mechanisms at play during formation and expansion of MS lesions [45, 51, 53]. Related to this, there is an increasing interest for normal-appearing white matter (NAWM) as starting point of MS pathogenesis, as NAWM in MS was found to contain numerous cellular and molecular changes in inflammatory and neuroprotective pathways [60, 82, 89]. Studying these pathways may lead to the identification of events that precede or determine permissiveness for development of new MS lesions .
Clinical and pathological heterogeneity is one of the most striking features of MS [5, 53, 54, 84]. A factor considered to be associated with heterogeneity of MS is activity of the hypothalamus–pituitary–adrenal (HPA) axis [25, 26, 32, 58, 75]. Recently, we found in a post-mortem study of 42 MS patients that high levels of cortisol in the cerebrospinal fluid (CSF) were associated with slow disease progression, in particular in females with secondary-progressive MS, whereas patients with low cortisol levels had fast disease progression, greater numbers of active lesions, and less remyelinated plaques . Moreover, NAWM of female secondary-progressive MS patients with high cortisol levels displayed elevated expression of glucocorticoid-responsive genes, such as CD163, and decreased expression of pro-inflammatory genes, such as tumor necrosis factor, when compared to NAWM of patients with low cortisol levels. These data strongly suggest that high HPA-axis activity, by means of cortisol hypersecretion, impacts on molecular mechanisms in the NAWM in such a way that it reduces permissiveness for lesion development. Defining the molecular pathways induced in NAWM by cortisol may therefore identify targets for the development of therapeutic strategies to slow down or halt MS progression.
The heterogeneity of MS has been established at the clinical, radiological, and neuropathological level, but not with regard to molecular mechanisms in NAWM [5, 53, 54]. Without addressing disease heterogeneity, an earlier study on genome-wide expression in MS NAWM showed an upregulation of genes involved in anti-inflammatory mechanisms in oligodendrocytes and increased mRNA levels of pro-inflammatory molecules . Thus, no study thus far addressed the question whether and to what extent gene expression profiles in NAWM of MS patients reflect disease severity and/or HPA-axis activity. Here, we applied an Agilent microarray to perform a genome-wide transcriptional analysis of post-mortem NAWM of 18 female subjects with secondary-progressive MS that displayed a strong heterogeneity in disease severity and HPA-axis activity. This allowed us to investigate gene expression in relation to disease severity and/or HPA-axis activity. No male subjects were included, as an earlier study by us indicate that the inverse correlation of HPA axis activity with both disease severity and pathology of MS was most prominent in females with secondary progressive MS . As an added benefit, this provided us with a more homogeneous population, excluding possible confounding effects of gender on our findings.
We analyzed our data using various approaches, including WGCNA and group-wise comparisons of control subjects and (subpopulations of) MS patients. The WGCNA enabled us to identify clusters of co-regulated genes that correlate with one or more indicators of MS-disease severity and HPA-axis activity. In addition, individual genes were also studied for their association with the same clinical and endocrinological parameters. This led us to uncover many novel genes positively or negatively associated with HPA-axis activity and/or severity of MS. In general, gene expression profiles associated with high cortisol production and mild MS were characterized by molecules that actively regulate inflammation, but also belong to pathways involved in proliferation of neural stem cells. The data indicate that HPA-axis activity strongly impacts on molecular mechanisms in NAWM of MS patients and strongly affects pathways associated with disease severity. In this way, our study identifies molecular targets that may be assessed for their potential to prevent MS pathology in NAWM or limit lesion formation.
Materials and methods
Snap-frozen tissue from control and MS-brain donors was obtained from the Netherlands Brain Bank, Amsterdam, The Netherlands (http://www.brainbank.nl). Informed consent was obtained for brain autopsy and the use of tissue and clinical information for research purposes. Exclusion criteria were death due to sepsis and glucocorticoid treatment within 8 weeks prior to death. Clinical diagnoses of MS were confirmed by a neurologist (Prof. C.H. Polman, VUmc, Amsterdam or Dr. S. Luchetti, NIN, Amsterdam, The Netherlands). An overview of the pathological and clinical data of the brain donors is depicted in Table 1.
Quantification of CRH-producing neurons and cortisol levels
Numbers of corticotropin-releasing hormone (CRH)-expressing neurons in the paraventricular nucleus (PVN) were quantified in fixed tissue as described previously [35, 58]. In short, serial 6-μm frontal sections were cut on a microtome. Delineation of the PVN was determined in thionine-stained sections. Each 100th section through the PVN was stained for CRH. Neurons that showed a nucleolus and expressed CRH were counted blinded. The total number of CRH-expressing neurons in the PVN was calculated on the basis of cell counts and the distance between the sections. Cortisol was measured by radioimmunoassay (Diagnostic Products Corporation, Los Angeles, CA, USA), using a radioactively labeled antibody that enables highly sensitive detection of cortisol levels in various types of fluid, including serum, CSF, and saliva.
Tissue processing and RNA isolation
Series of 10 cryostat sections (20 μm each) of subcortical NAWM were homogenized in Trizol (Invitrogen, Carlsbad, CA, USA). Sections preceding and following these series were stained by immunohistochemistry for proteolipid protein (PLP; Serotec, Oxford, UK) and HLA-DP, −DQ, −DR (DakoCytomation, Glostrup, Denmark) to confirm the absence of MS-lesion pathology, respectively, by ruling out demyelination and microglia/macrophage activation. RNA isolation and assessment of its quality by RNA integrity number (RIN) was performed as described previously . Along with this, RNA was extracted from snap-frozen tissue dissected from various anatomical regions of control and MS brains, including MS lesions and NAWM, as well from tonsil. This RNA which was pooled to create common reference complementary RNA (cRNA), which was co-hybridized to every microarray slide to enable accurate comparison of expression levels across different cDNA microarray experiments.
Labeling of isolated RNA was done using the Low Input Quick Amp Labeling kit (Agilent Technologies, Palo Alto, CA, USA), according to the manufacturer’s instructions. For whole-genome expression analysis, samples were hybridized to Agilent 4x44K v2 Whole Human Genome arrays (G4845A; Agilent Technologies), covering 27,958 genes. In brief, equal amounts of total RNA (50 ng) were amplified and labeled with either Cy3-CTP (experimental samples) or Cy5-CTP (reference material, obtained as described above) using the Low Input Quick Amp Labeling kit (Agilent Technologies). For hybridization, equal amounts (825 ng) of labeled samples were fragmented in Fragmentation Buffer (Agilent Technologies) for 30 min at 60 °C. Labeled and fragmented complementary RNA (cRNA) was hybridized to the array and incubated in a rotating hybridization chamber for 17 h at 60 °C. After hybridization, the array was washed subsequently for 5 min in 6 x saline sodium phosphate-EDTA (SSPE)/0.005% N-lauroylsarcosine, 1 min in 0.006 x SSPE/0.005% N-lauroylsarcosine, and 30 s in acetonitrile and dried quickly in nitrogen flow. The arrays were scanned at a resolution of 5 mm and at 5 and 100% photomultiplier tube settings using the Agilent DNA Microarray Scanner (Agilent Technologies). Scan data were extracted using Agilent Feature Extraction software (version 8.5.1; Agilent Technologies).
Normalization of gene expression and gene extraction
To allow for accurate comparison of expression levels across different cDNA microarray experiments, common reference cRNA was co-hybridized to every microarray. By using the reference cRNA, a ratio between the experimental and reference material could be calculated for every spot, and expression levels across different hybridizations were compared. These data were used for analysis by weighted gene co-expression network analysis (WGCNA) and by group-wise comparison between subgroups, which is further described below. For an overview of the experimental setup, see Fig. 1. Raw data from the extraction software was imported to the R statistical processing environment using linear models for microarray data (LIMMA)/Bioconductor (version 3.12.3) package. Feature and background non-uniformity outliers, as determined by the Agilent Feature Extraction software using default settings, were removed. The data were normalized using “between array normalization” with the “Gquantile” algorithm (http://www.bioconductor.org).
cDNA synthesis and quantitative real-time PCR
Using the same RNA from which cRNA was generated for hybridization, cDNA was synthesized with the Quantitect Reverse Transcription Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. Quantitative real-time polymerase chain reaction (PCR) was performed and analyzed as described elsewhere, with minor adaptations . The amount of cDNA used per reaction was based on an input of 5 ng original RNA in a final volume of 20 μl.
To obtain lists of genes associated with either disease severity or HPA-axis activity, a WGCNA was performed on the 2log intensities to find gene modules [44, 61]. WCGNA defines, in an unbiased way, clusters of genes that show a co-expression pattern amongst samples, creating (a network with) modules of genes that are co-regulated. A soft power threshold of 9 was determined in the recommended way , a signed network was used, and dendrogram cluster detection was performed using a minimum cluster size of 27 and a deepSplit value of 1.
In addition to the WCGNA, linear regression analysis was performed for group-wise comparisons, using the LIMMA/bioconductor package in R (version 3.12.3). This was carried out using the limma function “lmscFit” for single channel analysis on two-color arrays with correction for intraspot correlation. A false discovery rate of 0.05 was used by multiple testing correction with the Benjamini-Hochberg method. All models were pairwise comparisons of two conditions. These group-wise comparisons were made for gene expression between all MS patients and control subjects. Furthermore, the group of MS patients was separated into a high and a low cortisol subgroup by the median for this parameter for objective group-wise comparison to identify cortisol-associated gene expression. To perform unbiased group-wise comparisons to find gene expression profiles associated with severity of MS, the patient population was divided into two subgroups by the median disease duration. Correlation coefficients were interpreted as follows: r < 0.1: no correlation; 0.1 < r < 0.3: weak correlation; 0.3 < r < 0.5: moderate correlation; r > 0.5: strong correlation. The significance level was set to 0.05. For all analyses that did not include the microarray data, we used SPSS software version 24 (IBM, Armonk, NY, USA).
All modules identified by WGCNA were subjected to the online tool DAVID, version 6.7 (http://david.abcc.ncifcrf.gov) for functional gene ontology (GO) analysis, using the whole human genome as background. For enrichment of top GO classes, the significance level was set to 0.05 (unadjusted). The minimum number of genes in each GO term was 3.
Post-mortem NAWM of 18 female subjects with secondary-progressive MS that displayed a strong heterogeneity in disease severity and HPA-axis activity and white matter of 10 matched control subjects was analyzed for differences in gene expression using Agilent Human Gene Expression 4x44K v2 microarrays. Table 2 provides an overview of the brain donor characteristics, including age, post-mortem delay (PMD), pH value of the CSF, and quality of RNA from NAWM in control subjects and different subgroups of MS patients. This table also displays clinical characteristics of the MS patients, such as disease duration and cortisol levels in CSF. No differences were present between control subjects and MS patients in age, PMD, or RIN. However, pH was significantly lower in MS patients compared to controls (p = 0.001). At the same time, pH and PMD did not correlate with RIN, cortisol levels, or numbers of CRH-expressing neurons, indicating that they had no confounding effects. There was a significant correlation between cortisol and numbers of CRH-positive neurons in MS patients (r = 0.508, p = 0.031), corroborating the idea that indicators of HPA-axis activity measured after death reliably reflect the situation during life . In addition, cortisol levels in CSF correlated with duration of MS (r = 0.490, p = 0.039) and with a trend towards significance with time to a score of 6 on the expanded disability status scale (EDSS6) (r = 0.432, p = 0.095), confirming earlier findings that high HPA-axis activity coincides with slower disease progression and relatively mild MS [35, 58].
Weighted gene co-expression network analysis (WGCNA)
A first analysis on the microarray data was done by WGCNA, which in an unbiased way identifies clusters of genes that display a similar expression pattern across the different samples in the whole data set; clusters of genes that show the same expression pattern are referred to as modules. WGCNA can also be used to investigate whether such modules are related to traits of particular interest. Each module has an eigengene (the first principal component) that represents the average expression pattern of all genes within the module. Moreover, for each gene within a module, a correlation value to the module eigengene is calculated, which indicates to what extent the expression profile of an individual gene resembles that of the whole module [44, 61].
Two approaches were taken for WGCNA. The main analysis was performed on a dataset that included the samples from both control subjects and MS patients. Another analysis was performed on data from MS patients only, which was done to confirm findings made in the first approach and to identify additional gene expression profiles associated with HPA-axis activity and disease severity in MS. From both datasets, we selected out of a total of 27,958 genes covered by the microarray the 15,000 genes that showed the biggest variance across all samples for application of WGCNA. By hierarchical clustering using these data we found two outliers, one in the control group and one in the MS patient group (donor 01–069 and 98–158 respectively), which were excluded from further analyses. All data described below in the text, tables, and figures were generated by our main WGCNA analysis that included both control subjects and MS patients, unless stated otherwise. To find gene clusters associated with disease severity and HPA-axis activity in MS, we studied the correlation between gene modules identified by WCGNA and several clinical traits: disease duration and time to EDSS6 as parameters of disease severity, CSF cortisol levels and numbers of CRH-producing neurons as indicators of HPA-axis activity.
Modules generated by WCGNA correlate to HPA-axis activity and MS severity
Using the WGCNA method, we identified a total of 19 modules, of which the biggest module contained 2495 genes and the smallest 65 genes (Fig. 2). Of the identified modules, 9 showed an expression pattern that correlated with one or more of the traits of interest: cortisol levels in CSF, number of CRH-positive neurons, disease duration, and time to EDSS6. Notably, all modules that were correlated to CSF cortisol levels also showed a correlation in the same direction to numbers of CRH-expressing neurons, substantiating their association with activity of the HPA axis. In total, eight modules correlated only to indicators of HPA-axis activity, i.e. cortisol levels in CSF and numbers of CRH-expressing neurons in the PVN: the black module (r = − 0.47, p = 0.020; r = − 0.31, p = 0.100, respectively), yellow module (r = − 0.63, p < 0.001; r = − 0.57, p = 0.003, respectively), red module (r = 0.43, p = 0.030; r = 0.46, p = 0.020, respectively), salmon module (r = 0.40, p = 0.050; r = 0.52, p = 0.008, respectively), cyan module (r = 0.48, p = 0.010; r = 0.65, p < 0.001, respectively), pink module (r = 0.56, p = 0.003; r = 0.77, p < 0.001, respectively), tan module (r = − 0.63, p < 0.001; r = − 0.61, p = 0.001, respectively), and turquoise module (r = 0.25, p = 0.200; r = 0.46, p = 0.020, respectively).
Four modules stood out for their particular strong correlation to one or more traits. The first one is the lightgreen module, which positively correlated to cortisol levels, disease duration, and time to EDSS6. Secondly, the pink module showed a strong positive correlation to cortisol levels and to numbers of CRH-expressing neurons. In contrast, the tan and yellow module correlated negatively to cortisol levels and numbers of CRH neurons. The 10 genes most strongly connected to the expression pattern of the four modules (module eigengene) are depicted in Table 3. The GO classes enriched in these four modules according to functional annotation clustering analysis are depicted in Table 4. The results of the GO analysis of all the other modules are shown in Additional file 1.
A module strongly correlated to both severity of MS and HPA-axis activity
Among the gene clusters identified by the WCGNA, as depicted in Fig. 2, the lightgreen module had an eigengene that was positively correlated to cortisol levels (r = 0.56, p = 0.003), disease duration (r = 0.50, p = 0.010), and time to EDSS6 (r = 0.54, p = 0.006), indicating its positive association with both HPA-axis activity and disease severity. The module contained a total of 98 genes, encoding molecules, such as the S100 proteins S100A8, S100A9, and S100A12, which play a variety of roles in inflammation [3, 13, 21]. Functional annotation clustering indicated that the lightgreen module was enriched for genes involved in classes labeled by GO as ‘immune response’ (p = 6.0E-09), ‘negative regulation of cytokine biosynthetic process’ (p = 4.3E-03), and ‘regulation of IL-1β production’ (p = 4.3E-03), amongst others (Table 4). These GO classes contained several genes that are directly or indirectly implicated in MS, such as cystatin F (CST7), ghrelin (GHRL), IL-18 receptor accessory protein (IL18RAP), myeloperoxidase (MPO), matrix-metalloproteinase 9 (MMP9), and FMS-like tyrosine-3 (FLT3) [4, 11, 17, 24, 27, 29, 38, 43, 55, 70, 79]. The top 10 genes most strongly connected to the lightgreen module are involved in various molecular mechanisms, such as myelination (CST7), neurogenesis (ASPM), and inflammation (DEFA4, VSTM1, BPI, S100P, WISP3, NLRP12). Overall, the lightgreen module was characterized by molecules that actively regulate inflammation, though it also indicated changes in molecular pathways that affect multiple aspects of MS disease activity to slow down clinical progression.
A module positively correlated to HPA-axis activity, independent of disease severity
The pink module was positively correlated to cortisol levels (r = 0.56, p = 0.003) and numbers of CRH neurons (r = 0.77, p < 0.001). The module contained 310 genes and was enriched for molecules with the GO label ‘regulation of caspase activity’ (p = 9.4E-04), ‘heat shock protein binding’ (p = 8.0E-04), ‘regulation of lymphocyte activation’ (p = 2.0E-02), and ‘oxidation reduction’ (p = 1.1E-02), amongst others (Table 4). Moreover, the module was enriched for molecules involved in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway ‘antigen processing and binding’ (p = 1.7E-03). The molecules involved in regulation of caspase activity may be related to apoptosis of neurons and oligodendrocytes, but also play a role in microglia activation and neurotoxicity .
Interestingly, the pink module contained the interleukin-7 receptor (IL7R, CD127) and tumor necrosis factor ligand superfamily member 14 (TNFSF14), two molecules that are highly associated with genetic susceptibility to MS [28, 36]. In addition, the pink module contained various heat-shock proteins (HSP), which have been associated with MS pathogenesis, including HSP90, which functions as a chaperone protein for the glucocorticoid receptor (GR) [14, 66]. Also, sirtuin 7 (SRT7) was present in the GO class ‘cellular response to oxidative stress’, which has been described as a protective molecule in the brain and has been proposed as a therapeutic agent .
Similar to the observations regarding the lightgreen module, the pink module is also characterized by a combination of genes involved in active regulation of immunity and those involved in more neuroprotective mechanisms.
Modules negatively correlated to HPA-axis activity, independent of disease severity
The tan module had a module eigengene that was negatively correlated to cortisol levels (r = − 0.63, p = 7.0E-04) and numbers of CRH neurons (r = − 0.61, p = 1.0E-03). The module contained 253 genes and was enriched for molecules in the GO class ‘generation of precursor metabolites and energy’ (p = 7.6E-10) and ‘mitochondrion’ (p = 1.3E-10), which are likely related to changes in mitochondrial functioning and energy metabolism known to occur in MS (Table 4) .
The yellow module had a module eigengene that was strongly correlated negatively to cortisol levels (r = − 0.63, p = 7.0E-04) and numbers of CRH neurons (r = − 0.57, p = 3.0E-03). The module contained 922 genes and was enriched for molecules belonging to the GO class ‘lipid biosynthetic process’ (p = 1.7E-05) and ‘steroid metabolic process’ (p = 2.4E-02) (Table 4). One particularly interesting member of both enriched GO classes in the tan module was aconitase (ACO2). ACO2 is activated by iron and is involved in glutamate production , while both iron accumulation and glutamate excitotoxicity are implicated as pathogenic mechanisms in MS [30, 68].
Genes associated with severity of MS and HPA-axis activity, independent of WGCNA modules
The WGCNA also produces direct correlations between genes and traits of interest, independent of gene modules. These correlations can be used to select individual genes that may be strongly related to HPA-axis activity as well as to disease severity and may therefore offer potential in development of therapeutic strategies for MS. The 10 genes with the strongest positive or negative correlation with disease severity and cortisol are depicted in Table 5.
Several molecules among the genes positively correlated to cortisol levels are well-known for their prominent role in shaping adaptive immune responses. For example, retinoic acid receptor alpha (RXRA) showed the strongest positive correlation with cortisol levels. Interestingly, RXRA-mediated signaling enhances differentiation of functionally competent CD4+ regulatory T cells and potently inhibits the formation of CD4+ T helper 17 cells, the latter of which have been strongly implicated in various autoimmune pathologies, including MS [6, 15, 18, 19, 22, 39, 62, 78]. Indeed, retinoic acid and other vitamin-A derivatives were shown to be protective in animal models of autoimmune disease . Notably, IL7R was also present among the 10 genes most positively correlated to cortisol levels. In combination with its presence in the pink module, which showed the strongest positive correlation with cortisol levels in the CSF and numbers of CRH-expressing neurons on the PVN, this further confirms the positive association of IL7R expression with HPA-axis activity in MS. Another interesting members of the 10 genes most positively correlated to cortisol levels is gap junction alpha-4 protein (GJA4), as it is strongly expressed by endothelial cells and may affect leukocyte trafficking . Expression of ankyrin repeat domain-containing protein 16 (ANKRD16) showed the strongest negative correlation to levels of cortisol in CSF. However, not much is known about this protein. The gene most strongly correlated to disease duration is the transcription factor NF-E2 45 kDa subunit (NFE2), which has been shown to play a role in differentiation and maturation of erythroid cells . The strongest negative association to duration was seen for the gene coding for methionine sulphoxide reductase A (MSRA), which is involved in protection against oxidative stress .
Group-wise comparisons using linear regression (LIMMA)
The median cortisol level (220.5 nmol/l) was used to objectively divide the MS-patient population into two subgroups, allowing comparison of groups for gene expression profiles that might reveal associations with HPA-axis activity. Similarly, the median duration of MS was used to create an objective subdivision between patients with relatively severe or mild disease (i.e., having a disease course shorter or longer than 24 years, respectively), to compare these groups for identification of differences in gene expression profiles related to disease severity. Two patients in the study population had a disease duration of exactly 24 years. These were included in the group of patients with severe MS, as their disease duration was shorter than the median disease duration of 24.5 years in the entire MS-donor population of the Netherlands Brain Bank. The results of the different comparisons are shown in Fig. 3 and Additional file 2, with the 10 most strongly regulated genes in either direction highlighted. The GO classes enriched in the comparisons according to functional annotation clustering analysis are depicted in Table 6.
Comparison of MS patients and controls
Compared to control subjects, a total of 778 genes was significantly upregulated in MS patients, whereas 544 genes were downregulated. The four most strongly upregulated genes in NAWM of MS patients code for antibody subunits, which is possibly related to the synthesis of auto-antibodies in MS (Fig. 3a) . Furthermore, NAWM of MS patients showed an increased expression of heat-shock proteins, HSPA1A and HSPA6. Among the genes upregulated in NAWM of MS patients compared to that of control subjects, there was an enrichment for molecules in GO classes associated with induction of apoptosis, activation of caspases, regulation of T-cell activation, and cytokines binding (Table 6).
Compared to control NAWM, the gene showing the most strongly decreased expression in MS NAWM was transthyretin (TTR), which was downregulated 4-fold (Fig. 3a). TTR is an important carrier in serum and CSF for the thyroid hormone thyroxin (T4) and for retinol, a form of vitamin A. Importantly, both TTR and T4 have been implicated in MS. Oxidative modifications of TTR protein and decreased levels of T4 were present in the CSF, and not in the serum, of MS patients and were correlated with disease duration . Moreover, T4 was shown to play an important role in activating oligodendrocyte precursor activation and instill myelination .
GO classes overrepresented among the genes downregulated in NAWM of MS relative to control NAWM were related to neuron differentiation, cell projection, and regulation of lipid metabolism (Table 6).
Comparison of MS patients with high and low cortisol
In MS patients with high cortisol a total of 270 genes was upregulated, whereas 472 genes were downregulated compared to patients with low cortisol. Of note, the third most-strongly upregulated gene in MS patients with high cortisol by more than 4-fold was CD163 (Fig. 3b), which is a glucocorticoid-responsive gene that can be induced in myeloid immune cells, such as macrophages and microglia [59, 80].
Neuronal pentraxin-2 (NPTX2) was the gene most strongly increased in NAWM of MS patients with high cortisol, showing an almost 6-fold higher expression in comparison to NAWM of patients with low cortisol (Fig. 3b). NPTX2 is an immune-related molecule with structural similarities to several acute phase proteins that is thought to be essential in the compensatory synaptic response that occurs during prolonged neuronal inactivity . An interesting finding regarding the 10 most strongly downregulated genes in NAWM of MS patients with high cortisol was the presence of the purinergic receptor P2RY12 (Fig. 3b), a microglia signature gene shown to play a major role in microglial activation, synaptic plasticity, and closure of the injured blood-brain barrier [8, 52, 76].
Interestingly, analysis by functional annotation clustering on the genes that are upregulated in MS patients with high cortisol levels, in comparison to patients with low cortisol levels, are enriched for several GO classes associated with (regulation of) inflammation (Table 6). In contrast, there is an overrepresentation of GO classes associated with apoptosis and negative regulation of neurogenesis among the genes showing decreased expression in MS patients with high cortisol.
Comparison of MS patients with severe and mild disease
In total 202 genes were upregulated in patients with mild MS, whereas 154 genes were downregulated compared to those with severe MS. There is a clear overlap between genes upregulated in MS patients with mild disease and those with high cortisol, as for example S100A8 and solute carrier organic anion transporter family member 4A1 (SLCO4A1) are present in the top 10 most strongly upregulated genes in both groups. IL-1 receptor like 1 (IL1R1) is among the 10 most strongly upregulated genes in mild MS (Fig. 3c), and is known to be a receptor for IL-33 and has been shown to be involved in induction of Th2 responses during allergic inflammation . The oxidized low-density lipoprotein receptor 1 (OLR1) gene, also known as LOX-1, was the most downregulated gene in NAWM of MS patients mild disease (Fig. 3c). Interestingly, this gene was found to be strongly related to the extent of demyelination in white matter MS lesions . The upregulated genes in NAWM of MS patients with mild disease were enriched for GO classes involved in neuron development and differentiation as well as the inflammatory response.
Distinct expression profile of inflammatory and GC-related genes in patients with high cortisol or mild MS
For the next analysis, we selected genes included in GO classes related to either inflammation or glucocorticoid signaling, including: ‘regulation of acute inflammatory response’ (GO:0002673), ‘chronic inflammatory response’ (GO:0002544), ‘macrophage differentiation’ (GO:0030225), ‘microglial cell activation’ (GO:0001774), ‘glucocorticoid-receptor signaling pathway’ (GO:0042921), and ‘cellular response to glucocorticoids’ (GO:0071385) and ‘glucocorticoid biosynthetic process’ (GO:0006704). These GO classes were selected to specifically study to what extent HPA-axis activity impacts on molecular mechanisms in NAWM of MS and how this affects neuroinflammation and disease severity. To this end, the selected genes were compared for their expression between five groups: control subjects, MS patients with high and low cortisol levels, and MS patients with severe and mild disease. This was done by cluster analysis, to find out which groups most strongly resemble each other in expression profiles for the selected genes. Interestingly, a distinct expression profile was present for inflammatory and glucocorticoid-associated genes in MS patients with high cortisol and those with mild MS (Fig. 4a and b). These subgroups differed in this respect from MS patients with low cortisol or severe MS, and control subjects. Remarkably, the subgroups of MS patients with high cortisol and those with mild MS were especially similar with respect to higher average expression values for inflammatory genes, when compared to the other three subgroups.
Molecules selected for validation
Expression of in total 14 genes strongly associated with HPA-axis activity and/or disease duration was validated by qPCR. For this purpose, the same tissue was used as for generating the original microarray data. For several of the validated targets no alterations have been described in the context of MS before, such as leucine-rich repeat-containing protein 32 (LRRC32) and nidogen-1 (NID1). Others have been implicated in MS pathogenesis, but were never studied in human subjects for their association with disease severity or HPA-axis activity, such as IL7R, tissue transglutaminase (TGM2), MAC-inhibitory protein (CD59), and gap junction alpha-4 (GJA4) [10, 41, 77, 83]. For 11 out of 14 genes the expression profile was indeed confirmed by qPCR (Fig. 5). Genes for which the qPCR data did not confirm the expression profile detected by microarray were SLC4, CD44, and CST7. For CST7, the expression detected by the microarray as well as the qPCR was very low. Therefore, the expression differences for CST7 observed in the microarray data may have been a technical artifact caused by background variation. For SLC4 and CD44, quite some variation in expression between subgroups was detected by qPCR, but apparently those differences were not large enough to be significant.
This is the very first study, to our knowledge, that reports on genome-wide transcriptional changes in NAWM of MS in relation to disease severity as well as activity of the HPA axis. To this end, we purposefully included female MS brain donors that differed strongly in rate of disease progression and HPA-axis activity. This enabled us to define gene expression profiles associated with disease severity in NAWM that are cortisol-dependent and independent, thereby giving a novel and comprehensive insight into the molecular underpinnings of the clinical heterogeneity that is so characteristic for MS. Importantly, we identified a range of molecular changes associated with a more benign or aggressive disease course in MS that may be targeted for development of new therapeutic strategies. In addition, we uncovered disease-associated gene expression profiles by comparative analysis with white matter from control brain donors, yielding additional molecular targets that may be exploited therapeutically as well.
Data analysis was done using various approaches, including WGCNA and group-wise comparisons of control subjects and (subpopulations of) MS patients. By WGCNA, we were able to identify clusters of co-regulated genes that correlate with one or more indicators of MS-disease severity and HPA-axis activity. Irrespective of their module membership, genes were also studied individually for their association with the same clinical and endocrinological parameters. In this way, the data uncovered many novel genes positively or negatively associated with HPA-axis activity and/or severity of MS. In general, gene expression profiles associated with high cortisol production and mild MS were characterized by molecules that actively regulate inflammation, but also belong to pathways involved in proliferation of neural stem cells. Together, these data reveal that HPA-axis activity strongly impacts on molecular mechanisms/changes in NAWM of MS patients, but these changes are in part independent of the transcriptional changes associated with disease severity.
We show that gene expression profiles associated with high cortisol production and mild MS patients are characterized by molecules that negatively regulate inflammation and immunity, such as NLRP12, S1PR4, S100A8, S100A9, and S100A12 [1, 3, 13, 64, 74]. In this way, our study identifies various molecular targets that may be assessed for their potential to prevent MS pathology in NAWM or limit lesion formation.
The light green gene module is the most relevant module for unraveling the molecular mechanisms of cortisol-mediated suppression of MS-disease progression, as it was correlated to numbers of CRH neurons, cortisol levels, and duration of MS. Functional annotation clustering analysis pointed out that the lightgreen module contained several genes that may be associated with slower progression of MS by regulation of inflammation. Notable molecules present in the lightgreen cluster are S100A8 and S100A9, which are strongly expressed by myeloid immune cells, have well-established immunoregulatory properties and are implicated in protection against oxidative stress . The endogenous danger signal S100A9 plays a key role in immune escape of solid tumors, where its chronic expression in myeloid cells inhibits their maturation and thereby skews them to an immunosuppressive phenotype . Moreover, S100A8 and S100A9 serve as inflammatory biomarkers in several autoimmune disorders, such as systemic lupus erythematosus and inflammatory bowel disease [42, 46, 63]. The anti-inflammatory role of S100A8 and S100A9 are further indicated by the finding that glucocorticoids directly induce these proteins in human monocytes and dendritic cells, and that S100A8-positive macrophages are increased in synovial fluid after treatment of rheumatoid arthritis patients with high-dose methylprednisolone .
Another prominent immunoregulatory molecule in the lightgreen module was NLR family, pyrin domain-containing 12 (NLRP12), which was present in several enriched GO classes, such as ‘regulation of IL-1b production’ and ‘negative regulation of cytokine biosynthetic process’. Moreover, NLRP12 expression was highly correlated to cortisol levels. NLRP12 is strongly expressed by myeloid cells and has been found to suppress canonical and non-canonical NF-κB signalling [1, 2, 37, 50, 86]. As such, NLRP12 is able to inhibit Toll-like receptor-induced activation and chemokine production in monocytes and other myeloid cells  .
In our WGCNA analysis, we also assessed to what extent the expression pattern of a gene resembles that of the whole module by looking at the connectivity of genes to the module eigengene. Interestingly, GJA4 (gap-junction alpha-4) not only was among the 10 genes most strongly correlated to cortisol, but was also significantly upregulated in mild MS. The expression of GJA4 has been shown to protect against the formation of atherosclerotic plaques by decreased recruitment and local adhesion of monocytes, a mechanism that is thought to be crucially involved in formation of MS lesions . Therefore, high levels of cortisol in MS patients may increase expression of GJA4 on monocytes to diminished their recruitment into the CNS and thereby limit MS pathology.
By studying direct correlations between individual genes with cortisol levels and/or duration of MS, we found several genes that are implicated in (inhibition of) remyelination. Examples are ASPM (abnormal spindle-like microcephaly-associated protein) and TLR2 (Toll-like receptor 2) . Other genes may be more related to neuroprotection. One of these is NPTX2 (neuronal pentraxin-2), which showed a strong correlation to cortisol and was highly increased in MS patients with high cortisol compared to those with low cortisol. This gene is also known as neuronal activity-regulated pentraxin (NARP) and has been shown to be essential for long-term synaptic plasticity, in particular in formation and maintenance of excitatory synapses [69, 88]. As such, NPTX2 may be an important neuroprotective gene in MS that is induced by glucocorticoids.
The pink module showed the strongest positive correlation with HPA-axis activity, both for its correlation with cortisol levels in CSF and numbers of CRH-positive neurons in the PVN. The module was enriched for genes involved in several GO classes, such as ‘regulation of caspase activity’ and ‘heat shock protein binding’. In addition, also genes belonging to the GO class ‘regulation of lymphocyte activation’ were enriched. Among these, the presence of IL7R (IL-7 receptor) was most striking, as it shows a high allelic association with MS susceptibility [28, 41]. IL7R was also found to be induced by dexamethasone in human blood leukocytes, which may explain why its expression was strongly correlated to levels of CSF cortisol in MS patients . Ligation of IL7R by IL-7 was found to be required for autoimmune neuroinflammation in experimental autoimmune encephalomyelitis . Moreover, it was recently reported that downregulation of IL7R expression in oligodendrocytes contributes to CNS demyelination in zebrafish . Thus, IL7R may play opposing roles in MS pathogenesis depending on the cell type it is expressed on. This may also explain why no correlation between IL7R expression and disease severity was observed.
HSPA1A (heat shock protein family A member 1A, Hsp70) and SERPINA3 (serpin family A member 3, alpha1-antichymotrypsin) were higher expressed in NAWM of MS patients with high cortisol compared to those with low cortisol. Heat-shock proteins have been shown to play an important role in limiting T-cell mediated (chronic) inflammation [65, 81]. However, expression of these genes was increased in NAWM of all MS patients when compared to control subjects. This may indicate that upregulation of this pathway represents a general protective mechanism in NAWM of MS that is further enhanced under the influence of high HPA-axis activity.
Importantly, the expression profile of most genes selected for validation could be confirmed by qPCR. However, for the few genes selected for further analysis at the protein level, the signal generated by the antibodies used was not specific enough to enable us to draw clear conclusions about cellullar distribution and potential expression differences. Still, the qPCR data provide valuable biological evidence that observed expression differences in the microarray data are real. In this respect, LRRC32 is particularly notable gene, as qPCR analysis indicated clearly elevated levels of expression in MS patients with high cortisol or mild MS. LRRC32 is a key regulator of transforming growth factor beta (TGF-β), which is known to be involved maintaining the molecular and functional signature of microglia . Moreover, TGF-β is known to promote development of either CD4+ regulatory T cells . It is tempting to speculate that, in concert with increased IL7R expression in patients with high cortisol and mild MS, this may indicate the presence of mechanism that serves to maintain a population of T cells with a regulatory phenotype to limit autoimmunity. As the expression of LRRC32 was also very significantly increased in all MS patients compared to controls subjects (data not shown), this could be a general mechanism that is is at play in MS NAWM.
In summary, our data indicate that HPA-axis activity strongly impacts on molecular mechanisms in NAWM of MS patients and thereby has modulatory effects on MS-disease activity. This extends the findings we reported previously . When observing the molecular profiles associated with high HPA-axis activity and relatively mild MS, what is in general most striking is the enhanced expression of genes that actively regulate inflammation. At the same time, NAWM does not show signs of active inflammation. While this has also been indicated by previous studies, our data are the first to indicate more clearly how MS patients may benefit from approaches that promote specific physiological mechanisms for immunosuppression under conditions of chronic inflammation, such as induction of NLRP12, S100A8, and S100A9. At the same time, our data also indicate that high HPA-axis activity is associated with changes in molecular pathways that may affect multiple aspects of MS disease activity to slow down clinical progression. For example, there is a clear increase in the expression of genes involved in remyelination and pathways that likely exert a protective effect against for example oxidative stress, synaptic disintegration, and axonal damage. Most importantly, the study led us to uncover a range of molecular changes associated with a more benign or aggressive disease course in MS that may be targeted for development of new therapeutic strategies. This brings up the important question whether the current regimens for glucocorticoid treatment of MS patients actually promote the molecular mechanisms associated with slow disease progression detected in our study and, if not, how these regimens can be adjusted to ensure that they do. Additionally, it would be important to find out of those same molecular mechanisms associated with slow progression of MS in this study can be induced more effectively or directly without the use of glucocorticoids, especially in those MS patients that respond poorly to glucocorticoids treatment. By comparing MS patients to control subjects, we were able to define disease-associated gene expression profiles in NAWM, which revealed additional molecular targets that may be used for therapeutic exploitation. Since only females were included in this study, it would be important to verify to what extent our findings are also relevant for male patients. This strongly warrants in vitro and in vivo studies to validate which of the molecular targets identified here display the biggest therapeutic potential for the treatment of MS.
Allen IC, Wilson JE, Schneider M et al (2012) NLRP12 suppresses colon inflammation and tumorigenesis through the negative regulation of noncanonical NF-κB signaling. Immunity 36(5):742–754
Arthur JC, Lich JD, Aziz RK, Kotb M, Ting JP-Y (2007) Heat shock protein 90 associates with monarch-1 and regulates its ability to promote degradation of NF-kappaB-inducing kinase. J Immunol 179(9):6291–6296
Baillet A, Trocmé C, Berthier S et al (2010) Synovial fluid proteomic fingerprint: S100A8, S100A9 and S100A12 proteins discriminate rheumatoid arthritis from other inflammatory joint diseases. Rheumatology 49(4):671–682
Berilgen MS, Bulut S, Ustundag B, Tekatas A, Ayar A (2005) Patients with multiple sclerosis have higher levels of serum ghrelin. Neuro Endocrinol Lett 26(6):819–822
Bielekova B, Kadom N, Fisher E et al (2005) MRI as a marker for disease heterogeneity in multiple sclerosis. Neurology 65(7):1071–1076
Brucklacher-Waldert V, Stuerner K, Kolster M, Wolthausen J, Tolosa E (2009) Phenotypical and functional characterization of T helper 17 cells in multiple sclerosis. Brain 132(Pt 12):3329–3341
Burguillos MA, Deierborg T, Kavanagh E et al (2011) Caspase signalling controls microglia activation and neurotoxicity. Nature 472(7343):319–324
Butovsky O, Jedrychowski MP, Moore CS et al (2014) Identification of a unique TGF-β-dependent molecular and functional signature in microglia. Nat Neurosci 17(1):131–143
Calza L, Fernandez M, Giuliani A, Aloe L, Giardino L (2002) Thyroid hormone activates oligodendrocyte precursors and increases a myelin-forming protein and NGF content in the spinal cord during experimental allergic encephalomyelitis. Proc Natl Acad Sci U S A 99(5):3258–3263
Carrette F, Surh CD (2012) IL-7 signaling and CD127 receptor regulation in the control of T cell homeostasis. Semin Immunol 24(3):209–217
Chandler S, Coates R, Gearing A, Lury J, Wells G, Bone E (1995) Matrix metalloproteinases degrade myelin basic protein. Neurosci Lett 201(3):223–226
Chen W, Jin W, Hardegen N, Lei K-J, Li L, Marinos N, McGrady G, Wahl SM (2003) Conversion of peripheral CD4+CD25- naive T cells to CD4+CD25+ regulatory T cells by TGF-beta induction of transcription factor Foxp3. J Exp Med 198(12):1875–1886
Cheng P, Corzo CA, Luetteke N et al (2008) Inhibition of dendritic cell differentiation and accumulation of myeloid-derived suppressor cells in cancer is regulated by S100A9 protein. J Exp Med 205(10):2235–2249
Cid C, Alvarez-Cermeño JC, Camafeita E, Salinas M, Alcázar A (2004) Antibodies reactive to heat shock protein 90 induce oligodendrocyte precursor cell death in culture. Implications for demyelination in multiple sclerosis. FASEB J 18(2):409–411
Coombes JL, Siddiqui KRR, Arancibia-Cárcamo CV, Hall J, Sun C-M, Belkaid Y, Powrie F (2007) A functionally specialized population of mucosal CD103+ DCs induces Foxp3+ regulatory T cells via a TGF-beta and retinoic acid-dependent mechanism. J Exp Med 204(8):1757–1764
de Groot M, Verhaaren BFJ, de Boer R, Klein S, Hofman A, van der Lugt A, Ikram MA, Niessen WJ, Vernooij MW (2013) Changes in normal-appearing white matter precede development of white matter lesions. Stroke 44(4):1037–1042
DeBoy CA, Rus H, Tegla C, Cudrici C, Jones MV, Pardo CA, Small D, Whartenby KA, Calabresi PA (2010) FLT-3 expression and function on microglia in multiple sclerosis. Exp Mol Pathol 89(2):109–116
Denning TL, Wang Y-C, Patel SR, Williams IR, Pulendran B (2007) Lamina propria macrophages and dendritic cells differentially induce regulatory and interleukin 17-producing T cell responses. Nat Immunol 8(10):1086–1094
Elias KM, Laurence A, Davidson TS, Stephens G, Kanno Y, Shevach EM, O'Shea JJ (2008) Retinoic acid inhibits Th17 polarization and enhances FoxP3 expression through a Stat-3/Stat-5 independent signaling pathway. Blood 111(3):1013–1020
Elliott C, Lindner M, Arthur A et al (2012) Functional identification of pathogenic autoantibody responses in patients with multiple sclerosis. Brain 135(Pt 6):1819–1833
Foell D, Wittkowski H, Vogl T, Roth J (2007) S100 proteins expressed in phagocytes: a novel group of damage-associated molecular pattern molecules. J Leukoc Biol 81(1):28–37
Fujino S, Andoh A, Bamba S, Ogawa A, Hata K, Araki Y, Bamba T, Fujiyama Y (2003) Increased expression of interleukin 17 in inflammatory bowel disease. Gut 52(1):65–70
Galon J, Franchimont D, Hiroi N, Frey G, Boettner A, Ehrhart-Bornstein M, O'Shea JJ, Chrousos GP, Bornstein SR (2002) Gene profiling reveals unknown enhancing and suppressive actions of glucocorticoids on immune cells. FASEB J 16(1):61–71
Gillett A, Thessen Hedreul M, Khademi M, Espinosa A, Beyeen AD, Jagodic M, Kockum I, Harris RA, Olsson T (2010) Interleukin 18 receptor 1 expression distinguishes patients with multiple sclerosis. Mult Scler 16(9):1056–1065
Gold SM, Raji A, Huitinga I, Wiedemann K, Schulz K-H, Heesen C (2005) Hypothalamo-pituitary-adrenal axis activity predicts disease progression in multiple sclerosis. J Neuroimmunol 165(1–2):186–191
Grasser A, Möller A, Backmund H, Yassouridis A, Holsboer F (1996) Heterogeneity of hypothalamic-pituitary-adrenal system response to a combined dexamethasone-CRH test in multiple sclerosis. Exp Clin Endocrinol Diabetes 104(1):31–37
Gray E, Thomas TL, Betmouni S, Scolding N, Love S (2008) Elevated myeloperoxidase activity in white matter in multiple sclerosis. Neurosci Lett 444(2):195–198
Gregory SG, Schmidt S, Seth P et al (2007) Interleukin 7 receptor alpha chain (IL7R) shows allelic and functional association with multiple sclerosis. Nat Genet 39(9):1083–1091
Gutcher I, Urich E, Wolter K, Prinz M, Becher B (2006) Interleukin 18-independent engagement of interleukin 18 receptor-alpha is required for autoimmune inflammation. Nat Immunol 7(9):946–953
Hametner S, Wimmer I, Haider L, Pfeifenbring S, Brück W, Lassmann H (2013) Iron and neurodegeneration in the multiple sclerosis brain. Ann Neurol 74(6):848–861
Hanafy KA, Sloane JA (2011) Regulation of remyelination in multiple sclerosis. FEBS Lett 585(23):3821–3828
Heesen C, Gold SM, Raji A, Wiedemann K, Schulz K-H (2002) Cognitive impairment correlates with hypothalamo-pituitary-adrenal axis dysregulation in multiple sclerosis. Psychoneuroendocrinology 27(4):505–517
Hendrickx DAE, Koning N, Schuurman KG, van Strien ME, van Eden CG, Hamann J, Huitinga I (2013) Selective upregulation of scavenger receptors in and around demyelinating areas in multiple sclerosis. J Neuropathol Exp Neurol 72(2):106–118
Hsu K, Passey RJ, Endoh Y, Rahimi F, Youssef P, Yen T, Geczy CL (2005) Regulation of S100A8 by glucocorticoids. J Immunol 174(4):2318–2326
Huitinga I, Erkut ZA, van Beurden D, Swaab DF (2004) Impaired hypothalamus-pituitary-adrenal axis activity and more severe multiple sclerosis with hypothalamic lesions. Ann Neurol 55(1):37–45
International Multiple Sclerosis Genetics Consortium, Wellcome Trust Case Control Consortium 2, Sawcer S et al (2011) Genetic risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis. Nature 476(7359):214–219
Jéru I, Duquesnoy P, Fernandes-Alnemri T et al (2008) Mutations in NALP12 cause hereditary periodic fever syndromes. Proc Natl Acad Sci U S A 105(5):1614–1619
Jha S, Srivastava SY, Brickey WJ, Iocca H, Toews A, Morrison JP, Chen VS, Gris D, Matsushima GK, Ting JP-Y (2010) The inflammasome sensor, NLRP3, regulates CNS inflammation and demyelination via caspase-1 and interleukin-18. J Neurosci 30(47):15811–15820
Kehlen A, Pachnio A, Thiele K, Langner J (2003) Gene expression induced by interleukin-17 in fibroblast-like synoviocytes of patients with rheumatoid arthritis: upregulation of hyaluronan-binding protein TSG-6. Arthritis Res Ther 5(4):R186–R192
Koning N, Bö L, Hoek RM, Huitinga I (2007) Downregulation of macrophage inhibitory molecules in multiple sclerosis lesions. Ann Neurol 62(5):504–514
Kreft KL, Verbraak E, Wierenga-Wolf AF, van Meurs M, Oostra BA, Laman JD, Hintzen RQ (2012) The IL-7Rα pathway is quantitatively and functionally altered in CD8 T cells in multiple sclerosis. J Immunol 188(4):1874–1883
Kuruto R, Nozawa R, Takeishi K, Arai K, Yokota T, Takasaki Y (1990) Myeloid calcium binding proteins: expression in the differentiated HL-60 cells and detection in sera of patients with connective tissue diseases. J Biochem 108(4):650–653
La Russa A, Cittadella R, De Marco EV, Valentino P, Andreoli V, Trecroci F, Latorre V, Gambardella A, Quattrone A (2010) Single nucleotide polymorphism in the MMP-9 gene is associated with susceptibility to develop multiple sclerosis in an Italian case-control study. J Neuroimmunol 225(1–2):175–179
Langfelder P, Horvath S (2008) WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9(1):559
Lassmann H, Ransohoff RM (2004) The CD4-Th1 model for multiple sclerosis: a critical [correction of crucial] re-appraisal. Trends Immunol 25(3):132–137
Lawrance IC, Fiocchi C, Chakravarti S (2001) Ulcerative colitis and Crohn's disease: distinctive gene expression profiles and novel susceptibility candidate genes. Hum Mol Genet 10(5):445–456
Lawson BR, Gonzalez-Quintial R, Eleftheriadis T, Farrar MA, Miller SD, Sauer K, McGavern DB, Kono DH, Baccala R, Theofilopoulos AN (2015) Interleukin-7 is required for CD4(+) T cell activation and autoimmune neuroinflammation. Clin Immunol 161(2):260–269
Lei X, Cai S, Chen Y, Cui J, Wang Y, Li Z, Li Y (2017) Down-regulation of interleukin 7 receptor (IL-7R) contributes to central nervous system demyelination. Oncotarget 8(17):28395–28407
Li YJ, Higgins RR, Pak BJ, Shivdasani RA, Ney PA, Archer M, Ben-David Y (2001) p45(NFE2) is a negative regulator of erythroid proliferation which contributes to the progression of friend virus-induced erythroleukemias. Mol Cell Biol 21(1):73–80
Lich JD, Williams KL, Moore CB, Arthur JC, Davis BK, Taxman DJ, Ting JP-Y (2007) Monarch-1 suppresses non-canonical NF-kappaB activation and p52-dependent chemokine expression in monocytes. J Immunol 178(3):1256–1260
Lock C, Hermans G, Pedotti R et al (2002) Gene-microarray analysis of multiple sclerosis lesions yields new targets validated in autoimmune encephalomyelitis. Nat Med 8(5):500–508
Lou N, Takano T, Pei Y, Xavier AL, Goldman SA, Nedergaard M (2016) Purinergic receptor P2RY12-dependent microglial closure of the injured blood-brain barrier. Proc Natl Acad Sci U S A 113(4):1074–1079
Lucchinetti C, Brück W, Parisi J, Scheithauer B, Rodriguez M, Lassmann H (2000) Heterogeneity of multiple sclerosis lesions: implications for the pathogenesis of demyelination. Ann Neurol 47(6):707–717
Luchetti S, Fransen NL, van Eden CG, Ramaglia V, Mason M, Huitinga I (2018) Progressive multiple sclerosis patients show substantial lesion activity that correlates with clinical disease severity and sex: a retrospective autopsy cohort analysis. Acta Neuropathol 135(4) 511–528.
Ma J, Tanaka KF, Shimizu T, Bernard CCA, Kakita A, Takahashi H, Pfeiffer SE, Ikenaka K (2011) Microglial cystatin F expression is a sensitive indicator for ongoing demyelination with concurrent remyelination. J Neurosci Res 89(5):639–649
Massacesi L, Abbamondi AL, Giorgi C, Sarlo F, Lolli F, Amaducci L (1987) Suppression of experimental allergic encephalomyelitis by retinoic acid. J Neurol Sci 80(1):55–64
McGahan MC, Harned J, Mukunnemkeril M, Goralska M, Fleisher L, Ferrell JB (2005) Iron alters glutamate secretion by regulating cytosolic aconitase activity. Am J Phys Cell Physiol 288(5):C1117–C1124
Melief J, de Wit SJ, van Eden CG, Teunissen C, Hamann J, Uitdehaag BM, Swaab D, Huitinga I (2013) HPA axis activity in multiple sclerosis correlates with disease severity, lesion type and gene expression in normal-appearing white matter. Acta Neuropathol 126(2):237–249
Melief J, Koning N, Schuurman KG, van de Garde MDB, Smolders J, Hoek RM, van Eijk M, Hamann J, Huitinga I (2012) Phenotyping primary human microglia: tight regulation of LPS responsiveness. Glia 60(10):1506–1517
Melief J, Schuurman KG, van de Garde MDB, Smolders J, van Eijk M, Hamann J, Huitinga I (2013) Microglia in normal appearing white matter of multiple sclerosis are alerted but immunosuppressed. Glia 61(11):1848–1861
Miller JA, Horvath S, Geschwind DH (2010) Divergence of human and mouse brain transcriptome highlights Alzheimer disease pathways. Proc Natl Acad Sci U S A 107(28):12698–12703
Mucida D, Park Y, Kim G, Turovskaya O, Scott I, Kronenberg M, Cheroutre H (2007) Reciprocal TH17 and regulatory T cell differentiation mediated by retinoic acid. Science 317(5835):256–260
Nacken W, Roth J, Sorg C, Kerkhoff C (2003) S100A9/S100A8: myeloid representatives of the S100 protein family as prominent players in innate immunity. Microsc Res Tech 60(6):569–580
Olesch C, Ringel C, Brüne B, Weigert A (2017) Beyond immune cell migration: the emerging role of the Sphingosine-1-phosphate receptor S1PR4 as a modulator of innate immune cell activation. Mediat Inflamm 2017:1–12
Peferoen LAN, Gerritsen WH, Breur M, Ummenthum KMD, Peferoen-Baert RMB, van der Valk P, van Noort JM, Amor S (2015) Small heat shock proteins are induced during multiple sclerosis lesion development in white but not grey matter. Acta Neuropathol Commun 3(1):132
Picard D, Yamamoto KR (1987) Two signals mediate hormone-dependent nuclear localization of the glucocorticoid receptor. EMBO J 6(11):3333–3340
Pieragostino D, Del Boccio P, Di Ioia M et al (2013) Oxidative modifications of cerebral transthyretin are associated with multiple sclerosis. Proteomics 13(6):1002–1009
Pitt D, Werner P, Raine CS (2000) Glutamate excitotoxicity in a model of multiple sclerosis. Nat Med 6(1):67–70
Pribiag H, Stellwagen D (2014) Neuroimmune regulation of homeostatic synaptic plasticity. Neuropharmacology 78:13–22
Proost P, Vandamme J, Opdenakker G (1993) Leukocyte gelatinase B cleavage releases Encephalitogens from human myelin basic protein. Biochem Biophys Res Commun 192(3):1175–1181
Reed KE, Westphale EM, Larson DM, Wang HZ, Veenstra RD, Beyer EC (1993) Molecular cloning and functional expression of human connexin37, an endothelial cell gap junction protein. J Clin Invest 91(3):997–1004
Rice CM, Sun M, Kemp K, Gray E, Wilkins A, Scolding NJ (2012) Mitochondrial sirtuins--a new therapeutic target for repair and protection in multiple sclerosis. Eur J Neurosci 35(12):1887–1893
Savinko T, Karisola P, Lehtimäki S, Lappeteläinen A-M, Haapakoski R, Wolff H, Lauerma A, Alenius H (2013) ST2 regulates allergic airway inflammation and T-cell polarization in epicutaneously sensitized mice. J Invest Dermatol 133(11):2522–2529
Schulze T, Golfier S, Tabeling C, Räbel K, Gräler MH, Witzenrath M, Lipp M (2011) Sphingosine-1-phospate receptor 4 (S1P4) deficiency profoundly affects dendritic cell function and TH17-cell differentiation in a murine model. FASEB J 25(11):4024–4036
Schumann EM, Kümpfel T, Then Bergh F, Trenkwalder C, Holsboer F, Auer DP (2002) Activity of the hypothalamic-pituitary-adrenal axis in multiple sclerosis: correlations with gadolinium-enhancing lesions and ventricular volume. Ann Neurol 51(6):763–767
Sipe GO, Lowery RL, Tremblay ME, Kelly EA, Lamantia CE, Majewska AK (2016) Microglial P2Y12 is necessary for synaptic plasticity in mouse visual cortex. Nat Commun 7:10905
Song W-C (2006) Complement regulatory proteins and autoimmunity. Autoimmunity 39(5):403–410
Sun C-M, Hall JA, Blank RB, Bouladoux N, Oukka M, Mora JR, Belkaid Y (2007) Small intestine lamina propria dendritic cells promote de novo generation of Foxp3 T reg cells via retinoic acid. J Exp Med 204(8):1775–1785
Theil M-M, Miyake S, Mizuno M et al (2009) Suppression of experimental autoimmune encephalomyelitis by ghrelin. J Immunol 183(4):2859–2866
Van de Garde MDB, Martinez FO, Melgert BN, Hylkema MN, Jonkers RE, Hamann J (2014) Chronic exposure to fluticasone propionate shapes gene expression and modulates innate and adaptive activation pathways in human macrophages. J Immunol 192(3):1196–1208
van Eden W, van der Zee R, Prakken B (2005) Heat-shock proteins induce T-cell regulation of chronic inflammation. Nat Rev Immunol 5(4):318–330
van Horssen J, Singh S, van der Pol S et al (2012) Clusters of activated microglia in normal-appearing white matter show signs of innate immune activation. J Neuroinflammation 9:156
van Strien ME, de Vries HE, Chrobok NL et al (2015) Tissue transglutaminase contributes to experimental multiple sclerosis pathogenesis and clinical outcome by promoting macrophage migration. Brain Behav Immun 50:141–154
Weiner HL (2009) The challenge of multiple sclerosis: how do we cure a chronic heterogeneous disease? Ann Neurol 65(3):239–248
Weissbach H, Resnick L, Brot N (2005) Methionine sulfoxide reductases: history and cellular role in protecting against oxidative damage. Biochim Biophys Acta 1703(2):203–212
Williams KL, Lich JD, Duncan JA et al (2005) The CATERPILLER protein monarch-1 is an antagonist of toll-like receptor-, tumor necrosis factor alpha-, and Mycobacterium tuberculosis-induced pro-inflammatory signals. J Biol Chem 280(48):39914–39924
Wong CW, Christen T, Roth I, Chadjichristos CE, Derouette J-P, Foglia BF, Chanson M, Goodenough DA, Kwak BR (2006) Connexin37 protects against atherosclerosis by regulating monocyte adhesion. Nat Med 12(8):950–954
Xu D, Hopf C, Reddy R, Cho RW, Guo L, Lanahan A, Petralia RS, Wenthold RJ, O'Brien RJ, Worley P (2003) Narp and NP1 form heterocomplexes that function in developmental and activity-dependent synaptic plasticity. Neuron 39(3):513–528
Zeis T, Graumann U, Reynolds R, Schaeren-Wiemers N (2008) Normal-appearing white matter in multiple sclerosis is in a subtle balance between inflammation and neuroprotection. Brain 131(Pt 1):288–303
We thank the Netherlands Brain Bank team for their excellent work and Anke Essing for technical support.
This study was financed by MS Research grant MS03-525ext and Vriendenloterij grant MS09–666.
Availability of data and materials
Microarray data have been uploaded in the Gene Expression Omnibus (GEO) database. The GEO accession number to access the complete dataset is: GSE126802.
Ethics approval and consent to participate
The Netherlands Brain Bank received permission to perform autopsies and to use tissue and medical records from the Ethical Committee of the VU University medical center (VUmc, Amsterdam, The Netherlands). All donors have given informed consent for autopsy and use of their brain tissue for research purposes.
Consent for publication
All authors have approved the manuscript for submission.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.