Next Article in Journal
Peripapillary Oxygenation and Retinal Vascular Responsiveness to Flicker Light in Primary Open Angle Glaucoma
Next Article in Special Issue
Whole Exome Sequencing Enhanced Imputation Identifies 85 Metabolite Associations in the Alpine CHRIS Cohort
Previous Article in Journal
PPARβ/δ Augments IL-1β-Induced COX-2 Expression and PGE2 Biosynthesis in Human Mesangial Cells via the Activation of SIRT1
Previous Article in Special Issue
Using Mendelian Randomisation to Prioritise Candidate Maternal Metabolic Traits Influencing Offspring Birthweight
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 12
Issue 7
10.3390/metabo12070596
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Architecture of Untargeted Lipidomics in Cardiometabolic-Disease Patients Combines Strong Polygenic Control and Pleiotropy

by
Francois Brial
1,2,
Lyamine Hedjazi
3,
Kazuhiro Sonomura
4,
Cynthia Al Hageh
5,
Pierre Zalloua
5,
Fumihiko Matsuda
1,6 and
Dominique Gauguier
1,2,6,*
1
Center for Genomic Medicine, Graduate School of Medicine Kyoto University, Kyoto 606-8501, Japan
2
INSERM UMR 1124, Université Paris Cité, 45 rue des Saint-Pères, 75006 Paris, France
3
Beemetrix SAS, 30 Avenue Carnot, 91300 Massy, France
4
Life Science Research Center, Technology Research Laboratory, Shimadzu Corporation, Kyoto 606-8501, Japan
5
College of Medicine and Health Sciences, Khalifa University, Abu Dhabi P.O. Box 17666, United Arab Emirates
6
McGill University and Genome Quebec Innovation Centre, 740 Doctor Penfield Avenue, Montreal, QC H3A 0G1, Canada
*
Author to whom correspondence should be addressed.
Metabolites 2022, 12(7), 596; https://doi.org/10.3390/metabo12070596
Submission received: 5 June 2022 / Revised: 20 June 2022 / Accepted: 23 June 2022 / Published: 27 June 2022
(This article belongs to the Special Issue The Intersection of Metabolomics and Genomics and Their Role in Human Health)
Browse Figures
Versions Notes

Abstract

:
Analysis of the genetic control of small metabolites provides powerful information on the regulation of the endpoints of genome expression. We carried out untargeted liquid chromatography–high-resolution mass spectrometry in 273 individuals characterized for pathophysiological elements of the cardiometabolic syndrome. We quantified 3013 serum lipidomic features, which we used in both genome-wide association studies (GWAS), using a panel of over 2.5 M imputed single-nucleotide polymorphisms (SNPs), and metabolome-wide association studies (MWAS) with phenotypes. Genetic analyses showed that 926 SNPs at 551 genetic loci significantly (q-value < 10−8) regulate the abundance of 74 lipidomic features in the group, with evidence of monogenic control for only 22 of these. In addition to this strong polygenic control of serum lipids, our results underscore instances of pleiotropy, when a single genetic locus controls the abundance of several distinct lipid features. Using the LIPID MAPS database, we assigned putative lipids, predominantly fatty acyls and sterol lipids, to 77% of the lipidome signals mapped to the genome. We identified significant correlations between lipids and clinical and biochemical phenotypes. These results demonstrate the power of untargeted lipidomic profiling for high-density quantitative molecular phenotyping in human-genetic studies and illustrate the complex genetic control of lipid metabolism.

Graphical Abstract

1. Introduction

Molecular-phenotyping tools based on transcriptome, proteome and metabolome technologies provide detailed information on the molecular pathways and biomarkers relevant to disease etiopathogenesis. Their application in the context of genome-wide association studies (GWAS) of complex disorders can enhance our understanding of the genetic control of genome expression and to dissect out disease variables into multiple, intermediate disease traits and molecular phenotypes [1,2]. Metabolomics, which analyses the multivariate data representing a range of small metabolites in a biological sample, has already been used in humans to map the genetic determinants of the quantitative variations of metabolites [3]. Owing to the role of altered plasma-lipid profiles in many chronic-disease manifestations, including chronic kidney disease, cardiovascular risk, dyslipidemia and neurological disorders, the detection and quantification of lipids in a biospecimen through lipidomics has emerged as a promising approach to correlate variations in blood lipids with these diseases [4,5,6].
Even though elevated blood LDL cholesterol is known to be a major risk factor for coronary heart disease and stroke, lipidomics enables a hypothesis-free strategy for broadening the search for the biomarkers associated with these diseases to a wide range of lipid species and to uncover novel targets beyond traditional lipids that can predict or reduce the risk of cardiovascular diseases [7,8]. Among examples of lipid classes that can be detected and quantified through lipidomic technologies, ceramides are involved in vascular inflammation and apoptosis and may have a higher potential to predict coronary heart disease than LDL cholesterol [9]. Ceramides, but more prominently the phospholipid species, alter the progression to ischemic cardiomyopathy [10]). Beyond associations between lipids and disease, combining genetics and lipidomics allows the identification of the genetic factors involved in the coordinated regulation of lipid species, thus inferring functional connections between different lipid species and causal relationships between lipid species and disease status or disease endophenotypes. The most robust GWAS studies of blood-lipid metabolism have focused on circulating total, LDL and HDL cholesterol and triglycerides, which are easily quantified using standard, clinical chemistry assays [11,12]. The extension of GWAS to deeper analyses of lipid species requires mass-spectrometry (MS) technologies and analytical methods that allow for the enhanced efficiency and coverage of lipidome profiling [13]. The application of MS-based lipidomics to GWAS was initially based on targeted analysis of blood sphingomyelins and ceramides [14] and was recently extended to increasing numbers of known lipids [15,16].
Here, we applied liquid chromatography–mass spectrometry (LC–MS) to a group of 273 individuals well-characterized for clinical and biochemical phenotypes relevant to cardiometabolic diseases, to analyse the genetic architecture of lipid metabolism in humans. We were able to identify evidence of the pleiotropy and strong polygenic control of lipids and proposed annotations for lipidomic signals mapped to the human genome. This study demonstrates the power of untargeted lipidomics for high-density quantitative molecular phenotyping in humans and illustrates the complex genetic control of blood-lipid metabolism.

2. Results

2.1. Clinical-Data Analysis

The study group has a mean age of 57.4 ± 0.7 years and 56.4% (n = 154) of the individuals were males (Table 1). All individuals in the cohort were devoid of evidence of coronary artery stenosis, as assessed by an angiogram analysis. Analyses of the pathophysiological components of the cardiometabolic syndrome revealed that 132 individuals (49%) were obese (BMI > 30 kg/m2), 46 had type 2 diabetes (17%), 147 were hypertensive (54%) and 119 were hyperlipidemic (44%), with a similar proportion of affected males and females (Table 2).

2.2. General Features of Untargeted-Lipidome Data

Untargeted-lipidome profiling retrieved 3013 spectral features characterized by a specific mass-to-charge ratio (m/z) and retention time (RT) (1529 in the negative-ionization mode and 1484 in the positive-ionization mode) that met the acceptance criterion (i.e., Relative Standard Deviation (RSD) < 30%, also referred to as Coefficient of Variation CV) (Supplementary Table S1). Multivariate Principal Component Analysis (PCA) analysis showed the absence of strong technical drift during spectral-data acquisition in the cohort, as illustrated by the PCA scores’ 2D plot representation of the QC samples in the two ionization modes (Supplementary Figure S1). The QC samples were tightly clustered, which indicates an acceptable reproducibility of the retained set of metabolic features as well as a good stability of the LC–MS-profiling experiments.

2.3. General Features of Untargeted-Lipidome Data

Genome-wide association of untargeted-lipidome-profiling data identified 5501 statistically significant associations (FDR-adjusted p-value; q-value < 10−8) between SNPs and spectral features (1905 in the negative ionization mode and 3596 in the positive ionization mode). Further analyses of lipid features and their isotopes reduced the analyses to 926 significant associations, between 551 distinct SNP loci and apparently independent lipidome features (Figure 1) (Supplementary Table S2). Eventually, only 74 lipidome features showed evidence of statistical association (q-value < 10−8) to a genetic locus in the cohort (25 in the negative ionization mode and 49 in the positive ionization mode) (Table 3).
Evidence of polygenic control was observed for 52 lipidome features (Table 3), as illustrated with the compound detected, m/z: 277.22 (negative-ionization mode), which was controlled by genetic loci in chromosomes 6 (rs7749100 in DST, q-value = 1.903 × 10−13), 13 (rs1410818, q-value = 4.31 × 10−10) and 20 (rs11699738 in SOGA1, q-value = 4.75 × 10−9) (Figure 2, Supplementary Table S2). Such strong polygenic regulations of lipid metabolism are further illustrated in Figure 3A, with the associations of m/z 271.23, 345.24 and 828.58 (negative-ionization mode), with multiple distinct genetic loci. The compound characterized by an m/z of 345.24 was significantly associated with eight distinct genetic loci on chromosomes 2 (rs2005181 in BABAM2, q-value = 5.68 × 10−10), 4 (rs292037, q-value = 1.93 × 10−13 and rs12500579 in ANK2, q-value = 4.24 × 10−9), 5 (rs10076673 in PITX1, q-value = 7.40 × 10−12), 6 (rs7749100 in DST, q-value = 3.11 × 10−10), 7 (rs2069827 in STEAP1B, q-value = 1.23 × 10−11), 9 (rs7037093, q-value = 2.59 × 10−12) and 13 (rs1410818, q-value = 1.38 × 10−14) (Figure 3A, Supplementary Table S2).
The remaining 22 lipidomic features exhibited evidence of monogenic control. For example, several lipidomic signals acquired by the positive-ionization mode were controlled by a single marker locus on chromosomes 2 (rs12997234 in DPP10 with m/z 568.340 and 590.3213), 3 (rs2002218 in IQSEC1 with m/z 712.645), 5 (rs13362253 in MSX2 with m/z 766.574), 6 (rs7759479 in DST with m/z 279.232, rs6928180 in GRIK2 with m/z 344.279, 370.295, 398.326, 400.342 and 426.357, rs1009439 in RCAN2 with m/z 377.266 and m/z 379.282), 8 (rs6992234 with m/z 204.123), 15 (rs11855528 in CMIP with m/z 612.556 and rs11071737 in RAB8B with m/z 932.864), 16 (rs2292329 in NECAB2 with m/z 922.785) and 20 (rs2260930 in SEL1L2 with m/z 780.553) (Table 3).

2.4. Genetic Analysis of Lipid Metabolism Uncovers Evidence of Pleiotropy

We identified 44 SNP loci that control two or more metabolic features, indicating potential pleiotropic effects of genetic variants, as illustrated in Figure 3B, where closely linked SNPs on chromosomes 6 and 13 are associated with a different m/z. For example, the above-mentioned SNP rs6928180 in GRIK2 was associated with several lipidome features under monogenic control (m/z 344.279, q-value = 1.89 × 10−23; m/z 370.295, q-value = 1.14 × 10−32; m/z 398.326, q-value = 4.96 × 10−34; m/z 400.342, q-value = 3.68 × 10−28; m/z 426.357, q-value = 7.38 × 10−18) suggesting a pleiotropic effect of variants in GRIK2 on distinct but coordinately regulated lipids (Table 3). Along the same line, marker rs12997234 on chromosome 2 in an intron of DPP10 was exclusively associated with the monogenic control of m/z 568.34 (q-value = 1.73 × 10−11) and m/z 590.32 (q-value = 2.93 × 10−17) in the positive-ionization mode and with m/z 612.33 (q-value = 1.46 × 10−9) in the negative-ionization mode (Table 3). The most striking example of pleiotropy was detected on chromosome 13 at the locus rs1410818 and 11 distinct m/z values (Supplementary Table S2).

2.5. Assignment of Lipids to Lipidomic Features Mapped to the Human Genome

We next carried out the identification of candidate lipids for each of the 74 features showing evidence of genetic control. Using the LIPID MAPS database, we were able to annotate 26 lipidome signals with a single lipid, including 10 which were controlled by a single genetic locus (Table 3). Several lipid candidates could be proposed for the remaining 48 lipidome features, which prevented the unambiguous assignment of lipids. The vast majority of assigned lipids were fatty acyls (27), sterol lipids (23), triacylgycerols (9) and, to a lesser extent, a combination of phosphatidylcholines, phosphatidylethanolamine and phosphatidylserines (20).

2.6. Metabolome-Wide Association Studies Identify Metabolites Associated with Clinical and Biochemical Phenotypes

To test for evidence of association between clinical and variations in biochemical phenotypes and compounds from the lipidome dataset mapped to the human genome, linear regression was performed. Results from associations with a nominal p < 0.05 are given in Supplementary Table 3. Significant associations (q-value < 0.05) with multiple metabolic features were detected for cardiometabolic disease (Table 4). Fewer significant associations were identified for family history of hypertension (m/z 695.511 and 938.536) and for variations in body-mass index (m/z 774.543, 833.588, 834.591 and 832.584), total cholesterol (m/z 758.569 and 759.572) and HDL cholesterol (m/z 367.228 and 213.146) (Figure 4, Table 4). Family history of diabetes also showed evidence of marginal association to the feature m/z 695.511 (nominal p-value = 0.036) (Supplementary Table S3). Associations to family history of hypertension and diabetes independent to association to the diseases suggest that the underlying lipidomic feature may be a predictive marker of both diseases.
We did not identify statistically significant associations to LDL cholesterol or triacylglycerols. However, over 60 lipidomic features showed marginal evidence of co-association (nominal p-value < 0.05) to both LDL and HDL cholesterol (e.g., m/z 129.98 and 171.99) and five features (m/z 213.15, 367.23, 367.26, 369.27 and 722.50) were marginally associated to triacylglycerols and total, HDL and LDL cholesterol (Supplementary Table S3). No significant associations were found between spectral data and other phenotypes. We were able to assign one or several putative lipids to 14 lipidome signals, including ST 27:2;O;Hex and ST 28:1;O5, which were found to be regulated by multiple genetic loci (Table 4).

3. Discussion

We report results from the genome mapping of untargeted serum lipidomics in a group of individuals characterized for pathophysiological features of the cardiometabolic syndrome. We identified evidence of strong polygenic control of lipid features and instances of mechanisms of pleiotropy in the regulation of lipid metabolism. These observations illustrate the complex genetic architecture of serum lipid regulation and provide novel information beyond the genetic control of cholesterol metabolism.
Both proton nuclear magnetic resonance (1H NMR) and mass spectrometry (MS) have been successfully used to map the genetic control of predominantly serum metabolites in genome-wide association studies (GWAS) in humans [17]. Collectively, over 1800 metabolomic data (i.e., known and unknown metabolites and ratios) have been associated with over 40,000 unique SNPs [18]. Among these, MS-lipidomic data provide significant advances in our understanding of the etiopathogenesis of diseases characterized by anomalies in lipid metabolism [19]. Untargeted lipidomics, a hypothesis-free strategy that has the power of deepening quantitative lipid analyses to unassigned lipids, remains challenging due to the breadth and intrinsic complexity of known lipids, which differ in terms of physicochemical properties [13,20]. As a consequence, harmonization of sample preparation for such a heterogeneous group of molecules is a problematic issue that limits detection and quantification of the broad diversity of lipid species [21]. In addition, variations in MS-instrument stability affect repeatability within and between experiments. Finally, the unambiguous assignment of putative lipids to MS-spectral signals remains an important methodological consideration in the application of untargeted MS lipidomics in GWAS.
Polygenic control is a hallmark of GWAS of human chronic diseases and complex phenotypes, and the genetic regulation of metabolomic profiling data does not make any exceptions [22,23,24]. We show that serum-lipid abundance exhibits predominant polygenic control, when a single metabolite is associated with several unlinked SNPs. Results from lipidomic GWAS have shown that about 30% of lipids are associated with several genetic loci [16]. Specifically, loci on chromosomes 2 and 4 control triglyceride TAG(50:1;0), loci on chromosomes 8 and 11 are associated with triglyceride TAG(52:3;0) and loci on chromosomes 12 and 18 control lysophosphatidylcholine LPC(14:0;0) [15]. This pattern of polygenic control suggests either functional redundancy of proteins in the regulation of lipid metabolic pathways, or the involvement of distinct proteins each contributing in parallel or in concert to interconnected mechanisms of lipid sensing, synthesis, transport and degradation.
Our association results also suggest apparent pleiotropy when a single genetic locus controls multiple, different lipidomic features. It is expected to occur in metabolic processes, since altered regulation of an individual protein involved in an enzymatic reaction or metabolite binding or transport may result in changes in interconnected biological pathways affecting multiple metabolites. An excess of distinct lipid species associated with genomic regions in lipidomic GWAS suggests the widespread occurrence of this phenomenon in the regulation of lipid metabolism [15,23,24]. Harshfield et al. reported the genetic mapping of 181 lipids to only 24 genomic regions [16], and Tabassum et al. identified associations to 42 lipid species in 11 genomic regions [15], thus implying that one genomic region is associated with several lipids. One of the most striking examples of pleiotropy in lipidomic GWAS is the GCKR locus, which is associated with over 30 lipid species [25]. The eicosanoid metabolic network, which involves 28 proteins for the production of over 150 lipids, provides a further example of pleiotropy in the regulation of lipid biology [26]. These coordinately regulated lipid clusters suggest the existence of genetically-determined “lipidotypes”.
Combined with clinical data, lipidomic-based phenotyping allows the definition of disease-associated biomarkers as well as druggable-metabolite targets. Integrating genotyping data can identify instances of co-localization of disease-risk SNPs and loci associated with metabolomic features, which may represent disease-causative molecular biomarkers [15,16,27]. With the exception of SEL1L2 and SYT9, gene loci showing evidence of monogenic control of lipids in our study have been associated with disease-relevant phenotypes (e.g., body mass index), biochemical variables (e.g., creatinine) and behavioral traits in the GWAS repository (www.ebi.ac.uk/gwas/, accessed on 1 May 2022). Interestingly, multiple SNPs, the locus of the gene encoding pleckstrin and the Sec7 domain containing 3 (PSD3), which controls the level of a carnitine in our study, have been consistently associated with triglycerides and cholesterol levels as well as type 2 diabetes and obesity [28], and their downregulation results in reduced hepatic lipids in vitro and protects against fatty liver in vivo in mice [29].
Considering the breadth of circulating lipid species [7,21] and their roles in cardiovascular diseases [19], we were able to map the genetic control of several lipid species, mostly fatty acyls, phospholipids and triglycerides. On the other hand, we were unable to identify genetic loci associated with several important lipid species, including, for example, sphingomyelins and ceramides, which are involved in cardiovascular risk [9,30]. This may be caused by technical issues with data acquisition and the relatively modest sample size of the study but may also be accounted for by specific clinical features of the individuals selected in our study. Absence of coronary-artery stenosis in these individuals suggests reduced cardiovascular risk and, therefore, potentially limited quantitative variations in blood ceramides in cases and controls that may prevent genetic mapping. In support of this hypothesis, we did not identify statistically significant associations between lipidomic features and hypertension, which might nevertheless be improved with the use of intermediate, quantitative phenotypes, including measures of blood pressure. In addition, the fact that CMD patients may be under various medications, including lipid-lowering drugs (statins) or anti-diabetic treatments that result in improved control of blood pressure [31], may explain the absence of statistically significant associations between lipidomic features and hypertension in our study. However, our results suggest a role of lipids in the family history of hypertension, which may represent disease-predictive markers.

4. Materials and Methods

4.1. Study Subjects

The study group consisted of 273 subjects selected from a larger study recruited between 2006 and 2009 for inclusion in the FGENTCARD patient collection, primarily designed to map the genetic determinants of coronary artery stenosis [32]. Individuals from the FGENTCARD cohort were originally referred to a catheterization care unit for clinical evaluation. A 20 mL blood sample was collected in overnight fasted individuals from the peripheral femoral artery during the coronary angiography for serum preparation. Patients provided a written consent for the whole study including genomic analyses. The Institutional Review Board (IRB) at the Lebanese American University approved the study protocol.
Body weight, body-mass index (BMI) and blood chemistry (total, HDL and LDL cholesterol, triglycerides) were determined. Evidence of diabetes (fasting glucose > 125 mg/dl), hypertension (blood pressure > 10/14 mm Hg) and obesity (BMI > 30) was recorded in individuals’ medical charts. Evidence of cardiometabolic disease (CMD) was assessed by presence of at least three anomalies (diabetes, hypertension, BMI > 30 kg/m2 and HDL < 40 mg/dl). All 273 individuals selected for this genetic study were devoid of vessel stenosis, assessed through coronary angiography carried out at a single recruitment site. Family history of diabetes and hypertension, defined by presence of the disease in a sibling, parent or second-degree relative, was also recorded.
Statistical analysis of clinical and biochemical data was performed using two-way ANOVA. Differences were considered statistically significant with a p < 0.05.

4.2. Chemicals

Isopropanol, acetonitrile, formic acid and ammonium formate were LC–MS Chromasolv® Fluka and high-performance liquid chromatography (HPLC) quality and were purchased from Sigma-Aldrich (Sigma-Aldrich, Saint-Quentin Fallavier, France). Ultra-pure water (resistivity: 18 mΩ) was obtained with a Milli-Q Integral purification system (Millipore, Molsheim, France) fitted with a 0.22 µm filter. The mobile phase was prepared with a solvent containing 400 mL of water, 600 mL of acetonitrile, 0.1% formic acid and 0.630 g of ammonium formate, and a solvent containing 100 mL of acetonitrile, 900 mL of isopropanol, 0.1% formic acid and 0.630 g of ammonium formate.

4.3. Sample Preparation

Lipid extraction from serum was performed using isopropanol (1:6, v/v), as recommended by the MS-equipment supplier, which is the most robust solvent enabling a broad coverage and recovery of lipid species from serum [33]. Experiments were carried out with 50 µL serum aliquots. Samples were then centrifuged at 14,000 g, and supernatants were then transferred to vials for injection in the UPLC system.

4.4. UPLC analysis

A Waters Acquity UPLC® (Waters Corp, Saint-Quentin en Yvelines, France) fitted with a Acquity CSH C18 column (2.1 × 150 mm, 1.7 µm) and a corresponding guard column (Acquity CSH 1.7μM) (Waters Corp, Saint-Quentin en Yvelines, France) were used to analyse lipid compounds in serum samples as previously described [34]. The oven temperature was set at 55 °C. The flow rate used for these experiments was 400 µL/min and a volume of 5 µL of sample was injected. The total run time was 24 min. A binary gradient consisted of above-described mobile phases was used according to Waters’ recommendation. Mobile phase B was maintained at 99% during 4 min at the end of the gradient.

4.5. Mass Spectrometry

Mass spectrometry was carried out as previously [34]. The UPLC system was coupled with a Q-Exactive™ Hybrid Quadrupole-Orbitrap mass spectrometer (Thermo Fisher Scientific, Illkirch, France). Infusion of a calibration mixture (caffeine, MRFA and Ultramark® 1621) was used for calibration of the instrument. Parameters of the heated-electrospray (HESI-II, Thermo Fisher Scientific, Illkirch, France) interface were as follows: S-Lens 50 V, Sheat gas: 65, Auxiliary gas: 25 arbitrary units, capillary voltage 3 kV, capillary temperature 350 °C and vaporization temperature 60 °C. The maximum target capacity of the C-trap (autogain control, AGC) target was defined as 3e6 ions and the maximum injection time was 200 ms. Full scans were obtained in positive and negative ion modes simultaneously with a resolution of 70,000 full width at half maximum (FWHM), in the scan range of mass-to-charge ratio (m/z) of 85–1275.

4.6. Untargeted Lipidomic Data Analysis

Analysis of MS data derived from UPLC complied with standard protocols and food and drug administration (FDA) guidelines [35,36], as previously described (34). XCMS tools implemented in R statistical language (v 3.1.0) (http://www.bioconductor.org, accessed on 10 May 2020) were used for preprocessing steps of MS data analysis (peak picking, peak grouping, retention-time correction, annotation of isotopes and adducts). Profiles of positive and negative ionization modes were separately extracted and converted into mzXML format for preprocessing by the XCMS tools. Identification of Regions of Interest (ROI) used the wavelet-based peak-picking approach (centwave). MS-data preprocessing resulted in a peak table listing lipidomic features characterized by a retention time (RT), mass-to-charge ratio (m/z) and corresponding intensity for each serum sample.
A data matrix reduction was applied to retain spectral features consistently found in the individuals. Over 40% of missing values were withdrawn. Performance and reliability of the analytical process and compliance of data with FDA-acceptance criteria [37] were also verified through a quality assurance (QA) strategy, based on analysis of a pooled quality-control (QC) sample, which was injected every 10 samples throughout the analytical run. Median fold-change-normalization approach [38] was applied on the retained MS features, followed by a generalized log-transformation. A threshold of 30% calculated for each metabolic feature in the QC samples was set for relative standard deviation (RSD), which is an accepted standard to assess data reproducibility in metabolomic studies [35,36]. Four samples were identified as outliers and were discarded from the study. The resulting matrix was then used for multivariate and univariate statistical analyses (principal component analysis and linear regression).

4.7. Metabolome-Genome Wide Association Studies (mGWAS)

All individuals were genotyped by Illumina Human610-Quad BeadChip and Illumina Human660W-Quad BeadChip, respectively (552,510 overlapping SNPs), as part of the FGENTCARD consortium [32]. All SNPs with over 98% genotyping success rate, minor allele frequency above 1% and in Hardy-Weinberg equilibrium (p-value > 1 × 10−7) were included in the analysis. An imputation across the whole genome to CEU HapMap population as a reference was performed using the IMPUTE2 tool [39], which yielded 2,573,690 SNPs. The plink tool [40] was used to perform both association analyses based on an additive genetic model. An FDR adjusted p-value (q-value) < 1 × 10−8 was considered to be significant genome-wide. Plotting circles were generated using an in-house tool specifically developed to illustrate mGWAS associations.

4.8. Metabolome-Wide Association Studies (MWAS)

A linear-regression model was applied to carry out MWAS through the assessment of association, between each metabolic feature with clinical and biochemical continuous phenotypes (total, HDL and LDL cholesterol, triglycerides). Normality assumption of the residuals of each metabolic feature was investigated by Shapiro–Wilk test. The R statistical language was used to perform the linear regression and compute a p-value for each metabolic feature with a threshold of significance set to 0.05. Adjustment for age and sex was performed by including them as covariates in the statistical model. False discovery rates (FDR) were corrected using the Benjamini-Hochberg method to adjust P-values for false discovery involving multiple comparisons.

4.9. Assignment of Lipid Features

Annotation of lipid candidates corresponding to lipidome signals was carried out using the free resource LIPID MAPS (https://www.lipidmaps.org, accessed on 1 May 2022). We initially performed bulk-structure searches and subsequently refined our analysis by interrogating the LIPID MAPS Structure Database (LMSD) with a list of precursor ions. We entered the list of precursor ion m/z and chose appropriate polarity for the adduct ions. We defined a mass tolerance of ±0.001 m/z and sorted our data according to the delta between the input m/z and the m/z of candidate proposed in the database.

5. Conclusions

Results from our untargeted-lipidomic profiling provide information on fundamental mechanisms regulating serum lipids in humans. Replication of these findings in larger study populations and further analyses, such as MS/MS validation experiments designed to unambiguously assign lipids to lipidomic features, are required.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo12070596/s1. Figure S1: 2-D Principal component analysis of mass spectrometry data in the cohort representing the scores of the first components; Table S1: Lipidome spectral features acquired in serum samples from the study population; Table S2: List of all SNPs showing evidence of statistically significant association with lipidome fearures; Table S3: Associations between lipidomic features and clinical and biochemical data.

Author Contributions

P.Z., F.M. and D.G. conceived the study. D.G. wrote the manuscript. P.Z. provided patient serum samples. F.B., K.S., C.A.H. and L.H. performed metabolomic-data processing, statistical-data analyses, and genetic-association analyses. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge financial support of the European Commission for collection of the patient cohort (FGENTCARD, LSHGCT-2006-037683). F.M. and D.G. acknowledge the financial support from the Inserm “Projet de Recherche International” Diabetomarkers.

Institutional Review Board Statement

Patients provided a written consent for the whole study including genomic analyses. The Institutional Review Board (IRB) at the Lebanese American University approved the study protocol.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data is contained within the article or supplementary material. The data presented in this study are available in Supplementary Tables.

Conflicts of Interest

The authors declare no competing financial interests.

References

  1. Akiyama, M. Multi-omics study for interpretation of genome-wide association study. J. Hum. Genet. 2021, 66, 3–10. [Google Scholar] [CrossRef] [PubMed]
  2. Cookson, W.; Liang, L.; Abecasis, G.; Moffatt, M.; Lathrop, M. Mapping complex disease traits with global gene expression. Nat. Rev. Genet. 2009, 10, 184–194. [Google Scholar] [CrossRef] [PubMed]
  3. Suhre, K.; Raffler, J.; Kastenmüller, G. Biochemical insights from population studies with genetics and metabolomics. Arch. Biochem. Biophys. 2016, 589, 168–176. [Google Scholar] [CrossRef]
  4. Meikle, P.J.; Christopher, M.J. Lipidomics is providing new insight into the metabolic syndrome and its sequelae. Curr. Opin. Lipidol. 2011, 22, 210–215. [Google Scholar] [CrossRef] [PubMed]
  5. Baek, J.; He, C.; Afshinnia, F.; Michailidis, G.; Pennathur, S. Lipidomic approaches to dissect dysregulated lipid metabolism in kidney disease. Nat. Rev. Nephrol. 2022, 18, 38–55. [Google Scholar] [CrossRef]
  6. Xiao, C.; Rossignol, F.; Vaz, F.M.; Ferreira, C.R. Inherited disorders of complex lipid metabolism: A clinical review. J. Inherit. Metab. Dis. 2021, 44, 809–825. [Google Scholar] [CrossRef]
  7. Quehenberger, O.; Dennis, E.A. The human plasma lipidome. N. Engl. J. Med. 2011, 365, 1812–1823. [Google Scholar] [CrossRef]
  8. Meikle, T.G.; Huynh, K.; Giles, C.; Meikle, P.J. Clinical lipidomics: Realizing the potential of lipid profiling. J. Lipid Res. 2021, 62, 100127. [Google Scholar] [CrossRef]
  9. McGurk, K.A.; Keavney, B.D.; Nicolaou, A. Circulating ceramides as biomarkers of cardiovascular disease: Evidence from phenotypic and genomic studies. Atherosclerosis 2021, 327, 18–30. [Google Scholar] [CrossRef]
  10. Yang, L.; Wang, L.; Deng, Y.; Sun, L.; Lou, B.; Yuan, Z.; Wu, Y.; Zhou, B.; Liu, J.; She, J. Serum lipids profiling perturbances in patients with ischemic heart disease and ischemic cardiomyopathy. Lipids Health Dis. 2020, 19, 89. [Google Scholar] [CrossRef]
  11. Teslovich, T.M.; Musunuru, K.; Smith, A.V.; Edmondson, A.C.; Stylianou, I.M.; Koseki, M.; Pirruccello, J.P.; Ripatti, S.; Chasman, D.I.; Willer, C.J.; et al. Biological, clinical and population relevance of 95 loci for blood lipids. Nature 2010, 466, 707–713. [Google Scholar] [CrossRef] [PubMed]
  12. Willer, C.J.; Schmidt, E.M.; Sengupta, S.; Peloso, G.M.; Gustafsson, S.; Kanoni, S.; Ganna, A.; Chen, J.; Buchkovich, M.L.; Mora, S.; et al. Discovery and refinement of loci associated with lipid levels. Nat. Genet. 2013, 45, 1274–1283. [Google Scholar] [PubMed] [Green Version]
  13. Zhao, Y.Y.; Wu, S.P.; Liu, S.; Zhang, Y.; Lin, R.C. Ultra-performance liquid chromatography-mass spectrometry as a sensitive and powerful technology in lipidomic applications. Chem. Biol. Interact. 2014, 220, 181–192. [Google Scholar] [CrossRef] [PubMed]
  14. Demirkan, A.; van Duijn, C.M.; Ugocsai, P.; Isaacs, A.; Pramstaller, P.P.; Liebisch, G.; Wilson, J.F.; Johansson, Å.; Rudan, I.; Aulchenko, Y.S.; et al. Genome-wide association study identifies novel loci associated with circulating phospho- and sphingolipid concentrations. PLoS Genet. 2012, 8, e1002490. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Tabassum, R.; Rämö, J.T.; Ripatti, P.; Koskela, J.T.; Kurki, M.; Karjalainen, J.; Palta, P.; Hassan, S.; Nunez-Fontarnau, J.; Kiiskinen, T.T.J.; et al. Genetic architecture of human plasma lipidome and its link to cardiovascular disease. Nat. Commun. 2019, 10, 4329. [Google Scholar] [CrossRef] [Green Version]
  16. Harshfield, E.L.; Fauman, E.B.; Stacey, D.; Paul, D.S.; Ziemek, D.; Ong, R.M.Y.; Danesh, J.; Butterworth, A.S.; Rasheed, A.; Sattar, T.; et al. Genome-wide analysis of blood lipid metabolites in over 5000 South Asians reveals biological insights at cardiometabolic disease loci. BMC Med. 2021, 19, 232. [Google Scholar] [CrossRef]
  17. Kastenmüller, G.; Raffler, J.; Gieger, C.; Suhre, K. Genetics of human metabolism: An update. Hum. Mol. Genet. 2015, 24, R93–R101. [Google Scholar] [CrossRef] [Green Version]
  18. Hagenbeek, F.A.; Pool, R.; van Dongen, J.; Draisma, H.H.M.; Jan Hottenga, J.; Willemsen, G.; Abdellaoui, A.; Fedko, I.O.; den Braber, A.; Visser, P.J.; et al. Heritability estimates for 361 blood metabolites across 40 genome-wide association studies. Nat. Commun. 2020, 11, 39. [Google Scholar] [CrossRef] [Green Version]
  19. Tabassum, R.; Ripatti, S. Integrating lipidomics and genomics: Emerging tools to understand cardiovascular diseases. Cell. Mol. Life Sci. 2021, 78, 2565–2584. [Google Scholar] [CrossRef]
  20. Mundra, P.A.; Shaw, J.E.; Meikle, P.J. Lipidomic analyses in epidemiology. Int. J. Epidemiol. 2016, 45, 1329–1338. [Google Scholar] [CrossRef] [Green Version]
  21. Burla, B.; Arita, M.; Bendt, A.K.; Cazenave-Gassiot, A.; Dennis, E.A.; Ekroos, K.; Han, X.; Ikeda, K.; Liebisch, G.; Lin, M.K. MS-based lipidomics of human blood plasma: A community-initiated position paper to develop accepted guidelines. J. Lipid Res. 2018, 59, 2001–2017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Lotta, L.A.; Pietzner, M.; Stewart, I.D.; Wittemans, L.B.L.; Li, C.; Bonelli, R.; Raffler, J.; Biggs, E.K.; Oliver-Williams, C.; Auyeung, V.P.W.; et al. A cross-platform approach identifies genetic regulators of human metabolism and health. Nat. Genet. 2021, 53, 54–64. [Google Scholar] [CrossRef] [PubMed]
  23. Gallois, A.; Mefford, J.; Ko, A.; Vaysse, A.; Julienne, H.; Ala-Korpela, M.; Laakso, M.; Zaitlen, N.; Pajukanta, P.; Aschard, H. A comprehensive study of metabolite genetics reveals strong pleiotropy and heterogeneity across time and context. Nat. Commun. 2019, 10, 4788. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kanai, M.; Akiyama, M.; Takahashi, A.; Matoba, N.; Momozawa, Y.; Ikeda, M.; Iwata, N.; Ikegawa, S.; Hirata, M.; Matsuda, K.; et al. Genetic analysis of quantitative traits in the Japanese population links cell types to complex human diseases. Nat. Genet. 2018, 50, 390–400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Fernandes Silva, L.; Vangipurapu, J.; Kuulasmaa, T.; Laakso, M. An intronic variant in the GCKR gene is associated with multiple lipids. Sci. Rep. 2019, 9, 10240. [Google Scholar] [CrossRef] [Green Version]
  26. Buczynski, M.W.; Dumlao, D.S.; Dennis, E.A. Thematic Review Series: Proteomics. An integrated omics analysis of eicosanoid biology. J. Lipid Res. 2009, 50, 1015–1038. [Google Scholar] [CrossRef] [Green Version]
  27. Shin, S.Y.; Fauman, E.B.; Petersen, A.K.; Krumsiek, J.; Santos, R.; Huang, J.; Arnold, M.; Erte, I.; Forgetta, V.; Yang, T.; et al. An atlas of genetic influences on human blood metabolites. Nat. Genet. 2014, 46, 543–550. [Google Scholar] [CrossRef] [Green Version]
  28. Gong, S.; Xu, C.; Wang, L.; Liu, Y.; Owusu, D.; Bailey, B.A.; Lid, Y.; Wang, K. Genetic association analysis of polymorphisms in PSD3 gene with obesity, type 2 diabetes, and HDL cholesterol. Diabetes Res. Clin. Pract. 2017, 126, 105–114. [Google Scholar] [CrossRef]
  29. Mancina, R.M.; Sasidharan, K.; Lindblom, A.; Wei, Y.; Ciociola, E.; Jamialahmadi, O.; Pingitore, P.; Andréasson, A.; Pellegrini, G.; Baselli, G.; et al. PSD3 downregulation confers protection against fatty liver disease. Nat. Metab. 2022, 4, 60–75. [Google Scholar] [CrossRef]
  30. Wang, D.D.; Toledo, E.; Hruby, A.; Rosner, B.A.; Willett, W.C.; Sun, Q.; Razquin, C.; Zheng, Y.; Ruiz-Canela, M.; Guasch-Ferré, M.; et al. Plasma Ceramides, Mediterranean Diet, and Incident Cardiovascular Disease in the PREDIMED Trial (Prevención con Dieta Mediterránea). Circulation 2017, 135, 2028–2040. [Google Scholar] [CrossRef] [Green Version]
  31. Maloberti, A.; Bruno, R.M.; Facchetti, R.; Grassi, G.; Taddei, S.; Ghiadoni, L.; Giannattasio, C. The role of metabolic syndrome in blood pressure control and pulse wave velocity progression over a 3.5 years in treated hypertensive patients. Eur. J. Intern. Med. 2020, 76, 107–109. [Google Scholar] [CrossRef] [PubMed]
  32. Hager, J.; Kamatani, Y.; Cazier, J.B.; Youhanna, S.; Ghassibe-Sabbagh, M.; Platt, D.E.; Abchee, A.B.; Romanos, J.; Khazen, G.; Othman, R.; et al. Genome-wide association study in a Lebanese cohort confirms PHACTR1 as a major determinant of coronary artery stenosis. PLoS ONE 2012, 7, e38663. [Google Scholar]
  33. Sarafian, M.H.; Gaudin, M.; Lewis, M.R.; Martin, F.P.; Holmes, E.; Nicholson, J.K.; Dumas, M.E. Objective set of criteria for optimization of sample preparation procedures for ultra-high throughput untargeted blood plasma lipid profiling by ultra performance liquid chromatography-mass spectrometry. Anal. Chem. 2014, 86, 5766–5774. [Google Scholar] [CrossRef] [PubMed]
  34. Zalloua, P.; Kadar, H.; Hariri, E.; Abi Farraj, L.; Brial, F.; Hedjazi, L.; le Lay, A.; Colleu, A.; Dubus, J.; Touboul, D.; et al. Untargeted Mass Spectrometry Lipidomics identifies correlation between serum sphingomyelins and plasma cholesterol. Lipids Health Dis. 2019, 18, 38. [Google Scholar] [CrossRef]
  35. Dunn, W.B.; Broadhurst, D.; Begley, P.; Zelena, E.; Francis-McIntyre, S.; Anderson, N.; Brown, M.; Knowles, J.D.; Halsall, A.; Haselden, J.N.; et al. Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat. Protoc. 2011, 6, 1060–1083. [Google Scholar] [CrossRef]
  36. Want, E.J.; Wilson, I.D.; Gika, H.; Theodoridis, G.; Plumb, R.S.; Shockcor, J.; Holmes, E.; Nicholson, J.K. Global metabolic profiling procedures for urine using UPLC-MS. Nat. Protoc. 2010, 5, 1005–1018. [Google Scholar] [CrossRef]
  37. Dunn, W.B.; Wilson, I.D.; Nicholls, A.W.; Broadhurst, D. The importance of experimental design and QC samples in large-scale and MS-driven untargeted metabolomic studies of humans. Bioanalysis 2012, 4, 2249–2264. [Google Scholar] [CrossRef] [Green Version]
  38. Veselkov, K.A.; Vingara, L.K.; Masson, P.; Robinette, S.L.; Want, E.; Li, J.V.; Barton, R.H.; Boursier-Neyret, C.; Walther, B.; Ebbels, T.M.; et al. Optimized preprocessing of ultra-performance liquid chromatography/mass spectrometry urinary metabolic profiles for improved information recovery. Anal. Chem. 2011, 83, 5864–5872. [Google Scholar] [CrossRef]
  39. Howie, B.N.; Donnelly, P.; Marchini, J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet. 2009, 5, e1000529. [Google Scholar] [CrossRef] [Green Version]
  40. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
Metabolites 12 00596 g001 550
Figure 1. Genome-wide association study of metabolomic features (mGWAS) in the study group. Data are shown for metabolic features acquired in positive (A) and negative (B) ionization modes, showing evidence of significant association (LOD > 8) with an SNP locus. Chromosomes are color-coded on the circle. The colors of the lines indicate the chromosomal location of SNP loci showing evidence of significant association with metabolic features, characterized by a mass-to-charge ratio (horizontal axes). Details of genetic results are given in Supplementary Table S2.
Figure 1. Genome-wide association study of metabolomic features (mGWAS) in the study group. Data are shown for metabolic features acquired in positive (A) and negative (B) ionization modes, showing evidence of significant association (LOD > 8) with an SNP locus. Chromosomes are color-coded on the circle. The colors of the lines indicate the chromosomal location of SNP loci showing evidence of significant association with metabolic features, characterized by a mass-to-charge ratio (horizontal axes). Details of genetic results are given in Supplementary Table S2.
Metabolites 12 00596 g001
Metabolites 12 00596 g002 550
Figure 2. Manhattan plot illustrating the polygenic control of metabolic features. Genome-wide association study was carried out with over 2.5 M imputed SNPs, for the metabolomic feature characterized by a mass-to-charge ratio of 227.216 and a retention time of 67.49. Chromosomes are color-coded. Evidence of significant associations (LOD >8) with this metabolic feature were found on chromosomes 1, 5, 6, 13 and 20. The Y-axis corresponds to the significance of the association (−Log10 p-values). The X-axis represents the physical location of the variant colored by chromosome.
Figure 2. Manhattan plot illustrating the polygenic control of metabolic features. Genome-wide association study was carried out with over 2.5 M imputed SNPs, for the metabolomic feature characterized by a mass-to-charge ratio of 227.216 and a retention time of 67.49. Chromosomes are color-coded. Evidence of significant associations (LOD >8) with this metabolic feature were found on chromosomes 1, 5, 6, 13 and 20. The Y-axis corresponds to the significance of the association (−Log10 p-values). The X-axis represents the physical location of the variant colored by chromosome.
Metabolites 12 00596 g002
Metabolites 12 00596 g003 550
Figure 3. Architectural characteristics of genetic associations to metabolic features. Evidence of polygenic control of metabolites (A) and potential pleiotropic effects of genetic loci on metabolite abundance (B) were identified, following metabolomic analysis of serum samples of 273 individuals. The colours of the lines indicate the chromosomal location of SNP loci showing evidence of significant association (LOD > 8), with the abundance of a specific metabolic feature. Evidence of polygenic control of the abundance of metabolic features was found for compounds characterized by mass-to-charge ratios (horizontal axis) of 271.23 (red), 345.24 (blue) and 828.58 (purple) (A). Potential pleiotropic effects were detected for SNP loci on chromosomes 6 (red lines) and 13 (blue lines), significantly associated with metabolic features characterized by distinct mass-to-charge ratios on the horizontal axis (B). Details of genetic results are given in Supplementary Table S2.
Figure 3. Architectural characteristics of genetic associations to metabolic features. Evidence of polygenic control of metabolites (A) and potential pleiotropic effects of genetic loci on metabolite abundance (B) were identified, following metabolomic analysis of serum samples of 273 individuals. The colours of the lines indicate the chromosomal location of SNP loci showing evidence of significant association (LOD > 8), with the abundance of a specific metabolic feature. Evidence of polygenic control of the abundance of metabolic features was found for compounds characterized by mass-to-charge ratios (horizontal axis) of 271.23 (red), 345.24 (blue) and 828.58 (purple) (A). Potential pleiotropic effects were detected for SNP loci on chromosomes 6 (red lines) and 13 (blue lines), significantly associated with metabolic features characterized by distinct mass-to-charge ratios on the horizontal axis (B). Details of genetic results are given in Supplementary Table S2.
Metabolites 12 00596 g003
Metabolites 12 00596 g004 550
Figure 4. Metabolome-wide association studies (MWAS) in patients with cardiometabolic syndrome. Correlations were tested between clinical and biochemical phenotypes and serum metabolic features characterized by a mass-to-charge ratio (m/z) shown on the x-axes. Data are shown for body-mass index (A), family history of hypertension (B), total cholesterol (C) and HDL cholesterol (D,E). The Y-axis corresponds to the adjusted false-discovery rate (FDR). Regression analysis was adjusted for age and sex effects by including them as covariates in the model. pos, positive ionization mode; neg, negative ionization mode.
Figure 4. Metabolome-wide association studies (MWAS) in patients with cardiometabolic syndrome. Correlations were tested between clinical and biochemical phenotypes and serum metabolic features characterized by a mass-to-charge ratio (m/z) shown on the x-axes. Data are shown for body-mass index (A), family history of hypertension (B), total cholesterol (C) and HDL cholesterol (D,E). The Y-axis corresponds to the adjusted false-discovery rate (FDR). Regression analysis was adjusted for age and sex effects by including them as covariates in the model. pos, positive ionization mode; neg, negative ionization mode.
Metabolites 12 00596 g004
Table
Table 1. Clinical and biochemical features of individuals in the study group used for metabolomic profiling. Individuals were selected for absence of coronary stenosis. Data are given as means ± SEM. Number of cases are reported in parentheses. Gender differences were tested using two-way ANOVA.
Table 1. Clinical and biochemical features of individuals in the study group used for metabolomic profiling. Individuals were selected for absence of coronary stenosis. Data are given as means ± SEM. Number of cases are reported in parentheses. Gender differences were tested using two-way ANOVA.
AllFemalesMales
MeanRangeMeanRangeMeanRange
Age57.4 ± 0.7 (273)30–8361.4 ± 0.9 (119)38–8354.4 ± 0.9 (154)30–81
Body weight (kg)83.13 ± 0.99 (269)50–15077.69 ± 1.44 (118)52–15087.39 ± 1.26 (151)50–130
BMI (kg/m2)30.37 ± 0.33 (268)18.96–55.7731.36 ± 0.56 (118)20.34–55.7729.59 ± 0.37 (150)18.96–44.29
Glucose (mg/dL)107.95 ± 2.19 (219)60–299111.41 ± 3.98 (98)62–299105.14 ± 2.29 (121)60–255
Total cholesterol (mg/dL)187.89 ± 2.83 (266)71–357196.35 ± 4.12 (114)71–345181.55 ± 3.81 (152)76–357
HDL cholesterol (mg/dL)41.87 ± 0.80 (266)18–9046.10 ± 1.22 (115)18–8538.65 ± 0.98 (151)18–90
LDL cholesterol (mg/dL)113.90 ± 2.29 (261)24–254117.21 ± 3.22 (115)34–240111.29 ± 3.21 (146)24–254
Triglycerides (mg/dL)176.58 ± 7.03 (273)9–1215167.87 ± 8.12 (119)9–580183.30 ± 10.77 (154)9–1215
Table
Table 2. Pathophysiological components and risk factors of the cardiometabolic syndrome in individuals of the study group. Number of cases is reported and percentages are given in parentheses.
Table 2. Pathophysiological components and risk factors of the cardiometabolic syndrome in individuals of the study group. Number of cases is reported and percentages are given in parentheses.
AllMalesFemales
Body mass index > 30 (kg/m2)132 (49%)66 (44%)66 (56%)
HDL cholesterol < 40 (mg/dl)128 (48%)94 (62%)34 (30%)
Fasting glycemia > 125 mg/dl36 (16%)16 (13%)20 (20%)
Type 2 diabetes46 (17%)23 (15%)23 (19%)
Hypertension147 (54%)73 (47%)74 (62%)
Hyperlipidemia119 (44%)67 (44%)52 (44%)
Family history of hypertension187 (69%)99 (64%)88 (74%)
Family history of type 2 diabetes155 (57%)83 (54%)72 (61%)
Table
Table 3. Genetic control of lipidomic signals mapped to the genome and proposed lipid assignments. Lipidome data, acquired with a Q Exactive Hybrid Quadrupole-Orbitrap mass spectrometer fitted with a Waters Acquity CSH C18 column, were tested for genetic association with genotyped SNPs in the study group (n = 273). Features were characterized by their retention time (RT) and their mass-to-charge ratio (m/z). Details of SNPs and statistical significance of lipidome features under monogenic control are reported. Full list of genetically mapped LC–MS lipidomic features and details and distinct SNP markers associated with lipid features under polygenic control are given in Supplementary Table 2. Candidate lipids proposed for lipidome features were identified through the LIPID MAPS Structure Database (https://www.lipidmaps.org, accessed on 4 June 2022). CAR, Acyl carnitine; DG, Diacylglycerol; FA, Fatty acyl; FOH, Fatty alcohol; LPA, Lipophosphatydicacid; LPC, Lysophosphatidylcholine; MG, Monoradylglycerol; NAE, N-acyl ethanolamine; PA, Phosphatidic acid; PC, Phosphatidylcholine; PE, Phosphatidylethanolamine; PS, Phosphatidylserine; ST, Sterol lipid; TG, Triacylgycerol; WE, Wax ester.
Table 3. Genetic control of lipidomic signals mapped to the genome and proposed lipid assignments. Lipidome data, acquired with a Q Exactive Hybrid Quadrupole-Orbitrap mass spectrometer fitted with a Waters Acquity CSH C18 column, were tested for genetic association with genotyped SNPs in the study group (n = 273). Features were characterized by their retention time (RT) and their mass-to-charge ratio (m/z). Details of SNPs and statistical significance of lipidome features under monogenic control are reported. Full list of genetically mapped LC–MS lipidomic features and details and distinct SNP markers associated with lipid features under polygenic control are given in Supplementary Table 2. Candidate lipids proposed for lipidome features were identified through the LIPID MAPS Structure Database (https://www.lipidmaps.org, accessed on 4 June 2022). CAR, Acyl carnitine; DG, Diacylglycerol; FA, Fatty acyl; FOH, Fatty alcohol; LPA, Lipophosphatydicacid; LPC, Lysophosphatidylcholine; MG, Monoradylglycerol; NAE, N-acyl ethanolamine; PA, Phosphatidic acid; PC, Phosphatidylcholine; PE, Phosphatidylethanolamine; PS, Phosphatidylserine; ST, Sterol lipid; TG, Triacylgycerol; WE, Wax ester.
Positive-Ionization Mode
m/zRTGenetic ControlClosest MarkerClosest GenePutative Lipid
204.12337.098Monogenicrs6992234 (c8)PSD3CAR 2:0 (C9H17NO4)
277.21667.495Polygenic- FA 18:4 (C18H28O2), ST 18:1;O2 (C18H28O2), FA 18:3;O (C18H30O3)
279.23266.953Monogenicrs7759479 (c6)DSTFA 17:4 (C17H26O2)
295.22767.515Polygenic- FA 18:3;O (C18H30O3), FA 18:2;O2 (C18H32O4)
303.23272.294Polygenic- FA 20:5 (C20H30O2), ST 20:2;O2 (C20H30O2), FA 20:4;O (C20H32O3)
305.24774.887Polygenic- FA 20:4 (C20H32O2), ST 20:1;O2 (C20H32O2),FA 20:3;O (C20H34O3)
319.22666.276Polygenic- FA 20:5;O (C20H30O3), FA 20:4;O2 (C20H32O4)
343.22471.225Polygenic- FA 20:4;O (C20H32O3Na)
344.27952.103Monogenicrs6928180 (c6)GRIK2CAR 12:0 (C19H37NO4), FA 19:2;O2 (C19H34O4), FOH 19:3;O3 (C19H34O4)
356.38876.354Polygenic- -
370.29556.145Monogenicrs6928180 (c6)GRIK2CAR 14:1 (C21H39NO4), CAR 14:0;O (C21H41NO5), FA 21:3;O2 (C21H36O4)
377.266110.856Monogenicrs1009439 (c6)RCAN2FA 21:2;O2 (C21H38O4Na), MG 18:2 (C21H38O4Na)
379.282145.907Monogenicrs1009439 (c6)RCAN2FA 21:1;O2 (C21H40O4Na), MG 18:1 (C21H40O4Na), WE 21:1;O2 (C21H40O4Na)
398.32667.497Monogenicrs6928180 (c6)GRIK2-
400.34282.533Monogenicrs6928180 (c6)GRIK2CAR 16:0 (C23H45NO4), FA 23:2;O2 (C23H42O4)
426.35788.672Monogenicrs6928180 (c6)GRIK2CAR 18:1 (C25H47NO4), CAR 18:0;O (C25H49NO5)
429.373309.265Polygenic- ST 29:2;O2 (C29H48O2), ST 29:1;O3 (C29H50O3)
431.352314.575Polygenic- ST 28:2;O3 (C28H46O3), ST 28:1;O4 (C28H48O4)
447.347365.330Polygenic- ST 28:2;O4 (C28H46O4), ST 28:1;O5 (C28H48O5)
448.391309.387Polygenic- -
469.365309.438Polygenic- ST 29:1;O3 (C29H50O3Na)
518.32463.675Polygenic- LPC 18:3 (C26H48NO7P), PC 18:1 (C26H50NO8P)
563.551133.091Polygenic- -
568.34067.238Monogenicrs12997234 (c2)DPP10LPC 22:6 (C30H50NO7P)
590.32167.252Monogenicrs12997234 (c2)DPP10LPC 22:6 (C30H50NO7PNa)
612.556808.044Monogenicrs11855528 (c15)CEMIPDG 34:1 (C37H70O5), DG 35:2 (C37H70O5)
646.03158.383Polygenic- -
662.02562.334Polygenic- -
712.645897.105Monogenicrs2002218 (c3)IQSEC1TG 40:0 (C43H82O6)
738.660898.395Polygenic- TG 42:1 (C45H84O6)
756.553408.519Polygenic- PC 34:3 (C42H78NO8P),PE 37:3 (C42H78NO8P), PS O-36:2 (C42H80NO9P), PA 39:4 (C42H75O8P)
758.560408.446Polygenic- -
758.569457.168Polygenic- PC 34:2 (C42H80NO8P), PC 37:2 (C42H80NO8P), PS O-36:1 (C42H82NO9P), PA 39:3 (C42H77O8P)
766.574442.363Monogenicrs13362253 (c5)MSX2PC O-36:5 (C44H80NO7P), PC 36:3 (C44H82NO8P), PE 39:3 (C44H82NO8P)
780.553373.605Monogenicrs2260930 (c20)SEL1L2PC 36:5 (C44H78NO8P), PE 39:5 (C44H78NO8P), PC 36:4;O (C44H80NO9P), PS O-38:4 (C44H80NO9P), PA 41:6 (C44H75O8P)
784.584560.683Polygenic- PC 36:3 (C44H82NO8P), PE 39:3 (C44H82NO8P), PA 41:4 (C44H79O8P)
792.707921.958Polygenic- TG 46:2 (C49H90O6)
864.764887.193Polygenic- -
876.728841.945Polygenic- -
886.749911.605Polygenic- -
888.764928.842Polygenic- -
890.771929.103Polygenic- -
894.754922.854Polygenic- TG 54:7 (C57H96O6)
912.764912.510Polygenic- -
914.779929.523Polygenic- -
922.785939.142Monogenicrs2292329 (c16)NECAB2TG 56:7 (C59H100O6)
932.8641004.391Monogenicrs11071737 (c15)RAB8BTG 56:2 (C59H110O6)
946.785930.853Polygenic- TG 58:9 (C61H100O6)
948.800946.043Polygenic- TG 58:8 (C61H102O6)
Negative-Ionization Mode
187.00636.489Polygenic- -
271.228113.649Polygenic- FA 16:0;O (C16H32O3)
293.21364.408Polygenic- FA 18:3;O (C18H30O3)
295.22864.394Monogenicrs7760515 (c6)DSTFA 18:2;O (C18H32O3)
303.233129.783Polygenic- ST 20:1;O2 (C20H32O2)
311.22364.059Polygenic- FA 18:2;O2 (C18H32O4), FA 17:2 (C17H30O2), WE 17:2 (C17H30O2), WE 16:2 (C16H28O2), FA 16:2 (C16H28O2)
317.21262.651Monogenicrs7193436 (c16)MVDFA 20:5;O (C20H30O3), ST 19:2;O (C19H28O)
319.22870.158Polygenic- FA 20:4;O (C20H32O3), ST 19:1;O (C19H30O)
321.24371.306Polygenic- FA 20:3;O (C20H34O3), ST 19:0;O (C19H32O)
327.233118.705Polygenic- FA 22:6 (C22H32O2)
343.22865.947Polygenic- FA 22:6;O (C22H32O3), ST 22:3;O3 (C22H32O3), ST 20:3;O (C20H28O)
345.24468.352Polygenic- ST 21:2;O (C21H32O), ST 20:2;O (C20H30O)
409.23680.634Polygenic- LPA 16:0 (C19H39O7P)
433.23668.781Polygenic- LPA 18:2 (C21H39O7P)
437.29160.227Polygenic- ST 24:1;O4 (C24H40O4),FA 23:4;O2 (C23H38O4),FOH 23:5;O3 (C23H38O4),MG 20:4 (C23H38O4),ST 23:1;O4 (C23H38O4)
446.377287.415Polygenic- NAE 24:0 (C26H53NO2), TG 55:5 (C58H102O6)
448.30747.807Polygenic- ST 24:1;O4;G (C26H43NO5)
457.23666.170Polygenic- ST 24:2;O6 (C24H38O6)
591.391200.190Polygenic- ST 27:2;O;Hex (C33H54O6)
605.406223.252Monogenicrs1487842 (c11)SYT9ST 27:2;O;Hex (C33H54O6)
612.33164.327Monogenicrs12997234 (c2)DPP10LPC 22:6 (C30H50NO7P),LPE 24:6 (C29H48NO7P)
804.567435.379Monogenicrs2655474 (c9)ELAVL2PC O-36:3 (C44H84NO7P)
812.582530.577Polygenic- PC O-36:4 (C44H82NO7P), PC O-35:4 (C43H80NO7P), PE O-38:4 (C43H80NO7P)
828.577487.561Polygenic- -
828.577514.160Polygenic- -
Table
Table 4. Significant associations between lipidomic features and clinical and biochemical phenotypes in the study group. Lipidomic features were independently acquired in negative- and positive-ionization modes in serum samples from a study group of 273 individuals. Linear regression was used to compute a P-value statistic for each metabolic feature, which was corrected for multiple testing using the Benjamini-Hochberg method to calculate adjusted p-values. Significant evidence of association was obtained for cardiometabolic disease (CMD), family history (FH) of hypertension, body-mass index (BMI) and total and HDL cholesterol. CMD was assessed by presence of at least three anomalies (diabetes, hypertension, BMI > 30kg/m2, HDL < 40mg/dl). Results from association analysis for all phenotypes that did not reach statistical significance following correction for multiple testing (nominal p-value < 0.05) are shown in Supplementary Table S3. Mass-to-charge ratio (m/z) and retention time (RT) are reported for each lipidome feature. Assignment of lipid candidates for lipidome features was performed using LIPID MAPS (https://www.lipidmaps.org, accessed 1 May 2022). CAR, Acyl carnitine; FA, Fatty acyl; CL, Cardiolipin; NAT, N-acyl amide; PE, Phosphatidylethanolamine; PG, Phosphatidylglycerol; ST, Sterol lipid.
Table 4. Significant associations between lipidomic features and clinical and biochemical phenotypes in the study group. Lipidomic features were independently acquired in negative- and positive-ionization modes in serum samples from a study group of 273 individuals. Linear regression was used to compute a P-value statistic for each metabolic feature, which was corrected for multiple testing using the Benjamini-Hochberg method to calculate adjusted p-values. Significant evidence of association was obtained for cardiometabolic disease (CMD), family history (FH) of hypertension, body-mass index (BMI) and total and HDL cholesterol. CMD was assessed by presence of at least three anomalies (diabetes, hypertension, BMI > 30kg/m2, HDL < 40mg/dl). Results from association analysis for all phenotypes that did not reach statistical significance following correction for multiple testing (nominal p-value < 0.05) are shown in Supplementary Table S3. Mass-to-charge ratio (m/z) and retention time (RT) are reported for each lipidome feature. Assignment of lipid candidates for lipidome features was performed using LIPID MAPS (https://www.lipidmaps.org, accessed 1 May 2022). CAR, Acyl carnitine; FA, Fatty acyl; CL, Cardiolipin; NAT, N-acyl amide; PE, Phosphatidylethanolamine; PG, Phosphatidylglycerol; ST, Sterol lipid.
Ionization Modem/zRTPAdjusted PCorrelationR SquaredAdjusted R SquaredPutative Lipid
CMDNegative317.05948.7456.19 × 10−96.09 × 10−60.1050.1250.115-
Negative319.05648.7597.97 × 10−96.09 × 10−60.0610.1230.113-
Negative386.23759.8456.06 × 10−83.09 × 10−50.0580.1120.102NAT 18:2 (C20H37NO4S)
Negative466.308161.7818.74 × 10−72.74 × 10−40.0590.1020.092CAR 18:3 (C25H43NO4)
Negative465.305162.0101.02 × 10−62.74 × 10−40.0530.1030.093ST 27:1;O;S (C27H46O4S)
Negative497.12248.7071.07 × 10−62.74 × 10−40.1330.0930.083-
Negative231.02148.7307.22 × 10−60.0020.0150.0800.070FA 7:4;O4 (C7H6O6)
Negative233.01848.7598.94 × 10−60.0020.1500.0790.068-
Negative313.239115.0771.44 × 10−50.0020.1270.0840.073-
Negative463.344138.7129.16 × 10−50.0140.0160.0570.046ST 28:1;O5 (C28H48O5),ST 27:1;O3 (C27H46O3),ST 26:1;O3 (C26H44O3)
Negative551.359180.9072.40 × 10−40.0330.1400.0710.061-
Negative591.391200.1902.85 × 10−40.0360.1270.0560.046ST 27:2;O;He × (C33H54O6)
Negative592.394200.0093.79 × 10−40.0430.1240.0550.045PE 25:0 (C30H60NO8P)
Negative607.386200.3033.91 × 10−40.0430.1140.0470.036ST 27:1;O;GlcA (C33H54O7)
FH HypertensionNegative695.511336.9907.62 × 10−60.0120.0290.0930.083-
Negative938.536440.6933.61 × 10−50.0280.1040.0680.058-
BMIPositive774.543527.9851.80 × 10−50.0270.1820.0910.081-
Positive833.588430.1885.81 × 10−50.0370.1740.0700.060PG 40:4 (C46H83O10PLi)
Positive834.591429.7479.24 × 10−50.0370.1690.0680.057Hex 2Cer 32:1;O2 (C44H83NO13)
Positive832.584429.5129.85 × 10−50.0370.1610.0640.053PC 40:7 (C48H82NO8P), PS O-42:6 (C48H84NO9P)
Total CholesterolPositive758.569457.1681.26 × 10−60.002−0.0120.0850.075-
Positive759.572457.3702.35 × 10−60.0020.0220.0840.074CL 76:2 (C85H162O17P2)
HDL CholesterolNegative367.22884.9692.44 × 10−50.0370.0100.0780.068ST 24:5;O3 (C24H32O3)
Positive213.14649.5625.72 × 10−60.0080.0130.0910.081FA 13:4 (C13H18O2Li),WE 13:4 (C13H18O2Li)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brial, F.; Hedjazi, L.; Sonomura, K.; Al Hageh, C.; Zalloua, P.; Matsuda, F.; Gauguier, D. Genetic Architecture of Untargeted Lipidomics in Cardiometabolic-Disease Patients Combines Strong Polygenic Control and Pleiotropy. Metabolites 2022, 12, 596. https://doi.org/10.3390/metabo12070596

AMA Style

Brial F, Hedjazi L, Sonomura K, Al Hageh C, Zalloua P, Matsuda F, Gauguier D. Genetic Architecture of Untargeted Lipidomics in Cardiometabolic-Disease Patients Combines Strong Polygenic Control and Pleiotropy. Metabolites. 2022; 12(7):596. https://doi.org/10.3390/metabo12070596

Chicago/Turabian Style

Brial, Francois, Lyamine Hedjazi, Kazuhiro Sonomura, Cynthia Al Hageh, Pierre Zalloua, Fumihiko Matsuda, and Dominique Gauguier. 2022. "Genetic Architecture of Untargeted Lipidomics in Cardiometabolic-Disease Patients Combines Strong Polygenic Control and Pleiotropy" Metabolites 12, no. 7: 596. https://doi.org/10.3390/metabo12070596

APA Style

Brial, F., Hedjazi, L., Sonomura, K., Al Hageh, C., Zalloua, P., Matsuda, F., & Gauguier, D. (2022). Genetic Architecture of Untargeted Lipidomics in Cardiometabolic-Disease Patients Combines Strong Polygenic Control and Pleiotropy. Metabolites, 12(7), 596. https://doi.org/10.3390/metabo12070596

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop