- Research Article
- Open Access
Genetic variants and physical activity interact to affect bone density in Hispanic children
BMC Pediatrics volume 21, Article number: 79 (2021)
Our aim was to investigate if moderate to vigorous physical activity (MVPA), calcium intake interacts with bone mineral density (BMD)-related single nucleotide polymorphisms (SNPs) to influence BMD in 750 Hispanic children (4-19y) of the cross-sectional Viva La Familia Study.
Physical activity and dietary intake were measured by accelerometers and multiple-pass 24 h dietary recalls, respectively. Total body and lumbar spine BMD were measured by dual energy X-ray absorptiometry. A polygenic risk score (PRS) was computed based on SNPs identified in published literature. Regression analysis was conducted with PRSs, MVPA and calcium intake with total body and lumbar spine BMD.
We found evidence of statistically significant interaction effects between the PRS and MVPA on total body BMD and lumbar spine BMD (p < 0.05). Higher PRS was associated with a lower total body BMD (β = − 0.040 ± 0.009, p = 1.1 × 10− 5) and lumbar spine BMD (β = − 0.042 ± 0.013, p = 0.0016) in low MVPA group, as compared to high MVPA group (β = − 0.015 ± 0.006, p = 0.02; β = 0.008 ± 0.01, p = 0.4, respectively).
The study indicated that calcium intake does not modify the relationship between genetic variants and BMD, while it implied physical activity interacts with genetic variants to affect BMD in Hispanic children. Due to limited sample size of our study, future research on gene by environment interaction on bone health and functional studies to provide biological insights are needed.
Bone health in Hispanic children with high genetic risk for low BMD is benefitted more by MVPA than children with low genetic risk. Our results may be useful to predict disease risk and tailor dietary and physical activity advice delivery to people, especially children.
Osteoporosis and low bone mass are major public health threats with great economic burden in the U.S. Maximizing peak bone mass in early life is an important strategy to prevent osteoporosis in later life. Bonjour et al. (2009) reported that an increase of peak bone mass in childhood by one standard deviation could reduce fracture risk as much as 50% in adulthood . Physical activity and calcium intake are two of the most important modifiable factors for bone health and their benefits to bone health could persist into young adulthood [2, 3]. Ethnic-wise distribution shows that Hispanic adolescents are less physically active, and have lower calcium intake compared to their White counterparts [4, 5]. Despite their importance, the role of physical activity and calcium intake with respect to genetic susceptibility to low bone mass has not been sufficiently studied in Hispanic children.
Studies that account for the full complexity of genetic and environmental factors influencing BMD have been mostly conducted in adults typically focusing on few candidate genes, with very few studies being conducted in children or capturing the effects of multiple genes [6, 7]. Polygenic risk scores (PRSs), which could capture the genetic effects of multiple genes, could be a better approach to assess risk for osteoporosis. In one study of European children, significant interaction was found between physical activity and a PRS computed from selected SNPs identified from published genome-wide association studies (GWAS) . However, none of these studies have been conducted in a Hispanic population nor investigated genome-wide SNPs to identify significant genotypes that would interact with physical activity or calcium intake to affect BMD. Therefore, our aim was to examine if physical activity or calcium intake interact with BMD-associated PRS to affect BMD in an understudied and underserved population of Hispanic children aged 4–19 years old.
Study population: the viva La Familia study (VFS)
The VFS, a family study of 1030 Hispanic children (4-19y) from 319 families enrolled in Houston, TX, was designed to examine the genetic and environmental factors influencing obesity and its comorbidities, with the majority of families of Mexican-American descent [9, 10]. Each family was chosen from an overweight proband (4-19y) with the use of a bivariate procedure, ie, overweight ≥95th percentile for BMI and ≥ 85th percentile for fat mass. In addition, it was required to have ≥3 children between the ages of 4 and 19 y within the families [9, 10]. Related families were also recruited to expand the relative pair types. Interviews were conducted to obtain family pedigree, socio-demographics, lifestyle information, and medical histories at baseline. Puberty stage was self-reported in children. Anthropometrics were obtained from children and their parents . Dual energy X-ray absorptiometry (DXA) with a Delphi-A whole-body scanner (Hologic Inc., Waltham, MA) was used to obtain total body and lumbar spine BMD. Data analysis was conducted in about 750 children for whom complete baseline data (genetic data, dietary intake, physical activity and bone measurements) were available.
Written informed consent or assent was given by all enrolled children’ s parents. The protocol was approved by the Institutional Review Boards for Human Subject Research at Baylor College of Medicine and Affiliated Hospitals and Texas Biomedical Research Institute and the analysis was approved by the University of North Carolina at Chapel Hill. All methods were performed in accordance with the relevant guidelines and regulations.
Actiwatch accelerometers (Mini Mitter, Bend, OR) were used to monitor free-living physical activity continuously for 3 days, with details previously described , yielding activity counts at 1 minute intervals. Using thresholds corresponding to physical activity ratios of 1.5, 3.0, and 6.0 or energy expenditure of 0.01, 0.04, and 0.10 kcal/kg per minute, based on validated study in children . Moderate to vigorous activity (MVPA) is more likely to capture the weight-bearing physical activity than light or sedentary physical activity, therefore we chose MVPA (mins/day) as a key measure for physical activity in this study.
A multiple-pass 24 h dietary recall was recorded on 2 occasions, ~ 2–4 weeks apart, at the Children’s Nutrition Research Center by a registered dietitian using Nutrition Data System Software (database version, 2005, Nutrition Coordinating Center, University of Minnesota, Minneapolis) and food models and household measures/dishware. This system automates interviewing, editing, and coding of dietary intake data. The multiple-pass 24 h dietary recall method uses 3 distinct passes to garner information about a subject’s food intake during the preceding 24 h. The 24-h recalls were obtained without prior notice and children aged < 7 y were assisted by their mothers. Mean nutrient intakes were used in the analysis, computed from the two 24-h dietary recalls .
Genotyping of common variants and sequencing of exome variants
For common variants, genotyping for 1.1 million SNPs was conducted using marker assays included on the Illumina HumanOmni1-Quad v1.0 BeadChips, with details described previously . For exome sequencing, the custom NimbleGen VCRome 2.1 capture reagent was used to capture the entire exome for each DNA sample and details have been previously published .
Weight was measured with a digital balance to the nearest 0.1 kg, and height was measured to the nearest 1 mm with a stadiometer. BMI was calculated as kg/m2 and BMI z-score was calculated based on Centers for Disease Control and Prevention (CDC)  standardized sex and age growth charts on the basis of the L (normality) M (median) S (dispersion) formula:
Tanner stages of maturation were self-reported based on pubic hair and breast and penis development illustrated with drawings . Puberty stages were defined as pre-pubertal (Tanner stage = 1) and pubertal (Tanner stage > 1) .
Index SNPs were independent SNPs that were associated with total body BMD and lumbar spine BMD in a recent paper focusing on life-course genome-wide meta-analysis of bone health [17, 18]. Within a ± 500 kb region of each index SNP, SNPs with the lowest p value in VFS (“VFS_SNP”) were identified. Linkage disequilibrium (LD) between the index SNPs and the corresponding “VFS_SNP” were computed by LD link (https://ldlink.nci.nih.gov) in corresponding ethnic populations in which the index SNPs were originally reported. LD (r2) greater than 0.2 was considered as the threshold that SNPs in VFS might be a better signal of index SNPs. Proxy SNPs (r2 > 0.8 based on LD link, reference population: CEU) within the ±500 kb region of the index SNPs were also identified, when index SNPs found in literature were not available in VFS data. Within each 1 Mb region, the number of independent SNPs was calculated based on LD structure between SNPs in our cohort. The p value threshold for each 1 Mb region was determined using the Bonferroni adjustment (p value threshold).
Polygenic risk score (PRS) computation
SNPs of PRS include: 1) SNPs that are statistically significant or nominally significant per GWAS threshold (p < 1.0 × 10− 7 and p < 1.0 × 10− 6, respectively) and exome sequencing analysis (p < 9.2 × 10− 7 and p < 9.2 × 10− 6, respectively) ; 2) Index SNPs or Proxy SNPs in 1 Mb region that were statistically significantly and nominally significantly associated with bone health (p < p value threshold and p < 0.05, respectively). Nominally significantly associated-SNPs were included as these SNPs may have reached genome-wide significance if we had greater power, and more powerful PRS typically are constructed by incorporating both suggestive and significantly associated SNPs .; 3) “VFS_SNP” that have no strong LD (r2 < 0.2) with index SNPs but were statistically significantly associated with bone outcomes (p < p value threshold); 4) When “VFS_SNP” had strong LD (r2 > 0.2) with index SNPs and were statistically significantly associated with bone outcomes, we conducted conditional analyses including both proxy SNPs and “VFS_SNP” to determine the better signal in place of the Index SNPs, and the better one was included in PRS. As genetic variants related to BMD are site specific , separate PRS for total body BMD and lumbar spine BMD were computed, with details shown in Table 1. In addition, we conducted sensitivity analysis with PRS that included only statistically significant SNPs (secondary PRS), with details shown in Additional file 1. Minor allele frequencies of selected SNPs were above 0.01. Each SNP was recoded as 0, 1, or 2 based on the number of risk alleles (associated with lower BMD). PRS was constructed by summing up the risk alleles at each locus assuming equal effects at each locus: PRS = SNP1 + SNP2 + … + SNPn.
Total body BMD and lumbar spine BMD were inverse normalized to attain normal distribution in the analysis [21, 22]. PRS was used as continuous variables as we postulated a linear relationship, but we also used binary variables categorized based on medians (> = median vs below) to test sensitivity of results. Environmental variables: dietary intake (calcium, calcium to phosphorous ratio) and MVPA (mins/d) were included alternatively as continuous variables or binary variables categorized based on median (> = median vs below) . Covariates included age, sex, BMI Z-score, and puberty stage. As BMD tends to be lower for smaller bones compared to larger bones , we considered adding height into the current model for sensitivity analysis. However, as height was strongly correlated with age in our cohort (r = 0.91, p value < 0.001), we did not add height as a covariate, in addition to age, in the model to avoid issues of collinearity.
First, main effects were examined by conducting association analysis between environmental variables and BMD with linear regression. For the gene-environment interaction analysis, the modified generalized estimating equation was used and the regression model included the PRSs, the environment variables (calcium intake, calcium/phosphorus, MVPA, alternatively), the interactions term between the PRS and environment variables, and other covariates (age, sex BMI Z-score and puberty stage). For interaction analysis between PRS and environmental factors, the statistical significance level was considered as p < 0.05. Interactions between individual variants and environmental factors were considered as secondary analyses. Statistical analyses were performed using the software SUGEN (http://dlin.web.unc.edu/software/SUGEN/) .
The characteristics of VFS participants are shown in Table 2 and Additional file 1. The children in VFS have a mean age and BMI of 11.0 ± 4.1 y and 25.1 ± 7.6 kg/m2, respectively. The mean calcium intake was 881.5 ± 377.3 mg with a median of 841.5 mg, which is lower than recommended dietary allowances for calcium . The mean dietary calcium/phosphorous ratio was 0.75 ± 0.19 with a median of 0.75, which is lower than the dietary optimal ratio of 1.3 . The mean MVPA was 82.3 ± 51.8 mins/d with a median of 74 mins/d, which is higher than the recommended physical activity level for youth . The lower TBBMD and LSBM observed in high MVPA group could be potentially due to younger age and higher number of pre-pubertal children in high MVPA group. The characteristics of VFS children by PRS groups are shown in Additional file 1.
Moderate to vigorous physical activity (MVPA) – genetic interaction on BMD
No statistically significant associations were found between MVPA and BMD or genetic variants and MVPA. However, we found statistically significant effects of interactions between binary MVPA (mins/d) and total body BMD PRS on total body BMD (p = 0.005) (Table 3). Children with low MVPA (mins/d) on average have 0.040 g/cm2 total body BMD less in the low PRS group as compared to the high PRS group, whereas children with high MVPA (mins/d) on average have 0.015 g/cm2 total body BMD less in the low PRS group as compared to the high PRS group. We also found statistically significant interactions between binary MVPA (mins/d) and total body BMD PRS in relation to lumbar spine BMD (p = 0.001) (Table 3). Children with low MVPA (mins/d) on average have 0.042 g/cm2 lumbar spine BMD less in the low PRS group as compared to the high PRS group, whereas children with high MVPA (mins/d) on average have 0.008 g/cm2 lumbar spine BMD less in the low PRS group as compared to the high PRS group. In further analysis with individual SNPs included in the PRS, a statistically significant interaction was found between MVPA (mins/d) and rs10953178 at 26S Proteasome Complex Subunit (SEM1) (p = 0.046) and rs914153 (p = 0.028) to affect total body BMD, and rs452369 at Solute Carrier Family 8 Member A1 (SLC8A1) (p = 0.038), rs10953178 (SEM1) (p = 0.0074) and rs914153 at chromosome 21q21 (p = 0.0058) to affect lumbar spine BMD (Table 3). In the sensitivity analysis using secondary PRS, similar results were found (Additional file 1).
Calcium intake– genetic interaction on BMD
No statistically significant associations were found between calcium intake and BMD, but dietary calcium/phosphorous ratio was positively associated with total body BMD (p = 0.02). Also, no genome-wide or exome-wide significant genetic variants were identified for calcium intake and dietary calcium/phosphorous ratio. We also found no statistically significant interaction effects between PRS and calcium intake or dietary calcium/phosphorous ratio in association with total body BMD or lumbar spine BMD. However, in the secondary analysis with individual SNPs, we found significant interaction effects between calcium/phosphorous ratio and rs370055571 of Mitogen-Activated Protein Kinase Kinase Kinase Kinase 3 (MAP 4 K3) in relation to total body BMD (p = 0.02) and rs452369 of SLC8A1 in relation to lumbar spine BMD (p = 0.036). In sensitivity analysis with secondary PRS, we found statistically significant interaction effects of lumbar spine BMD PRS with binary dietary calcium/phosphorous ratio on lumbar spine BMD (p = 0.037), indicating that high dietary calcium/phosphorous ratio was more beneficial in children with higher genetic susceptibility to low lumbar spine BMD as compared with children with lower genetic risk (Additional file 1).
All individual SNPs interaction results were in Additional file 1.
Our study addresses the research gap to investigate the effects of interaction between PRSs and calcium intake and physical activity on BMD in Hispanic children. We found higher MVPA to be more beneficial for children with higher genetic predisposition to low bone density in comparison to children less prone to low bone density. No statistically significant interaction effects on BMD were found between genetic risk and calcium intake or dietary calcium/phosphorous ratio.
Most studies that investigated interaction between physical activity and genetic variants used a candidate gene approach in primarily adult populations . Studies have found that estrogen receptor-alpha (ESRα) polymorphisms may modulate the effect of exercise on BMD at weight bearing sites in children . The vitamin D receptor (VDR) genotypes were also attributed to bone metabolic response variance following resistance training in young adult males . One recent study in European children constructed a PRS based on GWAS-identified BMD loci in adults and found no evidence of effect of physical activity-PRS interactions on BMD, but did find beneficial effects on BMD of high impact physical activity on children overall with high genetic risk of low BMD . We, on the contrary, found significant evidence of interaction between MVPA and PRSs in relation to total body BMD and lumbar spine BMD. Moreover, our results suggest that higher MVPA may benefit the BMD of children genetically susceptible to low BMD, although it may not provide meaningful benefit to children less genetically susceptible to low BMD. It is possible that these genes influence the bones’ responses to physical activity or physical activity changes the way these genes are expressed. Larger sample sizes may be needed to detect the interaction effects between continuous MVPA variables and PRS on BMD, compared to binary MVPA. In individual SNPs analysis, significant effect of interaction between MVPA and SLC8A1, SEM1 loci and rs914153 were found to associate with BMD. Both SLC8A1 and SEM1 have been reported to associate with BMD or fractures in GWAS studies [17, 31]. Exercise training was found to reduce SLC8A1 expression in animals susceptible to heart diseases  and SLC8A1 was associated with hand grip strength and bone density [17, 33]. Exercise training was also found in mice to affect ubiquitin-proteasome systems , which has a role in skeletal muscle remodeling during exercise . As skeletal muscle and muscle strength could influence bone strength , exercise might influence SLC8A1 expression and interact with SEM1 to affect ubiquitin-proteasome systems to further affect bone density. The SNP rs914153 is a novel locus that has not been reported to be associated with bone health or physical activity, so the mechanism is largely unknown. Given the scarcity of transcriptomic and epigenomic data from bone tissue , further research is needed to determine the functional impact of these genetic variants.
We did not find any statistically significant interaction between PRS and calcium intake or dietary calcium/phosphorous ratio to affect total body BMD and lumbar spine BMD. Calcium intake and dietary calcium/phosphorous ratio are mainly below the recommendation levels in our study and the genetic variants in PRS may exert more influence on BMD with higher calcium intake. Also, calcium intake was found in randomized controlled trials in pre-pubertal children to mainly increase BMD at appendicular skeletal sites but not lumbar spine potentially due to different calcification proportion in different sites [38, 39], which may partly explain the null interaction results on lumbar spine BMD. As PRSs only included several SNPs and our sample size is small, more genetic variants and larger sample size may be needed to detect the interaction effects. In individual variants analysis, we observed statistical interactions between dietary calcium/phosphorous ratio with MAP 4 K3 and SLC8A1 loci in relation to total body BMD and lumbar spine BMD respectively. As far as we know, no previous literature has reported the association between MAP 4 K3 and bone health. However, the interaction is biologically plausible. MAP 4 K3 encodes a kinase that is involved in the MAPK signaling pathway, which is critical to skeletal developments and maintenance and regulated by calcium . SLC8A1, which encodes a sodium/calcium exchanger, was found to be associated with heel BMD and total body BMD . Not all individual SNPs included in the PRS were found to interact with MVPA, calcium intake or calcium/phosphorous ratio probably due to involvements in different biological pathways or small effect sizes.
Our study has a few limitations. First, DXA-derived BMD is areal BMD (g/cm2) rather than volumetric BMD; although it is commonly used and considered as the gold standard to assess pediatric bone health, areal BMD tends to be lower for smaller bones compared to larger bones, even when the volumetric BMD is the same . Also, even though total body excluding head is recommended, we used total body BMD as outcome due to SNPs associated with total body BMD may be more appropriate to assess variation from childhood to old age . Second, there are some challenges to measure dietary intake in children accurately including parents’ ability to report young children’s food intake, children’s reliability to recall food intake, and adolescents’ concerns about body image, which results in greater within-between subject variances in nutrient intake in children compared to adults . Although 2 days measurement of calcium intake were estimated to have good accuracy in adults , more days may be required to capture the usual intake of calcium. Moreover, we did not capture the specific type of physical activity in accelerometers, especially, weight-bearing physical activity, which is well-known to be beneficial to BMD. But MVPA is more likely to capture the weight-bearing physical activity compared to light or sedentary physical activity. Our accelerometer may also deviate from the current recommendations . Also, the lumbar spine was not analyzed from a spine densitometry. Furthermore, although we account for main factors including age, sex, puberty stage, and BMI Z-scores, our results may not be comparable to all children studies on BMD due to the inconsistent covariates adjustment across studies and there could still be residual confounding . Additionally, the sample size is small and may suffer from false positive in the results. Future studies are needed to replicate and validate our findings. Also, limited number of SNPs are included and the conclusions could be different if different SNPs are included in the PRS. Finally, as we used cross-sectional data in the analysis, the results may not reflect causal processes and future research in longitudinal studies is needed.
Our study also has many strengths. To our knowledge, this is the first study that investigated interactions between dietary calcium and physical activity with PRS affecting BMD in Hispanic children. First, we constructed the PRS based on literature and our association results from GWAS and exome sequencing data, which could account for joint effects of multiple genetic variants and improve power for detection of interactions . Second, we created site-specific PRSs for interaction analysis due to genetic specificity to skeletal site , which may provide insights on how environments associate with bone density in people with skeletal site-specific genetic risk. Third, we had physical activity data measured by accelerometers, which is an objective, nonreactive approach to estimate physical activity compared to self-reported questionnaires .
In summary, MVPA is potentially more beneficial to bone health in Hispanic children who are at greater genetic risk for lower bone density. Future research on gene by environment interaction on bone health and functional studies to provide biological insights are needed.
Availability of data and materials
Data described in the manuscript, code book and analytic code will be made available upon request pending approval with corresponding author and all sequencing data described in the manuscript is publically available without restriction at https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000616.v2.p2.
Body mineral density
Moderate to vigorous physical activity
Single nucleotide polymorphisms
Dual energy X-ray absorptiometry
Polygenic Risk Score
- MAP 4K3 :
Mitogen-Activated Protein Kinase Kinase Kinase Kinase 3
- SLC8A1 :
Solute Carrier Family 8 Member A1
- SEM1 :
26S Proteasome Complex Subunit
- VDR :
Vitamin D receptor
- IL6 :
- ESR1 :
Estrogen receptor α
The Viva La Familia Study
Bonjour J-P, Chevalley T, Ferrari S, Rizzoli R. The importance and relevance of peak bone mass in the prevalence of osteoporosis. Salud Publica Mex. 2009;51:s5–17.
Gunter KB, Almstedt HC, Janz KF. Physical activity in childhood may be the key to optimizing lifespan skeletal health. Exerc Sport Sci Rev. 2012;40:13.
Ilich JZ, Kerstetter JE. Nutrition in bone health revisited: a story beyond calcium. J Am Coll Nutr. 2000;19:715–37.
Duke J, Huhman M, Heitzler C. Physical activity levels among children aged 9-13 years-United States, 2002. JAMA. 2003;290:1308.
Jackson KA, Savaiano DA. Lactose Maldigestion, calcium intake and osteoporosis in African-, Asian-, and Hispanic-Americans. J Am Coll Nutr. 2001;20:198S–207S.
Bonjour J-P, Chevalley T, Rizzoli R, Ferrari S. Gene-environment interactions in the skeletal response to nutrition and exercise during growth. Med Sport Sci. 2007;51:64–80.
Ferrari S, Karasik D. Gene--diet interactions on bone. In: Holick MF, Nieves JW, editors. Nutrition and bone health. New York: Springer New York; 2015. p. 21–36.
Mitchell JA, Chesi A, Elci O, McCormack SE, Roy SM, Kalkwarf HJ, et al. Physical activity benefits the skeleton of children genetically predisposed to lower bone density in adulthood. J Bone Miner Res. 2016;31:1504–12.
Butte NF, Cai GW, Cole SA, Comuzzie AG. Viva la Familia study: genetic and environmental contributions to childhood obesity and its comorbidities in the Hispanic population. Am J Clin Nutr. 2006;84:646–54.
Voruganti VS, Laston S, Haack K, Mehta NR, Cole SA, Butte NF, et al. Serum uric acid concentrations and SLC2A9 genetic variation in Hispanic children: the viva La Familia study. Am J Clin Nutr April. 2015;101:725–32.
Cai G, Cole SAA, Butte N, Bacino C, Diego V, Tan K, et al. A quantitative trait locus on chromosome 18q for physical activity and dietary intake in Hispanic children. Obesity (Silver Spring). 2006;14:1596–604.
Butte NF, Cai G, Cole SA, Wilson TA, Fisher JO, Zakeri IF, et al. Metabolic and behavioral predictors of weight gain in Hispanic children: the viva la Familia study. Am J Clin Nutr. 2007;85:1478–85.
Wilson TA, Adolph AL, Butte NF. Nutrient adequacy and diet quality in non-overweight and overweight Hispanic children of low socioeconomic status: the viva la Familia study. J Am Diet Assoc. 2009;109:1012–21.
Hou R, Cole SA, Graff M, Haack K, Laston S, Comuzzie AG, et al. Genetic variants affecting bone mineral density and bone mineral content at multiple skeletal sites in Hispanic children. Bone. 2020;132:115175.
Kuczmarski RJ, Ogden CL, Grummer-Strawn LM, Flegal KM, Guo SS, Wei R, et al. CDC growth charts. Adv Data. 2000;2000(314):1–27.
Antiretroviral therapy of HIV infection in infants and children: towards universal access: recommendations for a public health approach-2010 revision. World Health Organization. 2010. pp 1-94. ISBN-978-82-4-15-9980-1.
Medina-Gomez C, Kemp JP, Trajanoska K, Luan J, Chesi A, Ahluwalia TS, et al. Life-course genome-wide association study meta-analysis of Total body BMD and assessment of age-specific effects. Am J Hum Genet. 2018;102:88–102.
Mitchell JA, Cousminer DL, Zemel BS, Grant SFA, Chesi A. Genetics of pediatric bone strength. Bonekey Rep. 2016;5:1–6.
Evans DM, Visscher PM, Wray NR. Harnessing the information contained within genome-wide association studies to improve individual prediction of complex disease risk. Hum Mol Genet. 2009;18:3525–31.
Kemp JP, Medina-Gomez C, Estrada K, St Pourcain B, Heppe DH, Warrington NM, et al. Phenotypic dissection of bone mineral density reveals skeletal site specificity and facilitates the identification of novel loci in the genetic regulation of bone mass attainment. PLoS Genet. 2014;10(6):e1004423.
Kocks J, Ward K, Mughal Z, Moncayo R, Adams J, Högler W. Z-score comparability of bone mineral density reference databases for children. J Clin Endocrinol Metab. 2010;95:4652–9.
Leonard MB, Propert KJ, Zemel BS, Stallings VA, Feldman HI. Discrepancies in pediatric bone mineral density reference data: potential for misdiagnosis of osteopenia. J Pediatr. 1999;135(2 Pt 1):182–8.
2018 Physical Activity Guidelines Advisory Committee. 2018 Physical Activity Guidelines Advisory Committee Scientific Report. 2018.
Binkovitz LA, Henwood MJ. Pediatric DXA: Technique and interpretation. Pediatr Radiol. 2007;37:21–31.
Lin DY, Tao R, Kalsbeek WD, Zeng D, Gonzalez F, Fernández-Rhodes L, et al. Genetic association analysis under complex survey sampling: the Hispanic community health study/study of Latinos. Am J Hum Genet. 2014;95:675–88.
Ross CA, Taylor CL, Yaktine AL, del Vale HB. Committee to review dietary reference intakes for vitamin D and calcium food and nutrition board; Institute of Medicine. Diet Ref Intakes Calcium Vitam D. The National Academies Press; 2010.
Burckhardt P, Dawson-Hughes B, Weaver C. Nutritional influences on bone health; 2010.
Piercy KL, Troiano RP, Ballard RM, Carlson SA, Fulton JE, Galuska DA, et al. The physical activity guidelines for Americans. JAMA. 2018;320:2020–8.
Suuriniemi M, Mahonen A, Kovanen V, Alén M, Lyytikäinen A, Wang Q, et al. Association between exercise and pubertal BMD is modulated by estrogen receptor α genotype. J Bone Miner Res. 2004;19:1758–65.
Tajima O, Ashizawa N, Ishii T, Amagai H, Mashimo T, Liu LJ, et al. Interaction of the effects between vitamin D receptor polymorphism and exercise training on bone metabolism. J Appl Physiol. 2000;88:1271–6.
Estrada K, Styrkarsdottir U, Evangelou E, Hsu Y-H, Duncan EL, Ntzani EE, et al. Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture. Nat Genet. 2012;44:491–501.
Kukielka M, Holycross B, Billman GE. Endurance exercise training reduces cardiac sodium/calcium exchanger expression in animals susceptible to ventricular fibrillation. Front Physiol. 2011;2:3.
Willems SM, Wright DJ, Day FR, Trajanoska K, Joshi PK, Morris JA, et al. Large-scale GWAS identifies multiple loci for hand grip strength providing biological insights into muscular fitness. Nat Commun. 2017;8:16015.
Cunha TF, Moreira JBN, Paixao NA, Campos JC, Monteiro AWA, Bacurau AVN, et al. Aerobic exercise training upregulates skeletal muscle calpain and ubiquitin-proteasome systems in healthy mice. J Appl Physiol. 2012;112:1839–46.
Taillandier D, Combaret L, Pouch M-N, Samuels SE, Béchet D, Attaix D. The role of ubiquitin–proteasome-dependent proteolysis in the remodelling of skeletal muscle. Proc Nutr Soc. 2004;63:357–61.
Schönau E, Werhahn E, Schiedermaier U, Mokow E, Schiessl H, Scheidhauer K, et al. Influence of muscle strength on bone strength during childhood and adolescence. Horm Res Paediatr. 1996;45(Suppl. 1):63–6.
Sabik OL, Farber CR. Using GWAS to identify novel therapeutic targets for osteoporosis. Transl Res. 2017;181:15–26.
Chevalley T, Bonjour JP, Ferrari S, Hans D, Rizzoli R. Skeletal site selectivity in the effects of calcium supplementation on areal bone mineral density gain: a randomized, double-blind, placebo-controlled trial in prepubertal boys. J Clin Endocrinol Metab. 2005;90:3342–9.
Bonjour J-P, Carrie A-L, Ferrari S, Clavien H, Slosman D, Theintz G, et al. Calcium-enriched foods and bone mass growth in prepubertal girls: a randomized, double-blind, placebo-controlled trial. J Clin Invest. 1997;99:1287–94.
Thouverey C, Caverzasio J. Focus on the p38 MAPK signaling pathway in bone development and maintenance. Bonekey Rep. 2015;4:1–8.
Livingstone MBE, Robson PJ. Measurement of dietary intake in children. Proc Nutr Soc. 2000;59:279–93.
Palaniappan U, Cue RI, Payette H, Gray-Donald K. Implications of day-to-day variability on measurements of usual food and nutrient intakes. J Nutr. 2003;133:232–5.
Migueles JH, Cadenas-Sanchez C, Ekelund U, Nyström CD, Mora-Gonzalez J, Löf M, et al. Accelerometer data collection and processing criteria to assess physical activity and other outcomes: a systematic review and practical considerations. Sport Med. 2017;47:1821–45.
Harvey NC, Cole ZA, Crozier SR, Kim M, Ntani G, Goodfellow L, et al. Physical activity, calcium intake and childhood bone mineral: a population-based cross-sectional study. Osteoporos Int. 2012;23:121–30.
Liu S, Song Y. Building genetic scores to predict risk of complex diseases in humans: is it possible? Diabetes. 2010;59:2729–31.
Trost SG. State of the art reviews: measurement of physical activity in children and adolescents. Am J Lifestyle Med. 2007;1:299–314.
The Viva La Familia Study was supported by National Institute of Health DK080457 and United States Department of Agriculture/Agricultural Research Service (Corporative agreement 625–51000-053). The development of SOLAR was supported by National Institute of Health grant MH59490.
Ethics approval and consent to participate
Written informed consent or assent was given by all enrolled children’s parents. The protocol was approved by the Institutional Review Boards for Human Subject Research at Baylor College of Medicine and Affiliated Hospitals and Texas Biomedical Research Institute and the analysis was approved by the University of North Carolina at Chapel Hill.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplemental Table 1. Genetic Risk Score Compositions with only statistically significant BMD-associated SNPs a, b. Supplemental Table 2. Baseline characteristics of children in the Viva La Familia Study by PRS. Supplemental Table 3. Interaction between PRS (with only statistically significant BMD-associated SNPs) and moderate to vigorous physical activity on TBBMD and LSBMD a, b, c, d. Supplemental Table 4. Interaction between genetic markers and calcium intake, dietary calcium/phosphorous ratio on TBBMD and LSBMD a, b, c, d.
About this article
Cite this article
Hou, R., Cole, S.A., Graff, M. et al. Genetic variants and physical activity interact to affect bone density in Hispanic children. BMC Pediatr 21, 79 (2021). https://doi.org/10.1186/s12887-021-02537-y