High throughput prediction of chylomicron triglycerides in human plasma by nuclear magnetic resonance and chemometrics
© Savorani et al; licensee BioMed Central Ltd. 2010
Received: 29 September 2009
Accepted: 14 May 2010
Published: 14 May 2010
The lipid content of the chylomicrons is a key biomarker and risk factor of cardiovascular diseases and for the understanding of obesity. A high throughput determination of chylomicrons in human blood plasma is outlined.
The new method, which uses a combination of Nuclear Magnetic Resonance (NMR) analysis and multivariate calibration analysis (chemometrics), is based on a correlation analysis towards the established standard method (ultracentrifugation and colorimetric test kit) and enables extraordinarily fast, inexpensive, and robust prediction of triglyceride (TG) content in chylomicrons. It is the position and shape of the complex lipid methylene resonance band that determines the chylomicron TG status and this information is extracted by the multivariate regression method.
The resulting method is a relatively simple multivariate model that facilitates parsimonious and accurate prediction of chylomicron lipids from NMR spectra of blood. The chemometric model predicts the chylomicron TG content with a correlation coefficient (R) of 0.96 when plotted against density gradient ultracentrifugation data.
The new rapid method facilitates large scale clinical and nutritional trials with inclusion of diagnostics of chylomicron status and thus creates new opportunities for research in lifestyle diseases and obesity.
Most of our current knowledge on the link between lipids, lipoprotein metabolism and cardiovascular disease (CVD) development relies on observations made on fasting individuals. Actually, we spend the vast majority of our time in a non-fasting postprandial state and recent studies indicate that plasma postprandial triglyceride (TG) concentrations are better at predicting risk of CVD than fasting levels . When fat is absorbed in the human body, it enters the circulation in the form of intestinally derived TG-rich lipoproteins, essentially chylomicrons, usually within about 15 minutes after finishing a meal. The chylomicrons are released by exocytosis in the villi of the small intestine and are then secreted into the bloodstream at the thoracic duct's connection with the left subclavian vein. Chylomicrons are large micellar lipoproteins, having a diameter of 75 to 1200 nm, and are primarily composed of TG (85%) and contain some cholesterol and cholesteryl esters. They are surrounded by a surface monolayer of phospholipids, unesterified cholesterol and specific apoproteins. Until now Low Density Lipoprotein (LDL) cholesterol has been the primary suspect in the development of coronary artery disease . However, it remains unclear why 40% of people who are highly vulnerable to suffering a stroke or heart attack have low or normal low density lipoprotein (LDL) levels. It is possible that chylomicron adhesion on arterial walls may be as important a risk factor as low density lipoprotein (LDL) cholesterol in causing strokes and heart attacks . To determine an individual's ability to clear postprandial lipids or the capacity of a food ingredient to diminish lipid uptake, TG-rich lipoproteins and, in particular, chylomicrons can be measured after an oral fat load.
The "gold standard" analytical method for measuring TG concentration in chylomicrons and lipoproteins is fractionation by density gradient ultracentrifugation followed by colorimetric determination of TG (UC-TG) . Unfortunately this is a time and labor-intensive method wherefore the number of samples that can simultaneously be analyzed is small and, as a consequence, the assessment of postprandial lipoprotein metabolism is rarely done in large clinical trials. However, as recently demonstrated [5, 6], proton NMR spectroscopy and curve-fitting represent an alternative method for the measurement of postprandial TG-rich lipoproteins in humans.
Proton NMR (1H-NMR) spectroscopy has recently become an indispensable analytical technique for characterization of complex biological samples such as metabolite mapping in tissue and body fluids . Lipoproteins have quite similar compositions and thus their subgroups give rise to very similar NMR signals. However, it has been shown that NMR signals of lipoproteins in blood are weakly shifted in frequency due to the different densities of the lipoproteins, giving rise to different local magnetic fields . Extraction of quantitative information from NMR spectra is classically performed by curve fitting and/or integration of peaks [9–13]. However, in this study it is demonstrated that a multivariate pattern recognition approach, such as that employed for lipoprotein subgroups in human blood , can make quantization more efficient and provide robust and precise determinations. The mathematical basis for the fundamental advances in high throughput metabonomics  has largely been due to the use of multivariate models such as principal component analysis (PCA)  and partial least squares (PLS) regression , and facilitated by the ongoing computerized revolution.
In this study we demonstrate how NMR spectroscopy in combination with multivariate chemometric regression can be used to efficiently and accurately determine TG content of chylomicrons in human plasma. The results are regressed towards the ultracentrifugation reference method which represents the limiting step in exploiting the full potential of the new method here presented; for example, if the ultracentrifugation method is unable to perfectly separate the chylomicron fraction from the very close very low density lipoprotein (VLDL) one. To demonstrate the performance of the new method, human blood samples were obtained after ingestion of test meals varying in fiber content, as dietary fibers are known to affect postprandial lipemia .
The presented samples are a subset of human plasma samples which was part of a larger double-blinded randomized crossover single-meal study investigating the postprandial lipoprotein response as a function of dietary fiber intake, in which 18 young apparently healthy normal to moderately overweight males participated (age: 27.2 ± 4.0; BMI: 25.4 ± 2.2 kg/m2). No exclusion criteria were applied related to lipemia, but results indicate that subjects had normal to high cholesterol and TG plasma concentrations. Baseline blood samples (0 minutes) were drawn after an overnight fast. Hereafter, the subjects consumed one of four iso-caloric test breakfast meals (differing in dietary fiber content) within 20 minutes and additional blood samples were drawn at 30, 60, 120, 180, 240, 300 and 420 minutes. The blood samples were collected in tubes containing EthyleneDiamineTetraacetic Acid (EDTA) and placed on ice. The samples were centrifuged at 2800 rpm for 15 minutes at 4°C. Plasma was extracted and stored overnight at 4°C before NMR analysis. The large study design included a total of 576 plasma samples out of which 192 samples, equally representing each sampling time point, were selected for the UC-TG measurements for calibrating the PLS regression model. The clinical part of the study was carried out at the Department of Human Nutrition (University of Copenhagen) and was approved by the Municipal Ethical Committee of Copenhagen in accordance with the Helsinki declaration (KF 01-309595). The large study is described in more detail elsewhere (Kristensen et al., under review).
Samples for NMR analysis were prepared by mixing 500 μl plasma and 60 μl D2O. 1H NMR spectra were acquired at 37°C on a Bruker Avance 500 spectrometer (Bruker Biospin GmbH, Rheinstetten, Germany) operating at a Larmor frequency of 500.13 MHz for 1H (11.75 Tesla) and using a 1H/13C Z-gradient 4 mm diameter, 120 μl active volume flow-probe. The spectra were recorded using the zgcppr Bruker standard pulse sequence. The relaxation delay was 5.0 s during which water presaturation was performed at 4.700 ppm. The 90 degree pulses had duration of 8.0 μs. The acquisition time for acquiring 32 k data point free induction decays was 1.573 s using a spectral width of 20.8278 ppm. A total of 64 scans were acquired and each free induction decay was zero-filled to 64 k points and apodized by a Lorentzian line-broadening of 0.3 Hz. Phase and baseline correction were performed manually. All spectra were referenced to the anomeric proton of -D-Glucose at 5.23 ppm. Processing was performed using Bruker Topspin™ 1.3 software. Due to occasional problems with the NMR flow probe, only 153 sample spectra out of the selected 192 were available for comparison with the UC-TG data.
Chylomicron fractions were isolated by ultracentrifugation as described in details elsewhere . In short, 3 ml of plasma for each sample was carefully overlaid with 2.5 ml salt solution of density 1.006 Kg/L in ultracentrifugation tubes which were then centrifuged for 23 minutes at 30000 rpm in an ultracentrifuge. The tubes were then sliced 45 mm from the bottom and the TG concentration was measured for the top-fraction, containing primarily chylomicrons and presumably minor amounts of very low density lipoproteins (VLDL), and for a corresponding plasma sample (total TG) by a colorimetric test kit (Roche TG, Roche Diagnostics GmbH, Mannheim, Germany). The intra- and inter-assay precisions were 0.6% and 1.4%, respectively. For reference purposes, total cholesterol was assessed using a colorimetric test kit (Roche TG, Roche Diagnostics GmbH, Mannheim, Germany), intra-assay precision is 0.9%, no inter-assay precision is provided as all samples for the same subject were analyzed in the same run. Analytical sensitivity for both TG and total cholesterol is 0.01 mmol/L.
Multivariate data analysis
Multivariate data analysis (chemometrics) is a necessity to study the full dimension of NMR metabonomics data sets. Chemometrics is based on latent variables and facilitates an inductive exploratory strategy which is of fundamental scientific importance as a complement to the ubiquitous deductive research strategies . For developing the new method, PLS regression  was used. This chemometric technique applies to the simultaneous analysis of two sets of variables on the same objects and is used in quantitative spectroscopy to correlate the spectroscopic data (X-block - fast multivariate measurements) to a univariate physico-chemical data (y-block - time consuming and laborious measurements). The main purpose of the regression analysis in this study is to build a model enabling the prediction of the reference UC-TG characteristic (y) from a measured NMR spectrum (x). In matrix notation the linear model y = Xb is obtained, where b contains the regression coefficients that are determined during the calibration step. The PLS-calculation starts with a covariance based calculation and, when using mean centered data, the elements of the first latent PLS-factor (calculated as X'y) will be proportional to the covariance of the spectral intensities of X with respect to the y vector of reference measurements. This illustrates that the information in the X matrix will be extracted guided, or supervised, by the information in y. PLS can be applied to the entire NMR spectra as it has been previously done to lipoprotein subgroups in human plasma  but, in order to develop a more parsimonious regression model, the interval partial least squares method (iPLS)  was applied. This approach has previously proven to be successful for the determination of cholesterol in rodent plasma lipoprotein fractions . The method is identical to PLS but divides the spectrum into a number of subintervals for each of which is developed a unique PLS regression model. In this study all reported PLS models are validated using full leave-one-sample-out-at-a-time cross validation in order to avoid overfitting and to present robust results. All PLS and iPLS models were performed in Matlab® (2007b, The Mathworks Inc., Natick, MA, USA) using the i toolbox available at http://www.models.life.ku.dk/source/.
NMR peak assignment
Partial Least Squares calibration
Fasted chylomicron TG (n = 20)
Postprandial chylomicron TG (n = 133)
Plasma chylomicron TG (n = 153)
Plasma total cholesterol
The PLS calibration model resulting from the selected 153 samples using the entire NMR spectral information to develop the PLS model yielded a resulting correlation coefficient towards the UC-TG concentration, R2, of 0.90 and a root-mean-square-error-of-cross-validation (RMSECV) of 0.174 mmol/L TG using just 3 principal components. While this result is exceptionally encouraging, an attempt to improve the model and its interpretation was made by application of iPLS regression. In this extension to PLS the NMR spectra are divided into a number of small regions or intervals for each of which a local PLS regression model is calculated.
Time course study
The new method for determination of chylomicron TG concentration offers a low sample volume, high throughput and less expensive method that can allow extensive measurements in large clinical metabonomics trials and in future individual diagnostics. The time that elapses from when the blood is sampled to when the result is produced is reduced from several hours (ultracentrifugation) to minutes with the new method (NMR). The results achieved are consistent with previous NMR based methods [5, 6, 8] even though the new method avoids the use of curve fitting and peak integration steps, keeping the analysis simple and robust.
The iPLS approach, which for the first time is applied to the determination of chylomicron TG, results in a simple model which focuses on the most informative region of the spectrum and will thus be less sensitive to variations occurring on the overall metabolic profile.
The presented method can, with some additional calibration efforts, be used for robust and automatic prediction of the chylomicron TG plasma content as a routine analysis in hospitals and as a research tool for studying large scale clinical investigations and detailed postprandial time trials.
When considering potential interferences for the new NMR-based method, it should be emphasized that the method is sensitive to methylene groups and their environment. This mean all long chain fatty acids present in the chylomicrons, independently if they are part of phospholipids or diglycerides etc. Nevertheless, the spectral features of cholesteryl esters and other rare fats in the chylomicrons, that potentially can contribute to the methylene resonances, are quantitatively outnumbered in the supervised analysis and the high correlation to UC-TG is the proof. The method's sensitivity to TGs in chylomicrons is due to the slow diffusion of the large chylomicron particles which affects the shape and position of the methylene resonance. This NMR information can be used to calculate the diameter of the chylomicron particles as done by Dyrby et al. , but in this supervised study this information is implicitly used.
Since UC-TG is used as response in the regression analysis, and since this method does not allow a perfect separation between chylomicrons and very low density lipoproteins (VLDL), the new method, which performs correlation analysis, will per definition also be affected by the presence of this last class of lipoproteins when assessing the chylomicron TG content. If the inclusion of minor fractions of very low density lipoproteins (VLDL) in the correlation analysis is a problem, it should result in fasted plasma samples with a significant content of chylomicron TG. However, this is clearly not the case as the fasted plasma samples are all well predicted by the new method and do not compromise the linearity of the relationship between NMR spectral features and TG content. Albeit the concentration of chylomicron TG in the fasted state should theoretically be 0 mmol/L, it is well known that some individuals do not clear lipids at the same high rate, and that chylomicrons may be present even after an overnight fast . Furthermore, a few subjects had both high fasting total TG (~2 mmol/L) (data not shown) and chylomicron TG (~0.9 mmol/L) indicating poor lipid clearance capacity, resulting in a relatively high mean fasting chylomicron TG concentration.
Finally, the time course study conducted on five subjects demonstrated good potential of the new method for monitoring variations in chylomicron TG concentration over time which is comparable or better to what can be obtained using the ultracentrifugation method. Most importantly, the low sample volume requirements and high throughput capacity of the new method will allow much more frequent time sampling for studying early events and finer details in the time courses for postprandial triglyceride (TG) concentrations after intake of standardized meals.
Due to the escalating obesity epidemic and its cardio-metabolic complications, in combination with the established relationship between postprandial elevation of lipids and CVD risk, the search for food ingredients and drugs that effectively suppress postprandial lipemia has increased in recent years. To determine an individual's ability to clear postprandial lipids or the capacity of a food ingredient to diminish lipid uptake, TG-rich lipoproteins and, in particular, chylomicrons must be measured after an oral fat load.
The new method facilitates rapid measurement of chylomicron TG and creates new opportunities for research in lifestyle diseases and obesity and may become a valuable tool in nutritional research for assessment of absorption of exogenous diet-derived lipids. Furthermore, non-fasting concentration of TG in TG-rich lipoproteins including chylomicrons may replace fasting levels in assessing CVD risk once standard reference values have been developed, and thus a high throughput NMR method may become a valuable tool in prediction of disease risk.
List of abbreviations used
Nuclear Magnetic Resonance
Partial Least Squares regression
Root Mean Square Error of Cross Validation.
The Danish Ministry of Food, Agriculture and Fisheries is acknowledged for sponsoring the project "Urinary biomarkers for eating patterns using NMR spectroscopy and Chemometrics" (3304-FVFP-060706-01) and the European Commission in the Communities 6th Framework Programme for sponsoring the project HEALTHGRAIN (FOOD-CT-2005-514008). Bruker Biospin is acknowledged for technical support and partial sponsoring of the 500 MHz NMR flow system.
- Bansal S, Buring JE, Rifai N, Mora S, Sacks FM, Ridker PM: Fasting compared with nonfasting triglycerides and risk of cardiovascular events in women. JAMA-Journal of the American Medical Association. 2007, 298: 309-316. 10.1001/jama.298.3.309.View ArticleGoogle Scholar
- Lamarche B, Tchernof A, Moorjani S, Cantin B, Dagenais GR, Lupien PJ, Despres JP: Small, dense low-density lipoprotein particles as a predictor of the risk of ischemic heart disease in men. Prospective results from the Quebec Cardiovascular Study. Circulation. 1997, 95: 69-75.View ArticleGoogle Scholar
- Kolovou GD, Anagnostopoulou KK, Daskalopoulou SS, Mikhailidis DP, Cokkinos DV: Clinical relevance of postprandial lipaemia. Curr Med Chem. 2005, 12: 1931-1945. 10.2174/0929867054546609.View ArticleGoogle Scholar
- Chapman MJ, Goldstein S, Lagrange D, Laplaud PM: A Density Gradient Ultra-Centrifugal Procedure for the Isolation of the Major Lipoprotein Classes from Human-Serum. Journal of Lipid Research. 1981, 22: 339-358.Google Scholar
- Ala-Korpela M: Potential role of body fluid 1H NMR metabonomics as a prognostic and diagnostic tool. Expert Rev Mol Diagn. 2007, 7: 761-773. 10.1586/14737220.127.116.111.View ArticleGoogle Scholar
- Tsai MY, Georgopoulos A, Otvos JD, Ordovas JM, Hanson NQ, Peacock JM, Arnett DK: Comparison of Ultracentrifugation and Nuclear Magnetic Resonance Spectroscopy in the Quantification of Triglyceride-Rich Lipoproteins after an Oral Fat Load. Clin Chem. 2004, 50: 1201-1204. 10.1373/clinchem.2004.032938.View ArticleGoogle Scholar
- Lindon JC, Nicholson JK, Holmes E, Everett JR: Metabonomics: Metabolic processes studied by NMR spectroscopy of biofluids. Concepts Magn Reson. 2000, 12: 289-320. 10.1002/1099-0534(2000)12:5<289::AID-CMR3>3.0.CO;2-W.View ArticleGoogle Scholar
- Otvos JD, Jeyarajah EJ, Bennett DW: Quantification of plasma-lipoproteins by proton nuclear-magnetic-resonance spectroscopy. Clin Chem. 1991, 37: 377-386.Google Scholar
- Burdge GC, Powell J, Dadd T, Talbot D, Civil J, Calder PC: Acute consumption of fish oil improves postprandial VLDL profiles in healthy men aged 50-65 years. British Journal of Nutrition. 2009, 102: 160-165. 10.1017/S0007114508143550.View ArticleGoogle Scholar
- Lai CQ, Arnett DK, Corella D, Straka RJ, Tsai MY, Peacock JM, Adiconis X, Parnell LD, Hixson JE, Province MA, Ordovas JM: Fenofibrate effect on triglyceride and postprandial response of apolipoprotein A5 variants - The GOLDN study. Arteriosclerosis Thrombosis and Vascular Biology. 2007, 27: 1417-1425. 10.1161/ATVBAHA.107.140103.View ArticleGoogle Scholar
- Rosenson RS, Huskin AL, Wolff DA, Helenowski IB, Rademaker AW: Fenofibrate reduces fasting and postprandial inflammatory responses among hypertriglyceridemia patients with the metabolic syndrome. Atherosclerosis. 2008, 198: 381-388. 10.1016/j.atherosclerosis.2007.12.007.View ArticleGoogle Scholar
- Shen J, Arnett DK, Perez-Martinez P, Parnell LD, Lai CQ, Peacock JM, Hixson JE, Tsai MY, Straka RJ, Hopkins PN, Ordovas JM: The effect of IL6-174C/G polymorphism on postprandial triglyceride metabolism in the GOLDN study. Journal of Lipid Research. 2008, 49: 1839-1845. 10.1194/jlr.P700033-JLR200.View ArticleGoogle Scholar
- Stein JH, Merwood MA, Bellehumeur JB, McBride PE, Wiebe DA, Sosman JM: Postprandial lipoprotein changes in patients taking antiretroviral therapy for HIV infection. Arteriosclerosis Thrombosis and Vascular Biology. 2005, 25: 399-405. 10.1161/01.ATV.0000152233.80082.9c.View ArticleGoogle Scholar
- Petersen M, Dyrby M, Toubro S, Engelsen SB, Nørgaard L, Pedersen HT, Dyerberg J: Quantification of lipoprotein subclasses by proton nuclear magnetic resonance-based partial least-squares regression models. Clin Chem. 2005, 51: 1457-1461. 10.1373/clinchem.2004.046748.View ArticleGoogle Scholar
- Nicholson JK, Connelly J, Lindon JC, Holmes E: Metabonomics: a platform for studying drug toxicity and gene function. Nature Reviews Drug Discovery. 2002, 1: 153-161. 10.1038/nrd728.View ArticleGoogle Scholar
- Hotelling H: Analysis of a complex of statistical variables into principal components. J Educ Psychol. 1933, 24: 417-441. 10.1037/h0071325.View ArticleGoogle Scholar
- Wold S, Martens H, Wold H: The multivariate calibration-problem in chemistry solved by the PLS method. Edited by: Ruhe A, Kågström B. 1983, Heidelberg: Springer, 286-293. 973Google Scholar
- Lairon D, Play B, Jourdheuil-Rahmani D: Digestible and indigestible carbohydrates: interactions with postprandial lipid metabolism. Journal of Nutritional Biochemistry. 2007, 18: 217-227. 10.1016/j.jnutbio.2006.08.001.View ArticleGoogle Scholar
- Tholstrup T, Sandstrom B, Bysted A, Holmer G: Effect of 6 dietary fatty acids on the postprandial lipid profile, plasma fatty acids, lipoprotein lipase, and cholesterol ester transfer activities in healthy young men. Am J Clin Nutr. 2001, 73: 198-208.Google Scholar
- Munck L, Nørgaard L, Engelsen SB, Bro R, Andersson CA: Chemometrics in food science - a demonstration of the feasibility of a highly exploratory, inductive evaluation strategy of fundamental scientific significance. Chemom Intell Lab Syst. 1998, 44: 31-60. 10.1016/S0169-7439(98)00074-4.View ArticleGoogle Scholar
- Nørgaard L, Saudland A, Wagner J, Nielsen JP, Munck L, Engelsen SB: Interval partial least squares regression (i PLS): A comparative chemometric study with an example from the near infrared spectroscopy. Appl Spectrosc. 2000, 54: 413-419. 10.1366/0003702001949500.View ArticleGoogle Scholar
- Kristensen M, Savorani F, Ravn-Haren G, Poulsen M, Markowski J, Larsen FH, Dragsted LO, Engelsen SB: NMR and interval PLS as reliable methods for determination of cholesterol in rodent lipoprotein fractions. Metabolomics. 2010, 6: 129-136. 10.1007/s11306-009-0181-3.View ArticleGoogle Scholar
- Dyrby M, Petersen M, Whittaker AK, Lambert L, Nørgaard L, Bro R, Engelsen SB: Analysis of lipoproteins using 2D diffusion-edited NMR spectroscopy and multi-way chemometrics. Analytica Chimica Acta. 2005, 531: 209-216. 10.1016/j.aca.2004.10.052.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.