Vitamin D deficiency (VDD) is a widespread public health problem, with substantial regional and seasonal variability. The effects of VDD may extend beyond the musculoskeletal system; thus, it is important to characterize the extent of the problem and its determinants at the population level. Few studies have investigated the prevalence of VDD in Mesoamerica. In a study of Guatemalan adolescents, 12·8 % had VDD, defined as serum 25-hydroxyvitamin D (25(OH)D) concentration <50 nmol/l( Reference Naqvi, Solomons and Campos 1 ). In clinic-based studies of adults in Puerto Rico( Reference Diaz, Gonzalez and Ramos-Trautmann 2 , Reference Gonzalez, Ramos-Trautmann and Diaz-Luquis 3 ) and among older adults from Mexico( Reference Carrillo-Vega, Garcia-Pena and Gutierrez-Robledo 4 ), the prevalence of VDD was 36·0 and 36·9 %, respectively.
The main source of vitamin D is skin synthesis from exposure to sunlight, although it can also be obtained from some foods including fatty fish, mushrooms, egg yolks and fortified milk. Although skin exposure and time spent outdoors have been investigated as correlates of vitamin D status( Reference Naqvi, Solomons and Campos 1 ), the relative contributions of both sunlight and diet to the vitamin D status of Mesoamericans have not been examined. The demographic and socio-economic patterning of vitamin D status in the region is also unknown. Identifying correlates of serum 25(OH)D concentrations in adults and children across the Mesoamerican region is an important initial step to determine which specific subgroups may benefit from public health interventions to prevent and decrease the burden of VDD.
The vitamin D status of family members may be correlated due to shared habits or genetics( Reference Orton, Morris and Herrera 5 , 6 ). Nevertheless, the level of aggregation of this measure between parents and children has not been estimated. It is relevant to understand whether shared environmental or genetic factors may influence vitamin D status because interventions to improve it might be more effective if they are targeted to families rather than individuals.
The aim of the present study was to examine the associations of sociodemographic characteristics, dietary habits and time spent on outdoor activities with serum 25(OH)D concentrations and the aggregation of 25(OH)D concentrations among mothers, fathers and their children in nine Mesoamerican countries.
Methods
Study population and field methods
The present study was conducted in the context of the Nine Mesoamerican Countries Metabolic Syndrome Study (NiMeCoMeS), a cross-sectional survey investigating dietary and other environmental correlates of metabolic syndrome among school-aged children and their two biological parents. Details on the study design have been described previously( Reference Villamor, Finan and Ramirez-Zea 7 ). Briefly, between July 2011 and November 2013, we identified all public primary schools in peri-urban areas of the capital cities of Guatemala, El Salvador, the Dominican Republic, Honduras, Nicaragua, Costa Rica, Panama and Belize, and from the city of Tuxtla Gutiérrez in Chiapas, Mexico, by contacting each country’s Ministry of Education. We considered schools with at least 400 students of both sexes and grouped them by geographic areas roughly representing the four cardinal directions. Next, one to three schools were randomly selected from each geographic area in each city to reach the target enrolment of thirty families per country, a sample size determined by the availability of funding. The study team identified students aged 7–12 years as potentially eligible for recruitment using the list of enrolled students in each school and invited their parents to attend a meeting. The study aims and procedures were explained to the families in attendance and additional eligibility criteria were confirmed. These criteria included living with both biological parents, not being pregnant or having a pregnant mother, and not having a sibling already identified for potential participation in the study. Researchers provided consent forms to families at the end of the meeting and afforded them one week to consider participation. Study personnel collected the signed consent forms at schools and verbally confirmed assent to participate among children of consenting parents. The final sample included 267 families (Guatemala, n 31; El Salvador, n 30; Dominican Republic, n 30; Honduras, n 30; Nicaragua, n 31; Costa Rica, n 27; Panama, n 26; Belize, 31; and Mexico, n 31). The sampling strategy employed did not intend to produce a representative sample of the population in these countries.
Sociodemographic characteristics, food intake and health habits of children and their parents were determined through questionnaires administered during a research visit to the participant’s home or at a health centre. The questionnaires inquired on parental age, education level, maternal parity, and socio-economic status indicators including home and household assets ownership. Household food insecurity was assessed with the use of the Latin American and Caribbean Food Security Scale (ELCSA)( Reference Perez-Escamilla, Melgar-Quiñonez and Nord 8 ), a survey validated for the region. We evaluated participants’ usual diet over the past 12 months using a semi-quantitative ninety-seven-item FFQ that was administered separately to mothers, fathers and children. The FFQ had been previously validated in Costa Rican adults( Reference Kabagambe, Baylin and Allan 9 , Reference Campos, Willett and Peterson 10 ) and included major sources of energy and micronutrients in food groups such as fruits, vegetables, animal protein (dairy, eggs and meat), carbohydrates (breads, flours, and cereals), beverages (dairy, sugar-sweetened beverages, coffee and alcohol), fast foods, desserts, sauces and dressings, cooking fats and oils, and dietary supplements. For each item, research assistants described reference portion sizes in natural units or standard measures for commonly consumed servings among this population. There were nine frequency of intake options: ≥6 times/d, 4–5 times/d, 2–3 times/d, 1 time/d, 5–6 times/week, 2–4 times/week, 1 time/week, 1–3 times/month or <1 time/month.
The questionnaires also included questions regarding the time spent on activities typically conducted outdoors as an indicator of sun exposure. For adults, these questions involved the usual frequency of gardening, performing agricultural work or walking outdoors for 1 h during the previous 12 months. Responses were coded into one of ten frequency options: ≥6 times/d, 4–5 times/d, 2–3 times/d, 1 time/d, 5–6 times/week, 2–4 times/week, 1 time/week, 1–3 times/month, <1 time/month or never. For children, mothers reported on their child’s average daily time spent playing outdoors on weekdays and weekends.
Trained research assistants obtained anthropometric measures from each participant. Weight was measured in light clothing to the nearest 0·1 kg on digital scales (Tanita, Arlington Heights, IL, USA) and height was measured without shoes to the nearest 1 mm with the use of portable stadiometers (Seca, Hanover, MD, USA). At the end of the visit, a fasting blood sample was collected by venepuncture from each family member.
Laboratory methods
Blood samples were transported on ice to each country’s collaborating institution where serum was separated and cryopreserved at −20°C until laboratory analysis. Analysis of 25(OH)D was performed at Heartland Assays, Inc. (Ames, IA, USA) with the DiaSorin LIAISON 25-OH Vitamin D Total assay kit, a direct, competitive chemiluminescence immunoassay. This validated method has a sensitivity of 6·25 nmol/l, and an inter- and intra-assay CV of 11·2 and 8·1 %, respectively( Reference Horst 11 ).
Definition of outcomes
The main outcome of interest was mean serum 25(OH)D concentration. A secondary outcome was the prevalence of serum 25(OH)D <50 nmol/l, a cut-off point recommended by the US Endocrine Society to define deficiency( Reference Holick, Binkley and Bischoff-Ferrari 12 ). A lower cut-off point of 30 nmol/l, proposed to define deficiency by the US Institute of Medicine( 13 ), was not considered because only one person had 25(OH)D concentration under this threshold.
Definition of exposures
Adults
Education level was categorized according to the number of completed years of schooling (incomplete elementary, 1–5 years; complete elementary, 6 years; incomplete secondary, 7–11 years; complete secondary, 12 years; or post-secondary, ≥13 years). We categorized height into quartiles separately for men and women. BMI was calculated as kg/m2 and categorized as <25·0, 25·0–29·9, 30·0–34·9 or ≥35·0( 14 ). A category of <18·5 kg/m2 was not considered because only six people were under this cut-off point. The number of home assets was determined as the total number of affirmative responses to owning a car, bicycle, refrigerator, gas stove, electric stove, blender, microwave, washing machine, colour television, stereo or computer, or having in-house access to the Internet. Household food insecurity was categorized according to the number of affirmative responses to the questions in the ELCSA survey (no insecurity, 0; mild insecurity, 1–5; moderate insecurity, 6–10; severe insecurity, ≥11). From the FFQ, we selected foods with a high vitamin D content as dietary exposures according to food composition tables. These included dairy, eggs, and canned tuna or sardines. Dairy intake frequency was calculated from the sum of weighted intake frequencies of milk, cheese, cream and ice cream. Egg intake represented the sum of weighted intake frequencies of one hard boiled or fried egg. Supplement use was not considered as a dietary correlate of vitamin D because fewer than five participants reported use of supplements that may contain vitamin D. Covariates were categorized as presented in Tables 1 and 2.
Serum 25-hydroxyvitamin D (25(OH)D) concentrations according to sociodemographic characteristics in adults from nine Mesoamerican countries, 2011–2013
* Totals may be less than 492 because of missing values.
† Test for linear trend from linear regression models with serum 25(OH)D as the continuous outcome and a variable representing ordinal categories of each characteristic as the continuous predictor. χ 2 score test for sex, country and month of measurement. An exchangeable correlation structure was specified in all models to account for intrafamily clustering.
‡ Test for linear trend from Poisson regression models with serum 25(OH)D <50 nmol/l as the dichotomous outcome and a variable representing ordinal categories of each characteristic as the continuous predictor. χ 2 score test for sex and month of measurement. An exchangeable correlation structure was specified in all models to account for intrafamily clustering. Fisher’s exact test for country.
§ From a list that included car, bicycle, refrigerator, gas stove, electric stove, blender, microwave, washing machine, colour television, stereo, computer, and having in-home access to the Internet.
Serum 25-hydroxyvitamin D (25(OH)D) concentrations according to food intake frequency and time spent on outdoor activities in adults from nine Mesoamerican countries, 2011–2013
* Totals may be less than 492 because of missing values·
† Test for linear trend from linear regression models with serum 25(OH)D as the continuous outcome and a variable representing ordinal categories of each characteristic as the continuous predictor. χ 2 score test for agricultural work. An exchangeable correlation structure was specified in all models to account for intrafamily clustering.
‡ Test for linear trend from Poisson regression models with serum 25(OH)D <50 nmol/l as the dichotomous outcome and a variable representing ordinal categories of each characteristic as the continuous predictor. χ 2 score test for agricultural work. An exchangeable correlation structure was specified in all models to account for intrafamily clustering.
§ Sum of weighted intake frequencies of milk (200 g), cheese (30 g), cream (12 g) and ice cream (60 g).
║ Sum of weighted intake frequencies of one hard boiled or fried egg (60 g).
¶ Intake of 60 g.
Children
Children’s height-for-age Z-score and BMI-for-age Z-score were calculated using sex-specific growth references from the WHO( Reference de Onis 15 ). The parent who had a higher level of educational attainment determined the highest level of parental education. Maternal height in centimetres was categorized into quartiles and maternal BMI was categorized as described for adults. Mother’s and father’s serum 25(OH)D concentrations were categorized as <50, ≥50 to <75 or ≥75 nmol/l. Covariates were categorized as presented in Tables 4 and 5.
Data analysis
Adults and children were analysed separately. Samples of children from Panama were not available, thus the analytic sample comprised 492 adults and 223 children.
We first compared mean and sd serum 25(OH)D concentrations by categories of sociodemographic characteristics, food intake frequency and time spent on outdoor activities. For ordinal predictors we conducted tests for linear trend using linear regression models, with 25(OH)D concentration as the continuous outcome and a variable representing ordinal categories of the predictor as a continuous covariate. For nominal predictors, we used χ 2 score statistics. Next, we estimated the prevalence of serum 25(OH)D <50 nmol/l by levels of sociodemographic, diet and outdoor activity characteristics. For ordinal predictors in adults, we performed tests for linear trend with the use of Poisson regression models, with serum 25(OH)D <50 nmol/l as the dichotomous outcome and a variable representing ordinal categories of the predictor as a continuous covariate. An exchangeable correlation structure was specified in all models of adults to account for intrafamily clustering. For nominal predictors, we used the χ 2 score test. In children, we used the Cochrane–Armitage exact test for trend for ordinal predictors. For nominal predictors, we used the Fisher’s exact test.
We estimated adjusted differences with 95 % CI in mean serum 25(OH)D concentrations by categories of exposure variables with the use of multivariable linear regression. We included predictors that were significantly associated with the outcome in bivariate analysis at P<0·10; predictors that remained significant (P<0·05) in multivariable analysis were retained in the final models. Robust estimates of the variance were calculated in all models to overcome potential deviations from the multivariate normality assumption( Reference White 16 ).
Family aggregation
We estimated Pearson correlations with 95 % CI of serum 25(OH)D concentrations between mothers and fathers, mothers and children, and fathers and children. Next, we used linear random-effects models with serum 25(OH)D as the continuous outcome and a random effect for family to partition the variance into within-family and between-family components. We also considered fixed effects for age and sex to account for extraneous between-person variation. The intraclass correlation coefficient (ICC) was calculated as the ratio of the between-family variance to the sum of the between- and within-family variances. 95 % CI were obtained using the delta method( Reference Hankinson, Manson and Spiegelman 17 ). Higher values for an ICC imply that people from the same family would have more similar serum 25(OH)D concentrations than people from different families. We then estimated the ICC among mothers and fathers, mothers and children, and fathers and children, and compared them with the ICC for the whole family to examine which combinations of family membership were associated with greater aggregation. Higher mother–father than mother–child or father–child aggregation could indicate a greater role for shared environment rather than shared genetics in family clustering of 25(OH)D since mothers and fathers are typically not genetically related. In supplemental analysis, we considered whether potential clustering by country could affect the results on family aggregation.
Results
The mean (sd) age of mothers, fathers and children at enrolment was 37·0 (6·4), 40·6 (8·3) and 10·0 (1·7) years, respectively; 53·4 % of children were girls. Sociodemographic characteristics, food intake frequencies and time spent on outdoor activities varied by country (see online supplementary material, Supplemental Table 1).
Vitamin D correlates in adults
Among adults, mean (sd) serum 25(OH)D concentration was 81·3 (21·1) nmol/l; 82·5 (21·8) nmol/l in mothers and 80·1 (20·4) nmol/l in fathers. The prevalence of serum 25(OH)D <30, <50 and <75 nmol/l was, respectively, 0·2, 3·9 and 39·8 %. Mean serum 25(OH)D concentrations were inversely related to education level and number of home assets. Participants from Nicaragua and those with measurements obtained in May or June had the highest concentrations (Table 1). The prevalence of VDD was higher in adults with more household assets compared with that in adults with fewer home assets. Prevalence of VDD in adults varied across countries, ranging from 0·0 % in Mexico to 9·4 % in Costa Rica. Mean serum 25(OH)D concentrations were positively associated with time spent gardening (Table 2). In multivariable analysis, mean serum 25(OH)D concentrations varied significantly by country (P<0·0001); the highest concentrations were in Nicaragua, Mexico and Belize whereas the lowest were in Panama, Guatemala and the Dominican Republic (Table 3). Compared with people who never gardened, adjusted mean serum 25(OH)D concentration was 4·8 nmol/l higher (95 % CI 0·0, 9·6 nmol/l) in people who gardened for 1–6 h/week and 6·4 nmol/l higher (95 % CI −0·6, 13·3 nmol/l) in people who gardened for ≥1 h/d (P-trend=0·03). Month of measurement was not associated with mean serum 25(OH)D concentrations in multivariable analysis.
Multivariable-adjusted correlates of serum 25-hydroxyvitamin D (25(OH)D) concentrations in adults from nine Mesoamerican countries, 2011–2013
* From multivariable linear regression with serum 25(OH)D as the continuous outcome. Predictors included country and gardening. An exchangeable correlation structure was specified to account for intrafamily clustering.
† For country, χ 2 score test. For gardening, test for linear trend when a variable representing ordinal categories of gardening time was included as a continuous predictor.
Vitamin D correlates in children
Among children, mean (sd) serum 25(OH)D concentration was 79·5 (18·1) nmol/l. The prevalence of serum 25(OH)D <30, <50 and <75 nmol/l was, respectively, 0·0, 3·6 and 44·4 %. Mean serum 25(OH)D concentrations were positively related to mother’s height, parental 25(OH)D concentrations and Nicaraguan origin (Table 4). The prevalence of VDD was higher in children whose mothers had higher BMI compared with children whose mothers had lower BMI. In addition, the prevalence of VDD was higher in children whose mothers had serum 25(OH)D concentration of <50 or 50–<75 nmol/l compared with children whose mothers’ serum 25(OH)D concentration was ≥75 nmol/l. Prevalence of VDD ranged from 0·0 % in Guatemala, the Dominican Republic, Honduras and Mexico to 11·5 % in Costa Rica. Mean serum 25(OH)D concentrations were positively associated with frequency of dairy intake (Table 5). In multivariable analysis, mean serum 25(OH)D concentrations were positively associated with male sex, mother’s 25(OH)D concentrations and dairy intake frequency (Table 6). Every 1 nmol/l of maternal serum 25(OH)D was related to an adjusted 0·3 nmol/l higher 25(OH)D concentration in the children (95 % CI 0·2, 0·5 nmol/l; P<0·0001). Mean 25(OH)D concentration among children with dairy intake ≥2·5 times/d was 6·8 nmol/l higher than that of children with intake <5 times/week (95 % CI 0·7, 12·9 nmol/l; P=0·03). Finally, mean serum 25(OH)D concentrations were inversely related to mother’s BMI and to having >4 home assets. Maternal height and country of origin were not associated with the children’s serum 25(OH)D concentrations in multivariable analysis.
Serum 25-hydroxyvitamin D (25(OH)D) concentrations according to sociodemographic characteristics in children from eight Mesoamerican countries, 2011–2013
* Totals may be less than 223 because of missing values.
† Test for trend for ordinal variables. χ 2 score test for sex, country and month of measurement.
‡ Cochrane–Armitage exact test for trend for ordinal variables. Fisher’s exact test for sex, country and month of measurement.
§ According to the WHO reference( Reference de Onis 15 ).
║ From a list that included car, bicycle, refrigerator, gas stove, electric stove, blender, microwave, washing machine, colour television, stereo, computer, and having in-home access to the Internet.
Serum 25-hydroxyvitamin D (25(OH)D) concentrations according to food intake frequency and time spent on outdoor activities in children from eight Mesoamerican countries, 2011–2013
* Totals may be less than 223 because of missing values.
† Test for trend for ordinal variables.
‡ Cochrane–Armitage exact test for trend for ordinal variables.
§ Sum of weighted intake frequencies of milk (200 g), cheese (20 g), cream (6 g) and ice cream (60 g).
║ Sum of weighted intake frequencies of one hard boiled or fried egg (60 g).
¶ Intake of 25 g.
Multivariable-adjusted correlates of serum 25-hydroxyvitamin D (25(OH)D) concentrations in children from eight Mesoamerican countries, 2011–2013
* From a multivariable linear regression model with serum 25(OH)D as the continuous outcome. Predictors included all variables presented. Empirical variances were specified in the model.
† For height and frequency of dairy intake, test for linear trend when a variable representing ordinal categories of the characteristic was introduced as a continuous predictor. Wald test for sex, mother’s BMI, mother’s serum 25(OH)D concentrations and home assets. χ 2 score test for country.
‡ From a list that included car, bicycle, refrigerator, gas stove, electric stove, blender, microwave, washing machine, colour television, stereo, computer, and having in-home access to the Internet.
§ Sum of weighted intake frequencies of milk (200 g), cheese (20 g), cream (6 g) and ice cream (60 g).
Family aggregation
Serum 25(OH)D concentrations were positively correlated between family members (Table 7). The mother–child correlation (r=0·49, P<0·0001) was higher than the father–child (r=0·24, P=0·0004) or mother–father (r=0·23, P=0·0004) correlation. Thirty-one per cent of the total variance in serum 25(OH)D was explained by family membership. Adjustment for age and sex did not change the results (adjusted ICC=0·32; 95 % CI 0·24, 0·41). Family membership explained more variability in 25(OH)D concentration than country membership, which only explained 12 %. Mother–child aggregation of serum 25(OH)D concentrations was higher than overall family, father–child or mother–father aggregation.
Family aggregation of serum 25-hydroxyvitamin D (25(OH)D) in eight Mesoamerican countries, 2011–2013
* Pearson correlation coefficient.
† Variance partitioned from a linear random-effects model with serum 25(OH)D as the continuous outcome and a random effect for family. Separate models were fitted to partition the variance for each combination of two family members.
‡ Intraclass correlation coefficient (ICC)=
$$\sigma _{{{\rm between{\hbox-}family}}}^{{\rm 2}} \,/\,\left( {\sigma _{{{\rm between{\hbox-}family}}}^{2} {\plus}\sigma _{{{\rm within{\hbox-}family}}}^{2} } \right).$$
. 95 % CI were calculated using the delta method(
Reference Hankinson, Manson and Spiegelman
17
).
Discussion
In the current study of Mesoamerican families, the overall prevalence of serum 25(OH)D concentration <50 nmol/l was relatively low for both adults (3·9 %) and children (3·6 %). Only one person had a serum 25(OH)D concentration <30 nmol/l, the threshold for deficiency proposed by the US Institute of Medicine( 13 ). Although the burden of VDD in Latin America is not well documented, the prevalence of VDD in our study is among the lowest reported in the region. For example, the prevalence of VDD in the National Survey of Nutrition and Health in Mexican adults was 9·8 % in 2006( Reference Flores, Sanchez Romero and Macias 18 ). No other country in the Mesoamerican region has conducted a nationwide survey to estimate vitamin D serostatus. Studies of Puerto Rican adults with different underlying health conditions found higher prevalence of VDD, ranging from 14 to 43 %( Reference Diaz, Gonzalez and Ramos-Trautmann 2 , Reference Gonzalez, Ramos-Trautmann and Diaz-Luquis 3 , Reference Hernandez Davila, Rodriguez Rivera and Lopez Valentin 19 – Reference Suarez-Martinez, Perez and Cruz 21 ). Nevertheless, neither these investigations nor our study were representative and the results may not be directly comparable.
Prevalence estimates of VDD among children and adolescents in the Mesoamerican region are variable. In asthmatic school-aged children from Costa Rica( Reference Brehm, Celedon and Soto-Quiros 22 ) and among Afro-Caribe Guatemalan adolescents( Reference Naqvi, Solomons and Campos 1 ), the prevalence of VDD was 3·4 and 5·0 %, respectively. These estimates are comparable to ours; however, among children from Southern Mexico (5–11 years old) in 2012 and adolescents (13–19 years old) in 2006 the prevalence of VDD was 29·0 %( Reference Flores, Flores and Macias 23 ) and 11·5 %( Reference Flores, Sanchez Romero and Macias 18 ), respectively. These values are substantially higher than the estimates in our study. Similarly, VDD was much more prevalent in indigenous Mayan adolescents living near the coast (21 %)( Reference Naqvi, Solomons and Campos 1 ) than in our study from Mesoamerica. These differences suggest that there is high within- and between-country variation in vitamin D status in children from the Mesoamerican region.
In our investigation vitamin D status varied substantially by country, although the country-specific sample sizes were too small to allow for robust statistical inference. Adults and children from Nicaragua had the highest mean serum 25(OH)D concentrations. Reasons to explain this difference are speculative but could be related to socio-economic or dietary differences since Nicaraguans in the present study had fewer home assets and higher dairy intake than participants from other countries and these variables were associated with higher 25(OH)D concentrations.
We found that mean serum 25(OH)D concentrations in adults were positively related to time spent gardening but not to other outdoor activities. Time spent on different outdoor activities may exert distinct effects on serum 25(OH)D levels because they could vary in the factors that ultimately influence the amount of vitamin D produced in the skin, including the amount of skin exposed, the time of day of sun exposure, the use of sunscreen and the absolute amount of time exposed to the sun. For example, among elderly adults in Europe 25(OH)D concentrations were positively associated with gardening and cycling but not with walking outdoors( Reference van den Heuvel, van Schoor and de Jongh 24 ). The body postures during gardening may expose a larger skin surface to the sun compared with walking. In addition, the majority of ambient UV radiation occurs between 10.00 and 16.00 hours( Reference Holick 25 ). Perhaps people tend to garden within these hours, whereas other outdoor activities may occur outside this interval. In a study of Puerto Rican medical residents, time spent outdoors between 10.00 and 16.00 hours was positively associated with serum 25(OH)D( Reference Hernandez Davila, Rodriguez Rivera and Lopez Valentin 19 ).
We did not find associations between age, sex or dietary sources of vitamin D and serum 25(OH)D concentrations among adults. A relationship of vitamin D status with age or sex has not been consistently demonstrated in previous studies of Mesoamerican adults( Reference Diaz, Gonzalez and Ramos-Trautmann 2 , Reference Gonzalez, Ramos-Trautmann and Diaz-Luquis 3 , Reference Flores, Sanchez Romero and Macias 18 – Reference Suarez-Martinez, Perez and Cruz 21 , Reference Caro, Negron and Palacios 26 ). Contrary to our findings, in a study of Puerto Rican adults with BMI ≥25·0 kg/m2, a composite measure representing intake of vitamin D supplements, fish, margarine and milk( Reference Palacios, Gil and Perez 20 ) was positively related to serum 25(OH)D concentrations. The main source of dietary vitamin D in that study was supplement intake, whereas in our study population supplement use was negligible. Discrepancies in the results might be explained by differences in supplement or food intake between populations.
In children, mean serum 25(OH)D concentrations were positively associated with male sex, frequency of dairy intake, and both maternal height and serum 25(OH)D concentration. On the other hand, maternal BMI and the number of home assets were inversely related to 25(OH)D concentrations. Consistent with previous studies, boys had higher concentrations of serum 25(OH)D than girls( Reference Naqvi, Solomons and Campos 1 , Reference Brehm, Celedon and Soto-Quiros 22 , Reference Flores, Flores and Macias 23 , Reference Brehm, Acosta-Perez and Klei 27 ). Explanations for this difference are speculative; they could be due to sex-related variation in intake of foods with high vitamin D content or in time spent playing outdoors. Boys in the current study may have been less likely to wear sunscreen or may have had greater body surface exposure to sunlight compared with girls. We did not find a relationship between age and mean serum 25(OH)D in children, possibly because the age range in our study was limited. Mother’s height was positively related to children’s mean serum 25(OH)D concentrations, perhaps reflecting an association with socio-economic status since adult height can be an indicator of socio-economic status( Reference Subramanian, Ozaltin and Finlay 28 ). There was an inverse association between mothers’ BMI and children’s mean serum 25(OH)D concentrations. An association between mothers’ pre-pregnancy BMI and infants’ vitamin D status has been found previously( Reference Josefson, Feinglass and Rademaker 29 ). The mechanisms by which maternal BMI may influence children’s vitamin D status are largely unknown. Obese women transfer vitamin D to their offspring less effectively at birth( Reference Josefson, Feinglass and Rademaker 29 ), but this would not necessarily explain an association during middle childhood. On the other hand, the mother’s BMI may be positively correlated with the child’s and adiposity is inversely associated with circulating 25(OH)D concentrations( Reference Holick 25 ). Nevertheless, we did not find an association between the children’s BMI and vitamin D status. Mothers’ and children’s serum 25(OH)D concentrations were strongly correlated, consistent with findings from a cross-sectional study of Korean adolescents( Reference Kim, Oh and Namgung 30 ). This possibly reflects the effect of genetic and environmental exposures that are shared between mothers and their children. Serum 25(OH)D concentrations were lower in children from households with >4 home assets compared with those from households with ≤4 assets. This inverse association between a marker of socio-economic status and 25(OH)D is consistent with results from other studies in the region( Reference Flores, Sanchez Romero and Macias 18 , Reference Flores, Flores and Macias 23 , Reference Brehm, Acosta-Perez and Klei 27 ). Finally, in children, frequency of dairy intake was positively related to 25(OH)D concentrations. Fortification of milk with vitamin D is required in Guatemala, Honduras and Mexico( Reference Bowley 31 ) and in other countries some fortified milk is available. Randomized controlled trials of vitamin D-fortified foods demonstrate that intake of fortified foodstuffs results in higher 25(OH)D concentrations( Reference Black, Seamans and Cashman 32 ). Of note, dairy intake in our study was related to 25(OH)D concentrations in children but not in adults. This could be because the absolute consumption of dairy is usually higher and more variable in children than in adults( Reference Quann, Fulgoni and Auestad 33 ).
In our study, about a third of the variability in serum 25(OH)D concentration was explained by family membership. Family aggregation of serum 25(OH)D concentrations suggests that it may be effective to target families, rather than individuals, for screening of VDD and to implement interventions designed to improve vitamin D status. Family aggregation can be due to both shared genes and shared environment. Estimates of heritability of circulating 25(OH)D from twin studies range from 37 to 86 %( Reference Orton, Morris and Herrera 5 , Reference Karohl, Su and Kumari 6 , Reference Arguelles, Langman and Ariza 34 – Reference Snellman, Melhus and Gedeborg 37 ). This could be due to sharing of gene polymorphisms that are associated with 25(OH)D concentration( Reference Dastani, Li and Richards 38 ). We found that correlation and aggregation between mothers’ and fathers’ 25(OH)D concentrations were similar to those between fathers and children, but lower than those between mothers and children. Because mothers and fathers are expected to be genetically unrelated, these comparisons could indicate that shared environmental factors play an important role in the vitamin D serostatus of this population, over that played by shared genetics. For example, sun exposure behaviours may aggregate by family. In Denmark, children’s sun exposure-related behaviours were correlated with their mothers’ but not with their fathers’( Reference Bodekaer Larsen, Petersen and Philipsen 39 ).
The present study has several strengths. Serum 25(OH)D concentrations are widely used in epidemiological studies as an indicator of vitamin D status since they integrate dietary and sunlight sources of vitamin D. The design involving families allowed for the examination of serum 25(OH)D correlates in adults and children, as well as aggregation within families. The use of standardized methods allowed for comparisons across countries where data were not previously available. There are limitations as well. First, the cross-sectional nature of the study limits the ability to make causal inference. Second, the samples were not representative of the underlying populations. Selection bias might have occurred if the probability of inclusion in the study was not independent of the exposures and outcomes examined. Also, the relatively small within-country samples prevented country-specific analysis. The FFQ that we used has not been validated in children. In addition, information on skin pigmentation, total time spent outdoors or sunscreen use was not available. Assuming that some adult participants use sunscreen, the magnitude of the underlying effect of gardening on mean serum 25(OH)D concentration could have been underestimated. Finally, the DiaSorin LIAISON 25-OH Vitamin D Total assay we used to measure 25(OH)D may underestimate its concentration compared with the gold standard of liquid chromatography–tandem mass spectrometry. This may lead to an overestimation of VDD at the population level( Reference Roth, Schmidt-Gayk and Weber 40 ).
In conclusion, the prevalence of VDD was low in a group of adults and children from the Mesoamerican region. Serum 25(OH)D concentrations vary by country in both adults and children and exhibit a moderate degree of aggregation by family. Nationally representative surveys of vitamin D status are needed to determine individual- and country-level predictors of VDD in the region, including the effects of skin pigmentation, the use of sunscreen, clothing habits, environmental toxicants and air pollution.
Acknowledgements
Financial support: This study was funded by the US National Heart, Lung, and Blood Institute (contract number BAA-NHLBI-HV-09-12). The funder did not play any role in the design, conduct or interpretation of the study. Conflicts of interest: None of the authors have any conflicts of interest or financial disclosures. Authorship: M.R.-Z., A.V.R. and E.V. designed the research. S.L.R. and E.V. performed the statistical analyses, wrote the paper and had primary responsibility for the final content. All authors have read and approved the final version of the manuscript. Ethics of human subject participation: The research was conducted in accordance with guidelines laid down by the Declaration of Helsinki. The study was approved by the Health Sciences and Behavioral Sciences Institutional Review Board at the University of Michigan and by the ethics committees of all nine collaborating institutions.
Participants in the Nine Mesoamerican Countries Metabolic Syndrome Study (NiMeCoMeS) Group. Universidad de Ciencias y Artes, Tuxtla Gutiérrez, Chiapas, México: Erika López, Liz Peña, Alejandra Maldonado, Aldeni Vásquez, Aldrin López. University of Belize: Lilly Mahung, Diomar Salazar. Instituto de Nutrición de Centroamérica y Panamá (INCAP), Guatemala: Fernanda Kroker, María Alejandra Córdova, Regina García, Lilian Navas. Universidad Nacional de El Salvador: Josefina Sibrián, Mauricio Flores, Noel Avalos. Universidad de Honduras: Astarté Alegría, Jorge A. Sierra, Héctor Murillo. Universidad Nacional Autónoma de Nicaragua: Ana María Gutiérrez, Carmen María Flores, Mario Romero. Universidad de Costa Rica: Emilce Ulate, Natalia Valverde, Andrea Fiatt, Juan Manuel Valverde. Universidad de Panamá: Flavia Fontes, Raisa Rodríguez, Emérita Pons, Lino Chue, Elka González. Universidad de República Dominicana: Rafael Montero, Francisco Torres, Amarilis Then, Melvi Pérez.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/S1368980017001616
