Skip to main content

Sex-specific differences in plasma lipid profiles are associated with Gulf War Illness



Nearly 250,000 veterans from the 1990–1991 Gulf War have Gulf War Illness (GWI), a condition with heterogeneous pathobiology that remains difficult to diagnose. As such, availability of blood biomarkers that reflect the underlying biology of GWI would help clinicians provide appropriate care to ill veterans. In this study, we measured blood lipids to examine the influence of sex on the association between blood lipids and GWI diagnosis.


Plasma lipid extracts from GWI (n = 100) and control (n = 45) participants were subjected to reversed-phase nano-flow liquid chromatography-mass spectrometry analysis.


An influence of sex and GWI case status on plasma neutral lipid and phospholipid species was observed. Among male participants, triglycerides, diglycerides, and phosphatidylcholines were increased while cholesterol esters were decreased in GWI cases compared to controls. In female participants, ceramides were increased in GWI cases compared to controls. Among male participants, unsaturated triglycerides, phosphatidylcholine and diglycerides were increased while unsaturated cholesterol esters were lower in GWI cases compared to controls. The ratio of arachidonic acid- to docosahexaenoic acid-containing triglyceride species was increased in female and male GWI cases as compared to their sex-matched controls.


Differential modulation of neutral lipids and ratios of arachidonic acid to docosahexaenoic acid in male veterans with GWI suggest metabolic dysfunction and inflammation. Increases in ceramides among female veterans with GWI also suggest activation of inflammatory pathways. Future research should characterize how these lipids and their associated pathways relate to GWI pathology to identify biomarkers of the disorder.


Approximately 30% of the 700,000 veterans of the 1991 Gulf War (GW) suffer from Gulf War Illness (GWI), a chronic multi-symptom condition characterized by cognitive impairment, fatigue, and pain [1]. Veterans’ exposure to chemical warfare agents, pesticides, and pyridostigmine bromide are among the major contributors to GWI etiology [1,2,3,4]. Currently, there are no biomarkers for indicating the underlying pathology of this condition. Advances in GWI research suggest lipid dysfunction, bioenergetics deficits, and inflammation as key contributors to its ongoing pathology [5, 6]. Blood biomarkers of these biological processes would be a cost-effective and minimally invasive avenue to assist clinicians with accurate detection of GWI.

Evidence suggests that GW pesticides affect cellular lipid metabolism and lipid homeostasis [2, 7,8,9,10,11]. Altered lipid homeostasis can result in an abnormal accumulation of lipids that promote cellular inflammation [12, 13]. Abnormal plasma and brain lipid profiles after GW pesticide exposure in rodents and in veterans with GWI have been shown to accompany inflammation, neurobehavioral and bioenergetic deficits [14,15,16,17,18]. As such, altered plasma lipid profiles may provide a biomarker signature for the bioenergetic deficits and inflammation associated with GWI.

It is expected that since the pesticides to which GW veterans were exposed tend to disturb lipid homeostasis, owing to their high affinity for lipids, neutral lipids classes, such as triglycerides (TG), diglycerides (DG), and cholesterol esters (CE), and phospholipid (PL) classes would be affected in GWI [2, 9,10,11, 19, 20]. The degree of unsaturation of the fatty acids within the PL classes is known to alter metabolic and inflammatory systems, with many saturated fatty acid (SFA)-containing PL positively associated with chronic metabolic conditions whereas polyunsaturated fatty acid (PUFA)-containing PL are associated with inflammation [21,22,23,24]. Among these, PUFA containing the omega-3 fatty acid, docosahexaenoic acid (DHA), generate anti-inflammatory or inflammation-resolving lipid metabolites and PUFA containing the omega-6 fatty acid, arachidonic acid (AA), can produce metabolites that promote pro-inflammatory responses [23,24,25].

Changes in lipid profiles were also expected in relation to the sex of GW veterans, as many classes of plasma lipids are significantly associated with sex [26]. Women are up to four times more likely than men to be diagnosed with myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS), a chronic multi-symptom disease clinically similar to GWI that has been shown to have sex-linked differences in blood lipid profiles [27, 28]. Female GW veterans have a higher prevalence of GWI diagnosis and severe GWI-related symptoms than male GW veterans [29]. These findings suggest potential sex-linked differences in pathophysiological responses to GW-era chemical exposures, necessitating sex-stratified analyses.

We hypothesized that blood lipid levels would be affected both by GWI status and the sex of GW veterans. We evaluated the levels of several major blood lipid classes and their molecular species in three cross-sectional cohorts of participants classified as GWI or control. This study will further our understanding of plasma lipids as potential biological markers to track the ongoing inflammation and metabolic dysfunction of GWI, which may help to provide better care and management of this illness.

Materials and methods

Cohort characteristics

All protocols were conducted in compliance with the relevant guidelines and regulations and were approved by institutional review boards (IRB) at each institution. Informed consent was obtained from all participants. All plasma utilized within this study were from the following case–control veteran cohorts: (1) the Roskamp Neurology Clinic (RNC) cohort (n = 36), (2) the Boston Gulf War Illness Consortium (GWIC), and (3) the Fort Devens cohort (n = 58) [3, 30, 31]. Additional healthy control samples were provided by Dr. Klimas to supplement the control population of the study. GW veteran controls and the healthy civilian controls were grouped together for analysis in the GWIC cohort (n = 51). In all cohorts, either the Center for Disease Control (CDC) chronic multi-symptom illness definition or the Kansas GWI criteria were used to diagnose GWI [32, 33]. These criteria required that GW veterans must show symptoms in at least 3 of 6 symptom domains (fatigue/sleep problems, somatic pain, neurological/cognitive/mood symptoms, gastrointestinal symptoms, respiratory symptoms, and skin abnormalities) for GWI diagnosis. Exclusionary criteria have been described previously, primarily that participants were excluded if they reported prior diagnosis with a confounding medical condition [17, 30, 31]. Controls were either veterans deployed to the 1991 GW or healthy civilian controls, who did not meet the Kansas GWI criteria or any of the exclusionary criteria. Additional information about the self-report of cardiovascular health status was collected on RNC and GWIC participants. Self-report of any heart conditions, diabetes, high blood pressure or high cholesterol was combined into a single variable to assess cardiovascular risk factors as a group. Reporting of any one or more of these risk factors was coded as presence and a lack of reporting on any cardiovascular risk factors was coded as absence.

Data on the sex of the participants were gathered via self-report, with the options being either male or female. As per the National Institute of Health, sex refers to biological differences between males and females [34]. As per the World Health Organization, gender refers to the characteristics of men and women that are socially constructed [35]. Therefore, this paper evaluating biological lipid levels uses exclusively male/female terminology for participants rather than men/women terminology, as data on the gender of the participants was not collected.

Participant recruitment, blood collection, and blood storage procedures have been described previously [17, 18, 31]. All samples were collected using the same written standard operating procedures for performing phlebotomy, plasma separation and aliquoting. All participants had fasted prior to blood collection. There was no requirement of time of collection.

Due to low numbers in ethnic groups other than Caucasian, we combined participants who were African American, Aleutian Eskimo, American Indian, Asian, Pacific Islander, multiracial or other races into one category, called Non-Caucasian. Only one African American veteran with GWI was present in the RNC cohort. Of the eleven African American participants in GWIC, three were control participants and eight were veterans diagnosed with GWI. Of the six African American participants in the Fort Devens cohort, four were control participants and two were veterans diagnosed with GWI.

Lipidomic analysis and quantification

Plasma samples were randomized prior to use via a group-block design, to ensure that the experimenters were blinded to participants’ group membership and other characteristics. The extraction protocol and liquid chromatography–mass spectrometry (LC–MS) run settings have been previously described by our laboratory [18, 36]. Briefly, lipid extracts from 10 µL of plasma spiked with class-specific internal standards (to control for variability in lipid extractions) were resuspended in 200 µL of 70:30 LC–MS mobile phase A:B [36, 37]. Mobile phase A was composed of 42% water, 31% acetonitrile, 27% isopropanol, 10 mM ammonium formate, and 0.1% formic acid. Mobile phase B was composed of 90% isopropanol, 10% acetonitrile, 10 mM ammonium formate, and 0.1% formic acid. For LC–MS analysis, an Easy-nanoLC 1000 chromatography pumping system was used with a nanoflex electrospray ionization (ESI) source on a Thermo LTQ/orbitrap mass spectrometer.

Resuspended lipid extracts (4 µL) were injected onto an Acclaim PepMapTM 100 (75 µm × 2 cm, nanoViper, C18, 3 µm, 100 Å) trapping column attached to an Acclaim PepMapTM RSLC (75 µm × 15 cm, nanoViper, C18, 2 µm, 100 Å) analytical column for reversed-phase chromatographic separation of lipid species. The chromatographic gradient started at 30% B, then increased linearly from 1 to 40 min from 50 to 98% B, plateauing at 98% B from 40 to 50 min, and then decreasing to 30% B for 50 to 65 min. All samples were run in triplicate, accompanied by a blank and a quality control sample in each run. Full-scan mass spectrometry (MS) data was acquired in positive ion mode, over a mass range of 300–2000 m/z at a resolution of 30,000. A maximum inject time of 200 ms was used with 6 microscans/acquired scan. Thermo Xcalibur Qual Browser was used to examine the individual chromatographic data of samples for lipid retention times. Lipid peak areas were integrated using Tracefinder™ software, using target compound lists of analytes of interest, identifying them based on m/z and retention time. A high mass accuracy window of 5 mmu was used for detecting each individual lipid species. We detected several major classes of lipids including the phospholipid classes phosphatidylcholine (PC), phosphatidylethanolamine (PE), the neutral lipid classes TG, DG, and CE and the sphingolipid classes sphingomyelin (SM), ceramides (Cer) and hexosylceramides (HexCer), which were quantified as previously described [28, 36]. The quality control (QC) samples were then used to normalize between batches of data ([sample] x [batch QC/normalizing QC]). Once the coefficients of variance (CV) were calculated per class for the whole cohort, lipid species with an average CV > 25 were excluded from further analysis. Lipid species with less than 90% coverage were also excluded from our analyses.

To categorize lipids based on their degree of unsaturation, totals were calculated by adding up values for all lipid species measured belonging to that lipid class. Saturated fatty acids (SFA) were calculated by totaling the lipids with no double bonds; monounsaturated fatty acids (MUFA) were calculated by totaling the lipids with one double bond. Polyunsaturated fatty acids (PUFA) were calculated by totaling all measured lipid species with three or more double bonds. For lipid classes known to contain AA or DHA, their respective categorizations were made by totaling individual species known to contain the relevant type of fatty acid, as based on previous MS/MS confirmations.

Statistical analysis

Chi-square and one-way analysis of variance (ANOVA) tests were used, as appropriate, to compare between the three cohorts and to compare GWI and healthy controls for descriptive purposes. Each lipid species was assessed for normality and those found to be skewed were log transformed for parametric analysis. All lipidomics data were analyzed using mixed linear modeling (MLM) as described previously to examine the independent effects of sex and GWI diagnosis as well as to account for random effects associated with repeated measurements in lipidomics and metabolomics datasets [17, 36,37,38,39]. The Benjamini-Hochberg (B-H) correction procedure was performed to correct for multiple testing and the false discovery rate (FDR) limit was set at 0.05. The webtool MetaboAnalyst 5.0 was used to perform a hierarchical clustering analysis of the individual lipids with a coverage superior to 90% [40]. Clusters were determined based on the similarity of lipids as expressed by the Euclidean distance. Lipids that were adjacent based on the Ward algorithm were combined in clusters and an ANOVA was applied to select the top 50 lipids significantly modulated by sex and diagnosis represented in the dendrogram and heatmap.


Basic demographics

General demographic characteristics of participants from the three cohorts at the time of the blood collection are presented in Table 1. The difference in ages between the RNC participants and the GWIC participants was found to be significantly different (p < 0.01). The Fort Devens and GWIC cohorts were found to be significantly different for race (p < 0.05). Overall, the GWIC cohort was younger and had a higher proportion of non-Caucasian participants than the other two cohorts (Table 1). The percentage of female participants ranged from 13.7% to 27.8% across the 3 cohorts (Table 1).

Table 1 Descriptive table of demographic variables of the three cohorts: RNC, GWIC, and Fort Devens

When all cohorts were combined (Table 2), there were no significant differences for age and race between GWI cases and controls. However, the sex ratio was found to be significantly different between GWI cases and controls, as 16% of GWI cases were female compared to 31% of controls (p < 0.05). In a subset of participants, there were no differences between GWI cases and controls for BMI, the presence of diabetes and the presence of cardiovascular risk factors. There were no differences between males and females for the presence of diabetes or cardiovascular risk factors (p > 0.05). Similarly there were no significant differences between males and females for BMI (males: 30.1 ± 5.3 SD and females: 29.4 ± 4.2SD). Given the known associations between sex and blood lipid levels and differences in the proportion of male and female sex between GWI and control participants, subsequent analyses were stratified by sex to account for its confounding effects.

Table 2 Descriptive table of demographic variables of GWI cases and controls from the RNC, GWIC, and Fort Devens cohorts combined

Sex-specific differences were observed for the degree of unsaturation and AA and DHA content of several major classes of lipids in veterans with GWI

Figure 1 shows heatmaps of plasma profiles of eight different lipid classes that are differentially affected by sex and diagnostic differences among male and female participants (see Additional file 1: Fig. S1 for individual values). The total amount of PC, DG and TG were significantly increased while the total amount of CE was significantly decreased in male GWI compared to male controls (Fig. 1A). Among female participants, only total Cer were significantly increased in GWI compared to controls. We also evaluated the degree of fatty acid unsaturation in individual lipid species since neutral lipids are considered reservoirs of PUFA lipids (Fig. 1B). Compared to male controls, male veterans with GWI had a significant increase in the unsaturated lipids, consisting of DG-PUFA, TG-MUFA, TG-PUFA and PC-MUFA, while CE-PUFA was significantly decreased. Among female participants, TG-SFA levels were significantly lower in GWI compared to controls. The PUFA-containing lipid species were further grouped by their AA or DHA-content for each lipid class (Fig. 1C). Levels of DG-AA, DG-DHA, TG-AA and TG-DHA were significantly increased while CE-AA and CE-DHA were significantly decreased in male GWI cases compared to male controls. The ratios of PE-AA/PE-DHA and TG-AA/TG-DHA were both significantly increased in male GWI compared to male controls. The TG-AA/TG-DHA ratio was also significantly increased in female GWI compared to female controls. Independent effects of the presence of cardiovascular risk factors on these lipid profiles are presented in Additional file 2: Fig. S2.

Fig. 1
figure 1

Heatmap visualization of phospholipids and neutral lipids in plasma of GWI cases and controls. The heatmaps represent the average concentrations relative to male controls of the A total level of each lipid classes, B lipids stratified by their degree of unsaturation and C lipids containing arachidonic acid (AA) or docosahexaenoic acid (DHA) as well as their AA-containing lipids/DHA-containing lipids ratio in each group. Black asterisks indicate significant differences between GWI cases and their sex-matched controls (p < 0.05). Red asterisks indicate significant differences between males and females within the same diagnostic group (p < 0.05)

Hierarchical clustering analysis identifies sex-specific lipid species signature of GWI

Hierarchical clustering analyses identified two major clusters (cluster 1 and cluster 2) within the top 50 lipid species significantly modulated by sex and diagnosis (Fig. 2). Lipids within both clusters were affected by sex differences, irrespective of GWI diagnosis. Cluster 1 was primarily composed of the neutral lipid class CE and the sphingolipid class SM and showed a phenotype of increased levels in female GWI and decreased levels in male GWI. Cluster 2 was primarily composed of the neutral lipid classes DG and TG, which were increased in male participants compared to female participants, irrespective of GWI diagnosis. Lipids associated with each cluster further separated into 3 subgroups. The CEs of group A show controls to be generally similar regardless of sex while the lipids differed in concentration in the GWI veterans based on sex. Lipids associated with groups B and E were both principally influenced by sex regardless of the diagnostic status. Group B was composed of an ether PC and SM lipids, while Group E was composed of DGs and TGs. Lipids in group C differed between female GWI and female controls. In group D, neutral lipids of the male and female GWI groups were both increased compared to their sex-matched controls. Both group D and E were composed of the neutral lipids, DG and TG, and notably differed in their range of degree of unsaturation, with lipids in group D consisting of more double bonds and those in group E consisting of fewer double bonds. Group F consisted of 2 lipid species that are atypical of any other group, but do not match each other either. PE34:2 is highly increased in the female controls and decreased in all other groups, while TG 46:0 is highly decreased in the female GWI group and increased in all other groups.

Fig. 2
figure 2

Hierarchical clustering of the top 50 lipids. The relative concentration value of each lipid from the top 50 lipids significantly modulated by sex and diagnosis identified by MetaboAnalyst has been projected onto the heatmap for each group. The color legend identifies in red the lipids with higher concentration and identifies in blue the lipids with lower concentration

Of the top 50 lipids significantly modulated by sex and diagnosis, more lipid species were significantly modulated in male GWI cases (Fig. 3A) than in female GWI cases (Fig. 3B) as compared to their respective sex-matched controls. Among male participants, TG and DG were significantly increased by 1.18-fold at minimum (TG 58:9, log2FC = 0.24) and all CEs were significantly decreased by 0.88-fold at minimum (CE 20:3, log2FC = − 0.19) in GWI veterans versus controls. A few TGs and SMs were significantly increased by a minimum of 1.18-fold change (TG 56:8, log2FC = 0.24) in female GWI participants compared to female controls, except for TG 46:0, which was significantly decreased by 0.37-fold (log2FC = − 1.45).

Fig. 3
figure 3

Fold Change of the top 50 lipids by sex. Relative changes of the top 50 lipids significantly modulated by sex and diagnosis in males (A) and females (B) between GWI participants and GW cases are shown by fold change (FC) on the volcano plots. Each dot represents one individual lipid from the list of top 50 lipids identified by MetaboAnalyst. The red dots indicate the significant changes relative with p-value threshold set at 0.05 (horizontal dotted line)


Previous animal studies have shown that exposure to GW chemicals, such as those to which GW veterans were exposed, disturbs lipid homeostasis and contributes to inflammatory and metabolic dysfunction that are seen in veterans with GWI [7,8,9, 14, 15, 41]. To evaluate the translational significance of prior mechanistic studies, the current study analyzed plasma lipid profiles of veterans with GWI and control participants. Most of the differences pertained to the neutral lipids (CE, TG and DG), where the degree of unsaturation and ratios of AA-to-DHA were changed in GWI compared to controls. Our study shows that male veterans with GWI had significant increases in total neutral lipid content compared to male controls. Female veterans with GWI had elevated levels of Cer and lower levels of SFA TG species as compared to female controls. These sex-specific plasma lipid compositions in GWI veterans require further attention as they may indicate progressive metabolic and inflammatory processes, particularly in male veterans with GWI.

In the plasma, neutral lipids primarily appear in lipoprotein compartments due to their hydrophobic nature. The neutral lipid class of TG are major dietary fats, but they are also synthesized in the liver [42]. Because dietary TGs are digested to free fatty acids and monoglycerides in the intestines, plasma TG levels within lipoprotein particles are reflective of TG synthesis in liver via acylation with DG [43]. In plasma, an excess of TG, such as seen in the male veterans with GWI, is thought to reflect an overproduction of very low-density lipoprotein (VLDL) corresponding with the production of CE-depleted low-density lipoprotein (LDL) particles [42]. This production of CE-depleted LDL may explain the significant decrease in CE in male veterans with GWI. In the current study, female participants, irrespective of their diagnosis, had lower levels of TG and DG compared to their male counterparts. This has also been reported by the Genetics of Lipid-Lowering Drugs and Diet Network (GOLDN) study of plasma lipidome, who attributed this to normal sex-linked lipid differences [26]. In terms of the clinical implications of elevated TG in male veterans with GWI, elevated plasma TG is thought to be a risk factor for cardiovascular disease [44]. Elevated TG and the generation of CE-depleted LDL particles are also the defining features of dyslipidemia associated with insulin resistance and type 2 diabetes mellitus [42].

Veterans with GWI may be more sedentary and exercise less than healthy controls due to their condition, so this TG increase could reflect heightened body weight in veterans with GWI. A positive association between high TG and elevated body mass index (BMI) has been previously reported [45, 46]. As early as 1995, it was observed that veterans with severe GWI had slightly higher BMI than non-cases [33]. However, BMI did not differ between GWI and controls (in a subsample) and therefore unlikely to explain the observed differences in lipid profiles. A recent health survey of GW veterans from the Fort Devens cohort suggested increases in the prevalence of high blood pressure, diabetes, and cardiovascular disease over the last 30 years [4]. The Fort Devens Cohort was also assessed to have significantly higher rates of diabetes in both male and female GW veterans than the general population [47]. However, differences between GWI and controls were not evaluated in this prior study. In our current study, a lack of difference in cardiovascular risk factors between controls and GWI suggested that observed differences in lipid profiles between the two groups were not simply a reflection of the presence of cardiovascular risk factors in this study cohort. As many of the neutral lipid changes may be indicative of subsequent risk of cardiovascular disease with age, further assessement of blood lipids by clinicians may be helpful for monitoring cardiovascular health as veterans with GWI age.

In the current study, SFA-containing TG species were reduced in female veterans with GWI. The role of SFA in promoting cardiovascular disease remains controversial as the source of these SFA-containing lipids and the individual identity of each SFA influences their potential role in promoting cardiovascular disease [21, 48]. As the SFA-TG species measured ranged from TG 46:0 to TG 54:0, and therefore consisted of 46 to 54 carbons split amongst 3 fatty acid chains, it is likely that the SFA-TG species measured contained long-chain fatty acids, which are over 12 carbons in length [48]. However, their significance in the clinical presentation of female GWI is unclear. Since these same TG-SFA species were elevated among those with cardiovascular risk factors, future studies aimed at different aspects of dietary modification in male vs. female GW veterans may be required to better understand the role of SFA-containing TG species in the health, particularly cardiovascular health, of veterans with GWI.

Our previous studies using subsets of samples from these cohorts showed increases in proinflammatory cytokines, interleukin-1 beta (IL-1β), interferon gamma (IFN-γ) and IL-6, and increases of peripheral immune cells among veterans with GWI compared to controls [18]. Others have shown that veterans with GWI have increases in other blood inflammatory markers consisting of C-reactive proteins and matrix-metalloproteinase-2 [49]. Several proinflammatory cytokines, including IL-1β and IL-15, are associated with an increase in GWI symptom severity [50]. As such, an examination of lipids related to inflammation may provide additional information about the underlying biology of GWI. Since omega-6 AA or omega-3 DHA modulate inflammation and comprise a large proportion of PUFA, we focused on lipids that contain AA and DHA as their fatty acid side chains. Our previous studies have shown that the AA and DHA content within certain PL classes were associated with GWI [17]. We now show that ratios of AA-to-DHA were increased in TG species among male and female veterans with GWI. The ratio of AA-to-DHA was also increased in PE among male GWI veterans. As AA and DHA compose 50% and 40% of brain PUFAs respectively, sufficient AA is necessary for the growth, repair, and maintenance of neurons and DHA is used in neurotransmission [51]. However, AA metabolites, such as prostaglandins, are known to activate pro-inflammatory pathways whereas DHA metabolites, such as resolvins, contribute to inflammation resolution pathways [24, 52]. In addition, studies have shown that patients with a low DHA/AA ratio, (i.e. a higher AA-to-DHA ratio) had a higher risk of acute coronary syndrome than those with a high DHA/AA ratio, and this was significant for men in particular [53]. These differences in AA and DHA composition could reflect a generalized pro-inflammatory state in veterans with GWI, which would be consistent with reports of markers of chronic, low-grade inflammation being associated with GWI diagnosis [49, 54]. Hence, future studies are needed that interrogate systems-level omic data along with biological parameters of inflammation to identify specific inflammatory pathways that are associated with GWI pathogenesis.

Another indication of pro-inflammatory states in GWI would be the increase in Cer in female GWI veterans, as ceramides are known to stimulate the production of pro-inflammatory cytokines, such as IL-6 and TNF-α, in macrophages and are correlated with IL-6 levels in serum [55, 56]. Production of long-chain ceramides is known to induce apoptosis [57]. The increase in the long chain Cer measured in this study (d32:0–44:2) may therefore indicate increased cell death in female GWI veterans. Ceramide levels are shown to be increased with obesity and metabolic disorders, particularly diabetes. For instance, Boon and colleagues also found that Cer elevation is another marker of the dyslipidemia of type 2 diabetes mellitus and promotes insulin resistance [55]. Others have shown that plasma Cer increases are associated with obesity and may be indicative of lipotoxicity, which may ultimately contribute to the development of diabetes and insulin resisistance [58, 59]. These studies suggest a different process in males and females in relation to Cer-mediated metabolic dysfunction that warrants further investigation.

This study provides key information about blood biomarkers associated with adverse metabolic and inflammation-related lipid profiles. The BMI in females was non-significantly lower than in males and therefore did not present as a confounding factors for the differences in lipids seen between the two sex groups. These findings are largely consistent with previous studies showing sex-specific differences in lipid profiles. Another limitation of our study includes unavailability of dietary data and medication data, particularly on the use of statins and omega-3/fish oil supplements. Since statins lower TG levels in blood and our findings suggest increases of TG in male veterans with GWI, we anticipate that it is unlikely that statin use affected our study results. Similarly, we have previously shown that AA-to-DHA ratios are decreased after DHA supplementation [39]. Since we observed an increase in AA-to-DHA ratios in male GWI veterans, it is similarly unlikely that omega-3 or fish oil supplementation was a confounder in our study. As CE levels were lower in male GWI veterans, measuring free cholesterol would be useful in future studies to determine if either the level of free cholesterol or the rate of cholesterol esterification had decreased. Unfortunately, we could not reliably detect the labelled free cholesterol in our assay but will correct this in future studies. While we have previously examined free fatty acids in a smaller cohort of GW veterans, which suggested that free forms of PUFA were indeed affected in GWI [18], a large study evaluating free fatty acids and their bioactive lipid metabolites will help better elucidate the role of these fatty acids in ongoing inflammatory processes associated with GWI.


This study identifies distinct sex-specific plasma lipid profiles associated with GWI diagnosis, with significantly elevated TG and DG and lowered CE as a distinctive pattern in male GW veterans with GWI. Triglyceride AA-to-DHA ratios were also increased in male GW veterans with GWI. Since these lipid changes are associated with metabolic syndrome and inflammation, our findings suggest veterans with GWI have maladaptive biological processes that involve lipid transport and metabolism. While additional mechanistic work is needed to better elucidate lipid dysfunction in GWI, these blood lipid signatures suggest that focused studies are needed into the role of dyslipidemia in the ongoing pathophysiology of GWI. The current findings suggest that blood lipid monitoring may be beneficial for veterans with GWI. As such, further research is needed to characterize how these lipids and their associated pathways relate to GWI pathology to identify potential biomarkers for tracking general health related outcomes in GWI.

Data availability

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.



Gulf War Illness


Gulf War






Cholesterol esters




Saturated Fatty Acids


Polyunsaturated fatty acids


Docosahexaenoic acid


Arachidonic acid


Myalgic encephalomyelitis/chronic fatigue syndrome


Roskamp Neurology Clinic


Gulf War Illness Consortium


Center for Disease Control


Liquid chromatography–mass spectrometry


Electrospray ionization


Mass spectrometry












Quality control


Coefficients of Variance


Monounsaturated fatty acids


Analysis of variance


Mixed linear modeling




Standard deviation


Very low-density lipoproteins


Low-density lipoproteins


Genetics of Lipid-Lowering Drugs and Diet Network


Body mass Index


  1. White RF, et al. Recent research on Gulf War illness and other health problems in veterans of the 1991 Gulf War: effects of toxicant exposures during deployment. Cortex. 2016;74:449–75.

    Article  PubMed  Google Scholar 

  2. Sullivan K, et al. Neuropsychological functioning in military pesticide applicators from the Gulf War: Effects on information processing speed, attention and visual memory. Neurotoxicol Teratol. 2018;65:1–13.

    Article  CAS  PubMed  Google Scholar 

  3. Yee MK, et al. Longitudinal assessment of health symptoms in relation to neurotoxicant exposures in 1991 gulf war veterans the ft. Devens Cohort. J Occup Environ Med. 2020;62(9):663–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Zundel CG, et al. Changes in Health Status in the Ft Devens Gulf War Veterans Cohort: 1997–2017. Neurosci Insights. 2020.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Whistler T, et al. Impaired immune function in Gulf War Illness. BMC Med Genomics. 2009;2:1–13.

    Article  CAS  Google Scholar 

  6. Alshelh Z, et al. In-vivo imaging of neuroinflammation in veterans with Gulf War illness. Brain Behav Immun. 2020;87:498–507.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Ross MK, Borazjani A, Wang R, AllenCrow J, Xie S. Examination of the carboxylesterase phenotype in human liver. Arch Biochem Biophys. 2012;522(1):44–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Winrow CJ, Hemming ML, Allen DM, Quistad GB, Casida JE, Barlow C. Loss of neuropathy target esterase in mice links organophosphate exposure to hyperactivity. Nat Genet. 2003;33(4):477–85.

    Article  CAS  PubMed  Google Scholar 

  9. Michalovicz LT, Kelly KA, Sullivan K, O’Callaghan JP. Acetylcholinesterase inhibitor exposures as an initiating factor in the development of Gulf War Illness, a chronic neuroimmune disorder in deployed veterans. Neuropharmacology. 2020;171: 108073.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Feriani A, et al. Permethrin induced arterial retention of native and oxidized LDL in rats by promoting inflammation, oxidative stress and affecting LDL receptors, and collagen genes. Ecotoxicol Environ Saf. 2021;207: 111269.

    Article  CAS  PubMed  Google Scholar 

  11. Zhu L, Wang P, Sun YJ, Xu MY, Wu YJ. Disturbed phospholipid homeostasis in endoplasmic reticulum initiates tri-o-cresyl phosphate-induced delayed neurotoxicity. Sci Rep. 2016;6:1–9.

    Article  CAS  Google Scholar 

  12. Kotas ME, Medzhitov R. Homeostasis, Inflammation, and Disease Susceptibility. Cell. 2015;160(5):816–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Buckley ML, Ramji DP. The influence of dysfunctional signaling and lipid homeostasis in mediating the inflammatory responses during atherosclerosis. Biochim Biophys Acta - Mol Basis Dis. 2015;1852(7):1498–510.

    Article  CAS  Google Scholar 

  14. Abdullah L, et al. Chronic elevation of phosphocholine containing lipids in mice exposed to Gulf War agents pyridostigmine bromide and permethrin. Neurotoxicol Teratol. 2013;40:74–84.

    Article  CAS  PubMed  Google Scholar 

  15. Abdullah L, et al. Lipidomic profiling of phosphocholine containing brain lipids in mice with sensorimotor deficits and anxiety-like features after exposure to gulf war agents. NeuroMolecular Med. 2012;14(4):349–61.

    Article  CAS  PubMed  Google Scholar 

  16. Abdullah L, et al. Translational potential of long-term decreases in mitochondrial lipids in a mouse model of Gulf War Illness. Toxicology. 2016;372:22–33.

    Article  CAS  PubMed  Google Scholar 

  17. Emmerich T, et al. Phospholipid profiling of plasma from GW veterans and rodent models to identify potential biomarkers of Gulf War Illness. PLoS ONE. 2017;12(4):1–24.

    Article  CAS  Google Scholar 

  18. Joshi U, et al. Author Correction: Oleoylethanolamide treatment reduces neurobehavioral deficits and brain pathology in a mouse model of Gulf War Illness (Scientific Reports. Sci Rep. 2019;9:1.

    Article  CAS  Google Scholar 

  19. Krahmer N, Farese RV, Walther TC. Balancing the fat: Lipid droplets and human disease. EMBO Mol Med. 2013;5(7):973–83.

    Article  CAS  PubMed  Google Scholar 

  20. Jackson E, Shoemaker R, Larian N, Cassis L. Adipose tissue as a site of toxin accumulation. Compr Physiol. 2017;7(4):1085–135.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Briggs MA, Petersen KS, Kris-Etherton PM. Saturated fatty acids and cardiovascular disease: Replacements for saturated fat to reduce cardiovascular risk. Healthc. 2017.

    Article  Google Scholar 

  22. Calder PC. Omega-3 polyunsaturated fatty acids and inflammatory processes: Nutrition or pharmacology? Br J Clin Pharmacol. 2013;75(3):645–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Marion-Letellier R, Savoye G, Ghosh S. Polyunsaturated fatty acids and inflammation. IUBMB Life. 2015;67(9):659–67.

    Article  CAS  PubMed  Google Scholar 

  24. Regulska M, Szuster-Głuszczak M, Trojan E, Leśkiewicz M, Basta-Kaim A. The Emerging Role of the Double-Edged Impact of Arachidonic Acid- Derived Eicosanoids in the Neuroinflammatory Background of Depression. Curr Neuropharmacol. 2020;19(2):278–93.

    Article  Google Scholar 

  25. Kumar NG, et al. Dietary bioactive fatty acids as modulators of immune function: Implications on human health. Nutrients. 2019;11(12):1–15.

    Article  Google Scholar 

  26. Slade E, et al. Age and sex are associated with the plasma lipidome: findings from the GOLDN study. Lipids Health Dis. 2021;20(1):1–12.

    Article  CAS  Google Scholar 

  27. Unger ER, et al. CDC Grand Rounds: Chronic Fatigue Syndrome — Advancing Research and Clinical Education. MMWR Morb Mortal Wkly Rep. 2016;65(5051):1434–8.

    Article  PubMed  Google Scholar 

  28. Nkiliza A, et al. Sex-specific plasma lipid profiles of ME/CFS patients and their association with pain, fatigue, and cognitive symptoms. J Transl Med. 2021;19(1):1–15.

    Article  CAS  Google Scholar 

  29. Heboyan V, et al. Sex Differences in Gulf War Illness: A Reanalysis of Data from the CDC Air Force Study Using CDC and Modified Kansas Case Definitions. J Occup Environ Med. 2019;61(7):610–6.

    Article  PubMed  Google Scholar 

  30. Joshi U, et al. A permethrin metabolite is associated with adaptive immune responses in Gulf War Illness. Brain Behav Immun. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Steele L, et al. Brain–immune interactions as the basis of gulf war illness: Clinical assessment and deployment profile of 1990–1991 gulf war veterans in the gulf war illness consortium (gwic) multisite case-control study. Brain Sci. 2021;11:9.

    Article  Google Scholar 

  32. Steele L. Prevalence and Patterns of Gulf War Illness Steele Prevalence and Patterns of Gulf War Illness in Kansas Veterans: Association of Symptoms with Characteristics of Person, Place, and Time of Military Service. 2000.

  33. Fukuda K, et al. Chronic multisymptom illness affecting Air Force veterans of the Gulf War. J Am Med Assoc. 1998;280(11):981–8.

    Article  CAS  Google Scholar 

  34. NIH Sex Gender webpage.

  35. WHO. Gender and health in disasters. Who. 2015.

  36. Huguenard CJC, et al. Plasma Lipidomic Analyses in Cohorts With mTBI and/or PTSD Reveal Lipids Differentially Associated With Diagnosis and APOE ε4 Carrier Status. Front Physiol. 2020.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Folch J, Lees M, Stanley GHS. A simple method for the isolation and purification of total lipides from animal tissues. J Biol Chem. 1957;226(1):497–509.

    Article  CAS  PubMed  Google Scholar 

  38. Ernest B, Gooding JR, Campagna SR, Saxton AM, Voy BH. MetabR: An R script for linear model analysis of quantitative metabolomic data. BMC Res Notes. 2012;5(1):1.

    Article  Google Scholar 

  39. Abdullah L, et al. APOE ε4 specific imbalance of arachidonic acid and docosahexaenoic acid in serum phospholipids identifies individuals with preclinical Mild Cognitive Impairment/Alzheimer’s Disease. Aging. 2017;9(3):964–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Pang Z, et al. MetaboAnalyst 5.0: Narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 2021;49(W1):W388–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Abdullah L, et al. Proteomic CNS Profile of Delayed Cognitive Impairment in Mice Exposed to Gulf War Agents. Gulf. 2011.

    Article  Google Scholar 

  42. Ginsberg HN, Zhang YL, Hernandez-Ono A. Regulation of plasma triglycerides in insulin resistance and diabetes. Arch Med Res. 2005;36(3):232–40.

    Article  CAS  PubMed  Google Scholar 

  43. Alves-Bezerra M, Cohen DE. Triglyceride metabolism in the liver. Compr Physiol. 2018;8(1):1–22.

    Article  Google Scholar 

  44. Talayero BG, Sacks FM. The role of triglycerides in atherosclerosis. Curr Cardiol Rep. 2011;13(6):544–52.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Shamai L, et al. Association of body mass index and lipid profiles: Evaluation of a broad spectrum of body mass index patients including the morbidly obese. Obes Surg. 2011;21(1):42–7.

    Article  PubMed  Google Scholar 

  46. Cirulli ET, et al. Profound perturbation of the metabolome in obesity is associated with health risk. Cell Metab. 2019;29(2):488-500.e2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Zundel CG, et al. Rates of chronic medical conditions in 1991 Gulf war veterans compared to the general population. Int J Environ Res Public Health. 2019;16(6):1–16.

    Article  Google Scholar 

  48. Flock MR, Kris-Etherton PM. Diverse physiological effects of long-chain saturated fatty acids: Implications for cardiovascular disease. Curr Opin Clin Nutr Metab Care. 2013;16(2):133–40.

    Article  CAS  PubMed  Google Scholar 

  49. Johnson GJ, Slater BCS, Leis LA, Rector TS, Bach RR. Blood biomarkers of chronic inflammation in Gulf War illness. PLoS ONE. 2016;11(6):1–14.

    Article  CAS  Google Scholar 

  50. Parkitny L, Middleton S, Baker K, Younger J. Evidence for abnormal cytokine expression in Gulf War Illness : A preliminary analysis of daily immune monitoring data. BMC Immunol. 2015.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Liu JJ, Green P, JohnMann J, Rapoport SI, Sublette ME. Pathways of polyunsaturated fatty acid utilization: Implications for brain function in neuropsychiatric health and disease. Brain Res. 2015;1597(646):220–46.

    Article  CAS  PubMed  Google Scholar 

  52. López-Vicario C, et al. Pro-resolving mediators produced from EPA and DHA: Overview of the pathways involved and their mechanisms in metabolic syndrome and related liver diseases. Eur J Pharmacol. 2016;785:133–43.

    Article  CAS  PubMed  Google Scholar 

  53. Nishizaki Y, et al. Association between the docosahexaenoic acid to arachidonic acid ratio and acute coronary syndrome: A multicenter observational study. BMC Cardiovasc Disord. 2016;16(1):1–10.

    Article  CAS  Google Scholar 

  54. Garza-Lombó C, et al. Circulating HMGB1 is elevated in veterans with Gulf War Illness and triggers the persistent pro-inflammatory microglia phenotype in male C57Bl/6J mice. Transl Psychiatry. 2021;11(1):1–12.

    Article  CAS  Google Scholar 

  55. Boon J, et al. Ceramides contained in LDL are elevated in type 2 diabetes and promote inflammation and skeletal muscle insulin resistance. Diabetes. 2013;62(2):401–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. De Mello VDF, et al. Link between plasma ceramides, inflammation and insulin resistance: Association with serum IL-6 concentration in patients with coronary heart disease. Diabetologia. 2009;52(12):2612–5.

    Article  CAS  PubMed  Google Scholar 

  57. Siskind LJ, et al. The BCL-2 protein BAK is required for long-chain ceramide generation during apoptosis. J Biol Chem. 2010;285(16):11818–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Green CD, Maceyka M, Cowart LA, Spiegel S. Review Sphingolipids in metabolic disease : The good, the bad, and the unknown. Cell Metab. 2021;33(7):1293–306.

    Article  CAS  PubMed  Google Scholar 

  59. Szamatowicz J, Juchnicka I, Kuz M. “Ceramides and Sphingosino-1. Phosphate Obesit. 2021;12:1–7.

    Article  Google Scholar 

Download references


We would like to thank David Patterson for his help in recruitment. We would like to thank Daniel Niedospial for his help in proof reading. We are grateful to the veterans and study subjects who agreed to participate in this study.


This material is based upon work supported by a Congressionally Directed Medical Research Program (CDMRP) award (GW150056) and a VA Merit award (5I01BX002084-04) from the US Department of Veterans Affairs to Dr. Fiona Crawford and by the Roskamp Foundation. Plasma samples from GW veterans were provided by a CDMRP GWI consortium award (GW120037) to Dr. Kimberly Sullivan and a VA merit and a CDMRP award to Dr. Nancy Klimas.

Author information

Authors and Affiliations



SO: Conceptualization, Investigation, Formal Analysis, Writing—Original Draft, Writing—Review and Editing, Visualization, Data Curation AN: Formal Analysis, Writing—Original Draft, Writing—Review and Editing, Visualization. MP: Formal Analysis, Writing—Original Draft, Writing—Review and Editing, Visualization. JEE: Methodology, Formal Analysis, Data Curation, Writing—Review and Editing, NK: Resources, Funding acquisition, Writing—Review and Editing, APK: Resources, Investigation, Writing—Review and Editing, KS: Resources, Funding acquisition, Writing—Review and Editing, MHK: Resources, Writing—Review and Editing, MM: Conceptualization, Writing—Review and Editing, FC: Conceptualization, Methodology, Funding acquisition, Investigation, Resources, Supervision, Project administration, Writing—Review and Editing, LA: Conceptualization, Methodology, Resources, Writing—Original Draft, Writing—Review and Editing, Formal Analysis, Supervision, Funding acquisition, Investigation, Project administration. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Laila Abdullah.

Ethics declarations

Ethics approval and consent to participate

An informed consent was collected for each participant and was reviewed by the IRB of each respective site.

Consent for publication

Not applicable.

Competing interests

Not applicable.


The contents of this poster do not represent the views of the Department of Veterans Affairs or the United States Government.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Individual levels differences in major blood classes between GWI and controls. Bar graphs indicate mean ± SE and individual values of total lipids (A), degree of unsaturation of each lipid class (B) and AA and DHA content of these lipids (C) that were significantly different between controls and GWI. Black asterisks (*) indicate significant differences between GWI cases and their sex-matched controls (p < 0.05).

Additional file 2: Figure S2.

Association of blood lipid profiles and the presence of cardiovascular risk factors among the study population. Major lipids were further evaluated to determine the impact of presence of reported cardiovascular risk factors (no = 8 and yes = 32). Heatmap generated from z-scores transformed from individual concentrations. Differences in AA and DHA containing PL and neutral lipids were affected the presence of cardiovascular risk factors. Total lipid content was unaffected by cardiovascular risk factors. Black asterisks (*) indicate significant differences between GWI cases and their sex-matched controls (p < 0.05).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Oberlin, S., Nkiliza, A., Parks, M. et al. Sex-specific differences in plasma lipid profiles are associated with Gulf War Illness. J Transl Med 20, 73 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: