Next Article in Journal
The Occurrence of a Negative Energy Balance in Holstein-Friesian and Simmental Cows and Its Association with the Time of Resumption of Reproductive Activity
Next Article in Special Issue
Identification of Serum Oxylipins Associated with the Development of Coronary Artery Disease: A Nested Case-Control Study
Previous Article in Journal
Effects of Acute and Chronic Resistance Exercise on the Skeletal Muscle Metabolome
Previous Article in Special Issue
Metabolomics and Inflammatory Mediator Profiling for the Differentiation of Life-Threatening and Non-Severe Appendicitis in the Pediatric Population
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 12
Issue 5
10.3390/metabo12050446
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Circulating Metabolites as Biomarkers of Disease in Patients with Mesial Temporal Lobe Epilepsy

by
Alexandre B. Godoi
1,2,
Amanda M. do Canto
1,2,
Amanda Donatti
1,2,
Douglas C. Rosa
1,2,
Danielle C. F. Bruno
1,2,
Marina K. Alvim
2,3,
Clarissa L. Yasuda
2,3,
Lucas G. Martins
4,
Melissa Quintero
4,
Ljubica Tasic
4,
Fernando Cendes
2,3 and
Iscia Lopes-Cendes
1,2,*
1
Department of Translational Medicine, School of Medical Sciences, University of Campinas (UNICAMP), Campinas 13083-888, Brazil
2
Brazilian Institute of Neuroscience and Neurotechnology (BRAINN), Campinas 13083-888, Brazil
3
Department of Neurology, School of Medical Sciences, University of Campinas (UNICAMP), Campinas 13083-888, Brazil
4
Department of Organic Chemistry, Institute of Chemistry, University of Campinas (UNICAMP), Campinas 13083-888, Brazil
*
Author to whom correspondence should be addressed.
Metabolites 2022, 12(5), 446; https://doi.org/10.3390/metabo12050446
Submission received: 20 March 2022 / Revised: 6 May 2022 / Accepted: 12 May 2022 / Published: 17 May 2022
(This article belongs to the Special Issue Omics Technologies in Evaluating New Insights into Molecular Aspects of Human Health)
Browse Figures
Review Reports Versions Notes

Abstract

:
A major challenge in the clinical management of patients with mesial temporal lobe epilepsy (MTLE) is identifying those who do not respond to antiseizure medication (ASM), allowing for the timely pursuit of alternative treatments such as epilepsy surgery. Here, we investigated changes in plasma metabolites as biomarkers of disease in patients with MTLE. Furthermore, we used the metabolomics data to gain insights into the mechanisms underlying MTLE and response to ASM. We performed an untargeted metabolomic method using magnetic resonance spectroscopy and multi- and univariate statistical analyses to compare data obtained from plasma samples of 28 patients with MTLE compared to 28 controls. The patients were further divided according to response to ASM for a supplementary and preliminary comparison: 20 patients were refractory to treatment, and eight were responsive to ASM. We only included patients using carbamazepine in combination with clobazam. We analyzed the group of patients and controls and found that the profiles of glucose (p = 0.01), saturated lipids (p = 0.0002), isoleucine (p = 0.0001), β-hydroxybutyrate (p = 0.0003), and proline (p = 0.02) were different in patients compared to controls (p < 0.05). In addition, we found some suggestive metabolites (without enough predictability) by multivariate analysis (VIP scores > 2), such as lipoproteins, lactate, glucose, unsaturated lipids, isoleucine, and proline, that might be relevant to the process of pharmacoresistance in the comparison between patients with refractory and responsive MTLE. The identified metabolites for the comparison between MTLE patients and controls were linked to different biological pathways related to cell-energy metabolism and pathways related to inflammatory processes and the modulation of neurotransmitter release and activity in MTLE. In conclusion, in addition to insights into the mechanisms underlying MTLE, our results suggest that plasma metabolites may be used as disease biomarkers. These findings warrant further studies exploring the clinical use of metabolites to assist in decision-making when treating patients with MTLE.

Graphical Abstract

1. Introduction

Epilepsy is a chronic neurological disorder characterized by persistent and long-lasting hyperactivity of groups of neurons, increasing one’s propensity to develop epileptic seizures [1,2]. In mesial temporal lobe epilepsy (MTLE), the epileptogenic focus is localized in the medial structures of the temporal lobe, mainly in the hippocampus [3,4,5]. MTLE is the most common type of focal epilepsy in adults. About 67–89% of patients with MTLE with hippocampal sclerosis have drug-resistant epilepsy and do not respond to currently available antiseizure medications (ASMs) [6,7]. In this sense, it has been determined that patients who do not reach seizure freedom when taking two ASMs (in monotherapy or combination therapy) may be defined as drug-resistant [8]. In addition to being exposed to the seizures’ intrinsic deleterious effects, pharmacoresistant patients experience the severe side effects of the ASM, often prescribed as polytherapy and in high doses [9]. Therefore, a significant challenge in the clinical management of patients with MTLE is identifying those who do not respond to ASM, allowing for the pursuit of other types of treatment, such as epilepsy surgery [10,11]. Indeed, recent reports have highlighted the importance of identifying biomarkers for response to therapy in epilepsies, including MTLE [12]. Several hypotheses explain the mechanism of ASM resistance. The current consensus is that this is a multifactorial condition where gene–gene and gene–environment interactions play important roles [12].
Metabolomics based on Nuclear Magnetic Resonance (1H-NMR) is used for complex sample analysis, principally because of the high reproducibility of NMR data, little or no sample preparation, and sample integrity preservation, which is particularly relevant when evaluating rare clinical samples [13,14]. In addition, NMR-based metabolomics may identify changes in biological samples linked to morphological and biochemical alterations associated with disease, thus assisting in early and precise diagnosis. Furthermore, identifying metabolites associated with specific phenotypes can contribute to a better understanding of disease mechanisms [15].
This study investigated the feasibility of using plasma metabolites as disease and pharmacoresistance biomarkers in patients with MTLE. In addition, we used the metabolomics data generated to gain insights into the mechanisms underlying MTLE and response to ASM.

2. Results

2.1. Characteristics of Study Population

We ascertained 28 patients with a mean age of 54 years (ranging from 26 to 70); eight patients were responsive to treatment with ASM, reaching seizure freedom when using the combination of CBZ + CLB, and 20 were classified as refractory to treatment with this therapy, as they did not acquire seizure freedom. The main clinical characteristics of the patients studied are reported in Table 1. We only included patients currently using CBZ + CLB. However, all patients considered here as pharmacoresistant have failed several other ASM regimens. None of the patients presented generalized seizures 24 h before blood collection. Most patients presented signs of hippocampal sclerosis on magnetic resonance image-decreased hippocampus volume in T1 images with increased signal in T2/FLAIR [3]. Age at onset of epilepsy varied from 1 to 30 years old. We also studied samples from 28 controls with a mean age of 49 years. These were adults between the ages of 29 and 63 years.

2.2. Metabolomic Analysis

Overall, we identified 27 metabolites in the samples studied (Appendix A, Table A1). Additionally, we found different groups of biomolecules, such as lipids, in the plasma of patients and controls. These biomolecules corresponded to different saturations’ fatty acids; amino acids such as alanine, isoleucine, leucine, valine, and glutamine; and aromatic amino acids such as tyrosine, histidine, phenylalanine, lactate, and glucose (Figure 1). See Appendix A, Table A1, and Figure A1 for the chemical shift assignments.
The largest variation in the acquired NMR data directions was visualized using the PLS-DA method. The PLS-DA scores plot (Figure 2A) clusters data for the patients and controls and the first principal component (1), explaining 39.1% of the data variation. The PLS-DA model goodness of fit (R2) value was 0.83, with a prediction (Q2) of 0.23. The metabolites that discriminate between the groups were identified based on the variables’ importance in the projection score (VIP score) from the PLS-DA analysis. The first 15 VIPs were combined with univariate analysis such as t-test and fold change (FC) (Table 2). We also evaluated the AUC of the ROC curve of the most important features of the model to calculate their predictability power (Figure 3). As a result, we identified five significantly altered metabolites when comparing patients with controls: glucose, saturated lipids, isoleucine, β-hydroxybutyrate, and proline (Table 2 and Figure 2C).
Furthermore, we also built a model to compare the patients based on their pharmacological responses. In this sense, a PLS-DA model was developed for comparing refractory with responsive MTLE patients (Figure A3), displaying scores plot with clustering for the two groups and the first principal component (1), which explains 44.6% of the data variation.
In addition, a constructed model alone was not enough to consider the model prediction. However, despite not observing considerable predictability for our method in PLS-DA parameters, we proposed some candidate variables according to their PLS regression coefficients. Those were classified based on the variables’ relative VIP. These metabolites were lipoproteins, lactate, glucose, proline, isoleucine, and unsaturated lipids (Table 3 and Figure A3C). Nonetheless, these metabolites did not present significant differences by t-test when comparing patients with responsive and refractory MTLE.
A comparative evaluation of the metabolome between these two groups of patients, as presented here, could potentially reveal metabolic features related to the process of pharmacoresistance. However, as we did not obtain sufficient significance and predictive power in the statistical model, the information found for such a comparison, having lactate, glucose, proline, isoleucine, and unsaturated lipids only suggestively support their importance in such a condition. Thus, for a more assertive elucidation of the discriminating metabolites between patients with refractory and responsive MTLE, such an assessment should be pursued in a larger sample size.
Additionally, we performed pathway enrichment analysis with the differentially abundant metabolites to compare patients and controls. For such an analysis, we identified as main altered pathways: lactose degradation, glucose-alanine cycle, lactose synthesis, transfer of acetyl groups into mitochondria, glycolysis, gluconeogenesis, fatty-acid biosynthesis, galactose metabolism, sphingolipid metabolism, arginine, and proline metabolism, Warburg effect, valine, leucine, and isoleucine degradation pathway as main altered pathways in patients (Figure 4A).
We also analyzed the diseases to which the differently abundant metabolites found in the comparison between patients and controls have been previously linked. They are diabetes mellitus (non-insulin-dependent), persistent hyperinsulinemic hypoglycemia, pyruvate carboxylase deficiency, respiratory chain deficiencies, and long-chain 3-hydroxyacyl-CoA dehydrogenase (LCHAD) deficiency—for a complete overview of the diseases found to be linked to the differentially abundant metabolites, refer to Appendix A, Figure A2.
Finally, we built networks correlating metabolites with their putative-related genes using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. We found 18 genes associated with the metabolites capable of discriminating patients from controls (Figure 4B). Then, this set of genes was used to search the human phenotype ontology gene setlist and the Gene Set Enrichment Analysis (GSEA) database. We found genes related to elevated hepatic transaminase, energetic metabolism, amino-acid metabolism, and inflammatory processes.

3. Discussion

We used 1HNMR-based metabolomics to investigate plasma samples of patients with MTLE and controls. In addition, we compared patients with different responses to ASM, i.e., those refractory and responsive to ASM. We searched for metabolic traits that could discriminate between the groups using a hypothesis-free design. Furthermore, we aimed to gain insights into the mechanisms underlying MTLE and response to treatment with ASM with the metabolomics data obtained. Because different epilepsy syndromes present different mechanisms, which may influence the biological processes leading to pharmacoresistance [12], we performed our study exclusively on patients with a single, well-defined epilepsy syndrome, MTLE. Furthermore, the use of ASM may lead to different metabolomics signatures [16]; thus, we selected patients using the same therapeutic regimen, CBZ + CLB, minimizing the source of bias. Moreover, the occurrence of seizures in proximity to the time of sample collection may temporarily affect the metabolic profile in plasma; hence, we only collected samples from patients who were seizure-free 24 h before sample collection.
We found five altered metabolites comparing patients and controls: glucose, saturated lipids, 3-hydroxybutyrate, isoleucine, and proline. These biomolecules are related to metabolic processes involved mainly in energetic and amino-acid metabolisms, especially those related to the cell’s anaplerotic reactions. Noteworthy is the presence of lipids as one of the discriminating elements between patients and controls.
In this study, we identified lower glucose levels in patients compared to controls. In this sense, changes in the energetic metabolism of neurons and glial cells have been frequently reported in neurological disorders. They are related to the preferential and intense usage of glucose by the nervous tissue. In addition, neurotransmitters, such as glutamate, acetylcholine, and gamma-aminobutyric acid, are dependent on energy metabolism [17]. The main molecule responsible for transferring and storing energy, the adenosine 5’-triphosphate (ATP), is produced through the glucose catabolism in the glycolysis, tricarboxylic acid (TCA) cycle, oxidative phosphorylation, and several other pathways related to the energetic metabolism, which were also present in the enriched pathways analysis performed in this study (Figure 4). Moreover, it is well-known that energy metabolism is indispensable for neurotransmitter production and release and the activity of ion channels. These are the essential players in the known mechanisms leading to epilepsy [18,19]. Therefore, an energetic failure due to abnormalities in glucose concentration may lead to increased seizure susceptibility [20]. Indeed, previous 1H-NMR studies revealed abnormal glucose concentrations in adult patients under different combinations of ASM [16] and in drug-free pediatric patients [21], presenting focal and generalized seizure types. Unfortunately, these studies presented important limitations since the diagnosis of epilepsy syndrome was not informed for both studies [16,21], and one of the studies included patients using different types of ASM [16]. Furthermore, none of the studies discussed above [16,21] informed if the patients presented seizures 24 h before sample collection, which may have introduced significant heterogeneity and bias toward identifying altered metabolites due to the ASM used and/or/or after a major seizure.
In addition, other energy sources become necessary in an energy-failure scenario since the decrease in ATP availability requires the recruitment of other biomolecules for energy production. In this context, lipidic support is crucial to cell maintenance [22]. In our study, some lipids corresponding to saturated fatty acids have been found elevated in patients and can be related to altered energetic metabolism since this group of molecules, abundantly found in lipoproteins, such as VLDL, HDL, and LDL, can be constantly displaced to many different tissues to supply the energy input [23].
Other important energy sources for neural cells are gluconeogenesis, glycogenolysis, and beta-oxidation. The ketone bodies can also perform this role, being synthesized depending on the metabolic disposition of the cell in prolonged periods of energetical deficit [24]. These molecules—acetoacetate, acetone, and β-hydroxybutyrate—contribute mainly to the brain-tissue energy input, as they can cross the blood–brain barrier, supplying about 60% of the central nervous system (CNS) [25]. Our study found increased levels of β-hydroxybutyrate in patients, even though signals of saturated lipids overlay the chemical shift identified for this molecule. This finding also indicates a state of energetic failure in the cells of patients with MTLE.
In addition to the well-known therapeutic strategy of the ketogenic diet in epilepsy, diets rich in BCAA have been shown to have positive effects on seizure control in an animal model [26] and combined with the ketogenic diet in children with epilepsy [27]. However, a recent study failed to demonstrate that chronic ingestion of BCAA is effective in the long-term reduction of seizures in an animal model of MTLE [28]. It has been proposed that an imbalance in the concentrations of amino acids such as valine, leucine, and isoleucine may lead to increased glutamate levels in the nervous tissue, and consequently, to increased hyperexcitability and excitotoxicity, resulting in neuroinflammation, which is known to occur in epilepsy [27,29,30]. Interestingly, we found increased levels of isoleucine in patients compared to controls, indicating a complex relationship between isoleucine concentration in MTLE and response to ASM treatment.
Furthermore, we found increased proline levels in patients compared to controls. This nonpolar amino acid has been linked to seizures and cognitive dysfunction [31], as well as to PYCR2 gene expression linked to microcephaly and hypomethylation [32], along with the PARS2 gene related to infantile-onset encephalopathy [33]. Furthermore, isoleucine is a metabolite linked to genes involved in neurological disorders, such as dementia and seizures: BCAT1, BCAT2, and IARS2 (26980008 and 30098844).
Interestingly, we found that the metabolic profile of patients with MTLE was also linked to non-insulin-dependent diabetes mellitus (Mody). This association has been reported previously, and it may be disease-induced or linked to the chronic use of ASM [34,35,36]. Indeed, a disruption of glucose metabolism in epilepsy leading to the unusual use of glucose as a source of energy [20,37,38] and pro-inflammatory processes occurring in epilepsy [38,39,40] may be risk factors for the development of diabetes in patients with MTLE.
However, the influence of pharmacological therapy on the group of patients cannot be discarded. This set of alterations, mainly linked to the amino acids and lipidic metabolism pathways found in the comparison between patients and controls, could be related to the ASM use, not exclusively the disease hallmarks. Despite the low evidence of metabolic alterations in humans under the use of the ASMs addressed in this study, assessments on the influence of CBZ in animal models have already been reported. In this context, alterations in amino acids and carnitine metabolisms (related to the lipid metabolism) in models of exposure to different doses of CBZ using Zebrafish [41] and in a bivalve mollusc Mytilus galloprovincialis [42] were already described.
Overall, our study found that the levels of metabolites such as glucose, saturated lipids, β-hydroxybutyrate, isoleucine, and proline, can distinguish patients with MTLE from controls. In addition, the computational integration of the data obtained in the metabolomic analysis suggests that pathways related to energetic metabolism, excitatory neurotransmission, and inflammatory processes are abnormally regulated in patients with MTLE. Notably, recent studies using tissue from animal models of MTLE reported similar altered biological pathways and processes [38,39,40]. Finally, we found that isoleucine and proline, amino acids present in abnormal levels in patients, may be involved in the increased susceptibility to seizures.

Strengths and Limitations

Strengths of this study include a well-characterized and homogeneous cohort since only patients with confirmed MTLE using the same type of ASM were included in the study. In addition, all patients were followed prospectively by the same group of neurologists in an epilepsy clinic at a university hospital. In addition, we only included patients without generalized seizures in the last 24 h of blood collection. The diagnostic criteria used to identify patients with MTLE and characterize the response to treatment with ASM followed the recommendations by the International League Against Epilepsy [8,43]. Furthermore, we used an agnostic approach to identify metabolites in the patients’ plasma and controls through a robust, highly reproducible analytical technique. To our knowledge, this is the first study to comply with all the elements listed above.
Our study also presents some limitations. First, because we recruited such a homogenous group of patients and controls, our study counted on a limited number of investigated individuals. Besides that, the access to these patients’ biofluids depends on their active follow-up in the epilepsy clinic, which increases the difficulty of obtaining them. This occurs mainly with the responsive patients, since they do not often return to the hospital due to their positive response to prescribed ASMs. Therefore, this reduced number of recruited individuals may have impacted the weak significance and predictability observed in our findings in comparing patients with responsive and refractory MTLE.
Another important limitation in our study is the inherent influence of the ASMs on the metabolic profile when comparing patients to controls. Thus, the correlation between some of the metabolites identified may not be due to the disease but to the recurrent use of the ASMs well.

4. Conclusions

Overall, plasma levels of glucose, saturated lipids, β-hydroxybutyrate, isoleucine, and proline can distinguish patients with MTLE from controls. Most interestingly, glucose levels, isoleucine, proline, lipoproteins, unsaturated lipids, and lactate are different in patients with pharmacoresistant MTLE compared to responsive MTLE. In addition, the computational integration of the data obtained in the metabolomic analysis suggests that pathways related to energetic metabolism, excitatory neurotransmission, and inflammatory processes are abnormally regulated in patients with MTLE. The identified metabolites were linked to biological pathways related to cell-energy metabolism, inflammatory processes, and modulation of neurotransmitter release and activity, which can potentially contribute to the mechanisms underlying MTLE. Thus, in addition to insights into the mechanisms underlying MTLE, we provided a supplementary analysis comparing patients with refractory and responsive MTLE, showing suggestive evidence that plasma metabolites may be used as biomarkers of response to ASM in patients with MTLE. These findings warrant further studies exploring the clinical use of metabolites to assist in medical decisions in treating patients with MTLE.

5. Materials and Methods

5.1. Study Population

We studied 28 patients with MTLE at the University of Campinas (UNICAMP) hospital’s outpatient epilepsy clinic. All patients had confirmed MTLE diagnoses according to International League Against Epilepsy (ILAE) criteria [43] and were investigated using a standard protocol, including serial interictal electroencephalogram high-resolution 3 Tesla magnetic resonance imaging [44,45]. We only included adults between the ages of 18 and 65. Patients were classified into two groups according to the patients’ responses to the treatment with ASM refractory or responsive following ILAE recommendations [8]. In addition, we only included patients using carbamazepine (CBZ) with clobazam (CLB), the therapy most frequently used by our epilepsy clinic patients. None presented generalized seizures 24 h before blood collection. Furthermore, we excluded patients using drugs for secondary clinical conditions that could have interacted with the ASM. We also studied samples from 28 control individuals. These were adults between 18 and 65 who were not in treatment for any neurological or chronic clinical conditions. We also excluded individuals under therapy with medications that could alter the function of key metabolism enzymes, such as cytochrome P450 isoforms.

5.2. Blood Collection and Plasma Preparation

Peripheral blood (4 mL) was collected in standard EDTA tubes (Vacutainer; Becton Dickinson, Franklin Lakes, NJ, USA). Blood samples were kept on ice for up to two hours until the plasma was separated and centrifuged at 2500 rpm for 10 min at 4 °C. Then, the obtained plasma was aliquoted and stored at −80 °C until further analysis.

5.3. 1H-NMR Spectroscopy Analyses

Samples were thawed in an ice bath (4 °C), diluted 1:1 (v/v) in deuterium oxide (D2O, 99.9% with 0.03% of trimethylsilyl propanoic acid, TSP, from Sigma Aldrich, St. Louis, MO, USA), and transferred into 5 mm NMR tubes. A Bruker AVANCE III 600 spectrometer (Bruker Biospin, Karlsruhe, Germany), equipped with a TBI (Triple Resonance Broadband Inverse) probe, was used. All 1H-NMR spectra were acquired at 25 °C, for example, Carr–Purcell–Meiboom–Gill (CPMG) pulse sequence (cpmgpr1d), with 128 scans, was explored further. The two-dimensional (2D) total correlation spectroscopy (TOCSY) experiments were recorded with mlevgpphw5 pulse sequence and 300 scans of randomly selected samples.

5.4. Data Analysis: NMR Data Processing

The obtained spectra were processed and standardized using MestreNova software (MestrelabResearch S.L.). First, the spectra were referenced using TSP (0.0 ppm), then aligned, and their phase was corrected. Next, peaks were normalized by the total area, and the matrix was prepared using spectral bins (0.005 ppm). A normalization with a constant sum (100) of the entire spectrum intensity was employed, reducing the differences in concentration between the plasma samples. In addition, that, water and EDTA spectral regions were removed.
Such a data matrix was processed using MetaboAnalyst 5.0 platform (Xia, McGill University) [46]. Data filtering, normalization by the sum, and Pareto scaling (mean-centered and divided by the square root of the standard deviation of each variable) were used in data preprocessing before the statistical tests, and 56 samples against 1620 variables (bins) for 1H-NMR CPMG were analyzed.
Furthermore, we performed multivariate principal component analysis (PCA), partial least-squares discriminant analysis (PLS-DA), orthogonal (o) PLS-DA, leave-one-out cross-validation (LOOCV), and variable importance in projection (VIP) scores. The first fifteen VIPs were correlated with univariate analyses, such as fold change (FC) and the t-test (p-value < 0.05). The area under the curve (AUC) of the receiver operator characteristics (ROC) curve of the selected features was also used to evaluate their prediction power by using the MetaboAnalyst 5.0 platform (Xia, McGill University) [46]. The chemical compounds were identified with the support of the Human Metabolome Database (HMDB) h [47] and Biological Magnetic Resonance Data Bank (BMRB) [48]. Metabolic and disease pathways analyses were performed using the Metaboanalyst function Metabolic Set Enrichment Analysis (MSEA) with an over-representation analysis (ORA) algorithm. We also built a network of compound-gene relations using the MetScape app from Cytoscape [49] to identify enriched pathways and visualize changes in metabolite data. All the above enrichment analyses were generated using the Kyoto Encyclopedia of Genes and Genomes (KEGG) [50] and the Small Molecule Pathway Database (SMPDB) [51].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo12050446/s1.

Author Contributions

A.B.G., A.M.d.C. and I.L.-C. designed the research; A.B.G., A.M.d.C., A.D., D.C.F.B., M.K.A., C.L.Y., F.C. and I.L.-C. selected the cases, collected and provided samples, revised and provided clinical data of patients; A.B.G., A.M.d.C., A.D., D.C.R., M.Q. and L.G.M. performed the NMR experiments; A.B.G. and M.Q. performed the chemometrics analysis; L.T. provided reagents and supervised NMR metabolomic experiments; I.L.-C. provided reagents and supervised the study; A.B.G., A.M.d.C., M.Q., L.T. and I.L.-C. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Brazil (Grant numbers: 2013/07559-3, and 2018/24069-3), and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Brazil (Grant number 001). A.B.G. and D.C.R. received scholarships from FAPESP, Grant numbers: 2019/00213-0, and 2019/00048-0. A.M.C., A.D., D.C.F.B., and M.D. received fellowships from FAPESP, Grant numbers, 2019/25948-3, 2021/04441-8, 2017/26167-0, and 2021/04451-3. I.L.-C. is supported by a grant from Conselho Nacional de Pesquisa (CNPq), Brazil (Grant number: 311923/2019-4).

Institutional Review Board Statement

The project was approved by the local medical research review board and ethics committee (file number 12112913.3.0000.5404/257.020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data analyzed in this study has been presented in the manuscript and in the Supplemental Material.

Acknowledgments

We thank all the staff of the epilepsy outpatient clinic at the University of Campinas Hospital (UNICAMP, Brazil) for the technical support during the assessment of patients for this study. We are very grateful to the patients and their families for their collaboration.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table
Table A1. 1H-NMR chemical shifts assignments of the metabolites found in the plasma of patients and controls. Metabolites were identified according to the variables of importance in projection (VIPs) by the PLS-DA model are shown underlined for the comparison between patients and controls and between patients with refractory MTLE and responsive MTLE.
Table A1. 1H-NMR chemical shifts assignments of the metabolites found in the plasma of patients and controls. Metabolites were identified according to the variables of importance in projection (VIPs) by the PLS-DA model are shown underlined for the comparison between patients and controls and between patients with refractory MTLE and responsive MTLE.
Sugars
1Glucose3.41 (H4′), 3.54 (H2′), 3.72 (H3′), 3.84 (H5′), 5.24 (H1′)
3.25 (H2′), 3.41 (H4′), 3.47 (H5′), 3.72 (H6′), 3.90 (H6′), 4.65 (H1′)
Amino acids
2Isoleucine0.94 (δ-CH3), 1.01 (γ-CH2), 1.42 (γ-CH2), 3.67 (α-CH)
3Leucine0.96 (δ-CH3), 1.71 (β-CH2)
4Valine0.99 (γ-CH2), 1.04 (γ-CH3), 2.25–2.31 (β-CH), 3.60 (α-CH)
5Alanine1.48 (β-CH3), 3.78 (α-CH)
6Glutamate2.04 (β-CH2), 2.36 (γ-CH2), 3.71 (α-CH)
7Glutamine2.13 (β-CH2), 2.45 (γ-CH2), 3.78 (α-CH)
8Arginine1.72 (γ-CH2), 1.89 (β-CH2), 3.23 (δ-CH2), 3.73 (α-CH)
9Proline3.36 (δ′δ′-CH2), 3.41 (δ-CH2), 4.14 (α-CH)
10Histidine7.06 (H4), 7.77 (H2)
11Tyrosine6.89 (H3, H5), 7.19 (H2, H6)
12Phenylalanine3.24 (β′β′-CH), 7.32,7.36 (H2, H6), 7.42 (H3, 5H)
Lipids
13Low-density lipoproteins (LDL)0.90 (CH3), 1.30 (CH2)n
14Very-low-density lipoproteins (VLDL)0.76–0.93 (CH3), 1.24–1.37 (CH2)n
15High-density lipoproteins (HDL)0.80–0.85 (CH3), 1.21–1.23 (CH2)n
16Fatty-acid chain1.53–1.65 (-CH2CH2CO)
17Unsaturated lipids and N-acetyl glycoproteins1.96–2.09 (-CH2-CH=)
18Lipids HC=CH5.29–5.43 (CH)
Organic acids
19Lactate1.33 (CH3), 4.11 (CH)
20Acetic acid1.92 (CH2)
21Acetone2.24 (CH2)
22Acetoacetic acid2.28 (CH2)
23Citrate2.54 (CH2), 2.68 (CH2)
24Formate8.47 (CH)
Other compounds
25Ethanol1.20 (CH3); 3.67 (CH2)
26Creatinine3.94 (CH2), 3.04 (CH2)
27Creatine4.06 (CH2), 3.05 (CH2)
Metabolites 12 00446 g0a1 550
Figure A1. 2D TOCSY 1H-1H NMR obtained for one plasma sample. Spectral region 0.50 to 5.50 ppm is shown on the upper panel, and the region from 6.6 to 7.9 ppm is given on the bottom panel. The following metabolites were identified: 1. Glucose; 3. Leucine; 4. Valine; 5. Alanine; 6. Glutamate; 7. Glutamine; 8. Arginine; 9. Proline; 10. Histidine; 11. Tyrosine; 12. Phenylalanine; 13. LDL: Low-Density Lipoproteins; 14, VLDL: Very Low-Density Lipoproteins; 16. Fatty-Acid Chain; 17. Unsaturated lipids and N-acetyl-glycoproteins; 20. Lactate; 23. Citrate. (Numbers refer to ones shown in Table A1).
Figure A1. 2D TOCSY 1H-1H NMR obtained for one plasma sample. Spectral region 0.50 to 5.50 ppm is shown on the upper panel, and the region from 6.6 to 7.9 ppm is given on the bottom panel. The following metabolites were identified: 1. Glucose; 3. Leucine; 4. Valine; 5. Alanine; 6. Glutamate; 7. Glutamine; 8. Arginine; 9. Proline; 10. Histidine; 11. Tyrosine; 12. Phenylalanine; 13. LDL: Low-Density Lipoproteins; 14, VLDL: Very Low-Density Lipoproteins; 16. Fatty-Acid Chain; 17. Unsaturated lipids and N-acetyl-glycoproteins; 20. Lactate; 23. Citrate. (Numbers refer to ones shown in Table A1).
Metabolites 12 00446 g0a1
Metabolites 12 00446 g0a2 550
Figure A2. List of diseases identified by the Metabolite Set Enrichment Analysis (MSEA) using the metabolites found to discriminate between patients and controls.
Figure A2. List of diseases identified by the Metabolite Set Enrichment Analysis (MSEA) using the metabolites found to discriminate between patients and controls.
Metabolites 12 00446 g0a2
Metabolites 12 00446 g0a3 550
Figure A3. Multivariate analysis of 1H-NMR (CPMG) plasma spectra of responsive and refractory MTLE patients. (A) PLS-DA analysis with the accuracy of 44.6%, R2 0.13, and Q2 −0.19. (B) O-PLS-DA analysis. The information about response to treatment with ASM in patients with MTLE was not implemented into the models. (C) Box plots representing the variations of the relative concentrations (measured as peak intensities) of metabolites whose VIP scores > 2 according to PLS-DA results. The black dots represent the metabolite levels in all samples, and the yellow diamond represents the average value. Responsive MTLE (blue) and refractory MTLE (red).
Figure A3. Multivariate analysis of 1H-NMR (CPMG) plasma spectra of responsive and refractory MTLE patients. (A) PLS-DA analysis with the accuracy of 44.6%, R2 0.13, and Q2 −0.19. (B) O-PLS-DA analysis. The information about response to treatment with ASM in patients with MTLE was not implemented into the models. (C) Box plots representing the variations of the relative concentrations (measured as peak intensities) of metabolites whose VIP scores > 2 according to PLS-DA results. The black dots represent the metabolite levels in all samples, and the yellow diamond represents the average value. Responsive MTLE (blue) and refractory MTLE (red).
Metabolites 12 00446 g0a3

References

  1. Beghi, E. The epidemiology of epilepsy. Neuroepidemiology 2020, 54, 185–191. [Google Scholar] [CrossRef] [PubMed]
  2. Scheffer, I.E.; Berkovic, S.; Capovilla, G.; Connolly, M.B.; French, J.; Guilhoto, L.; Hirsch, E.; Jain, S.; Mathern, G.W.; Moshe, S.L.; et al. ILAE classification of the epilepsies: Position paper of the ILAE Commission for Classification and Terminology. Epilepsia 2017, 58, 512–521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Cendes, F.; Sakamoto, A.C.; Spreafico, R.; Bingaman, W.; Becker, A.J. Epilepsies associated with hippocampal sclerosis. Acta Neuropathol. 2014, 128, 21–37. [Google Scholar] [CrossRef]
  4. Engel, J. Mesial temporal lobe epilepsy: What have we learned? Neuroscientist 2001, 7, 340–352. [Google Scholar] [CrossRef] [PubMed]
  5. Fisher, R.S.; Acevedo, C.; Arzimanoglou, A.; Bogacz, A.; Cross, J.H.; Elger, C.E.; Engel, J.; Forsgren, L.; French, J.A.; Glynn, M.; et al. ILAE Official Report: A practical clinical definition of epilepsy. Epilepsia 2014, 55, 475–482. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Semah, F.; Picot, M.C.; Adam, C.; Broglin, D.; Arzimanoglou, A.; Bazin, B.; Cavalcanti, D.; Baulac, M. Is the underlying cause of epilepsy a major prognostic factor for recurrence? Neurology 1998, 51, 1256–1262. [Google Scholar] [CrossRef] [PubMed]
  7. Chipaux, M.; Szurhaj, W.; Vercueil, L.; Milh, M.; Villeneuve, N.; Cances, C.; Auvin, S.; Chassagnon, S.; Napuri, S.; Allaire, C.; et al. Epilepsy diagnostic and treatment needs identified with a collaborative database involving tertiary centers in France. Epilepsia 2016, 57, 757–769. [Google Scholar] [CrossRef] [Green Version]
  8. Kwan, P.; Arzimanoglou, A.; Berg, A.T.; Brodie, M.J.; Allen Hauser, W.; Mathern, G.; Moshé, S.L.; Perucca, E.; Wiebe, S.; French, J. Definition of drug resistant epilepsy: Consensus proposal by the ad hoc Task Force of the ILAE Commission on Therapeutic Strategies. Epilepsia 2010, 51, 1069–1077. [Google Scholar] [CrossRef]
  9. Laxer, K.D.; Trinka, E.; Hirsch, L.J.; Cendes, F.; Langfitt, J.; Delaney, N.; Resnick, T.; Benbadis, S.R. The consequences of refractory epilepsy and its treatment. Epilepsy Behav. 2014, 37, 59–70. [Google Scholar] [CrossRef] [Green Version]
  10. Lee, A.T.; Burke, J.F.; Chunduru, P.; Molinaro, A.M.; Knowlton, R.; Chang, E.F. A historical cohort of temporal lobe surgery for medically refractory epilepsy: A systematic review and meta-analysis to guide future nonrandomized controlled trial studies. J. Neurosurg. 2019, 133, 71–78. [Google Scholar] [CrossRef]
  11. De Souza, J.P.S.A.S.; Ayub, G.; Nogueira, M.; Zanao, T.; Lopes, T.M.; Pimentel-Silva, L.R.; Domene, V.; Marquez, G.; Yasuda, C.L.; Ribeiro, L.F.; et al. Temporopolar amygdalohippocampectomy: Seizure control and postoperative outcomes. J. Neurosurg. 2020, 134, 1044–1053. [Google Scholar] [CrossRef]
  12. Bruxel, E.M.; Canto, A.M.D.; Bruno, D.C.F.; Geraldis, J.C.; Lopes-Cendes, I. Multi-omic strategies applied to the study of pharmacoresistance in mesial temporal lobe epilepsy. Epilepsia Open 2021. [Google Scholar] [CrossRef] [PubMed]
  13. Stanisic, D.; Martins, L.G.; Tasic, L. Nuclear magnetic resonance spectroscopy in analyses of biological samples. In Tools and Trends in Bioanalytical Chemistry; Springer: Cham, Switzerland, 2022; ISBN 978-3-030-82381-8. [Google Scholar]
  14. Emwas, A.-H.M.; Salek, R.M.; Griffin, J.L.; Merzaban, J. NMR-based metabolomics in human disease diagnosis: Applications, limitations, and recommendations. Metabolomics 2013, 9, 1048–1072. [Google Scholar] [CrossRef]
  15. Gebregiworgis, T.; Powers, R. Application of NMR metabolomics to search for human disease biomarkers. Comb. Chem. High Throughput Screen. 2012, 15, 595–610. [Google Scholar] [CrossRef] [PubMed]
  16. Murgia, F.; Muroni, A.; Puligheddu, M.; Polizzi, L.; Barberini, L.; Orofino, G.; Solla, P.; Poddighe, S.; Del Carratore, F.; Griffin, J.L.; et al. Metabolomics as a tool for the characterization of drug-resistant epilepsy. Front. Neurol. 2017, 8, 459. [Google Scholar] [CrossRef] [PubMed]
  17. Deutch, A.Y.; Roth, R.H. Pharmacology and biochemistry of synaptic transmission. In From Molecules to Networks; Academic Press: Cambridge, MA, USA, 2014; pp. 207–237. ISBN 978-0-12-397179-1. [Google Scholar]
  18. Dhar-Chowdhury, P.; Malester, B.; Rajacic, P.; Coetzee, W.A. The regulation of ion channels and transporters by glycolytically derived ATP. Cell. Mol. Life Sci. 2007, 64, 3069–3083. [Google Scholar] [CrossRef]
  19. Boison, D.; Steinhäuser, C. Epilepsy and astrocyte energy metabolism. Glia 2017, 66, 1235–1243. [Google Scholar] [CrossRef] [PubMed]
  20. Reid, C.A.; Mullen, S.; Kim, T.H.; Petrou, S. Epilepsy, energy deficiency and new therapeutic approaches including diet. Pharmacol. Ther. 2014, 144, 192–201. [Google Scholar] [CrossRef]
  21. Boguszewicz, Ł.; Jamroz, E.; Ciszek, M.; Emich-Widera, E.; Kijonka, M.; Banasik, T.; Skorupa, A.; Sokół, M. NMR-based metabolomics in pediatric drug resistant epilepsy—Preliminary results. Sci. Rep. 2019, 9, 15035. [Google Scholar] [CrossRef] [Green Version]
  22. McGuire, K.; Beerman, K.A. Nutritional Sciences: From Fundamentals to Food; Thomson/Wadsworth: Belmont, CA, USA, 2007; ISBN 049512656X. [Google Scholar]
  23. Kersten, S.; Wahli, W. Peroxisome proliferator activated receptor agonists. In New Approaches to Drug Development; Jollès, P., Ed.; Experientia Supplementum; Birkhäuser Basel: Basel, Switzerland, 2000; Volume 89, pp. 141–151. ISBN 978-3-0348-9547-7. [Google Scholar]
  24. Cotter, D.G.; Schugar, R.C.; Crawford, P.A. Ketone body metabolism and cardiovascular disease. Am. J. Physiol. Heart Circ. Physiol. 2013, 304, H1060–H1076. [Google Scholar] [CrossRef] [Green Version]
  25. Zhang, Y.; Kuang, Y.; LaManna, J.C.; Puchowicz, M.A. Contribution of brain glucose and ketone bodies to oxidative metabolism. In Oxygen Transport to Tissue XXXIV; Welch, W.J., Palm, F., Bruley, D.F., Harrison, D.K., Eds.; Advances in Experimental Medicine and Biology; Springer: New York, NY, USA, 2013; Volume 765, pp. 365–370. ISBN 978-1-4614-4771-9. [Google Scholar]
  26. Hartman, A.L.; Santos, P.; O’Riordan, K.J.; Stafstrom, C.E.; Hardwick, M. Potent anti-seizure effects of D-leucine. Neurobiol. Dis. 2015, 82, 46–53. [Google Scholar] [CrossRef] [Green Version]
  27. Evangeliou, A.; Spilioti, M.; Doulioglou, V.; Kalaidopoulou, P.; Ilias, A.; Skarpalezou, A.; Katsanika, I.; Kalamitsou, S.; Vasilaki, K.; Chatziioanidis, I.; et al. Branched chain amino acids as adjunctive therapy to ketogenic diet in epilepsy: Pilot study and hypothesis. J. Child Neurol. 2009, 24, 1268–1272. [Google Scholar] [CrossRef] [PubMed]
  28. Gruenbaum, S.E.; Dhaher, R.; Rapuano, A.; Zaveri, H.P.; Tang, A.; De Lanerolle, N.; Eid, T. Effects of branched-chain amino acid supplementation on spontaneous seizures and neuronal viability in a model of mesial temporal lobe epilepsy. J. Neurosurg. Anesthesiol. 2019, 31, 247–256. [Google Scholar] [CrossRef] [PubMed]
  29. Dufour, F.; Nalecz, K.A.; Nalecz, M.J.; Nehlig, A. Modulation of absence seizures by branched-chain amino acids: Correlation with brain amino acid concentrations. Neurosci. Res. 2001, 40, 255–263. [Google Scholar] [CrossRef]
  30. Skeie, B.; Petersen, A.J.; Manner, T.; Askanazi, J.; Steen, P.A. Effects of valine, leucine, isoleucine, and a balanced amino acid solution on the seizure threshold to picrotoxin in rats. Pharmacol. Biochem. Behav. 1994, 48, 101–103. [Google Scholar] [CrossRef]
  31. Roussos, P.; Giakoumaki, S.G.; Bitsios, P. A risk PRODH haplotype affects sensorimotor gating, memory, schizotypy, and anxiety in healthy male subjects. Biol. Psychiatry 2009, 65, 1063–1070. [Google Scholar] [CrossRef]
  32. Nakayama, T.; Al-Maawali, A.; El-Quessny, M.; Rajab, A.; Khalil, S.; Stoler, J.M.; Tan, W.-H.; Nasir, R.; Schmitz-Abe, K.; Hill, R.S.; et al. Mutations in PYCR2, encoding pyrroline-5-carboxylate reductase 2, cause microcephaly and hypomyelination. Am. J. Hum. Genet. 2015, 96, 709–719. [Google Scholar] [CrossRef] [Green Version]
  33. Yin, X.; Tang, B.; Mao, X.; Peng, J.; Zeng, S.; Wang, Y.; Jiang, H.; Li, N. The genotypic and phenotypic spectrum of PARS2-related infantile-onset encephalopathy. J. Hum. Genet. 2018, 63, 971–980. [Google Scholar] [CrossRef]
  34. Lu, C.-L.; Chang, Y.-H.; Sun, Y.; Li, C.-Y. A population-based study of epilepsy incidence in association with type 2 diabetes and severe hypoglycaemia. Diabetes Res. Clin. Pract. 2018, 140, 97–106. [Google Scholar] [CrossRef]
  35. Nair, S.S.; Harikrishnan, S.; Sarma, P.S.; Thomas, S.V. Metabolic syndrome in young adults with epilepsy. Seizure 2016, 37, 61–64. [Google Scholar] [CrossRef] [Green Version]
  36. Pylvänen, V.; Pakarinen, A.; Knip, M.; Isojärvi, J. Insulin-related metabolic changes during treatment with valproate in patients with epilepsy. Epilepsy Behav. 2006, 8, 643–648. [Google Scholar] [CrossRef]
  37. McCorry, D.; Nicolson, A.; Smith, D.; Marson, A.; Feltbower, R.G.; Chadwick, D.W. An association between type 1 diabetes and idiopathic generalized epilepsy. Ann. Neurol. 2006, 59, 204–206. [Google Scholar] [CrossRef] [PubMed]
  38. Canto, A.M.; Vieira, A.S.; Matos, A.H.B.; Carvalho, B.S.; Henning, B.; Norwood, B.A.; Bauer, S.; Rosenow, F.; Gilioli, R.; Cendes, F.; et al. Laser microdissection-based microproteomics of the hippocampus of a rat epilepsy model reveals regional differences in protein abundances. Sci. Rep. 2020, 10, 4412. [Google Scholar] [CrossRef] [Green Version]
  39. Canto, A.M.; Matos, A.H.B.; Godoi, A.B.; Vieira, A.S.; Aoyama, B.B.; Rocha, C.S.; Henning, B.; Carvalho, B.S.; Pascoal, V.D.B.; Veiga, D.F.T.; et al. Multi-omics analysis suggests enhanced epileptogenesis in the Cornu ammonis 3 of the pilocarpine model of mesial temporal lobe epilepsy. Hippocampus 2020, 31, 122–139. [Google Scholar] [CrossRef] [PubMed]
  40. Vieira, A.S.; De Matos, A.H.; Canto, A.M.D.; Rocha, C.S.; Carvalho, B.S.; Pascoal, V.D.B.; Norwood, B.; Bauer, S.; Rosenow, F.; Gilioli, R.; et al. RNA sequencing reveals region-specific molecular mechanisms associated with epileptogenesis in a model of classical hippocampal sclerosis. Sci. Rep. 2016, 6, 22416. [Google Scholar] [CrossRef] [Green Version]
  41. Ribbenstedt, A.; Posselt, M.; Benskin, J.P. Toxicometabolomics and biotransformation product elucidation in single zebrafish embryos exposed to carbamazepine from environmentally-relevant to morphologically altering doses. Chem. Res. Toxicol. 2022, 35, 431–439. [Google Scholar] [CrossRef] [PubMed]
  42. Dumas, T.; Courant, F.; Almunia, C.; Boccard, J.; Rosain, D.; Duporté, G.; Armengaud, J.; Fenet, H.; Gomez, E. An integrated metabolomics and proteogenomics approach reveals molecular alterations following carbamazepine exposure in the male mussel Mytilus galloprovincialis. Chemosphere 2022, 286, 131793. [Google Scholar] [CrossRef]
  43. Blümcke, I.; Thom, M.; Aronica, E.; Armstrong, D.D.; Bartolomei, F.; Bernasconi, A.; Bernasconi, N.; Bien, C.G.; Cendes, F.; Coras, R.; et al. International consensus classification of hippocampal sclerosis in temporal lobe epilepsy: A Task Force report from the ILAE Commission on Diagnostic Methods. Epilepsia 2013, 54, 1315–1329. [Google Scholar] [CrossRef]
  44. Yasuda, C.L.; Morita, M.E.; Alessio, A.; Pereira, A.R.; Balthazar, M.L.F.; Saude, A.V.; Costa, A.L.F.; Cardoso, T.A.; Betting, L.E.; Guerreiro, C.A.M.; et al. Relationship between environmental factors and gray matter atrophy in refractory MTLE. Neurology 2010, 74, 1062–1068. [Google Scholar] [CrossRef]
  45. Yasuda, C.L.; Tedeschi, H.; Oliveira, E.L.P.; Ribas, G.C.; Costa, A.L.C.; Cardoso, T.A.M.O.; Montenegro, M.A.; Guerreiro, C.A.; Guerreiro, M.M.; Li, L.M.; et al. Comparison of short-term outcome between surgical and clinical treatment in temporal lobe epilepsy: A prospective study. Seizure 2006, 15, 35–40. [Google Scholar] [CrossRef] [Green Version]
  46. Pang, Z.; Chong, J.; Zhou, G.; de Lima Morais, D.A.; Chang, L.; Barrette, M.; Gauthier, C.; Jacques, P.-É.; Li, S.; Xia, J. MetaboAnalyst 5.0: Narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 2021, 49, W388–W396. [Google Scholar] [CrossRef]
  47. Wishart, D.S.; Feunang, Y.D.; Marcu, A.; Guo, A.C.; Liang, K.; Vázquez-Fresno, R.; Sajed, T.; Johnson, D.; Li, C.; Karu, N.; et al. HMDB 4.0: The Human Metabolome Database for 2018. Available online: http://www.hmdb.ca/ (accessed on 15 March 2022).
  48. Ulrich, E.L.; Akutsu, H.; Doreleijers, J.F.; Harano, Y.; Ioannidis, Y.E.; Lin, J.; Livny, M.; Mading, S.; Maziuk, D.; Miller, Z.; et al. BioMagResBank. Available online: http://www.bmrb.wisc.edu/ (accessed on 25 March 2022).
  49. Karnovsky, A.; Weymouth, T.; Hull, T.; Tarcea, V.G.; Scardoni, G.; Laudanna, C.; Sartor, M.A.; Stringer, K.A.; Jagadish, H.V.; Burant, C.; et al. Metscape 2 bioinformatics tool for the analysis and visualization of metabolomics and gene expression data. Bioinformatics 2012, 28, 373–380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef] [PubMed]
  51. Frolkis, A.; Knox, C.; Lim, E.; Jewison, T.; Law, V.; Hau, D.D.; Liu, P.; Gautam, B.; Ly, S.; Guo, A.C.; et al. SMPDB: The small molecule pathway database. Nucleic Acids Res. 2010, 38, D480–D487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Metabolites 12 00446 g001 550
Figure 1. The representative 1H-NMR spectra of plasma samples of control, responsive, and refractory MTLE subjects. Spectral regions from 0.5 to 9.0 ppm, acquired using CPMG (cpmgpr1d) pulse sequence, are shown. The 7.0 to 8.5 ppm regions were zoomed in (8×) for better visualization. The regions of D2O and EDTA were removed (//). The following metabolites were identified: lipoproteins; leucine, valine; isoleucine; lactate; alanine; acetate; N-acetyl-glycoproteins; O-acetyl-glycoproteins; glutamine; glucose; tyrosine; histidine; phenylalanine; formate.
Figure 1. The representative 1H-NMR spectra of plasma samples of control, responsive, and refractory MTLE subjects. Spectral regions from 0.5 to 9.0 ppm, acquired using CPMG (cpmgpr1d) pulse sequence, are shown. The 7.0 to 8.5 ppm regions were zoomed in (8×) for better visualization. The regions of D2O and EDTA were removed (//). The following metabolites were identified: lipoproteins; leucine, valine; isoleucine; lactate; alanine; acetate; N-acetyl-glycoproteins; O-acetyl-glycoproteins; glutamine; glucose; tyrosine; histidine; phenylalanine; formate.
Metabolites 12 00446 g001
Metabolites 12 00446 g002 550
Figure 2. Multivariate analysis of 1H-NMR (CPMG) plasma spectra. (A) PLS-DA results with the accuracy of 76%, R2 0.83, and Q2 0.23. (B) O-PLS-DA results. The information about response to treatment with ASM in patients with MTLE was not implemented into the models. (C) Box plots representing the variations of the relative concentrations (measured as peak intensities) of metabolites whose VIP scores > 2 according to PLS-DA results. The black dots represent the metabolite levels in all samples, and the yellow diamond represents the average value. Patients with MTLE (red) and controls (blue). * = overlaid signals.
Figure 2. Multivariate analysis of 1H-NMR (CPMG) plasma spectra. (A) PLS-DA results with the accuracy of 76%, R2 0.83, and Q2 0.23. (B) O-PLS-DA results. The information about response to treatment with ASM in patients with MTLE was not implemented into the models. (C) Box plots representing the variations of the relative concentrations (measured as peak intensities) of metabolites whose VIP scores > 2 according to PLS-DA results. The black dots represent the metabolite levels in all samples, and the yellow diamond represents the average value. Patients with MTLE (red) and controls (blue). * = overlaid signals.
Metabolites 12 00446 g002
Metabolites 12 00446 g003 550
Figure 3. Receiver operator characteristics (ROC) curves and their respective areas under the curve (AUC) were calculated for the most important features determined by VIP values to compare MTLE patients and controls. * = overlaid signals.
Figure 3. Receiver operator characteristics (ROC) curves and their respective areas under the curve (AUC) were calculated for the most important features determined by VIP values to compare MTLE patients and controls. * = overlaid signals.
Metabolites 12 00446 g003
Metabolites 12 00446 g004 550
Figure 4. Metabolic pathways were identified by the Metabolite Set Enrichment Analysis (MSEA) of the discriminant metabolites that were identified when comparing (A) patients and controls. (B) Network compound genes built with the discriminant metabolites comparing patients and controls. Blue circles—annotated genes; pink hexagon—metabolite’s ligands; red hexagon—discriminant metabolites present in the KEGG database.
Figure 4. Metabolic pathways were identified by the Metabolite Set Enrichment Analysis (MSEA) of the discriminant metabolites that were identified when comparing (A) patients and controls. (B) Network compound genes built with the discriminant metabolites comparing patients and controls. Blue circles—annotated genes; pink hexagon—metabolite’s ligands; red hexagon—discriminant metabolites present in the KEGG database.
Metabolites 12 00446 g004
Table
Table 1. The main characteristics of the 28 patients with MTLE included in the study. All patients were using the ASM combination, CBZ + CLB, and most of them presented signs of hippocampal sclerosis on magnetic resonance imaging—decreased hippocampus volume in T1 images with increased signal in T2/FLAIR images.
Table 1. The main characteristics of the 28 patients with MTLE included in the study. All patients were using the ASM combination, CBZ + CLB, and most of them presented signs of hippocampal sclerosis on magnetic resonance imaging—decreased hippocampus volume in T1 images with increased signal in T2/FLAIR images.
IDSexAge (Years)Age at Onset of EpilepsyHippocampal AbnormalitiesGroupResponse to Treatment with ASM Combination, CBZ + CLB
1F6315LHSMTLERefractory
2M6026BilateralMTLERefractory
3F571LHSMTLERefractory
4M5915LHSMTLERefractory
5F7017LHSMTLERefractory
6M5814RHSMTLERefractory
7F502LHSMTLERefractory
8M5019LHSMTLERefractory
9F377LHSMTLERefractory
10M265LHSMTLERefractory
11M5410LHSMTLERefractory
12F6016RHSMTLERefractory
13F267RHSMTLERefractory
14F6223BilateralMTLERefractory
15M6020RHSMTLERefractory
16M6214LHSMTLERefractory
17F612RHSMTLERefractory
18M461LHSMTLERefractory
19F548LHSMTLERefractory
20F5130NoneMTLERefractory
21F437LHSMTLEResponsive
22M458RHSMTLEResponsive
23F653LHSMTLEResponsive
24M5531RHSMTLEResponsive
25F5817RHSMTLEResponsive
26M4719LHSMTLEResponsive
27F5620LHSMTLEResponsive
28F7018LHSMTLEResponsive
MTLE: mesial temporal lobe epilepsy; ID: patient identification; sex: male/female (M/F); age: the age at investigation. LHS: left hippocampal sclerosis; RHA: right hippocampal sclerosis; bilateral: bilateral hippocampal sclerosis.
Table
Table 2. Table showing the metabolites identified in different concentrations and their respective chemical shifts elucidated by the highest VIP values (VIP score). The p-values, calculated from the t-test, FC (fold change), and false discovery rate (FDR) are also shown.
Table 2. Table showing the metabolites identified in different concentrations and their respective chemical shifts elucidated by the highest VIP values (VIP score). The p-values, calculated from the t-test, FC (fold change), and false discovery rate (FDR) are also shown.
MTLE versus Control—VIP Score
MetabolitesChemical Shift (Multip; Assign.)Vip Scorep-ValueFCFDR
Glucose3.68–3.78 (m, CH)5.126.0000 × 10−31.300.196
Saturated Lipids 0.83–0.87 (m, CH3)5.020.0918 × 10−30.760.023
Saturated Lipids, Isoleucine *1.21–1.25 (m, -CH2-)3.800.0817 × 10−50.820.023
β-Hydroxybutyrate, Saturated Lipids *1.20–1.24 (m, -CH2-3.150.129 × 10−30.820.027
Unsaturated lipids, Isoleucine, Proline, and N-Acetyl-glycoproteins *1.96–2.09 (m, -CH2-CH=)2.160.0905 × 10−30.850.198
Multip = multiplicity, where s (singlet), d (doublet), t (triplet), dd (doublet of doublets), m (multiplet), l (broad); Assign. = assignment of these signals; * = overlaid signals.
Table
Table 3. List of metabolites identified in different concentrations of refractory and responsive MTLE and their respective chemical shifts elucidated by the highest VIP values (VIP score). The p-values, calculated from the t-test, FC (fold change) are also shown.
Table 3. List of metabolites identified in different concentrations of refractory and responsive MTLE and their respective chemical shifts elucidated by the highest VIP values (VIP score). The p-values, calculated from the t-test, FC (fold change) are also shown.
Refractory MTLE versus Responsive MTLE—VIP Scores
MetabolitesChemical Shift (Multip.; Assign.)VIP Scorep-ValueFC
Lipoproteins1.28 (m, CH)6.660.051.209
Lactate1.33 (d, CH3)5.410.051.159
Glucose3.41 (m, CH2)4.810.050.752
Exclusively unsaturated lipid2.06 (l, CH2CH=)1.710.050.857
Isoleucine0.94 (t, CH3)1.570.050.886
Proline3.36 (m, CH)1.130.050.658
Multip = multiplicity, s (singlet), d (doublet), t (triplet), dd (doublet of doublet), m (multiplet), l (broad); Assign. = assignment of these signals.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Godoi, A.B.; do Canto, A.M.; Donatti, A.; Rosa, D.C.; Bruno, D.C.F.; Alvim, M.K.; Yasuda, C.L.; Martins, L.G.; Quintero, M.; Tasic, L.; et al. Circulating Metabolites as Biomarkers of Disease in Patients with Mesial Temporal Lobe Epilepsy. Metabolites 2022, 12, 446. https://doi.org/10.3390/metabo12050446

AMA Style

Godoi AB, do Canto AM, Donatti A, Rosa DC, Bruno DCF, Alvim MK, Yasuda CL, Martins LG, Quintero M, Tasic L, et al. Circulating Metabolites as Biomarkers of Disease in Patients with Mesial Temporal Lobe Epilepsy. Metabolites. 2022; 12(5):446. https://doi.org/10.3390/metabo12050446

Chicago/Turabian Style

Godoi, Alexandre B., Amanda M. do Canto, Amanda Donatti, Douglas C. Rosa, Danielle C. F. Bruno, Marina K. Alvim, Clarissa L. Yasuda, Lucas G. Martins, Melissa Quintero, Ljubica Tasic, and et al. 2022. "Circulating Metabolites as Biomarkers of Disease in Patients with Mesial Temporal Lobe Epilepsy" Metabolites 12, no. 5: 446. https://doi.org/10.3390/metabo12050446

APA Style

Godoi, A. B., do Canto, A. M., Donatti, A., Rosa, D. C., Bruno, D. C. F., Alvim, M. K., Yasuda, C. L., Martins, L. G., Quintero, M., Tasic, L., Cendes, F., & Lopes-Cendes, I. (2022). Circulating Metabolites as Biomarkers of Disease in Patients with Mesial Temporal Lobe Epilepsy. Metabolites, 12(5), 446. https://doi.org/10.3390/metabo12050446

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop