- 1Jining First People’s Hospital, Jining Medical University, Jining, China
- 2Affiliated Hospital of Jining Medical University, Jining Medical University, Jining, China
- 3Department of Clinical Translational Medicine, Jining Life Science Center, Jining, China
- 4Department of Pharmacy, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
Metabolomics is increasingly used to observe metabolic patterns and disease-specific metabolic biomarkers. However, serum metabolite analysis of moyamoya disease (MMD) is rarely reported. We investigated serum metabolites in MMD and compared them with those of healthy controls (HCs) using a non-targeted gas chromatography–mass spectrometry (GC–MS) approach to identify metabolic biomarkers associated with MMD. Forty-one patients with MMD diagnosed by cerebral angiography and 58 HCs were recruited for our study. Comparative analyses (univariate, multivariate, correlation, heatmaps, receiver operating characteristi curves) were performed between MMD patients and HCs. Twenty-five discriminating serum metabolic biomarkers between MMD patients and HCs were identified. Compared with HCs, MMD patients had higher levels of phenol, 2-hydroxybutyric acid, L-isoleucine, L-serine, glycerol, pelargonic acid, L-methionine, myristic acid, pyroglutamic acid, palmitic acid, palmitoleic acid, stearic acid, octadecanamide, monoglyceride (MG) (16:0/0:0/0:0), and MG (0:0/18:0/0:0), and lower levels of L-alanine, L-valine, urea, succinic acid, L-phenylalanine, L-threonine, L-tyrosine, edetic acid, and oleamide. These metabolic biomarkers are involved in several pathways and are closely associated with the metabolism of amino acids, lipids, carbohydrates, and carbohydrate translation. A GC–MS-based metabolomics approach could be useful in the clinical diagnosis of MMD. The identified biomarkers may be helpful to develop an objective diagnostic method for MMD and improve our understanding of MMD pathogenesis.
Introduction
Moyamoya disease (MMD) is a chronic, rare cerebrovascular-occlusive disease with a poor prognosis and considerable regional and racial differences (Kuroda and Houkin, 2008). The occurrence and development of MMD is multifactorial, and the underlying etiological and pathogenic mechanisms remain largely unclear (Huang et al., 2017). Most studies have focused on geographical distributions, sex predispositions, clinical manifestations, and age-specific characteristics (Okazaki et al., 2019; Sato et al., 2019). Basic research, including genomic and proteomic approaches, has been extensively conducted in the past 60 years. However, the physiological and metabolic mechanisms remain unclear.
China, especially southwest Shandong, has a large number of affected patients. MMD incidence peaks in children aged 5–9 years and adults aged 40–50 years with a female/male ratio of ∼2 and familial occurrence of 15% (Scott and Smith, 2009). The most common symptoms of MMD are headache, transient ischemic attack, infarction, and intracerebral hemorrhage. MMD often results in death, and early diagnosis is important. Cerebral angiography is the “gold standard” for diagnosis (Lee et al., 2013). However, cerebral angiography has flaws such as an inherent risk of anaphylaxis and nephropathy due to the use of contrast medium. The identification of new biomarkers to assist physicians in timely MMD diagnosis is needed to decrease morbidity and mortality.
Metabolomics is a promising approach that can identify metabolic patterns and disease-specific metabolic markers. This pivotal tool for biomarker discovery has improved the diagnosis of psychiatric illnesses (Cai et al., 2012), atherosclerosis (Chen et al., 2015), coronary artery disease (Turer et al., 2009), cerebral infarctions (Jung et al., 2011), and intracranial tumors (Pandey et al., 2017). However, information on metabolite profiles for MMD is limited. Only one metabolomic study on MMD patients has been conducted, and it focused on identifying metabolites in cerebrospinal fluid (CSF) using nuclear magnetic resonance (NMR) spectroscopy (Jeon et al., 2015). Therefore, more metabolomics studies are urgently needed to better understand the physiological and metabolic mechanisms involved in MMD.
We sought to elucidate the metabolic mechanisms underlying the occurrence and development of MMD. A GC–MS-based metabolomics approach coupled with uni- and multivariate analyses was employed to identify metabolic biomarkers associated with MMD. These results could help guide the development of an objective diagnostic method for MMD and provide insights into MMD pathogenesis.
Materials and Methods
Participants
Forty-one patients diagnosed by cerebral angiography as MMD and 58 HCs treated at the Jining First People’s Hospital, Jining Medical University between January 2018 and June 2019 were recruited for our study. Patients with MMD are often accompanied by basic diseases, such as hypertension, diabetes, and coronary artery disease. In order to make our research more rigorous, MMD patients with other types of diseases were not included in our study. In addition, the patients with MMD that we recruited were hospitalized patients who had a light diet 3 days before the operation. Similarly, 3 days’ diet of healthy people was recorded before we recruited to minimize the errors caused by these factors. Blood samples from all MMD patients were collected and subjected to centrifugation at 5,000 rpm for 10 min at room temperature and then stored at −80°C until analyses.
Materials and Instruments
Heptadecanoic acid (purity: ≥98%; lot no. SLBX4162), as an internal standard (IS), was from Sigma-Aldrich (St Louis, MO, United States). Methanol was of chromatographic grade and purchased from Thermo Fisher Scientific (Waltham, MA, United States). Water was purchased from Hangzhou Wahaha Company (Hangzhou, China). Pyridine (lot no. C10486013) was from Shanghai Macklin Biochemical (Shanghai, China). O-Methyl hydroxylamine hydrochloride (purity, 98.0%; lot no. 542171) was obtained from J&K Scientific Industries (Ambala, India). N,O-Bis(trimethylsilyl)trifluoroacetamide with 1% of trimethylchlorosilane (BSTFA + 1% TMCS) (v/v; lot no. BCBZ4865) was purchased from Sigma-Aldrich.
Preparation and Derivatization of Samples for GC–MS
Serum samples were processed according to the following procedure. First, 350 μl of methanol (containing 100 μg/ml of IS) was added to 100 μl of serum, vortexed, and centrifuged at 14,000 rpm for 10 min at 4°C. The supernatant was transferred to a 2-ml tube and evaporated to dryness at 37°C under the gentle flow of nitrogen gas. After the extracts had been dried, 80 μl of O-methyl hydroxylamine hydrochloride (15 mg/ml in pyridine) was added and mixed. The solution was incubated for 90 min at 70°C. Subsequently, 100 μl of BSTFA + 1% TMCS was added to each sample, followed by incubation for 60 min at 70°C. Samples were then detected by GC–MS.
GC–MS
GC–MS was done on a 7890B GC system equipped with a 7000 C mass spectrometer. Separation was conducted on an HP-5MS fused-silica capillary column (30 m × 0.25 mm × 0.25 μm) with high-purity helium as the carrier gas at a constant flowrate of 1.0 ml/min. Each 1-μl aliquot of derivatized solution was run in split mode (50:1), with helium as the carrier gas and a front inlet purge flow of 3 ml/min; the gas flowrate was 1 ml/min. The GC temperature program was set to begin at 60°C for 4 min, increased to 300°C at 8°C/min, with a final 5-min maintenance at 300°C. The temperatures of the injection, transfer line, and ion source were 280, 250, and 230°C, respectively. Electron impact ionization (−70 eV) was used, with an acquisition rate of 20 spectra/s in the MS setting. MS detection was conducted by electrospray ionization (ESI) in full-scan mode from mass/charge (m/z) values of 50–800.
In general, 20 μl of each serum sample from MMD and HCs was obtained and then vortexed and mixed. Thus, this is the so-called quality control (QC) sample. Fifteen samples and five QCs were performed for 1 day (one batch), and 15 samples run randomized while the five QCs were distributed evenly among the 15 samples. At last, nine samples and five QCs were done. The peak area and retention time (RT) of the IS (heptadecanoic acid) were applied to evaluate the stability of the sample injection.
Multivariate Statistical Analyses
GC–MS equipped with Unknowns Analysis and Agilent MassHunter Quantitative Analysis (for GC–MS) was used to process the GC data. This process enabled deconvolution, alignment, and data reduction to produce a list of m/z and RT pairs, with the corresponding intensities for all detected peaks from each data file in the dataset. The resulting table was exported into ExcelTM (Microsoft, Redmond, WA, United States), and the normalized peak area percentages were used as the percentage of corresponding intensities of each peak/total peak area. The resulting three-dimensional dataset including peak index (RT–m/z pairs), sample names (observations), and normalized peak area percentages was imported into SIMCA-P 14.0 (Umetrics, Umeå, Sweden) for statistical analyses. Center scaling, unit variance scaling, and pareto scaling are commonly used to perform the normalization data. In our paper, we adopted the pareto scaling. Too many missing values will cause difficulties for downstream analysis. There are several different methods for this purpose, such as replace by a small values, mean/median, k-nearest neighbor (KNN), probabilistic principal components analysis (PPCA), Bayesian PCA (BPCA) method, and singular value decomposition (SVD) method to impute the missing values (Kumar et al., 2017; Do et al., 2018). In our work, the default method replaces all the missing values with small values (the half of the minimum positive values in the original data) assuming to be the detection limit, and the data were not transformed.
A modified multicriteria assessment strategy was used to select variables. The assessment was used to reduce the number of variables and explore those that were most sensitive to interventions. The statistically significant threshold of variable importance in projection (VIP) values from the orthogonal partial least squares discriminant analysis (OPLS-DA) model was >1.0, and two-tailed Student’s t test differences of p < 0.05 were considered significant. “Fold change” was defined as the average mass response (area) ratio between two groups. Analyses of correlations, heatmaps, receiver operating characteristic (ROC) curves, and pathways were done using MetaboAnalyst 4.0.1 All results were shown in the “metabolome” view.
Results
Basic Participant Characteristics
Overall, 41 patients with MMD and 58 HCs were included in the analyses. Detailed information on the demographical and clinical characteristics of these participants is summarized in Table 1.
MMD patients were not significantly different from HCs with regard to age, sex, body mass index (BMI), tobacco smoking, or alcohol consumption (all p > 0.05). The types of MMD onset were cerebral infarction (39 cases) and hemorrhage (2 cases). In addition, there were 34 cases of bilateral MMD and 7 cases of unilateral MMD.
Metabolomics Analyses
Representative GC–MS total ion current (TIC) chromatograms of the QC serum sample showed strong signals (Figure 1). In addition, the relative standard deviation (RSD) in intra- and interday of the peak area and retention time (RT) of the IS were <15%, indicating that the analytical instrument operated within acceptable standard variations.
Figure 1. A representative gas chromatography–mass spectrometry (GC–MS) total ion chromatogram (TIC) of the quality control (QC) serum sample.
After analyses of unknown compounds and quantitative analyses, 114 metabolites were identified in each serum sample and then used in the subsequent multivariate analysis. The PCA scores plot for HCs and MMD patients were R2X = 0.816, Q2 = 0.301. The pairwise PLS-DA score plots also suggested that the MMD patients were statistically different from the HCs: R2X = 0.579, R2Y = 0.836, and Q2 = 0.603. OPLS-DA analyses were carried out to maximize discrimination. The results suggested that this model was efficient and clearly separated the MMD patients and HCs (R2X = 0.823, R2Y = 0.913, and Q2 = 0.783). Values approaching 1.0 indicate a stable model with predictive reliability. Additionally, a permutation test with 200 iterations verified that the constructed OPLS-DA model was valid and not overfitted, as the original R2 and Q2 values to the right were significantly higher than the corresponding permutated values to the left: R2 = 0.289, Q2 = −0.556. They are shown in Figure 2.
Figure 2. Multivariate statistical analysis between the control group and moyamoya disease (MMD) group. (A) Principal components analysis (PCA) scores plot; (B) partial least squares discriminant analysis (PLS-DA) scores plot; (C) orthogonal PLS-DA (OPLS-DA) scores plot; and (D) statistical validation of the OPLS-DA model through 200 × permutation testing.
Identification of Potential Biomarkers
Twenty-four metabolites to distinguish between the HCs and MMD patients were identified (VIP > 1, p < 0.05) (Table 2). Compared with HCs, patients with MMD were characterized by higher levels of phenol, 2-hydroxybutyric acid, L-isoleucine, L-serine, glycerol, pelargonic acid, L-methionine, pyroglutamic acid, myristic acid, palmitoleic acid, palmitic acid, stearic acid, octadecanamide, monoglyceride (MG) (16:0/0:0/0:0), and MG (0:0/18:0/0:0), as well as lower levels of L-alanine, L-valine, urea, succinic acid, L-threonine, L-phenylalanine, L-tyrosine, edetic acid, and oleamide. Their relationships were revealed by correlation analyses (Figure 3A). To better understand the metabolic differences between patients with MMD and HCs, data on identified metabolites were analyzed using clustering heatmaps. Even though sample clusters overlapped slightly, as shown in Figure 3B, most samples clearly grouped into two differentiated clusters, in agreement with OPLS analyses.
Table 2. List of assigned statistically significant metabolites between moyamoya disease (MMD) and healthy control (HC) group.
Figure 3. (A) Correlation analysis of the differential metabolites in moyamoya disease (MMD) patients and healthy controls (HCs). (B) Heat map for identified metabolites in MMD patients and HCs. The color of each section is proportional to the significance of change of metabolites (red, upregulated; blue, downregulated). Rows, samples; columns, metabolites.
ROC Curve Analyses
Further selection of potential indicator was performed by ROC analysis. Using different models, the value of the sensitivity and the area under the ROC curve (AUC) of these biomarker panels were both ≥0.8 (Figure 4A). An area of 1 represents a “perfect” test, so we obtained “good” efficiency for a clinical diagnosis for this set of metabolite biomarkers.
Figure 4. (A) Receiver operating characteristic (ROC) curves were for the sets of biomarker metabolites. (B) Summary of pathway analysis with MetaboAnalyst 4.0. (a) Aminoacyl-tRNA biosynthesis; (b) valine, leucine, and isoleucine biosynthesis; (c) propanoate metabolism; (d) phenylalanine metabolism; (e) cysteine and methionine metabolism; (f) alanine, aspartate, and glutamate metabolism; (g) phenylalanine, tyrosine, and tryptophan biosynthesis; (h) tyrosine metabolism; and (i) glycerolipid metabolism.
Analyses of Metabolic Pathways
We identified several pathways that may be significant (raw p < 0.5, impact >0) (Table 3). Nine pathways had the greatest significance: aminoacyl-tRNA biosynthesis: valine, leucine, and isoleucine biosynthesis; propanoate metabolism; phenylalanine metabolism; cysteine and methionine metabolism; alanine, aspartate, and glutamate metabolism; phenylalanine, tyrosine, and tryptophan biosynthesis; tyrosine metabolism; and glycerolipid metabolism (Figure 4B). The detailed results of the pathway analyses are shown in Table 3, with a summary shown in Figure 5.
Figure 5. Schematic diagram of the proposed metabolic pathways in moyamoya disease (MMD) compared to healthy controls (HCs). Red and blue represent up- and downregulated metabolites, respectively.
Discussion
Moyamoya disease is a multifactorial disorder that likely presents unique pathophysiological profiles in each individual. Genetics, proteomics, and imaging have been used to discover markers for MMD (Araki et al., 2010; Maruwaka et al., 2015). However, there is no existing marker that could aid the diagnosis of MMD.
GC–MS-based metabolomics was applied to profile metabolic biomarkers in the serum of 41 MMD patients and 58 HCs. Our study is the first to identify serum biomarkers in MMD patients. Novel biomarkers may assist researchers in understanding MMD pathogenesis and provide new therapeutic strategies.
Owing to high separation efficiency, sensitivity, specificity, and throughput, as well as the development of various derivatization technologies, GC–MS has become important for the application of metabolomics. However, a single technology that can cover all metabolites in biological systems is lacking, and each technology has its own technical advantages and disadvantages. We discovered more metabolites than that in an NMR study of CSF samples of MMD patients (Jeon et al., 2015). There are several reasons for this diffidence, such as different matrix, different controls, and different instrument analysis. Twenty-four discriminating metabolites between HCs and MMD patients were identified using GC–MS analysis and used to establish a biomarker panel using logistic regression. The AUC was ≥0.8 for all samples (Figure 4A), indicating a “good” classifier of MMD patients and HCs. Furthermore, the mechanism of MMD could be obtained by assessing the pathways underlying these biomarkers.
As seen in Figure 3B, these discriminating metabolites were involved in nine significantly different pathways related to the metabolism of amino acids, lipids, carbohydrates, and translation of carbohydrates.
Overall, 20% of the human body is composed of amino acids and their metabolites that are basic substrates and regulators in many metabolic pathways. Amino acid levels in patients with various diseases often differ from those of healthy individuals. Alterations in plasma concentrations of amino acids have been reported in liver fibrosis, non-small-cell lung cancer, aortic dissection, first-episode psychosis, and type-2 diabetes mellitus in patients with coronary artery disease (Yang et al., 2013; Leppik et al., 2018).
We found that levels of a panel of amino acids (L-alanine, L-isoleucine, L-valine, L-serine, L-threonine, L-methionine, L-phenylalanine, and L-tyrosine) were significantly different in MMD patients relative to HCs. Levels of L-alanine, L-valine, L-threonine, L-phenylalanine, and L-tyrosine were decreased in MMD patients, whereas L-isoleucine, L-serine, and L-methionine levels were increased. Such alterations of several amino acids suggest that amino acid metabolism is disturbed in MMD.
The synthesis of L-alanine is derived from pyruvate by alanine aminotransferase directly involving in gluconeogenesis and the alanine–glucose cycle and regulates glucose metabolism (Chen et al., 2017). Like gamma-aminobutyric acid, taurine, and glycine, L-alanine is an inhibitory neurotransmitter in the brain and is involved in lymphocyte reproduction and immunity. L-Alanine is disturbed in serum of MMD patients, but the association between L-alanine and MMD has not been investigated. Our study did not elucidate how altered levels of L-alanine affect MMD development; more studies are needed. L-Isoleucine and L-valine are branched chain amino acids (BCCAs). They are critical to human life and are particularly involved in stress, energy generation, and muscle metabolism (Zheng et al., 2017; Zhenyukh et al., 2017). The biosynthesis of L-isoleucine is initiated by the L-threonine deaminase reaction, whereas the pathway toward L-valine starts from pyruvate. Acetohydroxyacid synthase, ketol acid reductoisomerase (KARI), dihydroxyacid dehydratase (DHAD), and branched chain amino transferase (BCAT) were the four enzymes operating in L-isoleucine and L-valine biosynthesis (Galili et al., 2016). Abnormal changes in L-isoleucine and L-valine levels have been documented in first-episode psychosis (Leppik et al., 2018). Our study suggested that abnormal changes in L-isoleucine and L-valine levels are also associated with MMD.
Like glutamate, aspartate, and glycine, L-serine can act as an excitatory and inhibitory neurotransmitter (Yu et al., 2017; Martin et al., 2018). As seen in Figure 5, phosphoenolpyruvate (PEP) is one of the precursors to the shikimate pathway; in addition, PEP is also derived from oxaloacetate by phosphoenolpyruvate carboxykinase (PEPCK) in glycolysis; thus, they are linked to each other. L-Serine may be derived from the biosynthesis of the glycolytic intermediate 3-phosphoglycerate, which participates in cell proliferation and is necessary for specific functions in the central nervous system. Altered levels of L-serine in patients with psychiatric disorders underscore the amino acid’s importance in brain development and function. Thus, L-serine level alterations may be involved in MMD pathogenesis by affecting brain development and function. L-Threonine is an essential amino acid in humans, and severe deficiency causes neurological dysfunction (Liu et al., 2010). In our study, L-threonine level was lower in serum of MMD patients compared to HCs, indicating that L-threonine was altered in serum of MMD patients. Alterations in the ratio of L-tyrosine and L-phenylalanine levels could be a sign of compromised function of the dopaminergic system, which may be partly associated with MMD pathogenesis.
Collectively, our findings indicate that altered levels of amino acids could be linked to MMD. The precise mechanism by which amino acid levels influence the genesis and development of MMD should be investigated further. Pyroglutamate is converted to urea over glutamate and allantoin reactions; additionally, L-arginine could be converted to urea via urea cycle. Altered levels of pyroglutamate and urea were found in MMD. A panel of fatty acids (succinic acid, myristic acid, palmitoleic acid, palmitic acid, stearic acid, oleamide, octadecanamide, and pelargonic acid) could be used to distinguish between MMD patients and HCs. Succinic acid is an intermediate of the tricarboxylic acid (TCA) cycle (Bechthold et al., 2010; Lauritzen et al., 2014). Levels of succinic acid were significantly decreased in MMD patients compared with HCs. Lower levels of TCA-cycle intermediates may indicate alterations in the cycle. The TCA cycle is the core of cellular respiratory machinery and produces energy to power manufacture of compounds needed to defend against oxidative stress. Oxidation of fatty acids in mitochondria is responsible for approximately half of the total amount of adenosine triphosphate generated. Fatty acids also serve as the “building blocks” of cellular membranes after their esterification into phospholipids and are involved in signal transduction. Altered levels of fatty acids would lead to decreased energy production. However, whether such lower energy levels are associated with an increased risk of MMD is unknown.
The OPLS-DA method was used to analyze metabolic features and maximize discrimination between classes of compounds (Westerhuis et al., 2010). This approach reduced the effects of variability of non-relevant metabolites and helped identify serum metabolites contributing to differences between MMD patients and HCs. The reliability of OPLS-DA results was confirmed by correlation analyses and heatmaps.
Our study had two main limitations. First, MMD is an uncommon disease. We recruited 41 MMD patients, which may have introduced selection bias because of interindividual differences. Second, we only used a metabolomics approach, and the possibility that single-omic data restricted interpretation of our results cannot be excluded. Proteomics and genomics analyses are needed to confirm our findings.
In summary, we used a GC–MS platform to characterize the metabolic profiles of serum from MMD patients. Our analysis revealed important candidate metabolic biomarkers for MMD. Verification and validation studies with larger independent samples are necessary to demonstrate the utility of these metabolites as potential disease markers and to elucidate the pathophysiological mechanisms underlying MMD.
Data Availability Statement
The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.
Ethics Statement
Our work was carried out according to the Declaration of Helsinki. All participants provided written informed consent to allow use of their samples in these analyses, and our study protocol was reviewed and approved by the Medical Ethics Committee of the Jining First People’s Hospital (Jining, China) (No. 20170021).
Author Contributions
PJ and CC conception and design. CW, JZ, and WH acquisition of the data. CG and PJ analysis and interpretation of the data. CG drafting the manuscript. PJ and FJ critically revising the manuscript. CG, PJ, CC, and FJ reviewed submitted version of manuscript. PJ approved the final version of the manuscript on behalf of all authors. CG, DC, and YG statistical analysis. PJ study supervision.
Funding
Our work was supported by the National Natural Science Foundation of China (No. 81602846 to PJ), the Taishan Scholar Project of Shandong Province (No. tsqn201812159 to PJ), Natural Science Foundation of Shandong Province (No. ZR2018BH016 to CC), the Key Research and Development Program of Jining Science and Technology (No. 2018SMNS005 to FJ), and Teacher Research Support Fund of Jining Medical University (No. JYFC2018FKJ156 to CC).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Footnotes
References
Araki, Y., Yoshikawa, K., Okamoto, S., Sumitomo, M., Maruwaka, M., and Wakabayashi, T. (2010). Identification of novel biomarker candidates by proteomic analysis of cerebrospinal fluid from patients with moyamoya disease using SELDI-TOF-MS. BMC Neurol. 10:112. doi: 10.1186/1471-2377-10-112
Bechthold, I., Bretz, K., Kabasci, S., Kopitzky, R., and Springer, A. (2010). Succinic acid: anew platform chemical for biobased polymers from renewable resources. Chem. Eng. Technol. 31, 647–654. doi: 10.1002/ceat.200800063
Cai, H. L., Li, H. D., Yan, X. Z., Sun, B., Zhang, Q., Yan, M., et al. (2012). Metabolomic analysis of biochemical changes in the plasma and urine of first-episode neuroleptic-naive schizophrenia patients after treatment with risperidone. J. Proteome Res. 11, 4338–4350. doi: 10.1021/pr300459d
Chen, H. Y., Fang, S., and Wang, L. (2017). Interactions of 1-butyl-2,3-dimethylimidazolium bromide ionic liquid with glycine, l -alanine and l -valine: a volumetric and NMR spectroscopic study. J. Mol. Liq. 225, 706–712. doi: 10.1016/j.molliq.2016.10.125
Chen, X., Liu, L., Palacios, G., Gao, J., Zhang, N., Li, G., et al. (2015). Plasma metabolomics reveals biomarkers of the atherosclerosis. J. Sci. 2776–2783. doi: 10.1002/jssc.201000395
Do, K. T., Wahl, S., Raffler, J., Molnos, S., Laimighofer, M., Adamski, J., et al. (2018). Characterization of missing values in untargeted MS-based metabolomics data and evaluation of missing data handling strategies. Metabolomics 14:128. doi: 10.1007/s11306-018-1420-2
Galili, G., Amir, R., and Fernie, A. R. (2016). The regulation of essential amino acid synthesis and accumulation in plants. Annu. Rev. Plant Biol. 67, 153–178. doi: 10.1146/annurev-arplant-043015-112213
Huang, S., Guo, Z. N., Shi, M., Yang, Y., and Rao, M. (2017). Etiology and pathogenesis of moyamoya disease: an update on disease prevalence. Int. J. Stroke 12, 246–253. doi: 10.1177/1747493017694393
Jeon, J. P., Yun, T., Jin, X., Cho, W. S., Son, Y. J., Bang, J. S., et al. (2015). 1H-NMR-based metabolomic analysis of cerebrospinal fluid from adult bilateral moyamoya disease: comparison with unilateral moyamoya disease and atherosclerotic stenosis. Medicine 94:e629. doi: 10.1097/MD.0000000000000629
Jung, J. Y., Lee, H. S., Kang, D. G., Kim, N. S., Cha, M. H., Bang, O. S., et al. (2011). 1H-NMR-based metabolomics study of cerebral infarction. Stroke 42, 1282–1288. doi: 10.1161/STROKEAHA.110.598789
Kumar, N., Hoque, M. A., Shahjaman, M., Islam, S. M., and Mollah, M. N. (2017). Metabolomic biomarker identification in presence of outliers and missing values. Biomed. Res. Int. 2017:2437608. doi: 10.1155/2017/2437608
Kuroda, S., and Houkin, K. (2008). Moyamoya disease: current concepts and future perspectives. Lancet Neurol. 7, 1056–1066. doi: 10.1016/S1474-4422(08)70240-0
Lauritzen, I., Blondeau, N., Heurteaux, C., Widmann, C., Romey, G., and Lazdunski, M. (2014). Polyunsaturated fatty acids are potent neuroprotectors. EMBO J. 19, 1784–1793. doi: 10.1093/emboj/19.8.1784
Lee, J. H., Youn, T. J., Yoon, Y. E., Park, J. J., Hong, S. J., Chun, E. J., et al. (2013). Coronary artery stenosis in moyamoya disease: tissue characterization by 256-slice multi-detector CT and virtual histology. Circulation 127, 2063–2065. doi: 10.1161/CIRCULATIONAHA.112.136473
Leppik, L., Kriisa, K., Koido, K., Koch, K., Kajalaid, K., Haring, L., et al. (2018). Profiling of amino acids and their derivatives biogenic amines before and after antipsychotic treatment in first-episode psychosis. Front. Psychiatry 9:155. doi: 10.3389/fpsyt.2018.00155
Liu, J. Q., Tohru, D., Nobuya, I., Michihiko, K., Sakayu, S., and Hideaki, Y. (2010). Gene cloning, biochemical characterization and physiological role of a thermostable low-specificity L-threonine aldolase from Escherichia coli. Eur. J. Biochem. 255, 220–226. doi: 10.1046/j.1432-1327.1998.2550220.x
Martin, S. A., Souder, D. C., Miller, K. N., Clark, J. P., Sagar, A. K., Eliceiri, K. W., et al. (2018). GSK3β regulates brain energy metabolism. Cell Rep. 23, 1922–1931. doi: 10.1016/j.celrep.2018.04.045
Maruwaka, M., Yoshikawa, K., Okamoto, S., Araki, Y., Sumitomo, M., Kawamura, A., et al. (2015). Biomarker research for moyamoya disease in cerebrospinal fluid using surface-enhanced laser desorption/ionization time-of-flight mass spectrometry. J. Stroke Cerebrovasc. Dis. 24, 104–111. doi: 10.1016/j.jstrokecerebrovasdis.2014.07.028
Okazaki, S., Morimoto, T., Kamatani, Y., Kamimura, T., Kobayashi, H., Harada, K., et al. (2019). Moyamoya disease susceptibility variant RNF213 p.R4810K increases the risk of ischemic stroke attributable to large-artery atherosclerosis. Circulation 139, 295–298. doi: 10.1161/CIRCULATIONAHA.118.038439
Pandey, R., Caflisch, L., Lodi, A., Brenner, A. J., and Tiziani, S. (2017). Metabolomic signature of brain cancer. Mol. Carcinog. 56, 2355–2371. doi: 10.1002/mc.22694
Sato, Y., Kazumata, K., Nakatani, E., Houkin, K., and Kanatani, Y. (2019). Characteristics of moyamoya disease based on national registry data in Japan. Stroke 50, 1973–1980. doi: 10.1161/STROKEAHA.119.024689
Scott, R. M., and Smith, E. R. (2009). Moyamoya disease and moyamoya syndrome. N. Engl. J. Med. 360, 1226–1237. doi: 10.1056/NEJMra0804622
Turer, A., Stevens, R. Jr., Muehlbauer, M., Van-Der-Westhuizen, J., Mathew, J., Schwinn, D., et al. (2009). Metabolomic profiling reveals distinct patterns of myocardial substrate use in humans with coronary artery disease or left ventricular dysfunction during surgical ischemia/reperfusion. Circulation 119, 1736–1746. doi: 10.1161/CIRCULATIONAHA.108.816116
Westerhuis, J. A., van Velzen, E. J., Hoefsloot, H. C., and Smilde, A. K. (2010). Multivariate paired data analysis: multilevel PLSDA versus OPLSDA. Metabolomics 6, 119–128. doi: 10.1007/s11306-009-0185-z
Yang, J., Chen, T., Sun, L., Zhao, Z., Qi, X., Zhou, K., et al. (2013). Potential metabolite markers of schizophrenia. Mol. Psychiatry 18, 67–78. doi: 10.1038/mp.2011.131
Yu, J., Qin, B., Wu, F., Qin, S., Nowsheen, S., Shan, S., et al. (2017). Regulation of serine-threonine kinase akt activation by NAD(+)-dependent deacetylase SIRT7. Cell Rep. 18, 1229–1240. doi: 10.1016/j.celrep.2017.01.009
Zheng, L., Wei, H., He, P., Zhao, S., Xiang, Q., Pang, J., et al. (2017). Effects of supplementation of branched-chain amino acids to reduced-protein diet on skeletal muscle protein synthesis and degradation in the fed and fasted states in a piglet model. Nutrients 9:17. doi: 10.3390/nu9010017
Zhenyukh, O., Civantos, E., Ruizortega, M., Sánchez, M. S., Vázquez, C., Peiró, C., et al. (2017). High concentration of branched-chain amino acids promotes oxidative stress, inflammation and migration of human peripheral blood mononuclear cells via mTORC1 activation. Free Radic Biol. Med. 104, 165–177. doi: 10.1016/j.freeradbiomed.2017.01.009
Keywords: moyamoya disease, gas chromatography–mass spectrometry, biomarkers, serum, metabolomics
Citation: Geng C, Cui C, Guo Y, Wang C, Zhang J, Han W, Jin F, Chen D and Jiang P (2020) Metabolomic Profiling Revealed Potential Biomarkers in Patients With Moyamoya Disease. Front. Neurosci. 14:308. doi: 10.3389/fnins.2020.00308
Received: 12 November 2019; Accepted: 16 March 2020;
Published: 21 April 2020.
Edited by:
Antonio Paoli, University of Padua, ItalyReviewed by:
Shayne Mason, North-West University, South AfricaSung Ung Kang, Johns Hopkins University, United States
Copyright © 2020 Geng, Cui, Guo, Wang, Zhang, Han, Jin, Chen and Jiang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Feng Jin, amluZmVuZ3Nkam5AMTYzLmNvbQ==; Pei Jiang, amlhbmdwZWljc3VAc2luYS5jb20=
†These authors have contributed equally to this work