Atherosclerosis is an arterial disease that induces lumen narrowing and obstruction, and is the most common cause of heart attacks and strokes worldwide. Several lifestyle and nutritional factors have been identified to influence the formation of atherosclerotic plaque(Reference Chan, Lau and Yiu1–Reference He, Liu and Daviglus4), but the relationship between Zn intake and atherosclerosis has been inconclusive, although Zn has been regarded as an antioxidant(Reference Powell5). The protective effects of Zn on the process of atherosclerosis have been suggested through in vitro and in vivo studies(Reference Wilkins and Leake6–Reference Truong-Tran, Carter and Ruffin9). Zn has been shown to inhibit LDL oxidation by macrophages and endothelial cells in vitro (Reference Wilkins and Leake6) and to prevent lesion development in hypercholesterolaemic animals(Reference Witztum and Steinberg7). Zn has also been shown to inhibit apoptosis of endothelial cells by inhibiting the activity of caspases(Reference Perry, Smyth and Stennicke8, Reference Truong-Tran, Carter and Ruffin9). Nevertheless, limited epidemiological studies have been done on the relationship between Zn intake and CVD, and the associations of Zn intake with CVD were not consistent(Reference Singh, Niaz and Rastogi10–Reference Iskra and Majewski15). A cross-sectional survey showed an inverse relationship between dietary Zn intake and the prevalence of coronary artery disease(Reference Singh, Niaz and Rastogi10), and a cohort study reported that dietary Zn intake showed an inverse association with CVD mortality(Reference Lee, Folsom and Jacobs11). Two case–control studies reported that CVD patients had higher Zn intakes and lower serum Zn than the controls(Reference Alissa, Bahjri and Ahmed12, Reference Ghayour-Mobarhan, Taylor and Kazemi-Bajestani13). However, other case–control studies demonstrated that serum Zn in individuals with atherosclerosis obliterans was higher than that in the healthy control groups(Reference Iskra, Patelski and Majewski14, Reference Iskra and Majewski15). Therefore, the effects of dietary Zn on atherosclerosis remain uncertain.
Intima-media thickness (IMT) of large artery walls is a non-invasive marker of early arterial wall alteration. Carotid IMT measurements assessed by B-mode ultrasound have been used as surrogate markers for atherosclerosis in observational studies and clinical research(Reference de Groot, Hovingh and Wiegman16). Increase in carotid IMT was directly associated with an increased risk of CHD including myocardial infarction and stroke in cohort studies(Reference O'Leary, Polak and Kronmal17, Reference Bots, Hoes and Koudstaal18). To date, no epidemiological studies have been reported on the effect of dietary Zn intake on morphological changes in the vascular wall. Thus, we examined the relationship between dietary Zn intake and subclinical atherosclerosis as measured using carotid IMT among the middle-aged and elderly populations.
Subjects and methods
Study population
The Korean Multi-Rural Communities Cohort Study in Rural Communities is part of the Korean Genome Epidemiology Study. The Korean Multi-Rural Communities Cohort Study in Rural Communities was initiated in 2004 to construct a genomic cohort and identify risk factors for CVD. As baseline, 7818 subjects aged ≥ 40 years were recruited from the centres located in Yangpyeong, Namwon and Goryeong between January 2005 and November 2007. Yangpyeong is located 45 km east of Seoul, the capital of South Korea, and Namwon and Goryeong are located in southwestern and southeastern areas of South Korea, respectively. The majority of the subjects were farmers and housewives. Among 7818 subjects, we excluded subjects who had hypertension (n 2131), hyperlipidaemia (n 348), diabetes (n 734), heart disease (n 502), cerebrovascular disease (n 273) or cancer (n 173) at the baseline survey, and those who reported implausible dietary intake ( < 2091 or >16 726 kJ/d or 1000 g of alcohol/d or more than ten missing food items) (n 67), and those who reported that they had changed their usual diet during the previous year (n 277). A total of 4564 subjects aged 40–89 years were included in the final data analysis. The present study was conducted according to the guidelines laid down in the Declaration of Helsinki, and all procedures involving human subjects were approved by the Institutional Review Boards of Hanyang University, Chonnam National University and Keimyung University. Written informed consent was obtained from all subjects.
General characteristics, anthropometrics and biochemical variables
To overcome the limitations of the multicentre studies, we developed standardised protocols for a questionnaire and each examination procedure including measurements of height and weight and blood pressure and blood sampling. All interviewers and examiners were trained by the same personnel at the coordinating centre.
A structured questionnaire, including information on demographics, education, smoking, alcohol intake, exercise, medical history and female reproductive history, was administered by trained interviewers. Height was measured using a standard height scale to the nearest 0·1 cm, and weight was measured using a metric weight scale to the nearest 0·01 kg in light clothing without shoes. BMI was calculated as weight (kg)/height (m2). Waist circumference was measured half way between the lowest rib margin and the iliac crest. Blood pressure was measured from the right arm by auscultation using a standard sphygmomanometer and a standard cuff. Two consecutive measurements of blood pressure were done after each subject had been sitting for at least 5 min. Systolic blood pressure and diastolic blood pressure were recorded to the nearest 2 mmHg. If two systolic or diastolic blood pressures were more than 5 mmHg apart, an additional measurement was done, and the mean value of the closest two measurements was used for the subsequent analyses. Blood samples were collected in the morning after at least 8 h of fasting, and all biochemical markers were analysed on the same day. Plasma total cholesterol, TAG, glucose and HDL-cholesterol were measured using the ADVIA1650 Automatic Analyser (Siemens, New York, NY, USA). If plasma TAG were < 4000 mg/l, LDL-cholesterol was calculated as described by Friedewald et al. (Reference Friedewald, Levy and Fredrickson19).
Dietary assessment
Dietary data were collected by trained interviewers using a food-based FFQ in which each participant was asked to identify the usual frequency of consumption of 106 food items during the previous year and the average serving size consumed. The FFQ consisted of nine frequency categories ranging from ‘never or rare’ to ‘3 times/d’, and three serving sizes were specified. For food items with different seasonal availability, the participants were asked to mark how long they had eaten among four categories: 3, 6, 9 and 12 months. The validity and reproducibility of the FFQ have been reported in detail elsewhere(Reference Ahn, Kwon and Shim20). Briefly, the estimated nutrient intakes from the FFQ were compared with those from 12 d food records during the four seasons. Pearson's correlation coefficient for Zn was 0·33. Based on the cross-classification analysis, with respect to Zn intake, 71·8 % of the subjects were classified into the same or adjacent quartiles, and 6·5 % of the subjects were classified into opposite quartiles. Nutrient intake was calculated using a weighted frequency per day and a portion size per unit in each food item. The seventh edition of the Food Composition Table of Korea was used as the nutrient database(21). The major contributing food sources of Zn were cereal grains (56 %) followed by oysters (3 %) and milk (3 %) (data not shown). Nevertheless, oyster ingestion explained 54 % of the total variation in dietary Zn intake, which is attributed to the high Zn content of oysters and large variation in oyster consumption among the subjects.
Measurement of intima-media thickness
IMT of common carotid artery was evaluated with the subject in the supine position using high-resolution B-mode ultrasound (SonoAce-9900; Medison Company Limited, Seoul, South Korea) equipped with a 7·5 MHz linear-array transducer. From the longitudinal view of the carotid bifurcation at a point 10 mm proximal to the common carotid artery, the maximal value of the IMT in a region free of plaque was measured in both carotid arteries. End-diastolic images were captured and saved as files for offline analysis by five sonographers at the three centres. A single trained reader at the reading centre in the Department of Preventive Medicine at Chonnam National University analysed the still images using the Sigma Scan Pro 5.0 (SPSS Inc., Chicago, IL, USA). Between- and within-sonographer reliabilities were evaluated using 180 cases and thirty cases, respectively, and within-reader reliability was tested using thirty-six cases. CV between and within sonographers (9·4–9·8 and 4·7–4·9 %, respectively), intraclass correlation coefficients between and within sonographers (r 0·85–0·88 and r 0·90–0·91, respectively), a CV within readers (3·8–4·1 %) and a intraclass correlation coefficients within readers (r 0·92–0·94) all indicated that the IMT measures taken from the three centres by five sonographers and one reader yielded highly reproducible values of common carotid IMT.
Statistical analysis
Nutrient intakes were adjusted for total energy intake by the residual method. Subjects were categorised into quintiles by dietary Zn intake. Stratification analysis was performed according to age of the subjects. The subjects aged < 65 years were considered ‘the middle-aged’, and the subjects aged ≥ 65 years were considered ‘the elderly’. The general linear model and the Cochran–Mantel–Haenszel analysis with adjustment for age and sex were used to determine differences in means or distribution and the linear trends across the Zn intake groups. The trend tests were conducted by treating the median value of each group as continuous variable in a multivariate model. The common carotid IMT was calculated as the mean IMT measurement of both carotid arteries. The differences in the means of IMT between Zn intake groups were assessed using general linear model. When statistically significant effects were demonstrated, Tukey's post hoc comparison test was used to identify group differences at P < 0·05. The normal values of carotid IMT have been established based on the distribution of IMT values in a general healthy population(Reference Howard, Sharrett and Heiss22, Reference Denarie, Gariepy and Chironi23), so that normal values are different depending on the methodology used for the measurement, age and sex. Subclinical atherosclerosis was defined by the value >80th upper percentile of common carotid IMT (0·8232 mm) similar to the study done by He et al. (Reference He, Liu and Daviglus4). However, since IMT tends to increase with age, this upper normal limit (0·8232 mm) can be too strict for elderly subjects. When this limit was applied to elderly subjects, 36 % of the elderly subjects aged ≥ 65 years were categorised as having subclinical atherosclerosis. Therefore, another criterion (IMT = 1 mm) was also applied to determine the risk of subclinical atherosclerosis. A value of IMT ≥ 1 mm at any age is related to the significant increase in myocardial infarction and/or cerebrovascular disease(Reference Simon, Gariepy and Chironi24). Non-conditional logistic regression analysis was applied to obtain the OR and corresponding 95 % CI for subclinical atherosclerosis. Potential confounders were identified by examining the distributions of characteristics that can affect atherosclerosis by Zn intake groups. The variables that showed significantly different distributions among Zn intake groups or showed significant linear trends across Zn intake groups were considered as potential confounders and adjusted in the analyses. Three different models were applied to examine the associations of Zn intake with IMT. Age and sex were adjusted as confounders in the first model, and education, alcohol intake (g/d), smoking status (current, former or non-smoker), regular exercise (yes/no), waist circumference and serum HDL-cholesterol were further adjusted in the second model. The final model additionally adjusted dietary factors such as quintiles of total energy, protein, fat, carbohydrate, β-carotene, vitamins C and E, folate, dietary fibre and cholesterol intakes. All statistical analyses were performed using the SAS software (version 9.1; SAS Institute, Inc., Cary, NC, USA), and P values < 0·05 were considered significant.
Results
The present study subjects consisted of Koreans aged 40–89 years sampled in three rural areas. General characteristics of the subjects are given in Table 1. Men made up 40·2 % of the total number of subjects. The mean ages of men and women were 61·3 and 59·7 years, respectively. The proportion of the subjects aged ≥ 65 years was 38·5 % of the total number of subjects. The proportions of the subjects whose BMI was >25 kg/m2 were 31·9 % of the men and 36·1 % of the women, respectively. The proportion of postmenopausal women was 78·4 %. Men were more likely to have a higher education, and to be current smokers and current drinkers than women. Means of total daily energy intakes were 7261·8 kJ (1736·6 kcal) for men and 6370·7 kJ (1523·5 kcal) for women. Means of daily Zn intake for men and women were 6·1 and 6·0 mg/d after adjustment for total energy intake by the residual method.
* All nutrients were total energy adjusted by residual method after log transformation.
To determine the potential confounding factors, the distributions of selected characteristics were examined across five Zn intake groups, which are given in Table 2. The median Zn intakes of the first, second, third, fourth and fifth quintiles were 5·00, 5·48, 5·85, 6·29 and 7·19 mg/d, respectively. Age, total energy intake, carbohydrate and proportion of current smokers decreased across the Zn intake groups. The proportions of higher education, current drinkers, regular exercise, waist circumference, serum HDL-cholesterol, and intakes of alcohol, protein, fat, β-carotene, vitamin C, folate, vitamin E, fibre and cholesterol increased across Zn intake groups. Differences in the proportion of males were significant. All variables for which significant linear trends or significant differences in means were used as covariates are given in Tables 3 and 4.
Q, quintile.
* P values for differences across groups and P values for the trends were determined using general linear model for continuous variables after adjustment for age and sex, and using the Cochran–Mantel–Haenszel test for categorical variables after adjustment for age (men and premenopausal women) or age and sex (higher education, current smoker, current drinker and regular exercise).
† All continuous variables were adjusted for age and sex.
‡ All nutrients were total energy adjusted by residual method after log transformation.
§ The proportion among women subjects.
∥ ≥ High school graduation.
Q, quintile.
a,b Mean values with unlike superscript letters within a row were significantly different among the three groups by Tukey's multiple comparison test.
* P values for differences across groups were obtained using the general linear model.
† Adjusted for age and sex.
‡ Adjusted for age, sex, education, alcohol intake, smoking, regular exercise, waist circumference and HDL-cholesterol.
§ Adjusted for age, sex, education, alcohol intake, smoking, regular exercise, waist circumference, HDL-cholesterol and dietary factors (quintiles of total energy, protein, fat, carbohydrate, β-carotene, vitamins C and E, folate, dietary fibre and cholesterol).
Q, quintile; IMT, intima-media thickness.
* The OR were adjusted for age and sex.
† The OR were adjusted for age, sex, education, alcohol intake, smoking, regular exercise, waist circumference and HDL-cholesterol.
‡ The OR were adjusted for age, sex, education, alcohol intake, smoking, regular exercise, waist circumference, HDL-cholesterol and dietary factors (quintiles of total energy, protein, fat, carbohydrate, β-carotene, vitamins C and E, folate, dietary fibre and cholesterol).
Adjusted IMT means by Zn intake group are given in Table 3. The middle-aged ( < 65 years) and the elderly ( ≥ 65 years) populations showed similar median Zn intakes. Means of IMT in the elderly population were higher than those in the middle-aged population. After adjustment for age, sex, education, alcohol intake, smoking, regular exercise, waist circumference, HDL-cholesterol and dietary factors (quintiles of total energy, protein, fat, carbohydrate, β-carotene, vitamins C and E, folate, dietary fibre and cholesterol intakes), adjusted IMT mean in the first quintile of Zn intake group was significantly higher than that in the fifth quintile group among the total subjects. In the middle-aged population, the first quintile group showed higher adjusted IMT means than the fifth quintile group in the first and second models. With respect to the elderly population, after adjusting for potential confounders, adjusted IMT mean in the first quintile group was higher than that in the fifth quintile group, but the difference was NS. The interaction between quintile of Zn intake and age group was not statistically significant (data not shown). Adjusted IMT means by Zn intake from seafood and grain showed the same direction with the relationship between total dietary Zn and IMT (data not shown).
The relationships between Zn intake and the risk of subclinical atherosclerosis are given in Table 4. When subclinical atherosclerosis was defined by >80th percentile value of IMT (0·8232 mm), 10·4 % of the middle-aged population and 35·7 % of the elderly population were classified as having subclinical atherosclerosis. After adjustment for potential confounders in the third model, Zn intake was inversely related to subclinical atherosclerosis (5th v. 1st quintile, OR 0·64, 95 % CI 0·45, 0·90, P for trend = 0·069). In the middle-aged population, a significant inverse correlation between Zn intake and subclinical atherosclerosis was found in the third model (5th v. 1st quintile, OR 0·44, 95 % CI 0·26, 0·76, P for trend = 0·048). However, no significant relation was found in the elderly population.
When subclinical atherosclerosis was defined by ≥ 1 mm of IMT, 1·6 % of the middle-aged population and 7·9 % of the elderly population were classified as having subclinical atherosclerosis. Zn intake was inversely correlated with the risk of subclinical atherosclerosis after adjusting for potential confounders (5th v. 1st quintile, OR 0·34, 95 % CI 0·16, 0·70, P for trend = 0·005). The elderly subjects showed an inverse relationship between Zn intake and the risk of subclinical atherosclerosis (5th v. 1st quintile, OR 0·33, 95 % CI 0·14, 0·77, P for trend = 0·013). No significant correlations were found in middle-aged adult subjects.
When the subjects who were diagnosed as having hypertension, hyperlipidaemia, diabetes, heart disease, cerebrovascular disease or cancer were included and diagnosed diseases were adjusted in the model, similar trends were found between Zn intake and the risk of carotid atherosclerosis (results not shown).
Discussion
Carotid IMT is known as a surrogate marker for carotid atherosclerosis. In the present cross-sectional study, the low Zn intake group displayed a higher mean IMT than the high Zn intake group. An inverse relationship between dietary Zn intake and common carotid atherosclerosis was demonstrated in the middle-aged and elderly populations. The present findings suggest that adequate Zn intake may have a beneficial effect on the risk of atherosclerosis among the middle-aged and elderly populations.
In the present study, daily Zn intake was estimated using the FFQ. However, it is inappropriate to compare the absolute dietary Zn intake estimated using the FFQ to the estimated average requirements or dietary Zn intakes estimated by other dietary assessment methods. Therefore, we conducted 3 d 24 h recalls in a subgroup of the subjects (n 880) (data not shown). The mean energy intakes for men and women in the subgroup were 7439 (sd 2273) kJ and 5890 (sd 1750) kJ, which were similar to the energy intakes estimated by the FFQ. The means of daily Zn intakes for men (n 350) and women (n 530) were 8·7 and 7·4 mg/d, which were higher than the Korean estimated average requirements for Zn(25), but lower than Zn intakes estimated by multiple day 24 h recalls from a healthy Chinese sample (8·9–11·9 mg/d)(Reference Ma, Li and Jin26), an older European sample (11·0–11·3 mg/d) and an adult sample from the United States (9·7–14·2 mg/d)(Reference Andriollo-Sanchez, Hininger-Favier and Meunier27–Reference Moshfegh, Goldman and Cleveland29).
Carotid IMT has been used as an intermediate or proxy end point for CVD. Carotid IMT has been consistently associated with the incidence of CHD in cohort studies(Reference O'Leary, Polak and Kronmal17, Reference Bots, Hoes and Koudstaal18, Reference Chambless, Heiss and Folsom30). The association between carotid IMT and diet, such as consumption of fruit, olive oil, fat, alcohol, fish and isoflavone, has been determined(Reference Chan, Lau and Yiu1–Reference He, Liu and Daviglus4, Reference Buil-Cosiales, Irimia and Berrade31), but the relationship between carotid IMT and Zn intake has not been investigated.
In the present study, the risk of subclinical atherosclerosis in the high Zn intake group was significantly lower than that in the low Zn intake group. In particular, when 80th upper percentile values of common carotid IMT were used to determine the OR, after adjusting for potential confounders, the risks of subclinical atherosclerosis in the fifth quintile group were lower than those in the first quintile group among the middle-aged populations (5th v. 1st quintile, OR 0·44, 95 % CI 0·26, 0·76, P for trend = 0·048). The elderly population did not show any significant associations, which can be attributed to the low cut-off value for subclinical atherosclerosis for the elderly subjects. Since IMT increases with age, the low cut-off for subclinical atherosclerosis was not appropriate for detecting the effect of dietary Zn among the elderly subjects. When IMT ≥ 1 was used to determine the OR, a significant inverse relationship between Zn intake and subclinical atherosclerosis was found only among the elderly subjects aged ≥ 65 years (5th v. 1st quintile, OR 0·33, 95 % CI 0·14, 0·77, P for trend = 0·013). Since the number of subjects with IMT ≥ 1 was small among the middle-aged adult subjects, any significant association was not found in the middle-aged adult subjects. Taken together, these results suggest that adequate dietary Zn intake may have a beneficial effect on the risk of atherosclerosis among the middle-aged adult and elderly populations, and different criteria for subclinical atherosclerosis may be needed between age groups.
Limited epidemiological studies have examined the relationship between Zn intake and atherosclerosis. Dietary Zn intake was positively correlated with serum Zn levels in cross-sectional studies(Reference Singh, Niaz and Rastogi10, Reference Ghayour-Mobarhan, Taylor and New32), and lower dietary Zn intake and lower concentrations of serum Zn were associated with a higher prevalence of CHD(Reference Singh, Niaz and Rastogi10, Reference Alissa, Bahjri and Ahmed12) or an increased cardiovascular mortality(Reference Reunanen, Knekt and Marniemi33). Lower serum Zn concentrations have been reported in patients with various forms of coronary artery diseases in case–control studies(Reference Ghayour-Mobarhan, Taylor and Kazemi-Bajestani13, Reference Kazemi-Bajestani, Ghayour-Mobarhan and Ebrahimi34–Reference Kerkeni, Added and Ben Farhat36).
Several potential mechanisms may account for the cardiovascular benefits of Zn(Reference Wilkins and Leake6, Reference Perry, Smyth and Stennicke8, Reference Truong-Tran, Carter and Ruffin9, Reference Meerarani, Ramadass and Toborek37–Reference Beattie and Kwun43). Atherogenesis involves the oxidation of LDL and its uptake into macrophages, and apoptosis of endothelial cells. Zn ions inhibited LDL oxidation by macrophages and endothelial cells(Reference Wilkins and Leake6). The activity of caspase enzymes involved in apoptotic pathways was inhibited by Zn(Reference Perry, Smyth and Stennicke8, Reference Truong-Tran, Carter and Ruffin9). Zn deficiency increased apoptotic cell death induced by linoleic acid and TNF-α, and supplementation of Zn attenuated apoptosis(Reference Meerarani, Ramadass and Toborek37).
NO is an endogenous signalling molecule produced by NO synthase. NO plays a substantial role in platelet and leucocyte function in human subjects. Several studies have suggested that NO can inhibit leucocyte adhesion through various mechanisms(Reference Freedman and Loscalzo38–Reference Kubes40). Zn plays a key role in endothelial NO synthase function by stabilising a dimer of two subunits of NO synthase(Reference Zou, Shi and Cohen41). Zn deficiency can reduce the activity of NO synthase, which causes a change in NO production. Changes in NO production are associated with the development of atherosclerosis(Reference Anderson42). Indeed, NO is involved in Zn release from proteins by thiol nitrosylation(Reference Beattie and Kwun43).
However, in the case–control studies conducted by Iskra et al. (Reference Iskra, Patelski and Majewski14, Reference Iskra and Majewski15), serum concentrations of Zn in individuals with atherosclerosis obliterans were higher than those in the healthy control groups. Even if the present findings reflect a causal relationship between Zn in food and atherosclerosis, this would not imply that the same amount or more Zn in supplements would have the same effect. High dose of Zn supplements ( ≥ 50 mg) induced a decline in plasma HDL-cholesterol concentrations in healthy young males(Reference Chandra44–Reference Black, Medeiros and Brunett46), which demonstrates that a high dose of Zn supplementation may induce an atherogenic lipoprotein profile. The tolerable upper intake levels of Zn for Koreans and Americans are 35(25) and 40 mg/d(47), respectively. Therefore, the role of Zn in CHD needs to be investigated further.
There are several limitations to the present study which need to be considered when interpreting the findings. First, convenience sampling may preclude generalisation of the results in the present study. Secondly, since this was a cross-sectional study, we cannot conclude causality of Zn intake for subclinical atherosclerosis. Thirdly, there were no subjects who took Zn supplement on a regular basis, and a total of 269 subjects (5·9 %) answered that ‘they are taking multivitamin supplements’, but most of the subjects could not remember the brand names of the multivitamins. Therefore, Zn intake from supplement was not considered. Fourthly, during the assessment of dietary Zn status, the bioavailability of Zn should be taken into account(48). Phytate is a major inhibitory factor for Zn absorption, and phytate is present in whole grains, cereals and legumes. The phytate:Zn molar ratio has been used in several studies to predict the inhibitory effect of phytate on the bioavailability of Zn(Reference Ma, Li and Jin26, 48, Reference Turnlund, King and Keyes49). However, due to a lack of information on the phytate content of foods in the Korean Food Composition Table, the inhibitory effect of phytate on the bioavailability of Zn was not taken into account in the present study. Finally, despite the lack of a reliable and specific index, Zn status of the subjects needs to be assessed by biomarkers such as plasma and/or erythrocyte Zn concentrations as well as dietary intake, which was not done in the present study(Reference Neggers, Goldenberg and Tamura50).
To date, this is the first epidemiological study to determine the relationship between dietary Zn intake and IMT, an indicator of subclinical atherosclerosis. Among the rural middle-aged adult and elderly populations in Korea, inadequate intake of Zn was related to high IMT and high risk of subclinical atherosclerosis. The present findings suggest that adequate intake of Zn may have a beneficial effect in reducing the risk of atherosclerosis. Further studies are necessary to examine these findings in prospective studies.
Acknowledgements
The present study was supported by the Korea Centers for Disease Control and Prevention (grants 2004-347-6111-213, 2005-347-2400-2440-215, 2006-347-2400-2440-215 and 20070308455-00). The authors' responsibilities are as follows: B. Y. C., B.-Y. C., S.-S. K., Y.-H. L. and M. K. K. designed the overall study and oversaw the data collection; Y. J. Y. and M. K. K. conducted the statistical analysis and wrote the manuscript draft; P. S. P. and M. K. K. contributed to dietary data collection and B. Y. C., B.-Y. C., S.-S. K., Y.-H. L., P. S. P. and M. K. K. provided critical reviews. None of the authors has any personal or financial conflicts of interest.