- 1Western Human Nutrition Research Center, Agricultural Research Service, USDA, Davis, CA, United States
- 2Department of Nutrition, University of California, Davis, Davis, CA, United States
- 3West Coast Metabolomics Center, Genome Center, University of California, Davis, Davis, CA, United States
- 4Department of Surgery, Davis School of Medicine, University of California, Davis, Sacramento, CA, United States
- 5Center for Alimentary and Metabolic Science, Davis School of Medicine, University of California, Davis, Sacramento, CA, United States
The use of meal challenge tests to assess postprandial responses in carbohydrate and fat metabolism is well established in clinical nutrition research. However, challenge meal compositions and protocols remain a variable. Here, we validated a mixed macronutrient tolerance test (MMTT), containing 56-g palm oil, 59-g sucrose, and 26-g egg white protein for the parallel determination of insulin sensitivity and postprandial triglyceridemia in clinically healthy subjects. The MMTT was administered in two study populations. In one, women with overweight/obese BMIs (n = 43) involved in an 8-week dietary intervention were administered oral glucose tolerance tests (OGTTs) and MMTTs within 2 days of each other after 0, 2, and 8 weeks of the dietary intervention. In the other, 340 men and women between 18 and 64 years of age, with BMI from 18–40 kg/m2, completed the MMTT as part of a broad nutritional phenotyping effort. Postprandial blood collected at 0, 0.5, 3, and 6 h was used to measure glucose, insulin, and clinical lipid panels. The MMTT postprandial insulin-dependent glucose disposal was evaluated by using the Matsuda Index algorithm and the 0- and 3 h blood insulin and glucose measures. The resulting MMTT insulin sensitivity index (ISIMMTT) was strongly correlated (r = 0.77, p < 0.001) with the OGTT-dependent 2 h composite Matsuda index (ISIComposite), being related by the following equation: Log (ISIComposite) = [0.8751 x Log(ISIMMTT)] –0.2115. An area under the triglyceride excursion curve >11.15 mg/mL h–1 calculated from the 0, 3, and 6 h blood draws established mild-to-moderate triglyceridemia in agreement with ∼20% greater prevalence of hypertriglyceridemia than fasting indications. We also demonstrated that the product of the 0 to 3 h and 3 to 6 h triglyceride rate of change as a function of the triglyceride incremental area under the curve optimally stratified subjects by postprandial response patterns. Notably, ∼2% of the population showed minimal triglyceride appearance by 6 h, while ∼25% had increasing triglycerides through 6 h. Ultimately, using three blood draws, the MMTT allowed for the simultaneous determination of insulin sensitivity and postprandial triglyceridemia in individuals without clinically diagnosed disease.
Clinical Trial Registration: [https://clinicaltrials.gov/], identifier [NCT02298725; NCT02367287].
Introduction
Insulin sensitivity and plasma triglyceride levels are important factors in the evaluation of cardiometabolic risk, and various approaches to their assessment are available. In type 2 diabetes and its pre-clinical manifestations, a loss of insulin sensitivity leads to an exaggerated surge of blood insulin and glucose, following carbohydrate intake (1). In some individuals, insulin resistance is accompanied by hypertriglyceridemia, an independent cardiometabolic risk factor (2). While fasting triglycerides (TGs) have been classically used to establish triglyceridemic status, postprandial hypertriglyceridemia appears to provide a better predictor of cardiovascular disease risk in those without frank (i.e., pre-clinical) diabetes (2). Moreover, TGs measured 2 to 4 h postprandially, unlike fasting, have strong associations with cardiovascular events independent of both insulin resistance (IR) and high-density lipoprotein cholesterol (HDL-c) levels (3). Therefore, the simultaneous assessment of insulin sensitivity and postprandial triglyceridemic responses has value to clinical cardiovascular risk management and research, exploring the interindividual variability in this metabolic phenotype.
Approaches to assessing insulin sensitivity have understandably been glucose centric, with clinical indices based on either the homeostatic balance of insulin and glucose or how an individual regulates insulin in response to a standardized 75-g glucose challenge (4). While the glucose to insulin ratio provides potentially useful indication of insulin sensitivity in the absence of diabetes, this measure loses utility with elevated fasting glucose (5). Clinically relevant models validated against the gold standard euglycemic clamp include those assessing basal glucose and insulin homeostasis [e.g., the homeostasis model of insulin resistance (HOMA-IR) and quantitative insulin sensitivity check index (QUICKI)], and those incorporating the postprandial response to a glucose challenge [e.g., Matsuda insulin sensitivity index (ISIMatsuda), beta-cell disposition index], probing a glucose challenge response provides opportunities to assess pancreatic function and peripheral glucose disposal (6, 7). However, mixed macronutrient tolerance tests (MMTTs) are gaining popularity as they allow for a broader probe of the nutritional phenotype, including evaluations of metabolic flexibility (i.e., fuel switching), insulin sensitivity, and lipid tolerance (8–14).
Studying participants with and without type 2 diabetes, a standardized liquid MMTT was previously shown to effectively interrogate multiple metabolic parameters, including postprandial blood insulin, glucose, TGs, adipose lipolysis, amino acid metabolism, and more (15). In the current study, we characterized a similar MMTT but replacing a dairy-based protein powder with egg whites and dextrose with sucrose, reporting here its use for the assessments of surrogate measures of insulin sensitivity and lipid tolerance simultaneously in a non-diabetic population. We evaluated a liquid MMTT of 56-g palm oil, 59-g sucrose, and 26-g egg white protein to perturb both insulin and TG homeostasis. Ultimately, we show: 1) that the developed protocols provide robust measures of insulin sensitivity and postprandial triglyceridemia; 2) that the magnitude of the triglyceridemic response is variable across the population; 3) that postprandial triglycerides increased for up to 6 h in a large segment of this generally healthy study population; and 4) that ∼2% of the population shows a minimal postprandial triglyceride increases in this time frame.
Materials and Methods
Study Participants
Participants were from two independent clinical studies conducted at the United States Department of Agriculture - Agricultural Research Service - Western Human Nutrition Research Center (WHNRC) in Davis California. The individual Metabolism and Physiological Signatures Study (iMAPS; ClinicalTrials.gov: NCT02298725) recruited pre- and postmenopausal women with overweight to obese BMIs, who had <150 min/week of physical activity and ≥1 cardiometabolic risk factor (n = 44) in an 8-week feeding intervention to test the impact of diets, meeting the Dietary Guidelines for Americans on cardiometabolic risk factors. Inclusion criteria included age 20–65 years, BMI of 25–39.9 kg/m2, resting blood pressure of ≤140/90 mm Hg, and/or evidence of impaired glucose homeostasis and/or elevated fasting TGs as previously described (16). Subject body composition was determined at 0 and 8 weeks by dual-energy X-ray absorptiometry (DXA; Hologic Discovery QDR Series 84994; Hologic, Inc.). The participants were screened to ensure a sedentary lifestyle (i.e., <150 min of exercise per week) and asked to maintain their normal physical activity levels during the 8-week intervention. To confirm adherence, activity was monitored for 7-day periods 4 times over the course of the 8-week study using waist-worn accelerometers (Respironics® Actical™; Philips North America Co, Cambridge MA). The WHNRC Cross-Sectional Nutritional Phenotyping Study (Phenotyping Study; ClinicalTrials.gov: NCT02367287) recruited a cohort of 393 generally healthy individuals living near Davis, CA, starting in May of 2015. Population demographics were similar to the 2010 CA census: 61.1% (versus 60.6%), White; 15.4% (versus 13.7%), other race; 12.3% (versus 14.1%), Asian; 4.8% (versus 5.8%), Black; 4.8% (versus 4.7%), two or more races; 0.6% (versus 0.8%), native American; 1.1% (versus 0.4%), native Hawaiian or Pacific Islander. The study was designed to probe sex by age by BMI differences. Therefore, recruitment included both males and females between the ages of 18 and 65 and with BMI of 18.5–40 kg/m2 (17), and enrollment efforts strove for equal distributions across the three ages and three BMI bins by sex, and, therefore, do not constitute a true cross-section of the population. Body composition was determined by DXA and the participants underwent a variety of physiological and psychological tests to be presented in other reports. Physical activity was monitored between the two study visit days using waist-worn accelerometers (Respironics® Actical™; Philips North America Co). The final cohorts available for analysis and their sex-x-age-x-BMI-class distributions are presented in Figure 1.
Figure 1. Study recruitment and data restrictions for the mixed macronutrient tolerance test (MMTT) analysis. (A) The individual Metabolism and Physiological Signatures Study (iMAPS) recruited 52 female (F) individuals with BMIs in the overweight to the obese range. Of these, 44 completed the entire 8-week intervention. A single individual was lost to the triglyceride analysis due to difficulty in blood collection (B) The WHNRC nutritional phenotyping study recruited 393 male (M) and F participants with attempts to balance into 3 age and BMI ranges. Of these, 340 were retained for analyses.
Clinical Chemistry
Plasma glucose was measured by an enzyme-linked colorimetric assay on a Clinical Chemistry Analyzer (Cobas Integra 400 +; Roche Diagnostics Corporation). Serum insulin concentrations were measured by a competitive binding assay on an automatic analyzer (Cobas E 411; Roche Diagnostics). In the phenotyping study, TGs and cholesterol were measured on a Cobas Integra 400 +, while, in iMAPS, plasma TGs and cholesterol were run on an automatic analyzer (Beckman Coulter DXC800) at the UC Davis Health, Department of Pathology and Laboratory Medicine Clinical Laboratory. All assays met manufacturers’ recommendations with inter- and intraday variability of <2%.
Oral Glucose Tolerance Test
At weeks 0, 2, and 8 of iMAPS, 75 g of glucose was administered orally to volunteers after a ∼12 h overnight fast. Blood was collected by an antecubital vein indwelling catheter within 5 min of glucose ingestion, with four subsequent blood samples collected at 0.5, 1, 1.5, and 2 h as previously described (16).
Mixed Macronutrient Tolerance Test Pre-test Meals
All the participants had consumed a study-specific standardized pretest dinner the night before the MMTT. In iMAPS, the pretest dinner consisted of a chicken, cheese and bean burrito, corn chips, and lemonade, containing 34-g fat, 123 g of carbohydrates, and 32 g of protein (926 kcal). The fat consisted of ∼ 1:2:1 SFA/MUFA/PUFA, while the carbohydrates contained ∼25 g of simple sugars, 80 g of starch, and 20 g of fibers. In phenotyping, the pretest dinner consisted of stir-fried rice with vegetables, egg, and sweet and sour sauce, accompanied by roasted potato soup, a lemonade drink and a sorbet dessert. The meal contained ∼17-g fat, ∼160-g carbohydrates, and ∼20 g of protein (873 kcal). The fat consisted of ∼ 1:2:1 SFA/MUFA/PUFA, while the carbohydrates contained ∼90 g of simple sugars, 30 g of starch, and 9 g of fibers. Detailed compositions and links to meal recipes can be found in Supplemental Information to this manuscript.
Mixed Macronutrient Tolerance Test
The WHNRC MMTT was patterned after a similar meal developed as the “PhenFlex” challenge (15). The WHNRC MMTT contains palm oil, sucrose, and pasteurized liquid egg white as the main ingredients, with xanthan gum, gum celluloses, and maltodextrin as emulsifying agents, and vanilla, almond, and artificial butter flavorings to improve palatability. See Table 1 for the challenge meal recipe. We replaced the PhenFlex challenge dextrose with sucrose to better reflect the types of simple sugar regularly ingested, and to allow the potential impact of sucrose-derived fructose on lipid metabolism. The inclusion of a protein source also better reflects a standard meal, allows for protein-fat-carbohydrate interactions and provides an opportunity to evaluate variance in postprandial protein metabolism (15). The latter data have been collected as part of a metabolomics effort and will be the focus of a future manuscript. In iMAPS, the participants received the MMTT within 2 days, following the OGTT. In the phenotyping study, participant anthropometry, body composition, and physiological assessment data were collected; training for at-home dietary data collection and accelerometer placement and pre-test meals described below were provided on their 1st-study test day visit (17). The MMTT was conducted on the phenotyping 2nd-study visit test day.
Each participant was served the MMTT formula and 60 ml of water to cleanse the pallet after MMTT ingestion. A fasting blood sample was collected, the MMTT was consumed within 5 min, and postprandial blood collections occurred at 0.5, 3, and 6 h. Additional water intake was restricted until after the 0.5 h postprandial blood draw, followed by ad libitum deionized water for the remainder of the day. The first study participant was provided 350 g of the MMTT; however, material loss due to viscosity and sticking to the drinking vessel was noted, and so all other participants received 403 g to account for this (∼12 fluid oz. or 340 ml). The 403-g MMTT portion allowed for consumption of 380 ± 22 g (∼95%) of the provided material, controlling for the mixture viscosity. In the phenotyping study, of 357 normoglycemic participants completing the meal challenge test, 17 were excluded from the analysis as described in Figure 1. Thirteen of these consumed <85% of the 380-g average dose, with two <65%. The retained phenotyping (n = 340) and iMAPS participants consumed 99 ± 3% and 102 ± 4 of the average doses, respectively. For the participants used in the subsequent analysis, this equated to a delivered dose of 56 ± 2 g of palm oil, 59 ± 2 g of sucrose, and 26 ± 1 g of egg white protein. Based on compositional analyses performed by Covance Laboratories (Madison, WI), a dose of 380 g would provide an 840 kcal caloric load (fat– 60 cal%; carbohydrate – 28 cal%, protein– 12 cal%; moisture – 62.5 %). The fatty acid composition by weight of the total fat was 43% palmitate (16:0), 40% oleate (C18: 1n9), 9% linoleate (18:2n6), 4% stearate (18:0), and <1 % other detected fatty acid residues.
Indirect Calorimetry
The metabolic rate was assessed by trained physiologists using indirect calorimetry. Automated metabolic carts (TrueOne 2400, ParvoMedics) were used to measure resting and postprandial metabolic rates using an open circuit system. Measurement times closely coincided with blood collection times: 0, 0.75–1, 3, and 6 h postprandial. The participants had rested quietly for 5–10 min before beginning the assessment, and data were collected for ∼15–20 min with the participants in a semi-reclined position, wearing a securely fitting facemask, covering the nose and the mouth. Participant exhalation of inhaled room air was directed to the metabolic cart mixing chamber for volume and gas analyses. Respiratory exchange ratios (RERs) are calculated using observed volume of oxygen (V̇O2) consumed and volume of carbon dioxide (V̇CO2) produced using the equation V̇CO2/V̇O2. The resting and postprandial energy expenditure (EE) were estimated without urinary nitrogen correction using the Weir equation: EE = [(3.94 x V̇O2) + (1.1 x V̇CO2)] (18).
Insulin Sensitivity
Among the various reported procedures for estimating insulin sensitivity using plasma measures of glucose and insulin, the Matsuda index has wide acceptance and the strongest correlation with the euglycemic-hyperinsulinemic clamp (19). Insulin sensitivity was estimated from the fasting and postprandial insulin-glucose relationship as proposed by DeFronzo and Matsuda [DeFronzo and Matsuda, (20)]. Specifically, the 0 and 2 h glucose and insulin data from the OGTT were used to calculate the “composite” insulin sensitivity index (ISIComposite) (20, 21). An MMTT insulin sensitivity index (ISIMMTT) was calculated in a parallel fashion using the 0 and 3 h glucose concentration (mg/dL) and insulin concentrations (mU/L) with the following equation:
Importantly, a 2 h blood draw was not performed during the MMTT, and the 3 h blood draw represents a compromise, allowing the assessment of both insulin sensitivity and triglyceridemia while minimizing blood draws. For purposes of comparing the OGTT and MMTT, the participants were considered insulin resistant (IR) if the ISIComposite was <4.3, and insulin sensitive (IS) when above this cut-off (4, 22–24). It should be recognized that this ISI cut-off is not universally recognized and is somewhat arbitrary due to variability in enzyme-linked immunoassay antibody cross reactivity (25). However, its use here as a research tool allows binary segregation of subjects into distinct groups with altered postprandial glucose homeostasis for comparison. An individual was considered to have IR if indicated by the OGTT median category of the triplicate assessment, i.e., if two of the three determinations were IR, and then the participant was assigned as IR.
Postprandial Insulin Response Patterns
Mixed macronutrient tolerance test-dependent insulin response patterns were estimated based on secondary analysis of data reported for 75-g OGTT analyses with 0, 0.5, 1, 2, 3, and 4 h blood collections [30–32]. The following rules established the patterns based on the 0.5, and 3 blood collections in this study. Pattern I – normal = Peak insulin at 0.5 h, 3 h insulin <20% of 0.5 h insulin; Pattern II – delayed insulin decline = Peak insulin at 0.5 h, 3 h insulin <between 20 and 65% of 0.5 h insulin; Pattern III – delayed peak insulin = 3 h insulin >65% of 0.5 h insulin; Pattern IV – high-fasting insulin – 0 h insulin >50 μ Units/ml; Pattern V - low insulin = no insulin >15 μ Units/ml. The low-insulin cut-off was set at 50% of that suggested by Kraft et al. based on the 50% lower 0.5 h insulin in the MMTT vs. OGTT in the iMAPs cohort.
Fasting Triglyceridemia Assessments
The clinical practice guidelines of the Endocrine Society consider fasting TGs >150 mg/dL and <500 mg/dL as a clinical indication of mild-to-moderate hypertriglyceridemia, posing a risk for cardiovascular disease (26). Therefore, herein, the cutoffs were normal triglyceridemia (fasting TGs <150 mg/dL) or mild-to-moderate hypertriglyceridemia (fasting TGs 150–500 mg/dL) based on these criteria.
Postprandial Triglyceridemia Assessments
Postprandial triglyceridemia was estimated from areas under the curves of plasma TGs (AUCTG) by the trapezoidal rule using either the 0 and 3 h (3h AUCTG), or the 0, 3 and 6 h (6h AUCTG) TG measurements in mg/dL using the following equations, with results expressed in mg/ml h–1:
Postprandial Triglyceride Kinetic Response Phenotypes
To evaluate phenotypic variability in the rates of plasma TG change in the early and late postprandial periods, available pre-intervention iMAPS MMTT TG data (n = 43) were combined with the 340 MMTT-compliant phenotyping participants for a final cohort size of 383 individuals. The rate of TG change in the early period (kEP; i.e., 0 to 3 h) and the late period (kLP; i.e., 3 to 6 h) were calculated using the zero-order kinetic equations below, with results expressed in mg/dL min–1:
These rates were evaluated in relation to both the 6h AUCTG, and the 6 h incremental (i.e., the baseline corrected) area under the postprandial TG curve (incAUCTG). The subjects were further stratified into 5 intensity categories of either their 6h AUCTG or incAUCTG, using 20% cuts of the population-wide Log (AUCTG + 1) or Log (incAUCTG + 1) ranges.
Statistics
Data normality was assessed using the Anderson-Darling test, and, if necessary, transformations were optimized to obtain normal distributions prior to effect testing or variable inclusion in modeling efforts. Missing estimates of MMTT intake (n = 6 of 491) were replaced with the average intake. Sparse missing data for glucose, insulin, and TG measures were imputed using multivariate normal imputation, considering the entire dataset. Pearson’s correlations were used to evaluate all single-value correlations. Regressions of ISIComposite and ISIMMTT allowed transformation of the ISIMMTT into ISIComposite values. Triglyceridemia classification cutoffs were generated by logistic regression and receiver operator characteristic curve analyses (27). A least squares regression model with ISIMMTT as the outcome variable and BMI, age, sex as the fixed effects with BMI x age and BMI x sex interactions was used to evaluate relationships between these terms. Time-dependent changes in triglyceride levels within identified kinetic groups were evaluated using least squares regression mixed models with plasma triglyceride levels or fasting-corrected triglyceride levels as the outcome variables with time, the triglyceride kinetic pattern group, and AUCTG or incAUCTG intensity groups as fixed effects, with the participant as a random effect, followed by Tukey’s HSD post hoc testing. Stepwise linear regressions used to identify factors associated with postprandial phenotypes were used as a decrease in the Bayesian Information Criterion as the stopping function. All statistical analyses were conducted in Jmp Pro v 16 (SAS Institute, Cary NC, United States).
Results
Participant Characteristics
The physiological and general health status markers for the study participants are shown in Table 2. While phenotypically more homogeneous, the range of HOMA-IR observed in the iMAPS cohort was similar to that of the phenotyping study subjects and included both insulin-sensitive and insulin-resistant participants.
ISIMMTT calibration with ISIComposite– insulin sensitivity assessed by the OGTT-based ISIComposite at the baseline and after 2 and 8 weeks of feeding was unchanged in the iMAPS cohort (16). Therefore, these measures were considered replicate assessments of participant insulin sensitivity. The ISIComposite and ISIMMTT were calculated for each individual and had coefficients of variation (CVs) of 29 ± 19% and 25 ± 14%, respectively. The CVs for HOMA-IR calculated on the OGTT and MMTT test days were similar at 28 ± 24% and 28 ± 17%, respectively. Regression analysis allowed the determination of insulin resistance with the MMTT and a direct comparison of the two calculated ISIs. Based on the experimental data, we transformed the ISIMMTT into the ISIComposite scale using the following equation:
For the iMAPS cohort, the median assignment of the triplicate insulin sensitivity assessments was considered to represent the actual insulin sensitivity status of an individual. The ISIComposite identified 33 individuals with IR (i.e., ISI <4.3), 26 were indicated on all 3 study days. Similarly, the ISIMMTT identified 35 individuals with IR (i.e., ISI <4.3), with 29 on all 3 study days. In all, 39 of 43 (i.e., 90%) of the IR classifications were identical between the two protocols. Of the four that differed, all had borderline assignments. Three individuals were identified with IR on 2 MMTT test days, and 1 OGTT test day, while one individual was identified with IR on 1 MMTT test day, and 2 OGTT test days. By comparison, the HOMA-IR identified 31 individuals as having IR using the sex-specific cutoff of ≥2.05 for women (28). The average postprandial glucose and insulin responses for the iMAPS participants with median insulin-sensitive and insulin-resistant status determined by the OGTT and the MMTT are shown in Figures 2A–D. Five distinct patterns of postprandial insulin responses were identified analogous to those reported by Kraft et al., using OGTTs with 0, 0.5, 1, 2, 3, and 4 h blood draws (29–31). Pattern I represents a normal insulin response, Pattern II indicates a delayed return of postprandial insulin, Pattern III indicates a delayed peak insulin, Pattern IV indicates high-fasting insulin, and Pattern V indicates a low-insulin response. Based on the comparison of the OGTT and MMTT, a low-insulin response to the MMTT was set at 50% of the 30-μU/ml proposed by Kraft et al. As shown in Figure 2E, of the 43 iMAPS women with overweight to obese BMIs, 25 (i.e., ∼60%) showed a delayed insulin decline, and the remainder showed evidence of delayed peak insulin, with both patterns dominated by IR. In the phenotyping participants, all five patterns were detected (Figure 2F). The Insulin-sensitive participants primarily demonstrated Pattern I or V. Of the 123 of 340 individuals (i.e., ∼36%) estimated to have IR, 3 showed insulin Pattern I, while 67 showed insulin Pattern II, and 77 showed insulin Pattern III, while 6% showed Pattern V. A least squares regression model indicated that the Log (ISIMMTT) decreased with BMI in both males and females (p < 0.001), differing by sex when adjusted for BMI (p < 0.001), and while age itself was not a determinant (p = 0.2), age-x-BMI interactions indicated that the negative impact of BMI on ISI was elevated in the young (p = 0.0003). The prevalence of insulin resistance increased from 21 to 77% in the normal weight, overweight, and obese categories, respectively.
Figure 2. Glucose and insulin responses of insulin-sensitive and insulin-resistant iMAPS participants to a 75 g glucose oral glucose tolerance test (OGTT) and a 50 g sucrose containing mixed macronutrient meal challenge test (MMTT). The OGTT (A) glucose and (B) insulin responses of participants identified as insulin sensitive (n = 8) or insulin sensitive (n = 35) by the 2 h ISIComposite. The MMTT (C) glucose and (D) insulin responses of the participants identified as insulin sensitive (n = 10) or insulin sensitive (n = 33), the 3 h ISIMMTT. In panels (A–D): • = insulin sensitive; ○ = insulin resistance. (E) iMAPS participant MMTT-dependent postprandial insulin response patterns. (F) Phenotyping participant MMTT-dependent postprandial insulin response patterns. Postprandial response patterns were adapted from previously reported 4 h OGTT-dependent insulin response patterns (29–31). Pattern I – normal = peak insulin at 0.5 h, 3 h insulin <20% of 0.5 h insulin; Pattern II – delayed insulin decline = peak insulin at 0.5 h, 3 h insulin <between 20 and 65% of 0.5 h insulin; Pattern III – delayed peak insulin = 3 h insulin >65% of 0.5 h insulin; Pattern IV – high-fasting insulin – 0 h insulin >50-μ units/ml; Pattern V- low insulin = no insulin >15-μ units/ml. iMAPS results are the mean ± standard error of the means of triplicate measurements at 0, 2, and 8 weeks of intervention. Phenotyping results are means ± standard errors of the participants within each postprandial insulin pattern.
Triglyceridemia Assessments
Postprandial triglyceridemia following the MMTT was analyzed to assess repeatability of the measurement and to establish cutoffs for hypertriglyceridemia using this tool. The repeated measures in the iMAPS cohort were used to assess the reproducibility of the MMTT TG response, with the median triglyceridemia assignment of the triplicate assessments considered the actual triglyceridemia status for any given individual. The coefficient of variation was 17 ± 10% for the measurement of each time point across the three study visits (t0h = 17 ± 9%; t0.5h = 17 ± 11; t3h 17 ± 8%, t6h = 16 ± 10%), and 14 ± 7% and 24 ± 16% for the AUCTG and incAUCTG, respectively. Therefore, the 6 h MMTT-dependent TG excursions over the 8-week study were reasonably stable. In the phenotyping study, the 6 h postprandial AUCTG calculated using the four available data points was strongly correlated with the 6h AUCTG, not using the 0.5 h data (r2 = 1.), the 3h AUCTG (r2 = 0.95), and the fasting TG levels (r2 = 0.69).
To establish MMTT-based triglyceridemia cutoffs for normal triglyceridemia and mild-to-moderate hypertriglyceridemia, the phenotyping study postprandial TG excursions were compared to fasting triglyceridemia assignments using logistical regression analysis. Specifically, both the 3h and 6h AUCTGs were used to predict the average fasting triglyceridemia assignment in the 340 phenotyping participants, consuming at least 85% of the provided MMTT dose. Postprandial scoring was selected to maximize sensitivity at the expense of 1 specificity of the predicted outcome. As shown in Table 3, fasting TGs established the prevalence of mild-to-moderate triglyceridemia at 11% in the phenotyping participants, while the 3h AUCTG and 6h AUCTG showed 25% and 31 %, respectively. To maximize the mild-to-moderate hypertriglyceridemia identification, the 6h AUCTG was adopted going forward. Applying these estimators to the iMAPS participants, triglyceridemia assignments were generally reproducible across triplicate measures, with inconsistencies being highest in the fasting assessment and increasing as the degree of triglyceridemia approached the cutoffs. Specifically, regardless of the approach used, 37 of 43 (i.e., 86 %) iMAPS study participants showed identical triglyceridemic assignments on all test days. Of the remaining six individuals, fasting TGs indicated a single individual, while the 6h AUCTG indicated three individuals with median assignment of mild-to-moderate hypertriglyceridemia. The 6h AUCTG again estimated a higher degree of mild-to-moderate triglyceridemia (n = 22; ∼51%) relative to fasting estimates (n = 11; ∼25%) in this group of women with overweight to obese BMIs.
Table 3. Postprandial and fasting mild-to-moderate hypertriglyceridemia status agreement in the phenotyping study (n = 340).
Differences in triglyceridemia assignment between the fasting and postprandial estimators occurred near the assigned cutoffs of each approach. Considering the 71 phenotyping individuals with normal fasting but mild-to-moderate postprandial triglyceridemia, an intermediate normal-to-mild triglyceridemia group with 6h AUCTG between 11.15 and 14.47 mg/dL h–1 was established. As shown in Figure 3, this normal-to-mild triglyceridemic group had significantly higher fasting TG levels and a more pronounced postprandial response than the normal triglyceridemia group. Assessing the 383 available participants from the two studies combined, the ranges of measured 6h AUCTG (mg/ml h–1) in the normal, normal-to-mild, and mild-to-moderate postprandial triglyceridemia groups were 7.35 ± 1.97 mg/ml h–1 (n = 249), 14.5 ± 3.4 (n = 83), and 19.2 ± 5.8 (n = 51), respectively. Moreover, insulin sensitivity decreased as postprandial TGs increased, with the prevalence of IR being 40, 52, and 77% in the normal, normal-to-mild, and mild-to-moderate postprandial triglyceridemia groups, respectively.
Figure 3. Postprandial triglyceride response in phenotyping cohort individuals characterized with normal triglyceridemia, normal triglyceridemia to mild hypertriglyceridemia, mild-to-moderate hypertriglyceridemia. Cutoffs established for the normal, normal-mild, and mild-moderate postprandial triglyceridemia were <11.15, >11.15 and <14.47, and >14.47 mg/ml h–1 for the 6 h area under the triglyceride curve calculated from the 0, 3, and 6 h plasma measurements (i.e., 6h AUCTG). Measurements that do not share annotations differ at p < 0.05 by Tukey’s HSD post hoc analysis. The normal group 0 to 3 h slopes differ from other groups (p < 0.0001). The normal-mild group 3 to 6 h slopes differ from other groups (p < 0.0016). Results represent means ± standard deviations.
Postprandial Triglyceride Kinetic Response Phenotypes
To evaluate phenotypic variability, the TG kEP, and kLP, the 43 complete pre-intervention iMAPS study MMTT dataset was combined with the 340 MMTT-compliant phenotyping study dataset, yielding a cohort of 383 individuals. As shown in Figure 4, despite equivalent intake, considerable variability in postprandial serum TG behavior was observed. Grouping subjects based on their TG kEP and kLP as a function of the AUCTG defined four kinetic pattern groups: Group I - TG increase until 3 h and decreased to 6 h (n = 54; 16%); Group II – TG increased until 3 h and changed little between 3 and 6 h (n = 231; 68%); Group III – TG increased continuously through 6 h (n = 49, 14%); Group IV - no change in TG between 0 and 3 h, and a minimal increase at 6 h (n = 6, 2%). To characterize the TG concentration range distribution by the TG kinetic group, the population-wide AUCTG was also subdivided into 5 intensity categories of equal ranges [i.e., 20% cuts of the observed Log (AUCTG + 1) range]. Notably, the rates of TG change before and after the 3 h time point were independent of the fasting TG concentration. While kEP and kLP are linked, failure to obtain a return to baseline levels for most participants, not to mention the distribution of TGs between various lipoprotein particles with their own kinetic behaviors, prevents a true assessment of the meal TG absorption and elimination rates. However, as shown in Figure 5, the product of the kEP and kLP (i.e., kEP x kLP), when combined with the incAUCTG, provides a useful phenotypic descriptor of the postprandial TG behavior, with the quartiles of the Johnson normalized [(kEP x kLP)/incAUCTG], establishing four TG kinetic pattern groups, denoted A, B, C, and D. The distribution of the five incAUCTG intensity groups [i.e., 20% cuts of the observed Log (incAUCTG + 1) range] across these four incAUCTG kinetic groups is shown in Figure 6. While kLP was negatively correlated with kEP (p < 0.0001), considerable kEP-independent variability was observed. Translating this phenotypically, we propose that, after TGs appeared, they could disappear quickly, slowly, or continue to rise at the 6 h postprandial measurement. While the AUCTG intensity categories were distributed across all incAUCTG kinetic patterns, Group A had higher prevalence of mild-to-moderate hypertriglyceridemia than other kinetic phenotypes (A = 44%; B = 22%; C = 23%; D = 36%; χ2 p < 0.01).
Figure 4. The postprandial triglyceride area under the curve x kinetic response groups among 340 clinically healthy free-living individuals. Data shown are the mean ± SD for the concentrations of the participants identified within five equal intensity groups of the population-wide Log [AUCTG + 1] and one of 4 postprandial kinetic behaviors (Groups I–IV). (A) Group I plasma triglycerides appeared quickly and disappeared rapidly (n = 54; 16%); (B) Group II plasma triglycerides appeared moderately and disappeared slowly (n = 231; 68%); (C) Group III plasma triglycerides appeared continuously through 6 h (n = 49, 14%); (D) Group IV plasma triglycerides showed low or delayed postprandial appearance (n = 6, 1.8%). Postprandial responses were assigned using the following rules: Group I – [(ka*ke)/AUCTG] < 0.033 and [ka/AUCTG] >0.021; Group II –0.33 ≤[(ka*ke)/AUCTG] < 0.0056 and [ka/AUCTG] >0.01; Group III – [(ka*ke)/AUCTG] ≥0.0056 and (ka/AUCTG) >0.01; Group IV – (ka/AUCTG) ≤0.01]. Time-dependent changes in triglyceride levels within identified kinetic groups were evaluated using least squares regression mixed models with plasma triglyceride levels as the outcome variables with time, the AUCTG kinetic pattern group and the AUCTG intensity group as fixed effects, with the participant as a random effect, followed by Tukey’s HSD post hoc testing. Time points annotated with different letters within each TG kinetic group are different at p < 0.05.
Figure 5. Phenotyping cohort postprandial triglyceride kinetic analysis demonstrated significant phenotypic variance. Panels show: (A) the triglyceride rate of change in the 0 to 3 h early phase (kEP) as a function of the incremental area under the triglyceride curve (incAUCTG); (B) the triglyceride rate of change in the 3 to 6 h late phase (kLP) as a function of the incAUCTG; (C) A Johnson-normalized kEP x kLP product as a function of the Johnson-normalized incAUCTG. The quartiles of the Johnson Su [(kEP x kLP)/incAUCTG] defined four kinetic response groups (A-D): Group A – early-phase increase/substantial late-phase decrease (orange; n = 86); Group B – early-phase increase/minimal late-phase decrease (blue; n = 83); Group C – early-phase increase/no late-phase decrease (green; n = 85); Group D – early-phase increase/late-phase increase (purple; n = 84). Note, quartiles do not have the same number of participants due to a small percentage of individuals with identical values. The Johnson Su [(kEP x kLP)/incAUCTG] and kLP/incAUCTG differ between each kinetic group by one-way ANOVA with a Tukey post hoc analysis (p < 0.05). Symbols indicate the estimated postprandial triglyceridemia: • = Normal (6h AUCTG <11.15 mg/ml h–1); ○ = mild-moderate (6h AUCTG >11.15 mg/ml h– 1).
Figure 6. The postprandial triglyceride incremental area under the curve x kinetic response patterns of 340 clinically healthy free-living individuals. Triglyceride kinetic response types A, B, C, and D are defined by the quartiles of the Johnson Su [(kEP x kLP)/incAUCTG], which describes the product of the 0 to 3 h and 3 to 6 h triglyceride rates of change in relation to the incremental area under the triglyceride curve. The range of the population-wide Log [incAUCTG + 1] was further subdivided into five equal TG concentration intensity groups. (A) Plasma triglyceride kinetic response Group A showed a rapid early period increase and late-period decrease. (B) Plasma triglyceride kinetic response Group B showed modest early-period increase with minimal but significant late-period decrease. (C) Plasma triglyceride kinetic response Group C showed apparent modest early-period increase but insignificant TG change in the late period. (D) Plasma triglyceride kinetic response Group D showed elevating TG levels in both the early and late periods. The occurrence of mild-to-moderate triglyceridemia by the kinetic group was A (n = 38; 44%), B (n = 19; 22%), C (n = 20; 23%) and D (n = 30; 36%). Results are means ± stdev of each intensity group. Time-dependent changes in triglyceride levels within identified kinetic groups were evaluated using least squares regression mixed models with fasting-corrected plasma triglyceride levels as the outcome variables, time, the incAUCTG kinetic pattern group, and the incAUCTG intensity group as fixed effects, and the participant as a random effect, followed by Tukey’s HSD post hoc testing. Time points annotated with different letters within each TG kinetic group are different at p < 0.05.
Stability of Mixed Macronutrient Tolerance Test Triglyceride Kinetic Response
To assess the stability of the estimated TG kEP and kLP in responses to the MMTT, the kinetic group quartile cutoffs from the 383 measurements were applied to the 2-week and 8-week postprandial TG results of the iMAPS participants. As with the ISI and triglyceridemia assessments, the median categorical Johnson Su [(kEP x kLP)/incAUCTG] phenotype assignment of the 0-, 2-, and 8-week measurements was used to indicate the “true” phenotype of each participant. Substantial variability was observed in both the kEP and kLP (Figure 7) as well as their products. Considering the [(kEP x kLP)/incAUCTG], a 39 % misclassification rate was observed between the weekly and median kinetic group classification. Of the 43 participants, 6 (14 %) showed all measures in a single (kEP x kLP) group, and 16 (37%) had all measures assigned to one of the two adjacent groups (e.g., A-B-A). However, 12 participants (28 %) had at least two measures in one group, but a third of at least two groups away (e.g., A-A-D) and nine (21%) were assigned to different groups at each visit. The overall kEPCV was 30 ± 18%, and, for rates >0.5 mg/dL min–1 (n = 26), CVs were 26 ± 13%. In contrast, the overall kLP CV was 300 ± 765%, and, for rates >0.2 mg/dL min–1 (n = 18), CVs were 64 ± 38 %. Moreover, when controlling for the participant as a random effect, the 8-week dietary intervention was associated with a weak but significant decrease in kEP (p = 0.0038) and increase of kLP (p = 0.025). This relationship was reflected in strong correlations between the incAUCTG and kEP at 2 and 8 weeks, with the 0-week data (r2 ∼0.6, p < 0.001), but weaker relationships in kLP; (2-week vs. 0-week and 8-week kLP- r2 ∼0.2, p < 0.01; 0-week vs. 8-week kLP- p > 0.05). Therefore, the elimination rate appears to introduce the greatest degree of variability into this assessment of postprandial TG kinetic behavior.
Figure 7. The intra-individual variance in (A) early-phase postprandial triglyceride change (kEP) and (B) late-phase postprandial triglyceride change (kLP) from the plasma measured at 0, 2, and 8 weeks of the iMAPS dietary intervention. The 43 participants are ordered by decreasing kEP x kLP and colored by the median kEP x kLP kinetic quartile group (A, B, C, or D) of the three measurements. Considering all participants, and controlling for the participant as a random effect, EP increased (p = 0.012) and kLP tended to decrease (p = 0.06) over the 8-week intervention, but was not affected by diet type, age, or BMI.
Postprandial Triglyceride Kinetic Associations With Body Fat Distribution and Energetics
To identify factors associated with the rate of TG change in the early and late postprandial periods, associations between the kEP and kLP with physiological and metabolic factors were assessed. To simplify kinetic evaluations, the six subjects without minimal postprandial TG elevations (Figure 4) were removed from this analysis. Fasting TGs were strongly correlated with kEP (n = 375; RMSE = 0.92; r2 = 0.25; p < 0.0001). In stepwise linear regressions, considering body composition parameters, a model of kEP showing positive associations with the log of the android to the gynoid fat ratio [Log (And:Gyn); p < 0.0001] and negative associations with body mass (p = 0.0095) described ∼9% of the variance in this parameter (n = 375; RMSE = 0.96; r2 = 0.085; p < 0.0001). Including indirect calorimetry measures in the stepwise linear regression, kEP was best predicted by a model, including positive correlations with the Log [And:Gyn] and the baseline respiratory exchange ratio (RERt0; n = 375; RMSE = 0.86; r2 = 0.16; p < 0.0001). Similar results were found for AUCTG and incAUCTG (data not shown). In contrast, a linear model of kLP explained ∼5% of this factor (n = 375; RMSE = 0.39; r2 = 0.06; p < 0.0002) and included negative associations with the postprandial energy expenditure at 3 h (EE3h; p = 0.0001), the 0.75-1 h postprandial RER (RER0.75–1h; p = 0.016) and an interaction between these components (p = 0.037). We then ran a mixed model regression of EE, including the lean body mass, BMI, postprandial time, and either Johnson Su [(kEP x kLP)/incAUCTG] or the Johnson Su [(kEP x kLP)/incAUCTG]-defined kinetic group as main effects, with the participant as a random effect (Figure 8). The covariate-adjusted EE was negatively correlated with the Johnson Su [(kEP x kLP)/incAUCTG] (p = 0.00056), with Group A having higher EE than Groups C and D. Similar results were observed with the postprandial V̇O2 and V̇CO2, both of which were negatively correlated with kLP after adjusting for time, lean body mass, and BMI (p < 0.0005). While RER was negatively correlated with kLP when adjusted by time and lean body mass (p = 0.038), significant differences by the kinetic group were not observed.
Figure 8. Resting and postprandial energy metabolism shows subtle differences among MMTT triglyceride kinetic response groups. Mixed models of (A) resting and postprandial energy expenditure (EE) and (B) the respiratory exchange ratio (RER) were constructed using lean body mass, BMI, time, and the Johnson (kEP x kLP)/incAUCTG kinetic quartile group as fixed effects and the participant as random effect. Results are adjusted least square means ± standard errors. Results of contrast post-tests of Group A vs. other kinetic groups are shown. Constructing the same models with Johnson (kEP x kLP)/incAUCTG as a continuous variable indicated a negative correlation between this factor and EE (p = 0.0056) but not RER (p = 0.1).
Discussion
The metabolic dysregulation of glucose and TG homeostasis is linked to multiple adverse health outcomes, including type 2 diabetes, cardiovascular disease, stroke, and cognitive impairment (32). While these risk factors are associated with being overweight, variable risk and resilience to these metabolic perturbations exist in the general population and across BMI categories (33–35). Today, the assessments of glycated hemoglobin (HbA1c), fasting glucose, TGs, and cholesterol are commonly used to detect early signs of metabolic dysfunction and metabolic disease progression. If fasting serum glucose suggests the presence of type 2 diabetes, an oral OGTT may be prescribed to confirm diagnosis. Less appreciated is pre-diabetes, which is estimated to be present in ∼30% of the U.S. population. Lipemic risks are typically estimated from fasting plasma TGs and cholesterol levels (26, 36). While total cholesterol levels are minimally influenced by the postprandial state, remnant cholesterol [i.e., total cholesterol minus both high-density lipoprotein cholesterol (HDLc) and low-density lipoprotein cholesterol (LDLc)] is strongly correlated with postprandial triglycerides (37). Moreover, elevations in hepatic export and/or reductions in peripheral TG-rich lipoprotein clearance following a meal can result in prolonged hypertriglyceridemia despite normal fasting TG levels (38). Regardless of whether remnant cholesterol, triglycerides, or both are linked to disease risk and progression, non-fasting TG levels are clinically relevant cardiometabolic risk factors, independent of fasting TGs, LDLc, and HDLc (3, 36). Thus, while fasting blood analyte concentrations are useful for diagnosing cardiometabolic disease, probing the non-fasted state has clinical value and may aid in the identification of pre-emergent disease or differential disease risk across the population. While the clinical assessment of glucose tolerance is routine, standardized assessments of lipid tolerance are still being refined (14, 36, 39, 40). Mixed meal tests induce robust metabolic multi-organ responses that collectively reflect the adaptive responses to metabolic challenges, with limited metabolic flexibility indicating poor metabolic health (41, 42). While MMTT recommendations exist, the macronutrient sources and balance of reported MMTTs are quite variable, with broad ranges of fat (22–80 g), carbohydrate (11–75 g), and protein (3–36 g) being reported (36, 39, 42–47). In the current study, we evaluated a liquid MMTT, delivering ∼59-g palm oil, ∼59-g sucrose, and ∼29-g egg white protein to perturb both insulin and TG homeostasis. From our results in over 380 individuals, we harmonized cutoffs for postprandial insulin sensitivity and hypertriglyceridemia with reported cut-points using OGTTs and meal challenge tests, while evaluating the intra- and inter-individual variability in responses.
Previous studies have demonstrated that mixed meal challenges can accurately assess insulin sensitivity (9, 42, 44). Our results concur with these findings and extend them by allowing their transformation into an OGTT frame of reference, allowing seamless integration with the historical literature. Assessing insulin sensitivity by both a standard 2 h OGTT and the 6 h MMTT in a parallel crossover design provided nearly equivalent IR determinations and with similar precision. Specifically, for triplicate measures conducted over an 8-week intervention period, any single OGTT or MMTT measurement showed an ∼85% chance of identifying the median OGTT-assessed insulin sensitivity status for an individual. Moreover, inconsistencies in each test were highest for the participants with borderline values. It was also clear that, despite the lack of a 2 h and 4 h postprandial blood draw, sufficient resolution existed in the sampling design to identify phenotypic variation in postprandial insulin patterns that have been associated with the detection of occult diabetes (29–31).
Ingestion of a lipid-rich meal can be used to probe postprandial lipid handling, allowing segregation of individuals based on both their basal and dynamic lipid metabolism. Oral fat tolerance tests (OFTTs) are valuably clinically relevant tools for investigating postprandial lipid metabolism (14, 36, 39). Classically, OFTTs are conducted over 6 h with hourly sampling (39, 48). Reducing the OFTT to a 4 h duration can provide reliable postprandial lipemia assessments in most subjects but is less reliable in hypertriglyceridemic individuals (45, 47). Considering the ∼1.6 h lag in oral TG absorption (49), shorter time points have been deemed inappropriate. The current study supports this conjecture, as a significant proportion of the population had 6 h plasma TG levels above their 3 h time point, resulting in hypertriglyceridemia indications in more individuals. These findings suggest that a large segment of the population has plasma TG enrichment at these late postprandial time points, as previously reported with smaller cohorts (50). Optimizing postprandial TG cutoffs to identify fasting hypertriglyceridemic subjects (i.e., >150 mg/dL as defined by the Endocrine Society), we established both 3h and 6h AUCTG cutoffs, demarking normal, normal-to-mild, and mild-to-moderate hypertriglyceridemia. When applied to the 340 phenotyping participants, the 3 h and 6 h measures showed 37% and 6% misclassification rates compared to the 200-mg/dL cutoff for any post-OFTT TG measurement suggested by an expert panel of scientists and clinicians (36). Therefore, a 6h AUCTG cutoff of 11.15 mg/ml h–1, following the MMTT defined here, appears clinically relevant. However, while fasting TG measures indicated only 11% mild-to-moderate hypertriglyceridemia, the 6h AUCTG increased this to 31%. Moreover, ∼20% of the individuals showed 11.15 to 14.47 mg/ml h–1 of TGs, defining an intermediate group with apparently normal triglyceridemia to mild hypertriglyceridemia. Future studies should evaluate the cardiometabolic risk associated with postprandial TGs in this region relative to that >14.47 mg/ml h–1 of TGs, the apparent threshold for mild-to-moderate hypertriglyceridemia.
While the AUCTG integrates postprandial triglyceridemia into a single manageable value, the incAUCTG provides a better measure of the postprandial response to a high-fat meal (51). Regardless, both approaches mask the underlying kinetic behavior that can be used to phenotypically stratify subjects and provide insight into physiological mechanisms, driving postprandial lipemia (48, 52). It has become clear that post-ingestion, an early phase TG release from an enterocyte storage pool, occurs prior to the primary postprandial peak commonly occurring between 3 and 4 h (48). Using basic pharmacokinetic principles, the postprandial behavior can be segregated into pseudo-absorptive and pseudo-elimination phases, corresponding to the measurable appearance and disappearance of TGs in the blood stream. In the current study, we used the 0 to 3 h and 3 to 6 h periods to demark these early and late phases and found an array of patterns in both the magnitude and temporal kinetics of individual responses. While the average overall response and that of most subjects showed the highest measured TG concentration at the 3 h time point, roughly half of those individuals plateaued at that level and a quarter continued to rise through 6 h. Similarly, the patients with treatment-resistant cardiovascular disease and normal fasting TGs ingestion of a high-fat meal resulted in TG maximum at 4 h, where they plateaued and remained above fasting levels for up to 12 h after a high-fat meal (53). Another study compared the 6 h postprandial response to a high-fat test meal in normolipidemic lean and obese individuals, and mildly obese hyperlipidemic subjects with <20 individuals per group (52). In that study, obesity and fasting TG status were found to increase the timing of peak triglyceridemia, as well as the maximum concentration obtained, with substantial variability within groups. In the current study, while BMI and total body fat were not correlated with kEP, the And:Gyn ratio adjusted by total body mass showed a weak positive association with this rate. We also found that the kEP was positively correlated with fasting TG levels and higher in the hyperlipidemic than normolipidemic subjects. Previous studies have reported such associations between fasting and postprandial TG responses (39, 54, 55). Interestingly, kEP was also positively correlated with the RERt0, suggesting that the rates of triglyceride increases were higher when participant baseline fat oxidation rates were lower. Moreover, upon adjusting for RER, the And:Gyn ratio was a positive predictor of kEP regardless of sex, again consistent with other reported links between fat depot distribution and postprandial triglyceridemic responses (56). The TG kLP was considerably more variable, showing influences by both the magnitude of TGs achieved and postprandial whole-body energy metabolism. Specifically, results suggest that, when controlling for lean body mass, time, BMI, and the kEP, EE increased as the late phase kinetic rate decreased. The inverse relationship between the EE and kLP may reflect either subtle differences in fuel availability or utilization. Regardless, such inter- and intra-individual variability in TG clearance is not surprising, considering the complex physiological dynamics between the intestine, liver, muscle, and adipose tissue, along with genetic influences, that control postprandial TG levels in the plasma (39, 48). The infrequent sampling implemented here likely further contributes to the high variance in kLP. If plasma TG levels continued to rise after 3 h as expected in hypertriglyceridemic individuals (45, 47, 53), substantial error in the TG clearance estimate would be inherent in the collected data. Regardless, individuals with 6 h TGs that exceed their 3 h levels (i.e., those with a “D-response” type in our study) will have prolonged and underappreciated postprandial lipemia. In particular, AUCTG-based determinations of postprandial hypertriglyceridemia in such individuals would tend to undercount this condition, since those with borderline levels may exceed the 200 mg/dL at times after 6 h. Finally, if cardiometabolic risk is associated with the time that an individual maintains plasma TGs above a particular threshold, one would expect those with a “D-response type” would be at higher risk, as daily intake from multiple meals would be expected to exacerbate plasma TGs to a larger extent in these individuals. As a first step, future studies evaluating the multiple meal effects in individuals with these MMTT-defined postprandial TG phenotypes appear warranted.
Limitations
It should be appreciated that all postprandial indices of insulin sensitivity are influenced by other physiological factors, including beta-cell function and glucose absorption rates and must be interpreted carefully (22, 57). While the MMTT pretest dinner was controlled in both the iMAPS and phenotyping studies, these meals differed considerably with respect to the macronutrient balance when comparing the two studies. However, both dinners were relatively high carbohydrate meals relative to the MMTT. It has been reported that such a precursor high carbohydrate meal can blunt fat catabolism and may, therefore, exaggerate some assessments of postprandial triglyceridemia in morning meal challenges (48). Furthermore, a relatively small group of female participants (n = 43) were used to calibrate the MIMMTT cutoffs, and sex-specific cutoffs for HOMA-IR have been reported. In addition, sex-dependent differences in lipid metabolism are known, and the identified triglyceridemic cutoffs should not be considered an indication differential risk of cardiometabolic disease (58, 59). The lack of a 2 h blood draw does pose some difficulties in distinguishing between Pattern II and Pattern II postprandial responses, but the method is equivalent to an OGT for assessing insulin sensitivity. Finally, the limited number of postprandial blood draws likely increased the variability of postprandial TG kinetic assessments, particularly in the late postprandial phase. Therefore, measurements appearing to increase from 3 to 6 h may have plateaued prior to the terminal blood draw. Regardless, these individuals would appear to have prolonged postprandial triglyceridemia.
Conclusion
Using 0, 3, and 6 h blood draws following the ingestion of an MMTT comprised of a 840-kcal palm oil (60 cal%), sucrose (28 cal%), and egg white protein (12 cal%) liquid meal allowed for the simultaneous determination of insulin sensitivity and postprandial triglycerideimia status in clinically healthy individuals. The MMTT was acceptable to 95% of the participants. The 3 h ISIMMTT was transformed into a 75 g OGTT ISIComposite frame of reference and provided an equivalent indication of IR, with a cutoff of <4.3. This tool identified ∼36% of individuals in the phenotyping cohort with some impairment in carbohydrate metabolism, with prevalence increasing with BMI. Due to the recruitment strategy, this should not be interpreted as the prevalence of IR in the population. While IR was detected in individuals with normal, overweight, and obese BMIs, it became more prevalent at higher android fat distributions. As compared to fasting triglyceridemia assessments, a 6 h MMTT AUCTG of >11.15 mg/ml h–1increased the detection of hypertriglyceridemia from 11 to 31% in the phenotyping cohort. Moreover, this cut point provided equivalent stratification of normal triglyceridemia and mild-to-moderate hypertriglyceridemia indicated by any MMTT postprandial TG >200 mg/dL. Interestingly, ∼25% of the population showed rising plasma TGs through 6 h after intake, with rates of TG disappearance being weakly associated with the ability to metabolize fats. Therefore, the described procedures using an MMTT prepared from commonly available food materials provide results equivalent to an OGTT and OFTT in a single test, reporting on perturbations in both glucose homeostasis and daylong triglyceridemia.
Data Availability Statement
Requests for raw data supporting the conclusions of this article should be made by email to the corresponding author. Requests will be reviewed quarterly by a committee consisting of the study investigators.
Ethics Statement
The studies involving human participants were reviewed and approved by University of California, Davis Institutional Review Board. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JN, SA, CS, and NK conceived and developed the research plan. JN, SK, NK, and CS conducted the research. JN performed the statistical analysis. SK and KB provided statistical review. JN wrote the primary manuscript. SK, KB, NK, and SA provided significant editorial input of manuscript. All authors have primary responsibility for final content and read and approved the final manuscript.
Funding
This effort was supported by USDA Intramural Projects (2032-51530-022-00D, 2032-51530-025-00D, and 2032-51530-026-00-D) and the National Dairy Council. The USDA is an equal opportunity employer and provider. Additional support was provided by the National Center for Advancing Translational Sciences, National Institutes of Health, through UL1 TR001860 and TR000002. The authors also declare that the research received funding from Campbells Soup Company. The funder was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to thank Lindsay H. Allen and Kevin D. Laugero for their efforts in the study design and execution of the larger studies in which the current experiments were nested. We would like to acknowledge the technical support of several staff and volunteers at the WHNRC for making this research possible. In particular, Lacey Baldiviez, Eduardo Cervantes, Yasmine Bouzid, Joanne Arsenault, Dustin Burnett, Ellen Bonnell, Leslie Woodhouse, Joseph Domek, Tammy Freytag, Debra Standridge, and the BioAnalytical Laboratory support staff were instrumental in these efforts. Janet Peerson provided expert guidance and review of statistical analyses.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.877696/full#supplementary-material
Abbreviations
And, Gyn – android to gynoid fat ratio; AUC, area under the curve; BMI, body mass index; DXA, dual-energy x-ray absorptiometry; EE, energy expenditure; HDLc, high-density lipoprotein cholesterol; HOMA-IR, homeostasis model of insulin resistance; iMPAS, individual metabolism and physiological signatures study; incAUC, incremental area under the curve; IR, insulin resistance; ISI, insulin sensitivity index; kEP, rate of change in plasma triglycerides between 0 and 3 h; kLP, rate of change in plasma triglyceride between 3 and 6 h; LDLc, low-density lipoprotein cholesterol; MMTT, mixed macronutrient tolerance test; OFTT, oral fat tolerance test; OGTT, oral glucose tolerance test; RER, respiratory exchange ratio; TGs, triglycerides.
References
2. Enkhmaa B, Ozturk Z, Anuurad E, Berglund L. Postprandial lipoproteins and cardiovascular disease risk in diabetes mellitus. Curr Diab Rep. (2010) 10:61–9. doi: 10.1007/s11892-009-0088-4
3. 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. (2007) 298:309–16. doi: 10.1001/jama.298.3.309
4. Gutch M, Kumar S, Razi S, Gupta K, Gupta A. Assessment of insulin sensitivity/resistance. Indian J Endocrinol Metab. (2015) 19:160–4.
5. Quon MJ. Limitations of the fasting glucose to insulin ratio as an index of insulin sensitivity. J Clin Endocrinol Metab. (2001) 86:4615–7. doi: 10.1210/jcem.86.10.7952
6. Mohan V, Amutha A, Ranjani H, Unnikrishnan R, Datta M, Anjana RM, et al. Associations of beta-cell function and insulin resistance with youth-onset type 2 diabetes and prediabetes among Asian Indians. Diabetes Technol Ther. (2013) 15:315–22. doi: 10.1089/dia.2012.0259
7. Ishimwe MCS, Wentzel A, Shoup EM, Osei-Tutu NH, Hormenu T, Patterson AC, et al. Beta-cell failure rather than insulin resistance is the major cause of abnormal glucose tolerance in Africans: insight from the Africans in America study. BMJ Open Diabetes Res Care. (2021) 9:e002447. doi: 10.1136/bmjdrc-2021-002447
8. Maki KC, Mckenney JM, Farmer MV, Reeves MS, Dicklin MR. Indices of insulin sensitivity and secretion from a standard liquid meal test in subjects with type 2 diabetes, impaired or normal fasting glucose. Nutr J. (2009) 8:22. doi: 10.1186/1475-2891-8-22
9. Selimoglu H, Duran C, Kiyici S, Guclu M, Ersoy C, Ozkaya G, et al. Comparison of composite whole body insulin sensitivity index derived from mixed meal test and oral glucose tolerance test in insulin resistant obese subjects. Endocrine. (2009) 36:299–304. doi: 10.1007/s12020-009-9213-z
10. Phillips CM, Goumidi L, Bertrais S, Ferguson JF, Field MR, Kelly ED, et al. Additive effect of polymorphisms in the IL-6, LTA, and TNF-{alpha} genes and plasma fatty acid level modulate risk for the metabolic syndrome and its components. J Clin Endocrinol Metab. (2010) 95:1386–94. doi: 10.1210/jc.2009-1081
11. Pellis L, Van Erk MJ, Van Ommen B, Bakker GC, Hendriks HF, Cnubben NH, et al. Plasma metabolomics and proteomics profiling after a postprandial challenge reveal subtle diet effects on human metabolic status. Metabolomics. (2012) 8:347–59. doi: 10.1007/s11306-011-0320-5
12. Van Dijk SJ, Mensink M, Esser D, Feskens EJ, Muller M, Afman LA. Responses to high-fat challenges varying in fat type in subjects with different metabolic risk phenotypes: a randomized trial. PLoS One. (2012) 7:e41388. doi: 10.1371/journal.pone.0041388
13. Van Ommen B, Van Der Greef J, Ordovas JM, Daniel H. Phenotypic flexibility as key factor in the human nutrition and health relationship. Genes Nutr. (2014) 9:423. doi: 10.1007/s12263-014-0423-5
14. O’Doherty AF, Sathyapalan T, Rigby AS, Ingle L, Carroll S. The repeatability of the abbreviated (4-h) Oral Fat Tolerance Test and influence of prior acute aerobic exercise. Eur J Nutr. (2018) 57:309–18. doi: 10.1007/s00394-016-1320-z
15. Wopereis S, Wolvers D, Van Erk M, Gribnau M, Kremer B, Van Dorsten FA, et al. Assessment of inflammatory resilience in healthy subjects using dietary lipid and glucose challenges. BMC Med Genomics. (2013) 6:44. doi: 10.1186/1755-8794-6-44
16. Krishnan S, Adams SH, Allen LH, Laugero KD, Newman JW, Stephensen CB, et al. A randomized controlled-feeding trial based on the dietary guidelines for Americans on cardiometabolic health indexes. Am J Clin Nutr. (2018) 108:266–78. doi: 10.1093/ajcn/nqy113
17. Baldiviez LM, Keim NL, Laugero KD, Hwang DH, Huang L, Woodhouse LR, et al. Design and implementation of a cross-sectional nutritional phenotyping study in healthy US adults. BMC Nutr. (2017) 3:79. doi: 10.1186/s40795-017-0197-4
18. Weir JB. New methods for calculating metabolic rate with special reference to protein metabolism. J Physiol. (1949) 109:1–9. doi: 10.1113/jphysiol.1949.sp004363
19. Soonthornpun S, Setasuban W, Thamprasit A, Chayanunnukul W, Rattarasarn C, Geater A. Novel insulin sensitivity index derived from oral glucose tolerance test. J Clin Endocrinol Metab. (2003) 88:1019–23. doi: 10.1210/jc.2002-021127
20. DeFronzo RA, Matsuda M. Reduced time points to calculate the composite index. Diabetes Care. (2010) 33:e93–93. doi: 10.2337/dc10-0646
21. Matsuda M, DeFronzo RA. Insulin sensitivity indices obtained from oral glucose tolerance testing: comparison with the euglycemic insulin clamp. Diabetes Care. (1999) 22:1462–70. doi: 10.2337/diacare.22.9.1462
22. Hucking K, Watanabe RM, Stefanovski D, Bergman RN. OGTT-derived measures of insulin sensitivity are confounded by factors other than insulin sensitivity itself. Obesity (Silver Spring). (2008) 16:1938–45. doi: 10.1038/oby.2008.336
23. Takahara M, Katakami N, Kaneto H, Noguchi M, Shimomura I. Distribution of the matsuda index in Japanese healthy subjects. J Diabetes Investig. (2013) 4:369–71. doi: 10.1111/jdi.12056
24. Moshkovits Y, Rott D, Chetrit A, Dankner R. The association between insulin sensitivity indices, ECG findings and mortality: a 40-year cohort study. Cardiovasc Diabetol. (2021) 20:97. doi: 10.1186/s12933-021-01284-9
25. Chen Z, Caulfield MP, Mcphaul MJ, Reitz RE, Taylor SW, Clarke NJ. Quantitative insulin analysis using liquid chromatography-tandem mass spectrometry in a high-throughput clinical laboratory. Clin Chem. (2013) 59:1349–56. doi: 10.1373/clinchem.2012.199794
26. Berglund L, Brunzell JD, Goldberg AC, Goldberg IJ, Sacks F, Murad MH, et al. Evaluation and treatment of hypertriglyceridemia: an endocrine society clinical practice guideline. J Clin Endocrinol Metab. (2012) 97:2969–89. doi: 10.1210/jc.2011-3213
27. Lu HX, Du J, Wen DL, Sun JH, Chen MJ, Zhang AQ, et al. Development and validation of a novel predictive score for sepsis risk among trauma patients. World J Emerg Surg. (2019) 14:11. doi: 10.1186/s13017-019-0231-8
28. Gayoso-Diz P, Otero-Gonzalez A, Rodriguez-Alvarez MX, Gude F, Garcia F, De Francisco A, et al. Insulin resistance (HOMA-IR) cut-off values and the metabolic syndrome in a general adult population: effect of gender and age: EPIRCE cross-sectional study. BMC Endocr Disord. (2013) 13:47. doi: 10.1186/1472-6823-13-47
29. Kraft JR. Detection of diabetes mellitus in situ (Occult Diabetes). Lab Med. (1975) 6:10–22. doi: 10.1093/labmed/6.2.10
30. Crofts C, Schofield G, Zinn C, Wheldon M, Kraft J. Identifying hyperinsulinaemia in the absence of impaired glucose tolerance: an examination of the Kraft database. Diabetes Res Clin Pract. (2016) 118:50–7. doi: 10.1016/j.diabres.2016.06.007
31. Dinicolantonio JJ, Bhutani J, Jh OK, Crofts C. Postprandial insulin assay as the earliest biomarker for diagnosing pre-diabetes, type 2 diabetes and increased cardiovascular risk. Open Heart. (2017) 4:e000656. doi: 10.1136/openhrt-2017-000656
32. Barber TM, Kyrou I, Randeva HS, Weickert MO. Mechanisms of insulin resistance at the crossroad of obesity with associated metabolic abnormalities and cognitive dysfunction. Int J Mol Sci. (2021) 22:546. doi: 10.3390/ijms22020546
33. Meigs JB, D’agostino RB Sr., Wilson PW, Cupples LA, Nathan DM, Singer DE. Risk variable clustering in the insulin resistance syndrome. The Framingham Offspring Study. Diabetes. (1997) 46:1594–600. doi: 10.2337/diacare.46.10.1594
34. Nichols GA, Horberg M, Koebnick C, Young DR, Waitzfelder B, Sherwood NE, et al. Cardiometabolic Risk Factors Among 1.3 Million Adults With Overweight or Obesity, but Not Diabetes, in 10 Geographically Diverse Regions of the United States, 2012-2013. Prev Chronic Dis. (2017) 14:E22. doi: 10.5888/pcd14.160438
35. Kim H, Kim K, Shin S. Cardiometabolic risk factor in obese and normal weight individuals in community dwelling men. Int J Environ Res Public Health. (2020) 17:8925. doi: 10.3390/ijerph17238925
36. Kolovou GD, Mikhailidis DP, Kovar J, Lairon D, Nordestgaard BG, Ooi TC, et al. Assessment and clinical relevance of non-fasting and postprandial triglycerides: an expert panel statement. Curr Vasc Pharmacol. (2011) 9:258–70. doi: 10.2174/157016111795495549
37. Nordestgaard BG, Freiberg JJ. Clinical relevance of non-fasting and postprandial hypertriglyceridemia and remnant cholesterol. Curr Vasc Pharmacol. (2011) 9:281–6. doi: 10.2174/157016111795495585
38. Keirns BH, Sciarrillo CM, Koemel NA, Emerson SR. Fasting, non-fasting and postprandial triglycerides for screening cardiometabolic risk. J Nutr Sci. (2021) 10:e75. doi: 10.1017/jns.2021.73
39. Ryan MF, O’grada CM, Morris C, Segurado R, Walsh MC, Gibney ER, et al. Within-person variation in the postprandial lipemic response of healthy adults. Am J Clin Nutr. (2013) 97:261–7. doi: 10.3945/ajcn.112.047936
40. Tentolouris N, Kanellos PT, Siami E, Athanasopoulou E, Chaviaras N, Kolovou G, et al. Assessment of the validity and reproducibility of a novel standardized test meal for the study of postprandial triacylglycerol concentrations. Lipids. (2017) 52:675–86. doi: 10.1007/s11745-017-4275-9
41. Stroeve JH, Van Wietmarschen H, Kremer BH, Van Ommen B, Wopereis S. Phenotypic flexibility as a measure of health: the optimal nutritional stress response test. Genes Nutr. (2015) 10:459. doi: 10.1007/s12263-015-0459-1
42. Wopereis S, Stroeve JHM, Stafleu A, Bakker GCM, Burggraaf J, Van Erk MJ, et al. Multi-parameter comparison of a standardized mixed meal tolerance test in healthy and type 2 diabetic subjects: the PhenFlex challenge. Genes Nutr. (2017) 12:21. doi: 10.1186/s12263-017-0570-6
43. Zhang JQ, Ji LL, Fogt DL, Fretwell VS. Effect of exercise duration on postprandial hypertriglyceridemia in men with metabolic syndrome. J Appl Physiol (1985). (2007) 103:1339–45. doi: 10.1152/japplphysiol.00181.2007
44. Rijkelijkhuizen JM, Girman CJ, Mari A, Alssema M, Rhodes T, Nijpels G, et al. Classical and model-based estimates of beta-cell function during a mixed meal vs. an OGTT in a population-based cohort. Diabetes Res Clin Pract. (2009) 83:280–8. doi: 10.1016/j.diabres.2008.11.017
45. Maraki M, Aggelopoulou N, Christodoulou N, Katsarou C, Anapliotis P, Kavouras SA, et al. Validity of abbreviated oral fat tolerance tests for assessing postprandial lipemia. Clin Nutr. (2011) 30:852–7. doi: 10.1016/j.clnu.2011.05.003
46. Lee CJ, Brown T, Magnuson TH, Egan JM, Carlson O, Elahi D. Hormonal response to a mixed-meal challenge after reversal of gastric bypass for hypoglycemia. J Clin Endocrinol Metab. (2013) 98:E1208–12. doi: 10.1210/jc.2013-1151
47. O’Doherty AF, Sathyapalan T, Rigby AS, Ingle L, Carroll S. The repeatability of the abbreviated (4-h) Oral Fat Tolerance Test and influence of prior acute aerobic exercise. Eur J Nutr. (2016) 57:309–18.
48. Lambert JE, Parks EJ. Postprandial metabolism of meal triglyceride in humans. Biochim Biophys Acta. (2012) 1821:721–6. doi: 10.1016/j.bbalip.2012.01.006
49. Morentin Gutierrez P, Yates J, Nilsson C, Birtles S. Evolving data analysis of an oral lipid tolerance test toward the standard for the oral glucose tolerance test: cross species modeling effects of AZD7687 on plasma triacylglycerol. Pharmacol Res Perspect. (2019) 7:e00465. doi: 10.1002/prp2.465
50. Tushuizen ME, Pouwels PJ, Bontemps S, Rustemeijer C, Matikainen N, Heine RJ, et al. Postprandial lipid and apolipoprotein responses following three consecutive meals associate with liver fat content in type 2 diabetes and the metabolic syndrome. Atherosclerosis. (2010) 211:308–14. doi: 10.1016/j.atherosclerosis.2010.02.002
51. Carstensen M, Thomsen C, Hermansen K. Incremental area under response curve more accurately describes the triglyceride response to an oral fat load in both healthy and type 2 diabetic subjects. Metabolism. (2003) 52:1034–7. doi: 10.1016/s0026-0495(03)00155-0
52. Tiihonen K, Rautonen N, Alhoniemi E, Ahotupa M, Stowell J, Vasankari T. Postprandial triglyceride response in normolipidemic, hyperlipidemic and obese subjects - the influence of polydextrose, a non-digestible carbohydrate. Nutr J. (2015) 14:23. doi: 10.1186/s12937-015-0009-0
53. Samson CE, Galia AL, Llave KI, Zacarias MB, Mercado-Asis LB. Postprandial peaking and plateauing of triglycerides and VLDL in patients with underlying cardiovascular diseases despite treatment. Int J Endocrinol Metab. (2012) 10:587–93. doi: 10.5812/ijem.4783
54. Jackson KG, Knapper-Francis JM, Morgan LM, Webb DH, Zampelas A, Williams CM. Exaggerated postprandial lipaemia and lower post-heparin lipoprotein lipase activity in middle-aged men. Clin Sci (Lond). (2003) 105:457–66. doi: 10.1042/CS20030134
55. Van Oostrom AJ, Alipour A, Sijmonsma TP, Verseyden C, Dallinga-Thie GM, Plokker HW, et al. Comparison of different methods to investigate postprandial lipaemia. Neth J Med. (2009) 67:13–20.
56. Min KB, Min JY. Android and gynoid fat percentages and serum lipid levels in United States adults. Clin Endocrinol (Oxf). (2015) 82:377–87. doi: 10.1111/cen.12505
57. Pisprasert V, Ingram KH, Lopez-Davila MF, Munoz AJ, Garvey WT. Limitations in the use of indices using glucose and insulin levels to predict insulin sensitivity: impact of race and gender and superiority of the indices derived from oral glucose tolerance test in African Americans. Diabetes Care. (2013) 36:845–53. doi: 10.2337/dc12-0840
58. Olefsky J, Farquhar JW, Reaven GM. Sex difference in the kinetics of triglyceride metabolism in normal and hypertriglyceridaemic human subjects. Eur J Clin Invest. (1974) 4:121–7. doi: 10.1111/j.1365-2362.1974.tb00382.x
Keywords: insulin sensitivity, meal challenge test, phenotyping, triglyceridemia, postprandial triglyceridemia, fat tolerance test, insulin patterns
Citation: Newman JW, Krishnan S, Borkowski K, Adams SH, Stephensen CB and Keim NL (2022) Assessing Insulin Sensitivity and Postprandial Triglyceridemic Response Phenotypes With a Mixed Macronutrient Tolerance Test. Front. Nutr. 9:877696. doi: 10.3389/fnut.2022.877696
Received: 17 February 2022; Accepted: 25 March 2022;
Published: 11 May 2022.
Edited by:
Julie-Anne Nazare, Université Claude Bernard Lyon 1, FranceReviewed by:
Stephen F. Burns, Nanyang Technological University, SingaporeHidetaka Hamasaki, Hamasaki Clinic, Japan
Copyright © 2022 Newman, Krishnan, Borkowski, Adams, Stephensen and Keim. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: John W. Newman, am9obi5uZXdtYW4yQHVzZGEuZ292
†Present address: Sridevi Krishnan, Department of Pediatrics, School of Medicine, University of California, San Diego, San Diego, CA, United States