| Literature DB >> 27797709 |
Ágnes A Fekete1,2, Carlotta Giromini3, Yianna Chatzidiakou1, D Ian Givens2, Julie A Lovegrove4.
Abstract
BACKGROUND: Cardiovascular diseases (CVDs) are the greatest cause of death globally, and their reduction is a key public-health target. High blood pressure (BP) affects 1 in 3 people in the United Kingdom, and previous studies have shown that milk consumption is associated with lower BP.Entities:
Keywords: angiotensin-converting enzyme inhibition; augmentation index; blood pressure; cardiometabolic health; caseinate; dairy; flow-mediated dilatation; milk protein; vascular function; whey
Mesh:
Substances:
Year: 2016 PMID: 27797709 PMCID: PMC5118733 DOI: 10.3945/ajcn.116.137919
Source DB: PubMed Journal: Am J Clin Nutr ISSN: 0002-9165 Impact factor: 7.045
FIGURE 1Participant inclusion flowchart of the Whey2Go study.
Baseline characteristics of participants
| Value | |
| Total (men/women), | 38 (20/18) |
| Age, y | 52.9 ± 2.1 |
| BMI, kg/m2 | 27.1 ± 0.8 |
| SBP, mm Hg | 135.5 ± 2.2 |
| DBP, mm Hg | 84.0 ± 1.9 |
| Total cholesterol, mmol/L | 5.11 ± 0.2 |
| Triacylglycerol, mmol/L | 1.12 ± 0.1 |
| Glucose, mmol/L | 5.51 ± 0.1 |
DBP, diastolic blood pressure; SBP, systolic blood pressure.
Mean ± SEM (all such values).
Twenty-four-hour ambulatory and clinic blood pressure at baseline (week 0) and postintervention (week 8) for whey protein, calcium-caseinate, and control interventions
| Whey protein | Calcium caseinate | Control | ||||||||
| Baseline | Post | Δ | Baseline | Post | Δ | Baseline | Post | Δ | ||
| 24-h SBP, mm Hg | 130.0 ± 2.2 | 127.1 ± 2.3 | −2.9 ± 1.1a | 129.7 ± 2.0 | 130.3 ± 2.6 | 0.6 ± 1.7a,b | 129.5 ± 2.3 | 130.5 ± 2.3 | 1.0 ± 1.1b | 0.041 |
| Day SBP, mm Hg | 134.2 ± 2.2 | 131.3 ± 2.4 | −2.8 ± 1.2a | 133.7 ± 2.2 | 135.4 ± 2.7 | 1.7 ± 1.7b | 133.4 ± 2.2 | 134.8 ± 2.4 | 1.3 ± 1.2a,b | 0.027 |
| Night SBP, mm Hg | 115.2 ± 2.4 | 112.2 ± 2.5 | −3.3 ± 1.5 | 115.5 ± 2.4 | 113.6 ± 3.0 | −2.0 ± 2.5 | 114.7 ± 2.3 | 115.7 ± 2.4 | 0.9 ± 1.4 | 0.110 |
| 24-h DBP, mm Hg | 78.9 ± 1.5 | 76.9 ± 1.5 | −2.0 ± 0.7a | 78.6 ± 1.3 | 79.0 ± 1.6 | 0.3 ± 1.0a,b | 78.0 ± 1.5 | 78.4 ± 1.6 | 0.5 ± 0.6b | 0.039 |
| Day DBP, mm Hg | 81.9 ± 1.6 | 79.8 ± 1.6 | −2.1 ± 0.8a | 81.6 ± 1.5 | 82.8 ± 1.7 | 1.2 ± 1.1b | 80.8 ± 1.5 | 81.5 ± 1.6 | 0.7 ± 0.7a,b | 0.027 |
| Night DBP, mm Hg | 68.0 ± 1.5 | 66.2 ± 1.5 | −1.8 ± 0.8 | 67.4 ± 1.6 | 66.1 ± 1.8 | −1.3 ± 1.5 | 67.2 ± 1.5 | 67.7 ± 1.5 | 0.5 ± 1.0 | 0.086 |
| Clinic SBP, mm Hg | 132.1 ± 2.1 | 127.9 ± 2.0 | −4.2 ± 1.3a | 129.3 ± 2.0 | 128.4 ± 2.0 | −0.9 ± 1.5a,b | 130.2 ± 2.2 | 130.9 ± 2.2 | 0.7 ± 1.4b | 0.035 |
| Clinic DBP, mm Hg | 76.7 ± 1.6 | 74.6 ± 1.5 | −2.1 ± 0.8 | 74.9 ± 1.6 | 73.6 ± 1.7 | −1.3 ± 0.8 | 74.8 ± 1.5 | 75.2 ± 1.8 | 0.5 ± 0.9 | 0.100 |
| Clinic PP, mm Hg | 55.5 ± 1.6 | 53.3 ± 1.3 | −2.2 ± 1.1 | 54.4 ± 1.9 | 54.8 ± 1.3 | 0.4 ± 1.1 | 55.4 ± 1.6 | 55.7 ± 1.5 | 0.2 ± 1.0 | 0.061 |
All values are means ± SEMs. n = 38. Baselines were significantly different (P ≤ 0.05) from one another except for 24-h SBP and DBP and day DBP. Different superscript letters within a row refer to treatment groups different from one another (P ≤ 0.05). DBP, diastolic blood pressure; PP, pulse pressure; SBP, systolic blood pressure; Δ change from baseline.
Overall between-group treatment effects for each Δ were obtained with the use of a linear mixed-model ANOVA with baseline values for the variable of interest and prognostic values such as age, sex, and BMI. Tukey-Kramer correction was used for the post hoc analysis to adjust for multiple testing.
Vascular measurements at baseline (week 0) and postintervention (week 8) for whey protein, calcium-caseinate, and control interventions
| Whey protein | Calcium caseinate | Control | ||||||||
| Baseline | Post | Δ | Baseline | Post | Δ | Baseline | Post | Δ | ||
| Endothelial function | ||||||||||
| %FMD | 4.79 ± 0.3 | 5.38 ± 0.4 | 0.59 ± 0.2a | 4.83 ± 0.3 | 4.94 ± 0.3 | 0.11 ± 0.2a | 4.79 ± 0.3 | 4.07 ± 0.3 | −0.72 ± 0.2b | <0.001 |
| Arterial stiffness | ||||||||||
| Assessed by PWA | ||||||||||
| Augmentation index at 75 beats/min, % | 17.6 ± 2.2 | 17.0 ± 2.2 | −0.6 ± 0.8 | 18.2 ± 2.3 | 17.1 ± 2.4 | −1.0 ± 0.8 | 18.2 ± 2.3 | 17.2 ± 2.1 | −1.0 ± 0.8 | 0.946 |
| Peripheral SP, mm Hg | 132.3 ± 2.2 | 127.5 ± 2.0 | −4.8 ± 1.4a | 129.1 ± 2.0 | 129.0 ± 2.0 | −0.1 ± 1.4a,b | 130.1 ± 2.2 | 131.0 ± 2.2 | 0.9 ± 1.4b | 0.007 |
| Peripheral DP, mm Hg | 76.6 ± 1.5 | 74.7 ± 1.5 | −1.9 ± 0.8 | 74.7 ± 1.6 | 73.7 ± 1.6 | −1.0 ± 0.8 | 74.7 ± 1.5 | 75.3 ± 1.8 | 0.6 ± 0.9 | 0.096 |
| Peripheral MP, mm Hg | 95.9 ± 2.0 | 92.5 ± 2.0 | −3.4 ± 0.9a | 94.1 ± 1.8 | 92.9 ± 1.9 | −1.2 ± 1.0a,b | 94.3 ± 1.8 | 94.6 ± 2.0 | 0.3 ± 0.9b | 0.032 |
| Central SP, mm Hg | 121.7 ± 2.5 | 116.7 ± 2.3 | −5.0 ± 1.3a | 118.6 ± 2.3 | 117.9 ± 2.3 | −0.7 ± 1.3a,b | 119.2 ± 2.4 | 119.6 ± 2.6 | 0.4 ± 1.3b | 0.011 |
| Central DP, mm Hg | 77.7 ± 1.6 | 75.6 ± 1.5 | −2.1 ± 0.8 | 75.8 ± 1.7 | 74.8 ± 1.7 | −1.0 ± 0.8 | 75.8 ± 1.5 | 76.4 ± 1.8 | 0.6 ± 0.9 | 0.094 |
| Central MP, mm Hg | 96.6 ± 1.8 | 93.3 ± 1.7 | −3.3 ± 0.8a | 94.1 ± 1.8 | 93.4 ± 1.8 | −0.7 ± 0.9a,b | 94.2 ± 1.8 | 94.9 ± 2.0 | 0.7 ± 1.0b | 0.023 |
| Assessed by DVP | ||||||||||
| Stiffness index, m/s | 7.53 ± 0.4 | 7.42 ± 0.4 | −0.11 ± 0.3 | 8.01 ± 0.4 | 7.92 ± 0.4 | −0.09 ± 0.2 | 7.59 ± 0.4 | 7.55 ± 0.4 | −0.04 ± 0.3 | 0.542 |
| Reflection index, % | 62.1 ± 2.9 | 66.2 ± 2.6 | 4.1 ± 2.3 | 65.0 ± 2.3 | 65.3 ± 2.7 | 0.3 ± 1.9 | 64.7 ± 2.4 | 64.7 ± 2.5 | 0.01 ± 1.6 | 0.515 |
| PPT, m/s | 246.5 ± 10.0 | 249.1 ± 9.9 | 2.6 ± 9.3 | 233.1 ± 10.5 | 235.2 ± 10.0 | 2.1 ± 5.2 | 244.2 ± 10.4 | 242.9 ± 9.5 | −1.3 ± 8.6 | 0.463 |
All values are means ± SEMs. n = 38 (except for FMD and DVP, for which n = 36). Baselines were significantly different (P ≤ 0.05) from one another except for peripheral DP and MP and central DP. Different superscript letters within a row refer to treatment groups differing from one another, P ≤ 0.05. DP, diastolic pressure; DVP, digital volume pulse; FMD; flow-mediated dilation; MP, mean pressure; PPT, peak to peak; PWA, pulse-wave analysis; SP, systolic pressure; Δ change from baseline; %FMD, flow-mediated dilatation percentage.
Overall between-group treatment effects for each Δ were obtained with the use of linear mixed-model ANOVA with baseline values for the variable of interest and prognostic values such as age, sex, BMI. Tukey-Kramer correction was used for the post hoc analysis to adjust for multiple testing.
Fasting lipid profile, indexes of insulin resistance, and vascular biomarkers at baseline (week 0) and postintervention (week 8) for whey protein, calcium-caseinate, and control interventions
| Whey protein | Calcium caseinate | Control | ||||||||
| Baseline | Post | Δ | Baseline | Post | Δ | Baseline | Post | Δ | ||
| Fasting lipid profile | ||||||||||
| Total cholesterol, mmol/L | 5.05 ± 0.2 | 4.88 ± 0.1 | −0.18 ± 0.1a | 5.01 ± 0.2 | 4.89 ± 0.2 | −0.12 ± 0.1a | 5.05 ± 0.2 | 5.12 ± 0.2 | 0.08 ± 0.1b | 0.010 |
| LDL cholesterol, mmol/L | 3.09 ± 0.1 | 3.00 ± 0.1 | −0.09 ± 0.1 | 3.08 ± 0.1 | 3.05 ± 0.1 | −0.03 ± 0.1 | 3.13 ± 0.1 | 3.15 ± 0.1 | 0.01 ± 0.1 | 0.337 |
| HDL cholesterol, mmol/L | 1.36 ± 0.1 | 1.34 ± 0.1 | −0.02 ± 0.02 | 1.35 ± 0.1 | 1.31 ± 0.1 | −0.04 ± 0.03 | 1.34 ± 0.1 | 1.37 ± 0.1 | 0.02 ± 0.01 | 0.160 |
| Total cholesterol:HDL cholesterol, mmol/L | 4.00 ± 0.2 | 3.81 ± 0.2 | −0.19 ± 0.1 | 3.97 ± 0.2 | 3.95 ± 0.2 | −0.01 ± 0.1 | 4.02 ± 0.2 | 4.02 ± 0.2 | 0.00 ± 0.1 | 0.094 |
| LDL cholesterol:HDL cholesterol, mmol/L | 2.48 ± 0.2 | 2.38 ± 0.1 | −0.10 ± 0.1 | 2.47 ± 0.2 | 2.50 ± 0.2 | 0.03 ± 0.1 | 2.53 ± 0.2 | 2.49 ± 0.2 | −0.04 ± 0.1 | 0.168 |
| Triacylglycerol, mmol/L | 1.31 ± 0.1 | 1.16 ± 0.1 | −0.15 ± 0.0a | 1.28 ± 0.1 | 1.16 ± 0.1 | −0.12 ± 0.1a,b | 1.25 ± 0.1 | 1.34 ± 0.1 | 0.08 ± 0.1b | 0.018 |
| Indexes of insulin resistance | ||||||||||
| Glucose, mmol/L | 5.04 ± 0.1 | 5.02 ± 0.1 | −0.03 ± 0.01 | 5.03 ± 0.1 | 5.10 ± 0.1 | 0.07 ± 0.1 | 5.00 ± 0.1 | 5.12 ± 0.1 | 0.12 ± 0.1 | 0.111 |
| Insulin, pmol/L | 52.4 ± 5.3 | 51.4 ± 5.0 | −1.17 ± 4.1 | 46.1 ± 4.3 | 51.7 ± 4.2 | 5.53 ± 2.7 | 47.6 ± 4.8 | 52.8 ± 5.0 | 5.25 ± 3.2 | 0.498 |
| NEFAs, μmol/L | 487.7 ± 35.0 | 414.4 ± 29.1 | −73.3 ± 30.7 | 528.7 ± 34.1 | 408.6 ± 28.4 | −120.1 ± 22.0 | 533.1 ± 32.8 | 464.0 ± 29.2 | −69.1 ± 25.8 | 0.091 |
| HOMA-IR | 1.99 ± 0.2 | 1.94 ± 0.2 | −0.05 ± 0.1 | 1.75 ± 0.2 | 1.98 ± 0.2 | 0.23 ± 0.1 | 1.79 ± 0.2 | 2.03 ± 0.2 | 0.24 ± 0.1 | 0.378 |
| QUICKI | 0.501 ± 0.0 | 0.502 ± 0.0 | 0.001 ± 0.0 | 0.503 ± 0.0 | 0.500 ± 0.0 | −0.003 ± 0.0 | 0.503 ± 0.0 | 0.500 ± 0.0 | −0.003 ± 0.0 | 0.135 |
| rQUICKI | 0.510 ± 0.0 | 0.511 ± 0.0 | 0.001 ± 0.0 | 0.511 ± 0.0 | 0.509 ± 0.0 | −0.002 ± 0.0 | 0.512 ± 0.0 | 0.509 ± 0.0 | −0.003 ± 0.0 | 0.102 |
| Vascular biomarkers | ||||||||||
| sVCAM-1, ng/mL | 637.3 ± 21.0 | 639.6 ± 22.5 | 2.3 ± 11.1a,b | 641.6 ± 22.5 | 618.5 ± 21.4 | −23.1 ± 11.7a | 627.7 ± 20.3 | 649.4 ± 23.2 | 21.7 ± 9.9b | 0.011 |
| sICAM-1, ng/mL | 224.7 ± 10.6 | 213 ± 9.0 | −11.7 ± 5.0a | 221.7 ± 10.4 | 218.4 ± 9.0 | −3.2 ± 5.6a,b | 223.6 ± 8.9 | 228.9 ± 10.7 | 5.3 ± 4.9b | 0.039 |
| E-selectin, ng/mL | 33.7 ± 2.1 | 34.15 ± 2.0 | 0.45 ± 0.8 | 34.46 ± 2.2 | 34.90 ± 2.1 | 0.44 ± 0.7 | 34.43 ± 2.2 | 36.32 ± 2.4 | 1.89 ± 0.8 | 0.303 |
| P-selectin, ng/mL | 33.5 ± 0.4 | 33.6 ± 0.3 | 0.09 ± 0.1 | 33.9 ± 0.3 | 33.5 ± 0.4 | −0.41 ± 0.2 | 33.9 ± 0.3 | 33.6 ± 0.4 | −0.28 ± 0.2 | 0.143 |
| CRP, mg/L | 1.25 ± 0.2 | 1.14 ± 0.2 | −0.11 ± 0.2 | 1.13 ± 0.1 | 1.46 ± 0.3 | 0.33 ± 0.3 | 1.10 ± 0.2 | 1.28 ± 0.2 | 0.18 ± 0.2 | 0.906 |
| IL-6, pg/mL | 1.31 ± 0.1 | 1.20 ± 0.1 | −0.11 ± 0.1 | 1.28 ± 0.1 | 1.27 ± 0.1 | −0.02 ± 0.1 | 1.21 ± 0.1 | 1.25 ± 0.1 | 0.04 ± 0.1 | 0.926 |
| TNF-α, pg/mL | 1.23 ± 0.1 | 1.24 ± 0.1 | 0.01 ± 0.1 | 1.24 ± 0.1 | 1.25 ± 0.1 | 0.01 ± 0.1 | 1.22 ± 0.1 | 1.27 ± 0.1 | 0.05 ± 0.1 | 0.856 |
| NO2, μmol/L | 1.45 ± 0.1 | 1.51 ± 0.1 | 0.06 ± 0.2 | 1.41 ± 0.1 | 1.49 ± 0.1 | 0.08 ± 0.2 | 1.51 ± 0.1 | 1.35 ± 0.1 | −0.17 ± 0.1 | 0.842 |
| NO3, μmol/L | 12.8 ± 1.9 | 12.1 ± 1.6 | −0.75 ± 1.4 | 12.0 ± 1.2 | 10.7 ± 1.1 | −1.30 ± 1.1 | 10.8 ± 1.1 | 11.7 ± 1.7 | −0.98 ± 1.2 | 0.811 |
| ACE activity, U/L | 4.65 ± 0.2 | 4.83 ± 0.2 | 0.18 ± 0.1 | 4.76 ± 0.2 | 5.01 ± 0.2 | 0.25 ± 0.1 | 4.58 ± 0.2 | 5.03 ± 0.2 | 0.45 ± 0.1 | 0.503 |
All values are means ± SEMs. n = 38 (except for sICAM-1 and E-selectin, for which n = 36; for P-selectin and IL-6, for which n = 37; and for CRP, for which n = 35). Baselines were not significantly different (P ≤ 0.05) from one another except for E-selectin and IL-6. Data were log transformed for CRP, IL-6, and TNF-α. Different superscript letters within a row refer to treatment groups different from one another (P ≤ 0.05). ACE, angiotensin-converting enzyme; CRP, C-reactive protein; NEFA, nonesterified fatty acid; QUICKI, quantitative insulin sensitivity index, rQUICKI, revised quantitative insulin sensitivity index; sICAM-1, soluble intercellular adhesion molecule 1; sVCAM-1, soluble vascular cell adhesion molecule 1; Δ, change from baseline.
Overall between-group treatment effects for each Δ were obtained with the use of a linear mixed-model ANOVA with baseline values for the variable of interest and prognostic values such as age, sex, and BMI. Tukey-Kramer correction was used for the post hoc analysis to adjust for multiple testing.
Reported nutrient intake at baseline (week 0) and postintervention (week 8) for whey protein, calcium-caseinate and control interventions
| Whey protein | Calcium caseinate | Control | ||||||||
| Baseline | Post | Δ | Baseline | Post | Δ | Baseline | Post | Δ | ||
| Energy, MJ/d | 9.5 ± 0.5 | 9.4 ± 0.5 | −0.09 ± 0.4 | 9.6 ± 0.5 | 9.2 ± 0.4 | −0.46 ± 0.3 | 9.2 ± 0.5 | 9.7 ± 0.5 | 0.50 ± 0.4 | 0.149 |
| Total fat, % of TE | 34.6 ± 1.0 | 32.2 ± 1.2 | −2.4 ± 0.8 | 34.8 ± 1.2 | 31.8 ± 1.1 | −3.0 ± 0.9 | 35.2 ± 1.2 | 32.3 ± 1.2 | −2.9 ± 1.1 | 0.906 |
| SFA, % of TE | 12.9 ± 0.6 | 12.7 ± 0.8 | −0.2 ± 0.6 | 12.4 ± 0.7 | 11.7 ± 0.7 | −0.7 ± 0.4 | 13.4 ± 0.6 | 12.2 ± 0.7 | −1.2 ± 0.6 | 0.391 |
| MUFA, % of TE | 11.7 ± 0.5 | 10.8 ± 0.4 | −1.0 ± 0.3 | 12.0 ± 0.5 | 10.9 ± 0.4 | −1.1 ± 0.4 | 12.3 ± 0.6 | 11.0 ± 0.5 | −1.3 ± 0.5 | 0.981 |
| n–6 PUFA, % of TE | 4.81 ± 0.2 | 4.11 ± 0.2 | −0.70 ± 0.2 | 5.22 ± 0.4 | 4.23 ± 0.2 | −0.99 ± 0.3 | 4.58 ± 0.3 | 4.37 ± 0.3 | −0.21 ± 0.3 | 0.464 |
| n–3 PUFA, % of TE | 0.78 ± 0.1 | 0.79 ± 0.1 | 0.01 ± 0.6 | 0.84 ± 0.1 | 0.85 ± 0.1 | 0.01 ± 0.6 | 0.74 ± 0.1 | 0.71 ± 0.1 | −0.03 ± 0.6 | 0.495 |
| Total PUFA, % of TE | 5.59 ± 0.3 | 5.04 ± 0.3 | −0.55 ± 0.4 | 6.06 ± 0.4 | 4.81 ± 0.3 | −1.24 ± 0.4 | 5.32 ± 0.3 | 5.08 ± 0.3 | −0.24 ± 0.4 | 0.456 |
| 0.88 ± 0.1 | 0.90 ± 0.1 | 0.02 ± 0.1 | 0.97 ± 0.1 | 0.95 ± 0.1 | −0.02 ± 0.1 | 1.00 ± 0.1 | 0.97 ± 0.1 | −0.03 ± 0.1 | 0.972 | |
| Protein, % of TE | 16.0 ± 0.5 | 24.4 ± 0.9 | 8.4 ± 0.8a | 15.9 ± 0.5 | 24.4 ± 0.7 | 8.5 ± 0.6a | 16.2 ± 0.6 | 15.2 ± 0.6 | −1.00 ± 0.6b | <0.001 |
| Carbohydrates, % of TE | 49.4 ± 1.3 | 44.0 ± 1.2 | −5.5 ± 1.0a | 49.7 ± 1.4 | 43.7 ± 1.4 | −6.0 ± 1.3a | 48.5 ± 1.7 | 54.8 ± 1.6 | 6.3 ± 1.9b | <0.001 |
| Alcohol, % of TE | 3.02 ± 0.7 | 2.60 ± 0.5 | −0.42 ± 0.5 | 2.59 ± 0.6 | 2.76 ± 0.6 | 0.18 ± 0.4 | 3.20 ± 0.7 | 2.48 ± 0.7 | −0.72 ± 0.6 | 0.066 |
| Dietary fiber (AOAC), g/d | 22.6 ± 1.4 | 20.1 ± 1.2 | −2.49 ± 1.1 | 23.0 ± 1.4 | 20.5 ± 1.2 | −2.56 ± 0.8 | 21.3 ± 1.1 | 20.4 ± 1.3 | −0.92 ± 0.9 | 0.399 |
| Sodium, g/d | 3.18 ± 0.2 | 3.12 ± 0.2 | −0.07 ± 0.2 | 3.37 ± 0.2 | 2.81 ± 0.2 | −0.57 ± 0.2 | 3.57 ± 0.5 | 2.95 ± 0.2 | −0.62 ± 0.5 | 0.215 |
| Potassium, g/d | 5.22 ± 1.2 | 5.21 ± 1.2 | −0.01 ± 0.2 | 4.54 ± 0.7 | 4.71 ± 0.9 | 0.17 ± 0.4 | 4.10 ± 0.6 | 4.80 ± 0.9 | 0.70 ± 0.8 | 0.582 |
| Calcium, g/d | 1.06 ± 0.1 | 1.19 ± 0.1 | 0.13 ± 0.1b | 1.09 ± 0.1 | 1.57 ± 0.1 | 0.48 ± 0.1a | 0.97 ± 0.1 | 1.00 ± 0.1 | 0.03 ± 0.1b | <0.001 |
All values are means ± SEMs. n = 38. Data were log transformed for alcohol intake. Baselines were not significantly different (P ≤ 0.05) from one another. Different superscript letters within a row refer to treatment groups different from one another (P ≤ 0.05). AOAC, Association of Official Analytic Chemists; TE, total energy; Δ change from baseline.
Overall between-group treatment effects for each Δ were obtained with the use of a linear mixed-model ANOVA with baseline values for the variable of interest and prognostic values such as age, sex, and BMI. Tukey-Kramer correction was used for the post hoc analysis to adjust for multiple testing.