
==== Front
Cardiovasc Diabetol
Cardiovasc Diabetol
Cardiovascular Diabetology
1475-2840
BioMed Central London

39261922
2429
10.1186/s12933-024-02429-2
Research
Mendelian randomization analysis reveals causal effects of blood lipidome on gestational diabetes mellitus
Dong Yao 123
Hu An-qun 4
Han Bai-xue 123
Cao Meng-ting 123
Liu Hai-yan 4
Li Zong-guang 4
Li Qing 5
Zheng Ying-jie yjzheng@fudan.edu.cn

123
1 https://ror.org/013q1eq08 grid.8547.e 0000 0001 0125 2443 Department of Epidemiology, School of Public Health, Fudan University, 130 Dong-an Rd., Shanghai, 200032 China
2 grid.8547.e 0000 0001 0125 2443 Laboratory of Public Health Safety, Ministry of Education, School of Public Health, Fudan University, Shanghai, 200032 China
3 grid.8547.e 0000 0001 0125 2443 Key Laboratory for Health Technology Assessment, National Commission of Health and Family Planning, Fudan University, Shanghai, 200032 China
4 Department of Clinical Laboratory, Anqing Municipal Hospital, Anqing, 246003 China
5 Department of Obstetrics and Gynecology, Anqing Municipal Hospital, Anqing, 246003 China
11 9 2024
11 9 2024
2024
23 3358 7 2024
2 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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-nc-nd/4.0/.
Background

Observational studies have revealed associations between maternal lipid metabolites and gestational diabetes mellitus (GDM). However, whether these associations are causal remain uncertain.

Objective

To evaluate the causal relationship between lipid metabolites and GDM.

Methods

A two-sample Mendelian randomization (MR) analysis was performed based on summary statistics. Sensitivity analyses, validation analyses and reverse MR analyses were conducted to assess the robustness of the MR results. Additionally, a phenome-wide MR (Phe-MR) analysis was performed to evaluate potential side effects of the targeted lipid metabolites.

Results

A total of 295 lipid metabolites were included in this study, 29 of them had three or more instrumental variables (IVs) suitable for sensitivity analyses. The ratio of triglycerides to phosphoglycerides (TG_by_PG) was identified as a potential causal biomarker for GDM (inverse variance weighted (IVW) estimate: odds ratio (OR) = 2.147, 95% confidential interval (95% CI) 1.415–3.257, P = 3.26e−4), which was confirmed by validation and reverse MR results. Two other lipid metabolites, palmitoyl sphingomyelin (d18:1/16:0) (PSM(d18:1/16:0)) (IVW estimate: OR = 0.747, 95% CI 0.583–0.956, P = 0.021) and triglycerides in very small very low-density lipoprotein (XS_VLDL_TG) (IVW estimate: OR = 2.948, 95% CI 1.197–5.215, P = 0.015), were identified as suggestive potential biomarkers for GDM using a conventional cut-off P-value of 0.05. Phe-MR results indicated that lowering TG_by_PG had detrimental effects on two diseases but advantageous effects on the other 13 diseases.

Conclusion

Genetically predicted elevated TG_by_PG are causally associated with an increased risk of GDM. Side-effect profiles indicate that TG_by_PG might be a target for GDM prevention, though caution is advised due to potential adverse effects on other conditions.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12933-024-02429-2.

Keywords

Gestational diabetes mellitus
Lipid metabolites
Two-sample Mendelian randomization
Phenome-wide Mendelian randomization
National Natural Science Foundation of China82173582 issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcIntroduction

Gestational diabetes mellitus (GDM), characterized by hyperglycemia developing during pregnancy, is the most frequent pregnancy complication [1]. GDM increases the risk of multiple short-term and long-term adverse complications in both mother and offspring, including preterm birth, hypertension, obesity, impaired glucose metabolism, and cardiovascular disease [2, 3], but relatively little is known about the modifiable targets for prevention and treatment strategies of GDM. Randomized clinical trials have shown that treatment (diet control and exercise) for mild GDM (diagnosed between 24+0 ~ 30+6 weeks of pregnancy) does not significantly reduce adverse outcomes including type 2 diabetes mellitus (T2DM), metabolic syndrome, and obesity within seven years postpartum [4, 5]. Additional research is necessary to comprehend the causes of GDM and to identify potential intervention targets.

The onset of GDM is closely related to lipid metabolism. Metabolic changes, such as physiological insulin resistance [6] and altered lipid metabolism, [7] occur during pregnancy, but these changes are more pronounced in women with GDM [8, 9], suggesting potential metabolic dysfunction. Lipids play a crucial role in human metabolism, and understanding their relationship with metabolic diseases may help us elucidate the pathological mechanism and identify potential intervention targets. Advances in high-throughput methods have made lipidomics a powerful tool for systematically uncovering lipidomic biomarkers and understanding the pathogenesis of GDM [10].

An increasing number of observational studies have identified connections between certain lipid metabolites—such as fatty acyls (FAs), glycerolipids (GLs), glycerophospholipids (GPs), sphingolipids (SPs), and lipoproteins—and the risk of GDM [11–13]. Additionally, individual lipid species within certain lipid class may have opposing associations with GDM [14–16], highlighting the necessity of lipidomics research to understand the complex biological effects of individual lipid species on the development of GDM. However, due to the complex correlation and interaction between glucose metabolism and lipid metabolism [17], it is challenging to determine the causal relationship and the direction of the association between lipid profiles and GDM, especially in observational studies that are prone to confounding and reverse causation. The causal effect of each lipid metabolite on GDM remains unclear, and potential side effects of targeting lipidomic biomarkers for intervention have not been explored. With the rapid increase in genome-wide association studies (GWAS) focusing on blood lipidomes [18–21] and GDM [22], Mendelian randomization (MR) analysis can be utilized to evaluate the causal relationship between lipidomes and GDM using GWAS summary statistics.

In this context, we conducted a two-sample MR study to fulfill two objectives. First, we aimed to examine the causal relationship between lipid metabolites and GDM risk in order to identify potential lipid intervention targets and gain deeper insights into the pathological mechanisms of GDM. Second, prior to clinical trials, we applied a phenome-wide MR (Phe-MR) analysis aimed to assess unanticipated adverse effects of potential lipid biomarker intervention, offering a comprehensive appraisal of their clinical safety.

Methods

Study design

This study was conducted and reported according to the STROBE-MR guideline [19]. We first conducted a two-sample MR analysis to evaluate the causal relationships between 295 lipid metabolites and GDM, and then an extended Phe-MR analysis of 1,210 non-GDM diseases to predict a wide range of potential side effects associated with targeting identified lipid metabolites (Fig. 1).Fig. 1 Conceptual framework of Mendelian randomization study. The study consists of two parts. a Part 1, we assessed the causality for the association between 295 lipid metabolites and the risk of GDM. b Part 2, we investigated a broad spectrum of side effects associated with targeting identified lipid metabolites in 1210 non-GDM diseases. GWAS = genome-wide association study, MR = Mendelian randomization, N in a = number of lipid metabolites, N in b = number of phenotypes

GWAS summary statistic for lipid metabolites

Summary statistics of single nucleotide polymorphisms (SNPs) associated with the blood lipidome were taken from three large-scale GWASs of European individuals [18–21] (Table 1). (1) Surendran et al. analyzed 271 lipid metabolic traits in 14,296 participants (8,455 from the INTERVAL study, and 5,841 from the EPIC-Norfolk study) with up to 87,696,888 imputed autosomal SNPs [18]. (2) Cadby et al. conducted a GWAS using data from the Busselton Health Study (BHS) involving 596 lipid species (539 of them found to have statistically significant associations) and 4,492 individuals of predominantly European ancestry, with 39,117,105 SNPs available for analysis after imputation [19]. (3) Karjalainen et al. conducted a meta-analysis of 213 lipid and lipoprotein parameters across 33 cohorts, with up to 13,389,637 imputed autosomal SNPs included in the meta-analysis [20], and our study specifically utilized data from 92,664 non-Finnish European participants from 21 of these cohorts, with detailed cohort information provided in Table S1.Table 1 Characteristics of blood lipid metabolites GWAS studies

References	Cohort(s)	Population	Sample size	Mean Age (y)	Female (%)	No. SNP	Platform	Metabolites*	Sample type	
Surendran et al. (2022) [18]	INTERVAL	European ancestor	8,455	44.0	48.8	87,696,888	MS	271 (98)	Fasting plasma	
	EPIC-Norfolk		5,841	59.8	53.3					
Cadby et al. (2022) [19]	BHS	Australian	4,492	50.8	56.0	39,117,105	MS	539 (135)	Fasting serum	
Karjalainen et al. (2024) [20]	Meta-GWAS	Non-Finnish European	92,664 (21 cohorts)	46.3	56.3	13,389,637	NMR	213 (62)	Mostly fasting serum/plasma	
Total								1023 (295)		
BHS the Busselton Health Study, EPIC-Norfolk the EPIC-Norfolk study, GWAS genome-wide association study, INTERVAL the INTERVAL study, MS mass spectrometry, NMR  nuclear magnetic resonance

*The numbers in brackets are the number of lipid metabolites included in the following analyses

GWAS summary statistics for GDM

The GWAS summary statistics for GDM were obtained from the Release 10 results provided by the FinnGen consortium. This GWAS included data from 14,718 GDM cases (identified from registry data using ICD-9 and 10 codes O24.4) and 215,592 parous female controls. All participants in this study were of European (Finnish) ancestry [22].

Instrument variables and lipid metabolites selection

Three key assumptions of MR must be satisfied in order to ensure the validity of the instrument variables (IVs): (1) IVs are strongly associated with lipid metabolites; (2) IVs are independent of any observed or unobserved confounders of lipid metabolites-GDM associations; (3) IVs affect GDM only through their effect in lipid metabolites without any alternative pathways.

To meet the above assumptions, IVs were selected based on rigorous criteria. Specifically, (1) genetic variants were screened with genome-wide significance (P < 5e−8); (2) for each lipid metabolites, genetic variants were clumped to retain independent SNPs using a linkage disequilibrium threshold of r2 = 0.001 and a clumping window size of 500 kb; (3) The strength of the selected IVs was evaluated using the F-statistic, calculated according to Eq. 1, where r2 represents the proportion of variance explained by the IVs, N is the sample size, and k is the number of IVs. IVs of F-statistic < 10 were considered weak and were excluded in the subsequent analysis; (4) SNP associated with more than five lipid metabolites were excluded to avoid pleiotropy effect.1 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F-statistic=\frac{{r}^{2}\times (N-1-k)}{(1-{r}^{2})\times k}$$\end{document}

We also eliminated metabolites with less than 0.5% variance explained by IVs to ensure adequate statistical power for legitimate causal inferences. Lipid metabolites retained were categorized based on the number of IVs: those with three or more IVs were designated as tier 1, for which sensitivity analyses and robust causal association estimates were available [23]. Lipid metabolites with less than three IVs were classified as tier 2. The main results were based on tier 1 lipid metabolites, while tier 2 lipid metabolites were presented as suggestive findings.

Mendelian randomization analysis

A two-sample MR analysis was carried out to assess the causal effect of lipid metabolites on GDM, using the inverse variance weighted (IVW) method as the primary approach. The MR-Egger regression, weighted median, weighted mode methods and constrained maximum likelihood and model averaging (cML-MA) method were conducted to assess the robustness of the IVW estimate. The absence of horizontal pleiotropy is a prerequisite for the validity of IVW results [24]. Assuming the instrument strength is independent of direct effect (InSIDE), MR-Egger regression assesses pleiotropy via its intercept term [25]. The result of the MR-Egger regression is consistent with IVW if the intercept term is equal to zero, which suggests that horizontal pleiotropy does not exist [25]. The weighted mode estimate demonstrates higher power, less bias, and lower type I error rates compared to MR-Egger regression if the InSIDE hypothesis is violated [26, 27]. By removing significant outliers, the MR-PRESSO analysis detects and attempts to reduce horizontal pleiotropy, but it also depends on InSIDE assumptions and requires that at least 50% of the genetic variants are valid IVs [28]. The weighted median approach can yield a reliable causal estimate when up to 50% of genetic variants were invalid [29]. To control correlated and uncorrelated pleiotropic effects, we employed cML-MA, which does not rely on the InSIDE assumption [30].

Sensitivity analyses

A series of sensitivity analyses were carried out in order to assess the robustness of causal effects following the main MR analysis. (1) Cochran’s Q test was conducted to measure heterogeneity across IVs [31]. A fixed-effects model was employed when no heterogeneity was detected among IVs, while a random-effects IVW model was utilized otherwise. (2) As mentioned above, MR-Egger regression estimates an intercept term of the IVW test, with any deviation from zero indicating potential directional pleiotropy [24]. (3) Additionally, a leave-one-out analysis re-evaluates the MR association by sequentially excluding each individual IV [32].

Validation analyses and test for reverse causality

Significant associations were validated using summary data from UK Biobank (UKB) obtained from the IEU OpenGWAS platform (https://gwas.mrcieu.ac.uk/). We further performed a reverse MR analysis (i.e., GDM as exposure, and the identified lipid metabolites as outcome) utilizing SNPs that were linked to GDM as IVs to investigate if GDM had any causal impact on the identified lipid metabolites. As summary data released by Karjalainen et al. did not contain all lipid metabolites–SNP associations and was not qualified as an outcome dataset when conducting MR analysis, we used summary data from UKB as mentioned above to verify the reverse causation between significant lipid metabolites and GDM. All IV screening criteria, MR analysis strategies, and sensitivity analysis methods used in the validation and reverse MR analysis were consistent with those applied in the primary MR analysis.

Colocalization analysis

We performed colocalization analysis using validation dataset to detect whether the identified lipid metabolites and GDM share common causal variant(s). For each lipid metabolite-GDM pair, we set the window size to ± 250 kb centered on each IV. The evidence for colocalization was evaluated by calculating the posterior probability (PP) for hypothesis 4 (PP.H4), which posits that the same causal variant(s) are responsible for the associations observed between the lipid metabolites and GDM. A threshold of PP.H4 > 0.8 was applied to indicate strong evidence of colocalization.

Phenome-wide Mendelian randomization analyses

Phe-MR analysis was conducted to evaluate the potential side effects associated with hypothetical interventions aimed at reducing GDM by targeting specific lipid metabolites. Summary GWAS data of 2,447 disease traits with 412,181 Finnish participants were acquired from FinnGen consortium with 21,311,942 variants in the FinnGen release 10 results (https://www.finngen.fi/en/access_results). Disease traits were defined in terms of “Phenotype Description” in the original FinnGen data, and were organized into ICD code format to facilitate correspondence to specific diseases in different studies. Disease traits with fewer than 500 cases were excluded due to concerns about statistical power.

Additionally, we selected representative phenotypes to minimize inherent redundancy between PheCodes, thereby improving the interpretability of the results. Ultimately, 1,210 non-GDM disease traits were included in the Phe-MR analysis to further investigate the potential side effects of GDM-related metabolites (Fig. 1; Table S2). Genetic variants for GDM-related lipidomic biomarkers were derived from the same GWAS as in the main MR analysis. The final Phe-MR results were standardized based on the associations between lipid metabolites and GDM, reflecting a change in lipid metabolite levels equivalent to a 10% reduction in GDM risk. This standardization served three purposes: (1) to provide a clinical context (i.e., what side effects might arise if identified biomarkers are targeted therapeutically for GDM?); (2) to harmonize the directionality among identified lipid metabolites that could either increase or decrease GDM risk; and (3) to facilitate a direct comparison of the magnitude and directionality of the side effects.

Statistics and reproducibility

All analyses were conducted using R version 4.3.0. The strength of causal associations was evaluated using odds ratios (ORs) and corresponding 95% confidence intervals (CIs). An observed P-value < 1.92e−3 (Bonferroni-corrected significance threshold calculated as 0.05 divided by 26 [for 26 tier 1 lipid metabolites]) was considered as statistically significant for a potential causal association [33], and P-value < 0.05 was taken as suggestive significant threshold in this study. All P-values were two-sided. Tier 1 lipid metabolites that reach the significance threshold were considered to have potential causal relationships with GDM. Lipid metabolites from both tier 1 and tier 2 that met the suggestive threshold were considered as suggestive findings; however, these lipid metabolites lack robust estimation and therefore need to be treated with caution when reporting.

For the Phe-MR analysis, a significance threshold of P < 0.01 was employed, given the primary objective of screening for potential side effects of the target metabolites. This threshold was chosen to minimize false negatives, emphasizing the identification of potential associations for further investigation.

Phenotypic variance of specific lipid metabolite explained by relevant IVs were calculated using the “gtx” package (downloaded from https://www.rdocumentation.org/packages/gtx/versions/0.0.8). MR analysis and clumping were carried out using ‘TwoSampleMR’ package [23] and ‘MendelianRandomization’ package [34]. cML-MA method was conducted using ‘MRcML’ package [30]. Colocalization was carried out using the “coloc” package (https://github.com/chr1swallace/coloc).

Results

Strength of the genetic instruments for lipid metabolites

A total of 1023 lipid metabolites and 12,408 unique SNPs were identified from three GWASs. After filtering, 793 lipid metabolites were retained, and 409 SNPs qualified as IVs. After harmonizing effect alleles between the exposure and outcome datasets, 396 IVs were available for 790 lipid metabolites. To mitigate potential horizontal pleiotropy, SNPs associated with more than five lipid metabolites were excluded, resulting in 295 lipid metabolites and 239 SNPs included in the main MR analysis (Fig. 2). These lipid metabolites were categorized into eight groups: FA, GL, GP, SP, SL, as well as apolipoproteins (Apo), lipid subclasses (e.g., total ceramide), and lipoprotein subclasses (e.g., total lipids in very large HDL). Among these, 26 lipid metabolites belong to the tier 1 group, as they had three or more IVs, enabling sensitivity analysis and providing robust causal association estimates. The remaining lipid metabolites were categorized as tier 2 group. The main results are based on the tier 1 group lipid metabolites, while tier 2 lipid metabolites are presented as suggestive findings.Fig. 2 The flow chart of instrumental variables selection. GWAS genome-wide association study, IVW inverse variance weighted, MR Mendelian randomization; SNP = single nucleotide polymorphism

Detailed information of the SNP data used as IVs in this MR study is listed in Tables S3 and S4, and an overview of lipides metabolites is available in Table S5. The range of metabolite variation explained by the IVs was found to be 0.514% to 26.864% (Table S5). The F statistics for the genetic instruments associated with blood metabolites ranged from 27.123 to 5977.797 (Table S3), which indicating that they were sufficiently powered to detect MR associations.

Causal effect of lipid metabolites on GDM risk

Associations between tier 1 and tier 2 lipid metabolites and the risk of GDM in the main MR analysis are illustrated in Fig. 3 and Table S6. A genetically determined high ratio of triglycerides to phosphoglycerides (TG_by_PG) was potentially causally associated with an increased risk of GDM (IVW estimate: OR = 2.147, 95% CI: 1.415–3.257, P = 3.26e-4). Additionally, tier 1 lipid metabolite palmitoyl sphingomyelin (d18:1/16:0) (PSM(d18:1/16:0)) (IVW estimate: OR = 0.747, 95% CI: 0.583–0.956, P = 0.021), and tier 2 lipid triglycerides in very small very low-density lipoprotein (XS_VLDL_TG) (IVW estimate: OR = 2.498, 95% CI: 1.197–5.215, P = 0.015) exhibited suggestive causal associations with GDM (Fig. 3, Table 2). Moreover, although the significance threshold was not reached, certain lipid subclasses, such as fatty acids (FAs) and glycerolipids (GLs), appeared to have stronger associations with GDM compared to other lipid subclasses (Fig. 3).Fig. 3 Circular Manhattan plot displaying the associations between lipid metabolites and the risk of GDM in the Mendelian randomization analysis. The dashed red line represents the Bonferroni-corrected significance threshold (P < 1.92e−3), while the inner line indicates the suggestive significance threshold (P < 0.05). (Suggestive) Significant lipid metabolites are labeled accordingly. Blue bars denote tier 1 lipid metabolites while grey bars denote tier 2 lipid metabolites. The 295 blood lipid metabolites are categorized and color-coded by lipid classes as shown in Fig. 1 and Table S5. Detailed results for the associations between blood metabolites and GDM by Mendelian randomization analysis are presented in Table S6. GDM gestational diabetes mellitus, TG_by_PG ratio of triglycerides to phosphoglycerides, PSM (d18:1/16:0) palmitoyl sphingomyelin (d18:1/16:0), XS_VLDL_TG triglycerides in very small VLDL

Table 2 MR results of significant lipid metabolites and GDM

Lipid metabolites	nsnp	IVW (Fixed-effect)	Cochran’s Q test	MR-Egger	
P-value	OR (95% CI)	Q	Q_pval	intercept	P-value	
TG_by_PG	4	3.26e-4	2.147 (1.415–3.257)	0.085	0.994	0.011	0.818	
PSM(d18:1/16:0)*	3	0.021	0.747 (0.583–0.956)	1.616	0.446	0.036	0.587	
XS_VLDL_TG*	2	0.015	2.948 (1.197–5.215)	0.006	0.936	/	/	
CI confidence interval, IVW inverse variance weighted, nsnp number of single nucleotide polymorphisms, OR odds ratio; PSM (d18:1/16:0) = palmitoyl sphingomyelin (d18:1/16:0); TG_by_PG = ratio of triglycerides to phosphoglycerides; XS_VLDL_TG = triglycerides in very small VLDL

*: suggestive findings

/: Unattainable

Sensitivity analyses results

The results of Cochran’s IVW Q test indicated significant heterogeneity of IVs associated with five lipid metabolites (L_HDL_TG, N-palmitoyl-sphingosine (d18:1/16:0), 5alpha-androstan-3beta,17beta-diol disulfate, Total Cer, epiandrosterone sulfate), and random-effects IVW models were conducted to evaluate causal associations between these five lipid metabolites and GDM (Table S6, Table S7). For other lipid metabolites, no significant heterogeneity in IVs was observed (Table 2, Table S7). Additionally, MR-Egger regression intercept analysis showed no significant directional horizontal pleiotropy in the associations between lipid metabolites and GDM (Table 2, Table S7). Leave-one-out analysis indicated that all tests are driven by multiple IVs (Figure S1).

Validation and reverse MR results

The relationship between TG_by_PG and GDM was validated using summary data from 114,999 individuals from UKB after excluding participants of non-European descent (downloaded from https://gwas.mrcieu.ac.uk/datasets/met-d-TG_by_PG/). Similarly, the association betweenXS_VLDL_TG and GDM was validated using summary data from 115,078 individuals in the UKB (data accessed from https://gwas.mrcieu.ac.uk/datasets/met-d-XS_VLDL_TG/). However, we were unable to verify the association between PSM(d18:1/16:0) and GDM as no other GWAS summary dataset for PSM(d18:1/16:0) was available for MR analysis.

Causal associations between TG_by_PG and GDM was validated using validation data (IVW estimate: OR = 1.125, 95% CI: 1.045–1.211, P = 0.002), although causal association between XS_VLDL_TG was not verified by the validation dataset (IVW estimate: OR = 1.033, 95% CI: 0.974–1.095, P = 0.279) (Table S8). Cochran’s IVW Q test showed no significant heterogeneity in TG_by_PG IVs (Q_pval = 0.977). Though the MR-Egger regression intercept analysis indicated the presence of horizontal pleiotropy, further MR-PRESSO analysis and cML-MA analysis found no significant outliers, and the association remained significant when accounting for associated pleiotropy (cML-MA estimate: OR = 1.126, 95% CI: 1.046–1.212, P = 0.002) (Table S8).

In reverse MR analyses, no significant association was found between GDM and TG_by_PG (P = 0.770) or XS_VLDL_TG (P = 0.670) (Table S9), indicating no potential reverse causation. Cochran’s IVW Q test showed no significant heterogeneity in GDM IVs (Q_pval > 0.05), and the MR-Egger regression intercept analysis did not find significant horizontal pleiotropy (P-value > 0.05) (Table S9).

Besides, strong colocalization evidence was observed between TG_by_PG and GDM (PP.H4 = 0.954), as shown in Table S10.

Phe-MR results

In the Phe-MR analysis using the IVW method, a total of 45 associations reached a significance threshold of P = 0.01. Sensitivity analyses with the methods of Cochran’s IVW Q test and MR-Egger revealed no significant heterogeneity or pleiotropy in these 45 associations (Table S11).

In brief, lowering TG_by_PG showed detrimental effects on the risk of peripheral retinal degeneration and intrahepatic cholestasis of pregnancy (ICP), while being beneficial for 13 other diseases. The most common diseases related to TG_by_PG fell into the categories of endocrine / metabolic diseases and genitourinary system (Fig. 4, Table 3). Interestingly, lowering XS_VLDL_TG also had detrimental effects on the risk of ICP, suggesting that reducing TG_by_PG and XS_VLDL_TG as intervention targets for GDM may increase the risk of ICP. Raising PSM(d18:1/16:0) demonstrated deleterious effects on enthesopathies of limb (excluding foot) and osteoporosis with pathological fracture, but was beneficial for perforation of the tympanic membrane (Fig. S2, Table S11).Fig. 4 Potential on-target side effects associated with TG_by_PG, TG and PG intervention revealed by phenome–wide Mendelian randomization analysis. Odds ratios with their 95% confidence intervals represent the effect estimates on the risk of various non-GDM diseases of per 10% reduction in GDM risk by targeting a TG_by_PG, b triglycerides and c phosphoglycerides, respectively. Associations that are above the horizontal line with a black dash indicate deleterious side effects, while those below represent beneficial side effects. The upper red dashed horizontal line (odds ratio = 1.10) represents the point at which decreased GDM risk is counterbalanced by an equal increase in disease risk. CI   confidence interval, OR  odds ratio, TG_by_PG ratio of triglycerides to phosphoglycerides

Table 3 Descriptive summary of significant phenome-wide Mendelian randomization findings representing on-target side effects of biomarker intervention

Exposures	Outcome	IVW	Cochran’s Q test	MR-Egger	
OR (95% CI)	P-value	Q	Q_pval	intercept	P-value	
TG_by_PG	Disorders of lipoprotein metabolism	0.980 (0.967–0.993)	0.003	5.464	0.792	0.001	0.791	
	Metabolic disorders	0.980 (0.968–0.993)	0.003	6.553	0.767	0.003	0.572	
	Hyperlipidaemia, other/unspecified	0.975 (0.958–0.993)	0.006	7.568	0.911	-0.004	0.553	
	Familial hypercholesterolemia, with ischemic heart disease	0.957 (0.930–0.984)	0.002	15.038	0.920	0.008	0.495	
	Migraine	0.991 (0.984–0.997)	0.006	20.735	0.707	0.002	0.421	
	Psoriasis	0.984 (0.972–0.996)	0.010	25.414	0.551	0.005	0.395	
	Endometriosis of ovary	0.977 (0.962–0.992)	0.003	18.821	0.926	0.006	0.365	
	Endometriosis of pelvic peritoneum	0.973 (0.955–0.990)	0.002	26.359	0.606	0.002	0.776	
	Pilonidal cyst	0.966 (0.945–0.988)	0.003	22.325	0.807	0.005	0.604	
	Peripheral retinal degeneration	1.061 (1.014–1.109)	0.010	23.896	0.777	− 0.001	0.961	
	Asthma and allergy	0.976 (0.969–0.993)	0.005	22.455	0.868	− 0.001	0.897	
	Haemorrhage in early pregnancy	0.982 (0.971–0.994)	0.002	21.860	0.911	0.004	0.433	
	Intrahepatic Cholestasis of Pregnancy	1.040 (1.015–1.065)	0.001	17.679	0.913	− 0.004	0.655	
	Vestibular neuronitis	0.969 (0.947–0.992)	0.009	22.956	0.851	0.007	0.472	
	Inflammatory diseases of female pelvic organs	0.989 (0.981–0.996)	0.004	33.116	0.462	0.002	0.528	
Total triglycerides	Intrahepatic Cholestasis of Pregnancy*	1.607 (1.048–2.462)	0.029	14.097	0.442	0.021	0.177	
	Carcinoma in situ of cervix uteri	0.519 (0.332–0.812)	0.004	13.652	0.552	− 0.020	0.213	
	Familial hypercholesterolemia	0.155 (0.048–0.492)	0.002	0.804	0.938	0.006	0.834	
	Familial hypercholesterolemia, with ischemic heart disease	0.388 (0.201–0.752)	0.005	4.035	0.854	0.026	0.271	
	Calcific aortic valvular stenosis, operated	0.602 (0.429–0.842)	0.003	3.664	0.961	0.007	0.552	
	Coronary artery bypass grafting	0.518 (0.325–0.822)	0.006	5.215	0.516	0.025	0.142	
	Non-rheumatic valve diseases	0.726 (0.576–0.912)	0.007	9.753	0.638	− 0.006	0.377	
	Nasal polyp	1.583 (1.219–2.052)	0.001	7.531	0.912	0.001	0.954	
	Asthma and allergy	0.638 (0.474–0.862)	0.003	11.543	0.827	0.004	0.734	
	Impacted teeth	0.677 (0.507–0.902)	0.008	7.562	0.961	0.010	0.331	
	TMD muscular pain linked with fibromyalgia	0.755 (0.612–0.932)	0.009	12.357	0.778	0.002	0.797	
	Endometriosis of pelvic peritoneum	0.585 (0.421–0.812)	0.001	11.732	0.699	0.007	0.586	
	Endometriosis ASRM stages 1,2	0.672 (0.512–0.882)	0.004	9.029	0.876	0.007	0.496	
Total phosphoglycerides	Intrahepatic Cholestasis of Pregnancy*	0.945 (0.528–1.691)	0.848	9.796	0.200	− 0.104	0.116	
	Cerebral cysts	2.170 (1.270–3.708)	0.005	3.168	0.957	0.023	0.630	
	Peritonsillar abscess	1.518 (1.137–2.027)	0.005	4.416	0.731	− 0.005	0.885	
	Endometriosis ASRM stages 1,2	1.597 (1.196–2.132)	0.002	5.070	0.750	− 0.004	0.862	
	Glomerular disorders	0.569 (0.390–0.830)	0.003	9.195	0.419	− 0.012	0.739	
CI confidence interval, IVW inverse variance weighted, OR odds ratio, TG_by_PG ratio of triglycerides to phosphoglycerides

*P > 0.01

Since the indicator TG_by_PG is composed of two lipids–triglyceride (TG) and phosphoglycerides (PG)–adjusting TG_by_PG involves modifying the levels of these two lipids. Therefore, we conducted Phe-MR analysis on TG and PG separately. Lowering TG showed a detrimental effect on the risk of nasal polyps and beneficial effect on 11 other diseases. Raising PG had detrimental effects on the risk of cerebral cysts, peritonsillar abscess and endometriosis (ASRM stages 1,2), while demonstrating a beneficial effect on glomerular disorders. Lowering TG also showed detrimental effect on the risk of ICP, although this association did not reach the significance threshold of 0.01 (P = 0.029). No significant association was observed between PG levels and ICP (P = 0.848) (Fig. 4, Table 3).

Discussion

Among the 26 tier 1 lipid metabolites included, we identified TG_by_PG as a potential causal biomarker for GDM (IVW OR = 2.147, 95% CI: 1.415–3.257, P = 3.26e-4), and this association was confirmed by validation and reverse MR results. Additionally, two other lipid metabolites, PSM (d18:1/16:0) (IVW OR = 0.747, 95% CI: 0.583–0.956, P = 0.021) and XS_VLDL_TG (IVW OR = 2.948, 95% CI: 1.197–5.215, P = 0.015), were identified as suggestive potential causal biomarkers for GDM. Furthermore, Phe-MR analysis was used to assess potential on-target side effects related to GDM prevention via reducing TG_by_PG and XS_VLDL_TG levels as well as raising PSM (d18:1/16:0) levels. Beyond GDM, lowering TG_by_PG had detrimental effects on two diseases and beneficial effects on 13 diseases. Our findings highlight TG_by_PG as a lipid metabolite of interest in the pathogenesis of GDM, while the associations for XS_VLDL_TG and PSM(d18:1/16:0) are suggestive and warrant further investigation to confirm their potential roles.

Our study is the first to indicate a positive correlation between TG_by_PG level and GDM risk. However, previous studies have linked TG_by_PG to overweight/obese BMI [35] and immune response [36], and both overweight/obesity [37] and immune dysregulation [38] are associated with GDM, suggesting that TG_by_PG may mediate the development of GDM. Furthermore, Sequeiros et al. found that maternal TG_by_PG levels in the second trimester and in cord blood positively correlated with offspring ascending growth profiles (which is associated with higher risk of childhood adiposity) in the first two years post-delivery [39], highlighting TG_by_PG as a potential target for GDM intervention.

The association of TG_by_PG with GDM may be independently linked to its constituent lipid indices, triglycerides (TGs) and phosphoglycerides (PGs). TGs is the primary energy storage in mammals [40, 41], while PGs or GPs, which are glycerol-based phospholipids, constitute the main structural components of biological membranes in eukaryotic cells [42] and serve as sources of signaling molecules [43]. Observational studies have independently associated both TGs and PGs with GDM. Elevated TG levels have been linked to an increased risk of GDM [16, 44–46], whereas reduced PG levels are correlated with higher risk of GDM across different ethnic groups [16, 47–49]. The links between TGs or PGs with GDM can be explained by a higher degree of insulin resistance, which is not adequately compensated by increased β-cell secretion [44]. Evidence showed that PGs with ester-linkages and ether-linkages have protective roles against oxidative damage [50, 51], and it is postulated that oxidative stress, secondary to TG-related dyslipidemia, may lead to pancreatic β-cells dysfunction and suppression of insulin gene expression [52]. Elevated TG levels or a deficiency in PGs with ester-linkages and ether-linkages could reduce antioxidant protection, making pancreatic β cells more susceptible to oxidative stress [50, 51]. Given the low antioxidant capacity of β cells [53], increased oxidative damage may lead to β cell dysfunction, impairing insulin secretion and exacerbating physical insulin resistance during pregnancy, both key factors in the pathogenesis of GDM [51].

TG and PG are both important lipid classes influence GDM and are interconnected, as their synthesis pathways share common precursors such as phosphatidic acid (PA) and diacylglycerol (DAG) [54]. Therefore, the synthesis of PG might compete with TG synthesis for these shared intermediates. And targeting key enzymes involved in TG synthesis, such as phosphatidate phosphatase (PAP) [54], could be an effect strategy to reduce TG_by_PG level. Lifestyle modifications and several pharmacological agents, including fibrates, niacin, and long-chain omega-3 fatty acids, are known to reduce TG levels [55]. However, there are currently no clinical drugs or interventions specifically targeting PG, nor do these existing treatments address the shared synthesis pathway of TG and PG. Thus, TG appears to be a more practical intervention target at present. Nonetheless, Phe-MR results suggested that lowing TG_by_PG levels might increase the risk of ICP, and this side effect is mainly caused by reduced TG levels. Which indicated that TG_by_PG may be a potential target for preventing GDM, however, in pregnant women at high risk of ICP, a strategy that increases PG levels rather than reducing TG levels may be more appropriate.

Our results indicated that higher levels of TGs in very small very-low-density lipoprotein (VLDL) and lower levels of PSM (d18:1/16:0) is suggestively associated with an increased risk of GDM. A study in a Western European population found an association between XS_VLDL_TG and 2-h OGTT results (OR = 2.319, 95% CI: 1.103–4.874, Wu-Hausman P-value = 0.042), [56], aligning with our findings. However, although some studies have reported the association between certain sphingomyelin (SM) species and GDM [16, 57–59], PSM has rarely been reported [60], and our study is the first to identify the association between PSM (d18:1/16:0) and GDM.

Although further research is necessary to confirm these suggestive findings, they are supported by relevant biological mechanisms. The liver constantly synthesizes TGs by utilizing free fatty acids and carbohydrates, secreting them into circulation in the core of VLDL [61]. However, we did not find associations between other TG-rich lipoproteins (TRLs) (e.g., TG in small high-density lipoprotein) and GDM, as reported by other observational lipidomics studies [16, 46]. This may be explained by the relationship between the size of VLDL particles and their susceptibility to lipid oxidation, with smaller dense VLDL being more susceptible to oxidative modification [62]. The connection between PSM (d18:1/16:0) and GDM may be explained by sphingomyelin (SM) metabolism pathway, which is the most abundant SLs in mammalian tissues [63] and have been experimentally linked to pancreatic β-cell failure, insulin resistance, and T2DM [64].

To our knowledge, this is the first systematic MR study utilizing lipid metabolites as exposures to assess their causal effects on GDM risk, highlighting several strengths of our study. Firstly, this two-sample MR study was conducted based on three large-scale GWASs. The large sample size provided sufficient statistical power for making valid causal inferences. Secondly, our results were robust by means of stringent IVs selection criteria and a series of sensitivity analyses. The results were also confirmed by validation analysis and reverse MR analysis. Thirdly, we screened potential intervention targets using Phe-MR analysis in order to comprehensively forecast the on-target side effects of identified lipid metabolites.

However, there are several limitations to note. Firstly, our outcomes are pregnancy-specific and limited to women of childbearing age, while the exposure data come from the general population. Although all three GWASs of lipid metabolites were adjusted for sex and age, this study may not have identified sex- or age-specific instrumental variables (IVs) or IVs specifically related to lipid metabolites during pregnancy. Consequently, pregnancy lipidome GWAS data from pregnant women are needed for verification. Secondly, although the stringent IV selection criteria, the number of IVs available for each lipid metabolite was limited. While this approach helps reduce pleiotropy effects and guaranteed the robustness of our causal estimates, it may also inadvertently remove SNPs that influence multiple lipid metabolites through shared biological pathways, thereby potentially diminishing our power to detect causal relationships between lipid metabolites and GDM. Thirdly, the GWAS data used in this MR study came only from European populations (mostly whites), which may have limited the generalizability of our findings when applied to non-European populations and other racial groups. Nevertheless, this restriction minimized population and race stratification bias, though further studies are needed to confirm our findings in populations with different ethnic backgrounds. Moreover, due to strict selection criteria, only 295 of 1,023 lipid metabolites from three large GWAS were included in this study, which represent only a small fraction of the blood lipidome. Further research is required to investigate the associations between additional lipid metabolites and GDM. Lastly, not all the cohorts included in this study were measured using the exactly same blood sample collection methods, such as serum or plasma, fasting or non-fasting blood, which may introduce bias across studies. To validate our findings, further research using larger GWAS with fasting blood samples or individual-level data which allows stratification analysis is necessary.

The public health and clinical implications of our study lies in below: by integrating lipidomics with genomics, we provided novel insights that may contribute to the search for promising and safe intervention targets for GDM. Our findings showed a potential causal association between TG_by_PG and GDM, indicating that TG_by_PG could serve as a target for preventive interventions, and in pregnant women at high risk of ICP, a more nuanced approach may involve increasing PG levels rather than lowering TG levels to minimize adverse outcomes. Additionally, XS_VLDL_TG and PSM(d18:1/16:0) were found to be suggestively associated with GDM though further research is necessary to confirm the results. Clinical trials will be essential to confirm the feasibility and safety of these lipid metabolites in the prevention of GDM. Validating these findings will enhance the precision of prevention strategies for GDM.

Conclusions

This systematic MR study identified a potential causal association between TG_by_PG and GDM, as well as suggestive significant associations between XS_VLDL_TG, PSM (d18:1/16:0) and GDM. Side-effect profiles were characterized to inform drug target prioritization, suggesting that TG_by_PG might be a potential target for the prevention and treatment of GDM in caution of the increased risk of ICP.

Supplementary Information

Supplementary Material 1.

Supplementary Material 2.

Abbreviations

BHS The Busselton Health Study

Cer Ceramide

CI Confidence interval

cML-MA Constrained maximum likelihood and model averaging

EPIC The EPIC-Norfolk study

FA Fatty acyl

GDM Gestational diabetes mellitus

GL Glycerolipid

GP Glycerophospholipid

GWAS Genome-wide association study

HDL High-density lipoprotein

ICP Intrahepatic cholestasis of pregnancy

InSIDE Instrument strength independent of direct effect

INTERVAL The INTERVAL study

IV Instrument variable

IVW Inverse variance-weighted

LDL Low-density lipoprotein

MR Mendelian randomization

OR Odds ratio

PG Phosphoglyceride

Phe-MR Phenome-wide MR

PK Polyketide

PR Prenol lipid

PSM(d18:1/16:0) Palmitoyl sphingomyelin (d18:1/16:0)

SL Saccharolipid

SM Sphingomyelin

SMases Sphingomyelinases

SMSs Sphingomyelin synthases

SNP Single nucleotide polymorphism

SP Sphingolipid

STROBE-MR Strengthening the reporting of observational studies in epidemiology using mendelian randomization

ST Sterol lipid

T2DM Type 2 diabetes mellitus

TG Triglycerides

TG_by_PG Ratio of triglycerides to phosphoglycerides

TRL TG-rich lipoprotein

UKB UK Biobank

VLDL Very-low-density lipoprotein

XS_VLDL_TG Triglycerides in very small very low-density lipoprotein

Acknowledgements

We thank all participants and investigators who contribute to the GWASs from which the summary data in our study were derived. We thank Dr. Qiang Xia from New York City Department of Health and Mental Hygiene for his assistance in reviewing the language and grammar. This research was supported by the Medical Science Data Center of Fudan University.

Author contributions

YD and AH performed statistical analyses, interpreted the results, and drafted the manuscript. BH, MC, HL, ZL and QL revised the manuscript. Study design and supervision were conducted by YZ. All authors read and approved the final manuscript.

Funding

This study was funded by the National Natural Science Foundation of China [Grant No.82173582].

Availability of data and materials

All data used in the present study are described in the “Methods” section and Table S1, and each of the original GWAS publications. The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Declarations

Ethics approval and consent to participants

This study only used publicly available summary statistics from relevant GWASs and FinnGen consortium, thus no ethics approval is required. Respective ethics approvals have been obtained by the GWAS investigator from all participating studies and the FinnGen consortium investigators.

Consent of publication

Not required.

Competing interests

The authors declare no competing interests.

Publisher's Note

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

Yao Dong and An-qun Hu have contributed equally to this work.
==== Refs
References

1. Shah NS Wang MC Freaney PM Perak AM Carnethon MR Kandula NR Trends in gestational diabetes at first live birth by race and ethnicity in the US, 2011–2019 JAMA 2021 326 7 660 669 10.1001/jama.2021.7217 34402831
Shah NS, Wang MC, Freaney PM, Perak AM, Carnethon MR, Kandula NR, et al. Trends in gestational diabetes at first live birth by race and ethnicity in the US, 2011–2019. JAMA. 2021;326(7):660–9.34402831 10.1001/jama.2021.7217
2. Ye W Luo C Huang J Li C Liu Z Liu F Gestational diabetes mellitus and adverse pregnancy outcomes: systematic review and meta-analysis BMJ 2022 377 e067946 10.1136/bmj-2021-067946 35613728
Ye W, Luo C, Huang J, Li C, Liu Z, Liu F. Gestational diabetes mellitus and adverse pregnancy outcomes: systematic review and meta-analysis. BMJ. 2022;377: e067946.35613728 10.1136/bmj-2021-067946
3. Group HSCR The hyperglycemia and adverse pregnancy outcome (HAPO) study Int J Gynaecol Obstet 2002 78 1 69 77 10.1016/S0020-7292(02)00092-9 12113977
Group HSCR. The hyperglycemia and adverse pregnancy outcome (HAPO) study. Int J Gynaecol Obstet. 2002;78(1):69–77.12113977 10.1016/S0020-7292(02)00092-9
4. Landon MB Spong CY Thom E Carpenter MW Ramin SM Casey B A multicenter, randomized trial of treatment for mild gestational diabetes N Engl J Med 2009 361 14 1339 1348 10.1056/NEJMoa0902430 19797280
Landon MB, Spong CY, Thom E, Carpenter MW, Ramin SM, Casey B, et al. A multicenter, randomized trial of treatment for mild gestational diabetes. N Engl J Med. 2009;361(14):1339–48.19797280 10.1056/NEJMoa0902430
5. Casey BM Rice MM Landon MB Varner MW Reddy UM Wapner RJ Effect of treatment of mild gestational diabetes on long-term maternal outcomes Am J Perinatol 2020 37 5 475 482 10.1055/s-0039-1681058 30866027
Casey BM, Rice MM, Landon MB, Varner MW, Reddy UM, Wapner RJ, et al. Effect of treatment of mild gestational diabetes on long-term maternal outcomes. Am J Perinatol. 2020;37(5):475–82.30866027 10.1055/s-0039-1681058
6. Catalano PM Huston L Amini SB Kalhan SC Longitudinal changes in glucose metabolism during pregnancy in obese women with normal glucose tolerance and gestational diabetes mellitus Am J Obstet Gynecol 1999 180 4 903 916 10.1016/S0002-9378(99)70662-9 10203659
Catalano PM, Huston L, Amini SB, Kalhan SC. Longitudinal changes in glucose metabolism during pregnancy in obese women with normal glucose tolerance and gestational diabetes mellitus. Am J Obstet Gynecol. 1999;180(4):903–16.10203659 10.1016/S0002-9378(99)70662-9
7. Wiznitzer A Mayer A Novack V Sheiner E Gilutz H Malhotra A Association of lipid levels during gestation with preeclampsia and gestational diabetes mellitus: a population-based study Am J Obstet Gynecol 2009 201 5 482.e481 488 10.1016/j.ajog.2009.05.032
Wiznitzer A, Mayer A, Novack V, Sheiner E, Gilutz H, Malhotra A, et al. Association of lipid levels during gestation with preeclampsia and gestational diabetes mellitus: a population-based study. Am J Obstet Gynecol. 2009;201(5):482.e481-488.10.1016/j.ajog.2009.05.032
8. Catalano PM Trying to understand gestational diabetes Diabet Med 2014 31 3 273 281 10.1111/dme.12381 24341419
Catalano PM. Trying to understand gestational diabetes. Diabet Med. 2014;31(3):273–81.24341419 10.1111/dme.12381
9. Plows JF Stanley JL Baker PN Reynolds CM Vickers MH The pathophysiology of gestational diabetes mellitus Int J Mol Sci 2018 19 11 3342 10.3390/ijms19113342 30373146
Plows JF, Stanley JL, Baker PN, Reynolds CM, Vickers MH. The pathophysiology of gestational diabetes mellitus. Int J Mol Sci. 2018;19(11):3342.30373146 10.3390/ijms19113342
10. Wang Y Pan XF Pan A Lipidomics in gestational diabetes mellitus Curr Opin Lipidol 2023 34 1 1 11 10.1097/MOL.0000000000000858 36637074
Wang Y, Pan XF, Pan A. Lipidomics in gestational diabetes mellitus. Curr Opin Lipidol. 2023;34(1):1–11.36637074 10.1097/MOL.0000000000000858
11. Liu J Li J Li S Leng J Li W Yang W Circulating lysophosphatidylcholines in early pregnancy and risk of gestational diabetes in chinese women J Clin Endocrinol Metab 2020 10.1210/clinem/dgaa058 32979271
Liu J, Li J, Li S, Leng J, Li W, Yang W, et al. Circulating lysophosphatidylcholines in early pregnancy and risk of gestational diabetes in chinese women. J Clin Endocrinol Metab. 2020. 10.1210/clinem/dgaa058.32979271 10.1210/clinem/dgaa058
12. Liu J Li J Li W Li N Huo X Wang H Predictive values of serum metabolites in early pregnancy and their possible pathways for gestational diabetes: a nested case-control study in Tianjin, China J Diabetes Complications 2021 35 12 108048 10.1016/j.jdiacomp.2021.108048 34563440
Liu J, Li J, Li W, Li N, Huo X, Wang H, et al. Predictive values of serum metabolites in early pregnancy and their possible pathways for gestational diabetes: a nested case-control study in Tianjin, China. J Diabetes Complications. 2021;35(12): 108048.34563440 10.1016/j.jdiacomp.2021.108048
13. Han BX Yan SS Yu H Xu Q Zhao QG Ma XL Causal effects of plasma proteome on osteoporosis and osteoarthritis Calcif Tissue Int 2023 112 3 350 358 10.1007/s00223-022-01049-w 36576504
Han BX, Yan SS, Yu H, Xu Q, Zhao QG, Ma XL, et al. Causal effects of plasma proteome on osteoporosis and osteoarthritis. Calcif Tissue Int. 2023;112(3):350–8.36576504 10.1007/s00223-022-01049-w
14. Zhu Y Tsai MY Sun Q Hinkle SN Rawal S Mendola P A prospective and longitudinal study of plasma phospholipid saturated fatty acid profile in relation to cardiometabolic biomarkers and the risk of gestational diabetes Am J Clin Nutr 2018 107 6 1017 1026 10.1093/ajcn/nqy051 29868913
Zhu Y, Tsai MY, Sun Q, Hinkle SN, Rawal S, Mendola P, et al. A prospective and longitudinal study of plasma phospholipid saturated fatty acid profile in relation to cardiometabolic biomarkers and the risk of gestational diabetes. Am J Clin Nutr. 2018;107(6):1017–26.29868913 10.1093/ajcn/nqy051
15. Zhu Y Li M Rahman ML Hinkle SN Wu J Weir NL Plasma phospholipid n-3 and n-6 polyunsaturated fatty acids in relation to cardiometabolic markers and gestational diabetes: A longitudinal study within the prospective NICHD Fetal Growth Studies PLoS Med 2019 16 9 e1002910 10.1371/journal.pmed.1002910 31518348
Zhu Y, Li M, Rahman ML, Hinkle SN, Wu J, Weir NL, et al. Plasma phospholipid n-3 and n-6 polyunsaturated fatty acids in relation to cardiometabolic markers and gestational diabetes: A longitudinal study within the prospective NICHD Fetal Growth Studies. PLoS Med. 2019;16(9): e1002910.31518348 10.1371/journal.pmed.1002910
16. Rahman ML Feng YA Fiehn O Albert PS Tsai MY Zhu Y Plasma lipidomics profile in pregnancy and gestational diabetes risk: a prospective study in a multiracial/ethnic cohort BMJ Open Diabetes Res Care 2021 9 1 e001551 10.1136/bmjdrc-2020-001551 33674279
Rahman ML, Feng YA, Fiehn O, Albert PS, Tsai MY, Zhu Y, et al. Plasma lipidomics profile in pregnancy and gestational diabetes risk: a prospective study in a multiracial/ethnic cohort. BMJ Open Diabetes Res Care. 2021;9(1):e001551.33674279 10.1136/bmjdrc-2020-001551
17. Parhofer KG Interaction between glucose and lipid metabolism: more than diabetic dyslipidemia Diabetes Metab J 2015 39 5 353 362 10.4093/dmj.2015.39.5.353 26566492
Parhofer KG. Interaction between glucose and lipid metabolism: more than diabetic dyslipidemia. Diabetes Metab J. 2015;39(5):353–62.26566492 10.4093/dmj.2015.39.5.353
18. Surendran P Stewart ID Au Yeung VPW Pietzner M Raffler J Wörheide MA Rare and common genetic determinants of metabolic individuality and their effects on human health Nat Med 2022 28 11 2321 2332 10.1038/s41591-022-02046-0 36357675
Surendran P, Stewart ID, Au Yeung VPW, Pietzner M, Raffler J, Wörheide MA, et al. Rare and common genetic determinants of metabolic individuality and their effects on human health. Nat Med. 2022;28(11):2321–32.36357675 10.1038/s41591-022-02046-0
19. Cadby G Giles C Melton PE Huynh K Mellett NA Duong T Comprehensive genetic analysis of the human lipidome identifies loci associated with lipid homeostasis with links to coronary artery disease Nat Commun 2022 13 1 3124 10.1038/s41467-022-30875-7 35668104
Cadby G, Giles C, Melton PE, Huynh K, Mellett NA, Duong T, et al. Comprehensive genetic analysis of the human lipidome identifies loci associated with lipid homeostasis with links to coronary artery disease. Nat Commun. 2022;13(1):3124.35668104 10.1038/s41467-022-30875-7
20. Karjalainen MK Karthikeyan S Oliver-Williams C Sliz E Allara E Fung WT Genome-wide characterization of circulating metabolic biomarkers Nature 2024 628 8006 130 138 10.1038/s41586-024-07148-y 38448586
Karjalainen MK, Karthikeyan S, Oliver-Williams C, Sliz E, Allara E, Fung WT, et al. Genome-wide characterization of circulating metabolic biomarkers. Nature. 2024;628(8006):130–8.38448586 10.1038/s41586-024-07148-y
21. Kettunen J Demirkan A Würtz P Draisma HH Haller T Rawal R Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA Nat Commun 2016 7 11122 10.1038/ncomms11122 27005778
Kettunen J, Demirkan A, Würtz P, Draisma HH, Haller T, Rawal R, et al. Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA. Nat Commun. 2016;7:11122.27005778 10.1038/ncomms11122
22. Kurki MI Karjalainen J Palta P Sipilä TP Kristiansson K Donner KM FinnGen provides genetic insights from a well-phenotyped isolated population Nature 2023 613 7944 508 518 10.1038/s41586-022-05473-8 36653562
Kurki MI, Karjalainen J, Palta P, Sipilä TP, Kristiansson K, Donner KM, et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature. 2023;613(7944):508–18.36653562 10.1038/s41586-022-05473-8
23. Hemani G Zheng J Elsworth B Wade KH Haberland V Baird D The MR-Base platform supports systematic causal inference across the human phenome Elife 2018 10.7554/eLife.34408 29846171
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018. 10.7554/eLife.34408.29846171 10.7554/eLife.34408
24. Burgess S Dudbridge F Thompson SG Combining information on multiple instrumental variables in Mendelian randomization: comparison of allele score and summarized data methods Stat Med 2016 35 11 1880 1906 10.1002/sim.6835 26661904
Burgess S, Dudbridge F, Thompson SG. Combining information on multiple instrumental variables in Mendelian randomization: comparison of allele score and summarized data methods. Stat Med. 2016;35(11):1880–906.26661904 10.1002/sim.6835
25. Bowden J Davey Smith G Burgess S Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression Int J Epidemiol 2015 44 2 512 525 10.1093/ije/dyv080 26050253
Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25.26050253 10.1093/ije/dyv080
26. Hartwig FP Davey Smith G Bowden J Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption Int J Epidemiol 2017 46 6 1985 1998 10.1093/ije/dyx102 29040600
Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. 2017;46(6):1985–98.29040600 10.1093/ije/dyx102
27. Burgess S Thompson SG Interpreting findings from Mendelian randomization using the MR-Egger method Eur J Epidemiol 2017 32 5 377 389 10.1007/s10654-017-0255-x 28527048
Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32(5):377–89.28527048 10.1007/s10654-017-0255-x
28. Verbanck M Chen CY Neale B Do R Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases Nat Genet 2018 50 5 693 698 10.1038/s41588-018-0099-7 29686387
Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693–8.29686387 10.1038/s41588-018-0099-7
29. Bowden J Davey Smith G Haycock PC Burgess S Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator Genet Epidemiol 2016 40 4 304 314 10.1002/gepi.21965 27061298
Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14.27061298 10.1002/gepi.21965
30. Xue H Shen X Pan W Constrained maximum likelihood-based Mendelian randomization robust to both correlated and uncorrelated pleiotropic effects Am J Hum Genet 2021 108 7 1251 1269 10.1016/j.ajhg.2021.05.014 34214446
Xue H, Shen X, Pan W. Constrained maximum likelihood-based Mendelian randomization robust to both correlated and uncorrelated pleiotropic effects. Am J Hum Genet. 2021;108(7):1251–69.34214446 10.1016/j.ajhg.2021.05.014
31. Bowden J Del Greco MF Minelli C Davey Smith G Sheehan N Thompson J A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization Stat Med 2017 36 11 1783 1802 10.1002/sim.7221 28114746
Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N, Thompson J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med. 2017;36(11):1783–802.28114746 10.1002/sim.7221
32. Cheng H Garrick DJ Fernando RL Efficient strategies for leave-one-out cross validation for genomic best linear unbiased prediction J Anim Sci Biotechnol 2017 8 38 10.1186/s40104-017-0164-6 28469846
Cheng H, Garrick DJ, Fernando RL. Efficient strategies for leave-one-out cross validation for genomic best linear unbiased prediction. J Anim Sci Biotechnol. 2017;8:38.28469846 10.1186/s40104-017-0164-6
33. Curtin F Schulz P Multiple correlations and Bonferroni's correction Biol Psychiatry 1998 44 8 775 777 10.1016/S0006-3223(98)00043-2 9798082
Curtin F, Schulz P. Multiple correlations and Bonferroni’s correction. Biol Psychiatry. 1998;44(8):775–7.9798082 10.1016/S0006-3223(98)00043-2
34. Yavorska OO Burgess S MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data Int J Epidemiol 2017 46 6 1734 1739 10.1093/ije/dyx034 28398548
Yavorska OO, Burgess S. MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data. Int J Epidemiol. 2017;46(6):1734–9.28398548 10.1093/ije/dyx034
35. Mensah GA Fuster V Murray CJL Roth GA Global Burden of Cardiovascular Diseases and Risks, 1990–2022 J Am Coll Cardiol 2023 82 25 2350 2473 10.1016/j.jacc.2023.11.007 38092509
Mensah GA, Fuster V, Murray CJL, Roth GA. Global Burden of Cardiovascular Diseases and Risks, 1990–2022. J Am Coll Cardiol. 2023;82(25):2350–473.38092509 10.1016/j.jacc.2023.11.007
36. Inouye M Kettunen J Soininen P Silander K Ripatti S Kumpula LS Metabonomic, transcriptomic, and genomic variation of a population cohort Mol Syst Biol 2010 6 441 10.1038/msb.2010.93 21179014
Inouye M, Kettunen J, Soininen P, Silander K, Ripatti S, Kumpula LS, et al. Metabonomic, transcriptomic, and genomic variation of a population cohort. Mol Syst Biol. 2010;6:441.21179014 10.1038/msb.2010.93
37. Sweeting A Hannah W Backman H Catalano P Feghali M Herman WH Epidemiology and management of gestational diabetes Lancet 2024 10.1016/S0140-6736(24)00825-0 38909620
Sweeting A, Hannah W, Backman H, Catalano P, Feghali M, Herman WH, et al. Epidemiology and management of gestational diabetes. Lancet. 2024. 10.1016/S0140-6736(24)00825-0.38909620 10.1016/S0140-6736(24)00825-0
38. McElwain CJ McCarthy FP McCarthy CM Gestational diabetes mellitus and maternal immune dysregulation: what we know so far Int J Mol Sci 2021 22 8 4261 10.3390/ijms22084261 33923959
McElwain CJ, McCarthy FP, McCarthy CM. Gestational diabetes mellitus and maternal immune dysregulation: what we know so far. Int J Mol Sci. 2021;22(8):4261.33923959 10.3390/ijms22084261
39. Blanco Sequeiros E Tuomaala AK Tabassum R Bergman PH Koivusalo SB Huvinen E Early ascending growth is associated with maternal lipoprotein profile during mid and late pregnancy and in cord blood Int J Obes (Lond) 2023 47 11 1081 1087 10.1038/s41366-023-01361-x 37592059
Blanco Sequeiros E, Tuomaala AK, Tabassum R, Bergman PH, Koivusalo SB, Huvinen E. Early ascending growth is associated with maternal lipoprotein profile during mid and late pregnancy and in cord blood. Int J Obes (Lond). 2023;47(11):1081–7.37592059 10.1038/s41366-023-01361-x
40. Coleman RA Lee DP Enzymes of triacylglycerol synthesis and their regulation Prog Lipid Res 2004 43 2 134 176 10.1016/S0163-7827(03)00051-1 14654091
Coleman RA, Lee DP. Enzymes of triacylglycerol synthesis and their regulation. Prog Lipid Res. 2004;43(2):134–76.14654091 10.1016/S0163-7827(03)00051-1
41. McLelland GL Lopez-Osias M Verzijl CRC Ellenbroek BD Oliveira RA Boon NJ Identification of an alternative triglyceride biosynthesis pathway Nature 2023 621 7977 171 178 10.1038/s41586-023-06497-4 37648867
McLelland GL, Lopez-Osias M, Verzijl CRC, Ellenbroek BD, Oliveira RA, Boon NJ, et al. Identification of an alternative triglyceride biosynthesis pathway. Nature. 2023;621(7977):171–8.37648867 10.1038/s41586-023-06497-4
42. Osawa T Fujikawa K Shimamoto K Structures, functions, and syntheses of glycero-glycophospholipids Front Chem 2024 12 1353688 10.3389/fchem.2024.1353688 38389730
Osawa T, Fujikawa K, Shimamoto K. Structures, functions, and syntheses of glycero-glycophospholipids. Front Chem. 2024;12:1353688.38389730 10.3389/fchem.2024.1353688
43. Casares D Escribá PV Rosselló CA Membrane Lipid Composition: Effect on Membrane and Organelle Structure, Function and Compartmentalization and Therapeutic Avenues Int J Mol Sci 2019 10.3390/ijms20092167 31052427
Casares D, Escribá PV, Rosselló CA. Membrane Lipid Composition: Effect on Membrane and Organelle Structure, Function and Compartmentalization and Therapeutic Avenues. Int J Mol Sci. 2019. 10.3390/ijms20092167.31052427 10.3390/ijms20092167
44. Eppel D Feichtinger M Lindner T Kotzaeridi G Rosicky I Yerlikaya-Schatten G Association between maternal triglycerides and disturbed glucose metabolism in pregnancy Acta Diabetol 2021 58 4 459 465 10.1007/s00592-020-01644-z 33387029
Eppel D, Feichtinger M, Lindner T, Kotzaeridi G, Rosicky I, Yerlikaya-Schatten G, et al. Association between maternal triglycerides and disturbed glucose metabolism in pregnancy. Acta Diabetol. 2021;58(4):459–65.33387029 10.1007/s00592-020-01644-z
45. Bao W Dar S Zhu Y Wu J Rawal S Li S Plasma concentrations of lipids during pregnancy and the risk of gestational diabetes mellitus: a longitudinal study J Diabetes 2018 10 6 487 495 10.1111/1753-0407.12563 28436169
Bao W, Dar S, Zhu Y, Wu J, Rawal S, Li S, et al. Plasma concentrations of lipids during pregnancy and the risk of gestational diabetes mellitus: a longitudinal study. J Diabetes. 2018;10(6):487–95.28436169 10.1111/1753-0407.12563
46. Lu L Koulman A Petry CJ Jenkins B Matthews L Hughes IA An unbiased lipidomics approach identifies early second trimester lipids predictive of maternal glycemic traits and gestational diabetes mellitus Diabetes Care 2016 39 12 2232 2239 10.2337/dc16-0863 27703025
Lu L, Koulman A, Petry CJ, Jenkins B, Matthews L, Hughes IA, et al. An unbiased lipidomics approach identifies early second trimester lipids predictive of maternal glycemic traits and gestational diabetes mellitus. Diabetes Care. 2016;39(12):2232–9.27703025 10.2337/dc16-0863
47. Law KP Mao X Han TL Zhang H Unsaturated plasma phospholipids are consistently lower in the patients diagnosed with gestational diabetes mellitus throughout pregnancy: a longitudinal metabolomics study of Chinese pregnant women part 1 Clin Chim Acta 2017 465 53 71 10.1016/j.cca.2016.12.010 27988319
Law KP, Mao X, Han TL, Zhang H. Unsaturated plasma phospholipids are consistently lower in the patients diagnosed with gestational diabetes mellitus throughout pregnancy: a longitudinal metabolomics study of Chinese pregnant women part 1. Clin Chim Acta. 2017;465:53–71.27988319 10.1016/j.cca.2016.12.010
48. Dudzik D Zorawski M Skotnicki M Zarzycki W Kozlowska G Bibik-Malinowska K Metabolic fingerprint of gestational diabetes mellitus J Proteomics 2014 103 57 71 10.1016/j.jprot.2014.03.025 24698665
Dudzik D, Zorawski M, Skotnicki M, Zarzycki W, Kozlowska G, Bibik-Malinowska K, et al. Metabolic fingerprint of gestational diabetes mellitus. J Proteomics. 2014;103:57–71.24698665 10.1016/j.jprot.2014.03.025
49. Lehmann R Friedrich T Krebiehl G Sonntag D Häring HU Fritsche A Metabolic profiles during an oral glucose tolerance test in pregnant women with and without gestational diabetes Exp Clin Endocrinol Diabetes 2015 123 7 483 438 10.1055/s-0035-1549887 26171623
Lehmann R, Friedrich T, Krebiehl G, Sonntag D, Häring HU, Fritsche A, et al. Metabolic profiles during an oral glucose tolerance test in pregnant women with and without gestational diabetes. Exp Clin Endocrinol Diabetes. 2015;123(7):483–438.26171623 10.1055/s-0035-1549887
50. Dean JM Lodhi IJ Structural and functional roles of ether lipids Protein Cell 2018 9 2 196 206 10.1007/s13238-017-0423-5 28523433
Dean JM, Lodhi IJ. Structural and functional roles of ether lipids. Protein Cell. 2018;9(2):196–206.28523433 10.1007/s13238-017-0423-5
51. Razquin C Toledo E Clish CB Ruiz-Canela M Dennis C Corella D Plasma lipidomic profiling and risk of type 2 diabetes in the PREDIMED trial Diabetes Care 2018 41 12 2617 2624 10.2337/dc18-0840 30327364
Razquin C, Toledo E, Clish CB, Ruiz-Canela M, Dennis C, Corella D, et al. Plasma lipidomic profiling and risk of type 2 diabetes in the PREDIMED trial. Diabetes Care. 2018;41(12):2617–24.30327364 10.2337/dc18-0840
52. Eguchi N Vaziri ND Dafoe DC Ichii H The role of oxidative stress in pancreatic β cell dysfunction in diabetes Int J Mol Sci 2021 22 4 1509 10.3390/ijms22041509 33546200
Eguchi N, Vaziri ND, Dafoe DC, Ichii H. The role of oxidative stress in pancreatic β cell dysfunction in diabetes. Int J Mol Sci. 2021;22(4):1509.33546200 10.3390/ijms22041509
53. Lenzen S Drinkgern J Tiedge M Low antioxidant enzyme gene expression in pancreatic islets compared with various other mouse tissues Free Radic Biol Med 1996 20 3 463 466 10.1016/0891-5849(96)02051-5 8720919
Lenzen S, Drinkgern J, Tiedge M. Low antioxidant enzyme gene expression in pancreatic islets compared with various other mouse tissues. Free Radic Biol Med. 1996;20(3):463–6.8720919 10.1016/0891-5849(96)02051-5
54. Lutkewitte AJ Finck BN Regulation of signaling and metabolism by lipin-mediated phosphatidic acid phosphohydrolase activity Biomolecules 2020 10 10 1386 10.3390/biom10101386 33003344
Lutkewitte AJ, Finck BN. Regulation of signaling and metabolism by lipin-mediated phosphatidic acid phosphohydrolase activity. Biomolecules. 2020;10(10):1386.33003344 10.3390/biom10101386
55. Marston NA Giugliano RP Im K Silverman MG O'Donoghue ML Wiviott SD Association between triglyceride lowering and reduction of cardiovascular risk across multiple lipid-lowering therapeutic classes: a systematic review and meta-regression analysis of randomized controlled trials Circulation 2019 140 16 1308 1317 10.1161/CIRCULATIONAHA.119.041998 31530008
Marston NA, Giugliano RP, Im K, Silverman MG, O’Donoghue ML, Wiviott SD, et al. Association between triglyceride lowering and reduction of cardiovascular risk across multiple lipid-lowering therapeutic classes: a systematic review and meta-regression analysis of randomized controlled trials. Circulation. 2019;140(16):1308–17.31530008 10.1161/CIRCULATIONAHA.119.041998
56. Fuller H Iles MM Moore JB Zulyniak MA Metabolic drivers of dysglycemia in pregnancy: ethnic-specific GWAS of 146 metabolites and 1-sample Mendelian randomization analyses in a UK multi-ethnic birth cohort Front Endocrinol (Lausanne) 2023 14 1157416 10.3389/fendo.2023.1157416 37255970
Fuller H, Iles MM, Moore JB, Zulyniak MA. Metabolic drivers of dysglycemia in pregnancy: ethnic-specific GWAS of 146 metabolites and 1-sample Mendelian randomization analyses in a UK multi-ethnic birth cohort. Front Endocrinol (Lausanne). 2023;14:1157416.37255970 10.3389/fendo.2023.1157416
57. McMichael LE Heath H Johnson CM Fanter R Alarcon N Quintana-Diaz A Metabolites involved in purine degradation, insulin resistance, and fatty acid oxidation are associated with prediction of Gestational diabetes in plasma Metabolomics 2021 17 12 105 10.1007/s11306-021-01857-5 34837546
McMichael LE, Heath H, Johnson CM, Fanter R, Alarcon N, Quintana-Diaz A, et al. Metabolites involved in purine degradation, insulin resistance, and fatty acid oxidation are associated with prediction of Gestational diabetes in plasma. Metabolomics. 2021;17(12):105.34837546 10.1007/s11306-021-01857-5
58. Wang Y Huang Y Wu P Ye Y Sun F Yang X Plasma lipidomics in early pregnancy and risk of gestational diabetes mellitus: a prospective nested case-control study in Chinese women Am J Clin Nutr 2021 114 5 1763 1773 10.1093/ajcn/nqab242 34477820
Wang Y, Huang Y, Wu P, Ye Y, Sun F, Yang X, et al. Plasma lipidomics in early pregnancy and risk of gestational diabetes mellitus: a prospective nested case-control study in Chinese women. Am J Clin Nutr. 2021;114(5):1763–73.34477820 10.1093/ajcn/nqab242
59. Wu P Wang Y Ye Y Yang X Huang Y Ye Y Liver biomarkers, lipid metabolites, and risk of gestational diabetes mellitus in a prospective study among Chinese pregnant women BMC Med 2023 21 1 150 10.1186/s12916-023-02818-6 37069659
Wu P, Wang Y, Ye Y, Yang X, Huang Y, Ye Y, et al. Liver biomarkers, lipid metabolites, and risk of gestational diabetes mellitus in a prospective study among Chinese pregnant women. BMC Med. 2023;21(1):150.37069659 10.1186/s12916-023-02818-6
60. Chen Y Jia H Qian X Wang J Yu M Gong Q Circulating palmitoyl sphingomyelin is associated with cardiovascular disease in individuals with type 2 diabetes: findings from the china da qing diabetes study Diabetes Care 2022 45 3 666 673 10.2337/dc21-1520 35165706
Chen Y, Jia H, Qian X, Wang J, Yu M, Gong Q, et al. Circulating palmitoyl sphingomyelin is associated with cardiovascular disease in individuals with type 2 diabetes: findings from the china da qing diabetes study. Diabetes Care. 2022;45(3):666–73.35165706 10.2337/dc21-1520
61. Walker HK HW, Hurst JW: Clinical methods: the history, physical, and laboratory examinations. 3rd edition; 1990.
62. McEneny J McMaster C Trimble ER Young IS Rapid isolation of VLDL subfractions: assessment of composition and susceptibility to copper-mediated oxidation J Lipid Res 2002 43 5 824 831 10.1016/S0022-2275(20)30125-5 11971954
McEneny J, McMaster C, Trimble ER, Young IS. Rapid isolation of VLDL subfractions: assessment of composition and susceptibility to copper-mediated oxidation. J Lipid Res. 2002;43(5):824–31.11971954 10.1016/S0022-2275(20)30125-5
63. Bienias K Fiedorowicz A Sadowska A Prokopiuk S Car H Regulation of sphingomyelin metabolism Pharmacol Rep 2016 68 3 570 581 10.1016/j.pharep.2015.12.008 26940196
Bienias K, Fiedorowicz A, Sadowska A, Prokopiuk S, Car H. Regulation of sphingomyelin metabolism. Pharmacol Rep. 2016;68(3):570–81.26940196 10.1016/j.pharep.2015.12.008
64. Holland WL Summers SA Sphingolipids, insulin resistance, and metabolic disease: new insights from in vivo manipulation of sphingolipid metabolism Endocr Rev 2008 29 4 381 402 10.1210/er.2007-0025 18451260
Holland WL, Summers SA. Sphingolipids, insulin resistance, and metabolic disease: new insights from in vivo manipulation of sphingolipid metabolism. Endocr Rev. 2008;29(4):381–402.18451260 10.1210/er.2007-0025
