Vascular lipidomics analysis reveales increased levels of phosphocholine and lysophosphocholine in atherosclerotic mice
Nutrition & Metabolism volume 20, Article number: 1 (2023)
Atherosclerosis (AS) is the major cause of cardiovascular disease, and dyslipidemia is a principal determinant of the initiation and progression of AS. Numerous works have analyzed the lipid signature of blood, but scarce information on the lipidome of vascular tissue is available. This study investigated the lipid profile in the aorta of ApoE−/− mice.
ApoE−/− mice were randomly divided into two groups: (1) the normal diet (ND) group and (2) the high-fat diet (HFD) group. After feeding for 8 weeks, the plasma low-density lipoprotein (LDL), total cholesterol (TC), and triglyceride (TGs) levels were measured. UHPLC-Q Exactive plus MS was used to assess the lipid profile using both positive and negative ionization modes.
LDL and TC levels were significantly increased in HFD mice, and lipid deposition, plaque area and collagen fiber levels were increased in HFD group. In addition, a total of 131 differential lipids were characterized, including 57 lipids with levels that were increased in the HFD group and 74 with levels that were decreased. Further analysis revealed that the levels of several differentially expressed phosphocholines (PCs) and lysophosphocholines (LPCs) were significantly increased. These PCs included PC (38:3), PC (36:4), PC (36:3), PC (36:2), PC (36:1), PC (34:1e), PC (34:1), PC (32:1), PC (18:0/18:1), and PC (38:5), and the LPCs included LPC (18:1), LPC (18:0) and LPC (16:0).
Our findings indicate the presence of a comprehensive lipid profile in the vascular tissue of atherosclerotic mice, particularly involving PC and LPC, which exhibited significantly increased levels in AS.
Cardiovascular and cerebrovascular diseases caused by atherosclerosis (AS) remain the leading cause of morbidity and mortality worldwide [1, 2]. AS is characterized by lipid disturbance, which leads to the formation of fibrofatty lesions in the arterial wall, subsequent hardening and narrowing of the arteries, and eventually the formation of plaques in the vessel wall [3, 4]. Lipid disturbance, vascular dysfunction, smooth muscle cell proliferation, oxidative stress and inflammation exacerbate the progression of atherosclerosis . A large number of studies show that lipid disorder is a principal determinant of the initiation and progression of AS . Endogenous lipids are the major constituents of cell membranes, which are extensively involved in the structural compartmentalization of signal transduction. Moreover, these lipids act as potential pathophysiological mediators and participate in regulating several intracellular functions, such as proliferation, apoptosis, oxidative stress and inflammation . Atherosclerosis is closely linked to lipid overload, including increased plasma levels of triglycerides (TG), low-density lipoprotein (LDL), intermediate-density lipoprotein (IDL) and VLDL [7, 8]. For example, LDL, a particle that carries cholesterol through the blood, is not only the initiator but also an accelerator of AS [9, 10]. Lipoprotein(a), a low-density lipoprotein (LDL) cholesterol-like particle that binds to apolipoprotein (A), has a potential causal relationship with atherosclerotic cardiovascular disease .
Lipidomics is a powerful analytical tool for analyzing the lipid composition in biological samples under physiological and pathological conditions to elucidate disease pathogenesis, identify biomarkers or study the effect and mechanism of drugs [12, 13]. Liquid chromatography coupled to MS (LC‒MS) is undoubtedly considered the current gold standard for the quantification of lipids . Lipidomics has been widely used in screening and mechanistic research on lipid biomarkers in atherosclerotic diseases in recent years. For example, the plasma lipid levels in subjects with TG/HDL-C levels ≥ 3.287 were found to be significantly different from those in subjects with TG/HDL-C levels < 1.391. The downregulated lipid including sphingomyelins (SM) (36:2), SM (38:2), SM (40:2), SM (41:2), SM (41:1), SM (42:4), SM (42:3), SM (42:1), electrophoresis (CE) (16:0), CE (18:0), CE (18:1), PC (38:2), and PC (36:2e)/PC (36:1p), and PE (38:1), to the contrary, all TG species and phosphatidylserine (PS) (36:1) were increased . In patients with large artery atherosclerotic cerebrovascular disease, 276 sphingolipids were detected in the plasma via lipidomic analysis, including SM (n = 162), glycosphingolipids (n = 72), Cer (n = 39), and ceramide phosphates (n = 3) . Moreover, Wang et al.  identified nine types of 732 lipid components in the plasma of ApoE−/− mice, including 11 lysophosphatidylethanolamine (LPE), 12 SM, 15 LPC, 27 cholesteryl esters (CE), 35 ceramides (Cer), 38 diacylglycerol (DAG), 60 phosphatidylethanolamine (PE), 64 PC, and 470 triglyceride (TAG). Therefore, further comprehensive exploration of the lipid characteristics related to atherosclerosis, including in serum and vascular tissue, is conducive to enabling an in-depth understanding of the pathogenesis of atherosclerosis and is also expected to provide new insights for the diagnosis, prevention and cure of atherosclerosis.
In summary, omics approaches are used to identify lipid profiles in plasma, serum, and plaques to establish a disease prediction model, elucidate disease pathogenesis, improve diagnostic efficiency and provide drug targets. To date, due to the limited sample size, although Jihan Talib et al. have developed a new Multi-ABLE method that can be used to assess small tissue samples (< 5 mg) and have performed aortic lipidomics in ApoE−/− mice , there has been no literature describing the use of classic LC‒MS techniques to assess the profile of vascular lipids in ApoE−/− mice. Therefore, we used classic LC‒MS/MS lipidomics technology to provide a basis for elucidating the lipid profile of atherosclerotic thoracic aortas in ApoE−/− mice.
Materials and methods
Reagents and materials
The chemicals methanol (No. A452-4), acetonitrile (No. A998-4), formic acid (No. A117-50), and isopropanol (No. A451-4) were purchased from ThermoFisher (ThermoFisher, Waltham, MA, USA), ammonium acetate (10001218) was purchased from Sinopharm (Sinopharm, Beijing, China), Lyso PC17:0 (No. 855676P) was purchased from Avanti (Avanti, Alabama, USA), chloroform (No. G75915B) was purchased from Greagent (Greagent, Shanghai). Distilled water was purchased from Watsons (Watsons, Guangzhou, China). All aqueous solutions were prepared using purified water at a Milli-Q grade.
Animal model establishment and specimen acquisition
This study was approved by the Chengdu Women’s and Children’s Central Hospital Ethics Committee. ApoE−/− mice (4 weeks old, male) were purchased from Jiangsu Jicui Yaokang Biotechnology Co., Ltd. Animals were kept in a specific pathogen-free (SPF) environment. The feeding conditions included a light–dark cycle for 12 h, temperature 22–25 °C, and a relative humidity of 60%. The mice were divided into two groups: (1) mice were fed a normal diet, and (2) mice were fed a high-fat diet (Yangzhou Promoter Biotechnology Co., Ltd., D12108C). After feeding for 8 weeks, the mice were fasted overnight and anesthetized with sodium pentobarbital solution (1%, 10 mL/kg) by intraperitoneal injection. The plasma was separated from the blood, which was collected from the abdominal aorta, after standing for 2 h and centrifugation at 8000 rpm for 8 min. The plasma was transferred to an EP tube and preserved at − 80 °C refrigerator before being used for analysis. The thoracic aorta was peeled off and preserved at − 80 °C refrigerator.
The measurement of serum TC, TG and LDL levels
The cryopreserved plasma (100 μL) obtained following fasting was thawed and then centrifuged. The biochemical markers associated with lipid metabolism, including the levels of LDL, TC, and TG, were detected using a commercial kit according to the manufacturer’s instructions (Changchun Huili Biotech Co., Ltd.).
The heart, including the aortic root, was fixed in paraformaldehyde. Then, the tissues were embedded in OCT compound and sliced into sections (7 μm) . Then, HE and Oil Red O staining were performed to observe the plaque area and lipid deposition, respectively. Furthermore, Masson’s trichrome staining was performed to observe the collagen area.
All images were photographed using a panoramic MIDI scanning microscope (3DHISTECH, Hungary) and then quantified using Image-Pro Plus 6.0 software. The plaque area was measured as the total area of plaque. Lipid deposition is presented as the percentage of positive cells in total plaque. The collagen area was measured as the collagen-positive area.
As previously described, the expression of CD68 in the aortic root was assessed by immunostaining . Briefly, the frozen sections were permeabilized with Triton X-100 (0.2%, 15 min) at room temperature. After incubation with rabbit anti-CD68 (1:100, Servicebio) at 4 °C overnight, the sections were washed, a secondary antibody was added, and the samples were incubated. Finally, DAPI was added to stain the nuclei. Images were acquired using a panoramic MIDI scanning microscope (3DHISTECH, Hungary) and then quantified using Image-Pro Plus 6.0 software.
Due to the limitation of sample weight, the thoracic aorta vessels of every two mice were combined in the analysis. After quickly freezing with liquid nitrogen, the samples were placed in and stored at − 80 °C refrigerator until analysis. Lipids were extracted as follows: (1) 20 μL Lyso PC-17:0 (0.1 mg/mL, dissolved in methanol) was added to the samples as internal standard to reflect the repeatability of the extraction. And then mixed the samples with methanol: water (1/1, vol/vol, 300 μL). Subsequently, the samples were subjected to grinding with two small steel balls (60 Hz, 2 min) in grinder (Tissuelyser-48, Shanghai) and then added 300 μL chloroform for ultrasonication for 10 min, and stored at − 20 °C for 20 min. After centrifugation for 10 min, 200 μL of the chloroform layer was transferred into a new centrifuge tube. (2) The residual samples were added to chloroform:methanol (2/1, vol/vol) and then vortexed for 30 s. The samples were ultrasonicated for 10 min, placed at − 20 °C for 20 min, and centrifuged for 10 min, and then taken 150 μL of the chloroform layer into the centrifuge tube. (3) The two chloroform layers were combined and then dried under a nitrogen stream. The samples were re-dissolved in 400 μL isopropanol:methanol (1/1, vol/vol), vortexed for 30 s and sonicated for 3 min. Before detection, the solution was filtered through a 0.22 μm organic phase pinhole filter. All sonication was performed in an ice-water bath, and centrifugation was performed at 13,000 rpm at 4 °C. QC samples were prepared by mixing aliquots of all samples, and then were injected at regular intervals (every 6 samples) into the samples throughout the analysis procedure to evaluate the stabilityof the equipment. All chemicals and solvents were analytical or HPLC grade.
Lipidomic profiling using LC‒MS
Lipid profiling was performed using an UPLC system (Waters, Japan) coupled to a Q-Exactive plus mass spectrometer equipped with a heated electrospray ionization (ESI) source (Thermo Fisher Scientific, Waltham, MA, USA). In detail, the samples were separated using an ACQUITY UPLC BEH C8 column (100 mm × 2.1 mm, 1.7 um, Waters) and detected using ESI positive and ESI negative ion modes. During the analysis, the samples were kept at 4 °C. The gradient elution system consisted of two mobile phases, (A) acetonitrile: water (6:4, v:v, containing 10 mmol ammonium formate, 0.1% formic acid) and (B) isopropanol: acetonitrile (9:1, v:v, containing 10 mmol ammonium formate, 0.1% formic acid). The flow rate was 0.26 mL/min, and the column temperature was 55 °C. The injection volumes were 2 μL and 4 μL in positive and negative modes, respectively. The separation was achieved using the following gradient: 0 min, 32% B; 1.5 min, 32% B; 15.5 min, 85% B; 15.6 min, 97% B; 18.0 min, 97% B; 18.1 min, 32% B; 20.0 min, 32% B.
The mass spectrometer was operated in both positive electrospray ionization (ESI+) mode and negative electrospray ionization (ESI−) mode. The acquisition parameters: the spray voltage value of positive ions is 3.5 kv, and 3.0 kv negative ions; Capillary Temp 300 °C; Mass range (m/z): 150–1500; Full ms resolution: 70,000; MS/MS resolution: 17,500.
Ingenuity pathway analysis (IPA) is a widely used bioinformatics software for predicting the correlations among components, targets, pathways, and diseases [21, 22]. The HMDB number of differential lipids was downloaded from the HMDB database (https://hmdb.ca/) . In canonical pathway analysis and analyses of disease and function, we took − log (P value) > 2 as the threshold and then defined Z score > 2 and Z score < − 2 as the threshold of significant activation and inhibition, respectively. The consistency scores were calculated by assessing regulatory effects and molecular networks, and the higher the consistency score was, the more accurate the regulatory effect analysis results. In addition, we used a P value of overlap < 0.05 as the inclusion criteria for upstream regulators.
The original data obtained from Q Exactive plus (LC‒MS/MS) analysis were processed by Lipid Search software. The molecular structure of lipids and the additive mode of its positive and negative ions were confirmed according to the parent ions and multistage mass spectrometry data of each individual sample. To sort the original data matrix, results were aligned according to a certain retention time range. The peak signals in each sample were normalized, by converting its signal intensity into the relative intensity in the spectrum, and multiplying it by 100,000. Lipids which the relative standard deviation (RSD) of QC greater than 30% were discarded. These highly variable ions should be rejected as unacceptable for the analysis. Following the extraction of data, peaks with missing values (ion intensity = 0) in more than 50% of groups were removed and replaced the zero value with half of the minimum value. Positive and negative ion data were combined into a data matrix. Principal component analysis (PCA) was carried out to monitor the overall sample distribution and the stability of the whole analysis process. Orthogonal partial least squares-discriminant analysis (OPLS-DA) and partial least squares-discriminant analysis (PLS-DA) were used to authenticate the results regarding the differential lipids. Variable importance of projection (VIP) value was used to assess the variable importance. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database was used to analyze metabolic pathway enrichment. The statistical graphs were generated by GraphPad Prism 8 software, and a two-tailed Student’s T test was used for statistical analysis. All data are presented as the means ± standard errors of the means, SEM). Differential lipids were selected with VIP values > 1.0 and p values < 0.05.
Effects of a high-fat diet on the lipid profile in ApoE −/− mice
We used male ApoE−/− mice fed a high-fat diet to establish an animal model of atherosclerosis and then collected blood vessel tissue to extract lipids for mass spectrometry analysis. The study design is shown in Fig. 1A. In brief, ApoE−/− mice were fed a high-fat diet (HFD) for 8 weeks to establish an atherosclerosis model, while ApoE−/− mice fed a normal diet (ND) were used as a control. The lipid species in thoracic aorta were assessed via HPLC‒MS/MS. The body weight was not significantly different between the ND group and the HFD group (Fig. 1B). The plasma levels of TC, LDL-C, and TG were detected to investigate the effects of HFD on total lipids. The levels of TC and LDL-C in the HFD group were significantly higher than those in the ND group, while no significant differences in TG levels were observed; Fig. 1C–E illustrates the details.
Assessment of atherosclerotic plaque formation in ApoE −/− mice
To verify that the model of atherosclerotic mice was successful, HE and Oil Red O staining were performed in frozen sections of the aortic root of ApoE−/− mice. The results of HE staining showed that 8 weeks of HFD treatment significantly increased plaque formation (Fig. 2A, B). Oil Red O staining showed that lipid deposition was also significantly increased in HFD-treated ApoE−/− mice (Fig. 2C, D). Masson's trichrome staining showed that the collagen content was significantly increased in HFD-treated ApoE−/− mice (Fig. 2E, F). In addition, the immunofluorescence staining of CD68 expression revealed that free monocyte/macrophage levels in the aortic root were increased (Fig. 2G, H). The above results revealed that ApoE−/− mice treated with HFD did indeed undergo the pathological changes associated with atherosclerosis, which suggested the success of animal model establishment.
UHPLC-Q Exactive MS-based lipidomics analysis
UHPLC-Q Exactive Plus-MS was used to lipidomic analysis of the thoracic aorta samples. The molecular structure of lipids and the additive mode of the positive and negative ions were identified according to the parent ions and multistage mass spectrometry data obtained from each individual sample. The base peak chromatograms (BPCs) of positive ions and negative ions in the ND group and HFD group are shown in Additional file 1: Fig. S1, respectively.
To evaluate the stability of the equipment, the lipid intensity distribution of the QC samples is shown in Fig. 3A, where the Y coordinate is the log10 value of the mass spectrum intensity. PCA was used first to monitor the overall distribution of the samples and the stability of the entire analysis process, as shown in (Additional file 2: Fig. S2). The distance between the two coordinate points on the score chart was far on the score map, indicating that there was a significant difference between the two samples. PLS-DA and OPLS-DA were used to identify differential lipids between groups (Fig. 3B, Additional file 2: Fig. S2). A random ranking method, the response ranking test, was used to evaluate the accuracy of the OPLS model. The value (green) of all Q2 on the left was lower than the initial value on the right, which shows that the model was effective (Fig. 3C). An S-plot showed that all the points were distributed in the first and third quadrants, and the closer the lipid molecules were to the lower left or upper right of the plot, the more significant the differences were (Fig. 3D).
Analysis of lipidomic profiling and screening of differential lipids
Based on the PLS-DA model, a total of 131 differential lipids were identified according to the conditions of VIP > 1 and P < 0.05. A total of 57 lipids exhibited increased levels, and 74 exhibited decreased levels in the HFD group (Fig. 4A). The details of the differential lipids are shown in Table 1. As shown in Fig. 4B, C, the expression levels of the 131 differential lipids were used to draw a heatmap and volcano plots for comparisons between the HFD and ND groups. Unlike in the ND group, the top 20 lipid components included PE (18:0/20:4), PE (18:1/20:4), TG (18:1/18:1/18:1), TG (18:1/18:2/18:2), PS (37:1), PE (16:0p/22:6), TG (18:0/18:1/18:1), TG (16:0/18:1/18:1), TG (16:1/16:1/18:2), TG (20:0/18:1/18:1), PC (34:1), TG (18:0/18:0/18:1), PA (24:2/20:4), PC (36:1), LPE (20:4), PA (26:2/20:3), PE (16:0p/20:4), TG (16:1/14:0/18:2), TG (16:0/18:2/18:2), and LPE (22:6). The details of the differential lipids that were analyzed in the top 20 significantly different metabolites sorted by VIP are shown in Table 2. We performed Z score analysis on the top 20 significantly different lipids sorted by VIP and transformed the data of different magnitudes into a unified measure of Z score for comparison (Fig. 4D). To study the degree of correlation between lipid levels, we performed correlation analysis using the top 20 significantly different lipids sorted by VIP and used the Pearson correlation coefficient to measure the degree of linear correlation between lipid levels. A positive correlation is indicated in red, while a negative correlation is indicated in blue. The larger the dot is, the larger the correlation coefficient between the two variables (Fig. 4E). In addition, we selected the relationship pair with p value < 0.05 and correlation > 0.95 to draw a correlation network diagram (Fig. 4F).
Metabolic pathway enrichment analysis
KEGG pathway enrichment analysis (http://www.kegg.jp/orhttp://www.genome.jp/kegg/) was used to predict metabolic pathways based on the lipids with increased and decreased levels. As shown in Fig. 5A, the lipids with increased levels were mainly associated with glycerophospholipid metabolism, choline metabolism, retrograde endocannabinoid signaling, autophagy-other, systemic lupus erythematosus, autophagy-animal and cholesterol metabolism. Interestingly, the lipids with decreased levels were mainly associated with fat digestion and absorption, regulation of lipolysis in adipocytes, retrograde endocannabinoid signaling, insulin resistance, thermogenesis, glycerophospholipid metabolism and Th1 and Th2 cells differentiation (Fig. 5B). These results indicate that lipid metabolism imbalance, autophagy and inflammatory responses may be involved in the mechanism underlying atherosclerosis. The above results showed that glycerophospholipid metabolism was significantly enriched in both the activated and inhibited KEGG pathways, which attracted our attention. Further analysis of individual glycerophospholipid pathways revealed that some lipids, including phosphatidyl-L-serine, phosphatidyl-ethanolamine, phosphatidylcholine (Lecithin) and 1-acyl-sn-glycero-3-phosphocholine, exhibited increased levels, which may play a critical function in the occurrence or progression of AS (Fig. 5C). Moreover, IPA network pathway analysis showed that differentially expressed lipids were closely related to NF-κB signaling, PI3K/AKT signaling and JAK/STAT signaling (Additional file 3: Fig. S3).
Lipid composition analysis
The types of the 131 lipid components are shown in Additional file 4: Fig. S4, including TG (n = 53), PE (n = 24), (PC (n = 10), PS (n = 10), SM (n = 7), phosphatidic acid (PA, n = 7), and other (n = 20) lipids. In addition, we showed changes in the expression levels of various lipids, including TG, PE, PC, PS, PA, SM and other lipids, in the ND and HFD groups in the form of a chromatic diagram (Additional file 4: Fig. S4). Interestingly, the levels of all lipids in the PC and LPC classes were significantly increased in the HDF group, which may play a nonnegligible effect in the onset and process of atherosclerosis.
As reported, PCs [24, 25] and LPCs  are closely related to atherosclerosis. To explore the role of PC and LPC in atherosclerosis, we further analyzed the expression levels of 10 PCs and 3 LPCs, which exhibited significantly increased levels in the HFD group. The levels of PC, including PC (38:3), PC (36:4), PC (36:3), PC (36:2), PC (36:1), PC (34:1e), PC (34:1), PC (32:1), PC (18:0/18:1), and PC (38:5), are shown in Fig. 6A–J, and those of LPC, including LPC (18:1), LPC(18:0) and LPC (16:0), are shown in Fig. 6K–M.
Dyslipidemia is a major risk factor for atherosclerosis and atherosclerosis-related cardiovascular disease . A consensus has been reached on the conclusion that elevated plasma LDL-cholesterol levels promote AS. LDL lipoproteins penetrate damaged endothelial cells, accumulate in the intima after oxidation, and subsequently recruit monocytes to endothelial cells, causing monocytes to differentiate into macrophages and form foam cells, initiating and accelerating the process of atherosclerosis . Unlike cholesterol, triglycerides are also associated with subclinical atherosclerosis and vascular inflammation . In this study, the lipid profile of the blood vessels in AS mice was investigated through lipidomic technology, and 131 differential lipids were identified in the thoracic aorta of ND- and HFD-treated ApoE−/− mice, including TG, PE, and PC. Unlike in the ND group, the top 20 lipid components included PE (18:0/20:4), PE (18:1/20:4), TG (18:1/18:1/18:1), TG (18:1/18:2/18:2), PS (37:1), PE (16:0p/22:6), TG (18:0/18:1/18:1), TG (16:0/18:1/18:1), TG (16:1/16:1/18:2), TG (20:0/18:1/18:1), PC (34:1), TG (18:0/18:0/18:1), PA (24:2/20:4), PC (36:1), LPE (20:4), PA (26:2/20:3), PE (16:0p/20:4), TG (16:1/14:0/18:2), TG (16:0/18:2/18:2), and LPE (22:6) in the AS group.
Recently, omics approaches have been used to improve the understanding of the pathogenesis of AS. For example, Wang et al.  used ultraperformance liquid chromatography coupled with quadrupole time-of-flight high-definition mass spectrometry (UPLC-Q-TOF/MS) to identify 17 lipid components closely related to the occurrence and development of AS in the serum of ApoE−/− mice. The lipids included DAG (14:0/18:2), PC (18:0/20:1), LPC (18:2), LPC (20:2), PC (18:2/22:6), PE (18:0/18:2), SM (D18:1/16:0), LPE (20:1), TAG48:1-FA12:0, TAG48:1-FA14:0, TAG49:2-FA17:0, TAG52:6-FA16:0, TAG53:3FA16:0, TAG53:3-FA17:0, TAG54:6-FA20:5, TAG54:2-FA16:0, and TAG54:0-FA18:0. Jihan et al.  developed a new Multi-ABLE method to analyze small tissue samples (< 5 mg) and determined the aortic lipid profile in ApoE−/− mice. The results showed that 52 lipid molecules exhibited increased abundance in lesion-containing tissue and included SM (n = 7), phospholipids (n = 24) and CL (n = 3). The results regarding the lipids identified by lipidomics in thoracic aorta of ApoE−/− mice in this article and those in previous studies are inconsistent, which suggested the data from vascular lipidomics should be further investigated to provide an in-depth understanding of the complex pathogenesis of AS.
Recent studies have shown that glycerol phospholipid metabolism plays a critical role in the pathogenesis of AS, and disordered glycerol phospholipid metabolism directly affects the process of AS [29, 30]. In our study, the KEGG enrichment analysis of the lipids that were identified to exhibit increased and decreased levels revealed enrichment in glycerol phospholipid metabolism, indicating that glycerol phospholipid metabolism may be an initiating factor; these results are conducive to further elucidating the pathogenesis of AS. As reported, the major classes of glycosphingolipids, comprising phosphatidylcholine (PC), PE, and PS, are closely related to atherosclerosis. In all nucleated mammalian cells, PC and PE are the most abundant phospholipids in cell membranes [31, 32]. PC is mainly synthesized by the CDP-choline pathway . In an additional pathway, under the condition of S-adenosylmethionine (AdoMet) as a methyl donor, PE can synthesize PC through three sequential methylation reactions . Impaired phosphatidylcholine biosynthesis has been reported to reduce atherosclerosis in ApoE−/− mice . Exteriorized PS and PE residues on the cell membrane in apoptotic cells trigger rapid phagocytosis and contribute to plaque vulnerability [35,36,37]. Among patients with CAD, SMs are directly related to serum levels of total cholesterol and lipoproteins, which suggests that sphingolipids play a critical role in modulating lipid serum levels. The results above highlight the importance of establishing new goals in the glycerol metabolic pathway for antiatherogenic therapies .
In addition, we found that the levels of 10 PCs from HFD-treated ApoE−/− mice were significantly higher than those in the ND group, including PC (38:3), PC (36:4), PC (36:3), PC (36:2), PC (36:1), PC (34:1e), PC (34:1), PC (32:1), PC (18:0/18:1), and PC (38:5). Phosphatidylcholine (PC) is a subclass of glycerophospholipids with choline as the head group. The physiological functions of PC include membrane integrity/stability, cell proliferation and survival, fuel and energy storage/source, glycerophospholipid metabolism and cell signaling [31, 39]. It has been reported that the levels of PCs (16:2/22:6, 16:2/20:5, 15:0/18:2) were decreased in stroke patients . In contrast, PC (16:0/20:4) levels were significantly higher in the serum of patients with severe calcific coronary artery disease than in patients with noncalcific coronary artery disease, while PC(18:2/18:2) and PC(36:3) levels had the opposite trend .
Lysophosphatidylcholine (LPC) is the main component of oxidized low-density lipoprotein (oxLDL). In cells, saturated LPC (16:0) and PC can be interconverted. Specifically, PC can be cleaved by phospholipase A2 (PLA2) to form LPC [42, 43], and LPC can be converted to PC by lysophosphatidylcholine acyltransferase (LPCAT) . In addition, lecithin-cholesterol acyltransferase (LCAT) can transfer fatty acids to free cholesterol to form LPC . LPC damages endothelial cells, induces smooth muscle cell proliferation and lymphocyte and macrophage migration, increases inflammatory mediator production and oxidative stress, and induces cell apoptosis and necrosis, which accelerate the progression of atherosclerosis and atherosclerosis-related cardiovascular diseases [26, 45, 46]. In patients with hypertension, LPCs can be considered potential predictors for further evaluation and validation of early vascular aging (EVA) . As reported, the levels of lysoPC 16:0, 18:0, and 18:1 in human plaques are markedly related to inflammatory mediators, such as monocyte chemoattractant protein-1 (MCP1), interleukin-1β (IL-1β), interleukin-6 (IL-6), and macrophage inflammatory protein-1β (MIP-1β) , which suggests that lysoPCs play a key role in plaque inflammation and stability. In our study, we also found that the levels of these three lipids were elevated in the vascular tissue of mice with atherosclerosis, which was consistent with the above research report. This finding indicates that LPC (16:0), LPC (18:0) and LPC (18:1) may participate in the process of atherosclerosis and be closely related to vascular injury.
Although we revealed the lipidomic profiling of vascular tissue of atherosclerotic mice by LC‒MS/MS technology and identified some differentially expressed lipid molecules, there are still some limitations. First, due to the continuous development of LC‒MS/MS technology, lipidomics methods are widely used to comprehensively analyze lipids and their interactions with other molecules. However, it is not clear whether all lipids in some biological samples can be comprehensively detected and accurately identified due to the complex composition of lipids. It is also difficult to conduct accurate localization analysis on the detected lipids. Therefore, the development of new instruments and new bioinformatics schemes, such as ion mobility mass spectrometry (IM-MS) , NanoLC‒MS , AFADESI-MSI and spatial metabolomics analysis , may bring new insights into the application of metabolomics and lipidomics in the future. Second, our study suggested that PCs and LPCs may have a strong association with atherosclerosis, but the causal relationship between them needs to be further clarified. In addition, we observed an interesting phenomenon. While confirming the successful establishment of the atherosclerosis model, there was no significant difference in the TG levels in the serum between mice in the HFD group and those in the ND group. More interestingly, our lipidomics spectrum analysis also revealed that the levels of TG lipid molecules were not increased in the HFD group, and the levels were in fact lower. It is possible that the specificity of TG as a marker of AS needs further reference.
Our findings indicate that in the comprehensive lipid profile in the vascular tissue of atherosclerotic mice, a total of 131 differential lipids were identified, including 57 lipids with levels that were increased in the HFD group and 74 with levels that were decreased. Further analysis revealed that the identified differentially expressed PCs and LPCs exhibited significantly increased levels. We particularly analyzed the details of PCs and LPCs, which exhibited significantly increased levels in AS. Therefore, targeting PCs, LPCs and lipid metabolism might be a novel and effective therapeutic strategy for atherosclerosis.
Availability of data and materials
- UHPLC-Q Exactive MS:
Ultrahigh-performance liquid chromatography plus Q-Exactive mass spectrometry
Liquid chromatography coupled to MS
Ingenuity pathway analysis
Principal component analysis
Orthogonal partial least squares-discriminant analysis
Partial least squares-discriminant analysis
Variable importance of projection
Kyoto Encyclopedia of Genes and Genomes
Base peak chromatograms
Oxidized low-density lipoprotein
Early vascular aging
Monocyte chemoattractant protein-1
Macrophage inflammatory protein-1β
Libby P, Buring JE, Badimon L, Hansson GK, Deanfield J, Bittencourt MS, Tokgozoglu L, Lewis EF. Atherosclerosis. Nat Rev Dis Primers. 2019;5:56.
World Health Organization. Cardiovascular diseases (CVDs). https://www.who.int/news-room/fact-sheets/detail/cardiovascular-diseases-(cvds) (2020). Accessed 27 Nov 2020.
Gimbrone MA Jr, Garcia-Cardena G. Endothelial cell dysfunction and the pathobiology of atherosclerosis. Circ Res. 2016;118:620–36.
Libby P, Kobold S. Inflammation: a common contributor to cancer, aging, and cardiovascular diseases-expanding the concept of cardio-oncology. Cardiovasc Res. 2019;115:824–9.
Bjorkegren JLM, Lusis AJ. Atherosclerosis: recent developments. Cell. 2022;185:1630–45.
Casares D, Escriba PV, Rossello CA. Membrane lipid composition: effect on membrane and organelle structure, function and compartmentalization and therapeutic avenues. Int J Mol Sci. 2019;20:20167.
Bobryshev YV, Nikiforov NG, Elizova NV, Orekhov AN. Macrophages and their contribution to the development of atherosclerosis. Results Probl Cell Differ. 2017;62:273–98.
Mineo C. Lipoprotein receptor signalling in atherosclerosis. Cardiovasc Res. 2020;116:1254–74.
Mehta A, Shapiro MD. Apolipoproteins in vascular biology and atherosclerotic disease. Nat Rev Cardiol. 2022;19:168–79.
Ference BA, Ginsberg HN, Graham I, Ray KK, Packard CJ, Bruckert E, Hegele RA, Krauss RM, Raal FJ, Schunkert H, et al. Low-density lipoproteins cause atherosclerotic cardiovascular disease. 1. Evidence from genetic, epidemiologic, and clinical studies. A consensus statement from the European Atherosclerosis Society Consensus Panel. Eur Heart J. 2017;38:2459–72.
Duarte Lau F, Giugliano RP. Lipoprotein(a) and its significance in cardiovascular disease: a review. JAMA Cardiol. 2022;7:760–9.
Han X, Gross RW. The foundations and development of lipidomics. J Lipid Res. 2022;63:100164.
Ahluwalia K, Ebright B, Chow K, Dave P, Mead A, Poblete R, Louie SG, Asante I. Lipidomics in understanding pathophysiology and pharmacologic effects in inflammatory diseases: considerations for drug development. Metabolites. 2022;12:333.
Han X, Gross RW. Electrospray ionization mass spectroscopic analysis of human erythrocyte plasma membrane phospholipids. Proc Natl Acad Sci U S A. 1994;91:10635–9.
Di Giorgi N, Michelucci E, Smit JM, Scholte A, El Mahdiui M, Knuuti J, Buechel RR, Teresinska A, Pizzi MN, Roque A, et al. A specific plasma lipid signature associated with high triglycerides and low HDL cholesterol identifies residual CAD risk in patients with chronic coronary syndrome. Atherosclerosis. 2021;339:1–11.
You Q, Peng Q, Yu Z, Jin H, Zhang J, Sun W, Huang Y. Plasma lipidomic analysis of sphingolipids in patients with large artery atherosclerosis cerebrovascular disease and cerebral small vessel disease. 2020. Biosci Rep. https://doi.org/10.1042/BSR20201519.
Wang H, Zhang L, Zhang X, Song J, Guo Q, Zhang X, Bai D. Prediction model for different progressions of Atherosclerosis in ApoE(-/-) mice based on lipidomics. J Pharm Biomed Anal. 2022;214:114734.
Talib J, Hains PG, Tumanov S, Hodson MP, Robinson PJ, Stocker R. Barocycler-based concurrent multiomics method to assess molecular changes associated with atherosclerosis using small amounts of arterial tissue from a single mouse. Anal Chem. 2019;91:12670–9.
Xing SS, Yang J, Li WJ, Li J, Chen L, Yang YT, Lei X, Li J, Wang K, Liu X. Salidroside decreases atherosclerosis plaque formation via inhibiting endothelial cell pyroptosis. Inflammation. 2020;43:433–40.
Xu S, Chen H, Ni H, Dai Q. Targeting HDAC6 attenuates nicotine-induced macrophage pyroptosis via NF-kappaB/NLRP3 pathway. Atherosclerosis. 2021;317:1–9.
Shuey MM, Xiang RR, Moss ME, Carvajal BV, Wang Y, Camarda N, Fabbri D, Rahman P, Ramsey J, Stepanian A, et al. Systems approach to integrating preclinical apolipoprotein E-Knockout investigations reveals novel etiologic pathways and master atherosclerosis network in humans. Arterioscler Thromb Vasc Biol. 2022;42:35–48.
Chen C, Cui S, Li W, Jin H, Fan J, Sun Y, Cui Z. Ingenuity pathway analysis of human facet joint tissues: insight into facet joint osteoarthritis. Exp Ther Med. 2020;19:2997–3008.
Wishart DS, Guo A, Oler E, Wang F, Anjum A, Peters H, Dizon R, Sayeeda Z, Tian S, Lee BL, et al. HMDB 5.0: the Human Metabolome Database for 2022. Nucleic Acids Res. 2022;50:D622–31.
Cole LK, Dolinsky VW, Dyck JR, Vance DE. Impaired phosphatidylcholine biosynthesis reduces atherosclerosis and prevents lipotoxic cardiac dysfunction in ApoE-/- Mice. Circ Res. 2011;108:686–94.
Aldana-Hernandez P, Azarcoya-Barrera J, van der Veen JN, Leonard KA, Zhao YY, Nelson R, Goruk S, Field CJ, Curtis JM, Richard C, Jacobs RL. Dietary phosphatidylcholine supplementation reduces atherosclerosis in Ldlr(-/-) male mice(2). J Nutr Biochem. 2021;92:108617.
Liu P, Zhu W, Chen C, Yan B, Zhu L, Chen X, Peng C. The mechanisms of lysophosphatidylcholine in the development of diseases. Life Sci. 2020;247:117443.
Pirillo A, Casula M, Olmastroni E, Norata GD, Catapano AL. Global epidemiology of dyslipidaemias. Nat Rev Cardiol. 2021;18:689–700.
Raposeiras-Roubin S, Rossello X, Oliva B, Fernandez-Friera L, Mendiguren JM, Andres V, Bueno H, Sanz J, Martinez de Vega V, Abu-Assi E, et al. Triglycerides and residual atherosclerotic risk. J Am Coll Cardiol. 2021;77:3031–41.
Su D, Liao L, Zeng Q, Liao Z, Liu Y, Jin C, Zhu G, Chen C, Yang M, Ai Z, Song Y. Study on the new anti-atherosclerosis activity of different Herba patriniae through down-regulating lysophosphatidylcholine of the glycerophospholipid metabolism pathway. Phytomedicine. 2022;94:153833.
Dang VT, Zhong LH, Huang A, Deng A, Werstuck GH. Glycosphingolipids promote pro-atherogenic pathways in the pathogenesis of hyperglycemia-induced accelerated atherosclerosis. Metabolomics. 2018;14:92.
Furse S, de Kroon AI. Phosphatidylcholine’s functions beyond that of a membrane brick. Mol Membr Biol. 2015;32:117–9.
van Meer G, Voelker DR, Feigenson GW. Membrane lipids: where they are and how they behave. Nat Rev Mol Cell Biol. 2008;9:112–24.
Kennedy EP, Weiss SB. The function of cytidine coenzymes in the biosynthesis of phospholipides. J Biol Chem. 1956;222:193–214.
Sundler R, Akesson B. Biosynthesis of phosphatidylethanolamines and phosphatidylcholines from ethanolamine and choline in rat liver. Biochem J. 1975;146:309–15.
Maiseyeu A, Mihai G, Kampfrath T, Simonetti OP, Sen CK, Roy S, Rajagopalan S, Parthasarathy S. Gadolinium-containing phosphatidylserine liposomes for molecular imaging of atherosclerosis. J Lipid Res. 2009;50:2157–63.
Chaudhry F, Kawai H, Johnson KW, Narula N, Shekhar A, Chaudhry F, Nakahara T, Tanimoto T, Kim D, Adapoe M, et al. Molecular imaging of apoptosis in atherosclerosis by targeting cell membrane phospholipid asymmetry. J Am Coll Cardiol. 2020;76:1862–74.
Schutters K, Kusters DH, Chatrou ML, Montero-Melendez T, Donners M, Deckers NM, Krysko DV, Vandenabeele P, Perretti M, Schurgers LJ, Reutelingsperger CP. Cell surface-expressed phosphatidylserine as therapeutic target to enhance phagocytosis of apoptotic cells. Cell Death Differ. 2013;20:49–56.
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.
Ridgway ND. The role of phosphatidylcholine and choline metabolites to cell proliferation and survival. Crit Rev Biochem Mol Biol. 2013;48:20–38.
Liu J, Yuan J, Zhao J, Zhang L, Wang Q, Wang G. Serum metabolomic patterns in young patients with ischemic stroke: a case study. Metabolomics. 2021;17:24.
Djekic D, Pinto R, Repsilber D, Hyotylainen T, Henein M. Serum untargeted lipidomic profiling reveals dysfunction of phospholipid metabolism in subclinical coronary artery disease. Vasc Health Risk Manag. 2019;15:123–35.
McMaster CR. From yeast to humans—roles of the Kennedy pathway for phosphatidylcholine synthesis. FEBS Lett. 2018;592:1256–72.
Mannheim D, Herrmann J, Versari D, Gossl M, Meyer FB, McConnell JP, Lerman LO, Lerman A. Enhanced expression of Lp-PLA2 and lysophosphatidylcholine in symptomatic carotid atherosclerotic plaques. Stroke. 2008;39:1448–55.
Vitali C, Cuchel M. Controversial role of lecithin: cholesterol acyltransferase in the development of atherosclerosis: new insights from an LCAT activator. Arterioscler Thromb Vasc Biol. 2021;41:377–9.
Aiyar N, Disa J, Ao Z, Ju H, Nerurkar S, Willette RN, Macphee CH, Johns DG, Douglas SA. Lysophosphatidylcholine induces inflammatory activation of human coronary artery smooth muscle cells. Mol Cell Biochem. 2007;295:113–20.
Schmitz G, Ruebsaamen K. Metabolism and atherogenic disease association of lysophosphatidylcholine. Atherosclerosis. 2010;208:10–8.
Polonis K, Wawrzyniak R, Daghir-Wojtkowiak E, Szyndler A, Chrostowska M, Melander O, Hoffmann M, Kordalewska M, Raczak-Gutknecht J, Bartosinska E, et al. Metabolomic signature of early vascular aging (EVA) in hypertension. Front Mol Biosci. 2020;7:12.
Goncalves I, Edsfeldt A, Ko NY, Grufman H, Berg K, Bjorkbacka H, Nitulescu M, Persson A, Nilsson M, Prehn C, et al. Evidence supporting a key role of Lp-PLA2-generated lysophosphatidylcholine in human atherosclerotic plaque inflammation. Arterioscler Thromb Vasc Biol. 2012;32:1505–12.
Eldrid C, Thalassinos K. Developments in tandem ion mobility mass spectrometry. Biochem Soc Trans. 2020;48:2457–66.
Zardini Buzatto A, Kwon BK, Li L. Development of a NanoLC-MS workflow for high-sensitivity global lipidomic analysis. Anal Chim Acta. 2020;1139:88–99.
Sun C, Li T, Song X, Huang L, Zang Q, Xu J, Bi N, Jiao G, Hao Y, Chen Y, et al. Spatially resolved metabolomics to discover tumor-associated metabolic alterations. Proc Natl Acad Sci U S A. 2019;116:52–7.
Alexandrov T. Spatial metabolomics and imaging mass spectrometry in the age of artificial intelligence. Annu Rev Biomed Data Sci. 2020;3:61–87.
This work was supported by Grants from the National Natural Science Foundation of China (Grant No. 82100475, No. 82201864) and Chengdu Medical Research Project (Grant No. 2022065).
Ethics approval and consent to participate
This study was approved by the Chengdu Women’s and Children’s Central Hospital Ethics Committee.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. The base peak chromatograms (BPCs) of positive ions and negative ions. A The BPCs of positive ions in the ND group. B The BPCs of negative ions in the ND group. C BPCs of positive ions in the HFD group. D The BPCs of negative ions in the HFD group. BPC, base peak chromatogram; ND, normal diet; HFD, high-fat diet.
. The multivariate statistical analysis. A PCA analysis. B OPLS-DA. PCA, principal component analysis; OPLS-DA, orthogonal partial least squares discriminant analysis.
. IPA network pathway analysis. The IPA network pathway analysis revealed that these differentially expressed lipids were related to NF-κB signaling, PI3K/AKT signaling, and JAK/STAT signaling. IPA, ingenuity pathway analysis; CP, classical signaling pathway.
. Chromatogram analysis of various lipid components. A A total of 131 differential lipids, including 53 TG (40.46%), 24 PE (18.32%), 10 PC (7.63%), 10 PS (7.63%), 7 SM (5.34%), 7 PA (5.34%) and 20 other lipids (15.27%), were identified. B TG, C PE, D PC, E PS, F PA, G SM, H Other lipids. TG, triglyceride; PE, phosphatidylethanolamine; PC, phosphatidylcholine; PS, phosphatidylserine; PA, phosphatidic acid; SM, sphingomyelin; LPE, lysophosphatidylethanolamine; LPC, lysophosphatidylcholine; DG, diglyceride; CL, cardiolipin; PG, phosphatidylglycerol; PEt, phosphatidylethanol; dMePE, dimethyl phosphatidylethanolamine.
About this article
Cite this article
Zhang, L., Xiong, L., Fan, L. et al. Vascular lipidomics analysis reveales increased levels of phosphocholine and lysophosphocholine in atherosclerotic mice. Nutr Metab (Lond) 20, 1 (2023). https://doi.org/10.1186/s12986-022-00723-y