|
|
||||||||
ORIGINAL RESEARCH COMMUNICATION |
1 From the Division of Epidemiology and Community Health, School of Public Health, University of Minnesota, Minneapolis, MN (JAN, LMS, and DRJ); the Department of Nutrition, University of Oslo, Oslo, Norway (DRJ); the Department of Epidemiology, German Institute of Human Nutrition Potsdam-Rehbruecke, Nuthetal, Germany (MBS); the Department of Pathology, University of Vermont, Burlington, VT (NSJ); the Department of Medicine, College of Physicians and Surgeons, Columbia University, New York, NY (RGB); and the Department of Public Health Sciences, Wake Forest University, Winston-Salem, NC (AGB)
2 Supported by training grant T32 HL07779 (to JAN) and contracts N01-HC-95159 through N01-HC-95166 from the National Heart, Lung, and Blood Institute, and General Clinical Research Center grant M01-RR00645 from the National Center for Research Resources. 3 Reprints not available. Address correspondence to JA Nettleton, University of Minnesota, Division of Epidemiology and Community Health, School of Public Health, 1300 South 2nd Street, Suite 300, Minneapolis, MN 55454. E-mail: nettleton{at}epi.umn.edu.
| ABSTRACT |
|---|
|
|
|---|
Objective: We investigated whether empirically derived dietary patterns are associated with coronary artery calcium (CAC) and common and internal carotid artery intima media thickness (IMT) and whether prior information about inflammatory processes would increase the strength of the associations.
Design: At baseline, dietary patterns were derived with the use of a food-frequency questionnaire, and inflammatory biomarkers, CAC, and IMT were measured in 5089 participants aged 4584 y, who had no clinical CVD or diabetes, in the Multi-Ethnic Study of Atherosclerosis. Dietary patterns based on variations in C-reactive protein, interleukin-6, homocysteine, and fibrinogen concentrations were created with reduced rank regression (RRR). Dietary patterns based on variations in food group intake were created with principal components analysis (PCA).
Results: The primary RRR (RRR 1) and PCA (PCA factor 1) dietary patterns were high in total and saturated fat and low in fiber and micronutrients. However, the food sources of these nutrients differed between the dietary patterns. RRR 1 was positively associated with CAC [Agatston score >0: OR (95% CI) for quartile 5 compared with quartile 1 = 1.34 (1.05, 1.71); ln(Agatston score = 1): P for trend = 0.023] and with common carotid IMT [
1.0 mm: OR (95% CI) for quartile 5 compared with quartile 1 = 1.33 (0.99, 1.79); ln(common carotid IMT): P for trend = 0.006]. PCA 1 was not associated with CAC or IMT.
Conclusion: The results suggest that subtle differences in dietary pattern composition, realized by incorporating measures of inflammatory processes, affect associations with markers of subclinical atherosclerosis.
Key Words: Dietary patterns principal components analysis reduced rank regression carotid artery intima media thickness coronary artery calcium
| INTRODUCTION |
|---|
|
|
|---|
Dietary pattern analysis facilitates the examination of food and nutrient synergies in relation to disease, whereas traditional reductive approaches may not capture such synergies and associations (30, 31). Recently, Hoffmann et al (32) proposed the application of reduced rank regression (RRR) in dietary pattern analysis. In contrast with the commonly used principal components analysis (PCA), which produces a linear combination of food groups (a factor) that maximally explains variation in food group intake, RRR produces a linear combination of food groups that maximally explains variation in investigator-specified intermediate response variables. Thus, the technique allows hypotheses regarding pathways (represented by the response variables) between diet and disease to be evaluated (33). Although PCA has been successful, it may not achieve the full potential of dietary pattern analysis because it relies solely on intercorrelations among dietary variables, which may not optimally represent the diet qualities most relevant to specific disease etiology (32).
In the present study we used RRR to derive dietary patterns on the basis of their ability to maximally explain variation in biochemical markers related to inflammation (ie, CRP, IL-6, fibrinogen, and homocysteine). A comparison of the food group composition of behavior-based PCA patterns with the food group composition of inflammation-based RRR patterns will provide an interesting perspective on the value of the 2 methods, ie how much information is lost if patterns are based on behavior alone and do not take into account inflammatory processes? On the basis of previous successes of PCA in this (26) and other (34) cohorts, we hypothesized that both PCA and RRR techniques would produce dietary patterns associated with CAC and carotid artery IMT. However, by incorporating inflammatory intermediates in the RRR analysis, we hypothesized that the RRR method would produce dietary patterns more strongly associated with CAC and IMT than dietary patterns derived by PCA.
| SUBJECTS AND METHODS |
|---|
|
|
|---|
The current cross-sectional investigation includes data from 5089 participants enrolled in MESA, including 2407 men and 2682 women aged 4584 y. This sample was chosen after the exclusion of individuals who were currently taking oral steroid or antiinflammatory asthma medications (n = 134), who had diabetes mellitus defined according to the American Diabetes Association 2003 criteria (n = 919) (36), and who provided insufficient or implausible dietary information (n = 630: 167 because of extreme energy intake and 463 because of implausible responses or improperly completed questionnaires; 26).
Carotid intima media thickness and coronary calcium
Thickness of the intima media layers of the common and internal carotid arteries was determined from images obtained by high-resolution B-mode ultrasonography (Logiq 700 ultrasound machine; GE Medical Systems, Waukesha, WI). Calculation of IMT was performed at the central MESA ultrasound reading center (TuftsNew England Medical Center, Boston, MA). In the present study, high IMT was defined as measures above the 75th percentile of the sample distribution, separately for each internal carotid IMT (1.3 mm) and common carotid IMT (1.0 mm).
Computed tomography (CT) of the coronary arteries was performed with cardiac-gated (at 80% of the R-R interval) electron beam scanners at 3 centers (Imatron C-150; Imatron Inc, San Francisco, CA) (37) and with a prospective electrocardiogram-triggered scan acquisition at 50% of the R-R interval with multidetector scanners (38) at the remaining 3 centers. These scanners are comparable in their ability to measure calcium (39). All scans were read at the central MESA CT reading center (Harbor-University of California Medical Center; Los Angeles, CA). CAC scores (Agatston scores) were determined by blinded CT analysts. Presence of CAC was defined as an Agatston score >0.
CRP, fibrinogen, interleukin-6, and homocysteine
CRP, fibrinogen, IL-6, and homocysteine concentrations were measured in blood samples collected at baseline, processed (40), and stored at 80 °C until analyzed. CRP and fibrinogen antigen were measured in plasma with a BNII nephelometer (N High Sensitivity CRP and N Antiserum to Human Fibrinogen; Dade Behring Inc, Deerfield, IL). IL-6 concentrations were measured by ultrasensitive enzyme-linked immunosorbent assay (Qantikine HS Human IL-6 Immunoassay; R&D Systems, Minneapolis, MN). Plasma homocysteine was measured by polarization immunoassay with an IMx Analyzer (IMx Homocysteine Assay; Axis Biochemicals ASA; Oslo, Norway). Analytic CVs were 3.6%, 2.6%, 6.3%, and 4.5% for CRP, fibrinogen antigen, IL-6, and homocysteine, respectively.
Dietary assessment
At the baseline examination, each participant completed a self-administered, modified-Block style food-frequency questionnaire (FFQ) (26, 41, 42). The FFQ was provided in English, Spanish, or Chinese with staff assistance available on request. Responses were converted to servings/d and categorized into 47 food groups.
Dietary pattern analysis
Dietary patterns were derived from 47 food groups with the use of 2 different empirical methods: PCA (26) and RRR (32) (SAS PROC FACTOR with varimax rotation and SAS PROC PLS with the RRR method option, respectively). In the RRR analysis, the specified response variables were CRP, IL-6, fibrinogen, and homocysteine (natural log values due to non-normal distribution of original values). For the PCA analysis, a 4-pattern solution was chosen after the evaluation of a scree plot of the eigenvalues and after the interpretability of the final solution was considered. For the RRR analysis, a 4-pattern solution was used in accordance with the number of response variables specified (32). The 4-factor PCA solution explained 26.2% of the variation in food group intake and 1.5% of the variation in response variables. The 4-factor RRR solution explained 12.6% of the variation in food group intake and 5.0% of the variation in response variables. To present data concisely and because other dietary patterns were not significantly associated with the outcomes, the results are presented for only the primary factor (dietary pattern) derived by the PCA and the primary factor derived by RRR.
A score was calculated for each participant for each PCA and RRR dietary pattern as a sum of the 47 food groups, each weighted according to the factor loadings. Participant scores were then categorized into quintiles separately for each dietary pattern. Thus, for each dietary pattern, quintile 5 was composed of persons whose diets conformed most closely to that particular pattern.
Assessment of other relevant variables
Information on demographics, education, medication use, smoking history, and physical activity were collected at the baseline examination with a combination of self-administered and interviewer-administered questionnaires. Waist circumference (cm) was measured at the umbilicus, and body mass index (BMI; in kg/m2) was calculated from measured height and weight.
Statistical analysis
All analyses were conducted with SAS version 9.1 (SAS Institute Inc, Cary, NC). The percentage variance in total dietary pattern score explained by each food group was calculated as follows (32):
![]() | (1) |
Odds ratios (ORs) for risk of high internal (
1.3 mm) and common (
1.0 mm) carotid IMT and CAC presence (Agatston score >0) were calculated for quintiles 25 of each dietary pattern by using quintile 1 as the reference category. The PCA-derived and RRR-derived dietary patterns were each entered separately into the multivariable logistic regression model, and the P for trend across quintiles was calculated by treating the categorical dietary pattern variable as a continuous term in all analyses. A linear regression analysis was also performed with the continuous natural log-transformed values of mean IMT of the common carotid artery, internal carotid artery, and CAC score (adding 1 to each before transformation to accommodate zeroes and very small values) as the dependent variables and the dietary patterns (quintile variables) as the independent variables. Again, each dietary pattern was entered into the regression model separately, with adjustment for the appropriate covariates (defined below).
In addition to age- and energy-adjusted analyses (model 1), we adjusted for the following confounding variables: field center (Baltimore City and County, MD; Chicago, IL; Forsyth County, NC; New York, NY; Los Angeles County, CA; and St Paul, MN), sex (male or female), race-ethnicity (white, black, Hispanic, and Chinese), education level (less than high school, high school, more than high school), active leisure activity (metabolic equivalents in minutes per week, includes moderate and vigorous activities), inactive leisure activity (metabolic equivalent in minutes per week, including activities such as television watching and automobile transportation), smoking (never versus current and pack-years), and dietary supplement use (at least weekly use versus nonuse) (model 2). We further investigated potential mediating pathways using the following 2 models: 1) main multivariable model listed above + waist circumference, and 2) main multivariable model + CVD risk factors with dietary origins (systolic blood pressure, HDL cholesterol, LDL cholesterol, triacylglycerols, insulin, and glucose). Additional adjustment for statin or hypertension medication use did not appreciably affect risk estimates; therefore, these variables were not included. Furthermore, exclusion of statin users (n = 656) before analysis did not materially affect results. We did not consider excluding users of hypertension medications because of the large numbers (n = 1646).
To assess homogeneity of risk estimates, stratified analyses were also conducted for BMI (<25 or
25) and waist circumference (<88.9 or
88.9 cm for women; <101.6 or
101.6 cm for men). Formal tests of interaction between dietary patterns and race-ethnicity, sex, current smoking status, BMI, or waist circumference were assessed by adding cross product terms to the main multivariable model.
To assess whether the associations ascribed to a dietary pattern were driven by single food group components, we performed analyses similar to those described above with the individual food groups contributing most to dietary pattern score variance as the independent variables (categorical and continuous analyses).
| RESULTS |
|---|
|
|
|---|
|
|
|
Dietary patterns and CRP, IL-6, fibrinogen, and homocysteine concentrations
Both PCA 1 and RRR 1 were significantly and positively correlated with CRP, IL-6, fibrinogen, and homocysteine (Table 1
), although the magnitudes of the correlations with RRR 1 were approximately twice as large as those with PCA 1. The percentages of total biomarker variance (RRR response variables) explained by each PCA 1 (0.23%) and RRR 1 (3.0%) also differed to a similar degree. PCA 1 explained 0.55%, 0.36%, 0.002%, and 0.01% of the variation in CRP, IL-6, fibrinogen, and homocysteine, respectively, whereas for RRR 1, these values were considerably higher (5.3%, 4.0%, 2.4%, and 0.34%, respectively).
Dietary patterns and participant characteristics
The distribution of participant characteristics across quintiles of PCA 1 was different from that of RRR 1 (Table 2
). For example, educational attainment was positively associated with PCA 1 scores, whereas educational attainment was negatively associated with RRR 1 scores. The directions of the associations for age, sex, and percentage Hispanic also differed between the 2 patterns (P for trend < 0.001 for all). However, both patterns were associated positively with inactive leisure, waist circumference, BMI, and smoking prevalence and negatively with supplement use (P for trend < 0.001 for all).
|
34% higher than for persons in the first quintile (OR: 1.34; 95% CI:.05, 1.71). Formal tests of interaction between dietary patterns and race, sex, current smoking status, waist circumference, or BMI were not statistically significant. The estimate for risk of high common carotid IMT based on RRR 1 scores in all participants was similar in magnitude to that of CAC > 0 (OR: 1.33), but the 95% CIs were wide (0.99, 1.79). If potential mechanistic variables were added to the main multivariable model, risk estimates for CAC > 0 were attenuatedmodestly after adjustment for waist circumference [quintile 5 versus quintile 1 (OR: 1.28; 95% CI: 1.00, 1.64)] and more notably after adjustment for systolic blood pressure, HDL cholesterol, LDL cholesterol, triacylglycerols, insulin, and glucose [quintile 5 versus quintile 1 (OR: 1.19; 95% CI: 0.93, 1.53)], which suggests that any relations between this particular dietary pattern and CAC are mediated through the effects of diet on traditional CVD risk factors other than inflammatory processes. Results of regression analysis of the continuous variables [mean internal IMT (mm), common IMT (mm), and CAC (Agatston score)] across dietary pattern quintiles were similar to those for the dichotomized analyses. Mean CAC (Agatston score) and mean common carotid IMT increased across quintiles of RRR 1 (P for trend = 0.023 and 0.006 for CAC and common carotid IMT, respectively), although the range of values across quintiles was rather narrow (
2% increase in common carotid IMT and
29% increase in CAC; Table 4
|
|
| DISCUSSION |
|---|
|
|
|---|
The associations noted in relation to our RRR dietary pattern are similar to those reported in the one other study that investigated the relation between empirically derived dietary patterns and IMT (20). In a cluster analysis in a cohort of white women, Millen et al (20) reported that women falling into dietary pattern clusters that were high in total and saturated fats and low in micronutrients and fiber had an increased risk of carotid artery atherosclerosis relative to women in a low-fat, high-fiber, and micronutrient-rich "heart healthy" dietary pattern cluster. Because our RRR dietary pattern was low in fiber and antioxidant vitamins but high in saturated fat, our findings are also consistent with the studies reporting the beneficial (14, 17, 18) and deleterious (saturated fat) (16) effects of these nutrients in relation to IMT, CAC, or both. Although our PCA dietary pattern also shared similar characteristics, it was not as discriminating as was the RRR dietary pattern in terms of inflammatory biomarker concentrations, which possibly explained the absence of associations between PCA 1 and IMT or CAC.
RRR applied to the nutritional discipline is currently fairly novel. One of our goals was to add to this smaller body of research by comparing the associations between IMT or CAC and an RRR dietary pattern with the associations between IMT or CAC and a dietary pattern derived by a more common method, PCA. We hoped that such a comparison would increase knowledge of the RRR method and also further our understanding of diet-disease relations. We found the RRR dietary pattern to be associated with CAC and common carotid IMT, whereas the PCA dietary pattern was not significantly associated with these markers. A comparison of these 2 dietary patterns may aid in identifying specific food groups with biological effects relevant to inflammation and the development of atherosclerosis. Several nutrients showed similar correlations with RRR 1 and PCA 1 (saturated fat, trans fat, fiber, and antioxidant vitamins), and intakes of the food groups fats and oils and processed meats were important determinants of both RRR 1 and PCA 1 dietary pattern scores. However, the consumption of cruciferous and dark-yellow vegetables, soy foods, and soda notably differed between RRR 1 and PCA 1 (as well as other food groups not among the highest score contributors; Appendixes A and B). It is possible that such subtle differences in dietary pattern composition are biologically relevant to the atherosclerotic process. For example, Adams et al (43) recently showed that diets high in green and yellow vegetables (freeze-dried peas, green beans, broccoli, corn, and carrots) inhibited aortic atherosclerosis in a mouse model of atherosclerosis. Soy protein-feeding studies in monkeys fed otherwise atherogenic diets showed reduced atherosclerosis compared with milk proteins (44, 45), which may have been due partly to modified monocyte-endothelial cell interactions in the early stages of atherosclerosis (46). Although soda consumption has not been evaluated previously in relation to atherosclerosis, its association with obesity and type 2 diabetes is well known (47, 48); thus, a negative effect on atherosclerosis would not be unexpected. However, in the present study, individual food groups were not associated with IMT or CAC, which suggests that all facets of the dietary pattern act synergistically. One important role of dietary pattern analysis is to guide further investigations at the food and nutrient level (30). Although it is likely that the effects of single constituents are partly dependent on the overall dietary context, randomized clinical studies are needed to disentangle the effects of single components from those of the larger dietary pattern.
The limitations of this study deserve discussion. First, this cross-sectional analysis may be confounded by recent dietary changes that are not reflective of lifetime dietary practices, which are influential in the development of atherosclerosis. Similarly, these cross-sectional data do not allow for mechanistic hypotheses to be tested, ie, "does diet impact the atherosclerotic process by influencing inflammatory processes? " We can only conjecture that all factors (diet, inflammation, and atherosclerotic disease) may be interrelated. Second, as a result of the number of statistical tests conducted, noted associations may have been significant by chance alone. Indeed, the ORs were small but were robust to adjustment for multiple confounders and were consistent across strata of several characteristics. Third, we caution that ORs overestimate the relative prevalence when outcomes are common, as was the case for CAC presence (
48% of cohort). Fourth, it is possible that the collection of biomarkers used in our RRR procedure did not adequately represent inflammatory processes relevant in the development of atherosclerosis. It is possible that the use of different biomarkers may have produced larger risk estimates, such as those reported in relation to coronary artery disease when more traditional CVD risk factors were used (33). However, we chose these biomarkers because of general interest in novel disease pathways. Last, the relation between diet and atherosclerosis may be mediated through pathways other than inflammation, and likewise, the relation between diet and CVD events may be mediated through pathways other than those reflected by carotid artery IMT and CAC (49).
Because RRR patterns are not based on dietary behavior, the applicative value of the resulting patterns in real-world settings remains largely unknown. In contrast, PCA dietary patterns reflect real-world food usage patterns but do not take into account biological measures of health (eg, inflammatory markers). Nevertheless, PCA-derived dietary patterns have consistently shown relevant associations with many CVD-related outcomes (34), including subclinical biomarkers whose concentrations are typically unknown by participants (25, 26). The unresolved question then becomes whether it is more important to identify a dietary pattern that reflects how people really eat and how that pattern is associated with disease outcomes, or alternatively, to identify an artificial dietary pattern (collection of foods) that is associated with disease outcomes. Likely, a valuable role for each method exists, depending on the question at handeg, PCA for planning behavioral interventions and RRR for etiologic investigations.
Dietary pattern research continues to be an area of interest because of the potential for food and nutrient synergies (30), utility of dietary patterns for making public health recommendations (50), and values of dietary pattern change (51-54). Our study is unique in that we were able to study a large multiethnic cohort, with a well-characterized measure of overall dietary intake and validated measures of inflammation and subclinical atherosclerosis. We found that a dietary pattern derived by RRR and based on variations in CRP, IL-6, fibrinogen, and homocysteine was positively associated with CAC and common carotid IMT, but a dietary pattern derived by PCA was not significantly associated with CAC or carotid artery IMT. Food group differences between these patterns did not solely explain the difference in relations to outcomes, which suggests the importance of the overall dietary pattern and method of dietary pattern derivation.
| ACKNOWLEDGMENTS |
|---|
The authors' responsibilities were as followsJAN: analytic design, data analysis, and manuscript preparation; LMS: critical review of the manuscript; MBS: data analysis and critical review of the manuscript; NSJ, RGB, and AGB: data acquisition and manuscript review; DRJ: data acquisition, critical review of the manuscript, and data analysis. None of the authors had a conflict of interest to report.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Franco, A. V Diez-Roux, J. A Nettleton, M. Lazo, F. Brancati, B. Caballero, T. Glass, and L. V Moore Availability of healthy foods and dietary patterns: the Multi-Ethnic Study of Atherosclerosis Am. J. Clinical Nutrition, March 1, 2009; 89(3): 897 - 904. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. DiBello, P. Kraft, S. T. McGarvey, R. Goldberg, H. Campos, and A. Baylin Comparison of 3 Methods for Identifying Dietary Patterns Associated With Risk of Disease Am. J. Epidemiol., December 15, 2008; 168(12): 1433 - 1443. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. M Hendricks, D M. Mwamburi, P. Newby, and C. A Wanke Dietary patterns and health and nutrition outcomes in men living with HIV infection Am. J. Clinical Nutrition, December 1, 2008; 88(6): 1584 - 1592. [Abstract] [Full Text] [PDF] |
||||
![]() |
U. Nothlings, S. P Murphy, L. R Wilkens, H. Boeing, M. B Schulze, H B. Bueno-de-Mesquita, D. S Michaud, A. Roddam, S. Rohrmann, A. Tjonneland, et al. A food pattern that is predictive of flavonol intake and risk of pancreatic cancer Am. J. Clinical Nutrition, December 1, 2008; 88(6): 1653 - 1662. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |