Do dietary patterns explain high prevalence of cardiovascular risk factors among Pakistani urban adults? A cross-sectional study

Background

The prevalence of cardiovascular disease (CVD) in South Asia is higher than in any other developing countries. The diversity of diets in populations among developing countries may be one explanation for the differences in CVD. This study was carried out to explore the association between dietary patterns and the presence of cardiovascular risk factors among Pakistani low income urban adults.

Methods

Socio demographic, physical activity and dietary information was collected from 1546 Pakistani subjects aged ≥40 years from the Control of Hypertension and Risk Attenuation (COBRA) study. Anthropometric, clinical and biomarker measures were assessed. Cluster analysis was used to identify dietary patterns from a food frequency questionnaire (FFQ) and multinomial regression was employed to investigate the association between dietary clusters and cardiovascular risk factors, using one of the dietary clusters as a reference category.

Results

The most prevalent CVD risk factors among participants were elevated low density lipoprotein- cholesterol (LDL-C) (69.8 %) and systolic blood pressure (SBP) (68.2 %), followed by central obesity (57.1 %), low levels of high density lipoprotein-Cholesterol (HDL-C) (56.3 %), overall obesity (46.0 %), high total cholesterol (32.3 %), and elevated fasting blood sugar (FBS) (34.9 %). The cluster analysis generated 3 non-overlapping diet patterns. Cluster 1(Traditional Pakistani Diet), was dominated by fruits, vegetables, milk products and chicken, included participants with high mean body mass index (BMI), waist circumference (WC), HDL-C and low mean SBP. Cluster 2 (Moderate Diet) reflected a moderate intake of most food items and included participants with significantly higher mean SBP. Cluster 3 (Fatty Diet) was characterized by high intake of beef, whole milk, paratha and lentils and those following this cluster had a low mean HDL-C and high SBP. In analyses controlling for age, gender, tobacco use and physical activity, participants in the Traditional diet cluster were more likely to be overweight (OR 1.39, 95 % CI 1.08 to 1.78) and high central adiposity (1.33, 95 % CI 1.04 to 1.71) than participants in the Moderate diet cluster, though less likely to have elevated SBP (OR 0.67, 95 % CI 0.51 to 0.86). LDL-C levels were higher in both the traditional Pakistani diet and Fatty diet cluster compared to the Moderate diet cluster. Low HDL-C was also more prevalent among the Traditional Pakistani diet cluster.

Conclusion

Among Pakistani population discernable diet patterns can be derived using clusters analysis. CVD risk factors prevalence differed by cluster membership, though relations for specific CVD risk are not consistent across clusters.

Globally cardiovascular diseases (CVD) are a major contributor to chronic disease burden [1]. South Asians, in particular have a higher rate of acute myocardial infarction, and develop CVD at a younger age than individuals from other countries [2, 3]; CVD is becoming a leading cause of premature adult mortality in the region [4], mainly due to increased occurrence of CVD risk factors including central obesity, hypertension, and dyslipidemias among South Asians [3, 5]. Moreover, compared to other groups within the South Asian population, the prevalence of CVD has been found to be as high as 26.9 % among Pakistanis who are 40 years and older [2].

The quality of diet has been associated with the development of CVD [6]. Diet is a complex exposure variable and individual dietary factors (nutrients) are highly correlated with one another. In addition to studies evaluating individual nutrients, dietary patterns are now being employed to better understand diet and disease relationships. Dietary pattern analyses integrate consumption of individual foods and nutrients into meaningful groups or patterns [7–9]. Cluster analysis is a data-driven analytical method useful to identify dietary patterns within a homogenous group of people; which are then treated as independent variables for analysis of associations with markers of disease risk and health outcomes [10].

Studies using CVD biomarkers can provide important evidence on the relation of dietary patterns and disease risk [11, 12]. For example, risk of CVD can be largely predicted by measurement of plasma lipids. Large prospective studies from the United States [7, 13, 14] and Europe [15, 16] linking diet patterns and triglycerides, cholesterol, low-density lipoprotein, and high-density lipoprotein cholesterol have shown inconsistency in the effect of overall diet on disease risk.

Given the well-recognized link between diet and CVD, an important step is to identify dietary patterns that are associated with CVD in specific populations [17, 18]. As there are variations and differences in dietary intakes across different populations, it is likely that risk factors of CVD differ among Pakistani people as compared to other populations [5, 19]. Limited research has been conducted in Pakistan on dietary patterns and CVD risk factors [20, 21]. Previously we reported an associations between dietary pattern derived from factor analysis and socio-demographic, anthropometric, lifestyle characteristics and hypertension in Pakistani adults enrolled in the Control of Hypertension and Risk Attenuation (COBRA) study [22, 23]. In the current study, our aim is to derive dietary patterns using the cluster analysis method and investigate how these diet patterns are associated with CVD risk among Pakistani people aged 40 years and above.

Study design and participants

This study was a secondary analysis of the data from the COBRA trial of 2004–7 that tested the effectiveness of home-based health education to prevent and control high blood pressure in low income urban adults in Karachi, Pakistan [24]. The sampling technique and participants particulars of the COBRA trial have been described in greater detail elsewhere [25]. Briefly, the sampling involved randomly selecting 12 out of 4200 lower middle income (average household income less than $70/month) household groups. Each group comprised of approximately 250 houses with an average household size of 7 individuals. The study population was further divided into three sub groups; children aged 5–14 years, younger adults aged 15–39 years, and older adults aged 40 years and above. At enrollment, all individuals aged 40 years or older were invited to participate in the study and 3143 agreed. After excluding 1418 individuals with self-reported history of hypertension, heart disease, kidney disease, diabetes and stroke (because their condition may have altered their eating patterns) and those (n = 179) who refused biochemical testing, 1546 participants’ data were analyzed.

Assessment of diet

Baseline dietary intake was measured using a modified version of the Harvard food frequency questionnaire (FFQ) [26] by trained interviewers. Modification of the FFQ was based on the researchers’ knowledge about Pakistani dietary habits and their previous work experiences with the study population. In addition, they reviewed available literature and information gathered through dietary surveys, which reflected the common dietary habits of the low income, Pakistani population. Using preliminary information, a 33-item FFQ was formulated to assess the intake of specific food items over the previous month without specifying the quantity or portion size. The FFQ contained questions about the frequency of intake of different food items such as type of locally used bread, meat types, cooked and uncooked vegetables, fruit and juices, milk and local dairy products, and cooked food dishes normally consumed by low income population and that may pose a risk to hypertension. The response options in the FFQ ranged from none to less than once per month, number of times per month, number of times per week, or number of times per day. The questionnaire was pretested for input on content, language clarity, and layout and revisions were made accordingly for final data collection from the study cohort. All weekly food frequency variables were used to derive dietary patterns via cluster analysis.

Assessment of CVD risk factors

Cardiovascular risk factors including body weight, height, WC and blood pressure were measured and recorded by trained field workers and technicians using calibrated equipment and standardized protocols described in detail elsewhere [24]. Body mass index (BMI) was calculated as weight/height² (kg/m²). WC was measured at the mid-point between the iliac crest and the lower rib margin. SBP and diastolic blood pressure (DBP) were measured with a calibrated automated sphygmomanometer (Omron HEM-737 Intellisense Blood Pressure Monitor) by trained staff after the participants sat quietly for more than five minutes. Three consecutive readings were recorded with a five minute gap between each reading. The average of the last two blood pressure readings were recorded for each subject. Blood samples were collected after 10 h of overnight fast for the estimation of serum total cholesterol (TC), LDL-C, HDL-C, TG and FBS. Lab analyses were performed at the Aga Khan University Hospital using the CHOD-PAP enzymatic end point method with the Colorimetric/Selectra Hitachi 912 analyzer (Roche Diagnostics, Switzerland) [27]. Lifestyle factors including physical activity was assessed with International Physical Activity Questionnaire (IPAQ) [28]. Average weekly duration of walking, moderate and vigorous activity was reported and computed as a metabolic equivalent (MET) score using Ainsworth's compendium [29]. Smoking status was assessed by gathering information on participant use of different forms of tobacco (cigarettes/cigars/biddi, pipe/huqqah, niswar, and tobacco chewing with or without pan leaf). Participants were classified if they had never used tobacco in any form (never) versus those who used tobacco in the past but had stopped use in the last year (past) and those who were currently using tobacco in some form (current) user .

Each of the CVD risk factors were defined according to the evidence-based guidelines in the prevention of CVD and established criteria were used to evaluate each participant for all risk factors. Participants were classified as: obese, kg/m² (BMI ≥ 25, Asian criteria) [30]; having central obesity (WC of ≥90 cm in males or ≥80 cm in females, Asian-specific) [24]; having elevated SBP (≥120 mmHg) [31]; high cholesterol (≥200 mg/dl); high LDL-C (≥100 mg/dl); high TG (≥150 mg/dl); low HDL-C (<40 mg/dl) and elevated FBS (≥100 mg/dl) [32, 33]. Smoking and physical activity risk was determined if participants reported use of tobacco in any form and were physically inactive. Each cardiovascular risk factor was converted into dichotomous variables (categorical) indicating their presence or absence.

Statistical analysis

Data were analyzed using SPSS, version 16.0 (SPSS Inc, Chicago, IL). Gender differences in basic sociodemographic and other baseline variables were assessed using independent sample t-tests and the values were reported as means ± SD if the data were continuous (age, BMI, WC, lipid biomarkers, and blood pressure). Categorical variables (educational status, tobacco use, physical activity) are presented as absolute and relative frequencies and their significance across gender were tested by chi-square test.

We then performed cluster analysis. Since characteristics of data were exploratory therefore we decided to use K-means cluster analysis, which identifies distinct non-overlapping groups, based on Euclidian distances [34]. FFQ data were converted to a mean of total servings from 33 food items and entered into a cluster analysis procedure. The group of food items that were consumed in a similar manner was identified; SPSS used 10 iterations to centralize the mean level of clusters. Each cluster solution was evaluated according to the degree of separation between clusters using an F statistic. Cluster solutions of 4 were evaluated and the resultant 3 cluster solutions showed the best proportion of subjects in each cluster and separation between them. The clusters were labeled arbitrarily. Individuals in a given cluster were thus distinctly similar to one another and different from individuals in other clusters with respect to differences in mean food intakes [10]. To compare the means of CVD associated variables across clusters, one-way analysis of variance (ANOVA) was carried out followed by post-hoc analysis using Tukey HSD test and p-value <0.05 was considered significant. Multinomial logistic regression was conducted to assess association between cardiovascular risk factors and dietary patterns while controlling for age, gender, physical activity and tobacco use. All reported p-values are based on two-sided tests and compared to a significance level of 5 %.

Participant characteristics

Complete data on 1546 subjects were available, representing 89 % participation rate. Table 1 summarizes the physical characteristics of the subjects. Men represented 54.9 % of the total sample with no significant differences in age between men and women. Women were less likely to have formal education as compared to men. The mean BMI in women was significantly higher than that of men, whereas SBP was higher in men as compared to women. For biomarkers, significant differences between men and women were seen for triglycerides and HDL-cholesterol. No other gender differences in blood lipids were found in the sample. When we analyzed the proportion of participants having various CVD risk factors, there were none without any CVD risk factors. Risk of central and overall obesity as defined by WHO Asian criteria was higher in women than men. High SBP was more prevalent in men than women. For abnormal lipid biomarkers as a CVD risk factor, men had lower HDL-C than women whereas elevated FBS was more prevalent among women than men. Lifestyle related CVD risk factors were also found to be significantly different in males and females, with a higher proportion in men who used tobacco and higher proportion of women who were physically active (Table 1).

Characteristics of clusters & CVD risk

Three discernable clusters or diet patterns were identified, though there was no clear “healthy" or "unhealthy" diet cluster in this population. A total of 447 subjects were in diet cluster 1 consuming mostly fruits, vegetables, milk without cream, yogurt and chicken. This cluster was named arbitrarily as the ‘Traditional Pakistani diet pattern’. Cluster 2 diet pattern was followed by the largest group of subjects, consuming moderate amount of most food items and therefore we labeled the cluster as 'moderate diet pattern'. In cluster 3, people were mainly consuming foods high in saturated fats such as paratha, whole milk, beef etc. and therefore we labeled this cluster as 'fatty diet pattern'. The naming of clusters was done arbitrarily and for ease in the discussion. Table 2 shows the frequency of consumption of each food items in each of the three diet clusters. Each diet cluster was further examined for CVD associated risk factors. Significant mean differences were seen for BMI, waist circumference, SBP and HDL-C across the three diet clusters (Table 3), with highest mean BMI, WC, and HDL-C among cluster 1 subjects.

Association of diet clusters with CVD risk factors were assessed by multinomial logistic regression analysis. The probability of being obese, having central obesity, elevated SBP were assessed across diet clusters. The reference category for the regression model was the moderate diet pattern (cluster 2). The model was adjusted for age, gender, tobacco use and physical activity. Table 4 shows that if a subject was following ‘Traditional Pakistani’, then he or she was more likely to be obese as measured by BMI (OR 1.39, 95 % CI 1.08 to 1.78; P < 0.01) and waist circumference (OR 1.33, 95 % CI 1.04 to 1.71; P < 0.01) as compared to other dietary clusters. However, subjects in this diet pattern (cluster1) were less likely to have a high SBP (OR 0.67, 95 % CI 0.51to 0.86; P < 0.01) as compared to others. Those subjects following ‘Traditional Pakistani’ diet pattern (OR 1.43, 95 % CI 1.10 to 1.86 P < 0.01) and fatty diet pattern (OR 1.53, 95 % CI 1.15 to 2.03 P < 0.01) were more likely to have higher LDL-C as compared to reference diet pattern (cluster 2). Subjects were more likely to have a high TG (OR 1.30, 95 % CI 0.99 to 1.63 P 0.054) and low HDL (OR 1.34, 95 % CI 1.03 to 1.73 P < 0.03) in cluster 1 characterized by as compared to others. Subjects following fatty diet pattern (cluster 3) were 26 % (OR 0.74, 95 % CI 0.56 to 0.98 P < 0.03) less likely to have elevated FBS than those following other dietary patterns. No other significant associations were seen for CVD risk factors in our data set.

South Asians are at higher risk of cardiovascular disease due to an increasing prevalence of its risk factors amongst them [35]. This cross-sectional study was undertaken to examine the association of diet patterns derived by cluster analysis with the prevalence of CVD risk factors (mainly obesity, lipid biomarkers, SBP and FBS) in Pakistani adults aged 40 years and older. Identification of overall diet trends that incorporates the interaction of different food items in a specific population is a more appropriate method to examine the relationship between diet and disease risk than focusing on individual foods or nutrients [7, 36]. Dietary pattern studies tend to identify at least one pattern or cluster that can be defined a healthier eating pattern, however, in some populations a typical healthy or unhealthy diet pattern may not exist but dietary behaviors might be characterized by different types of foods categorized as healthy or unhealthy within the similar diet patterns.

In Pakistan, limited research has been done to identify the different type of diets that are normally consumed by Pakistani people. In addition, currently no dietary guidelines, recommendations, or evidence-based dietary indices have been established by health and nutrition authorities in Pakistan. Therefore, studies using data driven methods including factor and cluster analysis are essential to identify existing dietary patterns among Pakistani people and inform public health initiatives in Pakistan. In our study, men and women following traditional Pakistani diet pattern (cluster 1) were characterized by high intake of foods including milk and milk products, potatoes and other vegetables, fruits, Pakistani bread, and bakery products with highest risk for obesity as compared to other diet patterns. Though, this diet pattern has some foods that are relatively healthier but because the usual intake in this group was highest for most foods, thereby increasing the total energy intake that may have been resulted in higher prevalence of overweight and obese adults in this pattern. Similar results were also observed in previous research studies [34, 37].

The relationship between dietary patterns and CVD risk could be manifold and complex in different populations [12, 38]. A large body of scientific evidence from observational studies [39], clinical trials [40], and genetic studies [41] have shown that lipids and glucose biomarkers may have causal relationship with CVD risk, though results have been inconsistent [7, 15]. Dyslipidemia, manifested by abnormal lipids with increased TC, LDL-C, and TG accompanied by a decrease in HDL-C, is considered a traditional risk factor for CVD [42]. It was difficult to compare our findings with earlier studies conducted in other populations because we were unable to identify a very distinct healthy or unhealthy diet pattern with which we could noticeably associate biomarkers of CVD risk. However, those following traditional Pakistani and fatty diet patterns (cluster 1 and 3) were relatively at a higher risk for LDL- cholesterol as compared to cluster 2. This could have been due to more food in general and high fatty food intake as reflected by these clusters. Similarly the risk for high TG was also observed in traditional Pakistani diet pattern (cluster 1). Previously Jafar has shown that factors including obesity, caloric-dense diet and sedentary lifestyle are likely to be associated with increased prevalence of CVD risk in Pakistani women [43]. In the current analysis those following pattern 3 were not consuming high levels of fatty foods including eggs, milk desserts, ice cream, and butter but had more CVD risk factors as compared to reference. The FFQ used in the assessment of food intake included a limited list of fatty foods reflective in a typical Pakistani diet, as the intent of the COBRA study was to determine only those types of foods that may pose a risk to hypertension. In addition, the FFQ was unable to capture the type and amount of fat used in cooking, and previous studies from Pakistan has shown an increased intake of saturated fats in daily cooking that may perhaps have contributed to an increased risk for those following pattern 3 [2, 43, 44].

Most data linking dietary patterns with biomarkers of CVD risk comes from observational studies. The Lyon Diet Heart Study, a randomized control trial was unable to show any association between lipid biomarkers of CVD risk (HDL-C, LDL-C, TG, and lipoprotein-a) and intake of diets characterized by fruits, vegetables, olive oil, and hard bread [45]. Thus, these results suggest that risk factors other than serum lipids may mediate risk between dietary patterns and CVD risk. We recommend caution in interpreting the relationship between the dietary patterns and biomarkers as ours was a cross-sectional study.

In our study, the prevalence of subjects with elevated SBP was found to be relatively high (68.2 %), though subjects who self-reported hypertension were excluded. Distributions of these participants were highest (71.8 %), across moderate diet pattern, despite the fact that the frequency of intake of many food items was much less than in other two clusters. Also notable was the fact that certain Pakistani foods that are supposedly proposed to increase blood pressure for example paratha, tandori nan, salty lassi, biryani etc. were eaten much less frequently in this group as compared to other cluster groups. In contrast, subjects in cluster 1 showed a protective effect on SBP risk that cannot be explained by more frequent consumption of above food items. This unusual observation could be explained by the fact that intake of fruits and vegetables along with some unhealthy food were seen in this group (Traditional Pakistani diet, cluster 1). Nevertheless, similar results have also been observed by Sadakane et al.[46], who has reported that the Western diet pattern was not associated with higher systolic and diastolic blood pressures in the Japanese population. We recommend further exploration of these findings in order to improve our understanding about the association between western type of diet in a Pakistani population and blood pressure.

This study is notable because of its large sample size and a high response rate (88.38 %). The extensive information gathered from the study subjects makes our study findings more informative and may have substantial implications for other similar populations in South Asia.

However, the study has several limitations. First, the results cannot be generalized to other regions of Pakistan as the study group was from one major city only and therefore cannot be assumed to be representative of all Pakistani people. However, participants from the COBRA study reflected almost all major ethnic groups living in Pakistan. Second, it was a cross sectional study, therefore, causality of the relationship between the dietary patterns and CVD risk factors cannot be fully ascertained. In addition to the nature of the study, recall bias was another concern because abnormal lipids may have changed participants’ diet patterns or led to a biased reporting of dietary habits. This issue, however, was mitigated by the exclusion of patients with self-reported diabetes, hypertension, stroke and heart disease from the study analyses. Third, we used a modified FFQ with limited food lists which was not a validated questionnaire for Pakistani population; though it was pretested in a similar sample population prior to its administration in terms of being understood by the population, and for capturing the typical dietary intake of the study population was done. Also, the FFQ did not require the participants to report the portion size or amount of food; therefore it precluded us from adjusting our regression model for total energy intake. Though, an alternative adjustment was done by adjusting the regression model for physical activity. It was difficult to quantify misreporting based on the available data that could be due to several reasons that we were unable to capture due to methodological limitations of secondary analysis. Fourth, since this was a study in which biomarkers were assessed, we need to consider many factors that may affect the appropriateness and sensitivity of biological samples in epidemiologic studies of dietary intake. A single sample may not have been be a representation of long term intake, especially if there is lot of variation in an individual’s diet. In addition, there are collection and measurement issues that may have affected the results in the current study [12].

Cluster analysis is a useful exploratory data driven technique for deriving dietary patterns in populations where fundamental dietary patterns do not exist. While identifying data driven or posteriori dietary patterns, many arbitrary but important decisions are made by the researcher, such as consolidation of food items into food groups, the number of patterns to derive from the data, labeling of patterns etc. This can make it difficult to compare and interpret the findings with other studies. However the exposure groups derived from cluster analysis are homogenous and non-overlapping, which is helpful in studies explaining relationship between diet patterns and risk for diseases [10]. As dietary guidelines for Pakistan get established, future studies of guideline-based indices and CVD risk will be possible and well warranted.

The present study shows that some of the CVD risk factors including obesity and hypertension are highly prevalent among low income urban Pakistani adults. The study further provides some evidence that CVD risk factors could be linked with certain type of diet patterns but not consistently for all CVD risk factors as well as in all population groups. However, these results are important for Pakistan as people here are already susceptible to a high prevalence of obesity, hypertension and dyslipidemia [47]. Effective Public health strategies and multiple lifestyle change initiatives are required against obesity and hypertension among Pakistani adults to avert rise in CVDs.

ANOVA:

One-way analysis of variance

BMI:

Body Mass Index

COBRA:

Control of blood pressure and risk attenuation

CVD:

Cardiovascular disease

DBP:

Diastolic blood pressure

FFQ:

Food frequency questionnaire

HDL-C:

High density lipoprotein-cholesterol

IPAQ:

International physical activity questionnaire

LDL-C:

Low density lipoprotein-cholesterol

MET:

Metabolic equivalent

OR:

Odds ratio

SBP:

Systolic blood pressure

TC:

total cholesterol

TG:

Ttriglycerides

FBS:

Fasting blood glucose (FBS)

WC:

Waist circumference

WHR:

Waist hip ratio

The authors thank members of the Hypertension Research Group (HRG) who were responsible for the COBRA study (clinical trial NCT00327574). We also sincerely thank Adnan Ali for his contribution in the statistical analysis of this research paper.

Funding

This study received no financial support and was part of the doctoral dissertation of Nilofer F. Safdar at the Department of Nutrition, School of Public Health and Health Sciences, University of Massachusetts, Amherst USA.

Availability of data and materials

Most of the data supporting our findings is contained within the manuscript. However, detailed SPSS data spread sheets are available with the authors.

Authors’ contributions

NFS conceptualized the research idea, study design, analysis, interpretation of the data, and drafting of the manuscript. EBJ provided guidance for study design, analysis, interpretation of the data and review of the manuscript. LC provided input on early drafts. TJ designed and conducted the COBRA study, provided guidance in conceptualizing the idea for the manuscript. N.C was closely involved in the design, analyses, and review of the manuscript. All authors have read and approved the final version of the manuscript.

Competing interests

The authors declare that they have no competing interest.

Consent for publication

Not applicable.

Ethics approval and consent to participate

The COBRA study was conducted in accordance with the principles of the Declaration of Helsinki and was approved by the Aga Khan University Ethical Review Committee. Written consents were obtained from the COBRA participants after explaining the objectives of the study. A separate approval was further obtained from the Human Subjects Committee and Internal Review Board (IRB) of the University of Massachusetts, Amherst, USA to use deidentified COBRA data for the secondary analysis for academic research and writing of manuscripts.

Affiliations

1. School of Public Health, Dow University of Health Sciences, KDA Scheme 33, Gulzar-e-Hijri, OJHA Campus, Suparco Road, Karachi, 75270, Pakistan

   • Nilofer F. Safdar

2. Department of Biostatistics and Epidemiology, University of Massachusetts, Amherst, MA, 01003-9304, USA

   • Elizabeth Bertone-Johnson

3. Department of Nutrition, University of Massachusetts, Amherst, MA, 01003-9282, USA

   • Lorraine Cordeiro
   •  & Nancy L. Cohen

4. Department of Community Health Sciences, Aga Khan University, Karachi, Pakistan

   • Tazeen H. Jafar

Corresponding author

Correspondence to Nilofer F. Safdar.

Reprints and Permissions