Skip to main content
Download PDF
Download ePub

Age and sex are associated with the plasma lipidome: findings from the GOLDN study

Lipids in Health and Disease volume 20, Article number: 30 (2021) Cite this article

Abstract

Background

Developing an understanding of the biochemistry of aging in both sexes is critical for managing disease throughout the lifespan. Lipidomic associations with age and sex have been reported, but prior studies are limited by measurements in serum rather than plasma or by participants taking lipid-lowering medications.

Methods

Our study included lipidomic data from 980 participants aged 18–87 years old from the Genetics of Lipid-Lowering Drugs and Diet Network (GOLDN). Participants were off lipid-lowering medications for at least 4 weeks, and signal intensities of 413 known lipid species were measured in plasma. We examined linear age and sex associations with signal intensity of (a) 413 lipid species; (b) 6 lipid classes (glycerolipids, glycerophospholipids, sphingolipids, sterol lipids, fatty acids, and acylcarnitines); and (c) 15 lipid subclasses; as well as with the particle sizes of three lipoproteins.

Results

Significant age associations were identified in 4 classes, 11 subclasses, 147 species, and particle size of one lipoprotein while significant sex differences were identified in 5 classes, 12 subclasses, 248 species, and particle sizes of two lipoproteins. For many lipid species (n = 97), age-related associations were significantly different between males and females. Age*sex interaction effects were most prevalent among phosphatidylcholines, sphingomyelins, and triglycerides.

Conclusion

We identified several lipid species, subclasses, and classes that differ by age and sex; these lipid phenotypes may serve as useful biomarkers for lipid changes and associated cardiovascular risk with aging in the future. Future studies of age-related changes throughout the adult lifespan of both sexes are warranted.

Trial registration

ClinicalTrials.gov NCT00083369; May 21, 2004.

Introduction

Global life expectancy has increased [1] and is predicted to continue increasing in many countries [2]. In developed nations, age is the primary risk factor for many common diseases including cardiovascular disease, cancer, and Alzheimer’s disease [3]. Human aging is a complex phenomenon influenced by environmental, genetic and other endogenous factors. For example, the aging process differs between sexes, and there are well established associations between sex and the incidence, symptoms, severity, age of onset, and response to treatments for many diseases. Management of disease in aging populations may be improved by better understanding the physiology of aging in both sexes. Given the complex biochemistry of aging, metabolomic approaches may offer important insights into the aging process and age-related disease etiologies. Lipidomics, in particular, is well positioned to investigate those pathologies. Lipids function as both structural and signaling molecules and are known to play important roles in age-related disorders such as cardiovascular disease [4], metabolic syndrome [5], macular degeneration [6], neurodegenerative diseases [7], and at least some types of stroke [8]. However, before lipidomics can be used to dissect and characterize disease, it is critical that we understand how the lipidome changes during normal, healthy aging in both sexes.

A number of well-powered metabolomic studies of age and sex have included lipidomic species [9,10,11,12]; however, most of these studies measured metabolites in serum, while current recommendations favor plasma for lipidomics due to the generation and degeneration of lipid species during the coagulation process and the similarity of plasma generated from fresh blood to the plasma compartment in vivo [13]. Beyene et al. studied a large Australian cohort (n = 10,339) to examine age and sex associations with individual plasma lipid species that are independent of lipoprotein metabolism [12]. Wong et al. conducted a plasma lipidomics study in a modestly sized sample (n = 100) of participants age > 55 years [14]. Darst et al. conducted a longitudinal plasma lipidomics analysis of 1212 participants with up to three study visits per participant, with follow-up occurring four years after baseline and every two years thereafter [15, 16]. However, these studies were not restricted to participants free of lipid-lowering medication, somewhat complicating the interpretation of their findings.

The Genetics of Lipid-Lowering Drugs and Diet Network (GOLDN) offers a context to investigate lipidomic markers of age and sex. Data on GOLDN’s 980 participants includes robust characterization of health status, diet, anthropometry, and demographics. No GOLDN participants were taking lipid-lowering drugs, including fish oils, at the time specimens were collected, thus avoiding the difficulties of interpretation that have characterized some previous studies. GOLDN participants ranged in age from 18 to 87 years, making it better suited for a study across the adult lifespan than previous studies with a more age-restricted sample. In this study we used liquid chromatography/mass spectrometry to quantify 413 lipidomic species in GOLDN plasma specimens. The goal of this study is to examine associations between age and each lipid species and to assess whether these associations across the adult lifespan differ between sexes.

Methods

The GOLDN Study has been described in detail elsewhere [17]. GOLDN participants were re-recruited from three-generational pedigrees from two National Heart, Lung and Blood Institute Family Heart Study field centers (Minneapolis, MN and Salt Lake City, UT); genetically related pedigree members of these probands were recruited as well. Nearly all participants were of European ancestry. Earlier studies demonstrated that Caucasians in Utah and Minnesota were homogeneous, and pooling data across centers did not threaten the validity of this study [18].

Exclusion criteria included the following: recent history (past six months) of myocardial infarction, coronary bypass surgery, coronary angioplasty or percutaneous transluminal coronary angioplasty; current use of warfarin or insulin; self-report of a positive history of kidney, gall bladder, liver, or pancreatic disease, or a history of nutrient malabsorption; serum concentrations of alanine transaminase exceeding 66 U/L in males or 44 U/L in females; serum concentrations of aspartate aminotransferase exceeding 52 U/L in males or 42 U/L in females; glomerular filtration rate < 30 ml/min/1.73 m2 estimated from the Modification of Diet in Renal Disease (MDRD) Study equation; fasting triglycerides ≥1500 mg/dL; pregnant women or women of childbearing potential not using contraception; and women nursing a child. Individuals who reported current use of prescription and/or over-the-counter hypolipidemic drugs or dietary supplements known to influence lipid values (e.g., fish oil, niacin, flaxseed oil) were required to consult with their clinician for approval to discontinue these lipid-lowering agents for four weeks prior to study participation.

Specimens were collected at a pre-intervention, baseline visit. Participants were asked to fast for ≥12 h and abstain from using alcohol for ≥24 h before the clinic visit [19]. Blood samples were centrifuged within 20 min of collection at 2000 g for 15 min at 4 °C [20]. Plasma samples were stored at − 80 °C.

Unique, known lipid species were analyzed by liquid chromatography coupled-mass spectrometry (LC-MS/MS). The data set was obtained as an untargeted lipidomic profiling experiment. Differences in signal intensity reflect differences in abundance of any given lipid species based on instrumental settings and solvents and buffers used during analysis. The relative concentrations are comparable within the experiments obtained with the same protocol but, by design, cannot be easily compared with absolute (targeted) measurements. The signals of lipids designated as (A) and (B) are related to double bond stereoisomers (cis/trans) which can be resolved by the chromatographic protocol used. Due to missing authentic reference compounds, no further annotations can be made. Moreover, tetrahedral stereoisomers (R/S) and regioisomers were not distinguished. Raw lipidomic data were normalized using Systematic Error Removal using Random Forest (SERRF) [21]. SERRF is a machine learning algorithm for large metabolomic data sets that characterizes systematic error such as batch effects and day-to-day variation in quality control samples and uses this model to reduce systematic error in study samples [21]. For lipid species measured on both positive and negative modes, the mode with the smaller standard deviation across participants was retained for analysis.

For statistical analysis, individual lipid species were classified into six major classes: glycerolipids, glycerophospholipids, sphingolipids, sterol lipids, fatty acids, and acylcarnitines. Four of the six major classes were further subdivided into 15 subclasses: glycerolipids were subdivided into triglycerides (TG) and diacylglycerols (DG); glycerophospholipids were subdivided into phosphatidylcholines (PC), phospatidylethanolamines (PE), phosphatidylinositols (PI), phosphatidylglycerols (PG), lysophosphatidylethanolamines (LPE), and lysophosphatidylcholines (LPC); sphingolipids were subdivided into sphingomyelins (SM), ceramides, lactosylceramides (LCer), glucosylceramides (GlcCer), and galactosylgalactosylceramides (GalGalCer); sterol lipids were subdivided into cholesterol and cholesteryl esters (CE). Fatty acids and acylcarnitines were not further subdivided. Analyses were performed at the class level, subclass level, and individual lipid level. For class and subclass analyses, each participant’s observed class/subclass intensity was calculated by summing the total signal intensity of lipid species belonging to the class/subclass for that participant.

Three sets of analyses were performed to examine the associations of age, sex, and their interaction with signal intensities of (1) classes of lipids, (2) subclasses of lipids, and (3) individual lipid species. The class intensities, subclass intensities, and individual lipid intensities were standardized to have mean 0 and standard deviation 1 prior to analysis. This standardization of outcome data allows for a direct comparison of estimated regression coefficients between models. Further, it aids in interpretation as a one unit change in outcome reflects a change of one standard deviation in signal intensity for all lipid species, classes, and subclasses. No further transformation was performed to the classes, subclasses, and individual lipid species.

A fourth set of analyses was also performed to examine the associations of age, sex, and their interaction with lipoprotein particle size measured by nuclear magnetic resonance (NMR) spectroscopy. Lipoproteins measured via NMR spectroscopy include very-low-density lipoprotein (VLDL), low-density lipoprotein (LDL), and high-density lipoprotein (HDL). All particle sizes are measured in nanometers (nm).

All analyses employed linear regression modeling to examine the associations between the class, subclass, lipid species, or lipoprotein particle size outcome (each separately) with age, sex, and age*sex interaction, while including adjustment for body mass index (BMI). Other studies have shown BMI associations with the plasma lipidome [22, 23]. As such, BMI is included as a covariate in all analyses in order to examine age and sex associations with the plasma lipidome that are not driven by BMI-related associations. Analyses of classes, subclasses, and individual lipid species also included adjustment for batch effects. For analyses including an age*sex interaction, the age variable was centered prior to analysis to aid in interpretability of the interaction term. For classes, subclasses, individual lipid species, and lipoprotein particle sizes exhibiting significant age*sex interaction effects, age associations with the class, subclass, lipid species, or lipoprotein particle size outcome were further explored in sex-stratified analyses, i.e., separate linear regression models were fit for males and females to explore the relationship between age and the class, subclass, lipid species, or lipoprotein particle size in a more easily interpretable form.

To account for multiple testing, a Benjamini-Hochberg adjustment (false discovery rate correction) was applied to control the false discovery rate in each set of analyses at 0.05 [24]. All reported P-values are the Benjamini-Hochberg-adjusted P-values, and thus, a significance level of 0.05 was utilized for all analyses. All analyses were performed using R version 4.0.2 [25].

Results

Participant demographics

Fasting plasma samples were analyzed for 980 participants with 413 unique, known lipid species identified. The mean age of participants was 48.3 years (standard deviation = 16.4 years), and 52.3% were female (Table 1). As expected by design of the GOLDN study, participants identified primarily as white (99.9%) (Table 1).

Table 1 GOLDN cohort characteristics, mean (standard deviation) or n (%)
Full size table

Age and sex are associated with lipid classes

Among the six classes of lipid species, four had significant associations between age and total class intensity after adjustment for sex, BMI, and batch (Table 2). These include glycerophospholipids, sphingolipids, sterol lipids, and fatty acids, which all increase significantly with age. On average, total intensity of glycerolipids and acylcarnitines did not change significantly with age after adjustment for sex, BMI, and batch (Table 2). Five of the six classes of lipid species had a significant association with sex (Table 2). Total intensities of glycerophospholipids, sphingolipids, sterol lipids, and fatty acids were significantly higher in women than men, and total intensity of glycerolipids was significantly higher in men than women after adjustment for age, BMI, and batch. The average total intensity of acylcarnitines did not differ significantly between males and females after adjustment for age, BMI, and batch (Table 2).

Table 2 Associations between total class intensity with age, sex, and their interaction
Full size table

In two classes of lipids, glycerophospholipids (P = 0.007) and glycerolipids (P = 0.022), age and sex had a significant interaction effect, indicating that the association between age and the total class intensity was different among males and females, while controlling for BMI and batch (Table 2). In both classes, the expected increase in the total class intensity with increasing age was larger in women than in men (Fig. 1). For a ten-year increase in age, our model predicts the total intensity of glycerophospholipids to increase by 0.05 standard deviations in men and by 0.15 standard deviations in women after controlling for BMI and batch (Table 2, Fig. 1). For a ten-year increase in age, we expect the total intensity of glycerolipids to increase by 0.13 standard deviations in men and by 0.25 standard deviations in women after controlling for BMI and batch (Table 2, Fig. 1). Correlations between total intensities of the six lipid classes are shown in Additional File 2.

Fig. 1
figure 1

Relationship between age and standardized total class intensity for glycerophospholipids (left) and glycerolipids (right). Dots represent observed values for males (blue) and females (red). Lines represent a linear regression of class intensity on age, performed separately for males (blue) and females (red). Only lipid classes with significant age*sex interaction effects from the overall analyses (see Table 2) are analyzed in these stratified models

Full size image

Age and sex are associated with lipid subclasses

Among the 15 subclasses of lipid species, 11 had significant associations with age after adjustment for sex, BMI, and batch (Table 3). Of these, cholesterol was the only subclass that decreased significantly with age. The other 10 (PC, PE, PI, PG, LPE, SM, ceramide, LCer, GlcCer, CE) all increased with age. Of the 15 subclasses of lipid species, 12 had significant associations with sex after adjustment for age, BMI, and batch (Table 3). On average, women had higher total subclass intensity of PC, PE, PI, PG, SM, GlcCer, and CE, while men had higher total subclass intensity of TG, DG, LPC, ceramides, and cholesterol after adjustment for age, BMI, and batch (Table 3).

Table 3 Associations between total subclass intensity with age, sex, and their interaction
Full size table

Five subclasses of lipid species exhibited a significant interaction effect between age and sex (Table 3). These included GalGalCer (P = 0.006), PE (P = 0.008), PC (P = 0.012), TG (P = 0.028), and LPC (P = 0.033). In GalGalCer and LPC, the total intensity remained approximately constant across ages in men, but in women there was a significant decrease in GalGalCer with age and a significant increase in LPC with age (Table 3, Fig. 2). In the PE, PC, and TG subclasses, total subclass intensity increased with age in both men and women, but the rate of increase was larger in women than men after controlling for BMI and batch (Table 3, Fig. 2).

Fig. 2
figure 2

Relationship between age and standardized total subclass intensity. Dots represent observed values for males (blue) and females (red). Lines represent a linear regression of subclass intensity on age, performed separately for males (blue) and females (red). Only lipid subclasses with significant age*sex interaction effects from the overall analyses (see Table 3) are analyzed in these stratified models. These include galactosylgalactosylceramides (top left), phosphatidylethanolamines (top center), phosphatidylcholines (top right), triglycerides (bottom left), and lysophosphatidylcholines (bottom center)

Full size image

Age and sex are associated with individual lipid species

Among the 413 individual lipid species, 147 had significant associations with age after adjustment for sex, BMI, and batch (Table 4, Additional File 1); on average, 141 increased with age, and 6 decreased with age. Among the 413 individual lipid species, 248 had a significant association with sex after adjustment for age, BMI, and batch (Table 4, Additional File 1). Of these, 103 lipids had higher levels in males and 145 had higher levels in females, on average (Table 4, Additional File 1).

Table 4 Associations between lipid signal intensity with age, sex, and their interaction (includes those with top 10 most significant age*sex interactions)
Full size table

For 97 individual lipid species, the association between age and lipid intensity was significantly different between males and females after adjustment for batch and BMI (Table 4, Additional File 1). All but three of these significant age*sex interaction effects were in the form of larger age-related increases in lipid intensity in females as compared to males (Fig. 3). Only PC (p-36:1) or PC (o-36:2), GalGalCer/LCer (d18:1/16:0), and PC (p-42:2) or PC (o-42:3) deviated from this pattern; the first two exhibited age-related declines in females with very little change in males, and the last exhibited an age-related decline in females and an age-related increase in males.

Fig. 3
figure 3

Relationship between age and standardized individual lipid intensity, stratified by sex. Bars indicate the age coefficient, representing the expected change in standardized lipid intensity for a one-year increase in age, from separate linear regression models for males (blue) and females (red). Only lipid species with significant age*sex interaction effects from the overall analyses (see Table 4) are analyzed in these stratified models. Lipid species are grouped by class for ease of presentation (panel A: glycerolipids, panel B: glycerophospholipids, panel C: sphingolipids, panel D: sterol lipids); within each class, lipids are ordered by size of the age*sex interaction effect (see Table 4). Lipids denoted A and B designate resolved cis/trans stereoisomers

Full size image

Analysis of lipoprotein particle sizes

Of the three lipoproteins measured via NMR spectroscopy, only VLDL exhibited a significant association with age after adjustment for sex and BMI (Table 5). This association between age and VLDL particle size was also significantly different between males and females after adjustment for BMI (Table 5, Fig. 4). In males, VLDL particle size decreased with age by an average of 0.090 nm per year; and in females, VLDL particle size slightly increased with age by an average of 0.036 nm per year, after adjusting for BMI (Table 5).

Table 5 Associations between lipoprotein particle size with age, sex, and their interaction
Full size table
Fig. 4
figure 4

Relationship between age and lipoprotein particle size, stratified by sex. Dots represent observed values for males (blue) and females (red). Lines represent a linear regression of lipoprotein particle size on age, performed separately for males (blue) and females (red). Only particles with significant age*sex interaction effects from the overall analyses (see Table 5) are analyzed in these stratified models. This includes only VLDL

Full size image

There was not a significant association between LDL and HDL particle size with age, but both lipoprotein particle sizes were significantly different between males and females after adjustment for BMI (Table 5). On average, LDL particles were 0.587 nm larger in females than in males, and HDL particles were 0.359 nm larger in females than in males, after adjustment for BMI (Table 5).

Discussion

Findings from the GOLDN study indicated many age and sex-related associations in the plasma lipidome in a large (n = 980) cohort of participants. GOLDN reflects a unique cohort for lipidomic analyses because all participants were treatment-naïve at the time of plasma collection and therefore, results are not confounded by the effects of lipid-lowering drugs. Many significant differences in the lipidomic measures were identified in the lipid classes, lipid subclasses, and individual lipid species for age and sex. These findings are highly consistent with the results from the Australian diabetes (AusDiab) cohort study which evaluated 706 lipids species in plasma; of the 11 subclasses reported in GOLDN to differ by sex, 8 were consistent in the direction of the sex difference and significance in the AusDiab cohort [12]. Similarly, of the 11 subclasses reported in GOLDN to differ by age, 9 were significantly associated with age in the AusDiab cohort [12]. Most interesting, however, is that the relationship between age and the lipidomic measures in GOLDN differed by sex; in one-third (two out of six) of lipid classes, one-third (five out of 15) of lipid subclasses, and about one fourth (97 out of 413) of individual lipid species, age demonstrated either a steeper increase or decline of the lipidomic measure between men and women. Four of the five subclasses reported in GOLDN to exhibit differential age-related associations by sex were also identified in the AusDiab cohort [12].

It is well understood that circulating lipoproteins are affected by female sex steroids, in particular estrogens which demonstrate more robust age-related fluctuations associated with menopause. In the GOLDN study, females had greater age-related increases in both glycerolipids and glycerophospholipids (Table 2, Fig. 1). Among glycerolipids, this sex difference in the relationship of age was also observed in triglycerides but not diglycerides (Table 3). These species are almost exclusively associated with circulating lipoproteins. The present study evaluated the plasma lipidome; further analysis of the lipidome within lipoproteins (LDL, HDL) and their subclasses (HDL1, 2, etc.) may reveal additional insight into the origins of plasma lipids, their potential as biomarkers, and perhaps the impact on aging and disease risk within biological sex.

Galactosylgalactosylceramides (GalGalCer) had an age-related decline in females that was not observed in males (Table 2, Fig. 1). While the plasma source of GalGalCer cannot be discerned from our lipidomic measures, α-galactosidase, a key metabolic enzyme in ceramides, is located on the X chromosome [26]. In addition, messenger RNA for the enzyme has been shown to be significantly repressed in MCF7 breast cancer cells in which estrogen receptor-α had been silenced by siRNAs [27,28,29]. Consistent with this observation, treatment of breast cancer cell lines with 17β-estradiol has been shown to produce a time-dependent increase in mRNA expression [30, 31]. Collectively, these data suggest a positive correlation between estrogen signaling and α-galactosidase expression that may underlie the age-associated changes in enzyme substrates over the lifespan in females but not in males.

We also found a sex-dependent age-related increase in choline-containing phospholipids, including phosphatidylcholines (PCs) which are the most abundant phospholipids in cellular membranes. Like PCs, PEs also demonstrated an age-related increase in females but not males. PEs are significantly less abundant in biological membranes than their PC counterparts (~ 20–30% if membrane lipids in liver), but they play critical roles in membrane biology and cellular signaling. PEs are asymmetrically distributed in plasma membranes where 80% are confined to the inner leaflet. The loss of that asymmetry is an inflammatory signaling event in a number of tissues, including heart [32]. PEs are also highly enriched in mitochondrial membranes. Thus, changes in their abundance across the lifespan may reflect well-established age-related declines in basal metabolic rate [33]. The most significant PE demonstrating an age*sex interaction was 40:7, which contains an omega – 3 PUFA. Such lipids are well established to influence a number of cardiovascular functions including endothelial and cardiomyocyte function and hemostasis [34].

Study strengths and limitations

This study has several strengths. First, the large sample size (n = 980) includes a robust characterization of underlying health characteristics such as diet, anthropometry, and medical conditions. Additionally, this large cohort spans a wide range of ages (18 to 87 years), offering a comprehensive view across the adult lifespan. Furthermore, all participants were free of lipid-lowering medications. The measurement of participants’ lipid levels in plasma is also a strength of this study. Although serum is more widely used in clinical practice, the coagulation process used to isolate serum can lead to generation or degeneration of lipid species. Plasma generated from fresh blood is considered to be the most similar to the plasma compartment in vivo; therefore, current recommendations for lipidomic analysis favor plasma over serum [13].

Despite these strengths, some limitations do exist. Examination of age associations in this study is limited by the cross-sectional nature of data collection; as such, our results may be used to target specific lipid species or classes in future longitudinal studies examining intra-individual age-related changes. Further, the generalizability of results from this study is limited to similar populations, and these findings should be confirmed in other populations. In particular, participants in the GOLDN cohort predominantly identify as white, so results should only be generalized to populations of European descent.

Residual plots and histograms of residuals were visually examined to assess modeling assumptions of homoscedasticity and normality of residuals for class and subclass analyses. No strong departures from these assumptions were observed, but some minor departures do exist, noted in the glycerolipids class and its two subclasses, TG and DG, as shown in Additional File 3. Because this analysis is a large lipidomic-scale screening phase, we utilized straightforward linear modeling. The intensities of some lipid species, classes, and subclasses may plateau at a certain age, possibly around menopause in women. Our lipidomic screening phase analysis has identified two classes, five subclasses, and 97 individual lipid species that exhibit age- and sex-related associations. We recommend further exploration of these associations in this smaller set of lipids using more complex methodologies that can capture non-linear trends, such as spline-based or change point analyses.

Conclusion

The results of this study show evidence of age and sex-related associations in several lipid species, classes, subclasses, and lipoproteins. These findings underscore the need for sex stratification when examining age-related changes in the plasma lipidome. Future studies may seek to explore more complex age-related associations and/or associations with disease in a small set of lipid species in order to preserve statistical power; the lipid species or groupings identified in our study should be targeted for this purpose. Our findings add to the understanding of changes in the plasma lipidome across the adult lifespan in both sexes, a feature that can inform advances in precision medicine to manage disease in aging populations.

Availability of data and materials

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Abbreviations

CE:

Cholesteryl ester

DG:

Diacylglycerol

GalGalCer:

Galactosylgalactosylceramide

GlcCer:

Glucosylceramide

GOLDN:

Genetics of Lipid-Lowering Drugs and Diet Network (study)

HOMA-IR:

Homeostatic model assessment of insulin resistance

LC-MS/MS:

Liquid chromatography coupled-mass spectrometry

LCer:

Lactosylceramide

LPC:

Lysophosphatidylcholine

LPE:

Lysophosphatidylethanolamine

MDRD:

Modification of Diet in Renal Disease (study)

PC:

Phosphatidylcholine

PE:

Phospatidylethanolamine

PG:

Phosphatidylglycerol

PI:

Phosphatidylinositol

SE:

Standard error

SERRF:

Systematic error removal using random forest

TG:

Triglyceride

References

  1. Gulland A. Global life expectancy increases by five years. BMJ. 2016;353:i2883.

    Article  Google Scholar 

  2. Foreman KJ, Marquez N, Dolgert A, Fukutaki K, Fullman N, McGaughey M, et al. Forecasting life expectancy, years of life lost, and all-cause and cause-specific mortality for 250 causes of death: reference and alternative scenarios for 2016-40 for 195 countries and territories. Lancet. 2018;392(10159):2052–90. https://doi.org/10.1016/S0140-6736(18)31694-5.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Niccoli T, Partridge L. Ageing as a risk factor for disease. Curr Biol. 2012;22(17):R741–52. https://doi.org/10.1016/j.cub.2012.07.024.

    Article  PubMed  CAS  Google Scholar 

  4. Ference BA, Graham I, Tokgozoglu L, Catapano AL. Impact of lipids on cardiovascular health: Jacc health promotion series. J Am Coll Cardiol. 2018;72(10):1141–56. https://doi.org/10.1016/j.jacc.2018.06.046.

    Article  PubMed  CAS  Google Scholar 

  5. Roche HM. Fatty acids and the metabolic syndrome. Proc Nutr Soc. 2005;64(1):23–9. https://doi.org/10.1079/PNS2004405.

    Article  PubMed  CAS  Google Scholar 

  6. Shen J, He J, Wang F. Association of lipids with age-related macular degeneration. Discov Med. 2016;22(120):129–45.

    PubMed  Google Scholar 

  7. Wong MW, Braidy N, Poljak A, Pickford R, Thambisetty M, Sachdev PS. Dysregulation of lipids in Alzheimer's disease and their role as potential biomarkers. Alzheimers Dement. 2017;13(7):810–27. https://doi.org/10.1016/j.jalz.2017.01.008.

    Article  PubMed  Google Scholar 

  8. Holmes MV, Millwood IY, Kartsonaki C, Hill MR, Bennett DA, Boxall R, et al. Lipids, lipoproteins, and metabolites and risk of myocardial infarction and stroke. J Am Coll Cardiol. 2018;71(6):620–32. https://doi.org/10.1016/j.jacc.2017.12.006.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. Mittelstrass K, Ried JS, Yu Z, Krumsiek J, Gieger C, Prehn C, et al. Discovery of sexual dimorphisms in metabolic and genetic biomarkers. PLoS Genet. 2011;7(8):e1002215. https://doi.org/10.1371/journal.pgen.1002215.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Dunn WB, Lin W, Broadhurst D, Begley P, Brown M, Zelena E, et al. Molecular phenotyping of a UK population: defining the human serum metabolome. Metabolomics. 2015;11(1):9–26. https://doi.org/10.1007/s11306-014-0707-1.

    Article  PubMed  CAS  Google Scholar 

  11. Weir JM, Wong G, Barlow CK, Greeve MA, Kowalczyk A, Almasy L, et al. Plasma lipid profiling in a large population-based cohort. J Lipid Res. 2013;54(10):2898–908. https://doi.org/10.1194/jlr.P035808.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Beyene HB, Olshansky G, AATS, Giles C, Huynh K, Cinel M, et al. High-coverage plasma lipidomics reveals novel sex-specific lipidomic fingerprints of age and BMI evidence from two large population cohort studies. PLoS Biol. 2020;18(9):e3000870. https://doi.org/10.1371/journal.pbio.3000870.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Burla B, Arita M, Arita M, Bendt AK, Cazenave-Gassiot A, Dennis EA, et al. MS-based lipidomics of human blood plasma: a community-initiated position paper to develop accepted guidelines. J Lipid Res. 2018;59(10):2001–17. https://doi.org/10.1194/jlr.S087163.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Wong MWK, Braidy N, Pickford R, Vafaee F, Crawford J, Muenchhoff J, et al. Plasma lipidome variation during the second half of the human lifespan is associated with age and sex but minimally with BMI. PLoS One. 2019;14(3):e0214141. https://doi.org/10.1371/journal.pone.0214141.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Darst BF, Koscik RL, Hogan KJ, Johnson SC, Engelman CD. Longitudinal plasma metabolomics of aging and sex. Aging (Albany NY). 2019;11(4):1262–82. https://doi.org/10.18632/aging.101837.

    Article  CAS  Google Scholar 

  16. Johnson SC, Koscik RL, Jonaitis EM, Clark LR, Mueller KD, Berman SE, et al. The Wisconsin registry for Alzheimer's prevention: a review of findings and current directions. Alzheimers Dement (Amst). 2018;10(1):130–42. https://doi.org/10.1016/j.dadm.2017.11.007.

    Article  Google Scholar 

  17. Lai CQ, Arnett DK, Corella D, Straka RJ, Tsai MY, Peacock JM, et al. Fenofibrate effect on triglyceride and postprandial response of apolipoprotein A5 variants: the GOLDN study. Arterioscler Thromb Vasc Biol. 2007;27(6):1417–25. https://doi.org/10.1161/ATVBAHA.107.140103.

    Article  PubMed  CAS  Google Scholar 

  18. Pankow JS, Province MA, Hunt SC, Arnett DK. Regarding "testing for population subdivision and association in four case-control studies". Am J Hum Genet. 2002;71(6):1478–80. https://doi.org/10.1086/344582.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. Tsai MY, Georgopoulos A, Otvos JD, Ordovas JM, Hanson NQ, Peacock JM, et al. Comparison of ultracentrifugation and nuclear magnetic resonance spectroscopy in the quantification of triglyceride-rich lipoproteins after an oral fat load. Clin Chem. 2004;50(7):1201–4. https://doi.org/10.1373/clinchem.2004.032938.

    Article  PubMed  CAS  Google Scholar 

  20. Irvin MR, Zhi D, Aslibekyan S, Claas SA, Absher DM, Ordovas JM, et al. Genomics of post-prandial lipidomic phenotypes in the genetics of lipid lowering drugs and diet network (GOLDN) study. PLoS One. 2014;9(6):e99509. https://doi.org/10.1371/journal.pone.0099509.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  21. Fan S, Kind T, Cajka T, Hazen SL, Tang WHW, Kaddurah-Daouk R, et al. Systematic error removal using random forest for normalizing large-scale untargeted lipidomics data. Anal Chem. 2019;91(5):3590–6. https://doi.org/10.1021/acs.analchem.8b05592.

    Article  PubMed  CAS  Google Scholar 

  22. Bagheri M, Djazayery A, Farzadfar F, Qi L, Yekaninejad MS, Aslibekyan S, et al. Plasma metabolomic profiling of amino acids and polar lipids in Iranian obese adults. Lipids Health Dis. 2019;18(1):94. https://doi.org/10.1186/s12944-019-1037-0.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Cirulli ET, Guo L, Leon Swisher C, Shah N, Huang L, Napier LA, et al. Profound perturbation of the metabolome in obesity is associated with health risk. Cell Metab. 2019;29:488–500 e2.

    Article  CAS  Google Scholar 

  24. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc Ser B. 1995;57:289–300.

    Google Scholar 

  25. R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2017.

    Google Scholar 

  26. Rebourcet R, Weil D, Van Cong N, Frézal J. Alpha-galactosidase: a dimeric enzyme dependent on a structural locus on the X chromosome. Cytogenet Cell Genet. 1975;14(3-6):406–8. https://doi.org/10.1159/000130394.

    Article  PubMed  CAS  Google Scholar 

  27. Al Saleh S, Al Mulla F, Luqmani YA. Estrogen receptor silencing induces epithelial to mesenchymal transition in human breast cancer cells. PLoS One. 2011;6(6):e20610. https://doi.org/10.1371/journal.pone.0020610.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Al Saleh S, Al Mulla F, Luqmani YA. NAGA - estrogen receptor alpha-silenced MCF7 breast cancer cells (identifier 77629993). GEO Profiles. 2011. https://www.ncbi.nlm.nih.gov/geoprofiles/77629993. Accessed 11 Dec 2020.

  29. Al Saleh S, Al Mulla F, Luqmani YA. NAGA - estrogen receptor alpha-silenced MCF7 breast cancer cells (identifier 77619624). GEO Profiles 2011. https://www.ncbi.nlm.nih.gov/geoprofiles/?term=77619624. Accessed 11 Dec 2020.

  30. Carroll JS, Meyer CA, Song J, Li W, Geistlinger TR, Eeckhoute J, et al. Genome-wide analysis of estrogen receptor binding sites. Nat Genet. 2006;38(11):1289–97. https://doi.org/10.1038/ng1901.

    Article  PubMed  CAS  Google Scholar 

  31. Carroll JS, Meyer CA, Song J, Li W, Geistlinger TR, Eeckhoute J, et al. GLA - estrogen effect on breast cancer cell line: time course (identifier 53367630). GEO Profiles. 2006; https://www.ncbi.nlm.nih.gov/geoprofiles/?term=53367630. Accessed 11 Dec 2020.

  32. Doktorova M, Symons JL, Levental I. Structural and functional consequences of reversible lipid asymmetry in living membranes. Nat Chem Biol. 2020;16(12):1321–30. https://doi.org/10.1038/s41589-020-00688-0.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Tamura Y, Kawano S, Endo T. Lipid homeostasis in mitochondria. Biol Chem. 2020;401(6-7):821–33. https://doi.org/10.1515/hsz-2020-0121.

    Article  PubMed  CAS  Google Scholar 

  34. Shibabaw T. Omega-3 polyunsaturated fatty acids: anti-inflammatory and anti-hypertriglyceridemia mechanisms in cardiovascular disease. Mol Cell Biochem. 2020;476(2):993–1003.

Download references

Acknowledgments

Not applicable.

Funding

GOLDN biospecimens and baseline phenotype data were collected with funding from National Heart, Lung, and Blood Institute (NHLBI) grant U01HL072524. This work was also supported in part by grants from the National Center for Advancing Translational Sciences (UL1TR001998), the National Institute of Diabetes and Digestive and Kidney Diseases (1R01DK113625), the National Institute on Aging (1R56AG057589), the NHLBI (R01HL091357), the National Institute of General Medical Sciences (1P20GM130456), and the American Heart Association (15SDG25760020). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

Author information

Authors and Affiliations

  1. Department of Biostatistics, University of Kentucky, 725 Rose St, Multidisciplinary Science Building, Suite 205, Lexington, KY, 40536, USA

    Emily Slade, Kevin Xie & David W. Fardo

  2. Department of Epidemiology, University of Alabama at Birmingham, Birmingham, AL, USA

    Marguerite R. Irvin

  3. College of Public Health, University of Kentucky, Lexington, KY, USA

    Donna K. Arnett & Steven A. Claas

  4. West Coast Metabolomics Center, University of California, Davis, CA, USA

    Tobias Kind

  5. Department of Pharmaceutical Sciences, University of Kentucky, Lexington, KY, USA

    Gregory A. Graf

Authors
  1. Emily Slade
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Marguerite R. Irvin
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Kevin Xie
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Donna K. Arnett
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Steven A. Claas
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Tobias Kind
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. David W. Fardo
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Gregory A. Graf
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

DKA and MRI designed the GOLDN studies and oversaw collection of all data. ES, DKA, DWF, and GAG conceived and designed the analyses presented here. ES and KX conducted the analyses presented here. ES, GAG, and SAC drafted the manuscript. TK oversaw the generation and quality control of the lipidomics data. All authors contributed to the interpretation of the data and read and approved the final manuscript.

Corresponding author

Correspondence to Emily Slade.

Ethics declarations

Ethics approval and consent to participate

Ethics approval and consent to participate in the GOLDN study was approved by the University of Minnesota’s Human Research Protection Program Institutional Review Board, University of Alabama at Birmingham’s Institutional Review Board for Human Use, University of Utah’s Institutional Review Board, Washington University in St Louis’s Human Research Protection Office Institutional Review Board, and the University of Kentucky’s Office of Research Integrity Institutional Review Board. Written informed consent was obtained from each participant.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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.

Continuation of Table 4. Associations between lipid signal intensity with age, sex, and their interaction (includes all lipid species).

Additional file 2.

Correlation between total class intensities.

Additional file 3.

Diagnostic plots for regression models of lipid classes and subclasses.

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Slade, E., Irvin, M.R., Xie, K. et al. Age and sex are associated with the plasma lipidome: findings from the GOLDN study. Lipids Health Dis 20, 30 (2021). https://doi.org/10.1186/s12944-021-01456-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12944-021-01456-2

Keywords

Lipids in Health and Disease

ISSN: 1476-511X