Study subjects
To determine the differences in lipid species of VLDL between MetS and non-MetS subjects, and to determine the lipid-lowering drug, that is, statin on lipid species of VLDL, this study enrolled participants at a single medical center. Those with any of the followings were excluded: significant coronary heart disease, myocardial infarction, congenital heart diseases, heart failure, significant heart valve diseases, cerebrovascular diseases, cancers, insulin therapy, and pregnant or breastfeeding women. Those who met any 3 of the following 5 components were defined as having MetS: (1) central obesity (waist ≥ 80 cm for women and ≥ 90 cm for men); (2) raised blood pressure (BP) (systolic BP ≥130mmHg or diastolic BP ≥85 mmHg or treatment of previously diagnosed HTN); (2) raised fasting glucose (≥100 mg/dL or diagnosed type 2 DM); (3) elevated triglycerides (≥150mg/dL or on triglyceride-lowering treatment); and (4) reduced HDL-C (< 50 mg/dL for women and < 40 mg/dL for men). Among the participants, 12 non-MetS and 27 MetS subjects were selected based on age and sex-matched. Fourteen patients with MetS (the MetS–off statin group) were requested to discontinue all lipid lowering drugs, that are statins, 14 days prior to the study visit for sample and data collection, while the remaining 13 patients with MetS were informed to continue ordinary medicine including statins (the MetS–on statin group). The study protocol was approved by the Kaohsiung Medical University Hospital Institutional Review Board (IRB) (KMUHIRB-E(I)-20170256) and registered with trial registration number ISRCTN 69295295 (retrospectively registered on June 9, 2020). All subjects signed an informed consent form before participating in the study. The study adhered to the principles of the Declaration of Helsinki. Each participant also underwent measurements of height, body weight, abdominal and hip circumferences, blood pressure, and heart rate at the study visit. Medical records, if available, were reviewed, and medication use was recorded.
VLDL isolation and lipid profiling for lipidome analysis
All study subjects were instructed to fast before last midnight until 20 mL venous blood draws were completed and collected in BD VACUETTE® EDTA tubes (Becton, Dickinson and Company, Franklin Lakes, NJ, USA) for subsequent VLDL isolation as described previously [14]. Briefly, the blood samples were centrifuged, separated from blood cells, and plasma was maintained for high-speed centrifugation at 10,000 rpm for 1 h to remove the upper chylomicrons. The ultracentrifuge at 40,000 rpm at 4°C for 24 h resulted in separated VLDL (density range of 1.006–1.063 g/mL) on the top. The VLDL samples were then subjected to lipid profiling using Lipotype GmbH (Dresden, Germany) [27].
Nomenclature
Lipid names and abbreviations were used: Cer, ceramide; Chol, cholesterol; DAG, diacylglycerol; HexCer, glucosyl/galactosyl ceramide; PA, phosphatidic acid; PC, phosphatidylcholine; PE, phosphatidylethanolamine; PG, phosphatidylglycerol; PI, phosphatidylinositol; PS, phosphatidylserine; and their respective lysospecies lyso-PA (LPA), lyso-PC (LPC), lyso-PE (LPE), lyso-PI (LPI), and lyso-PS (LPS); and their ether derivatives PC O-, PE O-, LPC O-, and LPE O-; SE, sterol ester; SM, sphingomyelin; SLs, sphingolipids; and TAG, triacylglycerols (also named triglycerides). Lipid species were annotated according to their molecular composition as follows: (lipid class)-(sum of carbon atoms in fatty acids):(sum of double bonds in fatty acids);(sum of hydroxyl groups in the long-chain base and the fatty acid moiety) (for example, SM-32:2;1). Where available, individual fatty acid compositions following the same rules are given in brackets (for example, 18:1;0-24:2;0).
Lipid extraction for MS lipidomics
Lipids were extracted using a two-step chloroform–methanol procedure. Samples were spiked with an internal lipid standard mixture containing cardiolipin (CL), 16:1/15:0/15:0/15:0; Cer, 18:1;2/17:0; DAG, 17:0/17:0; HexCer, 18:1;2/12:0; LPA, 17:0; LPC, 12:0; LPE, 17:1; LPG, 17:1; LPI, 17:1; LPS, 17:1; PA, 17:0/17:0; PC, 17:0/17:0; PE, 17:0/17:0; PG, 17:0/17:0; PI, 16:0/16:0; PS, 17:0/17:0; cholesterol ester (CE), 20:0; SM, 18:1;2/12:0;0; TAG, 17:0/17:0/17:0; and Chol. After the extraction, the organic phase was transferred to an infusion plate and dried in a speed vacuum concentrator. Each first-step dry extract was resuspended in 7.5 mM ammonium acetate in chloroform/methanol/propanol (1:2:4, v/v/v), and each second-step dry extract was resuspended in 33% ethanol solution of methylamine/chloroform/methanol (0.003:5:1, v/v/v). All liquid handling steps were performed using the Hamilton Robotics STARlet robotic platform with an anti-droplet control feature for organic solvent pipetting.
MS data acquisition
The samples were analyzed by direct infusion on a Q-Exactive mass spectrometer (Thermo Fisher Scientific Inc., Waltham, MA, USA) equipped with a TriVersa NanoMate ion source (Advion BioSciences, Inc. Ithaca, NY, USA). Samples were analyzed in both positive and negative ion modes with a resolution of 280,000 at m/z = 200 for MS and 17,500 for MS/MS experiments in a single acquisition. MS/MS was triggered by an inclusion list encompassing the corresponding MS mass ranges scanned in 1-Da increments. Both MS and MS/MS data were combined to monitor CE, DAG, and TAG ions as ammonium adducts; PC and PC O- as acetate adducts; and CL, PA, PE, PE O-, PG, PI, and PS as deprotonated anions. Only MS was used to monitor LPA, LPE, LPE O-, LPI, and LPS as deprotonated anions; Cer, HexCer, SM, LPC, and LPC O- as acetate adducts; and Chol as an ammonium adduct of an acetylated derivative [28].
Data analysis and post-processing
Lipid identification using LipotypeXplorer (2) was performed on the unprocessed mass spectra. For the MS-only mode, lipid identification was based on the molecular masses of intact molecules. MS/MS mode included the collision-induced fragmentation of lipid molecules, and lipid identification was based on both the intact masses and the masses of the fragments. Prior to normalization and further statistical analysis, lipid identification was filtered according to mass accuracy, occupation threshold, noise, and background [29]. The List of identified lipids and their intensities were stored in a database optimized for the particular structure inherent to lipidomic datasets. The intensity of the lipid class-specific internal standards was used for lipid quantification [30]. The identified lipid molecules were quantified by normalization to a lipid class-specific internal standard. The amounts of individual lipid molecules (species or subspecies, in p-moles) of a given lipid class were summed to yield the total amount of the lipid class. The amounts of the lipid classes were normalized to the total lipid amount, yielding mol. % per total lipids.
Lipidomic data processing
The processed lipidomic data were analyzed using LipidSig [31]. In brief, the lipid profiling data of each sample were scale-normalized to the total amount of lipids. Lipids classes with at least a two-fold change between the groups were identified as significantly changed in MetS or by the use of lipid-lowering drugs. Fisher’s exact test was then conducted to test the enrichment of the significantly changed lipids in each lipid class (such as PC, PE, and LPC). In addition, Spearmen’s correlation was used to examine the correlations between the lipid amount and clinical factors.
Laboratory testing for biochemical data
Biochemical data were obtained from the Department of Laboratory Medicine at Kaohsiung Medical University Hospital according to standard laboratory procedures. Technicians who performed the tests were blinded to the participants’ identity and data, including glucose, hemoglobin A1c (HbA1c), total cholesterol, triglycerides, VLDL, LDL-cholesterol, high-density lipoprotein (HDL)-cholesterol, alanine aminotransferase (ALT), creatinine, and uric acid.
Echocardiographic assessment
Echocardiography was performed for the measurement of left atrium (LA) diameter, maximum volume, and minimum volume by an experienced cardiologist using a transthoracic cardiac probe (Vivid 9E; General Electric Medical Systems, Horten, Norway), according to the standards of the American Society of Echocardiography [32]. LA volumes and total emptying fraction (EF) of the LA were derived using the modified Simpson’s method. Raw data were assessed while the researchers were blinded to the other data.
Measurement of electrocardiographic (ECG) parameters
Twelve-lead ECG was performed by medical technicians during the study visit. One experienced technician who was blinded to clinical information and data measured all subjects’ ECG parameters including P wave durations, PR intervals, QRS width, QTc intervals in lead II, and the duration and terminal force of P waves in lead V1 [33]. Only regular sinus rhythm was included for measurement and any arrhythmias, and bundle branch block, obvious ST-T abnormalities, and atrioventricular blocks were discarded.
Statistical analysis
All continuous variables (mean ± standard deviation) were compared among the three groups: non-MetS, MetS–off statin, and MetS–on statin by ordinary one-way ANOVA test followed by multiple comparisons using Tukey’s test. The results were considered statistically significant based on a P-value of < 0.05. Statistical analyses were performed using the statistical package in GraphPad Prism (version 9; GraphPad Software, Inc., San Diego, CA, USA) software system and SPSS statistical software (version 22; IBM Corp., Armonk, NY, USA) and SAS 9.4 software (SAS Institute Inc., Cary, NC, USA). Spearmen’s correlation was performed to determine the correlation between specific lipids and each parameter of cardiac remodelling (for the atrium and the ventricle) and the MetS.