Association between Dietary Diversity and Food Insecurity in Urban Households: A Cross-Sectional Survey of Various Ethnic Populations of Karachi, Pakistan.

BACKGROUND AND
PURPOSE: The world is experiencing a nutritional transition, yet some regions of South East Asia are still living under the shadows of catastrophic nutritional indicators. Pakistan bears a high dual burden of both communicable and chronic diseases. However, a major contributing factor of both is poor diet and nutrition. The causal pathway of stunting, underweight, wasting, and micronutrient deficiencies has beginnings in less dietary diversity and food insecurity. Current literature on Pakistan and other South Asian countries regarding food insecurity and dietary diversity largely focus on women and children aged less than 5 years; however, ethnic and cultural dynamics of the populace concerning their food and dietary practices in Pakistan have not been well explored. It is essential to have a clear insight into the factors involving malnutrition in different clusters, especially different cultures and ethnicities to target specific areas of interventional strategies versus a uniform approach for all. This study aimed to explore the level of food insecurity and dietary diversity in all major ethnic groups of Karachi, the largest city of Pakistan.
METHODS: Using the multistage random sampling, 535 households from five major ethnic groups residing in Karachi were included in the study. To measure food insecurity and dietary diversity (primary and secondary outcome, respectively), we used a FAO standardized questionnaire (version 3, 2007) that was translated first into local languages, piloted, and employed to collect information.
RESULTS: Severe food insecurity levels had an inverse relationship with the household dietary diversity, which persisted even after adjusting for other variables (β=-0.31, 95% CI=-0.65, -0.07). There was no statistically significant association of age, marital status, or education level on food diversity. The most food secure cultural people were those whose ancestors had migrated from India and were Urdu speaking communities, while Sindhi speakers were the least food secure community.
CONCLUSION: It is imperative to investigate more on the cultural causal factors leading to food insecurity to address the root causes of malnutrition and design new cultural-specific interventions that should be employed in large urban centers where different communities reside together.

Introduction

Changing global demographics and rapid urbanization are key factors that shape modern dietary patterns.1 Population growth in urban zones demands greater food supplies and availability, hence these urban centers are prone to have an increased cost of food. Additionally, rising inflation can have an inverse impact on dietary diversity and food security.1 Ending all forms of malnutrition and providing access to safe, adequate, and healthy food to all people worldwide is one of the aims of the United Nations’ Sustainable Development Goals (SDGs). Despite global efforts made to reduce malnutrition, which have lowered overall deaths due to nutritional causes, child deaths owing to undernutrition are still around 45% worldwide, and the prevalence of undernutrition related mortality remains considerably high in South East Asia.2

Pakistan, a low-to-middle-income country, which is striving to deal with substantial food security challenges in addition to other effect multipliers such as spurt population growth and economic crisis, is also experiencing a mass urbanization crisis.3 Pakistan is not merely self-sufficient in most of the food sources, but is also a net exporter of several agricultural products, yet it faces food insecurity at high levels. The calamity is evident by the “World Food program” report for Pakistan, stating that 18% of the Pakistani population are undernourished and this prevalence of undernourishment is also ranked “Serious” by the Global Hunger Map threshold.4 In Pakistan, almost 40% of children aged less than 5 years suffer from stunting (a key and extremely sensitive indicator of chronic malnutrition) and every third child is underweight (National Nutritional survey 2018). This indicates an alarming situation and mandates an immediate evidence-based comprehensive nutritional policy.5

Dietary diversity is linked with the access dimension of household food security. It serves as a measure of the nutritional quality of the diet as a variety of food items in the diet that is required to ensure adequate consumption of essential nutrients.6 M’Kaibi et al7 showed a positive relationship of dietary diversity with normal growth patterns.

Studies have shown that, besides specific fiscal and nutritional variables, multiple social factors are also strongly associated with food insecurity in Pakistan. These factors include maternal education, paternal employment, and poor healthcare system, vulnerable population, knowledge, awareness at household levels, and inequitable food supply.8,9 Almost half of the Pakistani population (48%) suffers from food insecurity, and the prevalence is much higher in rural areas and suburbs of large cities. It is also a major concern that 42% and 27.3% of Pakistani women of reproductive age suffer from iron deficiency anemia and vitamin A deficiency, respectively.5,10

With an increased cost of healthcare, high inflation rate, treating complications, and loss of productive life years, more importance should be given to preventive health, including the indispensable role of nutrition in reducing the risk of mortality and morbidity.11 The Food and Agriculture Organization of the United Nations (FAO) defines food security as “When all people, at all times, have physical, social, and economic access to sufficient, safe, and nutritious food, which meets their dietary needs and food preferences for an active and healthy life”.12 Key facets of food security include availability, accessibility, affordability, and acceptability.13

Research has demonstrated that an average urban Pakistani diet consists of a carbohydrate rich staple based diet.14 Staple foods provide a greater level of energy, however may not improve nutritional outcomes if not supplemented with other sources such as vitamins and minerals.14 Multiple studies demonstrate that increasing dietary diversity and food variety can successfully decrease food insecurity and the choice of food consumption in households, and these are driven by cultural choices and family preferences such as taste, quality, health concerns, and monetary elements.15,16 DiSantis et al17,18 and Burns et al have demonstrated that family preferences and price of the food are among the most important key factors associated with the number of food items consumed in a household diet in a low income country.

In Pakistan, during the last 60 years, there has not been a significant change seen in the indicators of malnutrition despite huge investments in different nutrition interventions. This has been reflected in numerous studies including various national nutrition surveys.19 Hence, there is a need to find the critical “missing link” which has been left unaccounted for so far. This critical link could be the ethnic dimension of nutrition. Most of the nutrition-specific and sensitive programs have been designed and relied on numbers of the malnourished population without understanding the cultural dynamics of the populace concerning their food and dietary practices. However, ethnic subgroup differences that could influence food insecurity and dietary diversity in different populations have not been explored so much, specifically among South Asians. The population of Pakistan includes five major ethnic subgroups coming from different parts of Central and South Asia: the Urdu speakers, the Punjabis, the Sindhis, the Baluchis, and the Pashtuns. These groups have distinct differences in cultural practices and values, dietary habits, health beliefs, and behaviors that could make them susceptible to food insecurity and dietary diversity, even within the same geographic location. It is also known that the lack of intermingling of these ethnic groups can further contribute to the relatively well-preserved sociocultural and dietary practices of each ethnic subgroup that could either protect or prevent food security within each ethnic group.

The main objective of this study was to conduct a preliminary quantitative survey of the population residing in Karachi to identify food insecurity (primary outcome) and dietary diversity (secondary outcome) among five ethnic groups in Pakistan to capture cultural relevance to the nutritional situation of households.

Materials and Methods

Study Design, Setting, and Population

This is a community-based, cross-sectional study that was conducted between August 2019 and July 2020 in five major municipal areas of Karachi. The five ethnic groups mentioned above represent more than 95% of the country’s population and they reside in clusters in different regions of Karachi, hence besides being the largest, it is the most culturally diverse city of Pakistan, also known as Mini Pakistan, with a population over 16 million.20 The study population was drawn from various catchment areas of Karachi including; Patel Para, Malir, Liyari, Sheerin Jinnah Colony, and North Nazimabad, that have a concentrated population of Punjabi, Sindhi, Balochi, Pushto, and Urdu speaking communities, respectively.

Sample Size and Sampling Technique

The study sample size was estimated using Openepi online calculator using a 60% proportion of food insecurity,21 a 4.5% margin of error with 95% confidence intervals, and a design effect of 1.0 for random sampling. The sample size came out to be 456. We sampled 535 households after inflating the sample size by 15% to account for the non-response rate and missing information.

Within the selected communities/cultures, a two-stage household sampling procedure was employed allowing for a partially convenient and randomized selection process. In the first phase, five districts of Karachi were chosen and, within those five districts, five areas concentrated with each ethnicity were purposely selected for the advantage of an improved sampling frame, followed by second phase sampling, which involved selecting households. A total of 535 households were approached from the five selected areas, with a systematic sampling technique starting with a random start.

Data Collection Procedure

The first house in “Union council” (smallest administrative unit of a district in Pakistan) within communities was selected by computer generated house or apartment number from the available registry of addresses from population department, the government of Sindh, Pakistan. The interviews were conducted by 15 trained professionals, who were residents of their selected areas (five from each area). However, team supervisors were present to monitor the activities during data collection. In case of denial or non-availability of the member of a household that had been selected, the next house on the right side (while standing facing the door of the first house) was selected and in rare instances where the next house was also not accessible to survey (denial or non-availability), the house after that was included.22 Forty-six households denied access or were not available and we moved to the second house, however, only three houses were required to be approached in the case of double denial. Hence, the first house response and participation rate was 91.4%.

Inclusion Criteria

This consisted of the person responsible for planning meals in the household (decision-maker). However, a person who prepared or cooked the meal if he or she was not the person who planned the meal was excluded. All study participants were females, as they are considered the decision-makers for food preparation in their homes in Pakistan; however they represented their entire households.

Data Collection Tools

To identify the prevalence of selected household’s food insecurity situation and dietary diversity and food consumption patterns, the FAO standardized questionnaires were used. These were Household Dietary Diversity Score23 and Household Food Insecurity and Access Scale (HFIAS).24 Both tools were standardized and have been used in many low-to-middle income countries in settings where malnutrition is a major concern.25–27

Face Validity

The questionnaire was translated into the local languages and all questions were discussed independently with 20 households (four from each ethnicity) from houses other than the selected five areas aiming to inquire regarding the clarity and understanding level of the questions and their relevance to the objectives of the study. Only two participants responded with multiple interpretations. Modifications were made in the questionnaire for cultural validation, by five experts of the field. Four questions were modified, and contextually appropriate examples were included.

HFIAS

The HFIAS questionnaire consists of a list of nine specific questions about apprehension and availability of and accessibility to foods for the household during the past month.24 The questionnaire was translated into five different languages including Urdu, Sindhi, Punjabi, Balochi, and Pashto for subjects included in the study. A standardized procedure for scoring was used, where zero was attributed if the event described by the question never occurred, 1 point if it occurred 1 or 2 times during the previous 30 days (rarely), 2 points if it occurred 3–10 times (sometimes), and 3 points if it occurred 10 times or more (often). For each household, the HFIAS score corresponded to the sum of these points and could range from 0 (food security) to 27 (maximum food insecurity)24 ().

Dietary Diversity Scale

Dietary diversity was estimated by organizing information on household food consumption from a set of 21 food items for 24 hours preceding the survey. These 21 items were clustered into six food groups: Pulses/grains; Fish/seafood; Fruits; Vegetables; Meat; and Dairy products. The six food group score, a measure of dietary diversity, is counted by adding up the value of the items reflecting the six food groups. For each household, the diet diversity can take any value from zero and six. A household would get the value zero if it does not consume any one of the 21 given items and it would get the value of six if it consumed items from all seven food groups.23 Food insecurity is the primary outcome, while dietary diversity is the secondary outcome. Study Covariates: Socio-demographic information about the respondents such as age, marital status, education level, ethnicity, employment status, number of family members and average household monthly income was collected. The number of family members sharing one kitchen was taken as a proxy (substitute) for household composition.

Data Analysis

We used SPSS version 22 for data analysis. Descriptive measures were reported as mean±SD, frequency, and percentages for all study variables. One-way ANOVA (for continuous variables) and Chi-square test (for categorical variables) were used to determine significant differences between the food insecurity groups.

Firstly, an association between food insecurity and diet diversity score was assessed using linear regression analysis. The household’s diet diversity score was taken as a continuous outcome variable and results were reported as a beta coefficient, standard error (SE), and p-value. Subsequently, binary logistic regression was analyzed to determine the association between diet diversity and food insecurity (dichotomous outcome variable: food secure (reference) and food insecure). Results were reported as odds ratio, 95% confidence interval (CI), and p-value. Univariate and multivariate regression analyses were used to estimate the association between food insecurity and diet diversity after adjusting for age, marital status, education level, ethnicity, employment status, and the number of family members. A p-value<0.05 was considered as significant.

Results

Socio-Demographic Characteristics of Respondents

Table 1 shows the basic characteristics of 535 households. The average age of the participants was 32.5 years, and the average household size was 5.8 persons. On average, more than 70% (n=445) of the representative participants had secondary school and higher education, while 16.8% did not attain any formal education. Overall, 63.7% (n=341) of the households were food secure. Approximately 20% (n=109) of households had severe food insecurity issues. Those who consumed more than three food groups were 23% (n=125) of household and almost an equal number reported consuming either three or less of food groups (Table 2). Figure 1 reflects consumption of the depicted number of food groups in an average meal by households.

Table 1

Socio-Demographic Characteristics of Respondents (n=535)

Characteristics	Mean±SD	n	%
Age in years (n=355)	35.7±12.2
Number of family members (n=408)	5.8±2.6
Average HH monthly income (n=352)	56,927.2±38,629.8
Age groups (n=358)
	≤35 years	201	56.1
	>35 years	157	43.9
Marital status
	Single	171	32
	Married	325	60.7
	Divorced/Separated	39	7.3
Education level
	No education	90	16.8
	Primary	63	11.8
	Secondary/Intermediate	178	33.3
	Graduation	138	25.8
	Masters	66	12.3
Ethnicity
	Urdu Speakers	114	21.3
	Sindhi	110	20.6
	Punjabi	106	19.8
	Balochi	103	19.3
	Pashto	102	19.1
Employment status
	Unemployed	176	32.9
	Employed	273	51
	Self-employed	86	16.1

Table 2

Food Security and Dietary Diversity Status among Household Participants (n=535) across Five Ethnic Populations at Baseline

	n	%
HFIAS Scores
Food secure	341	63.7
Mildly food insecure	29	5.4
Moderately food insecure	56	10.5
Severely food insecure	109	20.4
Diet Diversity
Food group score (0)	68	12.7
Food group score (1)	108	20.2
Food group score (2)	109	20.4
Food group score (3)	125	23.4
Food group score (4)	79	14.8
Food group score (5)	37	6.9
Food group score (6)	9	1.7

Figure 1

Dietary diversity across all ethnic groups (number of food groups consumed by the households).

Dietary Diversity

Figure 1 depicts the number of food groups in a household’s average meals, a little less than a quarter of households (23.4%, n=125) consume three groups and 23% (n=125) have more than three groups of food in daily meals. The remaining households (53%, n=285) were consuming meals that were of very low diversity (≤2 food groups as per FAO scale).

Table 3 shows household food insecurity regarding socio-demographic factors. Being single, married, or divorced or separated did not have a large effect on food insecurity status (67%, 61%, and 69%, respectively) with a p-value of 0.014.

Table 3

Socio-Demographic and Lifestyle Factors Associated with Adherence to Major Dietary Patterns among Pakistani Adults (n=535)

Characteristics	Household’s Food Insecurity
	Food Secure	Mildly Food Insecure	Moderately Food Insecure	Severely Food Insecure
	Mean±SD	Mean±SD	Mean±SD	Mean±SD	p-value*
HH Diet Diversity	2.5±1.6	2.7±1.2	2.3±1.0	1.9±1.5	0.006
Average HH monthly income (n=352)	62,059.5±31,496.1	60,448.2±477.5.6	49,742.8±43,138.1	51,250.5±42,997.6	0.063
Number of family members (n=408)	5.6±2.5	5.8±2.5	5.6±3.0	6.6±2.6	0.026
	n (%)	n (%)	n (%)	n (%)	p-value**
Age groups (n=358)
≤35 years	125 (62.2)	5 (2.5)	20 (10.0)	51 (25.4)	0.056
>35 years	106 (67.5)	8 (5.1)	20 (12.7)	23 (14.6)
Marital status
Single	115 (67.3)	4 (2.3)	10 (5.8)	42 (24.6)	0.014
Married	199 (1.2)	21 (6.5)	44 (13.5)	61 (18.8)
Divorced/Separated	27 (69.2)	4 (10.3)	2 (5.1)	6 (15.4)
Education level
No education	23 (25.6)	4 (4.4)	32 (35.6)	31 (34.4)	<0.001
Primary	34 (54.0)	5 (7.9)	13 (20.6)	11 (17.5)
Secondary/Intermediate	123 (9.1)	6 (3.9)	5 (2.8)	43 (24.2)
Graduation	109 (79.0)	6 (4.3)	4 (2.9)	19 (13.8)
Masters	52 (78.8)	7 (10.6)	2 (3.0)	5 (7.6)
Ethnicity
Urdu Speakers	88 (77.2)	3 (2.6)	7 (6.1)	16 (14.0)	0.031
Sindhi	60 (54.5)	6 (6.4)	17 (15.5)	26 (23.6)
Punjabi	66 (62.3)	5 (4.7)	13 (12.3)	22 (20.8)
Balochi	60 (58.3)	11 (10.7)	11 (10.7)	19 (20.4)
Pashto	67 (65.7)	3 (2.9)	8 (7.8)	24 (23.5)
Employment status
Unemployed	128 (72.7)	6 (3.4)	5 (2.8)	37 (21.0)	<0.001
Employed	148 (54.2)	21 (7.7)	47 (17.2)	57 (20.9)
Self-employed (Business)	65 (75.6)	2 (2.3)	4 (4.7)	15 (17.4)

Notes: *p-value calculated using ANOVA analysis. **p-value calculated using Chi-square analysis.

Another significant finding was that the Urdu-speaking households were the most food-secure group (77.2%) as compared to Pushto speaking households (67%), Punjabi population (66%), Sindhi (60%), and Balochi (60%) participants. Those who were unemployed, irrespective of their ethnicity, were significantly more (72.7%, p-value=0.001) food secure than those who had some kind of job (54.2%).

Severe food insecurity levels were found to have an inverse association with the household dietary diversity (β=−0.24, 95% CI=−0.39, −0.09), which persisted after adjusting for other variables (β=−0.31, 95% CI=−0.65, −0.07) (Table 4). There was no statistically significant association of age, marital status, or education level on food diversity.

Table 4

Association Between Household’s Diet Diversity Score and Household Food Insecurity

	Household’s Diet Diversity Score
	Model 1			Model 2
	Β Coefficient	SE	p-value	Β Coefficient	SE	p-value
Household’s food insecurity
Food secure	Ref.			Ref.
Mildly food insecure	0.12 (−0.11, 0.34)	0.11	0.309	0.22 (−0.23, 0.65)	0.22	0.306
Moderately food insecure	−0.04 (−0.22, 0.14)	0.09	0.684	−0.07 (−0.57, 0.19)	0.17	0.665
Severely food insecure	−0.24 (−0.39, −0.09)	0.07	0.002	−0.31 (−0.65, −0.07)	0.12	0.013
Age groups (n=358)
≤35 years	Ref.			Ref.
>35 years	0.06 (−0.07, 0.20)	0.06	0.344	0.04 (−0.19, 0.27)	0.12	0.739
Marital status
Single	Ref.			Ref.
Married	0.11 (−0.01, 0.23)	0.06	0.074	0.11 (−0.14, 0.37)	0.13	0.398
Divorced/Separated	0.16 (−0.06, 0.38)	0.11	0.158	−0.22 (−0.72, 0.27)	0.25	0.370
Education level
No education	Ref.			Ref.
Primary	0.03 (−0.17, 0.24)	0.11	0.733	0.01 (−0.34, 0.37)	0.18	0.940
Secondary/Intermediate	0.01 (−0.15, 0.17)	0.08	0.873	0.16(−0.14, 0.45)	0.15	0.303
Graduation	−0.02 (−0.20, 0.15)	0.08	0.779	−0.05 (−0.37, 0.27)	0.16	0.760
Masters	0.10 (−0.10, 0.30)	0.10	0.322	0.23 (−0.12, 0.58)	0.17	0.197
Ethnicity
Urdu Speaker	Ref.			Ref.
Sindhi	0.01 (−0.16, 0.17)	0.08	0.916	0.07 (−0.20, 0.34)	0.14	0.641
Punjabi	0.02 (−0.14, 0.19)	0.08	0.760	0.10 (−0.16, 0.36)	0.13	0.451
Balochi	−0.06 (−0.24, 0.11)	0.09	0.470	−0.04 (−0.35, 0.27)	0.16	0.818
Pashto	0.03 (−0.13, 0.21)	0.08	0.671	0.05 (−0.22, 0.33)	0.14	0.689
Employment status
Unemployed	Ref.			Ref.
Employed	0.14 (0.01, 0.26)	0.06	0.025	0.17 (−0.04, 0.38)	0.11	0.113
Self-employed	0.06 (−0.10, 0.23)	0.08	0.446	−0.07 (−0.36, 0.21)	0.14	0.608
Number of family members (n=408)	−0.01 (−0.02, 0.02)	0.01	0.874	−0.01 (−0.04, 0.02)	0.01	0.600

Notes: Model 1: Univariate analysis of HH diet diversity score with each study variable. Model 2: Multivariate analysis adjusted for each study variable (Household’s food insecurity, age, marital status, education level, ethnicity, employment status, and number of family members).

Ethnicity and Food Groups

The Urdu-speaking population were the group with the most diversified diet, having on average three or more than three groups in their diet. However the meat group was most commonly consumed by Sindhi speaking population (64.5%). The least abundant food group in all ethnicities was the seafood, which varied from 15.7–26.3%.

Household Food Insecurity

We analyzed the association of various factors on food insecurity through univariate and multivariate models. If the dietary diversity is increased by one unit, there is a likelihood that household food insecurity would reduce by 12% (cOR=0.88, 95% CI=0.78–0.99, p-value=0.042). However, after adjustment of other variables, there would be a reduction of 22% (aOR=0.78, 95% CI=0.66–0.92, p-value=0.002) in food insecurity. Being married or divorced/separated reduced the odds of food insecurity by 43% (aOR=0.57, 95% CI=0.33–0.99, p-value=0.049) and 80% (aOR= 0.20, 95% CI=0.06–0.65, p-value=0.008), respectively, as compared to being single. Education level of the mother was found to be significantly associated with food insecurity; the lower the educational level, the higher the food insecurity (No education: aOR=16.18, 95% CI=6.08–43.10, p-value<0.001, Primary: aOR=3.52, 95% CI=1.32–9.34, p-value=0.012). The Urdu-speaking population was the most food-secure group in the univariate model; however, after adjusting other factors, there was no significant difference between any ethnic groups (Table 5).

Table 5

Association of Household’s Diet Diversity on Household’s Food Insecurity

	Household’s Food Insecurity
	Model 1		Model 2
	cOR (95% CI)	p-value	aOR (95% CI)	p-value
Household’s diet diversity	0.88 (0.78–0.99)	0.042	0.78 (0.66–0.92)	0.002
Age groups (n=358)
≤35 years	Ref.		–
>35 years	0.79 (0.51–1.22)	0.296
Marital status
Single	Ref.		Ref.
Married	1.30 (0.88–1.92)	0.187	0.57 (0.33–0.99)	0.049
Divorced/Separated	0.91 (0.43–1.93)	0.812	0.20 (0.06–0.65)	0.008
Education level
Masters	Ref.		Ref.
Graduation	0.98 (0.48–2.02)	0.974	1.14 (0.44–2.86)	0.837
Secondary/Intermediate	1.66 (0.85–3.240)	0.138	2.00 (0.83–4.82)	0.120
Primary	3.16 (1.46–6.84)	0.003	3.52 (1.32–9.34)	0.012
No education	10.82 (5.07–23.06)	<0.001	16.18 (6.08–43.10)	<0.001
Ethnicity
Urdu Speakers	Ref.		Ref.
Sindhi	2.82 (1.58–5.02)	<0.001	1.92 (0.91–3.94)	0.083
Punjabi	2.05 (1.13–3.69)	0.017	1.85 (0.92–4.14)	0.114
Balochi	2.42 (1.34–4.36)	0.003	1.30 (0.59–2.77)	0.504
Pashto	1.76 (0.97–3.21)	0.062	1.34 (0.55–2.61)	0.472
Employment status
Unemployed	Ref.		Ref.
Employed	2.52 (1.49–3.38)	<0.001	2.44 (1.39–4.28)	0.002
Self-employed	0.86 (0.47–1.55)	0.623	1.29 (0.58–2.83)	0.524
Body Mass Index (n=371)
Underweight	Ref.		–
Normal	0.70 (0.23–2.10)	0.527
Overweight	0.55 (0.17–1.75)	0.318
Obese	0.55 (0.19–1.58)	0.268
Number of family members (n=408)	1.08 (1.00–1.16)	0.045	1.04 (0.95–1.13)	0.464

Notes: Model 1: Univariate analysis of HH food insecurity with each study variable. Model 2: Multivariate analysis adjusted for variables with p-values<0.250 in univariate analysis variable (Household’s diet diversity, marital status, education level, ethnicity, employment status, and number of family members).

Abbreviations: cOR, crude odds ratio; aOR, adjusted odds ratio; CI, confidence intervals.

Discussion

In this cross-sectional analysis, we have looked at the associations between food insecurity and dietary diversity across five ethnic groups in Pakistan. A positive association was seen between high food insecurity and low dietary diversity among major ethnic populations in Pakistan. The findings showed that households in all ethnic groups lacked dietary diversity. To our knowledge, this is the first investigation reporting a relationship between food insecurity and dietary diversity among different ethnic populations in Pakistan.

The majority of the households consumed less than three groups of food an average, being mainly dependent on carbohydrate-rich food items including wheat, rice, vegetables, and tubers. Similar findings were observed in other low-to-middle-income countries such as Ethiopia, Nigeria, and Tanzania.28–30 Pakistan is facing a vicious circle of the economic crisis resulting in a recent currency devaluation and the hike in prices of consumable goods including food items.31 Affordability is an important pillar of dietary diversity, thus compromise on diversity over quantity is comprehensible.

The results of this study showed that the Urdu-speaking communities are generally more food secure than other ethnicities. Given the fact that they are native residents of the city for decades and Balochi and Sindhi speakers have generally migrated from other cities and villages, a plausible reason could be that the families of individuals who have migrated from villages or small towns to big urban cities even after decades of residence may experience transitional or persistent food crises.32 That may be contributed to early non-employment status, switching to high accommodation prices in cities, and receiving loans.33 Haung et al34 showed that the children of the population who have moved from other areas were less food secure than their native peers, even after a long duration of residence and adjusting for some monetary variables. Another interesting and significant finding of this study was that the households with employed women were 2.5-times more food insecure than those households where the women were not employed. In the local context of Pakistan, the majority of women are housewives (usually men are bread earners) and women work outside of the home in cases of financial constraints such as low income or unemployment of husbands or fathers; however, this may not be true in all cases.35,36

It was found from our result that the educational level of the mother or food preparer endures an effect on food insecurity, the lower the educational level of the person who plans the meal, the higher the food insecurity. This may be secondary to an increase in nutritional education, better understanding of safe and hygienic practices by education.37,38 Numerous studies have found that lack of dietary diversity and scant utilization of locally produced nutritious food may exist across all wealth groups, highlighting the importance of nutrition literacy in addition to food affordability and accessibility in achieving household food security.39,40

Several limitations were identified in our study including concerns related to the study approach, sampling and cultural considerations. The cross-sectional study design provided a descriptive analysis of the study population at a given point in time, concerning several specific behavioral outcomes. However, due to design limitations, it was not possible to establish a clear temporal association between less dietary diversity and economic crisis, and increasing inflation faced by the Pakistani population. Whether less diversity occurred before the crisis or had a causal effect. Another concern is whether all the ethnicities carry the same cultural roots, beliefs, and practices that we assumed may have had influenced their diet or altered behaviors after moving to the urban center.

Recall bias is always a limiting factor, with issues related to accuracy or completeness of recollection of nutritional practices. However, communication limitations were minimized by data collectors who were equipped to be culturally sensitive and speak the same language of the selected participants.

In the ethnically diverse population of Pakistan our results have shown that there are clear ethnic differences. The most food-secure cultural group was those whose ancestors migrated from India and were called Urdu speakers, as compared to Sindhis who have been present in the Indus Valley (now in the province of Sind) for over 5,000 years. More research should be done to explore potential linguistic, religious, and sociopolitical variables that may strongly impact households’ diet in different cultures through a mixed methodology in both urban and rural areas. This may open up the avenues to new “culturally sensitive, cost-effective interventions” and a “culturally sensitive food index” for households that have high acceptability in all cultures.

Conclusion

This study highlighted the quantum of food insecurity and less dietary diversity in all ethnic groups, especially in the Sindhi and Balochi-speaking populations. Ethnic subgroup differences in susceptibility to food security and the reasons for such differences within the Pakistani population have provided additional information about lack of food security reasons when other factors like age, education, and income are similar. This information can potentially be of use by other low-income countries where ethnic differences exist, and food insecurity issues persist. Our results have given considerable hope for improvement in dietary behavior with sound communication and education strategies recognizing ethnic variations in the Pakistani population. However, these attributes are characteristic of many developing countries, to which our findings may be generalizable. We recommend that unmeasured ethnic variations in the food insecurity and dietary diversity deserve further study. There are cultural differences among food consumption patterns in households. The gravity of the nutritional situation in Pakistan mandates culturally sensitive approaches to tackle this mushrooming crisis.