ORIGINAL RESEARCH article

Front. Mol. Neurosci., 30 June 2020

Sec. Brain Disease Mechanisms

Volume 13 - 2020 | https://doi.org/10.3389/fnmol.2020.00080

Integrated Metabolomics and Proteomics Analysis Reveals Plasma Lipid Metabolic Disturbance in Patients With Parkinson’s Disease

  • LH

    Ling Hu 1,2

  • MD

    Mei-Xue Dong 2

  • YH

    Yan-Ling Huang 3

  • CL

    Chang-Qi Lu 1

  • QQ

    Qian Qian 1

  • CZ

    Chun-Cheng Zhang 4

  • XX

    Xiao-Min Xu 5

  • YL

    Yang Liu 1

  • GC

    Guang-Hui Chen 2

  • YW

    You-Dong Wei 1*

  • 1. Department of Neurology, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China

  • 2. Department of Neurology, Renmin Hospital of Wuhan University, Hubei General Hospital, Wuhan, China

  • 3. Department of Neurology, Chongqing University Central Hospital, Chongqing Emergency Medical Center, Chongqing, China

  • 4. Department of Neurology, The People’s Hospital of Tongliang District, Chongqing, China

  • 5. Department of Neurology, The Affiliated Hospital of Southwest Medical University, Luzhou, China

Article metrics

View details

72

Citations

9,8k

Views

3,5k

Downloads

Abstract

Parkinson’s disease (PD) is a common neurodegenerative disease in the elderly with a pathogenesis that remains unclear. We aimed to explore its pathogenesis through plasma integrated metabolomics and proteomics analysis. The clinical data of consecutively recruited PD patients and healthy controls were assessed. Fasting plasma samples were obtained and analyzed using metabolomics and proteomics methods. After that, differentially expressed metabolites and proteins were identified for further bioinformatics analysis. No significant difference was found in the clinical data between these two groups. Eighty-three metabolites were differentially expressed in PD patients identified by metabolomics analysis. These metabolites were predominately lipid and lipid-like molecules (63%), among which 25% were sphingolipids. The sphingolipid metabolism pathway was enriched and tended to be activated in the following KEGG pathway analysis. According to the proteomics analysis, forty proteins were identified to be differentially expressed, seven of which were apolipoproteins. Furthermore, five of the six top ranking Gene Ontology terms from cellular components and eleven of the other fourteen Gene Ontology terms from biological processes were directly associated with lipid metabolism. In KEGG pathway analysis, the five enriched pathways were also significantly related with lipid metabolism (p < 0.05). Overall, Parkinson’s disease is associated with plasma lipid metabolic disturbance, including an activated sphingolipid metabolism and decreased apolipoproteins.

Introduction

Parkinson’s disease (PD) is a common neurodegenerative disease and is expected to affect 14.2 million people worldwide by 2040 (Dorsey and Bloem, 2018). Its typical clinical features are motor symptoms, which result from accumulated Lewy bodies and the remarkable death of dopaminergic neurons in the substantia nigra pars compacta (SNpc). In addition, Parkinson’s disease is correlated with many non-motor symptoms, some of which even precede the motor symptoms by more than 10 years, including constipation, hyposmia, anxiety, depression, and rapid eye movement sleep behavior disorder (Emamzadeh and Surguchov, 2018). These non-motor symptoms are probably due to the disturbance of various neurotransmitters and multiple nerve systems in addition to the dopaminergic neurons in the SNpc (Lv and Ae, 2015). Most idiopathic PD patients are sporadic, and various environmental exposures have been identified as being correlated with the development of PD. Results of genome-wide association studies have also identified 90 independent mutations in more than 20 genes that increase risk factors for Parkinson’s disease (Blauwendraat et al., 2019). Therein, the strongest genetic risk factor is the Asn370Ser mutation of β-glucocerebrosidase, with an odds ratio greater than 5 (Sidransky et al., 2009). The interactions between environmental and genetic risk factors are still under investigation.

Drugs that increase intracerebral dopamine levels or directly stimulate dopamine receptors remain as the most common forms of treatment for PD patients (Connolly and Lang, 2014). These drugs are mainly for symptomatic treatment without neuroprotection or a so-called disease-modifying function. Researchers have found that mitochondrial oxidative stress and inflammation reaction play important roles in the apoptosis of dopaminergic neurons during the occurrence and development of Parkinson’s disease (Wei et al., 2018). However, drugs targeting oxidative stress have exhibited little therapeutic benefit (Trist et al., 2019). Research also indicates that Parkinson’s disease is related to the dysregulation of protein homeostasis, including intracellular and membrane protein trafficking, protein aggregation, and protein degradation by the ubiquitin-proteasome or lysosome-autophagy systems (Deng et al., 2018). Thus, the potential pathogenesis should be clarified and more effective disease-modifying treatments are urgently needed for PD patients.

Metabolomics is a comprehensive assessment of the total endogenous metabolites in a biological system, and proteomics can evaluate alterations at the protein level (Mostafa et al., 2016). These omics data can help researchers make many discoveries in the alteration of molecules and metabolic pathways following gene expression in PD patients. Plasma is an easily obtainable, non-invasive, and informative biofluid from patients, making it perfect to explore the pathogenesis of many neuropsychiatric disorders (Hu et al., 2016). Herein, we employed untargeted metabolomics and proteomics analysis to identify molecular changes in PD patients. Furthermore, an integrated bioinformatic analysis was performed to explore the probable pathogenesis of PD based on the above findings (Dong et al., 2018a).

Materials and Methods

Objects

Parkinson’s disease patients were consecutively recruited from April 2016 to February 2017 in the Department of Neurology, the First Affiliated Hospital of Chongqing Medical University, according to the criteria reported in a previous research (Dong et al., 2017b). The inclusion criteria was: (i) PD patients diagnosed using the European Federation of Neurological Societies and the International Parkinson Movement Disorder Society’s European Section (Berardelli et al., 2013), and (ii) patients only taking dopamine analogs or dopamine receptor agonists. The exclusion criteria was: (i) those with secondary Parkinson’s disease or Parkinson-plus syndrome; (ii) patients suffering from tumor, heart failure, chronic obstructive pulmonary disease, nephritis, infectious diseases, or any other severe chronic disease at the time of enrollment; and (iii) any history of stroke, brain surgery, head trauma, motor neuro disease, Alzheimer’s disease, or mental diseases.

A same amount of age- and sex-matched healthy controls (HCs) were included from the Department of Physical Examination. The controls were also without any history of illness in the central nervous system or suffering from any other severe disease. This study was approved by the ethics committee of the First Affiliated Hospital of Chongqing Medical University and performed in accordance with the Declaration of Helsinki. Informed consent was obtained from all individual participants included in the study. Clinical data, metabolomics analysis, and proteomics analysis were blindly assessed or performed (Dong et al., 2017a).

Clinical Data

The clinical data of all included participants were collected. Fasting plasma samples were obtained with an EDTA-K2 tube in the morning, and then stored at -80°C until experimental analysis. All the clinical scales were assessed by experienced neurologists. Statistical analyses were performed using Statistic Package for the Social Sciences 22.0 (IBM, Armonk, NY, United States). Categorical data are exhibited as absolute numbers while continuous data are exhibited as mean ± standard error. All the clinical data were compared between PD and HC groups using Pearson Chi-squared tests or Fisher exact tests for categorical data and Mann–Whitney U-tests with Bonferroni post hoc tests for continuous data, as appropriate (Dong et al., 2016).

Metabolomics Analysis

We adopted a Waters UPLC I-class system equipped with a binary solvent delivery manager (Waters Corporation, Milford, United States) to perform the untargeted liquid chromatography-mass spectrometry-based metabolomics (thirty-six objects per group); the detailed procedure was described in a previous research (Dong et al., 2018a). After that, data sets including m/z, peak retention time, and peak intensity of each ion were obtained. The m/z-peak retention time pairs were used to identify each ion based on Metlin1 and Human Metabolome Database (HMDB)2 while the peak intensity was deemed as the level of a metabolite. The data sets were further reduced by removing any peaks with a missing value in more than 60% of the total samples.

The positive and negative peak data were merged and multivariate statistical analyses were performed by the SIMCA-P 13.0 software package (Umetrics, Umeå, Sweden). The partial least squares-discriminant analysis (PLS-DA) model was constructed to show statistical differences and identify metabolites differentially expressed between PD patients and healthy controls, and this model was further validated by a permutation test with 200 iterations. Metabolites with variable influence on projection values (obtained from the PLS-DA model) of greater than 1.0, fold change values of greater than ±2, and p-values of less than 0.05 (Dong et al., 2018b) were identified to be differentially expressed. These differentially expressed metabolites were then classified by chemical taxonomy based on the HMDB database, and metabolic pathway analysis was further performed by MetaboAnalyst 4.03 (Chong et al., 2018).

Proteomics Analysis

The plasma of twenty-seven randomized healthy controls and thirty randomized PD patients were merged into three pooled samples for proteomics analysis. The details of the proteomics analysis have also been depicted previously (Dong et al., 2018a). Afterward, tandem mass spectrometry spectra were obtained and searched using the MASCOT engine 2.2 (Matrix Science, London, United Kingdom). The differentially expressed proteins were recognized by fold change values of greater than ±1.2 and p-values of less than 0.05. Gene ontology (GO) enrichment from cellular component, molecular function, and biological process, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed by Fisher’s exact test, considering the whole quantified protein annotations as the background data set. The obtained p-values from enrichment analysis was further converted into q-value using Benjamini-Hochberg correction and a q-value less than 0.05 was considered significant (Hochberg and Benjamini, 1990).

Results

Results of the Metabolomics Analysis

Thirty-six PD patients and healthy controls were separately included in the metabolomics analysis. The clinical data of these participants are shown in Table 1. No significant differences were observed between these two groups, including the levels of hemoglobin A1C (HbA1C) and blood lipid, or the incidences of diabetes mellitus and hypercholesterolemia. The mean UPDRS score and Hoehn–Yahr score of these PD patients were 41.11 ± 3.72 and 2.26 ± 0.17, respectively.

TABLE 1

Variable (SEM/%)HC (36)PD (36)p-valueVariable (SEM/%)HC (36)PD (36)p-value
Age (year)62.16 ± 1.7364.03 ± 0.950.348HbA1c (%)6.17 ± 0.206.08 ± 0.140.689
Gender, Male (%)18 (50%)18 (50%)1TC (mmol/L)4.42 ± 0.124.53 ± 0.130.516
Smoking history (%)10 (27.8%)9 (25.0%)0.789TG (mmol/L)1.53 ± 0.091.44 ± 0.130.6
Alcohol consumption (%)2 (5.6%)2 (5.6%)1HDL-C (mmol/L)1.30 ± 0.061.33 ± 0.060.794
Hypertension (%)15 (46.9%)17 (53.1%)0.617LDL-C (mmol/L)2.87 ± 0.122.91 ± 0.110.816
Diabetes mellitus (%)4 (12.5%)8 (25.0%)0.2Apo-A1 (g/L)1.37 ± 0.041.34 ± 0.040.658
Hypercholesterolemia (%)12 (37.5%)9 (28.1%)0.424Apo-B (g/L)0.91 ± 0.040.91 ± 0.040.992
BMI (kg/m2)24.21 ± 0.5723.00 ± 0.640.169Lpa (mg/L)190.76 ± 42.32279.53 ± 57.030.216
UPDRS score/41.11 ± 3.72/Hoehn-Yahr score/2.26 ± 0.17/

Clinical data of PD patients and healthy controls included in the untargeted liquid chromatography-mass spectrometry-based metabolomics analysis.

PD, Parkinson’s disease; SEM, standard error of the mean; HC, healthy control; HbA1c, hemoglobin A1c; TC, total cholesterol; TG, triglyceride; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; Apo-A1, apolipoprotein A1; Apo-B, apolipoprotein B; BMI, body mass index; Lpa, lipoprotein a; UPDRS, Unified Parkinson’s Disease Rating Scale.

After excluding internal standards, 6040 positive and 4363 negative peaks were detected in almost 98.81% of samples on average (Supplementary Table S1). Based on the above peaks, the PLS-DA score plot showed a significant difference between the PD and HC groups (R2X = 0.132, R2Y = 0.374, Q2 = 0.2) (Figure 1A). A permutation test with 200 iterations was performed and confirmed that the PLS-DA model was valid and not over-fitted because the original right R2 and Q2 values were observably higher than the corresponding permutated left values [R2 = (0.0, 0.746), Q2 = (0.0, -0.198)] (Figure 1B). According to the above results, eighty-three metabolites were differentially expressed between these two groups (Table 2). Of these, 63% were lipid and lipid-like molecules, 13% were organic acids and derivatives, 6% were phenylpropanoids and polyketides, 5% were organic oxygen compounds, 5% were organoheterocyclic compounds, 4% were benzenoids, 2% were nucleosides, nucleotides, and analogs, 1% was homogeneous non-metal compounds, and 1% was organosulfur compounds (Figure 1C). The lipid and lipid-like molecules were further categorized into sphingolipids (25%), glycerophospholipids (16%), glycerolipids (17%), fatty acyls (15%), prenol lipids (15%), and steroids and steroid derivatives (12%) (Figure 1D). These sphingolipids were ganglioside GD3 (d18:0/18:0), PS[15:0/20:4(8Z,11Z,14Z,17Z)], camellianin D, ceramide (d18:1/16:0), adrenorphin, ganglioside GD1a (d18:0/22:0), SM(d18:0/24:1(15Z)(OH)), thiamine(1+) diphosphate(1-), N-docosahexaenoyl phenylalanine, ganglioside GA2 (d18:1/12:0), SM(d18:0/22:0), glucosylceramide (d18:1/18:0), and AS 1–5.

TABLE 2

Compound NameCompound IDam/z valueRT (min)Ion modeFC valuebcVIP valuec-log(p-value)c
Coenzyme AHMDB0001423750.01177.6840Positive0.43232.46855.292
Ganglioside GD3 (d18:0/18:0)HMDB0011859737.84482.9037Positive2.49472.17126.520
PS(15:0/20:4(8Z,11Z,14Z,17Z))HMDB0012312752.48997.6840Positive0.41532.58905.192
4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanolHMDB0030032419.24013.8761Positive0.44342.00362.761
LupuloneHMDB0030041829.57192.8968Positive2.17492.44018.216
S-(Hydroxymethyl)glutathioneHMDB0004662113.37488.8773Positive2.52481.53403.835
Glycinoprenol 11HMDB0002003737.73292.8968Positive2.45382.27827.317
TG(24:0/15:0/o-18:0)HMDB47052919.92532.9037Positive2.11371.84445.140
Camellianin DHMDB0011920620.41297.7665Positive0.48532.52315.238
Dipiperamide CHMDB0034617539.26661.9177Positive3.86761.14031.496
DG(24:0/24:0/0:0)HMDB0007818758.04659.2399Positive2.00501.00101.418
IodideHMDB001223865.45238.6504Positive2.47511.42821.689
TG[24:1(15Z)/20:5(5Z,8Z,11Z,14Z,17Z)/ 18:4(6Z,9Z,12Z,15Z)]HMDB52343947.79352.9037Positive2.25862.30587.459
PC[DiMe(11,3)/DiMe(13,5)]HMDB0013456922.660410.4327Positive2.22551.15331.568
Lucidenic acid MHMDB0035973925.59139.7743Positive3.06681.68013.860
Pangamic acidHMDB0029949900.57749.5886Positive2.12751.45463.466
Ceramide (d18:1/16:0)HMDB0004949520.50858.7673Positive2.67571.83234.813
TyramineHMDB000030669.57518.6504Positive2.17231.40841.849
AdrenorphinHMDB0012081948.51082.9037Positive2.57862.49558.273
Phenylacetyl-CoAHMDB0006503886.13369.9806Positive2.43281.22232.087
TG[22:1(13Z)/22:2(13Z,16Z)/o-18:0]HMDB51764947.93732.8968Positive2.18082.42748.070
Ganglioside GD1a (d18:0/22:0)HMDB0011784948.08132.9037Positive2.24652.41187.750
KyotorphinHMDB0005768675.35258.1291Positive2.27651.45633.057
KerdlanHMDB0240277853.662410.4396Positive2.09731.17771.769
Threoninyl-GlutamineHMDB0029059495.24111.8902Positive4.78861.19711.612
TG[18:4(6Z,9Z,12Z,15Z)/22:6(4Z,7Z,10Z,13Z,16Z,19Z)/ 18:4(6Z,9Z,12Z,15Z)]HMDB0035640883.62349.9875Positive2.55441.24632.114
SM[d18:0/24:1(15Z)(OH)]HMDB0013469829.68159.7674Positive2.50751.51063.333
3-Methylglutaconyl-CoAHMDB0001057858.09959.5886Positive2.43451.19682.040
Dynorphin A (6-8)HMDB0012932887.62669.9737Positive2.16181.33472.575
Thiamine(1+) Diphosphate(1−)HMDB0011854875.99816.9452Positive2.08551.12722.154
LysoPC[16:1(9Z)/0:0]HMDB0010383987.61229.7674Positive2.81401.42392.992
DihydroceramideHMDB0006752686.604310.0268Positive2.68271.35792.580
D-UrobilinogenHMDB0004158591.31613.4380Positive2.43181.16141.396
N-Docosahexaenoyl phenylalanineHMDB0006482951.666810.4396Positive2.20631.35912.418
Alanyl-AsparagineHMDB0028682407.18861.8211Positive4.30231.26421.676
MG[0:0/22:4(7Z,10Z,13Z,16Z)/0:0]HMDB0011554840.64089.9875Positive2.74931.32652.430
IopromideHMDB0001493814.12019.5886Positive2.34651.33922.809
TaraxacosideHMDB0030055797.258310.1989Positive2.36521.57633.144
Elaidic carnitineHMDB0006464851.701010.5890Positive2.04741.09891.502
10-Hydroxy-2,8-decadiene-4,6-diynoic acidHMDB0031054177.05441.6216Positive0.00012.23884.311
PE-NMe[18:3(6Z,9Z,12Z)/22:6(4Z,7Z,10Z,13Z,16Z,19Z)]HMDB0009281400.77899.4806Positive2.10621.79564.829
Mandelic acidHMDB0000703153.05461.6216Positive0.00782.36864.793
Sesaminol 2-O-triglucosideHMDB0029556831.249410.1989Positive2.25061.63073.334
Polyporusterone CHMDB0000138953.622210.1989Positive2.02351.54233.039
Ganglioside GA2 (d18:1/12:0)HMDB0004888991.60739.7537Positive2.09881.32692.273
GlucoconringiinHMDB0001197810.15249.9806Positive2.17481.40052.739
Mannosyl-(N-acetylglucosaminyl)2-diphosphodolicholHMDB0012255919.63189.7605Positive2.70271.61243.791
3-phenylprop-2-enoic acidHMDB0000567149.05971.6216Positive0.00472.36434.806
SM(d18:0/22:0)HMDB0012091789.681710.1989Positive2.31981.22482.279
Alpha-linolenyl carnitineHMDB0006319843.64559.9806Positive2.22401.41522.813
Perulactone BHMDB0030119977.605010.1989Positive2.81901.52222.913
PE[DiMe(11,3)/DiMe(13,5)]HMDB0008682880.59949.8293Positive2.08441.70574.208
EbastineHMDB0035995939.653510.4327Positive2.23301.21811.702
TG[18:3(6Z,9Z,12Z)/18:4(6Z,9Z,12Z,15Z)/18:3(6Z,9Z,12Z)]HMDB0010470871.677410.4396Positive2.11821.21931.698
7-chloro-2-(3,4-dimethoxyphenyl)-3,5,6-trihydroxy-8-methoxy-4H-chromen-4-oneHMDB0001484816.11389.5955Positive2.02901.24752.113
(3beta,5alpha,9alpha,22E,24R)-5,9-Epidioxy-3-hydroxyergosta-7,22-dien-6-oneHMDB0032666885.632610.1989Positive2.18761.56693.119
Cholesterol sulfateHMDB0000653931.623810.2009Negative2.37741.58233.259
Methylimidazole acetaldehydeHMDB0004181371.18293.3773Negative0.45042.13792.501
Vanillactic acidHMDB0000913211.06081.8145Negative0.00272.42794.091
PC[14:1(9Z)/22:2(13Z,16Z)]HMDB0007921764.553510.0221Negative2.05901.36542.573
2-Octenedioic acidHMDB0000341343.14141.7113Negative0.02301.47652.370
1,2,3,4-Tetrahydro-beta-carbolineHMDB0012488515.28763.2529Negative2.01701.19221.332
MethyldopaHMDB0011754210.07681.6219Negative0.00002.57994.767
PE[20:2(11Z,14Z)/24:1(15Z)]HMDB0009311834.62329.8316Negative2.19001.83225.080
Glucosylceramide (d18:1/18:0)HMDB0004972708.615910.0221Negative2.98641.36082.665
cis-Hydroxy PerhexilineHMDB60644585.53059.2354Negative2.79831.37342.780
TG[18:4(6Z,9Z,12Z,15Z)/20:5(5Z,8Z,11Z,14Z,17Z)/ 18:4(6Z,9Z,12Z,15Z)]HMDB55513873.65699.7697Negative2.84001.67944.108
Cinncassiol D1 glucosideHMDB0034677513.27303.3567Negative2.47471.63982.277
AS 1-5HMDB0032843714.55358.7898Negative2.19151.45752.933
Glycochenodeoxycholic acid 3-glucuronideHMDB0002579606.32943.1567Negative2.56221.62892.650
TG[14:1(9Z)/14:0/14:1(9Z)]HMDB47724717.59259.4073Negative2.10381.33652.913
Asparaginyl-SerineHMDB0028740437.16283.6093Negative0.46501.80951.691
1-non-adecanoyl-glycero-3-phosphateHMDB62322903.59609.7766Negative2.61351.66103.764
Hexahydro-6,7-dihydroxy-5-(hydroxymethyl)-3-(2-hydroxyphenyl)-2H-pyrano[2,3-d]oxazol-2-oneHMDB0029234278.06481.6219Negative0.00012.57554.609
PE[DiMe(13,5)/MonoMe(13,5)]HMDB61496902.61039.8316Negative3.50421.66084.153
(2-hydroxy-2-{9-[(3-methylbut-2-enoyl)oxy]-2-oxo-2H,8H,9H-furo[2,3-h]chromen-8-yl}propoxy)sulfonic acidHMDB0001511421.05962.7668Negative0.44441.65911.564
(2S)-2-amino-3-[4-hydroxy-3-(sulfooxy)phenyl]-2-methylpropanoic acidHMDB0142153290.03421.6082Negative0.12252.11242.902
AsiaticosideHMDB36656979.59189.9829Negative2.27561.72614.518
Ponasteroside AHMDB0034091625.36004.2329Negative2.25311.08981.442
MG(12:0/0:0/0:0)HMDB72863821.63239.9829negative2.20341.63304.102
7-Chloro-6-demethylcepharadione BHMDB0031833340.03551.6219Negative0.00012.32573.726
NorpropoxypheneHMDB0011627974.59918.7760Negative2.39741.14481.437
Hordatine B glucosideHMDB0030460370.17824.0679Negative0.39561.71691.856

Key differentially expressed metabolites identified by untargeted liquid chromatography-mass spectrometry-based metabolomics analysis between PD patients and healthy control.

aCompound ID was mainly exhibited based on the Human Metabolome Database (www.hmdb.ca). bFC value was calculated as the ratio of the average mass response (area) between the two groups (FC value = PD/HC). Thus, FC values >1 indicate significantly higher levels in the PD group relative to the HC group while FC values <1 indicate significantly lower levels in the PD group. cOnly metabolites with FC values greater than ±2.0, VIP values greater than 1.0 and p-values less than 0.05 were deemed statistically significant. PD, Parkinson’s disease; RT, retention time; FC, fold change; VIP, variable influence on projection; HC, healthy control.

FIGURE 1

The related metabolic pathway is shown in Figure 2A, where only the sphingolipid metabolism pathway is significantly enriched (pathway impact = 0.473, q-value = 0.004). All six metabolites involved in sphingolipid metabolism had significantly increased, indicating that the pathway is activated in PD patients (Figure 2B).

FIGURE 2

Results of the Proteomics Analysis

The clinical data of all the participants included in the proteomics analysis are exhibited in Table 3. No statistical differences can be observed in these data, indicating that the following comparison of the two groups was reasonable. A total of 912 proteins were identified in plasma and protein ratio distributions between the two groups, which are shown in Figure 3A. Forty differentially expressed proteins were recognized, with fifteen proteins increased and twenty-five decreased (Figure 3B). The details of these proteins were shown in Table 4 and the seven differentially expressed apolipoproteins (apolipoprotein C-I, apolipoprotein C-III, protein APOC4-APOC2, apolipoprotein C-IV, apolipoprotein B variant, apolipoprotein B, and apolipoprotein M) were all significantly decreased.

TABLE 3

Variable (SEM/%)HC (27)PD (30)p-valueVariable (SEM/%)HC (27)PD (30)p-value
Age (year)66.59 ± 1.1469.17 ± 1.500.184HbA1c (%)5.99 ± 0.186.17 ± 0.190.515
Gender, Male (%)14 (51.9%)19 (63.3%)0.381TC (mmol/L)4.33 ± 0.174.32 ± 0.160.952
Smoking history (%)9 (33.3%)7 (23.3%)0.402TG (mmol/L)1.49 ± 0.151.23 ± 0.140.206
Alcohol consumption (%)3 (11.1%)1 (3.3%)0.53HDL-C (mmol/L)1.21 ± 0.081.24 ± 0.090.858
Hypertension (%)15 (55.6%)13 (43.3%)0.357LDL-C (mmol/L)2.81 ± 0.152.73 ± 0.130.679
Diabetes mellitus (%)2 (7.4%)5 (16.7%)0.51Apo-A1 (g/L)1.31 ± 0.051.30 ± 0.050.93
Hypercholesterolemia (%)4 (14.8%)6 (20.0%)0.869Apo-B (g/L)0.89 ± 0.040.85 ± 0.040.403
BMI (kg/m2)23.95 ± 0.6922.79 ± 0.580.197Lpa (mg/L)201.19 ± 48.58270.04 ± 66.490.407
UPDRS score/38.49 ± 3.98/Hoehn-Yahr score/2.40 ± 0.16/

Clinical data of PD patients and healthy controls included in the tandem mass tag-based proteomics analysis.

PD, Parkinson’s disease; SEM, standard error of the mean; HC, healthy control; HbA1c, hemoglobin A1c; TC, total cholesterol; TG, triglyceride; HDL-C, high-density lipoprotein cholesterol; LDL-C, low-density lipoprotein cholesterol; Apo-A1, apolipoprotein A1; Apo-B, apolipoprotein B; BMI, body mass index; Lpa, lipoprotein a; UPDRS, Unified Parkinson’s Disease Rating Scale.

TABLE 4

Uniprot IDProtein name (GeneName)MW [kDa]FC valueap-valueb
O75339Cartilage intermediate layer protein 1 (CILP)132.481.2470.003
A9UFC0Caspase 14 (CASP14)27.652.0770.010
Q5T619Zinc finger protein 648 (ZNF648)62.301.2760.010
Q04756Hepatocyte growth factor activator (HGFAC)70.641.6050.013
Q9UQ05Potassium voltage-gated channel subfamily H member 4 (KCNH4)111.621.2130.014
E5RIF9Carbonic anhydrase 1 (CA1)16.311.2300.021
P01023Alpha-2-macroglobulin (A2M)163.191.2510.025
B7Z544cDNA FLJ51742, highly similar to Inter-alpha-trypsin inhibitor heavy chain H498.291.2060.026
Q6UWP8Suprabasin (SBSN)60.501.2300.028
A0A087WYF1Laminin subunit alpha-2 (LAMA2)343.201.2540.030
P04275von Willebrand factor (VWF)309.061.2100.031
D3DQX7Serum amyloid A protein (SAA1)13.551.2640.032
A0A125U0U7MS-C1 heavy chain variable region13.091.2100.039
Q96K23cDNA FLJ14838 fis, clone OVARC1001726, weakly similar to APICAL-LIKE PROTEIN27.991.3800.041
D6RF20Vitamin D-binding protein(GC)16.071.2550.046
P02656Apolipoprotein C-III (APOC3)10.850.7180.001
Q1W658Follicle-stimulating hormone beta subunit (FSHB)8.450.3310.002
K7ER74Protein APOC4-APOC2 (APOC4-APOC2)20.040.6610.002
Q04695Keratin, type I cytoskeletal 17 (KRT17)48.080.4650.003
Q5NV68V4-1 protein (V4-1)11.230.3420.005
Q59HB3Apolipoprotein B variant183.460.8320.007
P04259Keratin, type II cytoskeletal 6B (KRT6B)60.030.6110.008
K7ERI9Apolipoprotein C-I (APOC1)8.640.7100.009
G5E968Chromogranin A (CHGA)34.250.7860.009
Q99592Zinc finger and BTB domain-containing protein 18 (ZBTB18)58.320.8300.009
P10720Platelet factor 4 variant (PF4V1)11.550.7780.010
O95445Apolipoprotein M (APOM)21.240.8310.011
A0A0G2JPR0Complement C4-A (C4A)192.750.6200.012
Q99857Tenascin-C9.760.8180.013
Q07507Dermatopontin (DPT)23.990.7240.015
O95576Pepsinogen8.710.8200.015
C0KRQ8Glycoprotein hormones alpha chain (CGA)10.200.5420.016
P55056Apolipoprotein C-IV (APOC4)14.540.8170.020
P35908Keratin, type II cytoskeletal 2 epidermal (KRT2)65.390.8270.033
E1B4S7Apolipoprotein B (APOB)25.220.8310.034
P80108Phosphatidylinositol-glycan-specific phospholipase D (GPLD1)92.280.7850.034
H7C0N0Inter-alpha-trypsin inhibitor heavy chain H1 (ITIH1)51.830.8110.039
H7BZ76Latent-transforming growth factor beta-binding protein 1 (LTBP1)28.800.8060.040
Q7RTS7Keratin, type II cytoskeletal 74 (KRT74)57.830.8070.040
B2RCB8cDNA, FLJ95971, highly similar to Homo sapiens protocadherin 12 (PCDH12)128.910.8190.044

Key differentially expressed proteins identified by tandem mass tag-based proteomics analysis between PD patients and healthy controls.

aFC value was calculated as the ratio of the average mass response (area) between the two groups (FC value = PD/HC). Thus, FC values >1 indicate significantly higher levels in the PD group relative to the HC group while FC values <1 indicate significantly lower levels in the PD group. bOnly proteins with FC values greater than ± 1.2 and p-values less than 0.05 were deemed statistically significant. PD, Parkinson’s disease; MW, molecular weight; FC, fold change; HC, healthy control.

FIGURE 3

There were 162, 25, and 22 GO terms significantly associated with biological processes, cellular components, and molecular function, respectively. The top 20 GO terms from GO enrichment analysis are shown in Figure 4 and Table 5. Five of the six top ranked GO terms from cellular components (very-low-density lipoprotein particle, triglyceride-rich lipoprotein particle, lipoprotein particle, plasma lipoprotein particle, and protein-lipid complex) were directly associated with lipids. The remaining top ranked GO term from cellular components was keratin filament. Eleven of the fourteen top ranked GO terms from biological processes were directly associated with lipid metabolism. These were regulation of very-low-density lipoprotein particle clearance, negative regulation of very-low-density lipoprotein particle clearance, regulation of the fatty acid biosynthetic process, regulation of lipoprotein particle clearance, regulation of the fatty acid metabolic process, negative regulation of lipoprotein particle clearance, the lipid metabolic process, protein–lipid complex remodeling, plasma lipoprotein particle remodeling, protein–lipid complex subunit organization, and plasma lipoprotein particle organization. The remaining three GO terms were negative regulation of receptor-mediated endocytosis, organic hydroxyl compound transport, and macromolecular complex remodeling. They were also indirectly associated with lipid metabolism because all of the involved proteins were lipoproteins.

TABLE 5

Go termsCategoryInvolved proteinsq-valueRich factor
Very-low-density lipoprotein particleCCO95445, K7ER74, P55056, Q59HB3, P02656, K7ERI9, E1B4S70.0210.292
Triglyceride-rich lipoprotein particleCCO95445, K7ER74, P55056, Q59HB3, P02656, K7ERI9, E1B4S70.0210.292
Keratin filamentCCP04259, P35908, Q04695, Q7RTS70.0300.500
Lipoprotein particleCCO95445, K7ER74, P55056, P02656, K7ERI9, Q59HB3, E1B4S70.0300.206
Plasma lipoprotein particleCCO95445, K7ER74, P55056, P02656, K7ERI9, Q59HB3, E1B4S70.0300.206
Protein-lipid complexCCO95445, K7ER74, P55056, P02656, K7ERI9, Q59HB3, E1B4S70.0300.206
Negative regulation of receptor-mediated endocytosisBPK7ER74, P02656, K7ERI90.0211.000
Regulation of very-low-density lipoprotein particle clearanceBPK7ER74, P02656, K7ERI90.0211.000
Negative regulation of very-low-density lipoprotein particle clearanceBPK7ER74, P02656, K7ERI90.0211.000
Regulation of fatty acid biosynthetic processBPK7ER74, P02656, K7ERI90.0300.750
Regulation of lipoprotein particle clearanceBPK7ER74, P02656, K7ERI90.0300.750
Regulation of fatty acid metabolic processBPK7ER74, P02656, K7ERI90.0300.750
Negative regulation of lipoprotein particle clearanceBPK7ER74, P02656, K7ERI90.0300.750
Organic hydroxy compound transportBPO95445, K7ER74, G5E968, Q59HB3, P02656, K7ERI9, E1B4S70.0300.212
Lipid metabolic processBPO95445, Q1W658, O75339, K7ER74, P55056, P02656, K7ERI9, Q59HB3, D6RF20, E1B4S7, P801080.0300.133
Protein-lipid complex remodelingBPO95445, K7ER74, P02656, K7ERI9, Q59HB3, E1B4S70.0300.240
Plasma lipoprotein particle remodelingBPO95445, K7ER74, P02656, K7ERI9, Q59HB3, E1B4S70.0300.240
Protein-lipid complex subunit organizationBPO95445, K7ER74, P02656, K7ERI9, Q59HB3, E1B4S70.0300.240
Plasma lipoprotein particle organizationBPO95445, K7ER74, P02656, K7ERI9, Q59HB3, E1B4S70.0300.240
Macromolecular complex remodelingBPO95445, K7ER74, P02656, K7ERI9, Q59HB3, E1B4S70.0300.240

The top 20 most enriched GO terms and involved proteins based on proteomics analysis between PD patients and healthy controls.

GO, gene ontology; PD, Parkinson’s disease; CC, cellular component; BP, biological process.

FIGURE 4

These differentially expressed proteins were also mapped to KEGG pathways and enriched in five pathways with p-values less than 0.05 (Table 6), including ovarian steroidogenesis, the GnRH signaling pathway, neuroactive ligand–receptor interaction, glycosylphosphatidylinositol (GPI)-anchor biosynthesis, and regulation of lipolysis in adipocytes. However, no significant differences can be observed after being adjusted by Benjamini-Hochberg correction (q-value > 0.05), indicating that these results were unstable.

TABLE 6

Enriched KEGG pathwayInvolved proteinp-valueq-valueRich factor
Ovarian steroidogenesisQ1W658 C0KRQ80.0050.1340.667
GnRH signaling pathwayQ1W658 C0KRQ80.0250.2170.333
Neuroactive ligand-receptor interactionQ1W658 C0KRQ80.0340.2170.286
Glycosylphosphatidylinositol (GPI)-anchor biosynthesisP801080.0430.2171
Regulation of lipolysis in adipocytesC0KRQ80.0430.2171

Enriched KEGG pathways and involved proteins based on proteomics analysis between PD patients and healthy control.

KEGG, Kyoto Encyclopedia of Genes and Genomes; PD, Parkinson’s disease.

Discussion

Systemic biology is the most useful way to overview the metabolic alterations of an organism in a pathological condition. Here, we have adopted an integrated metabolomics and proteomics analysis of plasma in PD patients to explore its potential pathogenesis for the first time. No significant differences were found in any of the clinical data between these two groups, suggesting the following omics analyses were reasonable.

Researchers from the United States performed the first metabolomics profiling of plasma in PD patients (Bogdanov et al., 2008). They identified altered levels of 8-hydroxy-2-deoxyguanosine, uric acid, and glutathione as potential blood biomarkers for PD patients without further pathway analysis. Other researchers also performed metabolomics analysis of cerebrospinal fluid (Lewitt et al., 2013) and urine (Luan et al., 2015) from PD patients, suggesting various perturbed metabolic pathways related to lipids, energy, fatty acids, bile acids, polyamine, and amino acids (Shao and Le, 2019). In this untargeted metabolomics analysis, significant metabolic differences were found between PD patients and healthy controls from the PLS-DA plot. Eighty-three differentially expressed metabolites were identified between the two groups, most of which were lipid and lipid-like molecules. The lipid and lipid-like molecules were further classified, and sphingolipids accounted for 25%. Furthermore, pathway analysis indicated that sphingolipid metabolism was the only significantly enriched metabolic pathway, and the pathway tended to be activated. Sphingolipids are a class of lipids containing a backbone of sphingoid bases and a set of aliphatic amino alcohols (e.g., ceramides, sphingomyelins, gangliosides, and cerebrosides). The involvement of sphingolipid metabolism was first reported in a yeast model of Parkinson’s disease (Lee et al., 2011). The Saccharomyces cerevisiae lipid elongase null mutants exhibited severe growth defects, accumulation of reactive oxygen species, aberrant protein trafficking, and a dramatic decrease in the survival of aged cells. α-Synuclein (α-syn) is a small soluble synaptic protein that is the major proteinaceous component of Lewy bodies, the pathological hallmark of Parkinson’s disease (Jo et al., 2004). It was reported that the toxicity of α-syn could be enhanced as a result of the disruption of ceramide-sphingolipid homeostasis in the endoplasmic reticulum (Lee et al., 2011). Whereas Gaucher disease-related sphingolipids (glucosylceramide, glucosylsphingosine, sphingosine, sphingosine-1-phosphate) are reported to induce α-syn aggregation, glucosylsphingosine can trigger the formation of oligomeric α-syn in human neurons (Taguchi et al., 2017). Autophagy can be impaired by the increased levels of glucosylceramide and sphingomyelin while also enhanced by the increased levels of ceramide (Indellicato and Trinchera, 2019). In total, the disruption of the sphingolipid metabolism leads to α-syn accumulation, influences the traffic of vesicles toward lysosomes, and affects autophagy. Mutations of glucocerebrosidase 1, an enzyme hydrolyzing glucocerebroside (a subclass of sphingolipids), are the most common genetic risk factor for PD, although the underlying biological processes are poorly understood (Gan-Or et al., 2008). Metabolomics analysis mainly focuses on small molecules (<1 KD), and in this research a proteomics analysis was further performed.

A lot of proteomics research has been conducted in PD patients, mainly from brain tissues (substantia nigra, thalamus, locus coeruleus, olfactory bulb, and cerebral cortex), biofluids (cerebrospinal fluid, tears, and blood) and subcellular structures (mitochondria) (Dixit et al., 2019). Most of this research referred to pathways involved in “mitochondrial dysfunction, oxidative stress, protein aggregation or degradation, autophagy, and inflammation”(Ren et al., 2015; Monti et al., 2018). Here, forty proteins were differentially expressed in PD patients from the plasma proteomics analysis and the levels of all the seven apolipoproteins (apolipoprotein C-I, apolipoprotein C-III, protein APOC4-APOC2, apolipoprotein C-IV, apolipoprotein B variant, apolipoprotein B, and apolipoprotein M) were significantly decreased. Another former study reported decreases of apolipoprotein AII and apolipoprotein E in Parkinson’s disease (Zhang et al., 2008). The decrease of apolipoprotein B-100 (Lehnert et al., 2012) and apolipoprotein A-I (Zhang et al., 2012) was also reported in some studies. There are a large number of lipids in the central nervous system, and about one-fourth of the total cholesterol has a function or is stored in the brain. These apolipoproteins can influence the deposition process of many proteins in neurodegenerative diseases, such as α-syn in PD (Zhang et al., 2012). According to the enriched GO analysis, five of the six top ranking GO terms from cellular components were directly associated with lipids, and eleven of the fourteen top ranking GO terms from biological processes were directly associated with lipid metabolism. Keratin filament from the cellular component is a part of the intermediate filament and constitutes the cytoskeleton while Lewy bodies are the accumulation of neurofilaments (intermediate filament) in Parkinson’s disease. The underlying pathogenesis needs further research to confirm. The remaining three GO terms from biological processes were also related with lipid metabolism, as all the involved proteins were lipoproteins. KEGG pathway analysis was performed with those differentially expressed proteins and it seemed to be enriched in five pathways. The GnRH signaling pathway and ovarian steroidogenesis directly regulate steroid hormones and participate in the production of steroids. These processes belong to neuroactive ligand–receptor interaction. Glycosylphosphatidylinositol (GPI)–anchor biosynthesis and regulation of lipolysis in adipocytes are direct detailed processes of lipid metabolism. However, the results were unstable after Benjamini-Hochberg correction, as no significant differences could be observed (q-value > 0.05).

Both sphingolipids and apolipoproteins participate in plasma lipid metabolism, and are mainly related with atherosclerotic diseases. A number of lipidomics studies have reported specific lipid alterations in the brain or plasma from PD patients, such as the composition alteration of lipid rafts in the frontal cortex (Alecu and Bennett, 2019). The interaction of α-syn and the synaptic membrane is thought to be critical in Parkinson’s disease and has been explored for many years (Kubo et al., 2005). It is reported that lipids and the ratio of lipids to proteins can regulate the aggregation propensity of α-syn (Galvagnion, 2017). α-Syn can interact with the lipid membrane, and that interaction further affects the oligomerization and aggregation of α-syn (Suzuki et al., 2018). The alteration of lipid chemical properties can also induce α-syn aggregation and affect the balance between functional and aberrant behavior of α-syn (Galvagnion et al., 2016). In addition, α-syn can bind to specific lipid molecules, and these lipid–protein conjugates can further help the transport of α-syn in the blood or across the blood-brain barrier (Emamzadeh and Allsop, 2017), enhance its interaction with synaptic membranes, and interfere with the catalytic activity of cytoplasmic and lysosomal lipases, thereby disrupting lipid metabolism (Alecu and Bennett, 2019). In the meantime, lipid dysregulation can also promote PD pathogenesis through oxidative stress and inflammation reaction. Apolipoprotein M can increase the level of circulating sphingosine 1-phosphate, activate the following signaling pathway, and induce inflammation reaction (Kurano et al., 2013). Platelet activating factors, a lipid proinflammatory mediator, play an important role in modulating progressive neurodegeneration in PD patients by intracellular trafficking of α-syn. Furthermore, there are a large number of genetic risk factors of PD involved in lipid metabolism, including PLA2G6 and SCARB2. These two genes are directly or indirectly involved in glycerophospholipid and sphingolipid metabolism (Alecu and Bennett, 2019). Mutations in the glucocerebrosidase 1 gene, which encodes a degrading enzyme for the glycolipid glucosylceramide, are regarded as strong risk factors for PD and dementia with Lewy bodies, and glucosylceramide has been confirmed to promote toxic conversion of α-syn. Moreover, pathological research has demonstrated the existence of α-syn in the brains of patients with lysosomal storage disorders, in which glycosphingolipids is found to be accumulated (Suzuki et al., 2018). In addition, α-syn can antagonize neurotrophic signaling of TrkB by repressing TrkB lipid raft distribution, decreasing its internalization, and reducing its axonal trafficking (Kang et al., 2017). The possible pathogenesis is that α-syn sequesters the early peroxidation products of fatty acids, thereby reducing the level of highly reactive lipid species (De Franceschi et al., 2017). All in all, dysregulated lipid metabolism is implicated in Parkinson’s disease, which might open a new therapeutic method for modifying Parkinson’s disease.

Conclusion

Integrated proteomics and metabolomics analysis reveals that Parkinson’s disease is associated with plasma lipid metabolic disturbance. The activated sphingolipid metabolism and decreased apolipoproteins are the probable potential pathogenesis for PD.

Statements

Data availability statement

The datasets presented in this study can be found in the supplementary materials (Supplementary Table S1) and online repository PeptideAtlas, with accession number PASS01568 (http://www.peptideatlas.org/PASS/PASS01568).

Ethics statement

The studies involving human participants were reviewed and approved by the Ethics committee of the First Affiliated Hospital of Chongqing Medical University. The patients/participants provided their written informed consent to participate in this study.

Author contributions

LH and Y-DW designed this study. LH, C-QL, QQ, C-CZ, and M-XD assessed the clinical data. LH, X-MX, and YL performed the experiments. M-XD, Y-LH, and G-HC wrote and revised the manuscript.

Funding

This study was funded by the Chongqing Health and Family Planning commission (No. 2017MSXM023) and the Fundamental Research Funds for the Central Universities (2042020kf0056).

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2020.00080/full#supplementary-material

TABLE S1

The matrix of metabolomics data included in the paper.

References

  • 1

    AlecuI.BennettS. A. L. (2019). Dysregulated lipid metabolism and its role in α-Synucleinopathy in Parkinson’s disease.Front. Neurosci.13:328. 10.3389/fnins.2019.00328

  • 2

    BerardelliA.WenningG. K.AntoniniA.BergD.BloemB. R.BonifatiV.et al (2013). EFNS/MDS-ES/ENS [corrected] recommendations for the diagnosis of Parkinson’s disease.Eur. J. Neurol.201634. 10.1111/ene.12022

  • 3

    BlauwendraatC.NallsM. A.SingletonA. B. (2019). The genetic architecture of Parkinson’s disease.Lancet19170178.

  • 4

    BogdanovM.MatsonW. R.WangL.MatsonT.Saunders-PullmanR.BressmanS. S.et al (2008). Metabolomic profiling to develop blood biomarkers for Parkinson’s disease.Brain131389396. 10.1093/brain/awm304

  • 5

    ChongJ.SoufanO.LiC.CarausI.LiS.BourqueG.et al (2018). MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis.Nucleic Acids Res.46W486W494.

  • 6

    ConnollyB. S.LangA. E. (2014). Pharmacological treatment of Parkinson disease: a review.J. Am. Med. Assoc.31116701683.

  • 7

    De FranceschiG.FecchioC.SharonR.SchapiraA. H. V.ProukakisC.BellottiV.et al (2017). α-Synuclein structural features inhibit harmful polyunsaturated fatty acid oxidation, suggesting roles in neuroprotection.J. Biol. Chem.29269276937. 10.1074/jbc.m116.765149

  • 8

    DengH.WangP.JankovicJ. (2018). The genetics of Parkinson disease.Ageing Res. Rev.427285.

  • 9

    DixitA.MehtaR.SinghA. K. (2019). Proteomics in human Parkinson’s disease: present scenario and future directions.Cell Mol. Neurobiol.39901915. 10.1007/s10571-019-00700-9

  • 10

    DongM. X.HuL.HuangY. J.XuX. M.LiuY.WeiY. D. (2017a). Cerebrovascular risk factors for patients with cerebral watershed infarction: a case-control study based on computed tomography angiography in a population from Southwest China.Medicine96:e7505. 10.1097/md.0000000000007505

  • 11

    DongM. X.XuX. M.HuL.LiuY.HuangY. J.WeiY. D. (2017b). Serum Butyrylcholinesterase activity: a biomarker for Parkinson’s disease and related dementia.Biomed. Res. Int.2017:1524107.

  • 12

    DongM. X.HuQ. C.ShenP.PanJ. X.WeiY. D.LiuY. Y.et al (2016). Recombinant tissue plasminogen activator induces neurological side effects independent on thrombolysis in mechanical animal models of focal cerebral infarction: a systematic review and meta-analysis.PLoS One11:e0158848. 10.1371/journal.pone.0158848

  • 13

    DongM. X.FengX.XuX. M.HuL.LiuY.JiaS. Y.et al (2018a). Integrated analysis reveals altered lipid and glucose metabolism and identifies NOTCH2 as a biomarker for Parkinson’s disease related depression.Front. Mol. Neurosci.11:257. 10.3389/fnmol.2018.00257

  • 14

    DongM.-X.LiC.-M.ShenP.HuQ.-C.WeiY.-D.RenY.-F.et al (2018b). Recombinant tissue plasminogen activator induces long-term anxiety-like behaviors via the ERK1/2-GAD1-GABA cascade in the hippocampus of a rat model.Neuropharmacology128119131. 10.1016/j.neuropharm.2017.09.039

  • 15

    DorseyE. R.BloemB. R. (2018). The parkinson pandemic-a call to action.JAMA Neurol.75910.

  • 16

    EmamzadehF. N.AllsopD. (2017). α-Synuclein interacts with lipoproteins in plasma.J. Mol. Neurosci.63165172. 10.1007/s12031-017-0967-0

  • 17

    EmamzadehF. N.SurguchovA. (2018). Parkinson’s disease: biomarkers, treatment, and risk factors.Front. Neurosci.12:612. 10.3389/fnins.2018.00612

  • 18

    GalvagnionC. (2017). The role of lipids interacting with α-Synuclein in the pathogenesis of Parkinson’s disease.J. Parkinson’s Dis.7433450. 10.3233/jpd-171103

  • 19

    GalvagnionC.BrownJ. W.OuberaiM. M.FlagmeierP.VendruscoloM.BuellA. K.et al (2016). Chemical properties of lipids strongly affect the kinetics of the membrane-induced aggregation of α-synuclein.Proc. Natl. Acad. Sci. U.S.A.11370657070. 10.1073/pnas.1601899113

  • 20

    Gan-OrZ.GiladiN.RozovskiU.ShifrinC.RosnerS.GurevichT.et al (2008). Genotype-phenotype correlations between GBA mutations and Parkinson disease risk and onset.Neurology7022772283. 10.1212/01.wnl.0000304039.11891.29

  • 21

    HochbergY.BenjaminiY. (1990). More powerful procedures for multiple significance testing.Stat. Med.9811818. 10.1002/sim.4780090710

  • 22

    HuL.DongM. X.ZhaoH.XuG. H.QinX. Y. (2016). Fibulin-5: a novel biomarker for evaluating severity and predicting prognosis in patients with acute intracerebral haemorrhage.Eur. J. Neurol.2311951201. 10.1111/ene.13013

  • 23

    IndellicatoR.TrincheraM. (2019). The link between Gaucher disease and Parkinson’s disease sheds light on old and novel disorders of sphingolipid metabolism.Int. J. Mol. Sci.20:E3304.

  • 24

    JoE.DarabieA. A.HanK.TandonA.FraserP. E.McLaurinJ. (2004). alpha-Synuclein-synaptosomal membrane interactions: implications for fibrillogenesis.Eur. J. Biochem.27131803189. 10.1111/j.1432-1033.2004.04250.x

  • 25

    KangS. S.ZhangZ.LiuX.ManfredssonF. P.BenskeyM. J.CaoX.et al (2017). TrkB neurotrophic activities are blocked by α-synuclein, triggering dopaminergic cell death in Parkinson’s disease.Proc. Natl. Acad. Sci. U.S.A.1141077310778. 10.1073/pnas.1713969114

  • 26

    KuboS.NemaniV. M.ChalkleyR. J.AnthonyM. D.HattoriN.MizunoY.et al (2005). A combinatorial code for the interaction of alpha-synuclein with membranes.J. Biol. Chem.2803166431672. 10.1074/jbc.m504894200

  • 27

    KuranoM.TsukamotoK.OhkawaR.HaraM.LinoJ.KageyamaY.et al (2013). Liver involvement in sphingosine 1-phosphate dynamism revealved by adenoviral hepatic overexpression of apolipoprotein M.Atherosclerosis229102109. 10.1016/j.atherosclerosis.2013.04.024

  • 28

    LeeY. J.WangS.SloneS. R.YacoubianT. A.WittS. N. (2011). Defects in very long chain fatty acid synthesis enhance alpha-synuclein toxicity in a yeast model of Parkinson’s disease.PLoS One6:e15946. 10.1371/journal.pone.0015946

  • 29

    LehnertS.JesseS.RistW.SteinackerP.SoininenH.HerukkaS. K.et al (2012). ITRAQ and multiple reaction monitoring as proteomic tools for biomarker search in cerebrospinal fluid of patients with Parkinson’s disease dementia.Exp. Neurol.234499505. 10.1016/j.expneurol.2012.01.024

  • 30

    LewittP. A.LiJ.LuM.BeachT. G.AdlerC. H.GuoL. (2013). 3-hydroxykynurenine and other Parkinson’s disease biomarkers discovered by metabolomic analysis.Mov. Disord.2816531660. 10.1002/mds.25555

  • 31

    LuanH.LiuL. F.MengN.TangZ.ChuaK. K.ChenL. L.et al (2015). LC-MS-based urinary metabolite signatures in idiopathic Parkinson’s disease.J. Proteome Res.14467478. 10.1021/pr500807t

  • 32

    LvK.AeL. (2015). Parkinson’s disease.Lancet386896912.

  • 33

    MontiC.ColugnatI.LopianoL.ChioA.AlberioT. (2018). Network analysis identifies disease-specific pathways for Parkinson’s disease.Mol. Neurobiol.55370381. 10.1007/s12035-016-0326-0

  • 34

    MostafaI.ZhuN.YooM. J.BalmantK. M.MisraB. B.DufresneC.et al (2016). New nodes and edges in the glucosinolate molecular network revealed by proteomics and metabolomics of Arabidopsis myb28/29 and cyp79B2/B3 glucosinolate mutants.J. Proteom.138119. 10.1016/j.jprot.2016.02.012

  • 35

    RenR.SunY.ZhaoX.PuX. (2015). Recent advances in biomarkers for Parkinson’s disease focusing on biochemicals, omics and neuroimaging.Clin. Chem. Lab. Med.5314951506.

  • 36

    ShaoY.LeW. (2019). Recent advances and perspectives of metabolomics-based investigations in Parkinson’s disease.Mol. Neurodegen.14:3.

  • 37

    SidranskyE.NallsM. A.AaslyJ. O.Aharon-PeretzJ.AnnesiG.BarbosaE. R.et al (2009). Multicenter analysis of glucocerebrosidase mutations in Parkinson’s disease.N. Engl. J. Med.36116511661.

  • 38

    SuzukiM.SangoK.WadaK.NagaiY. (2018). Pathological role of lipid interaction with α-synuclein in Parkinson’s disease.Neurochem. Int.11997106. 10.1016/j.neuint.2017.12.014

  • 39

    TaguchiY. V.LiuJ.RuanJ.PachecoJ.ZhangX.AbbasiJ.et al (2017). Glucosylsphingosine promotes α-Synuclein pathology in mutant GBA-associated Parkinson’s disease.J. Neurosci.3796179631. 10.1523/jneurosci.1525-17.2017

  • 40

    TristB. G.HareD. J.DoubleK. L. (2019). Oxidative stress in the aging substantia nigra and the etiology of Parkinson’s disease.Aging Cell18:e13031.

  • 41

    WeiZ.LiX.LiX.LiuQ.ChengY. (2018). Oxidative stress in parkinson’s disease: a systematic review and meta-analysis.Front. Mol. Neurosci.11:236. 10.3389/fnmol.2018.00236

  • 42

    ZhangJ.SokalI.PeskindE. R.QuinnJ. F.JankovicJ.KenneyC.et al (2008). CSF multianalyte profile distinguishes Alzheimer and Parkinson diseases.Am. J. Clin. Pathol.129526529. 10.1309/w01y0b808emeh12l

  • 43

    ZhangX.YinX.YuH.LiuX.YangF.YaoJ.et al (2012). Quantitative proteomic analysis of serum proteins in patients with Parkinson’s disease using an isobaric tag for relative and absolute quantification labeling, two-dimensional liquid chromatography, and tandem mass spectrometry.Analyst137490495. 10.1039/c1an15551b

Summary

Keywords

Parkinson’s disease, metabolomics, proteomics, lipid, sphingolipid metabolism, apolipoprotein

Citation

Hu L, Dong M-X, Huang Y-L, Lu C-Q, Qian Q, Zhang C-C, Xu X-M, Liu Y, Chen G-H and Wei Y-D (2020) Integrated Metabolomics and Proteomics Analysis Reveals Plasma Lipid Metabolic Disturbance in Patients With Parkinson’s Disease. Front. Mol. Neurosci. 13:80. doi: 10.3389/fnmol.2020.00080

Received

16 December 2019

Accepted

23 April 2020

Published

30 June 2020

Volume

13 - 2020

Edited by

Andrei Surguchov, University of Kansas Medical Center, United States

Reviewed by

Irina G. Sourgoutcheva, University of Kansas Medical Center, United States; Nawajes Mandal, The University of Tennessee Health Science Center (UTHSC), United States

Updates

Copyright

*Correspondence: You-Dong Wei,

These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics