Multidimensional Bone Density Phenotyping Reveals New Insights Into Genetic Regulation of the Pediatric Skeleton

Jonathan A. Mitchell, Alessandra Chesi, Diana L. Cousminer, Shana E. McCormack, Heidi J. Kalkwarf, Joan M. Lappe, Vicente Gilsanz, Sharon E. Oberfield, John A. Shepherd, Andrea Kelly, Babette S. Zemel, Struan F.A. Grant

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

Osteoporosis is a complex disease with developmental origins. It is therefore important to understand the genetic contribution to pediatric areal bone mineral density (aBMD). Individual skeletal site phenotyping has been primarily used to identify pediatric aBMD loci. However, this approach is limited because there is a degree of aBMD discordance across skeletal sites. We therefore applied a novel multidimensional phenotyping approach to further understand the genetic regulation of pediatric aBMD. Our sample comprised a prospective, longitudinal cohort of 1293 children of European ancestry (52% female; up to seven annual measurements). Principal components analysis was applied to dual-energy X-ray absorptiometry–derived aBMD Z-scores for total hip, femoral neck, spine, and distal radius to generate multidimensional aBMD phenotypes (ie, principal component scores). We tested the association between a genetic score (percentage of bone lowering alleles at 63 loci) and each principal component. We also performed a genomewide association study (GWAS) using the multiethnic baseline data (n = 1885) to identify novel loci associated with these principal components. The first component (PC1) reflected a concordant phenotypic model of the skeleton (eg, higher loading score = higher BMD across all sites). In contrast, PC2 was discordant for distal radius versus spine and hip aBMD, and PC3 was discordant for spine versus distal radius and hip aBMD. The genetic score was associated with PC1 (beta = –0.05, p = 3.9 × 10–10), but was not associated with discordant PC2 or PC3. Our GWAS discovered variation near CPED1 that associated with PC2 (rs67991850, p = 2.5 × 10–11) and near RAB11FIP5 (rs58649746, p = 4.8 × 10–9) that associated with PC3. In conclusion, an established bone fragility genetic summary score was associated with a concordant skeletal phenotype, but not discordant skeletal phenotypes. Novel associations were observed for the discordant multidimensional skeletal phenotypes that provide new biological insights into the developing skeleton.

Original languageEnglish (US)
Pages (from-to)812-821
Number of pages10
JournalJournal of Bone and Mineral Research
Volume33
Issue number5
DOIs
StatePublished - May 1 2018

Fingerprint

Skeleton
Bone Density
Pediatrics
Pelvic Bones
Phenotype
Spine
Bone and Bones
Femur Neck
Principal Component Analysis
Osteoporosis
Hip
Alleles
X-Rays

All Science Journal Classification (ASJC) codes

  • Endocrinology, Diabetes and Metabolism
  • Orthopedics and Sports Medicine

Cite this

Multidimensional Bone Density Phenotyping Reveals New Insights Into Genetic Regulation of the Pediatric Skeleton. / Mitchell, Jonathan A.; Chesi, Alessandra; Cousminer, Diana L.; McCormack, Shana E.; Kalkwarf, Heidi J.; Lappe, Joan M.; Gilsanz, Vicente; Oberfield, Sharon E.; Shepherd, John A.; Kelly, Andrea; Zemel, Babette S.; Grant, Struan F.A.

In: Journal of Bone and Mineral Research, Vol. 33, No. 5, 01.05.2018, p. 812-821.

Research output: Contribution to journalArticle

Mitchell, JA, Chesi, A, Cousminer, DL, McCormack, SE, Kalkwarf, HJ, Lappe, JM, Gilsanz, V, Oberfield, SE, Shepherd, JA, Kelly, A, Zemel, BS & Grant, SFA 2018, 'Multidimensional Bone Density Phenotyping Reveals New Insights Into Genetic Regulation of the Pediatric Skeleton', Journal of Bone and Mineral Research, vol. 33, no. 5, pp. 812-821. https://doi.org/10.1002/jbmr.3362
Mitchell, Jonathan A. ; Chesi, Alessandra ; Cousminer, Diana L. ; McCormack, Shana E. ; Kalkwarf, Heidi J. ; Lappe, Joan M. ; Gilsanz, Vicente ; Oberfield, Sharon E. ; Shepherd, John A. ; Kelly, Andrea ; Zemel, Babette S. ; Grant, Struan F.A. / Multidimensional Bone Density Phenotyping Reveals New Insights Into Genetic Regulation of the Pediatric Skeleton. In: Journal of Bone and Mineral Research. 2018 ; Vol. 33, No. 5. pp. 812-821.
@article{ecf6d927efc848918dfcf58a098304de,
title = "Multidimensional Bone Density Phenotyping Reveals New Insights Into Genetic Regulation of the Pediatric Skeleton",
abstract = "Osteoporosis is a complex disease with developmental origins. It is therefore important to understand the genetic contribution to pediatric areal bone mineral density (aBMD). Individual skeletal site phenotyping has been primarily used to identify pediatric aBMD loci. However, this approach is limited because there is a degree of aBMD discordance across skeletal sites. We therefore applied a novel multidimensional phenotyping approach to further understand the genetic regulation of pediatric aBMD. Our sample comprised a prospective, longitudinal cohort of 1293 children of European ancestry (52{\%} female; up to seven annual measurements). Principal components analysis was applied to dual-energy X-ray absorptiometry–derived aBMD Z-scores for total hip, femoral neck, spine, and distal radius to generate multidimensional aBMD phenotypes (ie, principal component scores). We tested the association between a genetic score (percentage of bone lowering alleles at 63 loci) and each principal component. We also performed a genomewide association study (GWAS) using the multiethnic baseline data (n = 1885) to identify novel loci associated with these principal components. The first component (PC1) reflected a concordant phenotypic model of the skeleton (eg, higher loading score = higher BMD across all sites). In contrast, PC2 was discordant for distal radius versus spine and hip aBMD, and PC3 was discordant for spine versus distal radius and hip aBMD. The genetic score was associated with PC1 (beta = –0.05, p = 3.9 × 10–10), but was not associated with discordant PC2 or PC3. Our GWAS discovered variation near CPED1 that associated with PC2 (rs67991850, p = 2.5 × 10–11) and near RAB11FIP5 (rs58649746, p = 4.8 × 10–9) that associated with PC3. In conclusion, an established bone fragility genetic summary score was associated with a concordant skeletal phenotype, but not discordant skeletal phenotypes. Novel associations were observed for the discordant multidimensional skeletal phenotypes that provide new biological insights into the developing skeleton.",
author = "Mitchell, {Jonathan A.} and Alessandra Chesi and Cousminer, {Diana L.} and McCormack, {Shana E.} and Kalkwarf, {Heidi J.} and Lappe, {Joan M.} and Vicente Gilsanz and Oberfield, {Sharon E.} and Shepherd, {John A.} and Andrea Kelly and Zemel, {Babette S.} and Grant, {Struan F.A.}",
year = "2018",
month = "5",
day = "1",
doi = "10.1002/jbmr.3362",
language = "English (US)",
volume = "33",
pages = "812--821",
journal = "Journal of Bone and Mineral Research",
issn = "0884-0431",
publisher = "Wiley-Blackwell",
number = "5",

}

TY - JOUR

T1 - Multidimensional Bone Density Phenotyping Reveals New Insights Into Genetic Regulation of the Pediatric Skeleton

AU - Mitchell, Jonathan A.

AU - Chesi, Alessandra

AU - Cousminer, Diana L.

AU - McCormack, Shana E.

AU - Kalkwarf, Heidi J.

AU - Lappe, Joan M.

AU - Gilsanz, Vicente

AU - Oberfield, Sharon E.

AU - Shepherd, John A.

AU - Kelly, Andrea

AU - Zemel, Babette S.

AU - Grant, Struan F.A.

PY - 2018/5/1

Y1 - 2018/5/1

N2 - Osteoporosis is a complex disease with developmental origins. It is therefore important to understand the genetic contribution to pediatric areal bone mineral density (aBMD). Individual skeletal site phenotyping has been primarily used to identify pediatric aBMD loci. However, this approach is limited because there is a degree of aBMD discordance across skeletal sites. We therefore applied a novel multidimensional phenotyping approach to further understand the genetic regulation of pediatric aBMD. Our sample comprised a prospective, longitudinal cohort of 1293 children of European ancestry (52% female; up to seven annual measurements). Principal components analysis was applied to dual-energy X-ray absorptiometry–derived aBMD Z-scores for total hip, femoral neck, spine, and distal radius to generate multidimensional aBMD phenotypes (ie, principal component scores). We tested the association between a genetic score (percentage of bone lowering alleles at 63 loci) and each principal component. We also performed a genomewide association study (GWAS) using the multiethnic baseline data (n = 1885) to identify novel loci associated with these principal components. The first component (PC1) reflected a concordant phenotypic model of the skeleton (eg, higher loading score = higher BMD across all sites). In contrast, PC2 was discordant for distal radius versus spine and hip aBMD, and PC3 was discordant for spine versus distal radius and hip aBMD. The genetic score was associated with PC1 (beta = –0.05, p = 3.9 × 10–10), but was not associated with discordant PC2 or PC3. Our GWAS discovered variation near CPED1 that associated with PC2 (rs67991850, p = 2.5 × 10–11) and near RAB11FIP5 (rs58649746, p = 4.8 × 10–9) that associated with PC3. In conclusion, an established bone fragility genetic summary score was associated with a concordant skeletal phenotype, but not discordant skeletal phenotypes. Novel associations were observed for the discordant multidimensional skeletal phenotypes that provide new biological insights into the developing skeleton.

AB - Osteoporosis is a complex disease with developmental origins. It is therefore important to understand the genetic contribution to pediatric areal bone mineral density (aBMD). Individual skeletal site phenotyping has been primarily used to identify pediatric aBMD loci. However, this approach is limited because there is a degree of aBMD discordance across skeletal sites. We therefore applied a novel multidimensional phenotyping approach to further understand the genetic regulation of pediatric aBMD. Our sample comprised a prospective, longitudinal cohort of 1293 children of European ancestry (52% female; up to seven annual measurements). Principal components analysis was applied to dual-energy X-ray absorptiometry–derived aBMD Z-scores for total hip, femoral neck, spine, and distal radius to generate multidimensional aBMD phenotypes (ie, principal component scores). We tested the association between a genetic score (percentage of bone lowering alleles at 63 loci) and each principal component. We also performed a genomewide association study (GWAS) using the multiethnic baseline data (n = 1885) to identify novel loci associated with these principal components. The first component (PC1) reflected a concordant phenotypic model of the skeleton (eg, higher loading score = higher BMD across all sites). In contrast, PC2 was discordant for distal radius versus spine and hip aBMD, and PC3 was discordant for spine versus distal radius and hip aBMD. The genetic score was associated with PC1 (beta = –0.05, p = 3.9 × 10–10), but was not associated with discordant PC2 or PC3. Our GWAS discovered variation near CPED1 that associated with PC2 (rs67991850, p = 2.5 × 10–11) and near RAB11FIP5 (rs58649746, p = 4.8 × 10–9) that associated with PC3. In conclusion, an established bone fragility genetic summary score was associated with a concordant skeletal phenotype, but not discordant skeletal phenotypes. Novel associations were observed for the discordant multidimensional skeletal phenotypes that provide new biological insights into the developing skeleton.

UR - http://www.scopus.com/inward/record.url?scp=85044508655&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85044508655&partnerID=8YFLogxK

U2 - 10.1002/jbmr.3362

DO - 10.1002/jbmr.3362

M3 - Article

VL - 33

SP - 812

EP - 821

JO - Journal of Bone and Mineral Research

JF - Journal of Bone and Mineral Research

SN - 0884-0431

IS - 5

ER -