|Home | About | Journals | Submit | Contact Us | Français|
Experimental Autoimmune Encephalomyelitis (EAE) is the most commonly used animal model for Multiple Sclerosis (MScl). CSF metabolomics in an acute EAE rat model was investigated using targetted LC–MS and GC–MS. Acute EAE in Lewis rats was induced by co-injection of Myelin Basic Protein with Complete Freund’s Adjuvant. CSF samples were collected at two time points: 10 days after inoculation, which was during the onset of the disease, and 14 days after inoculation, which was during the peak of the disease. The obtained metabolite profiles from the two time points of EAE development show profound differences between onset and the peak of the disease, suggesting significant changes in CNS metabolism over the course of MBP-induced neuroinflammation. Around the onset of EAE the metabolome profile shows significant decreases in arginine, alanine and branched amino acid levels, relative to controls. At the peak of the disease, significant increases in concentrations of multiple metabolites are observed, including glutamine, O-phosphoethanolamine, branched-chain amino acids and putrescine. Observed changes in metabolite levels suggest profound changes in CNS metabolism over the course of EAE. Affected pathways include nitric oxide synthesis, altered energy metabolism, polyamine synthesis and levels of endogenous antioxidants.
The online version of this article (doi:10.1007/s11306-011-0306-3) contains supplementary material, which is available to authorized users.
Multiple sclerosis (MScl) is a chronic demyelinating neurodegenerative central nervous system (CNS) disorder of autoimmune origin affecting over 1 million people worldwide. Development of the disease includes emergence of inflammatory lesions in white matter and dysfunction of neural conductivity. This leads to impairment of sensory and motor functions. Molecular processes associated with the onset and progression of MScl and its etiology are still unknown.
Because of obvious ethical and safety constraints, access to human brain and spinal cord system samples is very limited. Therefore in-depth investigation of molecular mechanisms of disorders within CNS is hampered. Biochemistry of CNS activity of affected patients can be monitored only indirectly, using cerebrospinal fluid (CSF) for diagnostic purposes. Also explorative studies of molecular mechanisms involved in the onset and progression of neurological diseases except for brain tumors are limited to CSF samples. The metabolites present in the CSF represent actual metabolism of CNS and the balance between blood and CSF. Its analysis can be helpful in the determination of markers for neurological disorders. However, sampling human CSF is an invasive procedure, hence access to such samples is much more limited compared to blood or urine. In order to partially avoid these limitations, animal models of neurological diseases have been introduced. An additional major advantage of animal models is less variation between individual animals compared to variation between individuals in clinical cohorts. Therefore the discovery of dominating factors of the onset and progression of the disease is more straightforward. On the other hand, biomarkers discovered by means of animal models may not be relevant for the human situation and such results need to be validated and confirmed using clinical data.
The most commonly used animal model of MScl is the experimental autoimmune encephalomyelitis (EAE). The course of EAE resembles the symptoms of MScl and can be used to study some aspects of MScl interconnected with neuroinflammation and neurodegeneration. In this model, disease is induced by inoculating animal with myelin-specific antigens. This paper describes a study investigating changes in metabolites caused by acute EAE; an autoimmune reaction was provoked using MBP co-injected with an adjuvant containing Mycobacterium tuberculosis proteins. In this model, gradual and reversible impairment of motor function is observed, however no clear myelin loss is observed.
Metabolomics, the comprehensive analysis of a wide range of metabolites, provides a novel perspective for the search of new disease biomarkers and drug targets, being an alternative and complementary approach to more established omics techniques such as genomics, transcriptomics or proteomics. Recent progress in metabolomics resulted in an increasing number of metabolomics applications in neurological research. In a recent review, Wishart et al. (2008) summarized the current status of CSF metabolomics reporting 308 metabolites together with their concentrations in CSF. Recent articles investigating CSF metabolome focused on untargeted approaches (Crews et al. 2009; Myint et al. 2009; Carrasco-Pancorbo et al. 2009; Wuolikainen et al. 2009) detecting high number of features. In this article we used a related approach emphasizing targeted metabolomics. We applied two platforms, allowing the fully validated analysis of 39 identified metabolites by LC–MS and 64 identified metabolites by GC–MS, with an overlap of 18 metabolites between both methods.
Male Lewis rats (Harlan Laboratories B.V., the Netherlands) kept under normal housing conditions with water and food ad libitum, weighing between 175 and 225 g at the start of the experiment, were inoculated on day 0 as previously described (Hendriks et al. 2004). Briefly, a 100 μl saline based emulsion containing 50 μl Complete Freund’s Adjuvant H37 RA (CFA, Difco Laboratories, Detroit, MI), 500 μl Mycobacterium tuberculosis type H37RA (Difco) and 20 μg guinea pig myelin basic protein (MBP) was injected subcutaneously in the pad of left hind paw of Isoflurane anaesthetized animals. Next to these MBP challenged rats, referred to as the EAE group, two control groups were included: a group of animals receiving the same emulsion without MBP (CFA group) and a healthy group undergoing anesthesia only (Healthy group, H). Each group consisted of 30 animals. Of each group, half of the animals were sacrificed to collect plasma and CSF on day 10 (day of onset of disease in EAE group) resulting in groups further referred as H10, CFA10 and EAE10, and the other half on day 14 (peak of disease in EAE group)—groups H14, CFA14 and EAE14.
Animals were grouped and housed three per cage and cages were randomized across treatments and disease duration. Disease symptoms and weights of all animals were recorded daily. The following scores for motor dysfunctions were used: 0, healthy animal with normal curling reflex at the tail; 1, paralysis of the tip of the tail; 2, loss of muscle tone at the base of the tail; 3, low posture of hind limbs; 4, instability at hips; 5, partial hind limb paralysis; 6, complete hind limb paralysis; 7, paralysis include midriff; 8, quadriplegia; 9, moribund; 10, death due to EAE. The animal experiments described were approved by the local Ethical Committee for Animal Experiments.
On day 10 and 14, animals were euthanized with CO2/O2 and the head of the rat was fixed in a holder. Terminal CSF samples were obtained by direct insertion of an insulin syringe needle (Myjector, 29G × 1/2″) via the arachnoid membrane into the Cisterna Magna. For this purpose a skin incision was made followed by a horizontal incision in the musculus trapezius pars descendens to reveal the Arachnoid membrane. A maximal volume of 60 μl was collected per animal. Each sample was centrifuged within 20 min after sampling, for 10 min at 2000×g at 4°C. After centrifugation, the supernatants were stored at −80°C for further analysis. Previous experiments have shown that collecting up to 60 μl using this technique and conditions provided hemoglobin-free CSF samples measured by ESI-Orbitrap (data not shown). As an additional check, fresh samples, supernatant and pellet size were visually scored for hemolysis and samples were discarded if positive (data not shown). After inspection, samples were aliquoted in volumes of 10 μl for analysis by different analytical platforms.
Due to the disturbed physiological state of animals of full EAE group at day 14 (EAE14), the success rate for CSF sampling for these animals was significantly lower, limiting the number of available samples. From 90 animals used in the experiment, 84 CSF samples were available for the LC–MS platform (11 in group EAE14) and 84 for the GC–MS platform (10 in group EAE14).
The LC–MS method employed is targeted to the analysis of primary and secondary amines using AccQ-Tag derivatization. 10 μl of CSF sample was spiked with 10 μl of an internal standard solution containing 13C15N-labeled amino acids, followed by addition of 100 μl of MeOH for deproteination. The mixture was vortexed for 10 s and centrifuged at 9400×g for 10 min at 10°C. The supernatant was transferred to a deactivated autosampler vial (Waters) and dried under N2. The residue was reconstituted in 80 μl of borate buffer (pH 8.5) vortexed for 10 s and treated with 20 μl of AccQ-TagAQC derivatization reagent (Waters). After keeping the mixture for 10 min at 55°C, the vials were transferred to an autosampler tray and cooled to 10°C until the injection. 1.0 μl of the reaction mixture was injected into the UPLC–MS/MS system.
An ACQUITY UPLC system with autosampler (Waters, Etten-Leur, The Netherlands) was coupled online with a Quattro Premier XE tandem quadrupole mass spectrometer (Waters) operated using Masslynx data acquisition software (version 4.1; Waters). The samples were analyzed by UPLC–MS/MS using an AccQ-Tag Ultra 100 mm × 2.1 mm (1.7 μm particle size) column (Waters). A binary gradient system of water–eluent A (10:1, v/v) (AccQ-Tag, Waters) and 100% eluent B (AccQ-Tag, Waters), was used. Elution of the analytes was achieved by ramping the percentage of eluent B from 0.1 to 90.0 in approx. 9.5 min using a combination of both linear and convex profiles. The flow-rate was 0.7 ml/min and the column temperature was maintained at 60°C. After each injection, the injection needle was washed with 200 μl 95:5 v/v ACN:water and subsequently with 600 μl 5:95 v/v ACN:water.
The Quattro Premier XE was used in the positive-ion electrospray mode and all analytes were monitored in Selective Reaction Monitoring (SRM) using nominal mass resolution; the width of the isolation/fragmentation window was set to 0.7 Da. All target compounds were monitored via the transition from the protonated molecule of the AccQ-Tag derivative to the common fragment at m/z 171. Collision energy and collision gas (Ar) pressure were 22 eV and 2.5 mbar, respectively. The identities of all compounds were confirmed by high resolution MS and identical retention times as authentic standards (data not shown). The acquisition run was divided into six time segments in order to decrease the number of simultaneous SRM transitions. All samples were analyzed in duplicate, and measurement order was randomized.
Acquired data were evaluated using Quanlynx software (Waters), by integration of assigned SRM peaks and normalization using proper internal standards. For analysis of amino acids, their 13C15N-labeled analogs were used. For other amines—the closest-eluting internal standard was employed. Resulting relative responses from duplicate injections were averaged after visual inspection.
The GC–MS platform used in this study requires 30 μl of rat CSF for sufficient sensitivity and coverage. Because only 10 μl CSF per rat was available, samples were pooled per three. This resulted in 30 μl pooled rat CSF samples and group sizes of 5 (groups H10, H14, CFA10, CFA14, EAE10) and 3 (group EAE14, limited number of samples available, see Section 2.2). The samples (30 μl) were deproteinized by adding 250 μl methanol and subsequently centrifuged for 10 min at 11800 rpm. The supernatant was dried under N2 followed by derivatization with methyl-N-(trimethylsilyl)-trifluoroacetamide (MSTFA) in pyridine similar to Koek et al. (2006). During the different steps in the sample work-up, i.e. prior to deproteinization, derivatization and injection, different (deuterated) internal standards were added at a level of approx. 20 ng/μl. The end volume was 45 and 1 μl aliquots of the derivatized samples were injected in splitless mode on a HP5-MS 30 m × 0.25 mm × 0.25 μm capillary column (Agilent Technologies, USA) using a temperature gradient from 70 to 320°C at a rate of 5°C/min. GC–MS analysis was performed using an Agilent 6890 gas chromatograph coupled to an Agilent 5973 mass selective detector (Agilent Technologies, USA). Detection was carried out using MS detection in electron impact mode and full scan monitoring mode (m/z 15–800). The electron impact for the generation of ions was 70 eV. Samples were analyzed in duplicate.
Data-pre-processing was carried out by composing target lists of peaks detected in the samples based on retention time and mass spectra and these peaks were integrated for all samples. All peak areas were subsequently normalized using internal standards. The resulting target lists were used for further statistical analysis. Identities were assigned based on the presence of identical mass spectra in an in-house database based on authentic standards (data not shown).
The significance of MBP treatment for all metabolites was assessed by linear modeling (Draper and Smith 1998). Three models were fitted for each metabolite. The first global model (Eq. 1) includes all data for each metabolite:
in which y = metabolite response; CFA = CFA injected (1 for CFA and EAE groups, otherwise 0); MBP = MBP injected (1 for EAE groups, otherwise 0); β0…β6 = model coefficients; time = number of days since start of the experiment (10 or 14); weightday0 = weight at day 0 (for LC–MS, individual rat LC–MS; for GC–MS, average of rats in pooled sample) [g].
To obtain a higher confidence, two additional single-day models were fitted for the samples collected on day 10 and day 14, respectively (Eq. 2):
in which (other symbols are explained above) ydayX = metabolite response at day X (X is 10 or 14) and α0…α3 = model coefficients.
For all three models, stepwise linear regression (Massart et al. 1997) was used to remove non-significant terms. The significance of the regression coefficients in the final models was obtained from a t test (H0: αi = 0, H1: αi ≠ 0 or H0: βi = 0, H1: βi ≠ 0).
A significant EAE effect in an individual single day model is present when there is a significant α3. However, here a distinction between a time dependent and a time independent effect cannot be made. For the combined model, a time-dependent effect is given by a significant β6 and a time-independent effect is given by a significant β3.
We distinguish between three effects: early stage EAE, full stage EAE and EAE progression. For an effect to be significant it must be present in both the combined model and at least one single day model, otherwise results are considered to be inconclusive.
A significant early stage EAE effect is reported when a significant α3 in the day 10 model is observed together with a significant β3 and/or β6 in the combined model. A significant α3 in the day 14 model together with a significant β3 and/or β6 in the combined model indicates full stage EAE. An EAE progression effect is reported when a significant β6 in the combined model is observed in correlation with a significant α3 in at least one of the single day models.
Fold changes are obtained by taking ratios between average model calculated metabolite responses per treatment group. The fold changes for the onset (EAE10/CFA10) and peak (EAE14/CFA14) of the disease are calculated from the respective single day model results (Eq. 2). Fold changes for EAE progression (EAE14/EAE10) are obtained from the global model results (Eq. 1).
Figure 1 shows the progression of animal growth (weight, left axis) for all experimental groups. For animals from EAE14, progression of EAE is also shown (neurological scores, right axis). The weight change during the course of the experiment is dependent on multiple factors and therefore difficult to use as quantitative measure in the data analysis. Nevertheless, it accounts for the overall well-being of the animal. Because of the fast growth of animals (4–6 g per day), the complex treatment and use of two sampling time-points, growth of the animals may affect the metabolic profile. In addition, after the onset of the disease, animals from EAE group started losing weight and this decrease is correlated with the increase of neurological scores (see Fig. 1). As a consequence, the weight of the animals from EAE14 group was significantly lower than from all the other groups. The exact cause for this effect is not certain. It is possible that progressing paralysis impaired the animals ability to feed themselves, however animals were able to eat and drink while lying. The degree to which weight loss influences the volume of CSF and concentration of metabolites is hard to access.
Principal component analysis (PCA) was applied for preliminary visual inspection of the metabolite profiles of the CSF samples belonging to the different groups (Figs. 2, ,3).3). Prior to PCA, each metabolite was scaled to zero mean and unit variance. Clear separation between the full-stage EAE group (EAE14) and the other groups is observed for both platforms. The GC–MS data (Fig. 3) also show partial separation of early stage EAE samples (group EAE10). For LC–MS (Fig. 2) this separation is not observed although some separation can be seen in higher principal components (see Figure S1 in supplementary material).
Table 1 summarizes the LC–MS linear modeling results. The table also indicates overlapping GC–MS targets and where the GC–MS linear modeling results confirm effect directionality. Although directionality is not always confirmed by linear modeling, visual inspection of the GC–MS data indicates agreement. There are no cases where LC–MS and GC–MS effects have opposite direction. The smaller groups in GC–MS as a result of sample pooling restrict the ability to draw unambiguously reliable biological conclusions based on that platform alone. Therefore, GC–MS linear modeling results are supplied as supplementary material.
Significant changes of metabolite levels detected by linear models were verified by pair wise comparison using double-sided t tests and visually inspected by boxplot per molecule. Figure 4 below shows a boxplot presenting concentration changes of putrescine—the metabolite that showed biggest fold-change between the experimental groups.
The design of the biological experiments, with three treatments groups—healthy controls, a CFA group with animals subjected to peripheral inflammation only and a full EAE treatment, allowed for separation of the concentration changes caused by EAE from possible variations due to adjuvant injection, including variations due to animal growth or peripheral inflammation.
Several significant changes in concentrations of metabolites both at the onset and on the peak of the disease were detected. Dissimilar patterns of changes observed on both time points may suggest activation of different metabolic pathways in the early and late stage of the disease. The early-stage EAE samples, taken on day 10 post inoculation resemble the moment of maximal CNS infiltration by blood monocytes and T-cells, however, at this moment no negative signs on health status are observed and animals do not have any neurological scores. Metabolite level changes in the early stage of the disease may be caused by processes related to the development of EAE and disruption of the blood-CSF barrier.
The full-stage EAE samples are taken on day 14, at the peak of the disease, when the paralysis of the animals is most severe, which is expressed in maximum neurological scores (see Fig. 1). The EAE model employed in this study involves spontaneous recovery of animals. In fact, due to the slightly heterogeneous course of the disease, 2 out of 15 animals from the EAE14 group had neurological scores which had already decreased relative to the previous day. Therefore, one can assume that some metabolites observed as significantly regulated in this sample series might be connected to pathological processes involved in EAE progression or recovery. The observed weight loss of EAE14 animals poses a difficult challenge in interpretation of observed changes. It was shown in the literature that weight loss and dietary restriction may cause changes in metabolite profiles in body fluids (Bollard et al. 2005). However, it is unknown how these changes affect CSF metabolic profile. For that reason the brief interpretation of biological relevance of observed changes presented below is mostly limited to compounds involved in processes previously associated with MScl and EAE.
Several metabolites targeted in this study and found significantly up/down-regulated were previously associated with molecular processes interconnected with development of MScl and EAE. Glutamic acid (glutamate), glutamine, gamma-aminobutyric acid (GABA), asparagine and taurine are involved in excitotoxicity and energy metabolism. Arginine and citrulline play a key role in NO synthesis. Arginine, ornithine and putrescine are involved in polyamine synthesis. Alanine and branched-chain amino acids (valine, leucine, isoleucine) are involved in energy metabolism. See Table 2 for overview.
l-glutamic acid (glutamate) is a ubiquitous neurotransmitter with several identified receptors and transport proteins both in neural and glial cells (Danbolt 2001). Over-stimulation of glutamate receptors may cause cell death due to excitotocity. Indeed, elevated levels of CSF glutamate were reported in clinical cases of multiple sclerosis (Sarchielli et al. 2003). Contrary to previous reports, as shown in Table 1, the level of glutamic acid in this study shows a statistically significant decrease in concentration at the peak of the disease, while no significant change is observed at the early stage. However, the glutamine, which is closely interconnected with glutamate through glutamine-glutamate cycle in neurons and astrocytes, shows down-regulation at the onset and up-regulation with full stage EAE. Conversion of glutamate into glutamine by glutamine synthetase (GS) enzyme localized in astrocytes is a key process protecting neurons from ammonia toxicity (Norenberg and Martinez-Hernandez 1979) and the only pathway of ammonia disposal due to the lack of complete urea cycle in CNS (Wiesinger 2007). Available evidence shows up-regulation of GS in astrocytes present in demyelinating lesions as assessed in a study of post mortem brain tissue of MS patients (Newcombe et al. 2008). However, a similar study performed in MBP-induced mouse EAE shows down-regulation of this enzyme (Hardin-Pouzet et al. 1997).
At the peak of the disease we observe a significant gamma-aminobutyric acid (GABA) increase. GABA is synthesized from glutamate by glutamate decarboxylase. It was shown that GABA, an inhibitory neurotransmitter, may attenuate glutamate-induced excitotoxicity (Ohkuma et al. 1994). In our study we observe significant increase in concentration of another inhibitory neurotransmitter, taurine, which was also shown to be protective against excitotoxicity (El Idrissi and Trenkner 1999).
It should also be emphasized that glutamate, besides having a profound role as a neurotransmitter in CNS, is one of the key metabolites in energy metabolism of neurons, glia and immune cells (Newsholme et al. 1999). It is possible that the global CSF levels of glutamine and glutamate involved in energy metabolism dominate over local high and excitotoxic concentrations in the most affected regions of brain and spinal cord. The observed discrepancy between available results on glutamate up- and down-regulation in EAE/MS suggests it may be involved in multiple different molecular processes associated with neuroinflammation and neurodegeneration.
Peroxynitrile anion (ONOO−) is a very reactive oxidizing agent, capable of inducing cell death through multiple pathways (Szabó et al. 2007). It emerges in the reaction between nitric oxide, produced mainly by inductive nitric oxide synthetase (iNOS) in activated immune cells and microglia, and the free radical superoxide O2−, coming from oxidative phosphorylation cycle in mitochondria. Currently available evidence suggests, besides promotion of peroxynitrile synthesis, that NO plays a complex regulatory role in mediating the immune response (van der Veen et al. 1997). Numerous reports suggest activity of peroxynitrile in active lesions in MScl (Smith and Lassmann 2002), as well as in EAE (van der Veen et al. 1997).
In our study we observe significant down-regulation of arginine, the main substrate for NO synthesis during the early stage of EAE. This effect may be interconnected with intense production of NO by cytotoxic T-cells, macrophages and activated microglia during development of neuroinflammation. However, the exact role of NO in EAE and MScl development remains unclear. Available reports show that inhibition of NO synthesis may attenuate (Ding et al. 1998) or enhance (Zielasek et al. 1995) EAE.
Interestingly, a second product of NO synthesis, citrulline, is also significantly down-regulated in the early stage of EAE, while one could expect increased production of citrulline linked to increased NO synthesis. However, immune and CNS cells are able to re-synthesize arginine in the so-called arginine-NO cycle (see Fig. 5, solid line), This cycle can be considered as a shortcut of the urea cycle (Fig. 5, dashed line), which due to the lack of ornithine carbamoyltransferase (OCT) is not complete in the CNS (Wiesinger 2001). It was demonstrated that arginine deficiency does not inhibit NO synthesis by T-cells, which under such circumstances are able to increase uptake of citrulline and up-regulate enzymes responsible for arginine synthesis (Bansal et al. 2004). Due to the involvement of arginine in other reactions, the citrulline–NO cycle does not operate in a stoichiometric manner, resulting in simultaneous decrease of arginine and citrulline levels.
It should be also noted that the most significant alternative reaction of arginine utilization is the synthesis of ornithine by the enzyme arginase (ARG). This reaction is known to take part in regulating NO synthesis by reducing the amount of available arginine (Bansal and Ochoa 2003). In the early stage of EAE we observe a decrease in the concentration of ornithine, supporting the hypothesis that arginine is probably being utilized by NO synthesis.
During the peak of EAE, arginine levels do not differ from controls, while citrulline levels continue to decrease. This may be explained by the reduction of NO synthesis interconnected with continued activity of other citrulline-NO cycle enzymes. At the same time we observe up-regulation of products of an alternative pathway of arginine utilization. The ornithine levels return to their normal states after significant decreases during the EAE onset, and this correlates with the significant 20-fold increase of putrescine levels, synthesised from ornithine by ornithine decarboxylase (see Fig. 4). These results may suggest a switch in arginine metabolism from nitric oxide production towards the synthesis of ornithine and further to polyamine synthesis, which could be associated with compensatory anti-inflammatory response and recovery. Interestingly, the rat with the highest observed level of putrescine is a clear outlier from the rest of the group (animal 76 on Fig. 4) and is one of two animals in the group that showed a neurological score decrease on the last day of experiment.
Putrescine was reported to be up-regulated in EAE (Bolton and Paul 2006), however this change was associated with disease progression rather than recovery mechanisms. On the other hand, the beneficial effect of increased polyamine synthesis on post-injury recovery of CNS was demonstrated in literature (Gilad et al. 1996; Cai et al. 2002). Available reports show the potential of putrescine to block NMDA glutamate receptors (Williams 1997) and therefore could possibly attenuate excitotoxic effects of glutamate.
It was postulated previously that over the course of inflammation, metabolism of arginine switches from NO synthesis in the early stage to arginine utilization in polyamine synthesis during tissue repair and wound healing (Nelin et al. 2007).
It was shown that antioxidants, such as uric acid, inhibit inflammation in EAE (Hooper et al. 2000). The inhibitory neurotransmitter taurine, observed as up-regulated in our study, is also known for its antioxidant and anti-inflammatory activity (Li et al. 2007). Levels of taurine increase over the course of EAE, suggesting that they may play a role in protection and recovery mechanisms. As the neurotoxicity of peroxynitrile might be dependent on the availability of superoxide rather than nitric oxide alone, it was suggested that MScl may be in fact connected with a malfunction of mitochondria (Mao and Reddy 2010) or imbalance of ATP synthesis and energy metabolism (Cadoux-Hudson et al. 1991; Amorini et al. 2009).
Besides changes discussed above, we also observe other changes in the metabolome that cannot be clearly connected to already established hypotheses of MScl/EAE mechanisms. In the early stage of the EAE we observe decreased levels of alanine and branched-chain amino-acids (BCAA: leucine, isoleucine and valine). Alanine and BCAAs are known to be utilized as a source for pyruvate for energy metabolism or de novo synthesis of macromolecules within neural (Hutson et al. 2007) and immune cells (Li et al. 2007). As the onset of EAE is associated with maximum infiltration of the CNS by blood monocytes and T-cells, the observed decrease of BCAAs and alanine may suggest these metabolites are utilized for energy metabolism by invading cells. Alternatively, pro-inflammatory cytokines released during EAE development may cause changes in energy metabolism of native CNS cells.
In the full stage EAE, among other changes, we observe a significant increase of O-phosphoethanolamine concentration, which is an important intermediate in phospholipid metabolism. Wheeler et al. (2008) reported changes in lipid contents in normal appearing white and grey matter of MScl patients, suggesting changes in lipid metabolism including increased sphingolipid turnover and phospholipid synthesis. The detected increase of O-phosphoethanolamine might be interconnected with these processes.
A targeted metabolomic approach of metabolites containing amino groups was applied to study changes in CNS metabolism in a MScl rat model. The obtained results showed significant changes in metabolic processes during the onset and progression of EAE. Excessive nitric oxide production along with changes in energy metabolism appear to play a key role in the development of the EAE, while a switch from arginine metabolism to polyamine synthesis may be associated with processes involved in spontaneous recovery of affected animals. Further research is required to confirm these hypotheses.
Results obtained in this study are partially in agreement with established knowledge on EAE and multiple sclerosis. However, in some cases we observe changes conflicting with previously published results. It should be noted that the vast majority of existing results were obtained in different experiments, executed in different paradigms, investigating similar, however still distinct biological conditions. Established EAE models differ from each other and still they do not cover all aspects of multiple sclerosis (Lassmann 2007). The ability to simultaneously measure a wide range of compounds is advantageous in systems biology approaches, like the one presented here. The use of a common experiment to probe multiple pathways in a biological system yields a chance for obtaining more comprehensive and coherent picture of molecular mechanisms involved.
Our results prove that metabolomics provides useful insights into the mechanisms underlying biological processes in the onset and development of multiple sclerosis using a rat model. High sensitivity and a wide range of metabolites were covered by the present analytical platforms, which allowed for a reliable screening of potential changes in CNS metabolism. In contrast to untargeted approaches, our LC/GC–MS analysis focuses on the measurement of well defined and confidently identified components. This allowed a more accurate description of the potential molecular mechanisms involved.
Authors would like to acknowledge Frans van der Kloet (LACDR) and Siem Heisterkamp (MSD) for batch design for LC–MS and GC–MS experiments, respectively. We also thank Amy Harms and Vanessa Gonzalez (both LACDR) for reading and commenting on this manuscript. This study was financially supported by Top Institute Pharma project D4-102. This study was supported by the research program of the Netherlands Metabolomics Centre (NMC), which is a part of The Netherlands Genomics Initiative/Netherlands Organization for Scientific Research.
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
Marek Noga and Adrie Dane contributed equally to this manuscript.