High-coverage plasma lipidomics reveals novel sex-specific lipidomic fingerprints of age and BMI: Evidence from two large population cohort studies.

Obesity and related metabolic diseases show clear sex-related differences. The growing burden of these diseases calls for better understanding of the age- and sex-related metabolic consequences. High-throughput lipidomic analyses of population-based cohorts offer an opportunity to identify disease-risk-associated biomarkers and to improve our understanding of lipid metabolism and biology at a population level. Here, we comprehensively examined the relationship between lipid classes/subclasses and molecular species with age, sex, and body mass index (BMI). Furthermore, we evaluated sex specificity in the association of the plasma lipidome with age and BMI. Some 747 targeted lipid measures, representing 706 molecular lipid species across 36 classes/subclasses, were measured using a high-performance liquid chromatography coupled mass spectrometer on a total of 10,339 participants from the Australian Diabetes, Obesity and Lifestyle Study (AusDiab), with 563 lipid species being validated externally on 4,207 participants of the Busselton Health Study (BHS). Heat maps were constructed to visualise the relative differences in lipidomic profile between men and women. Multivariable linear regression analyses, including sex-interaction terms, were performed to assess the associations of lipid species with cardiometabolic phenotypes. Associations with age and sex were found for 472 (66.9%) and 583 (82.6%) lipid species, respectively. We further demonstrated that age-associated lipidomic fingerprints differed by sex. Specific classes of ether-phospholipids and lysophospholipids (calculated as the sum composition of the species within the class) were inversely associated with age in men only. In analyses with women alone, higher triacylglycerol and lower lysoalkylphosphatidylcholine species were observed among postmenopausal women compared with premenopausal women. We also identified sex-specific associations of lipid species with obesity. Lysophospholipids were negatively associated with BMI in both sexes (with a larger effect size in men), whilst acylcarnitine species showed opposing associations based on sex (positive association in women and negative association in men). Finally, by utilising specific lipid ratios as a proxy for enzymatic activity, we identified stearoyl CoA desaturase (SCD-1), fatty acid desaturase 3 (FADS3), and plasmanylethanolamine Δ1-desaturase activities, as well as the sphingolipid metabolic pathway, as constituent perturbations of cardiometabolic phenotypes. Our analyses elucidate the effect of age and sex on lipid metabolism by offering a comprehensive view of the lipidomic profiles associated with common cardiometabolic risk factors. These findings have implications for age- and sex-dependent lipid metabolism in health and disease and suggest the need for sex stratification during lipid biomarker discovery, establishing biological reference intervals for assessment of disease risk.

Correction: High-coverage plasma lipidomics reveals novel sex-specific lipidomic fingerprints of age and BMI: Evidence from two large population cohort studies.

[This corrects the article DOI: 10.1371/journal.pbio.3000870.].

APOE ε2 resilience for Alzheimer's disease is mediated by plasma lipid species: Analysis of three independent cohort studies.

INTRODUCTION: The apolipoprotein E (APOE) genotype is the strongest genetic risk factor for late-onset Alzheimer's disease. However, its effect on lipid metabolic pathways, and their mediating effect on disease risk, is poorly understood. METHODS: We performed lipidomic analysis on three independent cohorts (the Australian Imaging, Biomarkers and Lifestyle [AIBL] flagship study, n = 1087; the Alzheimer's Disease Neuroimaging Initiative [ADNI] 1 study, n = 819; and the Busselton Health Study [BHS], n = 4384), and we defined associations between APOE ε2 and ε4 and 569 plasma/serum lipid species. Mediation analysis defined the proportion of the treatment effect of the APOE genotype mediated by plasma/serum lipid species. RESULTS: A total of 237 and 104 lipid species were associated with APOE ε2 and ε4, respectively. Of these 68 (ε2) and 24 (ε4) were associated with prevalent Alzheimer's disease. Individual lipid species or lipidomic models of APOE genotypes mediated up to 30% and 10% of APOE ε2 and ε4 treatment effect, respectively. DISCUSSION: Plasma lipid species mediate the treatment effect of APOE genotypes on Alzheimer's disease and as such represent a potential therapeutic target.

NormalizeMets: assessing, selecting and implementing statistical methods for normalizing metabolomics data.

INTRODUCTION: In metabolomics studies, unwanted variation inevitably arises from various sources. Normalization, that is the removal of unwanted variation, is an essential step in the statistical analysis of metabolomics data. However, metabolomics normalization is often considered an imprecise science due to the diverse sources of variation and the availability of a number of alternative strategies that may be implemented. OBJECTIVES: We highlight the need for comparative evaluation of different normalization methods and present software strategies to help ease this task for both data-oriented and biological researchers. METHODS: We present NormalizeMets-a joint graphical user interface within the familiar Microsoft Excel and freely-available R software for comparative evaluation of different normalization methods. The NormalizeMets R package along with the vignette describing the workflow can be downloaded from https://cran.r-project.org/web/packages/NormalizeMets/ . The Excel Interface and the Excel user guide are available on https://metabolomicstats.github.io/ExNormalizeMets . RESULTS: NormalizeMets allows for comparative evaluation of normalization methods using criteria that depend on the given dataset and the ultimate research question. Hence it guides researchers to assess, select and implement a suitable normalization method using either the familiar Microsoft Excel and/or freely-available R software. In addition, the package can be used for visualisation of metabolomics data using interactive graphical displays and to obtain end statistical results for clustering, classification, biomarker identification adjusting for confounding variables, and correlation analysis. CONCLUSION: NormalizeMets is designed for comparative evaluation of normalization methods, and can also be used to obtain end statistical results. The use of freely-available R software offers an attractive proposition for programming-oriented researchers, and the Excel interface offers a familiar alternative to most biological researchers. The package handles the data locally in the user's own computer allowing for reproducible code to be stored locally.

Development and validation of a plasmalogen score as an independent modifiable marker of metabolic health: population based observational studies and a placebo-controlled cross-over study.

BACKGROUND: Decreased levels of circulating ethanolamine plasmalogens [PE(P)], and a concurrent increase in phosphatidylethanolamine (PE) are consistently reported in various cardiometabolic conditions. Here we devised, a plasmalogen score (Pls Score) that mirrors a metabolic signal that encompasses the levels of PE(P) and PE and captures the natural variation in circulating plasmalogens and perturbations in their metabolism associated with disease, diet, and lifestyle. METHODS: We utilised, plasma lipidomes from the Australian Obesity, Diabetes and Lifestyle study (AusDiab; n = 10,339, 55% women) a nationwide cohort, to devise the Pls Score and validated this in the Busselton Health Study (BHS; n = 4,492, 56% women, serum lipidome) and in a placebo-controlled crossover trial involving Shark Liver Oil (SLO) supplementation (n = 10, 100% men). We examined the association of the Pls Score with cardiometabolic risk factors, type 2 diabetes mellitus (T2DM), cardiovascular disease and all-cause mortality (over 17 years). FINDINGS: In a model, adjusted for age, sex and BMI, individuals in the top quintile of the Pls Score (Q5) relative to Q1 had an OR of 0.31 (95% CI 0.21-0.43), 0.39 (95% CI 0.25-0.61) and 0.42 (95% CI 0.30-0.57) for prevalent T2DM, incident T2DM and prevalent cardiovascular disease respectively, and a 34% lower mortality risk (HR = 0.66; 95% CI 0.56-0.78). Significant associations between diet and lifestyle habits and Pls Score exist and these were validated through dietary supplementation of SLO that resulted in a marked change in the Pls Score. INTERPRETATION: The Pls Score as a measure that captures the natural variation in circulating plasmalogens, was not only inversely related to cardiometabolic risk and all-cause mortality but also associate with diet and lifestyle. Our results support the potential utility of the Pls Score as a biomarker for metabolic health and its responsiveness to dietary interventions. Further research is warranted to explore the underlying mechanisms and optimise the practical implementation of the Pls Score in clinical and population settings. FUNDING: National Health and Medical Research Council (NHMRC grant 233200), National Health and Medical Research Council of Australia (Project grant APP1101320), Health Promotion Foundation of Western Australia, and National Health and Medical Research Council of Australia Senior Research Fellowship (#1042095).

A lipidomic based metabolic age score captures cardiometabolic risk independent of chronological age.

BACKGROUND: Metabolic ageing biomarkers may capture the age-related shifts in metabolism, offering a precise representation of an individual's overall metabolic health. METHODS: Utilising comprehensive lipidomic datasets from two large independent population cohorts in Australia (n = 14,833, including 6630 males, 8203 females), we employed different machine learning models, to predict age, and calculated metabolic age scores (mAge). Furthermore, we defined the difference between mAge and age, termed mAgeΔ, which allow us to identify individuals sharing similar age but differing in their metabolic health status. FINDINGS: Upon stratification of the population into quintiles by mAgeΔ, we observed that participants in the top quintile group (Q5) were more likely to have cardiovascular disease (OR = 2.13, 95% CI = 1.62-2.83), had a 2.01-fold increased risk of 12-year incident cardiovascular events (HR = 2.01, 95% CI = 1.45-2.57), and a 1.56-fold increased risk of 17-year all-cause mortality (HR = 1.56, 95% CI = 1.34-1.79), relative to the individuals in the bottom quintile group (Q1). Survival analysis further revealed that men in the Q5 group faced the challenge of reaching a median survival rate due to cardiovascular events more than six years earlier and reaching a median survival rate due to all-cause mortality more than four years earlier than men in the Q1 group. INTERPRETATION: Our findings demonstrate that the mAge score captures age-related metabolic changes, predicts health outcomes, and has the potential to identify individuals at increased risk of metabolic diseases. FUNDING: The specific funding of this article is provided in the acknowledgements section.

Removing unwanted variation from large-scale RNA sequencing data with PRPS.

Accurate identification and effective removal of unwanted variation is essential to derive meaningful biological results from RNA sequencing (RNA-seq) data, especially when the data come from large and complex studies. Using RNA-seq data from The Cancer Genome Atlas (TCGA), we examined several sources of unwanted variation and demonstrate here how these can significantly compromise various downstream analyses, including cancer subtype identification, association between gene expression and survival outcomes and gene co-expression analysis. We propose a strategy, called pseudo-replicates of pseudo-samples (PRPS), for deploying our recently developed normalization method, called removing unwanted variation III (RUV-III), to remove the variation caused by library size, tumor purity and batch effects in TCGA RNA-seq data. We illustrate the value of our approach by comparing it to the standard TCGA normalizations on several TCGA RNA-seq datasets. RUV-III with PRPS can be used to integrate and normalize other large transcriptomic datasets coming from multiple laboratories or platforms.

Metabolic phenotyping of BMI to characterize cardiometabolic risk: evidence from large population-based cohorts.

Obesity is a risk factor for type 2 diabetes and cardiovascular disease. However, a substantial proportion of patients with these conditions have a seemingly normal body mass index (BMI). Conversely, not all obese individuals present with metabolic disorders giving rise to the concept of "metabolically healthy obese". We use lipidomic-based models for BMI to calculate a metabolic BMI score (mBMI) as a measure of metabolic dysregulation associated with obesity. Using the difference between mBMI and BMI (mBMIΔ), we identify individuals with a similar BMI but differing in their metabolic health and disease risk profiles. Exercise and diet associate with mBMIΔ suggesting the ability to modify mBMI with lifestyle intervention. Our findings show that, the mBMI score captures information on metabolic dysregulation that is independent of the measured BMI and so provides an opportunity to assess metabolic health to identify "at risk" individuals for targeted intervention and monitoring.

Challenges and opportunities for prevention and removal of unwanted variation in lipidomic studies.

Large 'omics studies are of particular interest to population and clinical research as they allow elucidation of biological pathways that are often out of reach of other methodologies. Typically, these information rich datasets are produced from multiple coordinated profiling studies that may include lipidomics, metabolomics, proteomics or other strategies to generate high dimensional data. In lipidomics, the generation of such data presents a series of unique technological and logistical challenges; to maximize the power (number of samples) and coverage (number of analytes) of the dataset while minimizing the sources of unwanted variation. Technological advances in analytical platforms, as well as computational approaches, have led to improvement of data quality - especially with regard to instrumental variation. In the small scale, it is possible to control systematic bias from beginning to end. However, as the size and complexity of datasets grow, it is inevitable that unwanted variation arises from multiple sources, some potentially unknown and out of the investigators control. Increases in cohort size and complexity have led to new challenges in sample collection, handling, storage, and preparation. If not considered and dealt with appropriately, this unwanted variation may undermine the quality of the data and reliability of any subsequent analysis. Here we review the various experimental phases where unwanted variation may be introduced and review general strategies and approaches to handle this variation, specifically addressing issues relevant to lipidomics studies.

Comprehensive genetic analysis of the human lipidome identifies loci associated with lipid homeostasis with links to coronary artery disease.

We integrated lipidomics and genomics to unravel the genetic architecture of lipid metabolism and identify genetic variants associated with lipid species putatively in the mechanistic pathway for coronary artery disease (CAD). We quantified 596 lipid species in serum from 4,492 individuals from the Busselton Health Study. The discovery GWAS identified 3,361 independent lipid-loci associations, involving 667 genomic regions (479 previously unreported), with validation in two independent cohorts. A meta-analysis revealed an additional 70 independent genomic regions associated with lipid species. We identified 134 lipid endophenotypes for CAD associated with 186 genomic loci. Associations between independent lipid-loci with coronary atherosclerosis were assessed in ∼456,000 individuals from the UK Biobank. Of the 53 lipid-loci that showed evidence of association (P < 1 × 10-3), 43 loci were associated with at least one lipid endophenotype. These findings illustrate the value of integrative biology to investigate the aetiology of atherosclerosis and CAD, with implications for other complex diseases.

Lipidomic Signatures of Changes in Adiposity: A Large Prospective Study of 5849 Adults from the Australian Diabetes, Obesity and Lifestyle Study.

Lipid metabolism is tightly linked to adiposity. Comprehensive lipidomic profiling offers new insights into the dysregulation of lipid metabolism in relation to weight gain. Here, we investigated the relationship of the human plasma lipidome and changes in waist circumference (WC) and body mass index (BMI). Adults (2653 men and 3196 women), 25-95 years old who attended the baseline survey of the Australian Diabetes, Obesity and Lifestyle Study (AusDiab) and the 5-year follow-up were enrolled. A targeted lipidomic approach was used to quantify 706 distinct molecular lipid species in the plasma samples. Multiple linear regression models were used to examine the relationship between the baseline lipidomic profile and changes in WC and BMI. Metabolic scores for change in WC were generated using a ridge regression model. Alkyl-diacylglycerol such as TG(O-50:2) [NL-18:1] displayed the strongest association with change in WC (ß-coefficient = 0.125 cm increment per SD increment in baseline lipid level, p = 2.78 × 10-11. Many lipid species containing linoleate (18:2) fatty acids were negatively associated with both WC and BMI gain. Compared to traditional models, multivariate models containing lipid species identify individuals at a greater risk of gaining WC: top quintile relative to bottom quintile (odds ratio, 95% CI = 5.4, 3.8-6.6 for women and 2.3, 1.7-3.0 for men). Our findings define metabolic profiles that characterize individuals at risk of weight gain or WC increase and provide important insight into the biological role of lipids in obesity.

Concordant peripheral lipidome signatures in two large clinical studies of Alzheimer's disease.

Changes to lipid metabolism are tightly associated with the onset and pathology of Alzheimer's disease (AD). Lipids are complex molecules comprising many isomeric and isobaric species, necessitating detailed analysis to enable interpretation of biological significance. Our expanded targeted lipidomics platform (569 species across 32 classes) allows for detailed lipid separation and characterisation. In this study we examined peripheral samples of two cohorts (AIBL, n = 1112 and ADNI, n = 800). We are able to identify concordant peripheral signatures associated with prevalent AD arising from lipid pathways including; ether lipids, sphingolipids (notably GM3 gangliosides) and lipid classes previously associated with cardiometabolic disease (phosphatidylethanolamine and triglycerides). We subsequently identified similar lipid signatures in both cohorts with future disease. Lastly, we developed multivariate lipid models that improved classification and prediction. Our results provide a holistic view between the lipidome and AD using a comprehensive approach, providing targets for further mechanistic investigation.