Characterizing the Urban Diet: Development of an urbanized diet index

doi:10.21203/rs.3.rs-929769/v1

Download PDF

Research

Characterizing the Urban Diet: Development of an urbanized diet index

https://doi.org/10.21203/rs.3.rs-929769/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Background

In recent decades China has experienced rapid urbanization leading to a major nutrition transition, with increased refined carbohydrates, added sweeteners, edible oils, and animal-source foods, and reduced legumes, vegetables, and fruits. These changes have accompanied increased prevalence of cardiometabolic disease (CMD). There is no single dietary measure that summarizes the distinct food changes across regions and levels of urbanization.

Objective

We aimed to develop a diet urbanization index that would represent a dietary pattern for urbanized diet across diverse regions.

Methods

Using a sample of adults (≥18 years) in the 2015 wave of the China Health and Nutrition Survey (CHNS; n = 14,024), we selected literature-based candidate dietary variables and tested their univariate associations with overall and within-region urbanization. Using iterative exclusion of select diet-related variables, we created six potential urbanized diet indices, which we examined relative to overall urbanization to select a final urbanized diet index based on a priori considerations, strength of association with urbanization, and minimal missingness. We tested stability of the final urbanized diet index across sociodemographic factors and used mixed effects logistic regression models to examine associations between the final urbanized diet index and CMD risk factors – hypertension (HTN), overweight, and type 2 diabetes mellitus (T2DM), adjusting for sociodemographics, overall urbanization, physical activity, and including random intercepts to account for correlation at community and household level.

Results

We identified a final urbanized diet index that captured dietary information unique to consumption of an urbanized diet and performed well across regions. We found a positive association (R² = 0.17, 0.01 SE) between the final urbanized diet index and overall urbanization in the fully adjusted model. We found no associations between the urbanized diet index with HTN, overweight, and T2DM in fully adjusted models (OR (95% CI) = 0.94 (0.89-1.0), 1.0 (0.95 - 1.05), 1.02 (0.94 - 1.10), respectively).

Conclusions

We derived an urbanized diet index that captured dietary urbanization that was distinct from overall urbanization and performed well across all regions of China. This urbanized diet index provides an alternative to measures of traditional versus Westernized diet that vary across regions due to different cultural dietary traditions.

Health Policy

Environmental Policy

China

diet

urbanization

cardiometabolic disease

Across the globe, the nutrition transition has shifted food consumption away from traditional diets high in vegetables, carbohydrates, fiber and whole foods, and low in fats to a diet comprising processed foods, saturated fats, added sugar, and little fiber (1). Through the process of urbanization, adoption of a “Western diet” is associated with a range of nutrition-related chronic diseases, such as obesity, diabetes, and cardiovascular diseases (2). However, characterizing a “Western diet” within and across diverse countries is challenging due to cultural differences in diet across diverse regions. For example, there are clear differences in urbanized diets in rice-based versus wheat-based culinary traditions across China (3).

Thus, there is a great need to develop urbanized diet measures that generalize across regions and cultural traditions. In addition, lower quality nutrition outcomes are more common in urban environments, particularly among low income residents who are sensitive to economic changes that impact access to healthy foods, yet, urban diets also tend to be more diverse, which is generally better, due to access to a wider variety of foods than in rural areas (4). However, this advantage is not equally distributed across socioeconomic strata, and this urban diet paradox presents another need to better understand the urbanized diet(4).

While measures of overall urbanization in China capture the multidimensionality of urbanization and move beyond the discrete urban-rural dichotomy (5), such measures do not specifically focus on dietary changes. Given the rapid and differential urbanization changes in dimensions of society with urbanization, there is considerable variation in diet across the full spectrum of urbanization (6). Yet, there is no single measure of urbanized diet to quantify changes in diet and diet-related behaviors that have occurred in recent years in China, a country with substantial heterogeneity in diet by geographic region (7) (8). Given the variation in diet across regions and dietary cultural traditions, the lack of a single urbanized diet measure leaves fundamental gaps in understanding how the dietary urbanization associates with diet-related chronic diseases over time.

China provides a unique opportunity to study changes in diet and its association with cardiometabolic disease (CMD) because of the rapid and differential pace of urbanization that has occurred in recent decades across its regions. With the relatively recent political and social reforms, China has transitioned from a country typically burdened with undernutrition to one with a rapid increase in obesity (9). Other non-communicable diseases, such as type 2 diabetes mellitus (T2DM), have become more prevalent and are now leading causes of morbidity, disability, and mortality in China (9). As a whole, diet-related non-communicable diseases have become common, which is of great concern (10). The China Health and Nutrition Survey (CHNS) is an extensive ongoing, prospective cohort study that began in 1989 with the intent to capture the health impacts of urbanization. The availability of large scale, high-quality CHNS dietary data provides an outstanding exemplar population in which to develop a measure of diet urbanization and investigate associations between diet and CMD.

We aimed to develop an urbanized diet index reflecting the total diet that would generalize across different geographic regions and across varied culinary traditions and consumption patterns. To this end, we sought to generate a single measure of urbanized diet using a wide variety of diet consumption and diet-related infrastructure variables. We capitalized upon the rapid and differential urbanization across China and the deeply characterized diet data of the CHNS to create an urbanized diet index. We then tested the association between this urbanized diet index with three key CMD risk factors – hypertension (HTN), overweight, and T2DM.

Study Participants

The China Health and Nutrition Survey (CHNS) is an ongoing household-based, cohort study. It began in 1989, with surveys completed every 2–4 years until 2015. Prior to 2011, the sample included data from nine provinces (Heilongjiang, Liaoning, Jiangsu, Shandong, Henan, Hubei, Hunan, Guangxi, Guizhou). In 2011 three megacities were added (Beijing, Shanghai and Chonqing) and in 2015 three additional provinces were added (Shaanxi, Yunnan, and Zhejiang). A stratified, multistage, clustered sampling design was used to select households and communities within each province or mega-city. The CHNS captures a variety of geographical areas, levels of economic development, and health indicators (11). Data were collected through a 7-day survey including collection of biomarker data, 3-day 24-hour diet recall data, eating behaviors, financial resources, and other relevant measures, making the CHNS a robust resource for studying diet urbanization and cardiometabolic disease. The study met the standards for the ethical treatment of participants and was approved by the Institutional Review Boards of the University of North Carolina at Chapel Hill and the National Institute for Nutrition and Health, Chinese Center for Disease Control and Prevention. Participants gave informed consent for participation (11).

To develop the urbanized diet index, we used data from all 12 provinces and three megacities in the 2015 wave of the CHNS for all adults, 18 years or older, who had diet data collected in 2015 (n = 17,191). All inclusion and exclusion criteria are shown in Fig. 1. We excluded participants who were pregnant, had implausible (< 500 kcals) or missing average daily energy intake (kcal), or had missing covariate or health outcome data.

Diet Data

Dietary intake was measured with 24-hour dietary recalls on three consecutive days. Household food consumption was determined by examining changes in inventory from the beginning to the end of each day, in combination with a weighing and measuring technique. Three-day inventory and dietary recall interviews were randomized to occur between Monday-Sunday and were conducted by trained investigators. Interviewers recorded the type of food, amount consumed, timing of consumption and location of foods consumed, which included food consumed outside the home, for each individual participant. Side dishes, condiments and spices were recorded and weighed at the beginning and end of the food inventory and allocated to household members based on individual 24-hour recall data. Food consumption data was converted to nutritional content using the Chinese Food Composition Table. Other variables were derived from the 24-hour recall data, including percent of daily calories consumed from a variety of foods or food groups (fruit, nuts and seeds, all snack foods, sweet snack foods, eggs, dairy, fried foods, fast food, instant noodles, high fat meat, carbohydrates, fat, animal products, processed foods), and daily averages of number of snacks eaten, number of food groups eaten, sodium intake, and fiber intake.

In addition, we used diet-related data from other parts of the CHNS questionnaire. For example, we used self-reported data on the amount of each type of alcohol consumed each week to generate an indicator for whether an individual consumed any wine. We used data from the household questionnaire to generate diet-related infrastructure indicators of urbanization, including household ownership of refrigerator and/or microwave.

Health Outcome Variables

Hypertension (HTN)

Blood pressure was measured three times on the right arm after 10 minutes of seated rest by physicians using standard mercury sphygmomanometers (12). A participant with systolic blood pressure ≥ 130 or diastolic blood pressure ≥ 80 was classified as having HTN based on the ACC/AHA blood pressure guidelines (13). Additionally, participants who reported a HTN diagnosis or who reported they were currently taking HTN medication were also defined as having hypertension.

Overweight

Height was measured without shoes to the nearest 0.1 cm using a portable stadiometer and weight was measured without shoes and in light clothing to the nearest 0.1 kg on a calibrated floor scale. Both were measured by trained anthropometrists and were used to calculate body mass index (BMI). Overweight was defined as having a BMI of 24 kg/m² or greater, based on the Chinese overweight BMI cut point (14).

Type II Diabetes Mellitus (T2DM)

Fasting blood glucose was measured through blood sample collection, after an overnight fast, by health workers, according to standard procedures (15). T2DM was determined based on individuals having either a fasting blood glucose level of 126 mg/dL or greater (16), the International Diabetes Federations (IDF) criteria for T2DM diagnosis (IDF, 2017), or a self-reported previous T2DM diagnosis or taking T2DM medications.

Urbanization Index

To define dietary consumption patterns relative to overall urbanization, we used the urbanization index, a validated multicomponent measure of urbanization in the CHNS that captures rapid and differential urbanization across China (5). This measure of overall urbanization, calculated at the community level, was derived using a multicomponent continuous scale based on very detailed direct measurement of community contextual measures. The 12 components that comprise this validated index include population density, economic activity, traditional markets, modern markets, transportation infrastructure, sanitation, communications, housing, education, diversity, health infrastructure and social services.

Covariates

Age, sex, educational level (highest attained), and smoking history were self-reported. Due to the high correlation between sex and smoking history, we derived a combined variable of sex and smoking status (female never smoker, male never smoker, female ever smoker, male former smoker, male current smoker) for use as a covariate. Region was categorized as North, Central, South, and Megacities. Per capita household income, in Yuan, was derived from individual and household questionnaires from time-use, asset, and economic activity. Average total daily energy intake, in calories, was collected across the repeated 24-hour recall data, which was validated using doubly labelled water. Physical activity was measured using a detailed weekly activity recall which captured occupational, domestic, travel and active leisure activity. Metabolic equivalent of task (MET) hours per week were calculated for each of these categories of physical activity using the Compendium of Physical Activity (18–20). A total physical activity, capturing average MET hours per week, was calculated as the sum across these four categories. More detail on these physical activity variables, including how the METs were calculated can be found elsewhere (21) (22).

Urbanized Diet Index Development

Select Dietary Variables of Interest

First, we selected a set of dietary variables with literature-based evidence for association with urbanization. Second, we examined consistency in the association between these variables with urbanization index, overall and within each region (North, Central, South) and within Megacities, by examining the mean and standard deviation for each continuous dietary measure (e.g., mg of sodium) and percentage for each dichotomous dietary measure (e.g., whether or not an individual consumed wine) by tertile of overall and within-region urbanization index. Third, we determined which variables to move forward, based on (1) association with overall urbanization, (2) consistent association with urbanization across regions, and (3) frequency of consumption > 5%. We used substantive differences in means or percentages across urbanization and region to determine inclusion/exclusion rather than formal statistical testing.

Score Individual Dietary Variables

We categorized each individual diet variable for scoring. For uncommonly consumed foods (< 80% consumers), we created a categorical variable for non-consumers and among consumers, quartiles of consumption. For commonly consumed foods (≥ 80% consumers), we created quintiles of consumption with non-consumers classified in the lowest quintile. Non-continuous variables were scored as dichotomous yes/no consumed or ownership.

To determine scoring of diet variables for the urbanized diet index, we examined the association between each of the candidate diet variables, as categorized above, with the overall urbanization index using logistic or multinomial logistic regressions in three sets of models: among commonly consumed foods (≥ 80% consumers); uncommonly consumed foods (< 80% consumers), and dichotomous ownership or consumption variables. For the non-dichotomous dietary variables, we present results as relative risk ratios (RRRs) for associations between each dietary variable and a one standard deviation change in urbanization index, using a referent of the lowest consumer or non-consumer group. For dichotomous variables, we present odds ratios (ORs) for associations between each variable and a one-unit change in overall urbanization index. Using these results, we scored the variables based on the strength and direction of their association with overall urbanization.

Analytical steps to create the urbanized diet index

First, we created six candidate urbanized diet indices for consideration based on inclusion and exclusion of specific food variables. Second, we tested associations for each of the six candidate diet urbanization indices with overall urbanization in unadjusted, age- and sex-adjusted, and fully adjusted (age, sex and smoking, average daily energy intake, region, educational attainment, per capita household income, and physical activity) mixed linear regression models, with random intercepts to account for correlation at the household and community level. We used these model results to select a final diet urbanization measure based on the strength of association with overall urbanization, and degree of missingness.

We then tested whether the final diet urbanization index was stable across varying sociodemographic characteristics using standardized residuals greater or less than two. Using a fully adjusted mixed model linear regression, we classified individuals into two groups, one with individuals who had less accurate prediction of final urbanized diet index (standardized residuals >|2|), and one with individuals who had more accurate prediction of final urbanized diet index (standardized residuals ≤ |2|). We generated average demographic data for each of these two groups and conducted ANOVA and chi-squared tests to determine statistical significance of differences across overall urbanization index, age, sex and smoking status, region, average daily energy intake, educational attainment, per capita household income, and physical activity.

Finally, we tested whether the final diet urbanization index was associated with three key CMD risk variables - HTN, overweight, and T2DM. We used mixed effects logistic regression models with urbanized diet index as the exposure, adjusted for age, sex/smoking, average daily energy intake, region, educational attainment, per capita household income, physical activity, and random intercepts to account for correlation at community and household levels. In a second set of models, we included all the previous covariates with the addition of the overall urbanization index (5), to compare associations with and without controlling for overall urbanization. We present these results, for a standardized diet urbanization index, as odds ratios for the association between each of the three CMD outcomes and a one standard deviation change in urbanized diet index.

Our sample had a slightly higher proportion of women than men (Table 1). Most participants had completed at least some post-primary schooling and lived in the Central or Southern regions of China.

Table 1

Descriptive statistics of sample (N = 14,037) for generation of Chinese urbanized diet indices, CHNS, 2015
Age (years) – mean (std)	51.45 (15.12)
Female – N (%)	8,050 (57.40)
Urbanization Index^*– mean (std)	72.52 (17.73)
Educational Attainment – N (%)
Completed no School	1,045 (7.45)
Completed some or all of Primary School	3,360 (23.96)
Completed some Post-Primary School	7,583 (54.07)
Completed College	2,036 (14.52)
Income ^† – mean (std)	24.19 (36.30)
Region – N (%)
North	1,587 (11.32)
Central	4,645 (33.12)
South	4,939 (35.22)
Megacities	2,853 (20.34)
^*Overall urbanization index is a validated multicomponent measure of urbanization in the CHNS (Jones-Smith & Popkin, 2010).
^†Income is gross per capita household income, in thousands of Yuan.

Select Dietary Variables of Interest

We identified 23 dietary variables (Table 2), as candidates for inclusion in the urbanized diet index and conducted preliminary analyses for each to determine patterns of association with overall urbanization index and region. Based on this analysis, we excluded alcohol consumption, percent of calories from fat, sodium intake, and fiber intake, variables which did not vary with overall urbanization. We excluded fast food consumption as it was infrequently consumed (< 5%). We excluded instant noodle consumption, which did not vary with urbanization and was infrequently consumed. Ultimately 17 variables were carried forward that were consistently associated with overall urbanization and region (Supplemental Table 1).

Table 2

Dietary responses across tertiles of overall urbanization index*, CHNS, 2015
	Urbanization Tertiles
	Low 29.2–61.2	Moderate 61.3–82.5	High 82.6–104.4	--
Urbanization Index – mean (SD)	50.7 (7.2)	72.9 (6.3)	91.2 (5.2)	--
^†Drink alcohol - % yes	27.3	27.9	27.9	Excluded - No association with urbanization
^†% kcal from fat – mean (SD)	32.4 (12.7)	35.4 (11.4)	35.1 (10.9)	Excluded - No association with urbanization
^†Sodium (mg) - mean (SD)	3896 (2264)	3974 (2221)	3780 (2092)	Excluded - No association with urbanization
^†Fiber (g) - mean (SD)	10.3 (7.2)	10.8 (7.3)	10.8 (6.9)	Excluded - No association with urbanization
^†Eat fast food – % yes	2.1	2.4	3.2	Excluded - Low consumption
^†Eat instant noodles - % yes	2.0	1.9	2.2	Excluded-No association with urbanization, low consumption
Drink wine - % yes	8.4	22.4	31.6	Included - Positive association with urbanization
Own fridge - % yes	85.2	94.0	96.9	Included - Positive association with urbanization
Own microwave - % yes	9.6	39.8	63.0	Included - Positive association with urbanization
Eat fruit - % yes	28.7	40.0	52.25	Included - Positive association with urbanization
Eat nuts/seeds - % yes	11.2	19.4	21.9	Included - Positive association with urbanization
Eat all snack foods - % yes	11.2	17.9	21.5	Included - Positive association with urbanization
Eat sweet snacks - % yes	6.4	12.5	15.5	Included - Positive association with urbanization
Eat eggs - % yes	55.4	70.8	77.2	Included - Positive association with urbanization
Eat dairy - % yes	9.42	19.1	35.8	Included - Positive association with urbanization
Eat away from home - % yes	15.7	30.2	32.0	Included - Positive association with urbanization
Eat fried foods - % yes	18.3	28.5	31.3	Included - Positive association with urbanization
Eat high fat meat - % yes	72.7	85.9	87.0	Included - Positive association with urbanization
Eat animal-source foods - % yes	77.7	90.3	92.6	Included - Positive association with urbanization
Eat processed foods - % yes	93.7	96.3	98.1	Included - Positive association with urbanization
# of snacks - mean (SD)	0.3 (0.5)	0.4 (0.6)	0.6 (0.8)	Included – Positive association with urbanization
# of food groups - mean (SD)	9.9 (2.9)	11.5 (3.2)	12.6 (3.4)	Included - Positive association with urbanization
% kcal from carbs – mean (SD)	55.6 (12.7)	51.1 (11.6)	50.7 (11.0)	Included – Negative association with urbanization
^* Overall urbanization index is a validated multicomponent measure of urbanization in the CHNS (Jones-Smith & Popkin, 2010).
^†Variables not included in the urbanized diet index development process.
23 variables considered for inclusion in urbanized diet index – drinking alcohol, % of calories from fat, sodium intake and fiber intake excluded for lack of variation with overall urbanization; eating fast food excluded for low consumption (< 5%); eating instant noodles excluded for lack of variation with urbanization and low consumption, carrying 17 variables forward for index development.

Score Individual Dietary Variables

Using logistic or multinomial logistic regressions, we examined the association between each of the 17 candidate variables in models based on frequency of consumption and/or ownership. (Fig. 2). Among the frequently consumed foods (Fig. 2A), most were positively associated with overall urbanization, with strongest association for number of food groups, and a negative association between percent of calories from carbohydrates with overall urbanization. Among the less frequently consumed foods (Fig. 2B), we found comparatively stronger positive associations for percent of calories from dairy, fruit and number of snacks with urbanization. All three dichotomous dietary variables (Fig. 2C) were positively associated with urbanization.

We assigned scores to each of the diet-related variables (Table 3), based on the magnitude and pattern of association with overall urbanization (Fig. 2). Scoring was based on consistency and direction of associations (Table 3).

Table 3

Scoring of variables based on magnitude and pattern of association with overall urbanization*†, CHNS, 2015
Commonly consumed foods (≥ 80% consumption)
	Quintile 1^‡	Quintile 2	Quintile 3	Quintile 4	Quintile 5
Number of Food Groups	0	1	2	3	4
% kcals from Animal Sources	0	1	2	3	2
% kcals from High Fat Meat	0	1	2	3	2
% kcals from Processed Foods	0	1	1	2	2
% kcals from Carbohydrates	4	3	2	1	0
Uncommonly consumed foods (< 80% consumption)
	Non-consumers	Quartile 1	Quartile 2	Quartile 3	Quartile 4
% kcals from Dairy Products	0	1	2	3	4
% kcals from Fruit	0	1	2	3	4
% kcals from Fried Food	0	1	2	3	2
% kcals from Eggs	0	1	1	2	2
Number of Snacks	0	1	2	3	4
Dichotomous dietary variables
	Non-consumers or non-owners	Consumers or owners	--	--	--
% kcals from Nuts & Seeds	0	2	--	--	--
% kcals from Away-from-Home Eating	0	2	--	--	--
% kcals from Sweet Snacks	0	2	--	--	--
% kcals from All Snack Foods	0	2	--	--	--
Own a microwave	0	1	--	--	--
Own a refrigerator	0	1	--	--	--
Drink wine	0	1	--	--	--
^* Using logistic or multinomial logistic regressions shown in Fig. 2.
^†Overall urbanization index is a validated multicomponent measure of urbanization in the CHNS (Jones-Smith & Popkin, 2010).
^‡Including non-consumers.

Create Candidate Urbanized Diet Indices

Each of the 17 candidate variables were carried forward to derive a set of candidate urbanized diet indices (Table 4). We excluded sets of variables for six candidate diet urbanization indices and tested each of these indices to derive a final urbanized diet index. Candidate indices were considered for exclusion if they included variables with considerable missingness or if they had a weaker magnitude of association with urbanization (e.g., RRRs and ORs close to one). Candidate indices were also excluded if they captured duplicate diet features (e.g., numbers of snacks and % of calories from snack foods).

Table 4

Urbanized diet index component variables, CHNS, 2015
	Urbanized Diet Index versions
Dietary Variables	1	2	3	4	5	6
Drinks Wine – yes/no	X
Owns a Refrigerator – yes/no	X	X
Owns a Microwave – yes/no	X	X
% of calories consumed from All Snack Foods	X	X	X		X
% of calories consumed from Sweet Snacks	X	X	X	X
Number of Snacks consumed	X	X	X	X	X	X
% of calories consumed from Fruit	X	X	X	X	X	X
% of calories consumed from Nuts & Seeds	X	X	X	X	X	X
% of calories consumed from Eggs	X	X	X	X	X	X
% of calories consumed from Dairy Products	X	X	X	X	X	X
% of calories consumed from Fried Foods	X	X	X	X	X	X
% of calories consumed from Processed Foods	X	X	X	X	X	X
% of calories consumed from High Fat Meat	X	X	X	X	X	X
% of calories consumed from Animal Source Food	X	X	X	X	X	X
% of calories consumed from Away-from-Home Eating	X	X	X	X	X	X
% of calories consumed from Carbohydrates	X	X	X	X	X	X
Number of Food Groups consumed	X	X	X	X	X	X
Percentages and number of foods consumed were determined based on 3-day averages.

Diet Index 1 included all variables, Diet Index 2 excluded wine consumption based on its high missingness, and Diet Index 3, in addition to excluding wine consumption as was done in Diet Index 2, dropped refrigerator and microwave ownership due to low magnitude of association with overall urbanization (Fig. 2). In Diet Indices 4, 5 and 6 we iteratively excluded percent of calories from sweet snack foods, all snack foods, and both snack variables, due to similarity in concept of snacking behavior, from Diet Index 3, which excluded wine consumption, refrigerator, and microwave ownership.

Select a Final Urbanized Diet Index

We compared the associations between each of the six candidate indices (Table 5) with overall urbanization, finding weak positive associations for each. Although the exclusion of the wine consumption variable (Urbanized Diet Index 2) made little difference in the R² relative to Urbanized Diet Index 1, its greater missingness (72%) was less desirable. Each sequential model resulted in minimal difference in R². Based on the strongest correlation with overall urbanization in the fully adjusted model (R² = 0.17, 0.01 SE) and low missingness (0.31%), we selected Urbanized Diet Index 2 as the final urbanized diet index.

Table 5

Associations between each candidate urbanized diet index and overall urbanization index*, CHNS, 2015

Diet Index	Description	N	Mean (SD)	Min	Max	% Missing	Model I^†	Model II ^‡	Model III^§
1	Wine consumption, owning a refrigerator, owning a microwave, daily average % of calories consumed from foods or food groups (fruit, nuts and seeds, all snack foods, sweet snack foods, eggs, dairy products, fried foods, away-from-home eating, high fat meat, carbohydrates, animal-source foods, and processed foods), daily average number of snacks consumed and food groups consumed.	3,866	15.76 (6.76)	0	38	72.40	0.18 (0.01)	0.17 (0.01)	0.16 (0.01)
2	Diet index 1, excluding wine consumption	13,981	15.36 (6.97)	0	39	0.31	0.18 (0.01)	0.18 (0.01)	0.17 (0.01)
3	Diet index 2, excluding owning a refrigerator and microwave	14,024	14.04 (6.72)	0	37	0.0	0.17 (0.01)	0.17 (0.01)	0.16 (0.01)
4	Diet index 3, excluding % of calories from all snack foods	14,024	13.72 (6.43)	0	35	0.0	0.16 (0.01)	0.16 (0.01)	0.15 (0.01)
5	Diet index 3, excluding % of calories from sweet snack foods	14,024	13.82 (6.49)	0	35	0.0	0.16 (0.01)	0.16 (0.01)	0.15 (0.01)
6	Diet index 3, excluding % of calories from all snack foods and sweet snack foods	14,024	13.50 (6.25)	0	33	0.0	0.16 (0.01)	0.16 (0.01)	0.15 (0.01)
Linear mixed model for association of candidate urbanized diet indices and overall urbanization index, accounting for correlations at community and household levels.
^*Urbanization index is a validated multicomponent measure of urbanization in the CHNS (Jones-Smith & Popkin, 2010).
^† Model I, unadjusted.
^‡ Model II, minimally adjusted for age and sex.
^§ Model III, fully adjusted for age, sex/ smoking status, average daily energy intake (kcals), region, education level, per capita gross household income, physical activity.

Test the Stability of the Final Urbanized Diet Index across Sociodemographic Factors

Using fully adjusted mixed model linear regressions, we tested stability of Urbanized Diet Index 2 across key sociodemographic factors with contrasts for standardized residuals ≤/>|2| using ANOVA and chi-squared tests. We found statistically significant differences in the accuracy of the prediction of urbanized diet index by overall urbanization, age, income, energy intake, region, and educational attainment (Table 6). Individuals with poorer accuracy in prediction of the urbanized diet index had on average a higher overall urbanization index, were younger, had a higher income, and had a lower caloric intake. Similarly, those with poorer accuracy were more likely to be from the Megacities or Central China, as compared to the North and South regions of China, and were more likely to have at least a college education.

Table 6

Differences in observed versus predicted* overall urbanization by urbanized diet, CHNS, 2015
	Standardized residuals ≤ \|2\|	Standardized residuals > \|2\|
N	13,626	355
^‡Overall Urbanization Index^† – mean (SD)	72.4 (17.8)	77.9 (15.1)
^‡Age (years) – mean (SD)	51.6 (15.1)	48.1 (16.4)
^‡Income^§ – mean (SD)	24.1 (36.3)	28.3 (37.2)
^‡Energy Intake (kcals) – mean (SD)	1844.1 (623.7)	1766.0 (607.6)
^‡Region – N(%)
North	1,559 (11.4)	24 (6.8)
Central	4,504 (33.1)	132 (37.2)
South	4,824 (35.4)	96 (27.0)
Megacities	2,739 (20.1)	103 (29.0)
^‡Educational Attainment – N (%)
Completed no School	1,016 (7.5)	23 (6.5)
Completed some or all of Primary School	3,286 (24.1)	64 (18.0)
Completed some Post-Primary School	7,385 (54.2)	179 (50.4)
Completed College	1,939 (14.2)	89 (25.1)
Sex and Smoking Status – N (%)
Female Never Smoker	7,475 (54.9)	204 (57.5)
Male Never Smoker	1,964 (14.4)	42 (11.8)
Female Ever Smoker	339 (2.5)	7 (2.0)
Male Former Smoker	431 (3.2)	15 (4.2)
Male Current Smoker	3,417 (25.1)	87 (24.5)
Physical Activity (METs) – mean (SD)	162.4 (172.3)	154.9 (140.0)
*Standardized residuals ≤ \|2\| versus > \|2\| from adjusted linear mixed model predicting overall urbanization index from urbanized diet index, controlling for age, sex and smoking status, daily energy intake, region, educational attainment, household income, daily physical activity, and accounting for community and household correlations.
^†Overall urbanization index is a validated multicomponent measure of urbanization in the CHNS (Jones-Smith & Popkin, 2010).
^‡ Statistically significant differences between the two groups.
^§Income is gross per capita household income, in thousands of Yuan.

Test the Association between the Final Urbanized Diet Index and CMD risk

We found a statistically significant association between Urbanized Diet Index 2 and T2DM in the fully adjusted model, which was attenuated with the addition of overall urbanization (Fig. 3). We found no statistically significant association between Urbanized Diet Index 2 and HTN or overweight in either model.

Using 2015 CHNS data we developed a measure that captures dietary urbanization and is positively associated with overall urbanization. Our analysis of associations between each of the individual dietary variables of interest and overall urbanization effectively informed our decisions about inclusion and exclusion of variables, and appropriate scoring of variables for generation of a set of urbanized diet indices. The set of six slightly different composite urbanized diet indices allowed for comparison of each measure’s association with overall urbanization. We determined that the urbanized diet index that was the best indicator of diet urbanization included the ownership of a refrigerator and microwave, percent of calories consumed from fruit, nuts and seeds, all snack foods, sweet snacks, eggs, dairy, fried food, away-from-home eating, high fat meat, carbohydrates, animal source foods and processed foods, as well as daily average number of snacks consumed, and daily average number of food groups consumed. The diet index excluded wine consumption. This decision was based on the little impact on the association with overall urbanization upon exclusion of the wine consumption variable, and a large reduction in percent of participants missing data. When adjusting for overall urbanization index, we saw no associations between the urbanized diet index and three CMDs - HTN, overweight and T2DM.

We saw a positive association between the urbanized diet index and overall urbanization index (R² = 0.17 (0.01 SE)). This suggests that the urbanized diet index does not solely represent urbanization and provides distinct information, specifically about diet. We assessed the stability of the final derived urbanized diet index across sociodemographic factors, finding comparatively lower predictive accuracy by overall urbanization, age, region, energy intake, educational attainment, and income, thus suggesting some variation in ability to discriminate according to sociodemographic factors.

Some of the results of our study were unanticipated. We expected that sodium consumption would be an important factor in diet urbanization, but our findings suggest otherwise. In preliminary analysis we saw only slight changes in sodium intake with changing urbanization. This led us to fully exclude the sodium variable from the index development process. This inconsistency may be related to differences in sources of sodium in China based on level of urbanization. While salt was traditionally used as a preservative, this method is being replaced by refrigeration, which becomes more available with increasing urbanization. Yet, high sodium intake persists mainly due to salt added in cooking and increased consumption of high sodium processed foods, especially in higher urbanization areas (23)(24). In addition, reduced fiber intake is considered a factor in the nutrition transition (1). However, in preliminary analyses we saw no trend between increasing urbanization and fiber intake, thus the fiber variable was not included in the index development. This may be due to increased consumption of foods high in fiber, like fruit, to make up for reduced consumption of whole grains in high urbanization communities (25).

Additionally, we saw inconsistent associations between each of the snack-related variables (number of snacks consumed, percent of calories from all snack foods, and percent of calories from sweet snack foods) and overall urbanization. While RRRs for associations with overall urbanization steadily increased with increasing number of snacks consumed, RRRs remained stable for percent of calories consumed from all snack foods and sweet snacks, indicating that the behavior of snacking may be a greater factor in diet urbanization than the proportion of calories consumed from snacks. This could be explained by differences in what data was used to construct these dietary variables. For the three-day dietary recall, individuals were asked to classify all foods they ate by meal, specifically breakfast, lunch, dinner or snack. Number of snacks was based on the number of times an individual reported they ate a snack, or a meal outside of breakfast, lunch or dinner. This snack could be comprised of one or more types of food groups. Earlier work in this cohort has suggested, for example, that fruits and beverages are often reported as snacks (26). However, percent of calories from snacks or sweet snacks were based on food group coding. Therefore, fruits and beverages, which constitute their own food group, would not be included in this category; however, foods like chips and snack cakes are included. This suggests while these measures might be correlated, they are capturing different components of an urbanized diet and therefore could easily result in differences in the association with urbanization.

We expected that HTN, overweight, and T2DM would be associated with urbanized diet, but we saw no association in models adjusted for overall urbanization. Though we found a positive association with T2DM in the fully adjusted model with no control for overall urbanization. It is important to consider that the index we generated does not measure quality of diet, thus a more urbanized diet does not necessarily indicate a less nutritious diet, which would suggest an association with CMD. For example, our index includes nutritious foods like fruit and nuts and seeds, along with adverse foods like fried foods and high fat meat, which can have opposing impacts on cardiometabolic health. However, in their examination using 2006 CHNS data, Wang et. al also found inconsistent associations between two diet quality measures and CMD odds; finding a negative association between the China Dietary Guideline Index and T2DM, but not with the tailored Alternative Eating Index, although both were negatively associated with abdominal obesity(27). The authors found that neither of the diet quality measures were associated with T2DM in women, and in both sexes neither were associated with high blood pressure or metabolic syndrome (27). It is possible that other factors, in addition to diet urbanization and diet quality, play a role in CMD outcomes amid a nutrition transition, such as macronutrient energy balance, the microbiome, or chronic stress response.

Our findings point to a complex dynamic related to urbanization, diet, and the nutrition transition, leaving more work to do to understand how exactly the nutrition transition leads to adverse cardiometabolic health. Better measurement of the urbanized diet may aid in better understanding of this complex dynamic, and our paper is one step in that direction. Traditional diets differ globally and within China, so there is no single set of factors that capture the nutrition transition. Thus, our paper presents a methodology approach that may be broadly applied to different populations, but likely with adaptation and tailoring to local cultural dietary traditions and regional variations.

A main strength of our study is the use of the large CHNS study sample, which includes extensive covariate data allowing us to adjust models for important confounding variables. Three 24-hour diet recalls allowed for the use of average diet data, which increased reliability of our analyses. In addition, the use of household inventories and questionnaires enabled us to use data gathered from multiple sources, further increasing reliability. The diet index includes variables capturing many aspects of an urbanized diet.

There were some limitations to our study and its results. Our findings are limited to the Chinese population and cannot be extrapolated to other populations, as only CHNS data was used. There is also great diversity within China, so we may not have accurately captured an urban diet for all people in China. In addition, the study only utilized 2015 data, thus preventing longitudinal analysis of changes in diet and health outcomes. As is true with all self-reported data, diet recalls may not be fully representative of each participant’s diet, and even three days of diet assessment may not be enough to capture some episodically consumed foods. Another minor limitation within the CHNS is that the CHNS survey questions did not specifically delineate between type I and type II DM. However, it can be presumed that most, if not all, cases of DM are type II, as type I DM has very low prevalence in China. As mentioned previously, the measure we developed captured the urbanized diet but does not provide insight into diet quality specifically, so further investigation should be done into the quality of an urbanized diet and associations with disease.

The methods used for this project, and the urbanized diet index itself, add to the resources available to study and better understand dietary urbanization. While there are existing measures that capture overall urbanization, we developed a measure that specifically captures dietary urbanization. The methodology presented can be used to further study changes in diet and its impacts in other urbanizing countries.

In our cross-sectional study of adults in the 2015 CHNS we developed a measure of diet urbanization that was positively associated with overall urbanization and provides information specifically about diet. However, our urbanized diet index was not associated with HTN, overweight, and T2DM prevalence in China, suggesting variation in diet quality with urbanization. The derived diet urbanization index can be used to expand upon our understanding of patterns of dietary urbanization and its relation to sociodemographic and health factors in China and potentially across low-middle income countries, and over time.

CMD

Cardiometabolic disease

T2DM

Type 2 diabetes mellitus

CHNS

China Health and Nutrition Survey

HTN

Hypertension

BMI

Body mass index

METs

Metabolic equivalent of task

Odds ratio

RRR

Relative risk ratio

Acknowledgements

The authors thank the National Institute for Nutrition and Health, Chinese Center for Disease Control and Prevention, the Carolina Population Center, the UNC-CH and all the people involved in the China Health and Nutrition Survey.

Ethics approval and consent to participate

The study met the standards for the ethical treatment of participants and was approved by the Institutional Review Boards of the University of North Carolina at Chapel Hill and the National Institute for Nutrition and Health, Chinese Center for Disease Control and Prevention. Participants gave informed consent for participation.

Consent for publication

Not applicable.

Availability of data and materials

Data available at: http://www.cpc.unc.edu/projects/china/

Competing interests

The authors declare that they have no competing interests.

Funding

This research uses data from China Health and Nutrition Survey (CHNS). We are grateful to research grant funding from the National Institute for Health (NIH), the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) for R01 HD30880, National Institute on Aging (NIA) for R01 AG065357, National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) for R01DK104371 and R01HL108427, the NIH Fogarty grant D43 TW009077 since 1989, and the China-Japan Friendship Hospital, Ministry of Health for support for CHNS 2009, Chinese National Human Genome Center at Shanghai since 2009, and Beijing Municipal Center for Disease Prevention and Control since 2011. We thank the National Institute for Nutrition and Health, China Center for Disease Control and Prevention, Beijing Municipal Center for Disease Control and Prevention, and the Chinese National Human Genome Center at Shanghai. We are also grateful for funding from the NICHD to Carolina Population Center (CPC) at the UNC-CH (NIH grant P2C HD050924). This project was also supported by the Gillian T. Cell Senior Thesis Research Award in the University of North Carolina College of Arts & Sciences administered by Honors Carolina

Authors' contributions

PGL, AGH and KM contributed to conception and design of the study, ACS and ES conducted data analyses, ACS, AGH, and PGL drafted the manuscript, KM, SD, BZ, and HW contributed to critical revision of the manuscript. All authors have read and approved the final manuscript.

Popkin BM. An overview on the nutrition transition and its health implications: the Bellagio meeting. Public Health Nutr. 2002 Feb;5(1a):93–103.
Popkin BM. Nutrition Transition and the Global Diabetes Epidemic. Curr Diab Rep. 2015 Sep;15(9):64.
Huang L, Wang H, Wang Z, Zhang J, Zhang B, Ding G. Regional Disparities in the Association between Cereal Consumption and Metabolic Syndrome: Results from the China Health and Nutrition Survey. Nutrients. 2019 Apr;11(4).
Vilar-Compte M, Burrola-Mendez S, Lozano-Marrufo A, Ferre-Eguiluz I, Flores D, Gaitan-Rossi P, et al. Urban poverty and nutrition challenges associated with accessibility to a healthy diet: a global systematic literature review. Int J Equity Health. 2021 Jan;20(40).
Jones-Smith JC, Popkin BM. Understanding community context and adult health changes in China: Development of an urbanicity scale. Soc Sci Med. 2010 Oct;71(8):1436–46.
Wang Z, Gordon-Larsen P, Siega-riz AM, Cai J, Wang H, Adair LS, et al. Sociodemographic Disparity in the Diet Quality Transition among Chinese adults, 1991 to 2011. Eur J Clin Nutr. 2017 Apr;71(4):486–93.
Li Y, He Y, Lai J, Zhang J, Fu P. Dietary patterns are associated with stroke in Chinese adults. J Nutr. 2011 Oct;141(10):1834–9.
Zhang X, Dagevos H, He Y, Van der Lans I, Zhai F. Consumption and corpulence in China: A consumer segmentation study based on the food perspective. Food Policy. 2008;33(1):37–47.
Gordon-Larsen P, Wang H, Popkin BM. Overweight dynamics in Chinese children and adults: Chinese obesity dynamics. Obes Rev. 2014 Jan;15:37–48.
Li M, Shi Z. Dietary Pattern during 1991–2011 and Its Association with Cardio Metabolic Risks in Chinese Adults: The China Health and Nutrition Survey. Nutrients. 2017 Nov 6;9(11):1218.
Popkin BM, Du S, Zhai F, Zhang B. Cohort Profile: The China Health and Nutrition Survey -- monitoring and understanding socio-economic and health change in China, 1989–2011. Int J Epidemiol. 2010;39(6):1435–40.
Wang Y, Howard AG, Adair LS, Wang H, Avery CL, Gordon-Larsen P. Waist Circumference Change is Associated with Blood Pressure Change, Independent of BMI Change. Obesity. 2019.
Whelton P, Carey R, Wilbert A, Casey D. 2017 ACC/AHA/AAPA/ABC/ACPM/AGS/APhA/ASH/ASPC/NMA/PCNA Guideline for the Prevention, Detection, Evaluation, and Management of High Blood Pressure in Adults. J Am Coll Cardiol. 2018 May;71(19).
Zhou B-F. Predictive values of body mass index and waist circumference for risk factors of certain related diseases in Chinese adults–study on optimal cut-off points of body mass index and waist circumference in Chinese adults. Biomed Environ Sci. 2002 Mar;15(1):83–96.
Li M, Shi Z. A Prospective Association of Nut Consumption with Cognitive Function in Chinese Adults Aged 55+_China Health and Nutrition Survey. J Nutr Health Aging. 2019;23(2):211–6.
Hao G, Liu K, Halbert JD, Chen H. Dietary sodium and potassium and risk of diabetes: A prospective study using data from the China Health and Nutrition Survey. Diabetes Metab [Internet]. 2019; Available from: https://www.sciencedirect.com/science/article/abs/pii/S1262363619301880?via%3Dihub.
IDF Clinical Practice Recommendations for managing Type 2 Diabetes in Primary Care. Int Diabetes Fed [Internet]. 2017; Available from: .
Ainsworth BE, Haskell WL, Leon AS, Jacobs DRJ, Montoye HJ, Sallis JF, et al. Compendium of physical activities: classification of energy costs of human physical activities. Med Sci Sports Exerc. 1993 Jan;25(1):71–80.
Ainsworth BE, Haskell WL, Whitt MC, Irwin ML, Swartz AM, Strath SJ, et al. Compendium of physical activities: an update of activity codes and MET intensities. Med Sci Sports Exerc. 2000 Sep;32(9 Suppl):498–504.
Ainsworth BE, Haskell WL, Herrmann SD, Meckes N, Bassett DRJ, Tudor-Locke C, et al. 2011 Compendium of Physical Activities: a second update of codes and MET values. Med Sci Sports Exerc. 2011 Aug;43(8):1575–81.
Ng S, Popkin B. Time use and physical activity: a shift away from movement across the globe. Obes Rev. 2012 Aug;13(8):659–80.
Ng S, Norton E, Popkin B. Why have physical activity levels declined among Chinese adults? Findings from the 1991–2006 China Health and Nutrition Surveys. Soc Sci Med. 2009 Apr;68(7):1305–14.
Du S, Wang H, Zhang B, Popkin B. Dietary Potassium intake remains low and sodium intake remains high, and most sodium is derived from home food preparation for Chinese Adults, 1991–2015 trends. J Nutr. 2020 May;150(5):1230–9.
Hipgrave D, Chang S, Li X. Salt and Sodium Intake in China. JAMA. 2016 Feb;315(7):703–5.
Wang HJ, Wang ZH, Zhang JG, Du WW, Su C, Zhang J, et al. Trends in dietary fiber intake in Chinese aged 45 years and above, 1991–2011. Eur J Clin Nutr. 2014;68:619–22.
Wang Z, Zhai F, Zhang B, Popkin BM. Trends in Chinese snacking behaviors and patterns and the social-demographic role between 1991 and 2009. Asia Pac J Clin Nutr. 2012;21(2):253–62.
Wang Z, Siega-Riz AM, Gordon-Larsen P, Cai J, Adair LS, Zhang B, et al. Diet quality and its association with type 2 diabetes and major cardiometabolic risk factors among adults in China. Nutr Metab Cardiovasc Dis. 2018 Oct 1;28(10):987–1001.

ACSDietPaperSupplement83.docx

Download PDF

Version 1

posted

You are reading this latest preprint version

Characterizing the Urban Diet: Development of an urbanized diet index

Status:

Version 1

Abstract

Figures

Background

Methods

Study Participants

Diet Data

Health Outcome Variables

Hypertension (HTN)

Overweight

Type II Diabetes Mellitus (T2DM)

Urbanization Index

Covariates

Urbanized Diet Index Development

Select Dietary Variables of Interest

Score Individual Dietary Variables

Analytical steps to create the urbanized diet index

Results

Select Dietary Variables of Interest

Score Individual Dietary Variables

Create Candidate Urbanized Diet Indices

Select a Final Urbanized Diet Index

Test the Stability of the Final Urbanized Diet Index across Sociodemographic Factors

Test the Association between the Final Urbanized Diet Index and CMD risk

Discussion

Conclusions

Abbreviations

Declarations

Acknowledgements

References

Supplementary Files

Status:

Version 1