Next Article in Journal
Aruncus dioicus var. kamtschaticus Extract Ameliorates Psoriasis-like Skin Inflammation via Akt/mTOR and JAK2/STAT3 Signaling Pathways in a Murine Model
Next Article in Special Issue
Policy to Decrease Low Birth Weight in Indonesia: Who Should Be the Target?
Previous Article in Journal
The Relation between Resistin (−420C/G) Single Nucleotide Variant, Resistin Serum Concentration, Carbohydrate, and Lipid Parameters and Fried Food Taste Preference in Patients with Hypertriglyceridemia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Role of Residence Area on Diet Diversity and Micronutrient Intake Adequacy in Urban and Rural Costa Rican Adolescents

by
Rafael Monge-Rojas
1,*,
Rulamán Vargas-Quesada
1 and
Georgina Gómez
2
1
Nutrition and Health Unit, Costa Rican Institute for Research and Education on Nutrition and Health (INCIENSA), Ministry of Health, Tres Ríos 4-2250, Costa Rica
2
Biochemistry Department, School of Medicine, University of Costa Rica, San José 11501-2060, Costa Rica
*
Author to whom correspondence should be addressed.
Nutrients 2022, 14(23), 5093; https://doi.org/10.3390/nu14235093
Submission received: 16 November 2022 / Accepted: 24 November 2022 / Published: 1 December 2022
(This article belongs to the Special Issue Nutrition Status and Policies in Low- and Middle-Income Countries)

Abstract

:
Dietary diversity might be essential to meet nutritional demands during adolescence. Diet diversity among 818 urban and rural Costa Rican adolescents aged 13–18 years was studied using the Minimum Dietary Diversity Score for Women. The Nutrient Adequacy Ratio (NAR) was calculated for 11 nutrients to estimate the nutrient adequacy of the diet. A NAR < 0.7 was considered inadequate for micronutrient intake. The optimal Diet Diversity Score (DDS) cut-off point for this study was 4, established using receiver-operating characteristic curves. The mean DDS for the overall sample was 4.17 ± 1.43, although DDS was significantly higher in adolescents from rural vs. urban areas (4.33 ± 1.43 vs. 4.00 ± 1.42, p-value = 0.001). The odds of having a diverse diet were 62% higher in rural vs. urban adolescents. Overall, 80–95% of adolescents reached a NAR ≥ 0.70 for 8 nutrients except for calcium, zinc, and vitamin A. The residence area plays a key role in adolescent dietary diversity. Although overall DDS was low, foods that make up the rural adolescent diet were nutritionally dense enough to satisfy the EAR for most micronutrients. A high DDS is not necessarily required for the diet to meet most micronutrient demands in adolescence. Improved dietary adequacy of vitamin A, zinc, and calcium is required due to the importance of these micronutrients in maintaining optimal health.

1. Introduction

The need for energy and nutrients increases significantly in adolescence to support the rapid rate of growth and development [1]. During this period of life, achieving an adequate nutrient intake is crucial; however, there is evidence that dietary quality declines from childhood to adolescence [2,3,4], and that adolescents, as a group, have poor eating habits that do not meet dietary recommendations for micronutrient intake [5]. Evidence suggests that healthy dietary habits established during adolescence continue into adulthood [6,7,8,9]; hence, adolescence has been identified as the best time to achieve dietary modifications aimed at enhancing health-conscious dietary habits [6,9,10]. Assessing the quality of the adolescent diet is necessary to prevent deficiencies that may limit optimal growth and development and increase the risk of developing cardiometabolic diseases in adulthood.
Several food-based and nutrient-based dietary quality indices have been developed in Europe and North America to assess the overall healthfulness of dietary intake in adolescents; however, these indices require intensive dietary evaluation, complex analysis, and a significant burden of resources [11]. Nevertheless, among dietary indices, the Dietary Diversity Score (DDS) has been identified at the individual level as a promising measurement tool, particularly for use in low-income countries (LICs) and middle-income countries (MICs) [12], because of its simplicity of implementation and low financial support required.
Food diversity has long been recognized as a key element of high-quality diets, based on the premise that no single food can provide the adequate amount of nutrients needed to maintain optimal health [12]; therefore, eating a variety of foods has been a longstanding recommendation in the dietary guidelines of various countries [13]. Lack of dietary diversity is usual among underserved adolescents from LICs and MICs, mainly because their diets are largely based on starchy staples, with little or no animal products, and few fresh fruits and vegetables [14,15,16].
Dietary diversity may be linked to four dietary archetypes arising from socio-ecological traps in the food systems prevalent in various countries, even in rural and urban areas of the same country [17,18]. The Traditional diet archetype is characterized by providing a healthy intake of various nutrients, whereas the Mixed diet archetype represents ample access to a large number of traditional and Western foods. The Undernourishing diet archetype is characterized by inadequate quality and potentially inadequate quantity of food, while the Overnourishing diet archetype tends to replace the traditional diet with a Western diet, which has an insufficient quantity of micronutrients rich and nutritious foods [17].
Due to the lack of a consensual definition of dietary diversity, several DDSs have been proposed. A recent systematic scoping review showed that many DDSs vary in terms of format (e.g., food item-based, food group-based, dietary guidelines-based, functional diversity-based, and others) and the dietary assessment of food intake (on a single 24-h recall or Food Frequency Questionnaire). However, most are consistent on the concept of counting foods or food groups consumed over a period [12].
Studies conducted on adolescents of LICs and MICs show that the DDS can range from 3.35 points [14] to 6.25 points [19]; nevertheless, these values are not comparable since the methodology to determine DDSs varies considerably from one study to another. What is consistent is that the higher the DDS, the more sufficient the diet is to meet micronutrient requirements, and that rural adolescents have a lower dietary diversity than urban ones [20,21], except in rural adolescents from India, where the DDS is equal to 8 points [22]. Nevertheless, this study has some relevant methodological aspects that make it notably different from other studies carried out in rural adolescents. First, foods were categorized into thirteen groups, and second, a simple count of food groups was made to calculate individual DDS, which ranged from 1 to 13.
To avoid so much variability in studying dietary diversity, it is necessary to establish a simple methodology that suits the context and food culture of a country or large geographic area. The Minimum Dietary Diversity for Women (MDD-W) indicator was developed by the Food and Agriculture Organization (FAO) of the United Nations and partners [23] to fill the need for a simple food-based indicator for measuring dietary diversity and micronutrient adequacy as key dimensions of diet quality in women of reproductive age. However, this indicator has also been used to assess the adequacy of micronutrient intake in children [24], adult men [25,26], and adolescent boys [25], evidencing its applicability in other population groups.
According to the MDD-W, women who consume 15 g or more from at least five out of ten defined food groups over a 24-h recall period are classified as having a minimally adequate DDS to meet the recommended intake of 11 important micronutrients (vitamin A, thiamine, riboflavin, niacin, vitamin B6, folate, vitamin B12, vitamin C, calcium, iron, and zinc). However, several studies [25,27,28] have pointed out the need to analyze the sensitivity and specificity of different DDS cut-off points to define which is most appropriate for reducing the risk of over-identifying population groups with nutritionally inadequate diets. Under this premise, other studies have established different DDS cut-off points; for example, ≥4 for Ugandan adolescents [15] and ≥5 for Mexican adolescents [25].
The aims of the current study were: (1) to determine the DDS cut-off point for the Costa Rican adolescent diet, (2) to compare dietary diversity between urban and rural adolescents, and (3) to estimate the proportion of adolescents that do not meet micronutrient dietary recommendations according to area of residence.

2. Materials and Methods

2.1. Study Population and Setting

Data came from a cross-sectional sample of adolescents (13–18 years old; 7–11th graders) enrolled in 18 schools (10 from urban areas, and 8 from rural areas) in the province of San José, Costa Rica, in 2017. Most Costa Rican adolescents (80%) are enrolled in school [29] and San José has the highest adolescent concentration (30%) in the country [30].
School selection and sample size were determined assuming a sampling error for a population proportion with correction for a finite population [31] following three phases: (1) Schools were selected using a proportional-size probability method to represent urban and rural municipalities within the province [32]. (2) Ten classrooms (two from each grade) were selected in each school using simple random sampling. In each classroom, all students were invited to participate in the study, provided an informed assent, and asked to obtain a signed informed consent from their parents if they were interested in participating in the study. (3) Among those who returned signed informed assent and informed consent forms, students were randomly selected to participate in the study. Close to 5% of the initial sample chose not to participate in the study before it started. The final study sample was 818 adolescents aged 13 to 18 years.
Data were collected during school hours in a private classroom. The Bioethics Committee of the Costa Rican Institute for Research and Education on Nutrition and Health (INCIENSA) approved all study protocols, and all guidelines for human subject research were strictly followed.

2.2. Sociodemographic Variables

A paper-based questionnaire was used to collect data on sex, age, residence area, socioeconomic status (SES), and nutritional status.

2.3. Anthropometric Assessment

Participants’ height and weight were determined by trained nutritionists following the methodology described by Preedy [33]. Body mass index (BMI) values were calculated, and adolescent nutritional status was determined using the BMI-for-age Z score recommended by WHO [34]: <−2: underweight; ≥−2 and <+1: healthy weight (eutrophy); ≥+1 and <+2: overweight, and ≥+2: obese. For data analysis, nutritional status was dichotomized as non-overweight (underweight and eutrophy) and overweight/obesity (overweight and obesity).

2.4. Dietary Intake Assessment

Dietary intake data were collected via 3-day food records completed by the participants in real time and reviewed by nutritionists. Participants were asked to complete the 3-day food records on two weekdays (Monday, Tuesday, Thursday, or Friday) and one weekend day (Saturday or Sunday). Half of the participants were randomly selected to record the foods and drinks that they consumed on Thursday, Friday, and Saturday, and the rest were asked to record their intake on Sunday, Monday, and Tuesday. Data were collected during nine months of the school year (February-November), reflecting seasonal variations for Costa Rica: rainy season (May-November) and dry season (December–April). The goal was to ensure that the data captured daily and seasonal variability in food consumption.
At each school, six trained nutritionists provided intake forms to the participants and taught them how to complete accurate food records by writing down detailed descriptions of what they ate and drank from the time they woke up in the morning to the time they went to bed at night, for three consecutive days. Adolescents were asked to write down the brand names of foods when appropriate, methods of preparation, and recipes for all the dishes and drinks whenever possible. The nutritionists taught the participants how to estimate the portion sizes of the foods and drinks they consumed using an established portion-size manual developed for Costa Rica [35]. The manual includes photographs and diagrams of typical local foods and their preparations, including 3–6 different portion sizes. Adolescents were also instructed to report portion size using measurements based on household utensils or volume and mass units.
Given the challenges with incomplete and inaccurate data when recording self-reported dietary data in young populations and specific demographic groups [36], the completed 3-day food records were thoroughly reviewed by the nutritionists by conducting one-on-one reviews with each participant during school hours. At this interview, the nutritionists inquired about commonly missed items or ingredients (i.e., added sweeteners, added fats, candies, beverages), added details about the kind of food or drinks that were consumed, verified or added any omitted portion sizes, and clarified any illegible items. The nutritionists used food models, fresh foods, and different utensils to verify serving and portion sizes. Only the food intake recorded by the adolescents on the second day (Monday or Friday) was used to determine the dietary diversity score. The second day of the 3-day food record was chosen because it was a weekday for the whole sample, making sure to select a day of habitual adolescent dietary intake. Meanwhile, the first day (Thursday or Sunday) or the third one (Saturday or Tuesday) would have been a weekend day at least for half the sample. This might induce bias in the dietary intake analysis since food consumption and diet quality differ between weekdays and weekends [37,38]. To corroborate this premise, a post-hoc sensitivity analysis was performed using the dietary data from the first day of the 3-day food records for all the analyses performed with data from day 2. This showed important differences in the results obtained when comparing the first and second days of the food record, which confirmed that including a weekend day modifies the results negatively (Table S1–S6, Figures S1 and S2).

2.5. Diet Diversity Score

The DDS was assessed using the MDD-W of reproductive age developed by FAO in 2016 as a valid indicator of micronutrient adequacy in women of childbearing age [23]. This is a proposal of a single indicator to assess dietary quality that has been used in populations different from the one for which it was originally designed, such as adult men [25,26] and adolescent boys [15,25].
The DDS was calculated at the individual level based on the Women’s Dietary Diversity Score Project food group classification [23]; the first step involved adding the total of food groups reported as consumed on the second day of the 3-day food records (Friday or Monday). This was done because the established methodology by FAO to analyze the MDD-W should be done based only on a single 24-h period over one day and night [23].
To determine dietary diversity, a score was created based on the consumption of at least 15 g/day of the following 10 food groups: (1) starchy staples (grains, roots and tubers, and plantains), (2) flesh foods (meat, poultry, and fish), (3) dark green leafy vegetables, (4) other vitamin A-rich fruits and vegetables (other vitamin A-rich F&V), (5) other vegetables, (6) other fruits, (7) pulses (beans, peas, and lentils), (8) milk and milk products, (9) eggs, and (10) nuts and seeds. One point was assigned if dietary intake of a food group was ≥15 g/day and 0 points if intake was <15 g/day; thus, the score ranged from 0 to 10 points. Higher scores would indicate a diet with greater diversity, as more food groups were eaten in amounts greater than or equal to 15 g/d. DDS analyses were performed by sex, residence area, socioeconomic status (SES), and nutritional status.

2.6. Nutrient Adequacy

To estimate the nutrient adequacy of the diet, a Nutrient Adequacy Ratio (NAR) was calculated for the intake of 11 nutrients. The NAR for a given nutrient is the ratio of a subject’s intake to the Estimated Average Requirement (EAR) for the adolescent’s sex and age category, according to the guidelines of the National Academy of Medicine of the United States [39].
The EAR is the appropriate dietary reference intake to use when assessing the adequacy of group intakes [40]. A NAR value less than 70% (NAR < 0.7) was considered indicative of inadequate nutrient intake, as suggested by Oldewage-Theron [41]. As an overall measure of nutrient adequacy, the Mean Adequacy Ratio (MAR) was calculated as has been previously described [42].
M e a n   A d e q u a c y   R a t i o   M A R = N A R   e a c h   t r u n c a t e d   a t   1 N u m b e r   o f   n u t r i e n t s
NAR was truncated at 1, so a nutrient with a high NAR could not compensate for a nutrient with a low NAR. MAR > 0.7 has been considered as the “gold standard” of nutritionally adequate intake [28].

2.7. Diet Diversity Score Cut-Off

According to the FAO methodology, the diet can be considered diverse if 5 or more food groups are consumed, and non-diverse if the consumption is less than 5 food groups [23]. However, the 5-food-group cut-off has been questioned and several studies have recommended determining the cut-off point according to the study population [25,27,28,41] to optimize balance among sensitivity, specificity, and positive predictive values for MAR [27]. In consequence, receiver-operating characteristic (ROC) curves were generated to obtain the adequate DDS cut-off for this study. As shown in Figure 1, different DDS cut-off points were tested for sensitivity and specificity versus different definitions of a nutritionally adequate diet, with MAR ranging from 0.65 to 0.75.
An optimal DDS cut-off point can identify as many nutritionally inadequate diets as possible (high sensitivity) without losing too much ability to identify nutritionally adequate diets (specificity). Therefore, nutritionally inadequate diets, defined here as MAR and DDS below a cut-off point, were defined as true positives. Nutritionally adequate diets, with MAR and DDS above a cut-off point, were defined as true negatives [28]. Figure 1 shows that the sensitivity and specificity were slightly influenced by changes in the MAR cut-off points, so we used a MAR of 0.70 as the cut-off point for a nutritionally adequate diet.
A MAR cut-off point equal to 0.70 and a DDS cut-off point equal to 4 resulted in 71.5% specificity and 78.0% sensitivity. If the DDS was increased to 5, specificity decreased to 43.5%, while sensitivity increased to 91.5%. Thus, a DDS cut-off point of 4 is, in this study, the one that resulted in higher than 70% values for both specificity and sensitivity. This cut-off point is considered robust, since the Youden J-index value (0.495) was the closest to 1, compared to the values observed while analyzing the other DDS cut-off points. Also, this DDS cut-off point offers a balance between a high negative predictive value and a modest positive predictive value for MAR (0.97 and 0.23, respectively), reducing the risk of over-identifying adolescents with nutritionally inadequate diets [43].

2.8. Statistical Analyses

Data were reported as means ± standard deviations (SD) for continuous variables, and frequencies (%) for categorical variables. The between-groups comparisons (age, sex, SES, and nutritional status) in the proportions of adolescents according to residence area were tested using the Wilcoxon test, or the chi-square test with Bonferroni post hoc test when required. The DDS was stratified by sex, SES, and nutritional status, and the difference between groups was tested with Student’s t or ANOVA tests as required. The odds of having a diverse diet according to sociodemographic and nutritional status variables were determined by a bivariate and multivariate logistic regression analysis. All tests were two-tailed, and p-values < 0.05 were considered statistically significant.
Consumption of food groups and nutrient intake were stratified by DDS (<4 or ≥4), and differences between groups were tested using the Wilcoxon test. Intakes were adjusted for 1000 kcal to assess the dietary diversity effect (DDS < 4 vs. DDS ≥ 4) on macronutrient and micronutrient intake, independently of diet quantity. NAR and MAR were stratified by diet diversity and comparisons were made using the Wilcoxon test. Spearman’s rank correlation coefficients were used to determine the association between the DDS and NAR/MAR scores. Finally, a chi-square test was used to estimate significant differences in the distribution of participants with NAR/MAR ≥ 0.70 according to DDS and residence area.
Data were analyzed using STATA software (version 14.1, 2015, College Station, TX, USA). Epi Info™ software (version 3.5.4, 2008) was used to process data from the 3-day food records, with information from the database of the School of Human Nutrition of the University of Costa Rica [44]. This database contains 1655 items from the following sources: 1307 foods from USDA, 80 foods from the food composition tables of the Institute of Nutrition for Central America and Panama (INCAP), 254 foods from recipes commonly consumed in Costa Rica (e.g., frescos, pastries, bread, cookies, and desserts), 13 foods from the Costa Rican mandatory fortification food group and updated information until 2017 (year in which the study data were collected).

3. Results

3.1. General Characteristics According to Residence Area

The mean age of the study sample was 15.0 (±1.7) years, with 63.6% girls, 50.2% urban, and 32.6% overweight/obese adolescents (Table 1). There was no significant difference between adolescents from urban and rural areas regarding age, sex, and overweight/obesity status (p > 0.05). Compared to urban adolescents, a higher proportion (66.9%) of rural adolescents was classified as low SES (p < 0.0001). Likewise, a higher proportion of urban adolescents was enrolled in private schools compared to rural adolescents.

3.2. Diet Diversity Score According to General Characteristics

The mean DDS for the overall sample was 4.17 ± 1.43 points out of 10.0 possible maximum score, with 66.5% of participants reaching 4 or more points, which is the criterion for a diverse diet (Table 2). DDS was higher in participants from rural vs. urban areas (4.33 ± 1.43 vs. 4.00 ± 1.42) (p = 0.001), and 72.0% of rural adolescents reached the criterion for a diverse diet, compared with 61.1% for urban adolescents (p = 0.001). No significant differences were found when analyzing the DDS by sex, SES, or nutritional status.

3.3. Consumption of Food Groups

As expected, adolescents with a diverse diet (DDS ≥ 4) showed a significantly higher consumption of all food groups than participants with a non-diverse diet (DDS < 4) (p < 0.05) (Table 3). Adolescents from rural areas had higher consumptions of starchy staples (Δ = 49.0 g/d; p < 0.0001); pulses (Δ = 31.0 g/d; p < 0.0001); other vegetables (Δ = 19.3 g/d; p <0.0001); eggs (Δ = 4.2 g/d; p = 0.029); and dark green leafy vegetables (Δ = 3.7 g/d; p = 0.016). In contrast, urban adolescents had a higher dietary intake of milk and milk products (Δ = 27.7 g/d; p = 0.010) (Table 3). There were no differences in the consumption of flesh foods, other fruits, other vitamin A-rich F&V, and nuts and seeds, between urban and rural adolescents (p > 0.05).
Figure 2 shows the proportion of adolescents who consumed the food groups that presented significant differences between urban and rural areas. The proportion of urban and rural adolescents that consumed the starchy food group was similar, possibly because the largest proportion (35%) of this food group corresponds to rice (data not shown), which is one of the main staple foods of the Costa Rican diet.
A higher proportion of rural vs. urban adolescents consumed at least 15 g of the following food groups: pulses (73.5% vs. 50.6%, p < 0.0001); other vegetables (54.1% vs. 39.2%, p < 0.0001); eggs (32.2% vs. 24.6%; p = 0.016); and dark green leafy vegetables (6.1% vs. 2.7%; p = 0.016). On the contrary, a higher proportion of urban vs. rural adolescents consumed 15 g or more of milk and milk products (57.4% vs. 47.7%; p = 0.005).
Table 4 shows a bivariate and multivariate logistic regression analysis with the sociodemographic variables related to having a diverse diet. According to the multivariate model, the odds of having a diverse diet are 62% higher in rural adolescents, after adjusting for sex, age, SES, and nutritional status (OR 1.62, 95% CI: 1.19–2.20, p-value = 0.002).

3.4. Nutrient Adequacy Ratio and Correlation with Diet Diversity Score

Table 5 shows NAR values for each nutrient. Mean NARs for calcium and vitamin A are below 70% of the EAR (NAR < 0.70), independently of diet diversity. Adolescents with a diverse diet (DDS ≥ 4) reported significantly higher NAR values compared to participants with a non-diverse diet (DDS < 4) for the 11 micronutrients assessed. MAR values were, as expected, significantly higher in the diverse diet group. For the whole sample, MAR and all NARs correlated positively and significantly with the DDS, showing an increase in the micronutrient adequacy ratio for all these vitamins and minerals as DDS increased. Higher correlations were obtained for zinc (0.354), calcium (0.383), and MAR (0.433).
Table 6 shows the proportion of participants reaching a NAR ≥ 0.70, according to diet diversity and residence area. Overall, less than 70% of the adolescents reached a NAR ≥ 0.70 for calcium, zinc, and vitamin A. As expected, for all nutrients the proportion of adolescents reaching a NAR ≥ 0.70 with a diverse diet (DDS ≥ 4) was significantly higher than for adolescents with a non-diverse diet (DDS < 4). Regarding residence area, a significantly higher proportion of rural vs. urban adolescents reached a NAR ≥ 0.70 for zinc (69.0% vs. 61.6%; p = 0.025), vitamin C (87.0% vs. 81.8%; p = 0.040), and vitamin B6 (87.0% vs. 79.3%; p = 0.003). In contrast, this proportion was higher in urban adolescents only for calcium (29.9% vs. 22.4%; p = 0.014). No significant differences were observed between rural and urban adolescents for other nutrients/MAR.

4. Discussion

This study shows that the DDS cut-off point for urban and rural Costa Rican adolescents was similar to the one for Ugandan adolescents (≥4) [15], and lower than the one for Mexican adolescents (≥5) [25]. The cut-off point of 4 was determined through sensitivity and specificity analyses, which are inversely related. For nutrition promotion and prevention of nutrient inadequacy, high sensitivity is desirable to accurately identify adolescents with inadequate micronutrient intake, as false positives pose no significant risk in this case [28]. However, the determination of a DDS cut-off point with a balance between sensitivity and specificity is essential to strengthening the analysis of dietary diversity, since failure to do so may increase the risk of overestimating the proportion of adolescents with inadequate micronutrient intakes. For example, in this study, a DDS cut-off point of 5 would have given a specificity of 43.5%, with a false positive rate of 56.5% (28.0% higher than with a DDS cut-off point of 4), which would predict a higher proportion of adolescents with inadequate micronutrient intakes.
In contrast to the evidence from LICs, where diet diversity is lower in rural areas, in Costa Rica (a MIC) diet diversity is higher in rural areas. Among the main causes of low diet diversity in LICs are dependence on home-grown staples and legumes, and limited access (economic and physical) to other key healthy foods that are available at the local markets—typical of the Undernourishing diet archetype [17].
In Costa Rica, food availability and access in rural areas (particularly in San José, where this study was conducted) are not very different from the urban areas of the province [45]. The conceptualization of ‘rural’ in Costa Rica may differ from what is established in LICs. Currently, the diversification of the rural productive structure no longer depends exclusively on agricultural activities. Instead, there has been a significant growth of service activities, which are increasingly not associated with agriculture. In addition, there is a more functional integration of rural and urban settings, which can be characterized as a process of counter-urbanization. This process has been facilitated by the development of new communication technologies and an increased income in rural populations, which help bridge the gap between the ‘modern’ and rural worlds [46].
However, rural areas still retain certain unique qualities such as family food patterns. It has been reported previously that the rural adolescent diet tends to be more traditional than the urban adolescent diet. The traditional diet in rural adolescents is characterized by adequate quantities, diversity, and quality of food to sustain health and nutrition [45]. This can help explain why rural adolescents have 62% odds of having a more diverse diet than their urban counterparts. The diet diversity observed in rural areas suggests that the Traditional diet archetype prevails in those areas.
Although the overall DDS was low (4.17 ± 1.43 out of 10 points) and similar to Uganda’s, there was a low proportion of rural and urban adolescents with inadequate intakes. Hence, even if a low DDS is an indicator of a monotonous diet [12], it does not necessarily lead to an inadequate micronutrient intake, as adequacy will depend on the foods that fall within that monotony. Beans and rice, which are the staple foods of the Costa Rican diet, provide a great variety of nutrients. Pulses are natural sources of thiamine, riboflavin, niacin, vitamin B6, folate, iron, and zinc [47], while fortified white rice provides thiamine, niacin, folic acid, and zinc [48].
In contrast to urban adolescents, a large proportion of rural adolescents consumes at least 15 g/day of pulses, dark green leafy vegetables, other vegetables, and affordable protein sources such as eggs, contributing to an improved DDS in rural areas —even if slightly above the minimum diet diversity set for this study. This finding suggests that neither a high DDS nor a high family income is necessarily required for the diet to meet the majority of micronutrient demands in adolescence.
The lower DDS found in urban areas may be reflecting a more established nutritional transition. It has been thoroughly demonstrated that rapid urbanization and high family income are critical drivers of this transition [49,50,51,52]. Nutritional transition is characterized by the progressive shift from a traditional diet to one rich in animal products and low in pulses [53], corresponding to the Overnourishing diet archetype [17], and evidenced in the urban adolescent diet. Therefore, increased consumption of costly foods (e.g., milk and milk products) at the expense of more affordable choices (pulses and eggs) does not guarantee a higher dietary diversity.
Although the prevalence of inadequate micronutrient intakes was low for most micronutrients, interventions need to be designed to increase the consumption of vitamin A-rich fruits and vegetables, eggs, and flesh foods—particularly in urban adolescents—to improve the dietary adequacy of vitamin A and zinc. Likewise, promoting the consumption of calcium-rich dairy products is urgently required, since the intake of this mineral is critical for both urban and rural adolescents.
We acknowledge the limitations of our study. Although the MDD-W was not designed for men, other studies conducted in adult and adolescent males [15,25,26] have demonstrated its usefulness in determining the adequacy of micronutrient intake; hence, its use in males should be considered and analyzed with greater deference. The MDD-W is a simple and affordable tool requiring very straightforward analysis, and the determination of food consumption does not involve complex dietary assessment methods.
Our study also has strengths that must be emphasized. First, we conducted the DDS analysis using a specific cutoff point for the population group, providing a balance between a high negative predictive value and a modest positive predictive value for MAR, and reducing the risk of over-identifying adolescents with nutritionally inadequate diets. Second, diet evaluation through intake records is considered the gold standard among other dietary assessment instruments [54]. Finally, while participants may not be representative of Costa Rica as a whole, the province of San Jose is where most of the national adolescent population (30%) is concentrated [30].

5. Conclusions

Our study shows that the area of residence plays a key role in adolescent dietary diversity. Although overall DDS was low, the foods that make up the rural adolescent diet were nutritionally dense enough to satisfy the EAR for most micronutrients. Therefore, our study suggests that a high DDS is not required to achieve ≥ 0.70 EAR intake for most micronutrients as long as the foods included in the diet have a high micronutrient density either naturally or through fortification. However, greater consumption of vitamin A-rich fruits and vegetables, eggs, and meats is required—particularly in urban adolescents—to improve the dietary adequacy of vitamin A and zinc. In addition, increased consumption of calcium-rich food sources is urgently required for all adolescents because inadequate intake of this mineral is notoriously high.

Supplementary Materials

The following are available online at: https://www.mdpi.com/article/10.3390/nu14235093/s1. Figure S1. Sensitivity and specificity of different DDS cut-off points with MAR cut-off points ranging from 0.65 to 0.75, performed as sensitivity analysis for day one of 3-day food record; Table S1. Diet Diversity Score (DDS) and proportion of participants reaching the minimum DDS 1 (n = 818), performed as sensitivity analysis for day one of 3-day food record; Table S2. Food group consumption of Costa Rican adolescents, according to Diet Diversity Score (DDS) and residence area, performed as sensitivity analysis for day one of 3-day food record; Figure S2. Proportion of participants consuming at least 15 g of each food group over one-day food record, according to residence area, performed as sensitivity analysis for day one of 3-day food record; Table S3. Bivariate and multivariate logistic regression analysis: Sociodemographic variables related to having a diverse diet, according to Diet Diversity Score (DDS) (n = 818), performed as sensitivity analysis for day one of 3-day food record; Table S4. Mean Nutrient Adequacy Ratio (NAR), according to Diet Diversity Score (DDS), performed as sensitivity analysis for day one of 3-day food record; Table S5. Proportion of participants reaching Nutrient Adequacy Ratio (NAR) ≥ 0.70, according to Diet Diversity Score (DDS) and residence area), performed as sensitivity analysis for day one of 3-day food record; Table S6. Mean Nutrient Adequacy Ratio (NAR) and Diet Diversity Score (DDS), comparing days of the 3-day food record, performed as sensitivity analysis.

Author Contributions

R.M.-R.: Conceived and designed the study, collected, and interpreted data, and wrote the manuscript. R.V.-Q.: Contributed to data analysis and interpretation and assisted in writing the manuscript. G.G.: Contributed to data analysis and assisted in writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Tobacco Control Program of the Department for Strategic Planning and Evaluation of Health Actions of the Costa Rican Ministry of Health (agreements DM-FG-4854-14 and DM-FG.1748-2018).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and in accordance with Law 9234, “Regulatory Law of Biomedical Research”, which regulates biomedical research in Costa Rica. The study protocol was approved by the Bio-ethics Committee of the Costa Rican Institute for Research and Education on Nutrition and Health (INCIENSA) under number IC-2007-01.

Informed Consent Statement

We obtained oral and written informed assents from all adolescents included in the study; likewise, we obtained written informed consents from their parents.

Data Availability Statement

The data presented in this study are available on request to the corresponding author.

Acknowledgments

The authors are grateful to Ana Leonor Rivera for her support in data collection.

Conflicts of Interest

The authors declare no conflict of interest. The funder had no role in the design of the study, collection, analyses, or interpretation of data, writing of the manuscript, or decision to publish the results.

References

  1. Norris, S.A.; Frongillo, E.A.; Black, M.M.; Dong, Y.; Fall, C.; Lampl, M.; Liese, A.D.; Naguib, M.; Prentice, A.; Rochat, T.; et al. Nutrition in Adolescent Growth and Development. Lancet 2022, 399, 172–184. [Google Scholar] [CrossRef]
  2. Demory-Luce, D.; Morales, M.; Nicklas, T.; Baranowski, T.; Zakeri, I.; Berenson, G. Changes in Food Group Consumption Patterns from Childhood to Young Adulthood: The Bogalusa Heart Study. J. Am. Diet. Assoc. 2004, 104, 1684–1691. [Google Scholar] [CrossRef]
  3. Story, M.; Neumark-Sztainer, D.; French, S. Individual and Environmental Influences on Adolescent Eating Behaviors. J. Am. Diet. Assoc. 2002, 102, S40–S51. [Google Scholar] [CrossRef] [PubMed]
  4. Birch, L.; Savage, J.S.; Ventura, A. Influences on the Development of Children’s Eating Behaviours: From Infancy to Adolescence. Can. J. Diet. Pr. Res. 2007, 68, s1–s56. [Google Scholar]
  5. Jenkins, S.; Horner, S.D. Barriers That Influence Eating Behaviors in Adolescents. J. Pediatr. Nurs. 2005, 20, 258–267. [Google Scholar] [CrossRef] [PubMed]
  6. Cruz, F.; Ramos, E.; Lopes, C.; Araújo, J. Tracking of Food and Nutrient Intake from Adolescence into Early Adulthood. Nutrition 2018, 55–56, 84–90. [Google Scholar] [CrossRef] [PubMed]
  7. Movassagh, E.; Baxter-Jones, A.; Kontulainen, S.; Whiting, S.; Vatanparast, H. Tracking Dietary Patterns over 20 Years from Childhood through Adolescence into Young Adulthood: The Saskatchewan Pediatric Bone Mineral Accrual Study. Nutrients 2017, 9, 990. [Google Scholar] [CrossRef] [PubMed]
  8. Craigie, A.M.; Lake, A.A.; Kelly, S.A.; Adamson, A.J.; Mathers, J.C. Tracking of Obesity-Related Behaviours from Childhood to Adulthood: A Systematic Review. Maturitas 2011, 70, 266–284. [Google Scholar] [CrossRef] [Green Version]
  9. Mikkilä, V.; Räsänen, L.; Raitakari, O.T.; Pietinen, P.; Viikari, J. Consistent Dietary Patterns Identified from Childhood to Adulthood: The Cardiovascular Risk in Young Finns Study. Br. J. Nutr. 2005, 93, 923–931. [Google Scholar] [CrossRef] [PubMed]
  10. Schneider, B.C.; Dumith, S.D.C.; Lopes, C.; Severo, M.; Assunção, M.C.F. How Do Tracking and Changes in Dietary Pattern during Adolescence Relate to the Amount of Body Fat in Early Adulthood? PLoS ONE 2016, 11, e0149299. [Google Scholar] [CrossRef] [PubMed]
  11. Marshall, S.; Burrows, T.; Collins, C.E. Systematic Review of Diet Quality Indices and Their Associations with Health-Related Outcomes in Children and Adolescents. J. Hum. Nutr. Diet. 2014, 27, 577–598. [Google Scholar] [CrossRef] [PubMed]
  12. Verger, E.O.; Le Port, A.; Borderon, A.; Bourbon, G.; Moursi, M.; Savy, M.; Mariotti, F.; Martin-Prevel, Y. Dietary Diversity Indicators and Their Associations with Dietary Adequacy and Health Outcomes: A Systematic Scoping Review. Adv. Nutr. 2021, 12, 1659–1672. [Google Scholar] [CrossRef]
  13. Herforth, A.; Arimond, M.; Álvarez-Sánchez, C.; Coates, J.; Christianson, K.; Muehlhoff, E. A Global Review of Food-Based Dietary Guidelines. Adv. Nutr. 2019, 10, 590–605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Baxter, J.B.; Wasan, Y.; Islam, M.; Cousens, S.; Soofi, S.B.; Ahmed, I.; Sellen, D.W.; Bhutta, Z.A. Dietary Diversity and Social Determinants of Nutrition among Late Adolescent Girls in Rural Pakistan. Matern. Child Nutr. 2022, 18, e13265. [Google Scholar] [CrossRef] [PubMed]
  15. Isabirye, N.; Bukenya, J.N.; Nakafeero, M.; Ssekamatte, T.; Guwatudde, D.; Fawzi, W. Dietary Diversity and Associated Factors among Adolescents in Eastern Uganda: A Cross-Sectional Study. BMC Public Health 2020, 20, 534. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Islam, M.R.; Rahman, S.M.; Tarafder, C.; Rahman, M.M.; Rahman, A.; Ekström, E.-C. Exploring Rural Adolescents’ Dietary Diversity and Its Socioeconomic Correlates: A Cross-Sectional Study from Matlab, Bangladesh. Nutrients 2020, 12, 2230. [Google Scholar] [CrossRef] [PubMed]
  17. Golden, C.D.; Gephart, J.A.; Eurich, J.G.; McCauley, D.J.; Sharp, M.K.; Andrew, N.L.; Seto, K.L. Social-Ecological Traps Link Food Systems to Nutritional Outcomes. Glob. Food Secur. 2021, 30, 100561. [Google Scholar] [CrossRef]
  18. Dornelles, A.Z.; Boonstra, W.J.; Delabre, I.; Denney, J.M.; Nunes, R.J.; Jentsch, A.; Nicholas, K.A.; Schröter, M.; Seppelt, R.; Settele, J.; et al. Transformation Archetypes in Global Food Systems. Sustain. Sci. 2022, 17, 1827–1840. [Google Scholar] [CrossRef]
  19. Mirmiran, P.; Azadbakht, L.; Esmaillzadeh, A.; Azizi, F. Dietary Diversity Score in Adolescents—A Good Indicator of the Nutritional Adequacy of Diets: Tehran Lipid and Glucose Study. Asia Pac. J. Clin. Nutr. 2004, 13, 56–60. [Google Scholar]
  20. Worku, M.; Hailemicael, G.; Asegedech, W. Dietary Diversity Score and Associated Factors among High School Adolescent Girls in Gurage Zone, Southwest Ethiopia. World J. Nutr. Health 2017, 5, 41–45. [Google Scholar] [CrossRef]
  21. Gyimah, L.A.; Annan, R.A.; Apprey, C.; Edusei, A.; Aduku, L.N.E.; Asamoah-Boakye, O.; Azanu, W.; Lutterodt, H. Dietary Diversity and Its Correlates among Pregnant Adolescent Girls in Ghana. PLoS ONE 2021, 16, e0247979. [Google Scholar] [CrossRef]
  22. Nithya, D.J.; Bhavani, R.V. Dietary Diversity and Its Relationship with Nutritional Status among Adolescents and Adults in Rural India. J. Biosoc. Sci. 2018, 50, 397–413. [Google Scholar] [CrossRef] [PubMed]
  23. FAO, FHI 360. Minimum Dietary Diversity for Women: A Guide for Measurement; FAO: Rome, Italy, 2016; ISBN 978-92-5-109153-1. [Google Scholar]
  24. Diop, L.; Becquey, E.; Turowska, Z.; Huybregts, L.; Ruel, M.T.; Gelli, A. Standard Minimum Dietary Diversity Indicators for Women or Infants and Young Children Are Good Predictors of Adequate Micronutrient Intakes in 24–59-Month-Old Children and Their Nonpregnant Nonbreastfeeding Mothers in Rural Burkina Faso. J. Nutr. 2021, 151, 412–422. [Google Scholar] [CrossRef] [PubMed]
  25. Rodríguez-Ramírez, S.; Sánchez-Pimienta, T.G.; Batis, C.; Cediel, G.; Marrón-Ponce, J.A. Minimum Dietary Diversity in Mexico: Establishment of Cutoff Point to Predict Micronutrients Adequacy. Eur. J. Clin. Nutr. 2022, 76, 739–745. [Google Scholar] [CrossRef]
  26. Gómez, G.; Fisberg, R.; Nogueira Previdelli, Á.; Hermes Sales, C.; Kovalskys, I.; Fisberg, M.; Herrera-Cuenca, M.; Cortés Sanabria, L.; García, M.; Pareja Torres, R.; et al. Diet Quality and Diet Diversity in Eight Latin American Countries: Results from the Latin American Study of Nutrition and Health (ELANS). Nutrients 2019, 11, 1605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Steyn, N.; Nel, J.; Nantel, G.; Kennedy, G.; Labadarios, D. Food Variety and Dietary Diversity Scores in Children: Are They Good Indicators of Dietary Adequacy? Public Health Nutr. 2006, 9, 644–650. [Google Scholar] [CrossRef] [Green Version]
  28. Hatløy, A.; Torheim, L.E.; Oshaug, A. Food Variety—A Good Indicator of Nutritional Adequacy of the Diet? A Case Study from an Urban Area in Mali, West Africa. Eur. J. Clin. Nutr. 1998, 52, 891–898. [Google Scholar] [CrossRef] [Green Version]
  29. Sistema de Información Estadística de Derechos de la Niñez y Adolescencia (SIEDNA). Personas Menores de Edad a La Luz Del Censo 2011/Underage Persons in Light of 2011 Census; UCR: San José, Costa Rica, 2013; Available online: https://accionsocial.ucr.ac.cr/sites/default/files/documentos/personas_menores_de_edad_a_la_luz_del_censo_2011.pdf (accessed on 10 October 2022).
  30. Programa Estado de la Nación. Sétimo Informe Estado de La Educación/Seventh State of Education Report; Programa Estado de la Nación: San José, Costa Rica, 2019; Available online: https://estadonacion.or.cr/wp-content/uploads/2019/08/Estado-Educaci%C3%B3n-RESUMEN-2019-WEB.pdf (accessed on 15 October 2022).
  31. Ryan, T.P. Sample Size Determination and Power. In Wiley Series in Probability and Statistics, 1st ed.; Wiley: New York, NY, USA, 2013; ISBN 978-1-118-43760-5. [Google Scholar]
  32. Alam, M. Selection of the Samples with Probability Proportional to Size. Sci. J. Appl. Math. Stat. 2015, 3, 230. [Google Scholar] [CrossRef]
  33. Sánchez-Muñoz, C.; Zabala, M.; Williams, K. Handbook of Anthropometry: Physical Measures of Human Form in Health and Disease; Preedy, V.R., Ed.; Springer: New York, NY, USA, 2012; ISBN 978-1-4419-1787-4. [Google Scholar]
  34. World Health Organization. WHO Child Growth Standards: Length/Height-for-Age, Weight-for-Age, Weight-for-Length, Weight-for-Height and Body Mass Index-for-Age: Methods and Development; De Onis, M., Ed.; WHO child growth standards; WHO Press: Geneva, Switzerland, 2006; ISBN 978-92-4-154693-5. [Google Scholar]
  35. Chinnock, A. Diario de Consumo de Alimentos/Food Consumption Diary; Instrumento para el Registro de Información/Instrument for Information Registration; UCR: San José, Costa Rica, 2007. [Google Scholar]
  36. Trevino, R.P.; Ravelo, A.V.; Senne-Duff, B.; Murad, M.; Diaz, J.F.; Birkenfeld, E. Poor Validity of Dietary Recall in Low-Income Hispanic Children Using Digital Food Imaging Analysis as the Reference. J. Food Nutr. Diet. 2016, 1, 107. [Google Scholar] [CrossRef] [Green Version]
  37. Dutch, D.C.; Golley, R.K.; Johnson, B.J. Diet Quality of Australian Children and Adolescents on Weekdays versus Weekend Days: A Secondary Analysis of the National Nutrition and Physical Activity Survey 2011–2012. Nutrients 2021, 13, 4128. [Google Scholar] [CrossRef]
  38. Monge-Rojas, R.; Vargas-Quesada, R.; Colón-Ramos, U.; Chinnock, A. Dietary Intake and Sources of Added Sugars in Various Food Environments in Costa Rican Adolescents. Nutrients 2022, 14, 959. [Google Scholar] [CrossRef]
  39. Institute of Medicine of the National Academies. Dietary Reference Intakes: The Essential Guide to Nutrient Requirements; National Academies Press: Washington, DC, USA, 2006; p. 11537. ISBN 978-0-309-15742-1. [Google Scholar]
  40. Institute of Medicine of the National Academies. Dietary Reference Intakes: Applications in Dietary Assessment; National Academies Press: Washington, DC, USA, 2000; p. 9956. ISBN 978-0-309-07183-3. [Google Scholar]
  41. Oldewage-Theron, W.H.; Kruger, R. Food Variety and Dietary Diversity as Indicators of the Dietary Adequacy and Health Status of an Elderly Population in Sharpeville, South Africa. J. Nutr. Elder. 2008, 27, 101–133. [Google Scholar] [CrossRef] [PubMed]
  42. Madden, J.P.; Goodman, S.J.; Guthrie, H.A. Validity of the 24-Hr. Recall. Analysis of Data Obtained from Elderly Subjects. J. Am. Diet. Assoc. 1976, 68, 143–147. [Google Scholar] [PubMed]
  43. McCarthy, H.; Dixon, M.; Crabtree, I.; Eaton-Evans, M.J.; McNulty, H. The Development and Evaluation of the Screening Tool for the Assessment of Malnutrition in Paediatrics (STAMP©) for Use by Healthcare Staff: Development and Evaluation of STAMP©. J. Hum. Nutr. Diet. 2012, 25, 311–318. [Google Scholar] [CrossRef]
  44. ValorNut UCR: Programa Para El Cálculo de Valor Nutritivo de Los Alimentos. Available online: https://www.nutricion.ucr.ac.cr/index.php/es/valornut (accessed on 2 October 2022).
  45. Arce, R. Lo Rural es Diverso: Evidencia para el Caso de Costa Rica; Rodríguez, A., Saborío, M., Eds.; Instituto Interamericano de Cooperación para la Agricultura (IICA): San José, Costa Rica, 2008; ISBN 978-92-9039-855-4. [Google Scholar]
  46. Valenciano Salazar, J.A.; Pagani Centeno, L.; Álvarez Madrigal, N. Seguridad y Soberanía Alimentaria en Costa Rica; Editorial UNA: Heredia, Costa Rica, 2020; ISBN 978-9930-588-01-7. [Google Scholar]
  47. Singh, N. Pulses: An Overview. J. Food Sci. 2017, 54, 853–857. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Ministerio de Salud de Costa Rica. Reforma Reglamento para el Enriquecimiento del Arroz. Decreto Ejecutivo N° 39985-S 2016. Available online: https://www.fao.org/faolex/results/details/es/c/LEX-FAOC162902 (accessed on 1 October 2022).
  49. Research Institute (IFPRI), I.F.P. Urbanization and the Nutrition Transition; International Food Policy Research Institute: Washington, DC, USA, 2017. [Google Scholar]
  50. Barría, R.M.; Amigo, H. Nutrition transition: A review of Latin American profile. Arch. Latinoam. Nutr. 2006, 56, 3–11. [Google Scholar] [PubMed]
  51. Popkin, B.M. Urbanization, Lifestyle Changes and the Nutrition Transition. World Dev. 1999, 27, 1905–1916. [Google Scholar] [CrossRef]
  52. Doak, C.M.; Adair, L.S.; Bentley, M.; Monteiro, C.; Popkin, B.M. The Dual Burden Household and the Nutrition Transition Paradox. Int. J. Obes. 2005, 29, 129–136. [Google Scholar] [CrossRef] [Green Version]
  53. Popkin, B.M. Nutritional Patterns and Transitions. Popul. Dev. Rev. 1993, 19, 138. [Google Scholar] [CrossRef]
  54. Willett, W. Nutritional Epidemiology. In Monographs in Epidemiology and Biostatistics, 3rd ed.; Oxford University Press: Oxford, UK; New York, NY, USA, 2013; ISBN 978-0-19-975403-8. [Google Scholar]
Figure 1. Sensitivity and specificity of different DDS cut-off points with MAR cut-off points ranging from 0.65 to 0.75. 1 Specificity: identifies nutritionally appropriate diets as adequate. 2 Sensitivity: identifies nutritionally inappropriate diets as inadequate. DDS: Diet Diversity Score; MAR: Mean Adequacy Ratio.
Figure 1. Sensitivity and specificity of different DDS cut-off points with MAR cut-off points ranging from 0.65 to 0.75. 1 Specificity: identifies nutritionally appropriate diets as adequate. 2 Sensitivity: identifies nutritionally inappropriate diets as inadequate. DDS: Diet Diversity Score; MAR: Mean Adequacy Ratio.
Nutrients 14 05093 g001
Figure 2. Proportion of participants consuming at least 15 g of each food group over one-day food record, according to residence area.
Figure 2. Proportion of participants consuming at least 15 g of each food group over one-day food record, according to residence area.
Nutrients 14 05093 g002
Table 1. General characteristics of Costa Rican adolescents, according to residence area.
Table 1. General characteristics of Costa Rican adolescents, according to residence area.
CharacteristicsOverall 1 (n = 818)Residence Area
Urban (n = 411)Rural (n = 407)p-Value 2
 Age (y)15.0 ± 1.714.9 ± 1.715.1 ± 1.70.208
Sex
 Female520 (63.6)259 (49.8) a261 (50.2) a0.741
 Male298 (36.4)152 (51.0) a146 (49.0) a
Socioeconomic status
 Low263 (32.2)87 (33.1) a179 (66.9) b<0.0001
 Middle325 (39.7)158 (48.6) a167 (51.4) a
 High230 (28.1)166 (72.2) a64 (27.8) b
Nutritional status
 Non-overweight551 (67.4)279 (50.6) a272 (49.4) a0.748
 Overweight/obesity267 (32.6)132 (49.4) a135 (50.6) a
1 Values are means ± SD or frequencies (%) unless otherwise indicated. a,b Labeled frequencies in a row without a common letter differ (p < 0.05). 2 p-values < 0.05 are statistically significant and were determined using the Wilcoxon test, or the chi-square test with Bonferroni post hoc test when required. p-values compare either means of age or proportions in categories of sex, socioeconomic status, and nutritional status between categories of residence area (urban/rural).
Table 2. Diet Diversity Score (DDS) and proportion of participants reaching the minimum DDS 1 (n = 818).
Table 2. Diet Diversity Score (DDS) and proportion of participants reaching the minimum DDS 1 (n = 818).
CharacteristicDiet Diversity Score (DDS)Participants Reaching the Minimum DDS
MeanSDp-Value 2n (%)p-Value 2
 Overall4.171.43 544 (66.5)
Sex
 Female4.181.430.679349 (67.1)0.624
 Male4.141.44 195 (65.4)
Residence area
 Urban4.001.420.001251 (61.1)0.001
 Rural4.331.43 293 (72.0)
Socioeconomic status
 Low4.281.440.303183 (69.6)0.420
 Middle4.111.39 213 (65.5)
 High4.111.48 148 (64.4)
Nutritional status
 Non overweight4.151.430.680367 (66.6)0.929
 Overweight/obesity4.191.43 177 (66.3)
1 Consumption of at least 4 out of 10 food groups. 2 p-values < 0.05 are statistically significant and were determined using the chi-square, Student’s t, or ANOVA tests.
Table 3. Food group consumption of Costa Rican adolescents, according to Diet Diversity Score (DDS) and residence area.
Table 3. Food group consumption of Costa Rican adolescents, according to Diet Diversity Score (DDS) and residence area.
Food Group 1 (g/d)Diet Diversity Residence Area
DDS < 4 (n = 274)DDS ≥ 4 (n = 544)p-Value 2Urban
(n = 411)
Rural
(n = 407)
p-Value 2
Starchy staples357.1 ± 184.3387.7 ± 197.90.031353.1 ± 190.0402.1 ± 196.8<0.0001
Milk and milk products91.1 ± 204.0166.5 ± 207.1<0.0001155.1 ± 220.7127.4 ± 195.80.010
Pulses55.1 ± 108.2106.9 ± 112.2<0.000174.1 ± 105.9105.1 ± 118.8<0.0001
Flesh foods56.5 ± 84.7102.3 ± 94.5<0.000191.3 ± 98.682.6 ± 88.60.245
Other fruits25.4 ± 78.191.8 ± 145.0<0.000173.4 ± 145.565.7 ± 109.20.816
Other vegetables11.8 ± 36.362.3 ± 86.8<0.000135.8 ± 69.455.1 ± 84.0<0.0001
Eggs7.0 ± 21.923.4 ± 37.9<0.000115.9 ± 32.720.1 ± 35.70.029
Other vitamin A-rich F&V5.9 ± 45.615.0 ± 57.3<0.000114.3 ± 64.19.7 ± 40.90.405
Dark green leafy vegetables0.9 ± 14.56.7 ± 31.4<0.00012.9 ± 21.36.6 ± 31.80.016
Nuts and seeds0.3 ± 3.42.8 ± 20.00.0031.8 ± 14.22.1 ± 18.50.467
1 Values are means ± SD. 2 p-values < 0.05 are statistically significant and were determined using the Wilcoxon test. DDS: Diet Diversity Score; F&V: Fruits and vegetables.
Table 4. Bivariate and multivariate logistic regression analysis: Sociodemographic variables related to having a diverse diet, according to Diet Diversity Score (DDS) (n = 818).
Table 4. Bivariate and multivariate logistic regression analysis: Sociodemographic variables related to having a diverse diet, according to Diet Diversity Score (DDS) (n = 818).
VariableBivariate AnalysisMultivariate Analysis 1
OR95% CIp-ValueAdj OR95% CIp-Value
Sex
 Women1.080.80–1.460.6241.080.80–1.460.622
Age1.091.00–1.190.0471.080.99–1.180.067
Residence area
 Rural1.641.22–2.200.0011.621.19–2.200.002
Socioeconomic status
 Middle0.830.59–1.180.2990.910.64–1.300.617
 High0.790.54–1.150.2170.980.66–1.460.922
Nutritional status
 Overweight/obesity0.990.72–1.340.9290.990.73–1.360.989
Constant---0.460.12–1.840.275
1 Overall model test p-value = 0.020. Goodness-of-fit test p-value = 0.470. Correct predictions: 66.5%.
Table 5. Mean Nutrient Adequacy Ratio (NAR), according to Diet Diversity Score (DDS).
Table 5. Mean Nutrient Adequacy Ratio (NAR), according to Diet Diversity Score (DDS).
NutrientsOverall (n = 818)DDS < 4 (n = 274)DDS ≥ 4 (n = 544) r1p-Value
MeanSDMeanSDMeanSDp-Value 2
Calcium0.4920.2920.3680.2710.5550.283<0.00010.383<0.0001
Iron0.9680.1310.9380.1900.9830.0840.0010.171<0.0001
Zinc0.7880.2490.6660.2840.8490.204<0.00010.354<0.0001
Vitamin C0.8850.2680.7780.3570.9390.187<0.00010.265<0.0001
Thiamin0.9650.1320.9230.1970.9850.073<0.00010.205<0.0001
Riboflavin0.9400.1730.8690.2480.9760.102<0.00010.276<0.0001
Niacin0.9540.1510.9000.2230.9810.084<0.00010.232<0.0001
Vitamin B60.8980.2130.8020.2840.9470.143<0.00010.311<0.0001
Folate equivalents0.9690.1280.9340.1890.9870.075<0.00010.163<0.0001
Cobalamin0.9570.1590.9080.2340.9820.093<0.00010.215<0.0001
Vitamin A0.6860.3100.5780.3380.7400.280<0.00010.309<0.0001
MAR0.8640.1350.7880.1760.9020.087<0.00010.433<0.0001
1 Spearman’s rank correlation coefficients (r) were calculated between each NAR/MAR value and the DDS for the whole sample. 2 p-values < 0.05 are statistically significant and were determined using the Wilcoxon test. NAR: Nutrient Adequacy Ratio; MAR: Mean Adequacy Ratio; DDS: Diet Diversity Score.
Table 6. Proportion of participants reaching Nutrient Adequacy Ratio (NAR) ≥ 0.70, according to Diet Diversity Score (DDS) and residence area.
Table 6. Proportion of participants reaching Nutrient Adequacy Ratio (NAR) ≥ 0.70, according to Diet Diversity Score (DDS) and residence area.
NutrientOverall 1 (n = 818)DDS < 4 (n = 274)DDS ≥ 4 (n = 544) Urban (n = 411)Rural (n = 407)
n (%)n (%)n (%)p-Value 2n (%)n (%)p-Value 2
Calcium214 (26.2)38 (13.9)176 (32.4)<0.0001123 (29.9)91 (22.4)0.014
Iron776 (94.9)249 (90.9)527 (96.9)<0.0001386 (93.9)390 (95.8)0.217
Zinc534 (65.3)127 (46.4)407 (74.8)<0.0001253 (61.6)281 (69.0)0.025
Vitamin C690 (84.4)195 (71.2)495 (91.0)<0.0001336 (81.8)354 (87.0)0.040
Thiamin774 (94.6)240 (87.6)534 (98.2)<0.0001386 (93.9)388 (95.3)0.370
Riboflavin743 (90.8)218 (79.6)525 (96.5)<0.0001368 (89.5)375 (92.1)0.198
Niacin752 (91.9)228 (83.2)524 (96.3)<0.0001378 (92.0)374 (92.0)0.967
Vitamin B6680 (83.1)184 (67.2)496 (91.2)<0.0001326 (79.3)354 (87.0)0.003
Folate equivalents779 (95.2)245 (89.4)534 (98.2)<0.0001389 (94.7)390 (95.8)0.430
Cobalamin762 (93.2)236 (86.1)526 (96.7)<0.0001386 (93.9)376 (92.4)0.385
Vitamin A430 (52.6)110 (40.2)320 (58.8)<0.0001215 (52.1)215 (52.8)0.883
MAR736 (90.0)210 (76.6)526 (96.7)<0.0001368 (89.5)368 (90.4)0.675
1 Values are frequencies (%) unless otherwise indicated. 2 p-values < 0.05 are statistically significant and were determined using the chi-square test. NAR: Nutrient Adequacy Ratio; MAR: Mean Adequacy Ratio; DDS: Diet Diversity Score.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Monge-Rojas, R.; Vargas-Quesada, R.; Gómez, G. Role of Residence Area on Diet Diversity and Micronutrient Intake Adequacy in Urban and Rural Costa Rican Adolescents. Nutrients 2022, 14, 5093. https://doi.org/10.3390/nu14235093

AMA Style

Monge-Rojas R, Vargas-Quesada R, Gómez G. Role of Residence Area on Diet Diversity and Micronutrient Intake Adequacy in Urban and Rural Costa Rican Adolescents. Nutrients. 2022; 14(23):5093. https://doi.org/10.3390/nu14235093

Chicago/Turabian Style

Monge-Rojas, Rafael, Rulamán Vargas-Quesada, and Georgina Gómez. 2022. "Role of Residence Area on Diet Diversity and Micronutrient Intake Adequacy in Urban and Rural Costa Rican Adolescents" Nutrients 14, no. 23: 5093. https://doi.org/10.3390/nu14235093

APA Style

Monge-Rojas, R., Vargas-Quesada, R., & Gómez, G. (2022). Role of Residence Area on Diet Diversity and Micronutrient Intake Adequacy in Urban and Rural Costa Rican Adolescents. Nutrients, 14(23), 5093. https://doi.org/10.3390/nu14235093

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop