Endothelial dysfunction and arterial stiffness are early predictors of atherosclerosis, hypertension and CVD( Reference Vlachopoulos, Aznaouridis and Stefanadis 1 , Reference Yeboah, Folsom and Burke 2 ). There is strong evidence that diet is related to endothelial dysfunction( Reference Brown and Hu 3 , Reference West 4 ) and, to a lesser degree, arterial stiffness( Reference Aatola, Koivistoinen and Hutri-Kahonen 5 , Reference Kesse-Guyot, Vergnaud and Fezeu 6 ). However, there is a significant variation in the methods used to quantify diet in previous studies, with the majority examining the intake of specific foods or nutrients rather than overall diet. Studies of chronic disease morbidity and mortality have indicated that the use of dietary indices, or diet quality scores, is a comprehensive approach that can provide a valuable insight into the relationship between diet and health( Reference McCullough, Feskanich and Stampfer 7 , Reference McNaughton, Bates and Mishra 8 ).
The Dietary Guidelines for Americans (DGA) are evidence-based recommendations that provide guidance for choosing an eating pattern that promotes health and prevents disease. The 2010 Guidelines emphasise greater intake of fruits, vegetables, low-fat dairy products, whole grains, and a variety of lean meats while maintaining appropriate weight through energy balance and physical activity( 9 ). The Dietary Guidelines Adherence Index (DGAI) is a tool that quantifies the degree to which key DGA recommendations are met. Developed in reference to the 2005 DGA( Reference Fogli-Cawley, Dwyer and Saltzman 10 ) and updated for the 2010 DGA( Reference Troy and Jacques 11 ), the DGAI provides an objective index of diet quality that is useful for standardising dietary assessments across studies. To our knowledge, no studies have evaluated whether overall diet quality is associated with measures of vascular function, particularly in a large community-based sample.
Vascular health declines with age despite the control of traditional risk factors. It is unclear whether age-related decline in vascular function is part of a normal physiological ageing process or a consequence of repeated exposure to lifestyle-related risk factors. Physiological changes with age probably interact with lifestyle risk factors to exacerbate arterial stiffness and endothelial dysfunction( Reference Lakatta and Levy 12 ). Given the burden of CVD on the USA's ageing population, there is a need for improved understanding of the interaction between age and lifestyle and its effect on vascular function.
Therefore, the purpose of the present study was to determine whether adherence to the 2010 DGA is associated with endothelial dysfunction and arterial stiffness in a cross-sectional sample of adults from the Framingham Heart Study. A secondary purpose was to determine whether age influences the association between diet quality and these measures of vascular health.
Experimental methods
Subjects
The Framingham Heart Study is a longitudinal, community-based study of risk factors for CVD. The present study includes dietary and vascular data collected during the seventh examination cycle of the Offspring cohort (1998–2001( Reference Feinleib, Kannel and Garrison 13 )) and the first examination cycle of the Third Generation cohort (2002–2005( Reference Splansky, Corey and Yang 14 )). The sample characteristics are presented in Table 1. The present analysis was approved by the Institutional Review Board at the Pennsylvania State University.
DGAI, Dietary Guidelines Adherence Index; bpm, beats per min.
* DGAI-2010 range was 0–100 possible points.
Dietary measurements
The Harvard semi-quantitative FFQ( Reference Fung, Rexrode and Mantzoros 15 ) was mailed to participants before the examination, and they were asked to bring the completed form to their appointment. The 126-item questionnaire assesses the consumption frequency of standard servings of foods and beverages during the last year with response selections ranging from ‘never or less than once per month’ to ‘6+ per d.’ The Harvard FFQ provides a space for participants to write-in up to three additional foods they frequently consumed that were not listed, and specifically asks for the type of breakfast cereal and cooking oil regularly used. Nutrient intakes are calculated by multiplying average intake with nutrient content of individual foods, based on the US Department of Agriculture food composition database and supplemented with other sources( Reference Rimm, Giovannucci and Stampfer 16 ).
The DGAI-2010 was applied to the FFQ data to determine the extent to which participants' diets are consistent with the 2010 DGA (see the online Supplementary material for further description and example calculation). The DGAI-2010 assessed the intake of fourteen food groups (fruit; dark green vegetables; orange and red vegetables; starchy vegetables; other vegetables; grains; milk; meat, protein, and eggs; seafood; nuts; legumes; sugar; variety in protein choices; and variety of fruits and vegetables) and eleven healthy choice or nutrient intake recommendations (amounts of total fat, saturated fat, trans-fat, cholesterol, Na, fibre and alcohol, and percentage of lean protein, low-fat milk, whole grains and whole fruits). Adherence to each DGAI-2010 item is scored on a continuous scale of 0–1, and the categories were summed and standardised to a range of 0–100 to create an overall score, with higher scores indicating greater adherence. An important component of the DGAI compared with other dietary quality assessment tools is the penalty assigned for overconsumption, which is in line with the 2005 and 2010 DGA emphasis on weight management. In other words, the DGAI avoids assigning a higher score to individuals who meet the recommended food intakes simply through eating more. Appropriate energy levels were calculated for each participant (based on height, weight, age, sex and physical activity estimates) and participants were penalised for consuming more than the recommended daily intake of energy-dense foods (e.g. starchy vegetables, specific protein sources, grains, meat and beans, and dairy products) for their energy intake.
Vascular measurements
Endothelial function was assessed by brachial artery flow-mediated dilation (FMD). Methodology and reproducibility data have been published previously( Reference Benjamin, Larson and Keyes 17 , Reference Hamburg, Palmisano and Larson 18 ). Briefly, brachial artery diameter (mm) was imaged in the supine position with high-resolution ultrasound at rest and 1 min after reactive hyperaemia that was induced by the 5 min cuff occlusion of forearm blood flow. Arterial diameter was measured offline using commercially available edge-detection software. Brachial FMD was calculated as the percentage change in brachial diameter during reactive hyperaemia from the resting state (%FMD), with lower values indicating greater endothelial dysfunction. Baseline and post-deflation hyperaemic flow velocity were assessed with Doppler imaging at baseline and for 15 s immediately post-deflation, as described previously( Reference Mitchell, Parise and Vita 19 ).
Central (aortic) arterial stiffness was assessed in the supine position with arterial tonometry, as described previously( Reference Mitchell, Guo and Benjamin 20 ). Briefly, blood pressure was obtained with an oscillometric (Offspring) or auscultatory (Third Generation) device, and mean arterial pressure was measured via brachial waveform planimetry. A tonometer recorded blood pulsations at the right carotid, brachial, radial and femoral arteries. Transit distances were measured from the suprasternal notch to each recording site. Tonometry waveforms were signal-averaged offline and calibrated using cuff pressures, as described previously. Carotid–femoral pulse wave velocity (PWV) was calculated from transit distances and tonometry waveforms, as described previously( Reference Mitchell, Parise and Benjamin 21 ), with greater PWV indicating greater arterial stiffness. The augmentation index was calculated from the carotid pressure waveform, as described previously( Reference Murgo, Westerhof and Giolma 22 ), with higher values reflecting greater relative wave reflection.
Covariates
Potential confounders of the relationship between diet and vascular health were considered in the present analysis in accordance with previous studies( Reference Benjamin, Larson and Keyes 17 , Reference Mitchell, Guo and Benjamin 20 ). All participants underwent routine medical examination at the time of vascular assessment to obtain the following characteristics: age; sex; race; BMI; heart rate; fasting glucose; total:HDL-cholesterol ratio; TAG; diabetes (defined as a fasting blood glucose of ≥ 7 mmol/l ( ≥ 126 mg/dl) or treatment with insulin or an oral hypoglycaemic agent); hypertension (defined as a systolic blood pressure of ≥ 140 mmHg and a diastolic blood pressure of ≥ 90 mmHg); or existing CVD (CHD, heart failure, stroke, transient ischaemic attack or intermittent claudication). Systolic and diastolic blood pressures were the average of two physician-measured readings at the Heart Study. Hormone replacement therapy, hypertension medication, lipid-lowering medication and cigarette smoking status (in the 6 h before vascular testing) were determined by self-report. A variable representing the timing of a walk test (performed concomitantly at Offspring Exam 7) in relation to the vascular assessments (before v. after or not done) was included. We also included variables denoting family relatedness (parent–child and sibling–sibling) and cohort.
Statistical analyses
Of the 7634 participants who attended the seventh Offspring exam (n 3539) or the first Third Generation exam (n 4095), 5887 had complete dietary and covariate data. Of these, brachial FMD data were available for 5521, flow data were available for 5067 and tonometry data were available for 5379. To maximise power, participants were included in the analyses for which complete data were available. To determine power for the present analysis, we reviewed an earlier Framingham Heart Study analysis of brachial FMD where a final model including eight predictors yielded a multiple R 2 of 0·16 for %FMD( Reference Benjamin, Larson and Keyes 17 ). In the present study, the sample size of 5521 (for brachial FMD data) provided >90 % power with an α of 0·05 to detect a change in the model R 2 of 0·01.
All analyses were conducted in SAS version 9.3 (SAS Institute, Inc.). DGAI-2010 scores were divided into equal quintiles according to the full sample (n 5887 total, n 1174 or 1175 per quintile). Means and 95 % CI of participant characteristics and potential covariates across quintile categories, adjusted for age and sex, were computed using general linear models. The statistical significance for trend was assessed using linear regression for continuous variables with the DGAI-2010 entered as a continuous score.
The DGAI-2010 score and all vascular outcome variables were assessed for normality; baseline flow velocity and PWV were positively skewed. A natural log transformation was applied to baseline flow velocity and an inverse transformation to PWV (1000/PWV). Quintile category means and 95 % CI of vascular characteristics, adjusted for clinical covariates (see below), were computed using general linear models. Analysis of the residual plots indicated that the assumption of linearity was met. The statistical significance for trend was assessed with the DGAI-2010 entered as a continuous score, and the generalised estimating equations approach was applied to account for the familial correlations in the present sample. First-order interactions between the DGAI-2010 and age were assessed for each of the vascular characteristics using model 2 (described below); variables with statistically significant interactions were stratified ( < 50 or ≥ 50 years) for further investigation.
For all vascular outcomes, two analyses were performed with family relatedness and cohort indicator variables included as covariates in all models. Model 1 adjusted for age and sex, and model 2 additionally adjusted for relevant clinical covariates (BMI, mean arterial pressure, heart rate and smoking status)( Reference Benjamin, Larson and Keyes 17 , Reference Mitchell, Guo and Benjamin 20 ). We explored the effect of further adjusting for total:HDL-cholesterol, TAG, diabetes, hypertension therapy, lipid therapy, hormone replacement therapy and prevalent CVD, and completing the walk test before vascular testing in a third model; however, this analysis yielded the same results as model 2 and is therefore not presented. For all analyses, P< 0·05 was considered statistically significant. Unless otherwise noted, results are presented as adjusted means and 95 % CI.
Results
The sample characteristics stratified by sex are presented in Table 1. The sample was 54 % women, with an average age of approximately 48 years for both men and women. The mean DGAI-2010 score was 55 for men and 61 for women. On average, both men and women were overweight, but men tended to have a worse metabolic profile and a higher use of anti-hypertensive and lipid-lowering medications. Increasing DGAI-2010 scores were significantly associated with increasing age (P< 0·001) and decreasing BMI (P< 0·001), heart rate (P< 0·001), total:HDL-cholesterol (P< 0·001), TAG (P< 0·001) and glucose (P< 0·001), and were significantly higher among women (P< 0·001) and non-smokers (P< 0·001) (data not shown).
The vascular characteristics according to the quintile categories of the DGAI-2010 are reported in Table 2 (model 1) and Table 3 (model 2). Baseline brachial artery diameter and FMD were not significantly associated with DGAI-2010 scores in model 1 or 2. Baseline mean flow velocity was lower with higher DGAI-2010 scores in both models. Surprisingly, hyperaemic mean flow velocity was lower with higher DGAI-2010 scores in model 1, though this association was blunted in the fully adjusted model. Further analysis indicated that concurrent adjustment for heart rate, BMI and smoking status (but not mean arterial pressure) attenuated the association between hyperaemic mean flow velocity and diet, with the greatest attenuation observed when smoking status was added to the model. Mean arterial pressure and carotid–femoral PWV were lower with higher dietary quintile scores in model 1, but the relationships were attenuated in model 2; further analysis indicated that adjustment for heart rate alone rendered the associations non-significant. The augmentation index was lower with increasing DGAI-2010 scores in both models.
* Derived with general linear models adjusted for age, sex, cohort and family relatedness.
† Derived from general estimating equations with the DGAI-2010 entered as a continuous score.
* Derived with general linear models adjusted for age, sex, cohort, family relatedness, BMI, mean arterial pressure, heart rate and smoking status.
† Derived from general estimating equations with the DGAI-2010 entered as a continuous score.
We tested the interactions between the DGAI-2010 and age for vascular characteristics using model 2, and found a significant interaction for mean arterial pressure, carotid–femoral PWV and augmentation index. Stratified analyses ( < 50 or ≥ 50 years; Table 4) indicated that mean arterial pressure is lower with higher DGAI-2010 scores in younger adults (β = − 0·03, P= 0·05), but not in older adults (β = 0·04, P= 0·09). Similarly, stratified analyses suggested that carotid–femoral PWV is lower with higher DGAI-2010 scores in younger adults (β = − 0·03, P= 0·01), but not in older adults (β = 0·001, P= 0·06), although neither association was statistically significant. The augmentation index in the younger group was significantly lower with higher DGAI-2010 scores (β = − 0·05, P= 0·01); although a similar association was indicated in the older group, it did not reach statistical significance (β = − 0·04, P= 0·06).
* Derived with general linear models adjusted for age (continuous), sex, BMI, mean arterial pressure, heart rate and smoking status.
† Derived from general estimating equations with the DGAI-2010 entered as a continuous score.
Discussion
In a large cross-sectional community-based cohort study, we have comprehensively evaluated the associations of adherence to the 2010 DGA with measures of vascular function. Vasodilator measures in both a conduit artery, assessed by brachial FMD, and the microvessels, assessed by reactive hyperaemia, were not associated with dietary adherence. Resting brachial flow velocity, but not diameter, was related to dietary adherence. The association of central aortic stiffness with diet in unadjusted models appeared to be related to concomitant risk factors. However, wave reflection assessed by the augmentation index was lower with greater dietary adherence, an association that was more pronounced in adults younger than 50 years.
The cross-sectional relationships between selected dietary components and FMD was previously examined in over 3000 adults in the Multi-Ethnic Study of Atherosclerosis cohort, and found that among women (but not men), regular fish intake was associated with higher FMD( Reference Anderson, Nettleton and Herrington 23 ); however, fish intake was the only component of diet reported. Numerous clinical trials have reported beneficial effects of dietary interventions on FMD, such as interventions low in fat( Reference Fuentes, Lopez-Miranda and Sanchez 24 – Reference Koh, Son and Ahn 27 ), rich in unsaturated fat( Reference Fuentes, Lopez-Miranda and Sanchez 24 , Reference West, Krick and Klein 28 , Reference Fuentes, Lopez-Miranda and Perez-Martinez 29 ), based on the Mediterranean diet( Reference Cuevas, Guasch and Castillo 30 – Reference Ryan, McInerney and Owens 33 ), or rich in protein( Reference Ferrara, Innelli and Palmieri 34 ). Additionally, a review of observational studies has concluded that diets rich in fruits and vegetables are inversely associated with biomarkers of endothelial dysfunction (such as cellular adhesion molecules and other pro-inflammatory markers), whereas Westernised diets rich in meat are positively associated with biomarkers of endothelial dysfunction( Reference Oude Griep, Wang and Chan 35 ). In the present study, we found that a dietary pattern in line with the 2010 DGA was not related to baseline brachial diameter and FMD. It is possible that the food groups and nutrients highlighted by the DGA are not those most important to vascular function, at least when assessed by brachial FMD, as these guidelines were meant to promote general health rather than prevent a specific condition such as vascular disease. Thus, the use of an overall index may be masking the effects of specific foods and nutrients, including those previously shown to modify endothelial function and arterial stiffness (e.g. nuts, chocolate, tea, red wine, n-3 fatty acids and Na)( Reference Brown and Hu 3 – Reference Kesse-Guyot, Vergnaud and Fezeu 6 , Reference van Trijp, Beulens and Bos 36 – Reference Crichton, Elias and Dore 42 ). Importantly, as the 2010 DGA index does not include a component specific to intake of fish rich in long-chain n-3 fatty acids or overall PUFA consumption, we are unable to compare our findings with those reported in the Multi-Ethnic Study of Atherosclerosis study described above. The differences between the present results and previous intervention trials may be explained by the limitations of cross-sectional observational studies and FFQ in assessing diet. Short-term intervention studies that provide food to participants can more accurately measure consumption of a particular food or dietary pattern, and thereby establish efficacy in modifying endothelial function.
Brachial flow velocities at rest and during hyperaemia reflect arterial properties in the microcirculation. In the present analysis, we have shown that increased adherence to the 2010 DGA is associated with lower baseline (resting) flow velocity. In the fully adjusted model, we observed a difference in mean baseline flow velocity between the bottom and top quintiles of diet scores of − 0·5 cm/s. Prior studies in the present cohort and others have demonstrated the associations between higher resting flow velocity and CVD risk factors (particularly metabolic risk factors)( Reference Hamburg, Larson and Vita 43 ), and there is evidence that higher resting flow may induce small-vessel damage( Reference Mitchell, Vita and Larson 44 ). The absolute difference in resting flow that we observed between quintile 1 and quintile 5 (0·5 cm/s) is similar in magnitude to the 0·39 cm/s increase predicted by every increase of 1·3 in total:HDL-cholesterol ratio and to the 0·75 cm/s increase predicted by every increase of 4·6 kg/m2 in BMI in a prior analysis of the Framingham Heart Study( Reference Mitchell, Vita and Larson 44 ). Taken together, the present results suggest that adherence to the 2010 DGA may be as important as other CVD risk factors in determining resting flow velocity.
Hyperaemic flow reflects small-vessel vasodilation in response to ischaemia, and also predicts CVD outcomes and correlates with CVD risk factors( Reference Anderson, Charbonneau and Title 45 – Reference Philpott, Lonn and Title 47 ). We found an unexpected trend towards a negative association for adherence to the DGA and hyperaemic flow velocity in the age- and sex-adjusted model that was blunted in the fully adjusted model. Further analysis indicated that heart rate, BMI and smoking status accounted for the association of DGAI-2010 scores and hyperaemic flow velocity.
Prior observational studies have indicated that diets rich in meat intake and high alcohol consumption are associated with greater arterial stiffness( Reference Kesse-Guyot, Vergnaud and Fezeu 6 , Reference van Trijp, Beulens and Bos 36 ), whereas diets with moderate alcohol consumption( Reference Sierksma, Lebrun and van der Schouw 37 – Reference van den Elzen, Sierksma and Oren 40 ), low Na intake( Reference Avolio, Clyde and Beard 41 ), greater fruit and vegetable consumption( Reference Aatola, Koivistoinen and Hutri-Kahonen 5 ), and greater consumption of dairy products( Reference Crichton, Elias and Dore 42 , Reference Livingstone, Lovegrove and Cockcroft 48 ) have been associated with lower arterial stiffness. In the present study, adherence to the DGA was related to mean arterial pressure, carotid–femoral PWV and augmentation index in the age- and sex-adjusted model, but only the augmentation index remained significantly associated with the DGAI-2010 after further adjustment for CVD risk factors. On average, the difference in the augmentation index (%) between the bottom and top quintiles of dietary scores was 1·3 %, which was similar to the increase of 0·93 % predicted by every 8·5 year increase in age within the present cohort( Reference Mitchell, Parise and Benjamin 21 ). This finding is consistent with reduced wave reflection and ventricular ejection( Reference Torjesen, Wang and Larson 49 ) with greater adherence to the DGA. Further analyses indicated that the relationship between PWV and diet observed in the age- and sex-adjusted model was no longer evident after adjustment for heart rate. Heart rate is an important potential confounder of associations with carotid–femoral PWV( Reference Mitchell, Guo and Benjamin 20 ), and researchers are encouraged to adjust for this in future studies.
There was a significant interaction between the DGAI-2010 and age for vascular stiffness measures that persisted after adjustment for heart rate and the other covariates in model 2. Age is the predominant risk factor for CVD( Reference Lakatta and Levy 12 ), and advancing age increases the risk despite the control of modifiable lifestyle factors( Reference Roger, Go and Lloyd-Jones 50 – Reference Wilson, D'Agostino and Levy 53 ). Stratified analyses indicated that mean arterial pressure and arterial wave reflection were lower with higher DGAI-2010 scores in adults younger than 50 years, but in those aged 50 years and older, the associations were not strong or were statistically non-significant. While stratified analyses for PWV were non-significant for both age groups, the trend towards lower PWV with higher DGAI-2010 scores was notably stronger in the younger group. Collectively, the present results indicate that for younger adults, following a diet that more closely resembles the 2010 DGA is associated with better vascular health. In contrast, for older adults, adherence to the 2010 DGA is unrelated to vascular health. Longitudinal studies and intervention studies with long-term follow-up are needed to understand the possible dietary contribution to vascular decline.
The goal of the present study was to examine the association between adherence to the DGA and vascular health. However, there may be limitations to this approach. The DGA are evidence-based recommendations that provide guidance for choosing an eating pattern that promotes health and prevents disease, but as noted above, these recommendations do not focus solely on vascular disease. Moreover, few individuals in this cohort consumed diets that closely adhered to the DGA, which may limit our ability to see benefits of this dietary pattern. Other limitations of the study include its cross-sectional nature that prevents us from drawing conclusions about causation and related mechanisms. The Framingham cohorts are overwhelmingly white; thus, generalisation to other races or ethnicities is limited. However, the use of this large well-characterised sample enables us to examine the relationship between diet and vascular health with consideration of CVD risk factors. In addition, the age range of this sample (19–89 years) allowed us to examine the relationship between diet quality and vascular health over a wide age range.
In conclusion, we have shown that adherence to the 2010 DGA is associated with measures of blood flow velocity and arterial wave reflection, but not related to brachial FMD. Importantly, we have demonstrated that diet may be particularly related to vascular health in adults younger than 50 years. Future studies should examine whether interventions that increase adherence to the DGA modify vascular health, especially among younger adults.
Supplementary material
To view supplementary material for this article, please visit http://dx.doi.org/10.1017/S0007114515000859
Acknowledgements
The present study was supported by the National Institutes of Health (to K. A. S., grant no. F31AG043224 and T32DK07658, Penn State grant no. UL1TR000127, and Framingham Heart Study grant no. HL076784, HL070100, HL060040, HL080124, HL071039, HL077447 and 2-K24-HL04334). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH. The Framingham Heart Study was conducted and supported by the National Heart, Lung, and Blood Institute in collaboration with Boston University (contract no. N01-HC-25915). This manuscript was reviewed by the participating Framingham Heart Study investigators for scientific content and consistency of data interpretation with previous Framingham Heart Study publications.
The authors' contributions are as follows: K. A. S. conducted the analysis and drafted the manuscript; D. N. P., M. C., P. F. J., L. M. T., N. W., N. M. H., J. A. V., R. S. V., G. F. M. and S. G. W. assisted in the creation, design, analysis and interpretation of the project. All authors critically revised and approved the final manuscript.
The authors declare that there are no conflicts of interest.