|
|
||||||||
ORIGINAL RESEARCH COMMUNICATION |
1 From the Jean Mayer US Department of Agriculture Human Nutrition Research Center on Aging at Tufts University, Boston (PKN and KLT), and the National Institute on Aging, National Institutes of Health, Baltimore (DM, JH, and RA)
2 Supported in part by the US Department of Agriculture, Agricultural Research Service, contract number 53-3K06-01; the National Institutes of Health, National Institute on Aging Intramural Program; and the General Mills Bell Institute of Health and Nutrition.
3 Reprints not available. Address correspondence to PK Newby, Jean Mayer USDA Human Nutrition Research Center on Aging at Tufts University, 711 Washington Street, 9th Floor, Boston, MA 02111. E-mail: pknewby{at}post.harvard.edu.
| ABSTRACT |
|---|
|
|
|---|
Objective: Our objective was to examine whether food patterns derived from exploratory factor analysis are related to anthropometric changes. We hypothesized that we would derive a healthy food pattern and that it would predict smaller changes in body mass index (BMI; in kg/m2) and waist circumference (in cm) than would other food patterns in models adjusted for baseline anthropometric measures.
Design: The subjects were 459 healthy men and women participating in the Baltimore Longitudinal Study of Aging. Diet was assessed by using 7-d dietary records, from which 40 food groups were formed and entered into a factor analysis.
Results: Six food patterns were derived. Factor 1 (reduced-fat dairy products, fruit, and fiber) was most strongly associated with fiber (r = 0.39) and loaded heavily on reduced-fat dairy products, cereal, and fruit and loaded moderately on fruit juice, nonwhite bread, nuts and seeds, whole grains, and beans and legumes. In a multivariate-adjusted model in which the highest and lowest quintiles were compared, factor 1 was inversely associated with annual change in BMI (ß = 0.51; 95% CI: 0.82, 0.20; P < 0.05; P for trend < 0.01) in women and inversely associated with annual change in waist circumference (ß = 1.06 cm; 95% CI: 1.88, 0.24 cm; P < 0.05; P for trend = 0.04) in both sexes.
Conclusion: Our results suggest that a pattern rich in reduced-fat dairy products and high-fiber foods may lead to smaller gains in BMI in women and smaller gains in waist circumference in both women and men.
Key Words: Food patterns dietary patterns factor analysis principal components obesity body composition diet assessment body mass index BMI waist circumference
| INTRODUCTION |
|---|
|
|
|---|
25], and 31% of these adults are obese (BMI
30) (1). Obesity results from weight gains over time; however, dietary determinants of weight gain remain controversial (2, 3). Inconsistent findings may be explained in part by limitations in the single-nutrient approach, which is commonly used in nutritional epidemiologic research (4, 5). Cluster and factor analyses have emerged as methods to empirically derive dietary (or food) patterns and possibly explain nutrient-disease relations. In our previous study, we found that patterns derived by using cluster analysis were significantly related to longitudinal changes in anthropometric measures (6). Specifically, the healthy dietary patternwhich was high in fruit, reduced-fat dairy products, high-fiber cereal, and whole grains and low in red and processed meats, fast food, and sodashowed smaller gains in both BMI and waist circumference than did the other dietary patterns (white bread, sweets, alcohol, and meat and potatoes).
In the current study, we continued to examine prospectively whether eating patterns are related to anthropometric changes in women and men participating in the Baltimore Longitudinal Study of Aging (BLSA) by using a distinctly different patterning methodfactor analysis. Specifically, we used principal components to extract factors, in which a given set of variables is transformed into a reduced set of correlated variables (7). Factor analysis is statistically quite different from cluster analysis. Whereas cluster analysis separates persons into mutually exclusive groups based on differences in mean food intakes (clusters), factor analysis separates foods into groups based on correlations between foods (factors), and persons receive a factor score for each of the derived factors. With the use of the factor analysis procedure, a persons dietary pattern would be best represented by looking at his or her factor scores for each of the derived vectors. Because factor analysis studies focus on individual food vectors when examining relations rather than the total diet for a given person, we refer in this article to the derived factors as food patterns and not as dietary patterns.
Other studies that have used factor or principal components analysis have empirically derived a "healthy" food pattern that is high in fruit, vegetables, seafood, fiber, and other "healthy" foods (8-11). None of these studies, nor other studies using factor analysis (12-14), examined whether food patterns derived by using factor analysis are related prospectively to anthropometric changes. Our objective was to examine whether food patterns derived from exploratory factor analysis, specifically, principal components analysis, are related to changes in BMI and waist circumference. We hypothesized that, using factor analysis, we would derive a healthy food pattern and that this pattern would lead to smaller changes in BMI and waist circumference than would other dietary patterns.
| SUBJECTS AND METHODS |
|---|
|
|
|---|
This analysis uses the same data set as described in our previous paper (6) to facilitate comparison of the factor and cluster analyses in deriving food and dietary patterns, respectively. In summary, our study population was limited to those entering the study on or after 1980 (n = 921) to avoid bias from changing dietary trends over time, and persons aged <30 y and >80 y were excluded (n = 138). Of the remaining 783 subjects, persons who had not completed
4 d of dietary records were omitted (n = 9), as were subjects whose food group intake appeared implausible (>6 SDs from the mean for each food group) (n = 63). An additional 134 subjects were omitted because they did not have
2 measures of height or weight (at the time the dietary record was collected and at follow-up). Finally, all subjects who received a diagnosis of cancer, diabetes, stroke, or heart disease either before or at baseline were omitted (n = 190) to create a disease-free cohort. After all exclusions, 459 subjects were available for the analysis. An additional 10 subjects were omitted from the waist circumference analysis because they did not have baseline or follow-up information on waist circumference; thus, data for 449 subjects were used in this analysis.
Dietary assessment
Dietary intake was assessed with the use of 7-d dietary records; reports detailing dietary collection methods and dietary intake in the BLSA population were published previously (16-18). In summary, study participants were instructed by trained dietitians to complete 7-d food records. Food records were completed at home by the participants and were sent back to the study center. Before 1992, the subjects were given food models and a booklet of food pictures to help them assess portion size. In 1994, the subjects were given a portable scale to weigh food portions. The participants were contacted by telephone with any questions about the diet records.
Dietary records from 1984 to 1991 were originally coded and entered into a nutrient database maintained by the BLSA, whereas diet records completed since 1994 were coded and entered into the Minnesota Nutrition Data System (Nutrition Coordinating Center, University of Minnesota, Minneapolis) at Tufts University; no dietary data were collected in 1992 and 1993. Dietary data from 1984 to 1991 were then reentered into the Minnesota Nutrition Data System, and nutrient intakes were back-adjusted to correct for changes in the food supply (eg, changes in nutrient content resulting from cereal fortification) by using data from the US Department of Agriculture to correspond to appropriate time intervals (19, 20).
Food pattern derivation
To perform a food pattern analysis with the use of principal components analysis, individual foods and food ingredients from the dietary records were first aggregated into groups. We formed 40 food groups, mainly according to macronutrient composition (eg, fat or fiber content) and culinary use; several foods (eg, pizza and eggs) composed their own groups (Appendix A
). Where possible, foods were separated into full- and reduced-fat groups (eg, high-fat dairy and reduced-fat dairy products).
|
The PROC FACTOR procedure in SAS (version 8.2; SAS Institute, Cary, NC) was used to perform the analysis. The procedure uses principal components analysis and orthogonal rotation (varimax option in SAS) to derive noncorrelated factors. As with other empirically derived pattern procedures, principal components analysis requires preselection of the number of factors. To decide which number of factors to retain, we ran factor solutions ranging from 2 to 15, including food-group components with an eigenvalue > 1 and examined both the scree plots and the factors themselves to see which set of factors most meaningfully described distinct food patterns. From these analyses, a 6-factor solution was selected. Factor loadings were calculated for each food group across the 6 factors (food patterns). A factor score was then calculated for each subject for each of the 6 factors, in which the standardized intakes of each of the 40 food groups were weighted by their factor loadings and summed (7); the sums were then standardized (
= 0, SD = 1). Food patterns were named according to both the foods that loaded most positively on the factor and according to how the factors correlated with nutrients (see Statistical analyses).
Some studies have found that factor solutions differ by sex (22). We examined factor solutions separately for men and women but chose to retain the factor solution that included both men and women to maximize statistical power. We also tested for an interaction between food pattern and sex in the regression analyses.
Anthropometric and covariate assessment
Anthropometric measurements were made following standardized procedures (23) and were fully described elsewhere (24). In summary, weight and height were measured for each subject at each visit, from which the BMI was calculated. Waist circumference was measured at each visit with an inelastic tape at the narrowest part of the torso at the end of expiration (24).
Demographic data were collected from each study participant at the first visit and were used to adjust for potential confounding in regression analyses. Race-ethnicity, physical activity, smoking habit, education, and vitamin supplement use were determined by questionnaire at the time that dietary records were collected. Physical activity was measured by an adapted version of the Harvard Alumni questionnaire, which asked participants about all daily activities (eg, activities at home, work, and during recreation or sports). The amount of time spent for each activity was summed across all activities to determine the daily energy output per body weight (kJ/kg) and was described previously (16, 25).
Statistical analyses
Factors were divided into quintiles, and sample means and frequencies were calculated at the baseline visit at which time dietary data were collected. The percentage of subjects who were overweight (BMI = 2529.99) or obese (BMI
30) was calculated by using recommended international cutoffs (26). Partial correlation coefficients (adjusted for age, sex, and energy intake) were calculated between each food pattern vector and selected macronutrients to further describe the patterns. Mean nutrient and food group intakes were calculated for each quintile of factor 1 (reduced-fat dairy products, fruit, and fiber) because it was the dominant food pattern vector.
Separate regression analyses were performed for each factor to test whether food patterns predicted changes in BMI and waist circumference. Factors were treated categorically (quintiles), in which the lowest quintile was the reference group; a test for trend was also performed by using a categorical variable that ranked the median factor score for each quintile. Annual change in BMI was calculated by subtracting BMI at visit 1 from BMI at visit 2, dividing by the time interval between visits (months), and then multiplying by 12 mo. Annual change in waist circumference was calculated by using the same algorithm. We tested the a priori hypotheses that the healthy food pattern we derived would predict smaller changes in BMI and waist circumference than would other food patterns in models adjusted for baseline anthropometric measures (baseline BMI and baseline waist circumference), age, sex, and sociodemographic covariates. Because total energy intake is the mechanism through which nutrients and foods may affect BMI, it is arguable whether energy should be included in the model; we fit a final model for each analysis in which we added total energy to the multivariate model to see whether adjustment for baseline energy would change our estimates. We also added terms to adjust for development of disease (stroke, cancer, diabetes, and heart disease) during the follow-up period. Because age may not be linearly related to BMI, we fit regression models using quadratic terms for age. We tested for interactions between food pattern and sex by using a cross-product term.
In secondary analyses, we performed 2 tests to compare the factor solution discussed in this study with the cluster solution reported in our previous study (6). In the first test, we calculated the mean factor scores for each cluster by using the same 40 food groups used in the current study. In the second test, we used Pitmans test to compare the SD of the residuals from the each of the regression analyses to examine which patterning method better predicted change in BMI and waist circumference. All analyses were performed by using the Statistical Analysis System (SAS) for WINDOWS (version 8.2; SAS Institute, Cary, NC).
| RESULTS |
|---|
|
|
|---|
|
|
|
|
Regression results associating change in BMI across quintiles for all 6 factors are presented in Table 5
. A test for interaction between food pattern and sex indicated that the effect of food pattern on change in BMI was significantly modified by sex for factor 1 (P = 0.008), but not for other factors (P > 0.05); therefore, the results are presented separately for women and men for factor 1 only. For women, the highest quintile of factor 1 was inversely associated with annual change in BMI in a comparison with the lowest quintile (ß = 0.51; 95% CI: 0.82, 0.20; P < 0.05) and the test for trend was significant (P < 0.01) in a multivariate-adjusted model. Factor 1 was inversely but not significantly related to annual change in BMI in men in a comparison of the fifth with the first quintile (ß = 0.14; 95% CI: 0.36, 0.08; P > 0.05; P for trend = 0.47). Factor 2 was directly associated with change in BMI in a comparison of the fifth with the first quintile (ß = 0.20; 95% CI: 0.04, 0.36; P < 0.05), and the test for trend was marginally significant (P = 0.05). Factors 46 were not significantly related to annual change in BMI (data not shown).
|
|
± SD: 0.98 ± 0.08), and the lowest mean score for factor 1 was seen for those subjects in the meat-and-potatoes cluster (
± SD: 0.48 ± 0.09) after adjustment for additional risk factors (data not shown). Pitmans test showed no significant difference between the factor solution and the cluster solution in predicting change in BMI (P = 0.26); however, the factor solution better predicted change in waist circumference than did the cluster solution (P = 0.03) (data not shown). | DISCUSSION |
|---|
|
|
|---|
The healthy pattern we observed is most similar to the cold-foods pattern derived in the study by Maskarinec et al (14), which had high factor loadings for cold breakfast cereals, fruit juice, and fruit. It also shares common elements with the cereals patterns derived in several other studies (21, 27, 28). Factor 1 is similar to the "healthy" cluster derived in our previous study, whereas the protein and alcohol and sweets factors are less similar to the alcohol and sweets clusters, respectively. Although these similarities point to some comparability in the methods, the extent to which factors and clusters yield different solutions is attributable to methodologic differences.
Our healthy pattern is less similar to the "prudent" pattern described in other studies (10, 29, 30), which is high in vegetables, fruit, legumes, and fish. Also, we did not observe a separate "Western" pattern, which is high in processed meat, refined grains, sweets and desserts, and French fries (10, 29). Rather, we derived a pattern, factor 1, that had high positive factor loadings for healthy foods and high negative factor loadings for less healthy foods. In our study, then, the "Western" pattern previously described (10, 29) appears embedded in factor 1, ie, persons with higher scores for factor 1 were consuming a healthy food pattern, whereas those with lower scores for factor 1 were consuming a Western food pattern.
Different effects of food patterns on BMI by sex were previously reported (10, 14, 27, 29, 31, 32), and an interaction between food patterns and sex was also observed in relation to coronary heart disease (33). In this study, of those in the top quintile of factor 1, 59% were vitamin users and 87% were never smokers, which suggests that a healthy diet is part of a healthy lifestyle (34). Although we adjusted for these factors in our analysis, there may be other potential confounders associated with dietary intake and change in BMI that may be different in men and women, for which we were not able to adjust. Our finding of an interaction between a fiber-rich food pattern and sex may reflect a true association or may be due to lack of power, lack of variation in factor 1 scores for men, different overlap in factor scores between men and women, uncontrolled confounding, measurement error, or some other reason. More research is needed to elucidate whether inconsistent findings between women and men are due to methodologic limitations or whether the biological effects of diet are truly modified by sex.
We are not aware of any studies that have examined the relation between food patterns derived by using factor analysis and waist circumference. The inverse relation seen here between factor 1 and anthropometric changes may be explained by the roles of fiber (35), glycemic index (36), and energy density (37-39) on food intake and appetite, as previously discussed (6); the roles of dairy products and calcium; or some othre reason. For example, the dietary fiber intake was >10 g greater in the highest quintile of factor 1 than in the lowest quintile of factor 1, with similar intakes in quintiles 24, which possibly explained why only the highest quintile was associated with changes in waist circumference and BMI. Understanding the reasons for the associations between the other food patterns we derived and anthropometric changes, however, is more difficult. Factor 2 was most strongly correlated with alcohol, and factor 3 was most strongly correlated with sucrose; these nutrients may help to explain the effects we observed. Because a persons total dietary pattern is represented by his or her scores on all of the derived factors, an examination of individual factor relations may obfuscate diet-disease relations within a person. Perhaps, because of this limitation, investigators have created mutually exclusive groups of persons based on their factor scores across several food patterns [eg, the healthiest pattern combined those in the highest quintile of the prudent pattern and the lowest quintile of the Western pattern (29). This approach, however, is not practical for a 6-factor solution. We therefore focused our analysis and discussion on the hypothesized healthy food pattern.
Our derived food patterns explained almost 30% of the variance in dietary intake, which is less than that observed in some studies (10, 29, 30) but similar in magnitude to that observed in other studies (22, 27). The amount of variance explained by individual factors may be expected to vary across studies, because it is a function of both the number of variables included in the analysis (7) and the correlation matrix itself (28). When we compared the factor solution from this study with the cluster solution from our previous study (6), no significant difference between solutions was observed in predicting change in BMI, but the factor solution presented here better predicted change in waist circumference. More research is needed to further examine whether significant differences in patterning solutions are consistently observed in explaining health outcomes.
Our study has several limitations. Our sample has limited variability in ethnicity, and education and food patterns differed according to these variables (4, 40). However, it is not clear whether the patterns we derived would have a different biologic effect in persons with different educational levels or ethnicities. Food patterns may also differ by sex (4, 22, 41, 42). Women tend to have healthier diets than do men (43, 44), which may explain why 60% of those in the top quintile of factor 1, our healthy pattern in this study, were women. In addition, we were unable to adjust for the role of genetics in our models, which may modify the relation between diet and adiposity (45).
The principal components method itself also has limitations that stem from several subjective decisions that investigators must make. First, one must decide how to treat the dietary input variables and whether these variables should be adjusted for total energy intake, because these decisions will affect the correlations between the variables and, hence, the factor solution. In this study, we chose to consider the dietary variables as a percentage of energy to obtain factors proportional to total energy intake. However, some reports have adjusted for intakes by using the nutrient residual approach (46) and others have used total grams (47) or servings (8, 10) unadjusted for total energy; other studies have transformed food group variables logarithmically to improve normality (14). Research is needed to examine how the treatment of dietary variables affects the factor solution, which may lead to insight on what is the best method for handling variables in nutritional studies using pattern analysis.
Another limitation of factor analysis is that the investigator is forced to prespecify the number of factors. Although eigenvalues, scree plots, and interpretability (7) are used to guide the investigator in determining the best factor solution, ultimately such a decision is subjective. The number of factors extracted by other studies has varied, ranging from 2 (48, 49) to 10 (50). The aforementioned reasons for the increasing complexity of analyzing >2 factors may explain why investigators using factor analysis are inclined to derive or report less, rather than more, food patterns. Our choice to retain 6 factors was based on eigenvalues, scree plots, interpretability, and similarity in range to other studies, in addition to our a priori belief that there are likely >2 meaningful food patterns in nature.
In summary, we derived a healthy food pattern and observed prospective associations with this pattern and BMI in women and with waist circumference in both sexes. Food patterns differ between populations, including by age (50), religion (51), and ethnicity (49), but the prudent, Western, and cereal-based patterns as well as the sweets and alcohol patterns seem to be reasonably reproducible among populations. Our results suggest that a fiber-rich food pattern that is high in reduced-fat dairy products, cereal, fruit, fruit juice, nonwhite bread, nuts and seed, whole grains, and beans and legumes may lead to smaller gains in BMI in women and smaller gains in waist circumference in both women and men.
| ACKNOWLEDGMENTS |
|---|
RA and JH contributed to the original design and to data collection for the BLSA. DM prepared the nondietary data for this analysis. PKN was responsible for the design and analysis of the study and drafted the manuscript. KLT contributed to the analysis and design of the study and oversaw the Tufts collaboration with the BLSA. All authors made critical comments during the preparation of the manuscript and fully accept responsibility for the work. No author had a financial interest or professional or personal affiliations that compromised the scientific integrity of this work.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
L. Johnson, A. P Mander, L. R Jones, P. M Emmett, and S. A Jebb Energy-dense, low-fiber, high-fat dietary pattern is associated with increased fatness in childhood Am. J. Clinical Nutrition, April 1, 2008; 87(4): 846 - 854. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. K. Good, N. Holschuh, A. M. Albertson, and A. L. Eldridge Whole Grain Consumption and Body Mass Index in Adult Women: An Analysis of NHANES 1999-2000 and the USDA Pyramid Servings Database J. Am. Coll. Nutr., February 1, 2008; 27(1): 80 - 87. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Esmaillzadeh and L. Azadbakht Major Dietary Patterns in Relation to General Obesity and Central Adiposity among Iranian Women J. Nutr., February 1, 2008; 138(2): 358 - 363. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Murakami, S. Sasaki, Y. Takahashi, K. Uenishi, M. Yamasaki, H. Hayabuchi, T. Goda, J. Oka, K. Baba, K. Ohki, et al. Hardness (difficulty of chewing) of the habitual diet in relation to body mass index and waist circumference in free-living Japanese women aged 18-22 y Am. J. Clinical Nutrition, July 1, 2007; 86(1): 206 - 213. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Esmaillzadeh, M. Kimiagar, Y. Mehrabi, L. Azadbakht, F. B Hu, and W. C Willett Dietary patterns, insulin resistance, and prevalence of the metabolic syndrome in women Am. J. Clinical Nutrition, March 1, 2007; 85(3): 910 - 918. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. McNaughton, G. D. Mishra, A. M. Stephen, and M. E. J. Wadsworth Dietary Patterns Throughout Adult Life Are Associated with Body Mass Index, Waist Circumference, Blood Pressure, and Red Cell Folate J. Nutr., January 1, 2007; 137(1): 99 - 105. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Mendez, B. M. Popkin, P. Jakszyn, A. Berenguer, M. J. Tormo, M. J. Sanchez, J. R. Quiros, G. Pera, C. Navarro, C. Martinez, et al. Adherence to a Mediterranean Diet Is Associated with Reduced 3-Year Incidence of Obesity J. Nutr., November 1, 2006; 136(11): 2934 - 2938. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. K. Newby, C. Weismayer, A. Akesson, K. L. Tucker, and A. Wolk Longitudinal Changes in Food Patterns Predict Changes in Weight and Body Mass Index and the Effects Are Greatest in Obese Women J. Nutr., October 1, 2006; 136(10): 2580 - 2587. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Okubo, S. Sasaki, H. Horiguchi, E. Oguma, K. Miyamoto, Y. Hosoi, M.-k. Kim, and F. Kayama Dietary patterns associated with bone mineral density in premenopausal Japanese farmwomen Am. J. Clinical Nutrition, May 1, 2006; 83(5): 1185 - 1192. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. K. Newby, C. Weismayer, A. Akesson, K. L. Tucker, and A. Wolk Long-Term Stability of Food Patterns Identified by Use of Factor Analysis among Swedish Women J. Nutr., March 1, 2006; 136(3): 626 - 633. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Schulz, U. Nothlings, K. Hoffmann, M. M. Bergmann, and H. Boeing Identification of a Food Pattern Characterized by High-Fiber and Low-Fat Food Choices Associated with Low Prospective Weight Change in the EPIC-Potsdam Cohort J. Nutr., May 1, 2005; 135(5): 1183 - 1189. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |