|
|
||||||||
ORIGINAL RESEARCH COMMUNICATION |
1 From the Department of Epidemiology, German Institute of Human Nutrition Potsdam, Rehbrücke, Germany (KH and HB), and the Department of Internal Medicine, University Hospital Hamburg-Eppendorf, Germany (BCZ and EW)
2 The funding for the CORA Study came primarily from the German Stifterverband für die Wissenschaft. 3 Address reprint requests and correspondence to K Hoffmann, Department of Epidemiology, German Institute of Human Nutrition, Arthur-Scheunert-Allee 114-116, 14558 Bergholz-Rehbrücke, Germany. E-mail: khoff{at}mail.dife.de.
| ABSTRACT |
|---|
|
|
|---|
Objective: The objective was to apply a new statistical method to construct a specific dietary pattern that is strongly associated with the risk of coronary artery disease (CAD).
Design: We applied reduced rank regression to a sample of 200 cases and 255 controls from the Coronary Risk Factors for Atherosclerosis in Women (CORA) Study. The CAD-specific dietary pattern was constructed by choosing intake data for 49 food groups as predictors and 5 established biomarkers for CAD as responses.
Results: A high score for the constructed dietary pattern was characterized by high intakes of meat, margarine, poultry, and sauce and low intakes of vegetarian dishes, wine, vegetables, and whole-grain cereals. After adjustment for known CAD risk factors, the relative risks from the lowest to the highest quintiles of the pattern score were 1.0, 1.1, 3.6, 6.2, and 12.3 (95% CI: 4.9, 30.9; P for trend < 0.0001). There was an approximate 4.5-fold difference in C-reactive protein and a 2-fold difference in C-peptide between the highest and lowest score quintiles of the study population. HDL-cholesterol concentrations ranged from 70 mg/dL in the lowest quintile to 49 mg/dL in the highest quintile of dietary pattern score.
Conclusion: The new statistical method, reduced rank regression, may be a useful tool for identifying dietary patterns that simultaneously affect the concentrations of known CAD biomarkers and the risk of developing CAD.
Key Words: Coronary artery disease risk factors nutrition dietary patterns biomarkers statistical methods
| INTRODUCTION |
|---|
|
|
|---|
A possible approach to derive dietary patterns consists of applying statistical dimension-reduction techniques to the full data set of all food intakes (13). Two well-known dimension-reduction techniques similar in their mathematical foundation are principal component analysis and factor analysis. They aid in the construction of linear functions of food intakes that explain as much variation of the original intake data as possible. Although both statistical methods have been widely applied to derive dietary patterns in epidemiologic studies of CAD (1417), neither of them aims to find out what combination of consumed foods mostly increases the risk of CAD. A dietary pattern explaining a high percentage of food consumption variation is not necessarily associated with the development of CAD, as can be seen from the nonsignificant results in the literature (16). Thus, another statistical method appropriate for constructing disease-specific dietary patterns is required.
Reduced rank regression (RRR) is a new dimension-reduction technique that can fill the gap. In contrast with the exploratory methods principal component analysis and factor analysis, RRR can allow for prior information about the pathway from diet to disease. Applying RRR to food groups and choosing biochemical markers for CAD as response variables, linear functions of food groups are derived by maximizing the proportion of explained biomarker variation. This approach ensures that a change in RRR pattern score results in a change in CAD biomarker concentrations, which is generally associated with a decrease or increase in CAD risk. We applied RRR to food consumption data collected through a food-frequency questionnaire (FFQ) in the Coronary Risk Factors for Atherosclerosis in Women (CORA) Study.
| SUBJECTS AND METHODS |
|---|
|
|
|---|
For each patient, 2 controls of the same age and from the same city district of Hamburg were randomly chosen on the basis of the population registry. They were invited by mail to participate in the study. The invitation was followed by a reminder or a telephone call after 14 d. Exclusion criteria for the controls were the same as for cases. The final control group comprised 255 women, which corresponded to a participation rate of 67% of those who were eligible.
Data collection
Dietary intake information was collected by a self-administered scanner-readable FFQ originally designed for and applied in the German part of the EPIC (European Prospective Investigation into Cancer and Nutrition) Study (18). The FFQ included questions on the frequency and portion size of 146 food items eaten during the previous year. Foods were classified into 49 food groups based on nutrient profiles or culinary usage (19). Total caloric intake was calculated by summing up energy intakes from all foods by using data from the German Food Code BLSII.2. A validation study among members of the EPIC-Potsdam cohort showed good correlation between nutrients assessed by this FFQ and multiple 24-h dietary recalls. Energy-adjusted correlation coefficients ranged from 0.54 for dietary fiber to 0.86 for alcohol (20).
Information about sociodemographic characteristics and suspected risk factors for CAD, including smoking habits, physical activity, hormone replacement therapy, medication use, and family history of cardiovascular disease were obtained by questionnaire and by conducting computer-assisted person-to-person interviews. Smokers were defined as cigarette smokers and were categorized into those subjects who smoked < and
20 cigarettes/d. Former smokers were categorized into those subjects who stopped smoking within the past
and >2y, because previous studies have indicated that much of the coronary disease risk attributable to smoking disappears within 2 y of quitting (21, 22). Women were defined as postmenopausal if they had no regular monthly period for >1 y or were receiving hormone replacement therapy.
Anthropometric measurements including body weight, body height, waist circumference, and hip circumference were obtained while the subjects wore no shoes and only light underwear. Body weight was measured to the nearest 0.5 kg. The waist was taken at the smallest circumference between the lower rib margin and the iliac crest, and the hip circumference was determined over the greater trochanters. All interviews and physical examinations were performed by the same trained investigator.
A fasting blood sample was collected as soon as possible, within at least 24 h in women who had experienced an acute myocardial infarction. Routine laboratory measurements were performed with the use of standard techniques. LDL cholesterol was calculated by using the formula of Friedewald. Blood pressure was measured 3 times after the interviews while the subjects were in a sitting position. The results of the second and third measurements were averaged (23). Subjects were considered to have hypertension if they were taking antihypertensive drugs or had a systolic blood pressure
140 mm Hg or a diastolic blood pressure of
90 mm Hg according to the definition of the World Health Organization (24). All participants obtained an evaluation of their dietary habits and their laboratory results. The study protocol was approved by the ethical committee of the Aertztekammer Hamburg.
Statistical methods
To derive a dietary pattern predictive of CAD risk we applied RRR. RRR is a statistical dimension-reduction technique similar to principal component analysis; however, it works with 2 different sets of variables called predictors and responses (25). Whereas principal component analysis determines linear functions of predictors by maximizing the explained variation of all predictor variables, RRR identifies linear functions of predictors that explain as much response variation as possible. Principal component analysis is a special type of RRR in which the 2 sets of variables are identical. In the subsequent application of RRR, we choose food groups as predictors and biochemical markers for CAD as responses. Thus, we calculated linear functions of food group intakes (dietary patterns) that explain as much variation in CAD biomarkers as possible.
The implementation of RRR is based on the determination of eigenvalues and corresponding eigenvectors of the covariance matrix of responses. Choosing the eigenvectors as coefficient vectors for responses results in different linear functions of responses called response scores. These response scores have to be projected on the space of predictors to have linear functions of predictors called factor scores. The successive extracted scores are uncorrelated because the eigenvectors are always orthogonal. To each factor score belongs an eigenvalue quantifying the proportion of response variation that is explained by this factor only. Therefore, factors should be sorted by decreasing eigenvalues. The first factor obtained by RRR is optimal in the sense that no other linear function of predictors explains more response variation than does this factor. The method RRR is implemented in the special procedure PLS of SAS for WINDOWS (release 8.02; SAS/STAT Users Guide 1999). A more detailed description of the method, including the SAS code and its application in nutritional epidemiology, can be found elsewhere (26).
Adjusted means were calculated, and the differences in adjusted means were tested to be significantly different from zero by using the statement LSMEANS of the general linear model procedure in SAS. Conditional logistic regression was used for estimating relative risks of CAD adjusted for potential confounding variables. Trends in CAD risk across a continuous risk factor were assessed in logistic models by categorizing the risk factor in quintiles and by using the quintile number as an independent variable.
Biomarkers for CAD
For this RRR analysis we exclusively chose CAD biomarkers as response variables. We confined ourselves to primary biomarkers measured in serum, which can provide the link between the consumption of specific food groups and the development of CAD. Those biomarkers allow simultaneously the measurement of dietary exposure and act as indicators of health effects before the final outcome (27). Conventional biomarkers for CAD are LDL and HDL cholesterol. It has been shown by clinical intervention that lowering LDL-cholesterol concentrations and raising HDL-cholesterol concentrations decrease the risk of CAD (28). In several studies the predictive utility of serum cholesterol has been shown (2932). However, despite its predictive utility, LDL cholesterol may not necessarily distinguish controls from cases at the time of manifestation of CAD (33).
Beside cholesterol we chose lipoprotein(a) as a further response variable for RRR because it seems to be a predictor of future myocardial infarction, coronary death, and angina pectoris (3439). Several prospective studies showed that elevated lipoprotein(a) further increased the risk of CAD in subjects with low HDL-cholesterol and high LDL-cholesterol concentrations (38, 39). Finally, we chose C-reactive protein (CRP)the classic acute-phase proteinand C-peptide as the fourth and fifth response variables in our statistical analysis. Recent studies of CRP showed a consistent positive association with future cardiovascular events in initially healthy subjects (4043). High serum concentrations of C-peptide, which indicate insulin resistance, significantly increased the risk of CAD in diabetic and nondiabetic patients (4446).
| RESULTS |
|---|
|
|
|---|
The major contributors to the dietary pattern are quantified in Table 1
. Intakes of meat, margarine, other vegetable fats and oils (except margarine and olive oil), poultry, and sauce were directly associated with the score perceivable by positive score parameters as well as positive correlations with the score. In contrast, consumption of vegetarian dishes, wine, cooked vegetables, raw vegetables, and whole-grain cereals and muesli decreased the pattern score. The 10 food groups listed in Table 1
explained together 71.9% of the variation in score; the most important contributors were vegetarian dishes, wine, and meat.
|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
What is the interpretation of a dietary pattern score derived by RRR? In contrast with a pattern obtained by principal component analysis or factor analysis, an RRR factor score does not necessarily represent a combination of foods and drinks that are often eaten together in the study population. An RRR dietary pattern is a combination of food intakes that affects biomarker concentrations if biomarkers were chosen as response variables. The aim of this RRR application was not to describe real-word dietary behaviors but rather to find out what variation in diet is important for the development of disease. Persons with a high RRR score have a higher disease risk than do those with a low score. Moreover, the score coefficient of a food group allows to evaluate the single effect on disease of this food group supposing that the intakes of all other food groups are unchanged.
Participants who ate more meat, margarine, other vegetable fats and oils (except margarine and olive oils), poultry, and sauce, and consumed less vegetarian dishes, wine, vegetables, and whole-grain cereals had a high pattern score associated with a high risk of CAD. The observed effects of the food groups are broadly consistent with findings of others. High consumption of red and processed meat and high-fat products have been identified as a main constituent of Western diet that is associated with an unfavorable biomarker profile (48) and elevated risk of CAD (14, 15). In the Seven Countries Study, animal foods (including meat) were highly correlated with CAD death rates (49). The role of margarine in the development of CAD has been intensively studied (50, 51). Margarine, especially hard margarine, is a major source of trans fatty acids, which increase LDL-cholesterol and decrease HDL-cholesterol concentrations, and in turn, increase the risk of heart disease (52, 53). Although the content of trans fatty acids in margarine has strongly decreased in Germany, some German sunflower margarines still contain these unfavorable fatty acids (54).
Consistent with our results, several studies have shown a significant inverse relation between vegetable intake and risk of CAD (4, 49, 55). In a pooled analysis of the Nurses Health Study and the Health Professionals Follow-up Study, the relative risk of CAD in a comparison of extreme quintiles of vegetable intake was 0.82 (95% CI: 0.71, 0.94) (6). However, although fruit consumption also has a protective effect on CAD incidence it does not contribute to the RRR dietary pattern. The beneficial effect of whole grain intake as a constituent of the CAD dietary pattern obtained in the CORA Study was previously found in 2 large prospective studies (2, 56, 57).
Many studies have shown that light-to-moderate alcohol consumption was associated with a lower CAD risk (5862). An experimental dose of 30 g ethanol/d increased the concentrations of HDL cholesterol by
4 mg/dL, which indicated that the effect of alcohol consumption on CAD risk may partially be mediated by HDL cholesterol (63). There is controversy concerning the type of alcoholic beverage responsible for the effect (62). Our finding suggests that wine has an effect on CAD incidence, which is mediated through the selected biomarkers.
A main limitation of the present study was its case-control design. Although we did not include inpatients who received dietary advice, we cannot rule out that some cases changed their dietary habits as a consequence that they perceived experiencing the first symptoms of CAD. Moreover, knowledge of disease status can cause recall bias. Furthermore, selection bias in the recruitment of controls from the general population was possible because only 67% of the control candidates participated in this study, and the participants were probably more health conscious than were the nonresponders. It is necessary to replicate this RRR analysis in large cohort studies. The preliminary results of applying RRR in a nested case-control study (26) suggest that the association between dietary pattern score and disease could be weaker in a study of prospective design, especially if the number of incident cases is small. A possible explanation for this suggestion is that the variation in RRR score will be low in studies with many controls, provided that the selected responses do not differ markedly between the controls. The stronger results found with this RRR application than with those found with the previous approach (26) can partly be attributed to the use of biomarkers instead of disease-related nutrients as response variables.
Choosing biomarkers as response variables in RRR, the optimal time to take blood samples for the measurement of biomarker concentrations is a critical issue. In a case-control study design, blood samples are always taken from the cases after the diagnosed onset of the disease. Thus, biomarker concentrations can be affected by disease and thus should not be considered indicators of health effects preceding the disease. In contrast, a cohort study would ensure that biomarker measurements always be made before the final outcome. However, the time period that elapses between the time of the collection of blood samples and the time of the onset of disease can be large and individually different, which decisively diminishes the predictive power of biomarkers. The best time point to predict disease by biomarkers is immediately before disease. Assuming that biomarker concentrations do not rapidly decrease or increase with the onset of disease, which should be satisfied for the selected biomarkers (except for CRP), blood sampling within 24 h of the onset of symptoms as accomplished in the CORA Study should be next to optimal.
In conclusion, the new statistical method RRR is a powerful tool to derive dietary patterns that predict the incidence of CAD. This study showed that a diet high in vegetarian dishes, wine, vegetables, and whole-grain cereals and low in meat, margarine, poultry, and sauce is associated with a more favorable biomarker profile and a reduced risk of CAD. Together with other available evidence, these results can be used when developing dietary recommendations to prevent premature cardiovascular events. The CAD-specific dietary pattern score may also be helpful in evaluating individual eating habits in dietary counsel.
| ACKNOWLEDGMENTS |
|---|
B-CZ and EW planned and conducted the CORA Study. KH conducted the statistical analysis. KH and HB wrote the first draft of the manuscript. All authors stated having no conflict of interest due to authorship of the submitted paper.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
F. Imamura, A. H Lichtenstein, G. E Dallal, J. B Meigs, and P. F Jacques Generalizability of dietary patterns associated with incidence of type 2 diabetes mellitus Am. J. Clinical Nutrition, October 1, 2009; 90(4): 1075 - 1083. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Heroux, I. Janssen, M. Lam, D.-c. Lee, J. R Hebert, X. Sui, and S. N Blair Dietary patterns and the risk of mortality: impact of cardiorespiratory fitness Int. J. Epidemiol., April 20, 2009; (2009) dyp191v1. [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] |
||||
![]() |
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] |
||||
![]() |
A. Esmaillzadeh and L. Azadbakht Food Intake Patterns May Explain the High Prevalence of Cardiovascular Risk Factors among Iranian Women J. Nutr., August 1, 2008; 138(8): 1469 - 1475. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. McNaughton, G. D. Mishra, and E. J. Brunner Dietary Patterns, Insulin Resistance, and Incidence of Type 2 Diabetes in the Whitehall II Study Diabetes Care, July 1, 2008; 31(7): 1343 - 1348. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Nanri, D. Yoshida, T. Yamaji, T. Mizoue, R. Takayanagi, and S. Kono Dietary patterns and C-reactive protein in Japanese men and women Am. J. Clinical Nutrition, May 1, 2008; 87(5): 1488 - 1496. [Abstract] [Full Text] [PDF] |
||||
![]() |
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] |
||||
![]() |
D. R Jacobs Jr, L. F. Andersen, and R. Blomhoff Whole-grain consumption is associated with a reduced risk of noncardiovascular, noncancer death attributed to inflammatory diseases in the Iowa Women's Health Study Am. J. Clinical Nutrition, June 1, 2007; 85(6): 1606 - 1614. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A Nettleton, L. M Steffen, M. B Schulze, N. S Jenny, R G. Barr, A. G Bertoni, and D. R Jacobs Jr Associations between markers of subclinical atherosclerosis and dietary patterns derived by principal components analysis and reduced rank regression in the Multi-Ethnic Study of Atherosclerosis (MESA) Am. J. Clinical Nutrition, June 1, 2007; 85(6): 1615 - 1625. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Gao, P. E. Wilde, A. H. Lichtenstein, and K. L. Tucker The 2005 USDA Food Guide Pyramid Is Associated with More Adequate Nutrient Intakes within Energy Constraints than the 1992 Pyramid J. Nutr., May 1, 2006; 136(5): 1341 - 1346. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. B Schulze, K. Hoffmann, J. E Manson, W. C Willett, J. B Meigs, C. Weikert, C. Heidemann, G. A Colditz, and F. B Hu Dietary pattern, inflammation, and incidence of type 2 diabetes in women Am. J. Clinical Nutrition, September 1, 2005; 82(3): 675 - 684. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Weikert, K. Hoffmann, J. Dierkes, B.-C. Zyriax, K. Klipstein-Grobusch, M. B. Schulze, R. Jung, E. Windler, and H. Boeing A Homocysteine Metabolism-Related Dietary Pattern and the Risk of Coronary Heart Disease in Two Independent German Study Populations J. Nutr., August 1, 2005; 135(8): 1981 - 1988. [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] |
||||
![]() |
J. Dierkes, K. Hoffmann, K. Klipstein-Grobusch, C. Weikert, H. Boeing, B.-C. Zyriax, E. Windler, and J. Kratzsch Low plasma pyridoxal-5'phosphate and cardiovascular disease risk in women: results from the Coronary Risk Factors for Atherosclerosis in Women Study Am. J. Clinical Nutrition, March 1, 2005; 81(3): 725 - 727. [Full Text] [PDF] |
||||
![]() |
S. Friso, D. Girelli, N. Martinelli, O. Olivieri, and R. Corrocher Reply to J Dierkes et al Am. J. Clinical Nutrition, March 1, 2005; 81(3): 727 - 728. [Full Text] [PDF] |
||||
![]() |
K. Hoffmann, M. B. Schulze, and H. Boeing THE AUTHORS REPLY Am. J. Epidemiol., December 1, 2004; 160(11): 1132 - 1133. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |