CVD is the leading cause of death worldwide, although it has been posited that up to 90 % of cases may be preventable simply by modifying diet and lifestyle( Reference Yusuf, Reddy and Ounpuu 1 – Reference Perk, De Backer and Gohlke 3 ). While CVD contributes less, proportionally, to overall mortality in the developing world (~23 % compared with ~49 % in more economically developed countries), these areas are home to the vast majority of the global population as well as over 60 % of worldwide CVD mortality( Reference Reddy and Yusuf 4 ). Additionally, further industrialization stands to yield a decrease in infectious disease mortality and concurrent increase in the burden of CVD( Reference Yusuf, Reddy and Ounpuu 1 , Reference Omran 5 , Reference Islam 6 ). Hence, as CVD is typically influenced by modifiable lifestyle factors or exposures, it becomes advantageous to recognize risk factors and preclinical markers of disease in populations with varying socio-economic characteristics. This not only facilitates more effective preventive strategies, but also allows for earlier interventions.
Carotid artery intima-media thickness (IMT) is a validated surrogate marker of preclinical atherosclerosis that has been shown to predict cardiovascular morbidity and mortality( Reference Stein, Korcarz and Hurst 7 ). Many studies have sought to examine CVD risk factors in relation to carotid IMT as a mechanistic intermediate. Particularly useful would be a better understanding of how diet, a critical modifiable risk factor for CVD, affects carotid IMT. Previous observational studies have shown high intakes of plant foods and fibre as well as low intakes of saturated fats and processed foods to be inversely associated with carotid IMT( Reference Nettleton, Schulze and Jiang 8 – Reference Yang, Zhang and Sun 11 ), while additional experimental studies have tested the effect of dietary or lifestyle interventions on carotid IMT( Reference Wildman, Schott and Brockwell 12 – Reference Sala-Vila, Romero-Mamani and Gilabert 15 ). More recently, it has been suggested that adherence to a ‘Mediterranean’ diet may protect against a higher burden of carotid atherosclerotic plaque( Reference Gardener, Wright and Cabral 16 ). Data are scarce, however, with respect to comparing whole dietary patterns within a given population( Reference Mikkila, Rasanen and Laaksonen 17 ). Additionally, the vast majority of the aforementioned studies were conducted in Western populations, and little is known about the role of diet in IMT in low-income and low-BMI countries.
In addition to analysing intakes of specific foods or nutrients, evaluation of diet patterns may yield a more comprehensive understanding of how this complex factor affects the aetiology of disease( Reference Hu 18 ). Compared with the ‘single nutrient’ approach, this method of analysis more fully accounts for interactions between nutrients, as well as interrelationships between dietary components that cause difficulty in distinguishing individual food or nutrient effects( Reference Willett 19 ).
In the present study, we examined the association of major food groups, as well as dietary patterns, with carotid IMT in a rural Bangladeshi population.
Participants and methods
Study population
The parent study, Health Effects of Arsenic Longitudinal Study (HEALS), is an ongoing population-based prospective cohort study in Araihazar, Bangladesh( Reference Ahsan, Chen and Parvez 20 ). Briefly, between October 2000 and May 2002, 11 746 men and women (‘original cohort’) were recruited from a well-defined 25 km2 geographical area, under the criteria that all were married (to reduce loss to follow-up), between 18 and 75 years old and had resided in the study area for at least 5 years. From 2006 to 2008, HEALS was expanded to include an additional 8287 participants (‘expansion cohort’) following the same methodologies. The overall response rate was 97 %. Study participants underwent baseline clinical assessment and structured interviews. Detailed explanations of data collection for such variables as height, weight, smoking status and blood pressure have been published elsewhere( Reference Pierce, Kalra and Argos 21 – Reference Wu, Chen and Parvez 23 ). Informed consent was obtained from study participants; study procedures were approved by the Ethical Committee of the Bangladesh Medical Research Council and the Institutional Review Boards of Columbia University and the University of Chicago.
Assessment of food intake
Trained interviewers administered a thirty-nine-item FFQ during the baseline assessment. Details of this questionnaire and its validation are available elsewhere( Reference Chen, Ahsan and Parvez 24 ). Briefly, HEALS investigators worked with local nutrition experts to identify all food items available at the village market in the study area. This formed the basis for a preliminary version of the FFQ to be developed through a series of subsequent extensive discussions with ten focus groups. After pilot testing among 120 local individuals who were not part of the cohort study, trained interviewers completed the FFQ with HEALS participants via in-person interviews. For each food item on the FFQ, participants were queried regarding the amount eaten per meal, as well as the frequency of consumption per day, month and year.
Analysis of 189 respondents randomly selected to complete 7 d food diaries showed no major differences when compared with data previously obtained in the FFQ (correlation 0·30 to 0·76), including common food items, macronutrients and many major micronutrients. This food diary sampling was the largest number possible for the present study, given the financial and logistical constraints of this more intensive form of dietary monitoring. The US Department of Agriculture’s Nutrient Database for Standard Reference (abbreviated version)( 25 ) and an Indian food nutrient database( Reference Gopalan, Rama Sastri and Balasubramanian 26 ) were used to convert food intakes to nutrient intake values. Although there is no ideal measure of dietary intakes, food diaries are known to be superior to FFQ and have previously been frequently used to assess the validity of FFQ( Reference Willett, Sampson and Stampfer 27 , Reference Erkkola, Karppinen and Javanainen 28 ).
Dietary patterns
From our raw FFQ data, intakes of major food groups were calculated by summing the average daily intake of the relevant individual food items, as listed in Appendix Table 1. Additionally, we have previously described our delineation of three major dietary patterns among participants of in the original cohort( Reference Chen, Factor-Litvak and Howe 22 ). Briefly, we performed principal component analysis using the PROC FACTOR procedure in the SAS statistical software package to identify dietary patterns, with an orthogonal rotation (the varimax option in SAS) and conventional criteria including eigenvalue, the scree test, proportion of variance accounted for and the interpretation criterion( Reference Hatcher 29 ). To avoid any arbitrary decisions on food group definitions, food items were directly entered into the principal component analysis as average daily intakes in grams or millilitres. Each food item received a factor loading associated with each diet pattern, and the factor loading represents the correlation coefficient between the food item and the diet pattern. To indicate a participant’s relative standing on each diet pattern in the population, a factor score, which is a linear composite of the optimally weighted food items by factor loadings, was constructed for each dietary pattern. Each participant received a factor score for each identified diet pattern. The identification of dietary patterns was conducted separately for the original and expansion cohort, as expansion cohort members were recruited later. The dietary patterns were mostly similar between the two populations (Appendix Table 2). Three major dietary patterns were identified. The ‘balanced’ pattern was comprised of steamed rice, red meat, fish, fruit and vegetables. An ‘animal protein’ diet was more heavily weighted towards eggs, milk, red meat, poultry, bread and vegetables. The ‘gourd and root vegetable’ diet relied heavily on a variety of gourds, radishes, pumpkin, sweet potato and spinach.
We have previously evaluated the associations of these diet patterns in relation to blood pressure at baseline( Reference Chen, Factor-Litvak and Howe 22 ), as well as risk of skin lesions( Reference Pierce, Argos and Chen 30 ) and cardiovascular mortality( Reference Chen, McClintock and Segers 31 ).
Assessment of carotid intima-media thickness
Carotid IMT was measured between April 2010 and September 2011. Eight hundred participants were randomly selected for IMT measurement, from the 11 224 original cohort members who provided urine samples at baseline, as part of a previous study on urinary arsenic and IMT( Reference Chen, Hakim and Parvez 32 ). An additional 700 participants were randomly sampled for IMT measurement from the 5136 participants older than 30 years of age in the expansion cohort, as the expansion cohort was, on average, younger. In total, carotid IMT was measured for 1151 individuals, consisting of 575 from the original cohort and 576 from the expansion cohort. A total of 349 participants did not complete IMT measurements due to deaths, move, serious illness or time constraints. After two cases were removed due to incomplete baseline data, the final study population consisted of 1149 participants.
Measurements of carotid IMT were conducted using a SonoSite MicroMaxx ultrasound machine (SonoSite, Inc., Bothell, WA, USA) equipped with an L38e/10–5 MHz transducer. All carotid imaging and IMT measurements were performed by a single physician who was trained and certified to perform carotid ultrasound measurements according to the specific ultrasound imaging and reading protocols developed, implemented and validated in the Oral Infections and Vascular Disease Epidemiology Study (INVEST)( Reference Desvarieux, Demmer and Rundek 33 ). The carotid arteries were scanned longitudinally in three carotid segments, using the optimal angle and the lateral extent of each carotid segment as defined relative to the tip of the flow divider, which is normally the most clearly defined anatomic reference in the proximity of the carotid bifurcation. IMT measurements were analysed offline with Matlab (Mathworks, Natick, MA, USA), which automatically calculated the distances between boundaries and expressed the results as the mean and maximal value. In accordance with the Manheim consensus, we used the mean of the near and far walls of the maximum common carotid artery IMT from both sides of the neck (mean of four measurements) as the main outcome variable, similar to previous studies( Reference O’Leary, Polak and Kronmal 34 – Reference Touboul, Hennerici and Meairs 36 ).
Assessment of covariates
All covariates were derived from the HEALS baseline interview. Sociodemographic factors utilized in the current analysis included age, sex, formal education (years), land ownership (yes/no), television ownership (yes/no), smoking (never, former or current), diabetes history and religious affiliation (Muslim, Hindu, Christian, other). With respect to tobacco, details of cigarettes and bidis (filterless locally produced cigarettes) were asked together, including information on past or current use, duration of use, age at start, number of sticks per day and age at quitting. Trained study physicians measured height and weight using a locally manufactured tape measure and a Misaki (Okaka, Japan) scale (calibrated weekly), respectively. Both height and weight were measured three times at baseline and averaged( Reference Pierce, Kalra and Argos 21 ). Trained clinicians measured blood pressure with an automatic sphygmomanometer( Reference Chen, Factor-Litvak and Howe 22 ). Measurements were taken with participants in a seated position after 5 min of rest, with the cuff around the upper left arm, in accordance with recommended guidelines. BMI was calculated as [average weight (kg)]/[average height (m)]2. While other data were collected at baseline, we selected the aforementioned variables in particular as they were deemed most relevant to cardiovascular risk factor analysis. Data on physical activity, lipids and female menopause status were not available.
Statistical analyses
Descriptive analyses were conducted first to assess relationships of population characteristics and dietary intake with carotid IMT. Significance testing was conducted using the χ 2 test for categorical variables and ANOVA for continuous variables. Continuous variables included carotid IMT, age, BMI, education, blood pressure and energy intake. Sex, cigarette use, diabetes, food group intake and dietary pattern adherence were assessed as categorical variables. The primary analysis evaluated the relationship between IMT and major food groups as well as dietary patterns. For this, linear regression was used, with dietary elements and patterns in quartiles as independent variables and common carotid IMT as the dependent variable. Quartiles of food groups were formed with respect to the combined cohort, as intake was recorded on an absolute scale. Due to relatively low levels of red meat and poultry consumption, all those with zero intake were grouped into the reference (first) quartile and individuals with non-zero intakes were then divided into the remaining higher quartiles. Because dietary patterns were derived separately in the original and expansion cohorts, cohort-specific quartiles for dietary patterns were used in the analyses. In the first model, we adjusted for age (years), sex, daily energy intake and smoking status (ever/never). Additional adjustment for educational attainment did not change results appreciably (data not shown). Subsequent models adjusted for BMI, blood pressure and diabetes mellitus (models 2 and 3). Physical activity metrics were not available for integration into these final models. IMT differences associated with a 1 sd difference in food groups and dietary pattern were computed for all models. Sensitivity analyses were conducted after excluding thirty-five participants with history (prior to date of IMT) of CHD or stroke to minimize the potential for reverse causality (i.e. dietary changes resulted from life-changing atherosclerotic cardiovascular events). All analysis was conducted using the statistical software package IBM SPSS Statistics version 20.
Results
Table 1 shows the distribution of dietary intakes and study population characteristics by carotid IMT quartile. Participants with a greater IMT measurement were more likely to be older (P<0·01). A positive correlation was seen between BMI and IMT (P=0·05). It is noteworthy that this population displayed a relatively low BMI, with a group mean of 19·95 kg/m2 at baseline; diabetes differed significantly by quartile (P<0·01), with baseline prevalence in the highest quartile equal to 4·1 % (compared with 1·0 % in the lowest quartile). The only food group or dietary pattern to achieve significance in descriptive analysis was the animal protein diet, as the percentage above the median intake of that pattern increased steadily from the first to fourth quartile (42·9 % to 61·5 %, respectively).
IMT, intima-media thickness.
* Blood pressure measurements were not valid for two study participants.
† Blood pressure measurements were not valid for one study participant.
‡ The χ 2 test was used for categorical variables and ANOVA for continuous variables.
The animal protein diet was positively associated with a variety of markers of socio-economic status, as shown in Table 2. These included educational attainment (P<0·01), television ownership (P<0·01), land ownership (P<0·01) and BMI (P<0·01). The gourd/root vegetable pattern was inversely related to both smoking (P<0·01) and male sex (P<0·01).
Q, quartile; IMT, intima-media thickness.
* Trend tests were performed by using linear regression that included energy-adjusted factor scores as the dependent variable and age, length of education, BMI, sex, smoking status, marital status, religion, television ownership or land ownership as the independent variable.
There was a significant difference (P<0·01) in carotid IMT between men (mean IMT 822·06 μm) and women (769·36 μm; Appendix Table 3). As expected, blood pressure showed a strong positive correlation with IMT (P<0·01). Although the population was relatively lean, BMI positively correlated with IMT as well (P<0·01).
The balanced diet was associated with lower IMT in the fully adjusted model, with each 1 sd increase in balanced diet adherence associated with a 4·95 (95 % CI −9·78, −0·11) μm decrease in IMT (P=0·045; Table 3). In model 1, each 1 sd increase in adherence to the animal protein diet was associated with an increase of 7·70 (95 % CI 2·45, 12·95) μm in IMT. This association, however, did not persist in the fully adjusted model (P=0·45). In all models, the gourd/root vegetable diet was strongly associated with increased IMT. Compared with participants in the bottom quartile of the gourd/root vegetable diet, IMT in participants in the second, third and fourth quartile was higher by 4·41 (95 % CI −8·95, 17·77) μm, 14·77 (95 % CI 1·33, 28·22) μm and 17·66 (95 % CI 3·97, 31·35) μm, respectively (model 3). IMT increase per 1 sd equalled 7·74 (95 % CI 2·86, 12·62) μm in the fully adjusted model (P<0·01). Sensitivity analyses excluding those with a history of CHD or stroke yielded very similar results (data not shown). It is noteworthy that, in contrast to the relationships observed in dietary pattern analysis, our data did not suggest that intakes of individual food groups were significantly related to IMT in this study population.
IMT, intima-media thickness; Q, quartile.
* P values were estimated from linear regression with dietary variables as continuous variables.
† Model 1 adjusted for age (years), sex, daily energy intake and smoking status (ever/never).
‡ Model 2 adjusted for age (years), sex, daily energy intake, smoking status (ever/never) and BMI.
§ Model 3 adjusted for age (years), sex, daily energy intake, smoking status (ever/never), BMI, systolic blood pressure and current diabetes diagnosis (yes/no).
Shown in Table 4 are the Spearman correlation coefficients between nutrient intakes and factor scores for dietary patterns, which signify adherence to each respective diet. Among all of the patterns, the balanced diet was more positively correlated with K, fibre, folate and vitamin B6; conversely, it was least correlated with polyunsaturated fat. The animal protein diet was highest in association with protein, Na, total fat, vitamin B12, polyunsaturated fat, saturated fat, monounsaturated fat and riboflavin, while correlating least with carbohydrate, fibre, folate and Mg. The gourd/root vegetable pattern was related to higher intake of carbohydrate and Mg; of the diets, it was least associated with protein, K, Na, total fat, vitamin B6, vitamin B12, saturated fat, monounsaturated fat and riboflavin.
Discussion
Our data support a positive association of a gourd/root vegetable diet and IMT level (7·74 μm increase in model 3; P<0·01). Additionally, adherence to the balanced diet was associated with lower IMT (4·95 μm decrease in model 3; P=0·045). Although single food groups were also measured in the study (meat, poultry, fish, fruit, vegetables), none displayed a statistically significant relationship with IMT. Given the high prevalence of CVD and considering that prior studies have shown clinically relevant differences in morbidity and mortality corresponding to IMT increases of 8 to 40 μm (albeit in studies of varying design types)( Reference Zureik, Touboul and Bonithon-Kopp 37 – Reference Crouse, Byington and Bond 40 ), the associations observed here could translate to substantial alterations in disease burden.
Many studies have provided evidence of certain foods or nutrients either enhancing or mitigating CVD risk( Reference Mozaffarian, Appel and Van Horn 41 ); however, individuals do not consume these dietary components in isolation. Dietary pattern analysis, although reliant on unique characteristics of food consumption in a given population that may limit generalizability of subsequent findings, can yield valuable information that may have otherwise been clouded in studies of single foods or nutrients. Indeed this was the case in the present study, as individual food groups did not alone convey any significant impact on IMT. Previous research on diet patterns and IMT, however, is limited, with one prior study finding no association for any diets identified in a middle-aged French population( Reference Kesse-Guyot, Vergnaud and Fezeu 42 ) and another finding significance only among men (but not women) in comparing ‘traditional’ and ‘health-conscious’ diets among a population of Finns below the age of 40 years( Reference Mikkila, Rasanen and Laaksonen 17 ).
Atherosclerosis is an inflammatory phenomenon that is an important pathological process in the development of ischaemic CVD( Reference Ross 43 ). It has been hypothesized that certain diet choices can lead to a chronic inflammatory state, predisposing to atherosclerosis and increasing subsequent risk of cardiovascular morbidity and mortality( Reference De Caterina, Zampolli and Del Turco 44 ). Our findings regarding the balanced diet in this population are largely consistent with existing nutrition literature on atherosclerotic development, as ‘cardioprotective’ diets have been described as being high in vegetables, nuts and monounsaturated fats( Reference De Caterina, Zampolli and Del Turco 44 – Reference Millen, Quatromoni and Nam 46 ), while ‘harmful’ factors include trans-fatty acids, foods with a high glycaemic index, processed meats, refined grains and high-fat foods( Reference Mente, de Koning and Shannon 45 , Reference Nettleton, Steffen and Schulze 47 ). This is illustrated in Table 4 where the balanced pattern corresponds with a lower intake of both saturated fat and total fat, as well as in Appendix Table 1, where the balanced pattern is shown to de-emphasize animal protein and integrate a wider variety of vegetables when compared with the animal protein pattern.
The findings presented here are consistent with some aspects our previously published data on risk of cardiovascular mortality associated with these same three dietary patterns in the parent HEALS study population( Reference Chen, McClintock and Segers 31 ). For example, we found borderline significantly reduced risk of CVD and heart disease associated with the balanced diet pattern (hazard ratio=0·86; 95 % CI 0·73, 1·01) that was associated with decreased IMT in the present study (P=0·045). Likewise, although the observed positive association between the animal protein diet and IMT was not statistically significant, the data were consistent with our previous analyses that reported an increased risk of mortality found for the animal protein diet (hazard ratio=1·13; 95 % CI 1·00, 1·28). The gourd/root vegetable diet, although associated with higher IMT here (P<0·01), was not related to risk of mortality in our prior analysis (hazard ratio=1·03; 95 % CI 0·89, 1·18). It is possible that the lack of agreement between findings of the present study and some of the mortality results may be due to the fact that we focus on a single preclinical indicator of atherosclerosis and, therefore, may not capture the full complement of multifactorial, complex pathological processes that underlie the association of diet with cardiovascular morbidity and mortality.
The observation that a gourd and root vegetable diet in our population is related to higher levels of IMT contradicts what has been found in some recent studies of food patterns( Reference Yang, Li and Zhang 10 , Reference Yang, Zhang and Sun 11 , Reference Key, Fraser and Thorogood 48 ). For example, a recent series of studies comparing vegetarian and omnivore diets in Chinese men associated a vegetarian diet with decreased IMT( Reference Yang, Li and Zhang 10 , Reference Yang, Zhang and Sun 11 ). It is crucial to note, however, that specific nutritional deficiencies may be present in certain types of vegetable-based diets and can predispose adherents to adverse health outcomes, especially in a lean population( Reference Lee, Woo and Chen 49 , Reference Davis and Kris-Etherton 50 ). In particular, such diets may be lacking in certain n-3 fatty acids (namely EPA and DHA)( Reference Davis and Kris-Etherton 50 ) and vitamin B12 Reference Ingenbleek and McCully (51 , Reference Siri, Verhoef and Kok 52 ), which is obtained solely through consumption of animal products( Reference Green and Kinsella 53 ). As a result, the nutrient deficiencies in such diets may negate the coinciding decreased fat intake that otherwise would have conveyed cardiovascular benefits. B12 insufficiency has been shown to be most widespread in developing countries( Reference Ingenbleek and McCully 51 , Reference Yajnik, Deshpande and Lubree 54 ), with South Asian populations considered to be at uniquely high risk( Reference Gammon, von Hurst and Coad 55 ). Consistently, we have reported a high prevalence of hyperhomocysteinaemia in our study population( Reference Gamble, Ahsan and Liu 56 ). It should lastly be considered that the gourd/root vegetable diet was most highly associated with carbohydrate intake, which may increase atherosclerosis when such foods take the place of unsaturated fats( Reference Mozaffarian, Rimm and Herrington 57 ).
Study limitations and strengths
To our knowledge, the present study is the first of dietary patterns and IMT in a low-income country. Strengths of our study include its large, population-based sample size as well as use of a validated FFQ and standardized IMT protocols. Additionally, religious motivations lead to restricted alcohol consumption in Bangladesh and remove a potential confounder, as does the fact that a particularly small segment of the study population reported use of any prescription drugs at baseline. Excluding contraceptives, the most common prescriptions were ulcer/gastro-oesophageal reflux medications (4·31 % of participants), antibiotics (1·43 %) and beta blockers (0·81 %)( Reference Scannell, Argos and Pierce 58 ). Limitations include the cross-sectional study design, with a time gap between FFQ administration (during baseline enrolment periods: 2000–2002 or 2006–2008) and IMT measurement (2010–2011). Although the absence of IMT measurement at recruitment makes it impossible to firmly establish a temporal relationship, atherosclerosis is a chronic disease that is likely to develop over a long period of time as a result of past and persistent exposure. Our FFQ was validated to measure long-term usual levels of common foods and nutrients( Reference Chen, Ahsan and Parvez 24 ); as such, further adjustment of our primary outcome models for recruitment period did not affect the significance of association in either the gourd/root vegetable or balanced diet (P<0·01). An additional limitation is that although diet patterns can be argued to be a more effective, realistic method of assessing food consumption, their determination still involves some subjectivity. The findings presented here are also limited by our lack of data on physical activity and lipids, known determinants of CVD. We, likewise, have no data regarding the menopause status of female participants. Lastly, certain characteristics of the study population also limit generalizability, as the population possessed a mean age of 39 years at enrolment (range 17–65 years) and a mean BMI of 19·95 kg/m2 (range 11·70–35·58 kg/m2).
Conclusions
In conclusion, we found a gourd/root vegetable diet to correlate with increased carotid IMT, while a balanced diet was associated with a decrease in IMT. These findings underscore the utility of whole dietary pattern analysis in identifying nutrient deficiencies associated with population-specific diets, and could be useful in determining the most effective strategies for chronic disease prevention in this and similar populations.
Acknowledgements
Financial support: This research was supported by the National Institutes of Health (grant numbers R01 ES017541, P42 ES010349 and P30 ES000260) and in part by the National Center for the Advancement of Translational Science (NCATS), National Institutes of Health (grant number UL1 TR000038). The National Institutes of Health had no role in the design, analysis or writing of this article. Conflict of interest: None. Authorship: Authorship for all individuals listed on this manuscript was justified through participation in the following: conception and design (T.R.M., Y.C., F.P., T.I., A.A., R.R.P., I.S., G.S., H.A.), data collection (F.P., T.I., A.A., R.R.P., I.S., G.S.), data management (F.P., T.I., A.A., R.R.P., I.S., G.S.), analysis and interpretation of data (T.R.M., Y.C., F.W., T.R., R.T.D., M.D., H.A.), drafting of the manuscript (T.R.M., Y.C.), revising the manuscript critically for important intellectual content (all authors) and final approval of the manuscript (all authors). Ethics of human subject participation: Informed consent was obtained from study participants; study procedures were approved by the Ethical Committee of the Bangladesh Medical Research Council and the Institutional Review Boards of Columbia University and the University of Chicago.