Dysregulation of the Kennedy Pathway and Tricarboxylic Acid Cycle in Myalgic Encephalomyelitis/Chronic Fatigue Syndrome
=======================================================================================================================

* Xiaoyu Che
* Christopher R. Brydges
* Yuanzhi Yu
* Adam Price
* Shreyas Joshi
* Ayan Roy
* Bohyun Lee
* Dinesh K. Barupal
* Aaron Cheng
* Dana March Palmer
* Susan Levine
* Daniel L. Peterson
* Suzanne D. Vernon
* Lucinda Bateman
* Mady Hornig
* Jose G. Montoya
* Anthony L. Komaroff
* Oliver Fiehn
* W. Ian Lipkin

## Abstract

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a chronic and debilitating disease that is characterized by unexplained physical fatigue unrelieved by rest. Symptoms also include cognitive and sensory dysfunction, sleeping disturbances, orthostatic intolerance and gastrointestinal problems. The pathogenesis is not fully understood. Using regression, Bayesian and enrichment analyses, we conducted targeted and untargeted metabolomic analysis of 888 metabolic analytes in plasma samples of 106 ME/CFS cases and 91 frequency-matched healthy controls. In ME/CFS cases, the regression, Bayesian and enrichment analyses all revealed abnormal levels of several membrane lipids indicating dysregulation of the Kennedy pathway: decreased plasma levels of plasmalogens, phosphatidylcholines, phosphatidylethanolamines, sphingomyelins, and phospholipid ethers. Enrichment analyses revealed decreased levels of cholines, ceramides and carnitines, and increased levels of long chain triglycerides, dicarboxylic acids, hydroxy-eicosapentaenoic acid, and the tricarboxylic acid cycle intermediates alpha-ketoglutarate and succinate. Using machine learning algorithms with selected metabolites as predictors, we were able to differentiate female ME/CFS cases from female controls (highest AUC=0.794) and ME/CFS cases without self-reported irritable bowel syndrome (sr-IBS) from controls without sr-IBS (highest AUC=0.873). Our findings are consistent with earlier ME/CFS work indicating compromised energy metabolism and redox imbalance, and highlight specific abnormalities that may provide insights into the pathogenesis of ME/CFS.

## Introduction

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a disease of unknown cause that is defined by impairment from fatigue lasting longer than six months, unrefreshing sleep, post-exertional malaise, and either cognitive dysfunction or orthostatic intolerance1. People with ME/CFS often report additional symptoms, such as gastrointestinal disturbances, influenza-like symptoms, and chronic pain2. It is estimated that ME/CFS affects between 0.4% to 2.5% of the global population, and 1.5 to 2.5 million people in the United States alone1, 3. There are no approved diagnostic tests for ME/CFS; medical providers must assess medical history, conduct a physical examination and exclude other disorders for diagnosis4, 5.

Prior metabolomic studies of patients with ME/CFS have provided insights into the potential pathogenesis and course of the disease, demonstrating disturbances in energy, lipid, amino acid, and redox metabolism6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16. Studies continue in the ongoing pursuit of identifying a metabolomic signature as biomarker for ME/CFS. However, there is a lack of consistency in the metabolites that are altered in ME/CFS patients across these studies. This may be the result of the heterogeneity of this disease.

Metabolic dimensions of ME/CFS may be related to sex; women are disproportionately affected by ME/CFS1, 17. Naviaux et al. (2016) found differences in metabolic pathway disturbances and altered metabolite levels when stratifying ME/CFS cases by sex15. Other have also reported sex-specific differences in plasma biomarkers14, 18, 19.

Comorbid gastrointestinal (GI) symptoms constitute a potential subtype in ME/CFS 11, 12, 14, 15, 19, 20, 21, 22, 23. Among those with ME/CFS, the presence or absence of self-reported physician diagnosed irritable bowel syndrome (sr-IBS), in particular, has highlighted differences in the plasma proteome relating to immune dysregulation and altered levels of metabolites within the metabolome14, 19. In a fecal metagenomics study, Nagy-Szakal et al. (2017) identified eleven bacterial species delineating differences between ME/CFS patients with and without sr-IBS and found relations between bacterial taxa and symptoms relating to fatigue and pain23.

In this study, we report targeted and untargeted analyses of 888 metabolic analytes comprising of primary metabolites, biogenic amines, complex lipids, and oxylipins in plasma of ME/CFS cases and controls. We identified altered metabolomic profiles between ME/CFS patients, controls, and subgroups within ME/CFS patients based on sex and sr-IBS.

## Materials and Methods

### Study population

Our starting population comprised 177 ME/CFS cases and 177 controls in ME/CFS clinics in Incline Village, NV; Miami FL; New York, NY; Salt Lake City, UT; and Palo Alto, CA. All ME/CFS cases met the 1994 CDC Fukuda24 and Canadian consensus criteria for ME/CFS25. All ME/CFS cases completed standardized screening and assessment instruments including medical history and symptom rating scales as well as a physical examination. Controls were matched to cases on age, sex, race/ethnicity, geographic/clinical site, and date of sampling (±30 days). Based on screening criteria, we excluded 5 ME/CFS cases that met any exclusion criteria from the 1994 CDC Fukuda and/or Canadian consensus criteria for ME/CFS such as having chronic infections, rheumatic and chronic inflammatory diseases, neurological disorders, psychiatric conditions, or were taking any immunomodulatory medication. Controls underwent the same screening process as ME/CFS subjects and were excluded if they reported ME/CFS or other conditions deemed by the recruiting physician to be inconsistent with a healthy control population. Controls were also excluded if they had a history of substance abuse, psychiatric illness, antibiotics in the prior three months, immunomodulatory medications in the prior year, and clinically significant findings on physical exam or screening laboratory tests. One control was excluded after prescreening based on these criteria. Additionally, 21 participants were excluded prior to baseline due to withdrawal from the study (n=18), loss to follow-up (n=2), and enrollment capacity (n=1). The baseline questionnaire was completed by with 327 participants. During the study, an additional 63 participants were excluded for study protocol deviations (n=25), loss to follow-up (n=25), and withdrawal from the study (n=13), resulting in a total of 264 participants.

For the analysis reported here, a sub-cohort was established based on complete survey and biospecimen data (blood, saliva and stool) at the first and last time points of the study and key demographic characteristics were frequency-matched to ensure that the nested cohort was similar to the full cohort. This sub-cohort consisted of 106 ME/CFS cases and 91 controls; the derivation of the sub-cohort is summarized in Figure 1. All participants provided informed written consent in accordance with protocols approved by the Institutional Review Board at Columbia University Irving Medical Center.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/22/2021.06.14.21258895/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/06/22/2021.06.14.21258895/F1)

Figure 1. 
Pipeline for sample selection.

### Plasma collection

Blood samples were collected into BD VacutainerTM Cell Preparation Tubes (CPT) with ethylenediaminetetraacetic acid (EDTA) anticoagulant between January, 2016 and June, 2016, and centrifuged to pellet red blood cells. The plasma was shipped to Columbia University at 4°C. After aliquoting, samples were stored at -80°C until thawed for metabolomics analyses. All the samples were analyzed within two years of collection.

### Clinical assessment

Clinical symptoms and baseline health status were assessed on the day of physical examination and biological sample collection from both case and control subjects using the following instruments: the Short Form 36 Health Survey (SF-36), the Multidimensional Fatigue Inventory (MFI), DePaul Symptom Questionnaire (DSQ)26, and Pittsburgh Sleep Quality Index (PSQI)27. The SF-36 includes the following subject-reported evaluations about current health status: physical and social functioning, physical and emotional limitations, vitality, pain, and general and mental health28. The MFI comprises of a 20-item self-reported questionnaire focused on general, physical and mental fatigue, reduced activity, and reduced motivation29. Cognitive function was tested based on the DSQ questionnaire data and was scored using a standard cognitive disturbance definition as well as a modified definition based on a subset of questionnaire variables. Sleeping disturbances linked to ME/CFS were tested and scored based on DSQ and PSQI questionnaire items. Each instrument was transformed into a 0–100 scale to facilitate combination and comparison wherein a score of 100 is equivalent to maximum disability or severity and a score of zero is equivalent to no disability or disturbance.

A diagnosis of sr-IBS was based on answers in the medical history form. Subjects were asked if they had received a previous IBS diagnosis by a physician and the date of that diagnosis. Of the 106 subjects with ME/CFS, 35 (33.0%) had sr-IBS. Of the 91 control subjects, 3 (3.3%) had sr-IBS.

### Metabolomics analysis

Samples were stored at -80°C before analysis. Untargeted metabolomics data were acquired using three chromatography/mass spectrometry-based assays (MS): (1) Primary metabolites such as mono- and disaccharides, hydroxyl- and amino acids were measured by gas chromatography/time-of-flight mass spectrometry (GC-TOF MS30) including data alignment and compound annotation using the BinBase database algorithm31. (2) Biogenic amines including microbial compounds such as trimethylamine N-oxide (TMAO), methylated and acetylated amino acids and short di- and tripeptides were measured by hydrophilic interaction liquid chromatography/quadrupole time-of-flight mass spectrometry (HILIC-QTOF MS). (3) Complex lipids including phosphoglycerolipids, triacylglycerides, sphingolipids, and free fatty acids were analyzed by liquid chromatography (LC)/quadrupole time-of-flight mass spectrometry (CSH-QTOF MS32). Targeted bioactive oxylipin assay included thromboxanes, prostaglandins, and hdyroxy-, keto- and epoxy-lipins. All LC-MS/MS data included diverse sets of internal standards. LC-MS data were processed by MS-DIAL vs. 4.0 software33, and the compounds were annotated based on accurate mass, retention time and MS/MS fragment matching using LipidBlast34 and Massbank of North America libraries35. MS-FLO was used to remove erroneous peaks and reduce the false discovery rate in LC datasets36. A total of 821 known metabolites were annotated. Some complex lipids were annotated in both positive (ESI+) and negative (ESI-) ion modes, resulting in a total of 888 metabolic analytes that were included in our analysis. Data were normalized by SERRF37. Residual technical errors were assessed by coefficients of variation (CV) for known metabolites.

### Statistical analyses

For each metabolic analyte, zero values reflecting a measurement below the detection limit, were replaced with 50% of its smallest available value. In each of the four metabolomics panels, outliers were identified through principal component analysis (PCA). In primary metabolites (PM), 6 outliers (4 cases and 2 controls) were identified and removed; in complex lipids (CL), there were 5 outliers (3 cases and 2 controls); in oxylipins (OL), there was 1 outlier (1 case); in biogenic amines (BA), 4 outliers (3 cases and 1 control) were eliminated.

To compare the levels of each metabolite between ME/CFS cases and controls, we employed a variety of regression models with the metabolite level as the dependent variable and the binary case/control status as the independent variable, adjusting for all the matching variables (age, sex, race/ethnicity, geographic/clinical site, and season of sampling), body mass index (BMI) and sr-IBS. We considered two options for the dependent variable: 1) original metabolite levels, and 2) natural log-transformed metabolite levels. Before log-transformation, if necessary, all data points in metabolic analytes were multiplied by a minimal factor to keep the feature on a positive domain. Four regression models were considered: Gaussian regression with identity link, Gaussian regression with log link, lognormal regression and Gamma regression with log link. The Bayesian information criterion (BIC) was used to select the best fitting transformation/regression combination. We then calculated the estimated coefficient for the case/control status, together with its 95% confidence interval (95% CI) and p-value. Multiple comparisons over all metabolites were corrected using the Benjamini-Hochberg procedure38 controlling the false discovery rate (FDR) at the 0.15 level. Additionally, chemical enrichment analyses were performed using ChemRICH39 to determine chemical classes that were significantly altered between groups. ChemRICH does not rely upon background databases for statistical calculations and provides enrichment analysis based upon chemical structure, as opposed to defined pathways that can be inherently flawed39.

For each metabolite, we also conducted Bayesian analysis with the best fitting transformation/regression combination using R packages “rstanarm”40 and “bayestestR”41. Default (weakly informative) prior distributions from rstanarm were applied adjusting the scales of the priors internally. We then calculated the Bayes factors (BFs) and 95% highest density credible intervals (HDIs). The BF of a single parameter indicates the degree by which the mass of the posterior distribution has shifted further away from or closer to the null value (zero), relative to the prior distribution42. Hence, the BF measures the strength of evidence in favor of the alternative hypothesis (β≠0) over the null hypothesis (β=0). The 95% credible interval in the Bayesian framework is the range, within which the effect has 95% probability of falling, given the observed data. It has a different interpretation from the 95% confidence interval in the frequentist framework which instead signifies that with a large number of repeated samples, 95% of such calculated confidence intervals would include the true value of the parameter. We considered a metabolite significantly associated with ME/CFS if it satisfied the following criteria: 1) FDR adjusted p-value < 0.15, 2) BF > 3, and 3) 95% HDIs not covering 0. Jeffreys (1961)43 suggested that the strength of evidence for the alternative hypothesis compared to the null hypothesis is regarded as noteworthy if BFs are above 3.

Naviaux et al. (2016)15 showed that potential diagnostic metabolites for ME/CFS in targeted metabolomics are different between male and female subjects. Accordingly, we conducted sex-stratified analyses in addition to analyses with the whole cohort. In our previous work with a different cohort, sr-IBS comorbidity was identified as the strongest driving factor in the separation of topological networks based on fecal microbiome and plasma metabolic pathways14, 23. We subsequently found different patterns in the relationships between plasma proteomic profiling and ME/CFS when comparing ME/CFS with or without sr-IBS to healthy controls19. Given this precedent, we tested the hypothesis that sr-IBS subgroups in ME/CFS patients have altered metabolic profiles in a stratified analysis. As there were only 3 control subjects with sr-IBS, we focused on the comparison of ME/CFS subjects without sr-IBS versus controls without sr-IBS.

To explore the utility of the metabolomics assay as a biomarker tool for ME/CFS, we employed four machine learning algorithms: least absolute shrinkage and selection operator (Lasso)44, adaptive Lasso (AdaLasso)45, Random Forests (RF)46 and XGBoost47. AdaLasso is different from Lasso in that AdaLasso has the oracle property that leads to consistent variable selection whereas Lasso is only consistent for variable selection under certain conditions on the shrinkage parameters and correlations48. However, neither outperforms the other consistently in predictions. For each of the algorithm, three sets of predictors were considered: 1) all metabolites, 2) metabolites with BF>1, and 3) metabolites with BF>3. The predictive models were first trained in the 80% randomly-selected training set using 10-fold cross-validation; the remaining 20% of the study population was used as the independent test set to validate model performance. We also applied the Bayesian Model Averaging (BMA) method49 that combines the predictions of multiple models using weighted averages in which the weights are Bayesian posterior probabilities that the given model is the true model, conditional on the training data. The predictive performance of the 5 models (Lasso, AdaLasso, RF, XGBoost and Model Average) using the three sets of predictors in the test set was evaluated using Area under the Receiver Operating Characteristic curve (AUROC) values and Receiver Operating Characteristic (ROC) curves.

Data analyses were performed using MATLAB Statistics Toolbox R2013a (MathWorks, Inc., Natick, MA) and R version 3.6.3 (RStudio, Inc., Boston, MA). All p-values were 2-tailed.

## Results

### Study population characteristics

The study included plasma samples from 106 ME/CFS cases and 91 healthy controls recruited from five sites across the United States. Demographic and clinical characteristics of the study population are shown in Table 1. ME/CFS cases and controls were similar for all the frequency matching variables except season of collection (Chi-squared p = 0.004). We adjusted for all the matching variables, BMI and sr-IBS in our statistical analyses to account for confounding. All scales in SF-36 and MFI were significantly different between the two cohorts (Wilcoxon rank-sum p < 0.001). The study population is similar to the prescreened cohort that consisted of 177 ME/CFS cases and 177 controls in sex (Chi-squared p=0.60), race (Chi-squared p=0.66) and age (Wilcoxon rank-sum p=0.65).

View this table:
[Table 1.](http://medrxiv.org/content/early/2021/06/22/2021.06.14.21258895/T1)

Table 1. 
Subject characteristics.

### Metabolomic dataset

Targeted and untargeted mass spectrometry platforms yielded data for 888 metabolic analytes comprising 100 primary metabolites (PM), 237 biogenic amines (BA), 480 complex lipids (CL), and 71 bioactive oxylipins (OL). Supplementary Table S1 shows the sample mean and the standard deviation (SD) of levels of each metabolite within all ME/CFS cases, all controls, female ME/CFS cases, female controls, male ME/CFS cases, male controls, ME/CFS cases without sr-IBS and controls without sr-IBS.

### ME/CFS is associated with altered metabolomic profile

In PM, BA and CL panels, lognormal regression models with log-transformed metabolite levels as dependent variables had the lowest BIC values and best fit the data; the estimated coefficients can be interpreted as the differences in the mean values of log-log transformation of metabolite levels between cases and controls. In OL panel, a mixture of lognormal and log-link Gamma regression models with original metabolite levels as dependent variables best fit the data. For lognormal regression models, the estimated coefficients are interpreted as the mean differences of log transformation of metabolite levels between two groups. For log-link Gamma regression models, the estimated coefficients are interpreted as the log of fold change between two groups.

We did not identify any metabolite as significantly associated with ME/CFS in the PM panel. In the BA panel, levels of acetaminophen were increased in ME/CFS cases compared to controls. In the CL panel, we found decreased levels of plasmalogens, unsaturated phospholipid ethers, unsaturated phosphatidylcholines (PC), an unsaturated sphingomyelin (SM), and an unsaturated lysophosphatidylcholines (LPC) in ME/CFS cases compared to controls. In the OL panel, decreased levels of Resolvin D1 were observed in ME/CFS cases compared to controls. Table 2 shows the estimated coefficients in the regression models of these metabolites, their associated 95% CIs, p-values, FDR adjusted p-values and BFs. Because we used weakly informative priors in Bayesian analysis, the 95% HDIs were extremely similar to the 95% CIs. We report estimations of HDIs in Supplementary Table S2 where estimations for all metabolites are shown.

View this table:
[Table 2.](http://medrxiv.org/content/early/2021/06/22/2021.06.14.21258895/T2)

Table 2. 
Metabolites significantly associated with ME/CFS or ME/CFS subgroups.

Set enrichment analysis of the results from the regression models (Figure 2A) revealed that ME/CFS subjects had reduced levels of plasmalogens, sphingomyelins, unsaturated phospholipid ethers, unsaturated ceramides, carnitines, saturated lysophospholipids, unsaturated lysophosphoethanolamines, unsaturated lysophosphatidylcholines, saturated triglycerides and prostaglandins. The majority of unsaturated phosphatidylcholines were also down-regulated in ME/CFS cases. Increased levels of hydroxy-eicosapentaenoic acid (HEPE), dicarboxylic acids, and the majority of unsaturated long chain triglycerides were found in ME/CFS cases compared to controls. There were mixed directional alterations in the food exposome and epoxy fatty acids (EpODE). Complete data from ChemRICH enrichment analysis are provided in Supplementary Table S3. Data from compound-level enrichment analysis for the significantly altered metabolic clusters are illustrated in Supplementary Table S4. Levels of choline in food exposome were reduced in ME/CFS (estimated coefficient β=-0.009, p=0.004); levels of succinic acid (β=0.022, p=0.007) and alpha-ketoglutarate (β=0.016, p=0.048) in dicarboxylic acids were elevated in ME/CFS.

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/22/2021.06.14.21258895/F2/graphic-7.medium.gif)

[](http://medrxiv.org/content/early/2021/06/22/2021.06.14.21258895/F2/graphic-7)

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/22/2021.06.14.21258895/F2/graphic-8.medium.gif)

[](http://medrxiv.org/content/early/2021/06/22/2021.06.14.21258895/F2/graphic-8)

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/22/2021.06.14.21258895/F2/graphic-9.medium.gif)

[](http://medrxiv.org/content/early/2021/06/22/2021.06.14.21258895/F2/graphic-9)

Figure 2. Chemical enrichment analyses using ChemRICH
**a** All ME/CFS v. controls. **b** Female ME/CFS v. female controls. c ME/CFS without sr-IBS v. controls without sr-IBS. The length of the bar represents altered ratio for each metabolic cluster. A bar restricted to the left of the centered vertical line indicates a metabolic cluster that is lower in ME/CFS patients. A bar restricted to the right of the centered vertical line indicates a metabolic cluster that is higher in ME/CFS patients. A bar that crosses the vertical line indicates a metabolic cluster that is dysregulated in mixed directions. The color represents significance. EpODE: epoxy octadecadienoic acid. HEPE: hydroxy eicosapentaenoic acid. ME/CFS: myalgic encephalomyelitis/chronic fatigue syndrome. sr-IBS: self-reported physician diagnosed irritable bowl syndrome.

### Altered metabolomic profiles in female and male ME/CFS patients

Naviaux et al. (2016)15 discovered that women with ME/CFS, but not men, had disturbed fatty acid and endocannabinoid metabolism. Accordingly, we repeated separately the analyses in female and male cohorts in our study population.

In female subjects, regression and Bayesian analyses (Table 2) revealed that levels of phosphatidylcholines (PC), phosphatidylethanolamines (PE) and sphingomyelins (SM) in the CL panel were decreased in ME/CFS patients compared to controls. In the BA panel, levels of two drug metabolites, alprazolam and acyclovir, were up-regulated in ME/CFS patients. We did not find the elevated levels of acetaminophen in female subjects that were observed in the entire ME/CFS (male and female population), presumably due to loss of power. Enrichment analysis in female subjects (Figure 2B) identified dysregulations in the same metabolic clusters as in the overall population. Complete data from enrichment analysis in female subjects are shown in Supplementary Table S6. In contrast, we did not find any metabolites significantly associated with risk of ME/CFS in male subjects. This may be due to limited sample size. Supplementary Table S5 shows the regression and Bayesian estimations for all metabolites in male and female cohorts.

### Altered metabolomics profile in ME/CFS patients without sr-IBS

Due to the limited sample size of subjects with sr-IBS (35 ME/CFS cases and 3 controls), we only compared levels of metabolites between ME/CFS cases without sr-IBS and controls without sr-IBS. Levels of phosphatidylcholines (PC) and phosphatidylethanolmines (PE) were decreased in ME/CFS patients in this subgroup (Table 2). In the ChemRICH enrichment analysis, the dysregulations in metabolite clusters found to be dysregulated in the subgroup without sr-IBS (Figure 2C) were all identified in the overall population (Figure 2A). Complete data pertaining to the regression, Bayesian and enrichment analyses are shown in Supplementary Tables S7 and S8.

### Assessment of the metabolomics assay as a potential diagnostic tool for ME/CFS

We considered three sets of metabolites as predictors to distinguish ME/CFS cases from controls, including all metabolites, metabolites with BF>1 and metabolites with BF>3. Each set of predictors was fitted in five different machine learning classifiers: Lasso, adaptive Lasso (AdaLasso), Random Forests (RF), XGBoost and Model Average. The classifiers were first trained in the 80% randomly-selected training set and then validated in the remaining 20% test set. Figure 3A, Figure 3B, and Figure 3C show the ROC curves and the AUROC values differentiating all ME/CFS cases from all controls, female ME/CFS from female controls, and ME/CFS without sr-IBS from controls without sr-IBS, respectively, in the test set. Although classifiers did not differentiate all ME/CFS from all controls, Lasso with BF>1 metabolites as predictors distinguished female ME/CFS patients from female controls with an AUROC value of 0.794 (95% CI: 0.612-0.976) and Lasso with BF>3 metabolites distinguished ME/CFS without sr-IBS from controls without sr-IBS with an AUROC value of 0.873 (95% CI: 0.747-0.999). The AUROC values and their associated 95% CIs of all the classifiers are shown in Supplementary Table S9.

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/22/2021.06.14.21258895/F3/graphic-10.medium.gif)

[](http://medrxiv.org/content/early/2021/06/22/2021.06.14.21258895/F3/graphic-10)

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/22/2021.06.14.21258895/F3/graphic-11.medium.gif)

[](http://medrxiv.org/content/early/2021/06/22/2021.06.14.21258895/F3/graphic-11)

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/22/2021.06.14.21258895/F3/graphic-12.medium.gif)

[](http://medrxiv.org/content/early/2021/06/22/2021.06.14.21258895/F3/graphic-12)

Figure 3. ME/CFS predictive modeling.
**a** Overall population. **b** Women only. **c** No GI complaints. To differentiate ME/CFS cases from healthy controls, we employed five machine learning algorithms: Lasso, adaptive Lasso (AdaLasso), Random Forests (RF), XGBoost, and Bayesian Model Averaging (Model average). For each algorithm, three sets of predictors were considered: 1) all metabolites, 2) metabolites with BayesFactor > 1, and 3) metabolites with BayesFactor > 3. The predictive models were first trained in the 80% randomly selected training set using 10-fold cross-validation, and the remaining 20% of the study population was used as the independent test set to validate model performance. ME/CFS: myalgic encephalomyelitis/chronic fatigue syndrome. sr-IBS: self-reported physician diagnosed irritable bowel syndrome. BF: BayesFactor. AUC: area under the receiver operating characteristic curve.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/22/2021.06.14.21258895/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2021/06/22/2021.06.14.21258895/F4)

Figure 4. Functional interaction network of altered metabolic clusters in ME/CFS
Metabolite levels that are decreased (orange), increased (blue), or mixed in direction (green) in the enrichment assay have been associated with oxidative stress, mitochondrial dysfunction, and neurodegeneration.

### Metabolites significantly associated with ME/CFS or ME/CFS subgroups do not strongly correlate with duration of illness or symptom severity scores

We investigated whether the plasma levels of metabolites, defined by their associations with ME/CFS and ME/CFS subgroups (Table 2), correlated with duration of illness in years, Short Form 36 Health Survey (SF-36) scales and Multidimensional Fatigue Inventory (MFI) scales using Spearman’s correlation tests. None of the correlation coefficients exceeded the absolute magnitude of 0.5 (data not shown).

## Discussion

Since the first reports of large-scale metabolomic studies in people with ME/CFS were published in 2016 by Naviaux15, Yamano16, and Fluge9, several research teams, including our own14 have reported metabolomic analyses of plasma. The majority describe abnormalities in energy metabolism, with most reporting decreased levels of phospholipids and suggesting abnormalities in mitochondrial activity that could contribute to fatigue and cognitive dysfunction. Here we report confirmation of decreased levels of phospholipids, including phosphatidylcholine and sphingolipids, and provide evidence for dysregulation of the Kennedy pathway and the tricarboxylic acid cycle. Our findings provide potential insights into the pathobiology of clinical features of ME/CFS by providing a mechanistic framework for understanding compromised energy production, loss of integrity of cellular and mitochondrial membranes, inflammation, impaired cognition, dysregulated autonomic function, impediments to repair of tissue injury, and redox imbalance.

The Kennedy pathway is responsible for the biosynthesis of phosphatidylcholines (PC) and phosphatidylethanolamines (PE), the two most abundant phospholipids in mammalian cells50. PC, the most abundant phospholipid in the mitochondrial membranes51, 52, is sourced from endoplasmic reticulum. PE is synthesized in mitochondria by the decarboxylation of phosphatidylserine by phosphatidylserine decarboxylase 1 (Psd1) at the inner mitochondrial membrane53. PC and PE are essential to the formation of intermediate structures in membrane fusion and fission events, for stabilizing membrane proteins into their optimal conformations, and for actin-filament disassembly in the end stage of cytokinesis54, 55, 56. In people with ME/CFS, we found decreased levels of PC and PE and their downstream products: ceramides, sphingomyelins, lysophosphatidycholines, phospholipid ethers, prostaglandin D2 (PGD2) and prostaglandin F2α (PGF2α).

One critical functional implication of reduced levels of PC and PE is impaired oxidative phosphorylation. PC depletion specifically affects the function of inner membrane protein translocases of mitochondria, including the TIM23 complex57. PE synthesis is critical for cytochrome bc1 complex III function in the mitochondrial inner membrane53. Preprotein binding to the TIM/TOM complex, which translocates proteins produced from nuclear DNA through the mitochondrial membrane for use in oxidative phosphorylation, is disturbed in PE-deficient mitochondria58, 59. Cytochrome c oxidase activity in the respiratory chain complex is also decreased with PE-deficiency60, 61. Reduced import of PE into the mitochondria results in the formation of respiration deficient cells55, and in mitochondrial dysfunction. Finally, reduced levels of lysophosphatidycholines and phospholipid ethers, as well as of PC and PE, can impede mitochondrial respiration55. Reduced synthesis of PGF2α and PGD2 in phospholipase A2γ-deficient mice induces mitochondrial dysfunction as well as oxidative stress that can contribute to further mitochondrial damage62.

Because PE and PC, and downstream metabolites in the Kennedy pathway, are important components of the lipid bilayer, the reduction in their levels has implications for signaling. Alteration in the levels or conformation of membrane components can adversely affect the functionality of proteins embedded in the membranes such as G protein coupled receptors (GPCRs). Phospholipids can act as direct allosteric modulators of GPCR activity through the lipid head group that affect ligand binding (agonist and antagonist) and receptor activation. For example, PE favors antagonist binding and stabilizes the inactive state of the receptor, whereas phosphatidylglycerols favor agonist binding and activation63.

Both PE and PC are precursors to many biologically active molecules that can act as second messengers. Prominent among them are Di-acyl glycerol (DAG), fatty acids, phosphatidic acid, lysophosphatidic acid, N-arachidonylethanolamine, N-palmitoylethanolamine, N-steroylethanolamine and arachidonic acid64, 65, 66. Ceramides, are not only structural components of membranes, but can also act as second messengers in modulating a range of cellular signaling pathways67.

Metabolomic analyses also revealed reductions in levels of plasmalogens and of resolvin D1. Plasmalogens are compounds with antioxidant functions that are synthesized by peroxisomes. Resolvin D1, a derivative of docosahexanoic acid (DHA), may contribute to resolution of inflammation by targeting dead cells for clearance by macrophages.68 As previously reported, 14, 69 we found a significant reduction in levels of carnitine. Carnitine regulates the cellular to mitochondrial ratio of free CoA to Acyl-CoA, removes the unwanted acyl groups and plays a key role in the transport of long-chain fatty acids from cytoplasm to the mitochondrial matrix for oxidation70. Depletion of carnitine reduces the generation of ATP from fatty acids, and may adversely impact the integrity of cell and mitochondrial membranes, and responses to inflammation and oxidative stress70. Finally, reduced levels of carnitine threaten the integrity of cell and mitochondrial membranes, increase oxidative stress, and reduce the ability to counter inflammation71. We also observed increased levels of long-chain triglycerides in ME/CFS. Depletion of carnitine leads to the accumulation of long-chain triglycerides that become targets for lipid peroxidation by mitochondria72. The accumulation of toxic lipid peroxidation products can lead to mitochondrial membrane damage18.

The tricarboxylic acid (TCA) cycle is a conserved pathway in aerobic organisms through which the acetyl-CoA from carbohydrates, fats and proteins is converted into ATP. We observed elevated levels of two TCA cycle intermediates, the dicarboxylic acids α-ketoglutarate (α-KG) and succinate in ME/CFS. Increased levels of α-KG have been reported previously in ME/CFS patients11, although we are not aware of previous reports of elevated levels of succinate. Abnormal levels of TCA cycle intermediates suggest inefficiencies in ATP production that may contribute to the fatigue and post-exertional malaise reported in ME/CFS. Increases in α-KG levels have been reported to induce severe metabolic impairment of pyruvate oxidation in the tricarboxylic acid cycle, leading to cell death73. Succinate accumulation has been reported to induce HIF-1α stabilization as well as the transcriptional activation of the pro-inflammatory cytokine IL-1β73. Elevated succinate levels contribute to increased oxidative stress and neuronal degeneration in rat models74. Oxidative stress, in turn, augments nitrosative stress75. Nitrosative stress, which has been documented in people with ME/CFS76, can lead to the increased production of peroxynitrite and downregulate the function of both alpha-ketoglutarate dehydrogenase and succinate dehydrogenase75, 76, 77. Infection is a common cause of nitrosative stress, and many ME/CFS patients report symptoms consistent with system infection prior to diagnosis.

Choline depletion was another potentially important finding in our study. Choline is an essential nutrient; 95% of it is utilized in the synthesis of PC via the Kennedy pathway50. The remaining 5% exists as either free choline or is used in the synthesis of phosphocholine, glycerophosphocholine, CDP-choline, acetylcholine, and other choline-containing phospholipids like sphingomyelin, plasmalogens and lysophosphatidylcholine. Each of these compounds contributes to maintenance of the structure and signaling functionality of the plasma membrane50, 66. IgG autoantibodies that specifically target GPCRs have been reported, even in healthy individuals, but are more commonly found in ME/CFS78, 79, particularly to autonomic nervous system targets including the M3 Acetylcholine receptor (M3AChR) and β2 Adrenergic receptor (β2AdR). Agonists for each of these receptors have choline precursors, acetylcholine (AC) and epinephrine (adrenaline), respectively. Choline also plays a role in the production of epinephrine, by donating the methyl group. Thus, choline deficiency could potentially lead to the autonomic dysfunction that is found in many people with ME/CFS, with reduced tissue blood flow and oxygen supply, leading to hypoxia, ischemia and fatigue80.

Impairments in cognition have been reported in ME/CFS81, 82. Our metabolomic data revealed reductions in levels of sphingomyelin, ceramides, and plasmalogens that may contribute to central nervous system dysfunction. Reduced levels of sphingomyelins have been reported in neurological disorders such as Alzheimer’s disease (AD)83, 84, Parkinson’s disease (PD)85 and multiple sclerosis86. In these examples of neurodegeneration, levels of sphingomyelin were reduced, but ceramide was increased87. In contrast, our data show reduced levels of both sphingomyelin and ceramide. Reduced levels of ceramide are associated with decreased cell survival, aberrant Purkinje cell dendritic differentiation88, and neurons having shorter axon plexus and fewer axonal branches89. Reductions in mitochondrial ceramide levels have also been shown to result in neuronal degeneration and reduced mitochondrial respiratory function as manifest with decreased mitochondrial basal and maximal oxygen consumption rates and decreases in spare respiratory capacity89. Reduced plasmalogen levels have been observed AD and PD90, 91, 92. Plasmalogen deficiencies may increase vulnerability of neural membranes to oxidative stress, destabilize membranes, and impair muscarinic cholinergic signaling and abnormal amyloid precursor processing91, 93, 94, 95.

## Conclusion

Our findings indicate a series of interconnected metabolic alterations in people with ME/CFS, that are consistent with two central abnormalities that may contribute to the pathogenesis of ME/CFS: disturbances in the Kennedy Pathway leading to reductions in levels of PC and PE and structural and functional disturbances of cellular and mitochondrial membranes, and reductions in levels of α-ketoglutarate and succinate that are consistent with an impairment in the TCA cycle.

## Supporting information

Supplemental Table S1 [[supplements/258895_file08.xlsx]](pending:yes)

Supplemental Table S2 [[supplements/258895_file09.xlsx]](pending:yes)

Supplemental Table S3 [[supplements/258895_file10.xlsx]](pending:yes)

Supplemental Table S4 [[supplements/258895_file11.xlsx]](pending:yes)

Supplemental Table S5 [[supplements/258895_file12.xlsx]](pending:yes)

Supplemental Table S6 [[supplements/258895_file13.xlsx]](pending:yes)

Supplemental Table S7 [[supplements/258895_file14.xlsx]](pending:yes)

Supplemental Table S8 [[supplements/258895_file15.xlsx]](pending:yes)

Supplemental Table S9 [[supplements/258895_file16.xlsx]](pending:yes)

## Data Availability

Data will be uploaded upon acceptance into a peer-reviewed journal.

## Ethics Statement

### Conflict of Interest

The authors declare no conflict of interest.

## Author Contributions

Author contributions: X.C., B.L., D.K.B., O.F. and W.I.L. developed the experimental design; S.L., D.L.P., S.D.V., L.B. and J.G.M. contributed to case and control recruitment and characterization; X.C., C.R.B., Y.Y., B.L., A.C. and D.M.P. contributed to biostatistical analysis. X.C., C.R.B., B.L., D.K.B. and A.C. contributed to bioinformatics analysis; X.C., A.R.P., S.J., A.R., A.L.K. and W.I.L contributed to the data interpretation and functional analysis; X.C., C.R.B., A.R.P., Y.Y., S.J., A.R., B.L., D.K.B., A.C., D.M.P., S.L., D.L.P., S.D.V., L.B., M.H., J.G.M., A.L.K., O.F. and W.I.L. contributed to the writing and the review of manuscript before submission for publication. All authors read and approved the final manuscript.

## Acknowledgements

This study was funded by National Institutes of Health U54 AI138370 (Center for Solutions for ME/CFS). We are grateful to Kelly Magnus of the Center for Infection and Immunity at Columbia University for assistance with manuscript preparation, to Kelly Paglia of the West Coast Metabolomics Center at University of California, Davis for her support and coordination, to the Chronic Fatigue Initiative of the Hutchins Family Foundation and the ME/CFS patients who provided the samples and inspiration that enabled our work. We dedicate this paper to the memory of Bohyun Lee and her contributions to research in ME/CFS.

*   Received June 14, 2021.
*   Revision received June 14, 2021.
*   Accepted June 22, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Committee on the Diagnostic Criteria for Myalgic Encephalomyelitis/Chronic Fatigue Syndrome; Board on the Health of Select Populations; Institute of Medicine Beyond Myalgic Encephalomyelitis/Chronic Fatigue Syndrome: Redefining an Illness. National Academdies Press (2015).
    
    

2.  2.Carruthers BM, et al. Myalgic encephalomyelitis: International Consensus Criteria. J Intern Med 270, 327–338 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2796.2011.02428.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21777306&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

3.  3.Jason LA, Mirin AA. Updating the National Academy of Medicine ME/CFS prevalence and economic impact figures to account for population growth and inflation. Fatigue: Biomedicine, Health & Behavior, (2021).
    
    

4.  4.Haney E, et al. Diagnostic Methods for Myalgic Encephalomyelitis/Chronic Fatigue Syndrome: A Systematic Review for a National Institutes of Health Pathways to Prevention Workshop. Ann Intern Med 162, 834–840 (2015).
    
    

5.  5.Scheibenbogen C, et al. The European ME/CFS Biomarker Landscape project: an initiative of the European network EUROMENE. J Transl Med 15, 162 (2017).
    
    

6.  6.Armstrong CW, McGregor NR, Sheedy JR, Buttfield I, Butt HL, Gooley PR. NMR metabolic profiling of serum identifies amino acid disturbances in chronic fatigue syndrome. Clin Chim Acta 413, 1525–1531 (2012).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22728138&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

7.  7.Armstrong CW, McGregor NR, Lewis DP, Butt HL, Gooley PR. Metabolic profiling reveals anomalous energy metabolism and oxidative stress pathways in chronic fatigue syndrome patients. Metabolomics 11, 1626–1639 (2015).
    
    

8.  8.Armstrong CW, McGregor NR, Lewis DP, Butt HL, Gooley PR. The association of fecal microbiota and fecal, blood serum and urine metabolites in myalgic encephalomyelitis/chronic fatigue syndrome. Metabolomics 13, (2017).
    
    

9.  9.Fluge O, et al. Metabolic profiling indicates impaired pyruvate dehydrogenase function in myalgic encephalopathy/chronic fatigue syndrome. JCI Insight 1, e89376 (2016).
    
    

10. 10.Germain A, Ruppert D, Levine SM, Hanson MR. Metabolic profiling of a myalgic encephalomyelitis/chronic fatigue syndrome discovery cohort reveals disturbances in fatty acid and lipid metabolism. Mol Biosyst 13, 371–379 (2017).
    
    

11. 11.Germain A, Ruppert D, Levine SM, Hanson MR. Prospective Biomarkers from Plasma Metabolomics of Myalgic Encephalomyelitis/Chronic Fatigue Syndrome Implicate Redox Imbalance in Disease Symptomatology. Metabolites 8, (2018).
    
    

12. 12.Germain A, Barupal DK, Levine SM, Hanson MR. Comprehensive Circulatory Metabolomics in ME/CFS Reveals Disrupted Metabolism of Acyl Lipids and Steroids. Metabolites 10, (2020).
    
    

13. 13.McGregor NR, Armstrong CW, Lewis DP, Gooley PR. Post-Exertional Malaise Is Associated with Hypermetabolism, Hypoacetylation and Purine Metabolism Deregulation in ME/CFS Cases. Diagnostics (Basel) 9, (2019).
    
    

14. 14.Nagy-Szakal D, et al. Insights into myalgic encephalomyelitis/chronic fatigue syndrome phenotypes through comprehensive metabolomics. Sci Rep 8, 10056 (2018).
    
    

15. 15.Naviaux RK, et al. Metabolic features of chronic fatigue syndrome. Proc Natl Acad Sci U S A 113, E5472–5480 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTEzLzM3L0U1NDcyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDYvMjIvMjAyMS4wNi4xNC4yMTI1ODg5NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

16. 16.Yamano E, et al. Index markers of chronic fatigue syndrome with dysfunction of TCA and urea cycles. Sci Rep 6, 34990 (2016).
    
    

17. 17.Valdez AR, et al. Estimating Prevalence, Demographics, and Costs of ME/CFS Using Large Scale Medical Claims Data and Machine Learning. Front Pediatr 6, 412 (2018).
    
    

18. 18.Tomic S, Brkic S, Maric D, Mikic AN. Lipid and protein oxidation in female patients with chronic fatigue syndrome. Arch Med Sci 8, 886–891 (2012).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23185200&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

19. 19.Milivojevic M, et al. Plasma proteomic profiling suggests an association between antigen driven clonal B cell expansion and ME/CFS. PLoS One 15, e0236148 (2020).
    
    

20. 20.Aaron LA, et al. Comorbid clinical conditions in chronic fatigue: a co-twin control study. J Gen Intern Med 16, 24–31 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1525-1497.2001.03419.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11251747&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000166979600004&link_type=ISI) 

21. 21.Giloteaux L, Goodrich JK, Walters WA, Levine SM, Ley RE, Hanson MR. Reduced diversity and altered composition of the gut microbiome in individuals with myalgic encephalomyelitis/chronic fatigue syndrome. Microbiome 4, 30 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40168-016-0171-4&link_type=DOI) 

22. 22.Maes M, Bosmans E, Kubera M. Increased expression of activation antigens on CD8+ T lymphocytes in Myalgic Encephalomyelitis/chronic fatigue syndrome: inverse associations with lowered CD19+ expression and CD4+/CD8+ ratio, but no associations with (auto)immune, leaky gut, oxidative and nitrosative stress biomarkers. Neuro Endocrinol Lett 36, 439–446 (2015).
    
    

23. 23.Nagy-Szakal D, et al. Fecal metagenomic profiles in subgroups of patients with myalgic encephalomyelitis/chronic fatigue syndrome. Microbiome 5, 44 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40168-017-0261-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28095919&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

24. 24.Fukuda K, Straus SE, Hickie I, Sharpe MC, Dobbins JG, Komaroff A. The chronic fatigue syndrome: a comprehensive approach to its definition and study. International Chronic Fatigue Syndrome Study Group. Ann Intern Med 121, 953–959 (1994).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/0003-4819-121-12-199412150-00009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7978722&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1994PW08300009&link_type=ISI) 

25. 25.Carruthers BM, et al. Myalgic Encephalomyelitis/Chronic Fatigue Syndrome. Journal of Chronic Fatigue Syndrome 11, 7–115 (2003).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1300/J092v11n01_02&link_type=DOI) 

26. 26.Jason LA, et al. The Development of a Revised Canadian Myalgic Encephalomyelitis Chronic Fatigue Syndrome Case Definition. American Journal of Biochemistry and Biotechnology 6, 120–135 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3844/ajbbsp.2010.120.135&link_type=DOI) 

27. 27.Buysse DJ, Reynolds CF, 3rd, Monk TH, Berman SR, Kupfer DJ. The Pittsburgh Sleep Quality Index: a new instrument for psychiatric practice and research. Psychiatry Res 28, 193–213 (1989).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0165-1781(89)90047-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2748771&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1989AB69000008&link_type=ISI) 

28. 28.Ware JE, Jr.., Sherbourne CD. The MOS 36-item short-form health survey (SF-36). I. Conceptual framework and item selection. Med Care 30, 473–483 (1992).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00005650-199206000-00002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1593914&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992HX94800002&link_type=ISI) 

29. 29.Smets EM, Garssen B, Bonke B, De Haes JC. The Multidimensional Fatigue Inventory (MFI) psychometric qualities of an instrument to assess fatigue. J Psychosom Res 39, 315–325 (1995).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0022-3999(94)00125-O&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7636775&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QX47500009&link_type=ISI) 

30. 30.Fiehn O. Metabolomics by Gas Chromatography-Mass Spectrometry: Combined Targeted and Untargeted Profiling. Curr Protoc Mol Biol 114, 30 34 31–30 34 32 (2016).
    
    

31. 31.Kind T, et al. FiehnLib: mass spectral and retention index libraries for metabolomics based on quadrupole and time-of-flight gas chromatography/mass spectrometry. Anal Chem 81, 10038–10048 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/ac9019522&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19928838&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

32. 32.Cajka T, Smilowitz JT, Fiehn O. Validating Quantitative Untargeted Lipidomics Across Nine Liquid Chromatography-High-Resolution Mass Spectrometry Platforms. Anal Chem 89, 12360–12368 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/acs.analchem.7b03404&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29064229&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

33. 33.Tsugawa H, et al. MS-DIAL: data-independent MS/MS deconvolution for comprehensive metabolome analysis. Nat Methods 12, 523–526 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.3393&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25938372&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

34. 34.Kind T, Liu KH, Lee DY, DeFelice B, Meissen JK, Fiehn O. LipidBlast in silico tandem mass spectrometry database for lipid identification. Nat Methods 10, 755–758 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.2551&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23817071&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000322453600025&link_type=ISI) 

35. 35.Bakovic M, Fullerton MD, Michel V. Metabolic and molecular aspects of ethanolamine phospholipid biosynthesis: the role of CTP:phosphoethanolamine cytidylyltransferase (Pcyt2). Biochem Cell Biol 85, 283–300 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1139/O07-006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17612623&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

36. 36.DeFelice BC, et al. Mass Spectral Feature List Optimizer (MS-FLO): A Tool To Minimize False Positive Peak Reports in Untargeted Liquid Chromatography-Mass Spectroscopy (LC-MS) Data Processing. Anal Chem 89, 3250–3255 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/acs.analchem.6b04372&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28225594&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

37. 37.Fan S, et al. Systematic Error Removal Using Random Forest for Normalizing Large-Scale Untargeted Lipidomics Data. Anal Chem 91, 3590–3596 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/acs.analchem.8b05592&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30758187&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

38. 38.Benjamini Y, Hochberg Y. Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B-Methodological 57, 289–300 (1995).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=ISI:A1995QE4&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

39. 39.Barupal DK, Fiehn O. Chemical Similarity Enrichment Analysis (ChemRICH) as alternative to biochemical pathway mapping for metabolomic datasets. Sci Rep 7, 14567 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-017-15231-w&link_type=DOI) 

40. 40.Goodrich B, Gabry J, Ali I, Brilleman S. rstanarm: Bayesian applied regression modeling via Stan.) (2020).
    
    

41. 41.Makowski D, Ben-Shachar MS, Lüdecke D. bayestestR: Describing Effects and their Uncertainty, Existence and Significance within the Bayseian Framework. Journal of Open Source Software 4, 1541 (2019).
    
    

42. 42.Makowski D, Ben-Shachar MS, Chen SHA, Ludecke D. Indices of Effect Existence and Significance in the Bayesian Framework. Front Psychol 10, 2767 (2019).
    
    

43. 43.Jeffreys H. Theory of Probability, 3rd edn. Clarendon Press (1961).
    
    

44. 44.Tibshirani R. Regression Shrinkage and Selection Via the Lasso. Journal of the Royal Statistical Society: Series B (Methodological) 58, 267–288 (1996).
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996TU31400017&link_type=ISI) 

45. 45.Zou H. The Adaptive Lasso and Its Oracle Properties. Journal of the American Statistical Association 101, 1418–1429 (2012).
    
    

46. 46.Breiman L. Random Forests. Machine Learning 45, 5–32 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1023/A:1010933404324&link_type=DOI) 

47. 47.Chen T, Guestrin C. XGBoost: A Scalable Tree Boosting System. In: 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining) (2016).
    
    

48. 48.Fan J, Li R. Variable Selection via Nonconcave Penalized Likelihood and its Oracle Properties. Journal of the American Statistical Association 96, 1348–1360 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1198/016214501753382273&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000172728000028&link_type=ISI) 

49. 49.Hoeting JA, Madigan D, Raftery AE, Volinsky CT. Bayesian Model Averaging: A Tutorial. Statistical Science 14, 382–401 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1214/ss/1009212519&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00008609&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000086094500003&link_type=ISI) 

50. 50.Gibellini F, Smith TK. The Kennedy pathway--De novo synthesis of phosphatidylethanolamine and phosphatidylcholine. IUBMB Life 62, 414–428 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/iub.337&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20503434&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000278980400002&link_type=ISI) 

51. 51.Sperka-Gottlieb CD, Hermetter A, Paltauf F, Daum G. Lipid topology and physical properties of the outer mitochondrial membrane of the yeast, Saccharomyces cerevisiae. Biochim Biophys Acta 946, 227–234 (1988).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0005-2736(88)90397-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3061466&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

52. 52.Zinser E, Sperka-Gottlieb CD, Fasch EV, Kohlwein SD, Paltauf F, Daum G. Phospholipid synthesis and lipid composition of subcellular membranes in the unicellular eukaryote Saccharomyces cerevisiae. J Bacteriol 173, 2026–2034 (1991).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MjoiamIiO3M6NToicmVzaWQiO3M6MTA6IjE3My82LzIwMjYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNi8yMi8yMDIxLjA2LjE0LjIxMjU4ODk1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

53. 53.Calzada E, et al. Phosphatidylethanolamine made in the inner mitochondrial membrane is essential for yeast cytochrome bc1 complex function. Nat Commun 10, 1432 (2019).
    
    

54. 54.Furt F, Moreau P. Importance of lipid metabolism for intracellular and mitochondrial membrane fusion/fission processes. Int J Biochem Cell Biol 41, 1828–1836 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biocel.2009.02.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19703652&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

55. 55.Birner R, Burgermeister M, Schneiter R, Daum G. Roles of phosphatidylethanolamine and of its several biosynthetic pathways in Saccharomyces cerevisiae. Mol Biol Cell 12, 997–1007 (2001).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6Im1vbGJpb2xjZWxsIjtzOjU6InJlc2lkIjtzOjg6IjEyLzQvOTk3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDYvMjIvMjAyMS4wNi4xNC4yMTI1ODg5NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

56. 56.Dowhan W, Bogdanov M. Lipid-dependent membrane protein topogenesis. Annu Rev Biochem 78, 515–540 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.biochem.77.060806.091251&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19489728&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268069200019&link_type=ISI) 

57. 57.Schuler MH, Di Bartolomeo F, Martensson CU, Daum G, Becker T. Phosphatidylcholine Affects Inner Membrane Protein Translocases of Mitochondria. J Biol Chem 291, 18718–18729 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyOTEvMzYvMTg3MTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNi8yMi8yMDIxLjA2LjE0LjIxMjU4ODk1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

58. 58.Gebert N, et al. Mitochondrial cardiolipin involved in outer-membrane protein biogenesis: implications for Barth syndrome. Curr Biol 19, 2133–2139 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cub.2009.10.074&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19962311&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

59. 59.Becker T, Horvath SE, Bottinger L, Gebert N, Daum G, Pfanner N. Role of phosphatidylethanolamine in the biogenesis of mitochondrial outer membrane proteins. J Biol Chem 288, 16451–16459 (2013).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyODgvMjMvMTY0NTEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNi8yMi8yMDIxLjA2LjE0LjIxMjU4ODk1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

60. 60.Tasseva G, Bai HD, Davidescu M, Haromy A, Michelakis E, Vance JE. Phosphatidylethanolamine deficiency in Mammalian mitochondria impairs oxidative phosphorylation and alters mitochondrial morphology. J Biol Chem 288, 4158–4173 (2013).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEwOiIyODgvNi80MTU4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDYvMjIvMjAyMS4wNi4xNC4yMTI1ODg5NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

61. 61.Bottinger L, et al. Phosphatidylethanolamine and cardiolipin differentially affect the stability of mitochondrial respiratory chain supercomplexes. J Mol Biol 423, 677–686 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jmb.2012.09.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22971339&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

62. 62.Yoda E, et al. Mitochondrial dysfunction and reduced prostaglandin synthesis in skeletal muscle of Group VIB Ca2+-independent phospholipase A2gamma-deficient mice. J Lipid Res 51, 3003–3015 (2010).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamxyIjtzOjU6InJlc2lkIjtzOjEwOiI1MS8xMC8zMDAzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDYvMjIvMjAyMS4wNi4xNC4yMTI1ODg5NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

63. 63.Dawaliby R, et al. Allosteric regulation of G protein-coupled receptor activity by phospholipids. Nat Chem Biol 12, 35–39 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nchembio.1960&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26571351&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

64. 64.Momchilova A, Markovska T. Phosphatidylethanolamine and phosphatidylcholine are sources of diacylglycerol in ras-transformed NIH 3T3 fibroblasts. Int J Biochem Cell Biol 31, 311–318 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1357-2725(98)00111-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10216963&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000079119300007&link_type=ISI) 

65. 65.Okamoto Y, Morishita J, Tsuboi K, Tonai T, Ueda N. Molecular characterization of a phospholipase D generating anandamide and its congeners. J Biol Chem 279, 5298–5305 (2004).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEwOiIyNzkvNy81Mjk4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDYvMjIvMjAyMS4wNi4xNC4yMTI1ODg5NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

66. 66.Li Z, Vance DE. Phosphatidylcholine and choline homeostasis. J Lipid Res 49, 1187–1194 (2008).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamxyIjtzOjU6InJlc2lkIjtzOjk6IjQ5LzYvMTE4NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA2LzIyLzIwMjEuMDYuMTQuMjEyNTg4OTUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

67. 67.Mathias S, Kolesnick R. Ceramide: a novel second messenger. Adv Lipid Res 25, 65–90 (1993).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8368154&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1993BZ68D00004&link_type=ISI) 

68. 68.Gerlach BD, et al. Resolvin D1 promotes the targeting and clearance of necroptotic cells. Cell Death Differ 27, 525–539 (2020).
    
    

69. 69.Plioplys AV, Plioplys S. Serum levels of carnitine in chronic fatigue syndrome: clinical correlates. Neuropsychobiology 32, 132–138 (1995).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8544970&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995TB04400005&link_type=ISI) 

70. 70.Flanagan JL, Simmons PA, Vehige J, Willcox MD, Garrett Q. Role of carnitine in disease. Nutr Metab (Lond) 7, 30 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1743-7075-7-30&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20398344&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

71. 71.Li JL, Wang QY, Luan HY, Kang ZC, Wang CB. Effects of L-carnitine against oxidative stress in human hepatocytes: involvement of peroxisome proliferator-activated receptor alpha. J Biomed Sci 19, 32 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1423-0127-19-32&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22435679&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

72. 72.Vacha GM, Giorcelli G, Siliprandi N, Corsi M. Favorable effects of L-carnitine treatment on hypertriglyceridemia in hemodialysis patients: decisive role of low levels of high-density lipoprotein-cholesterol. Am J Clin Nutr 38, 532–540 (1983).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo4OiIzOC80LzUzMiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA2LzIyLzIwMjEuMDYuMTQuMjEyNTg4OTUuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

73. 73.Martinez-Reyes I, Chandel NS. Mitochondrial TCA cycle metabolites control physiology and disease. Nat Commun 11, 102 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-13668-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31900386&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

74. 74.Zhang Y, et al. Succinate accumulation induces mitochondrial reactive oxygen species generation and promotes status epilepticus in the kainic acid rat model. Redox Biol 28, 101365 (2020).
    
    

75. 75.Palmieri EM, et al. Nitric oxide orchestrates metabolic rewiring in M1 macrophages by targeting aconitase 2 and pyruvate dehydrogenase. Nat Commun 11, 698 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-14433-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32019928&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

76. 76.Morris G, Berk M, Klein H, Walder K, Galecki P, Maes M. Nitrosative Stress, Hypernitrosylation, and Autoimmune Responses to Nitrosylated Proteins: New Pathways in Neuroprogressive Disorders Including Depression and Chronic Fatigue Syndrome. Mol Neurobiol 54, 4271–4291 (2017).
    
    

77. 77.Morris G, Maes M. Mitochondrial dysfunctions in myalgic encephalomyelitis/chronic fatigue syndrome explained by activated immuno-inflammatory, oxidative and nitrosative stress pathways. Metab Brain Dis 29, 19–36 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11011-013-9435-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24557875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

78. 78.Cabral-Marques O, et al. GPCR-specific autoantibody signatures are associated with physiological and pathological immune homeostasis. Nat Commun 9, 5224 (2018).
    
    

79. 79.Loebel M, et al. Antibodies to beta adrenergic and muscarinic cholinergic receptors in patients with Chronic Fatigue Syndrome. Brain Behav Immun 52, 32–39 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbi.2015.09.013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26399744&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

80. 80.Wirth K, Scheibenbogen C. A Unifying Hypothesis of the Pathophysiology of Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS): Recognitions from the finding of autoantibodies against ss2-adrenergic receptors. Autoimmun Rev 19, 102527 (2020).
    
    

81. 81.Capuron L, et al. Cognitive dysfunction relates to subjective report of mental fatigue in patients with chronic fatigue syndrome. Neuropsychopharmacology 31, 1777–1784 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.npp.1301005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16395303&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000239341200022&link_type=ISI) 

82. 82.Cockshell SJ, Mathias JL. Cognitive functioning in chronic fatigue syndrome: a meta-analysis. Psychol Med 40, 1253–1267 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291709992054&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20047703&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

83. 83.Han X, et al. Metabolomics in early Alzheimer’s disease: identification of altered plasma sphingolipidome using shotgun lipidomics. PLoS One 6, e21643 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0021643&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21779331&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

84. 84.Torretta E, et al. Particular CSF sphingolipid patterns identify iNPH and AD patients. Sci Rep 8, 13639 (2018).
    
    

85. 85.Stoessel D, et al. Promising Metabolite Profiles in the Plasma and CSF of Early Clinical Parkinson’s Disease. Front Aging Neurosci 10, 51 (2018).
    
    

86. 86.Pieragostino D, et al. Enhanced release of acid sphingomyelinase-enriched exosomes generates a lipidomics signature in CSF of Multiple Sclerosis patients. Sci Rep 8, 3071 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-018-21497-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29449691&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 

87. 87.Pujol-Lereis LM. Alteration of Sphingolipids in Biofluids: Implications for Neurodegenerative Diseases. Int J Mol Sci 20, (2019).
    
    

88. 88.Furuya S, Mitoma J, Makino A, Hirabayashi Y. Ceramide and its interconvertible metabolite sphingosine function as indispensable lipid factors involved in survival and dendritic differentiation of cerebellar Purkinje cells. J Neurochem 71, 366–377 (1998).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9648886&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000074283600041&link_type=ISI) 

89. 89.Schwartz NU, et al. Decreased ceramide underlies mitochondrial dysfunction in Charcot-Marie-Tooth 2F. FASEB J 32, 1716–1728 (2018).
    
    

90. 90.Guan Z, Wang Y, Cairns NJ, Lantos PL, Dallner G, Sindelar PJ. Decrease and structural modifications of phosphatidylethanolamine plasmalogen in the brain with Alzheimer disease. J Neuropathol Exp Neurol 58, 740–747 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00005072-199907000-00008&link_type=DOI) 

91. 91.Ginsberg L, Xuereb JH, Gershfeld NL. Membrane instability, plasmalogen content, and Alzheimer’s disease. J Neurochem 70, 2533–2538 (1998).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9603219&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000073647800033&link_type=ISI) 

92. 92.Han X, Holtzman DM, McKeel DW, Jr.. Plasmalogen deficiency in early Alzheimer’s disease subjects and in animal models: molecular characterization using electrospray ionization mass spectrometry. J Neurochem 77, 1168–1180 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1471-4159.2001.00332.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11359882&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000168771600022&link_type=ISI) 

93. 93.Farooqui AA, Horrocks LA. Plasmalogens: workhorse lipids of membranes in normal and injured neurons and glia. Neuroscientist 7, 232–245 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/107385840100700308&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11499402&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000169011900017&link_type=ISI) 

94. 94.Perichon R, Moser AB, Wallace WC, Cunningham SC, Roth GS, Moser HW. Peroxisomal disease cell lines with cellular plasmalogen deficiency have impaired muscarinic cholinergic signal transduction activity and amyloid precursor protein secretion. Biochem Biophys Res Commun 248, 57–61 (1998).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1006/bbrc.1998.8909&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9675085&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000074758800011&link_type=ISI) 

95. 95.Ifuku M, et al. Anti-inflammatory/anti-amyloidogenic effects of plasmalogens in lipopolysaccharide-induced neuroinflammation in adult mice. J Neuroinflammation 9, 197 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1742-2094-9-197&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22889165&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F22%2F2021.06.14.21258895.atom)