Comprehensive Single-Platform Lipidomic/Metabolomic Analysis Using Supercritical Fluid Chromatography-Mass Spectrometry

Supercritical fluid chromatography (SFC) is a rapidly expanding technique in the analysis of nonpolar to moderately polar substances and, more recently, also in the analysis of compounds with higher polarity. Herein, we demonstrate a proof of concept for the application of a commercial SFC instrument with electrospray ionization-mass spectrometry (MS) detection as a platform for the comprehensive analysis of metabolites with the full range of polarities, from nonpolar lipids up to highly polar metabolites. The developed single-platform SFC-MS lipidomic/metabolomic method is based on two consecutive injections of lipid and polar metabolite extracts from biphase methyl tert-butyl ether extraction using a diol column and two different gradient programs of methanol-water-ammonium formate modifier. Detailed development of the method focused mainly on the pressure limits of the system, the long-term repeatability of results, and the chromatographic performance, including optimization of the flow rate program, modifier composition and gradient, and injection solvent selection. The developed method enabled fast and comprehensive analysis of lipids and polar metabolites from plasma within a 24 min cycle with two injections using a simple analytical platform based on a single instrument, column, and mobile phase. Finally, the results from SFC-MS analysis of polar metabolites were compared with widely established liquid chromatography MS analysis in metabolomics. The comparison showed different separation selectivity of metabolites using both methods and overall lower sensitivity of the SFC-MS due to the higher flow rate and worse chromatographic performance.

Lipidomic profiling of human serum enables detection of pancreatic cancer

Pancreatic cancer has the worst prognosis among all cancers. Cancer screening of body fluids may improve the survival time prognosis of patients, who are often diagnosed too late at an incurable stage. Several studies report the dysregulation of lipid metabolism in tumor cells, suggesting that changes in the blood lipidome may accompany tumor growth. Here we show that the comprehensive mass spectrometric determination of a wide range of serum lipids reveals statistically significant differences between pancreatic cancer patients and healthy controls, as visualized by multivariate data analysis. Three phases of biomarker discovery research (discovery, qualification, and verification) are applied for 830 samples in total, which shows the dysregulation of some very long chain sphingomyelins, ceramides, and (lyso)phosphatidylcholines. The sensitivity and specificity to diagnose pancreatic cancer are over 90%, which outperforms CA 19-9, especially at an early stage, and is comparable to established diagnostic imaging methods. Furthermore, selected lipid species indicate a potential as prognostic biomarkers.

Lipidomic and metabolomic analysis reveals changes in biochemical pathways for non-small cell lung cancer tissues

Lung cancer represents one of the leading worldwide causes of cancer death, but the pathobiochemistry of this disease is still not fully understood. Here we characterize the lipidomic and metabolomic profiles of the tumor and surrounding normal tissues for 23 patients with non-small cell lung cancer. In total, 500 molecular species were identified and quantified by a combination of the lipidomic shotgun tandem mass spectrometry (MS/MS) analysis and the targeted metabolomic approach using liquid chromatography (LC) - MS/MS. The statistical evaluation includes multivariate and univariate methods with the emphasis on paired statistical approaches. Our research revealed significant changes in several biochemical pathways related to the central carbon metabolism, acylcarnitines, dipeptides as well as the disruption in the lipid metabolism observed mainly for glycerophospholipids, sphingolipids, and cholesteryl esters.

LipidQuant 1.0: automated data processing in lipid class separation-mass spectrometry quantitative workflows

SUMMARY

We present the LipidQuant 1.0 tool for automated data processing workflows in lipidomic quantitation based on lipid class separation coupled with high-resolution mass spectrometry. Lipid class separation workflows, such as hydrophilic interaction liquid chromatography or supercritical fluid chromatography, should be preferred in lipidomic quantitation due to the coionization of lipid class internal standards with analytes from the same class. The individual steps in the LipidQuant workflow are explained, including lipid identification, quantitation, isotopic correction, and reporting results. We show the application of LipidQuant data processing to a small cohort of human serum samples.

AVAILABILITY AND IMPLEMENTATION

The LipidQuant 1.0 is freely available at Zenodo https://doi.org/10.5281/zenodo.5151201 and https://holcapek.upce.cz/lipidquant.php.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

HILIC/ESI-MS determination of gangliosides and other polar lipid classes in renal cell carcinoma and surrounding normal tissues

Negative-ion hydrophilic liquid chromatography-electrospray ionization mass spectrometry (HILIC/ESI-MS) method has been optimized for the quantitative analysis of ganglioside (GM3) and other polar lipid classes, such as sulfohexosylceramides (SulfoHexCer), sulfodihexosylceramides (SulfoHex2Cer), phosphatidylglycerols (PG), phosphatidylinositols (PI), lysophosphatidylinositols (LPI), and phosphatidylserines (PS). The method is fully validated for the quantitation of the studied lipids in kidney normal and tumor tissues of renal cell carcinoma (RCC) patients based on the lipid class separation and the coelution of lipid class internal standard with the species from the same lipid class. The raw data are semi-automatically processed using our software LipidQuant and statistically evaluated using multivariate data analysis (MDA) methods, which allows the complete differentiation of both groups with 100% specificity and sensitivity. In total, 21 GM3, 28 SulfoHexCer, 26 SulfoHex2Cer, 10 PG, 19 PI, 4 LPI, and 7 PS are determined in the aqueous phase of lipidomic extracts from kidney tumor tissue samples and surrounding normal tissue samples of 20 RCC patients. S-plots allow the identification of most upregulated (PI 40:5, PI 40:4, GM3 34:1, and GM3 42:2) and most downregulated (PI 32:0, PI 34:0, PS 36:4, and LPI 16:0) lipids, which are primarily responsible for the differentiation of tumor and normal groups. Another confirmation of most dysregulated lipids is performed by the calculation of fold changes together with T and p values to highlight their statistical significance. The comparison of HILIC/ESI-MS data and matrix-assisted laser desorption/ionization mass spectrometric imaging (MALDI-MSI) data confirms that lipid dysregulation patterns are similar for both methods. Graphical abstract ᅟ.

Hydrophilic Interaction Liquid Chromatography-Mass Spectrometry Characterization of Gangliosides in Biological Samples

The hydrophilic interaction liquid chromatography (HILIC) coupled to a negative-ion electrospray ionization tandem mass spectrometry (ESI-MS/MS) method has been developed for the identification of a wide range of gangliosides in biological samples. Gangliosides consist of a backbone of sphingoid base and a polar oligosaccharide chain containing at least one sialic acid. Gangliosides are extracted by chloroform-methanol-water mixture, where an upper aqueous layer containing gangliosides and other polar lipid subclasses is further purified by C18 solid-phase extraction. The optimization of chromatographic conditions includes the column selection, mobile-phase composition, pH value, buffer type, and concentration with the goal to achieve the best chromatographic resolution and MS sensitivity. The identification of gangliosides and other polar lipids is based on accurate m/z values of [M-H]^- ions and fragment ions as well measured by high-resolution MS. The detailed interpretation of MS/MS spectra enables the generalization of fragmentation pathways, which is then used for the differentiation of a, b, and c series of gangliosides. The structural assignment is further confirmed by agreement with the predicted retention behavior in HILIC mode on the basis of the correlation among the ganglioside retention, the number of saccharide units, and their sequence. The final HILIC/ESI-MS/MS method is applied for the analysis of porcine brain, human kidney, lungs, plasma, and erythrocytes resulting in unambiguous identification of 145 ganglioside species from 19 subclasses, which represents the highest number of reported gangliosides. Moreover, 71 sulfatides and 59 polar phospholipids (phosphatidylserines, phosphatidylinositols, lysophosphatidylinositols, and phosphatidylglycerols) are detected within a 15 min run.

Lipidomic analysis of biological samples: Comparison of liquid chromatography, supercritical fluid chromatography and direct infusion mass spectrometry methods

Lipidomic analysis of biological samples in a clinical research represents challenging task for analytical methods given by the large number of samples and their extreme complexity. In this work, we compare direct infusion (DI) and chromatography - mass spectrometry (MS) lipidomic approaches represented by three analytical methods in terms of comprehensiveness, sample throughput, and validation results for the lipidomic analysis of biological samples represented by tumor tissue, surrounding normal tissue, plasma, and erythrocytes of kidney cancer patients. Methods are compared in one laboratory using the identical analytical protocol to ensure comparable conditions. Ultrahigh-performance liquid chromatography/MS (UHPLC/MS) method in hydrophilic interaction liquid chromatography mode and DI-MS method are used for this comparison as the most widely used methods for the lipidomic analysis together with ultrahigh-performance supercritical fluid chromatography/MS (UHPSFC/MS) method showing promising results in metabolomics analyses. The nontargeted analysis of pooled samples is performed using all tested methods and 610 lipid species within 23 lipid classes are identified. DI method provides the most comprehensive results due to identification of some polar lipid classes, which are not identified by UHPLC and UHPSFC methods. On the other hand, UHPSFC method provides an excellent sensitivity for less polar lipid classes and the highest sample throughput within 10min method time. The sample consumption of DI method is 125 times higher than for other methods, while only 40μL of organic solvent is used for one sample analysis compared to 3.5mL and 4.9mL in case of UHPLC and UHPSFC methods, respectively. Methods are validated for the quantitative lipidomic analysis of plasma samples with one internal standard for each lipid class. Results show applicability of all tested methods for the lipidomic analysis of biological samples depending on the analysis requirements.

Correlation of lipidomic composition of cell lines and tissues of breast cancer patients using hydrophilic interaction liquid chromatography/electrospray ionization mass spectrometry and multivariate data analysis

RATIONALE

The goal of this work is the comparison of differences in the lipidomic compositions of human cell lines derived from normal and cancerous breast tissues, and tumor vs. normal tissues obtained after the surgery of breast cancer patients.

METHODS

Hydrophilic interaction liquid chromatography/electrospray ionization mass spectrometry (HILIC/ESI-MS) using the single internal standard approach and response factors is used for the determination of relative abundances of individual lipid species from five lipid classes in total lipid extracts of cell lines and tissues. The supplementary information on the fatty acyl composition is obtained by gas chromatography/mass spectrometry (GC/MS) of fatty acid methyl esters. Multivariate data analysis (MDA) methods, such as nonsupervised principal component analysis (PCA), hierarchical clustering analysis (HCA) and supervised orthogonal partial least-squares discriminant analysis (OPLS-DA), are used for the visualization of differences between normal and tumor samples and the correlation of similarity between cell lines and tissues either for tumor or normal samples.

RESULTS

MDA methods are used for differentiation of sample groups and also for identification of the most up- and downregulated lipids in tumor samples in comparison to normal samples. Observed changes are subsequently generalized and correlated with data from tumor and normal tissues of breast cancer patients. In total, 123 lipid species are identified based on their retention behavior in HILIC and observed ions in ESI mass spectra, and relative abundances are determined.

CONCLUSIONS

MDA methods are applied for a clear differentiation between tumor and normal samples both for cell lines and tissues. The most upregulated lipids are phospholipids (PL) with a low degree of unsaturation (e.g., 32:1 and 34:1) and also some highly polyunsaturated PL (e.g., 40:6), while the most downregulated lipids are PL containing polyunsaturated fatty acyls (e.g., 20:4), plasmalogens and ether lipids. Copyright © 2016 John Wiley & Sons, Ltd.

Retention behavior of lipids in reversed-phase ultrahigh-performance liquid chromatography-electrospray ionization mass spectrometry

Reversed-phase ultrahigh-performance liquid chromatography (RP-UHPLC) method using two 15cm sub-2μm particles octadecylsilica gel columns is developed with the goal to separate and unambiguously identify a large number of lipid species in biological samples. The identification is performed by the coupling with high-resolution tandem mass spectrometry (MS/MS) using quadrupole - time-of-flight (QTOF) instrument. Electrospray ionization (ESI) full scan and tandem mass spectra are measured in both polarity modes with the mass accuracy better than 5ppm, which provides a high confidence of lipid identification. Over 400 lipid species covering 14 polar and nonpolar lipid classes from 5 lipid categories are identified in total lipid extracts of human plasma, human urine and porcine brain. The general dependences of relative retention times on relative carbon number or relative double bond number are constructed and fit with the second degree polynomial regression. The regular retention patterns in homologous lipid series provide additional identification point for UHPLC/MS lipidomic analysis, which increases the confidence of lipid identification. The reprocessing of previously published data by our and other groups measured in the RP mode and ultrahigh-performance supercritical fluid chromatography on the silica column shows more generic applicability of the polynomial regression for the description of retention behavior and the prediction of retention times. The novelty of this work is the characterization of general trends in the retention behavior of lipids within logical series with constant fatty acyl length or double bond number, which may be used as an additional criterion to increase the confidence of lipid identification.

Hydrophilic interaction liquid chromatography-mass spectrometry of (lyso)phosphatidic acids, (lyso)phosphatidylserines and other lipid classes

The goal of this work is a systematic optimization of hydrophilic interaction liquid chromatography (HILIC) separation of acidic lipid classes (namely phosphatidic acids-PA, lysophosphatidic acids-LPA, phosphatidylserines-PS and lysophosphatidylserines-LPS) and other lipid classes under mass spectrometry (MS) compatible conditions. The main parameters included in this optimization are the type of stationary phases used in HILIC, pH of the mobile phase, the type and concentration of mobile phase additives. Nine HILIC columns with different chemistries (unmodified silica, modified silica using diol, 2-picolylamine, diethylamine and 1-aminoanthracene and hydride silica) are compared with the emphasis on peak shapes of acidic lipid classes. The optimization of pH is correlated with the theoretical calculation of acidobasic equilibria of studied lipid classes. The final method using the hydride column, pH 4 adjusted by formic acid and the gradient of acetonitrile and 40 mmol/L of aqueous ammonium formate provides good peak shapes for all analyzed lipid classes including acidic lipids. This method is applied for the identification of lipids in real samples of porcine brain and kidney extracts.

Lipidomic differentiation between human kidney tumors and surrounding normal tissues using HILIC-HPLC/ESI-MS and multivariate data analysis

The characterization of differences among polar lipid classes in tumors and surrounding normal tissues of 20 kidney cancer patients is performed by hydrophilic interaction liquid chromatography (HILIC) coupled to electrospray ionization mass spectrometry (ESI-MS). The detailed analysis of identified lipid classes using relative abundances of characteristic ions in negative- and positive-ion modes is used for the determination of more than 120 individual lipid species containing attached fatty acyls of different chain length and double bond number. Lipid species are described using relative abundances, providing a better visualization of lipidomic differences between tumor and normal tissues. The multivariate data analysis methods using unsupervised principal component analysis (PCA) and supervised orthogonal partial least square (OPLS) are used for the characterization of statistically significant differences in identified lipid species. Ten most significant up- and down-regulated lipids in OPLS score plots are also displayed by box plots. A notable increase of relative abundances of lipids containing four and more double bonds is detected in tumor compared to normal tissues.

Lipidomic analysis of plasma, erythrocytes and lipoprotein fractions of cardiovascular disease patients using UHPLC/MS, MALDI-MS and multivariate data analysis

Differences among lipidomic profiles of healthy volunteers, obese people and three groups of cardiovascular disease (CVD) patients are investigated with the goal to differentiate individual groups based on the multivariate data analysis (MDA) of lipidomic data from plasma, erythrocytes and lipoprotein fractions of more than 50 subjects. Hydrophilic interaction liquid chromatography on ultrahigh-performance liquid chromatography (HILIC-UHPLC) column coupled with electrospray ionization mass spectrometry (ESI-MS) is used for the quantitation of four classes of polar lipids (phosphatidylethanolamines, phosphatidylcholines, sphingomyelins and lysophosphatidylcholines), normal-phase UHPLC-atmospheric pressure chemical ionization MS (NP-UHPLC/APCI-MS) is applied for the quantitation of five classes of nonpolar lipids (cholesteryl esters, triacylglycerols, sterols, 1,3-diacylglycerols and 1,2-diacylglycerols) and the potential of matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) is tested for the fast screening of all lipids without a chromatographic separation. Obtained results are processed by unsupervised (principal component analysis) and supervised (orthogonal partial least squares) MDA approaches to highlight the largest differences among individual groups and to identify lipid molecules with the highest impact on the group differentiation.

Determination of nonpolar and polar lipid classes in human plasma, erythrocytes and plasma lipoprotein fractions using ultrahigh-performance liquid chromatography-mass spectrometry

A novel normal-phase (NP) ultrahigh-performance liquid chromatography-mass spectrometry (UHPLC/MS) method is developed for a separation and quantitation of nonpolar lipid classes occurring in human plasma, erythrocytes and plasma lipoprotein fractions. The baseline class separation of cholesteryl esters (CE), cholesterol, triacylglycerols (TG), regioisomers of 1,2- and 1,3-diacylglycerols (DG) and 1-monoacylglycerols (1-MG) is achieved using an optimized hexane - 2-propanol-acetonitrile mobile phase within 18min for all nonpolar lipid classes or only 9min excluding monoacylglycerols not detected in studied samples. The determination of individual nonpolar lipid classes is performed by the response factor approach and the use of dioleoyl ethylene glycol as a single internal standard. Polar lipid classes, such as phosphatidylglycerols (PG), phosphatidylethanolamines (PE), phosphatidylcholines (PC), sphingomyelins (SM) and lysophosphatidylcholines (LPC), are separated by hydrophilic interaction liquid chromatography (HILIC) using 5mmol/L aqueous ammonium acetate-methanol-acetonitrile gradient within 13minutes. The quantitation of polar lipid classes is done by a similar approach as for nonpolar lipid classes, but a different internal standard (sphingosyl PE d17:1/12:0) is used. The complementary information on fatty acyl profiles after the transesterification of the total lipid extract is obtained by gas chromatography with flame ionization detection (GC/FID). The applicability of developed methodology for fast and comprehensive characterization of blood lipidome is illustrated on samples of human plasma, erythrocytes, high-density lipoprotein (HDL), low-density lipoprotein (LDL) and very low-density lipoprotein (VLDL) fractions.

Determination of lipidomic differences between human breast cancer and surrounding normal tissues using HILIC-HPLC/ESI-MS and multivariate data analysis

The comprehensive approach for the lipidomic characterization of human breast cancer and surrounding normal tissues is based on hydrophilic interaction liquid chromatography (HILIC)-electrospray ionization mass spectrometry (ESI-MS) quantitation of polar lipid classes of total lipid extracts followed by multivariate data analysis using unsupervised principal component analysis (PCA) and supervised orthogonal partial least square (OPLS). This analytical methodology is applied for the detailed lipidomic characterization of ten patients with the goal to find the statistically relevant differences between tumor and normal tissues. This strategy is selected for better visualization of differences, because the breast cancer tissue is compared with the surrounding healthy tissue of the same patient, therefore changes in the lipidome are caused predominantly by the tumor growth. A large increase of total concentrations for several lipid classes is observed, including phosphatidylinositols, phosphatidylethanolamines, phosphatidylcholines, and lysophosphatidylcholines. Concentrations of individual lipid species inside the abovementioned classes are also changed, and in some cases, these differences are statistically significant. PCA and OPLS analyses enable a clear differentiation of tumor and normal tissues based on changes of their lipidome. A notable decrease of relative abundances of ether and vinylether (plasmalogen) lipid species is detected for phosphatidylethanolamines, but no difference is apparent for phosphatidylcholines.

Nontargeted quantitation of lipid classes using hydrophilic interaction liquid chromatography-electrospray ionization mass spectrometry with single internal standard and response factor approach

The identification and quantitation of a wide range of lipids in complex biological samples is an essential requirement for the lipidomic studies. High-performance liquid chromatography-mass spectrometry (HPLC/MS) has the highest potential to obtain detailed information on the whole lipidome, but the reliable quantitation of multiple lipid classes is still a challenging task. In this work, we describe a new method for the nontargeted quantitation of polar lipid classes separated by hydrophilic interaction liquid chromatography (HILIC) followed by positive-ion electrospray ionization mass spectrometry (ESI-MS) using a single internal lipid standard to which all class specific response factors (RFs) are related to. The developed method enables the nontargeted quantitation of lipid classes and molecules inside these classes in contrast to the conventional targeted quantitation, which is based on predefined selected reaction monitoring (SRM) transitions for selected lipids only. In the nontargeted quantitation method described here, concentrations of lipid classes are obtained by the peak integration in HILIC chromatograms multiplied by their RFs related to the single internal standard (i.e., sphingosyl PE, d17:1/12:0) used as common reference for all polar lipid classes. The accuracy, reproducibility and robustness of the method have been checked by various means: (1) the comparison with conventional lipidomic quantitation using SRM scans on a triple quadrupole (QqQ) mass analyzer, (2) (31)P nuclear magnetic resonance (NMR) quantitation of the total lipid extract, (3) method robustness test using subsequent measurements by three different persons, (4) method transfer to different HPLC/MS systems using different chromatographic conditions, and (5) comparison with previously published results for identical samples, especially human reference plasma from the National Institute of Standards and Technology (NIST human plasma). Results on human plasma, egg yolk and porcine liver extracts are presented and discussed.

Lipidomic profiling of biological tissues using off-line two-dimensional high-performance liquid chromatography-mass spectrometry

Lipids are important components in all biological tissues having many essential roles associated with the proper function of the organism. Their analysis in the biological tissues and body fluids is a challenging task due to the extreme sample complexity of polar lipids and to their amphiphilic character. In this work, we describe a new method for the characterization of the lipid composition in various tissues, using off-line two-dimensional coupling of hydrophilic interaction liquid chromatography (HILIC) and reversed-phase (RP) high-performance liquid chromatography coupled to electrospray ionization (ESI) and atmospheric pressure chemical ionization (APCI) mass spectrometry. In the first dimension the total lipid extracts are fractioned using HILIC into individual lipid classes. In total, 19 lipid classes (+3 regioisomeric pairs) that cover a wide range of polarities are separated in one analytical run, which is the highest number of analyzed lipid classes reported so far. The lysophospholipid regioisomers are also separated in HILIC mode followed by the identification based on the characteristic ESI mass spectra. The collected fractions of the various lipid classes are further separated in the RP mode, which offers an excellent resolution of the individual lipid species. Their ESI or APCI mass spectra give correct information on the fatty acid composition and on the individual regioisomeric positions on the glycerol skeleton. Off-line coupling of both modes enables the comprehensive analysis of plant and animal samples as illustrated on the analysis of egg yolk, soya and porcine brain tissues.