Genes with specificity for expression in the round cell layer of the growth plate are enriched in genomewide association study (GWAS) of human height
Human adult height reflects the outcome of childhood skeletal growth. Growth plate (epiphyseal) chondrocytes are key determinants of height. As epiphyseal chondrocytes mature and proliferate, they pass through three developmental stages, which are organized into three distinct layers in the growth plate: (i) resting (round), (ii) proliferative (flat), and (iii) hypertrophic. Recent genomewide association studies (GWASs) of human height identified numerous associated loci, which are enriched for genes expressed in growth plate chondrocytes. However, it remains unclear which specific genes expressed in which layers of the growth plate regulate skeletal growth and human height. To connect the genetics of height and growth plate biology, we analyzed GWAS data through the lens of gene expression in the three dissected layers of murine newborn tibial growth plate. For each gene, we derived a specificity score for each growth plate layer and regressed these scores against gene-level p values from recent height GWAS data. We found that specificity for expression in the round cell layer, which contains chondrocytes early in maturation, is significantly associated with height GWAS p values (p = 8.5?×?10-9); this association remains after conditioning on specificity for the other cell layers. The association also remains after conditioning on membership in an “Online Mendelian Inheritance in Man (OMIM) gene set” (genes known to cause monogenic skeletal growth disorders, p?<?9.7?×?10-6). We replicated the association in RNA-sequencing (RNA-seq) data from maturing chondrocytes sampled at early and late time points during differentiation in vitro: we found that expression early in differentiation is significantly associated with p values from height GWASs (p = 6.1?×?10-10) and that this association remains after conditioning on expression at 10?days in culture and on the OMIM gene set (p?<?0.006). These findings newly implicate genes highlighted by GWASs of height and specifically expressed in the round cell layer as being potentially important regulators of skeletal biology. © 2021 American Society for Bone and Mineral Research (ASBMR). Overall design: WT GPLCs were grown for 1, 3, 5, or 10 days in culture. RNAseq was performed to compare these timepoints.