Emerging Insights into the Endocrine Regulation of Bone Homeostasis by Gut Microbiome.

Gut microbiota plays an important role in the regulation of bone homeostasis and bone health. Recent studies showed that these effects could be mediated through microbial metabolites released by the microbiota like short-chain fatty acids, metabolism of endogenous molecules such as bile acids, or a complex interplay between microbiota, the endocrine system, and the immune system. Importantly, some studies showed a reciprocal relationship between the endocrine system and gut microbiota. For instance, postmenopausal estrogen deficiency could lead to dysbiosis of the gut microbiota, which could in turn affect various immune response and bone remodeling. In addition, evidence showed that shift in the indigenous gut microbiota caused by antibiotics treatment may also impact normal skeletal growth and maturation. In this mini-review, we describe recent findings on the role of microbiome in bone homeostasis, with a particular focus on molecular mechanisms and their interactions with the endocrine and immune system. We will also discuss the recent findings on estrogen deficiency and microbiota dysbiosis, and the clinical implications for the development of new therapeutic strategies for osteoporosis and other bone disorders.

Nutrition, Other Environmental Influences, and Genetics in the Determination of Human Stature.

Linear growth during three distinct stages of life determines attained stature in adulthood: namely, in utero, early postnatal life, and puberty and the adolescent period. Individual host factors, genetics, and the environment, including nutrition, influence attained human stature. Each period of physical growth has its specific biological and environmental considerations. Recent epidemiologic investigations reveal a strong influence of prenatal factors on linear size at birth that in turn influence the postnatal growth trajectory. Although average population height changes have been documented in high-income regions, stature as a complex human trait is not well understood or easily modified. This review summarizes the biology of linear growth and its major drivers including nutrition from a life-course perspective, the genetics of programmed growth patterns or height, and gene-environment interactions that determine human stature in toto over the life span. Implications for public health interventions and knowledge gaps are discussed.

Growth disorders caused by variants in epigenetic regulators: progress and prospects.

Epigenetic modifications play an important role in regulation of transcription and gene expression. The molecular machinery governing epigenetic modifications, also known as epigenetic regulators, include non-coding RNA, chromatin remodelers, and enzymes or proteins responsible for binding, reading, writing and erasing DNA and histone modifications. Recent advancement in human genetics and high throughput sequencing technology have allowed the identification of causative variants, many of which are epigenetic regulators, for a wide variety of childhood growth disorders that include skeletal dysplasias, idiopathic short stature, and generalized overgrowth syndromes. In this review, we highlight the connection between epigenetic modifications, genetic variants in epigenetic regulators and childhood growth disorders being established over the past decade, discuss their insights into skeletal biology, and the potential of epidrugs as a new type of therapeutic intervention.

Epigenetic Causes of Overgrowth Syndromes.

Human overgrowth disorders are characterized by excessive prenatal and/or postnatal growth of various tissues. These disorders often present with tall stature, macrocephaly, and/or abdominal organomegaly and are sometimes associated with additional phenotypic abnormalities such as intellectual disability and increased cancer risk. As the genetic etiology of these disorders have been elucidated, a surprising pattern has emerged. Multiple monogenic overgrowth syndromes result from variants in epigenetic regulators: variants in histone methyltransferases NSD1 and EZH2 cause Sotos syndrome and Weaver syndrome, respectively, variants in DNA methyltransferase DNMT3A cause Tatton-Brown-Rahman syndrome, and variants in chromatin remodeler CHD8 cause an autism spectrum disorder with overgrowth. In addition, very recently, a variant in histone reader protein SPIN4 was identified in a new X-linked overgrowth disorder. In this review, we discuss the genetics of these overgrowth disorders and explore possible common underlying mechanisms by which epigenetic pathways regulate human body size.

Gut microbiota in regulation of childhood bone growth.

Childhood stunting and wasting, or decreased linear and ponderal growth associated with undernutrition, continue to be a major global public health challenge. Although many of the current therapeutic and dietary interventions have significantly reduced childhood mortality caused by undernutrition, there remain great inefficacies in improving childhood stunting. Longitudinal bone growth in children is governed by different genetic, nutritional and other environmental factors acting systemically on the endocrine system and locally at the growth plate. Recent studies have shown that this intricate interplay between nutritional and hormonal regulation of the growth plate could involve the gut microbiota, highlighting the importance of a holistic approach in tackling childhood undernutrition. In this review, I focus on the mechanistic insights provided by these recent advances in gut microbiota research and discuss ongoing development of microbiota-based therapeutics in humans, which could be the missing link in solving undernutrition and childhood stunting.

Identification of novel genes including NAV2 associated with isolated tall stature.

Very tall people attract much attention and represent a clinically and genetically heterogenous group of individuals. Identifying the genetic etiology can provide important insights into the molecular mechanisms regulating linear growth. We studied a three-generation pedigree with five isolated (non-syndromic) tall members and one individual with normal stature by whole exome sequencing; the tallest man had a height of 211 cm. Six heterozygous gene variants predicted as damaging were shared among the four genetically related tall individuals and not present in a family member with normal height. To gain insight into the putative role of these candidate genes in bone growth, we assessed the transcriptome of murine growth plate by microarray and RNA Seq. Two (Ift140, Nav2) of the six genes were well-expressed in the growth plate. Nav2 (p-value 1.91E-62) as well as Ift140 (p-value of 2.98E-06) showed significant downregulation of gene expression between the proliferative and hypertrophic zone, suggesting that these genes may be involved in the regulation of chondrocyte proliferation and/or hypertrophic differentiation. IFT140, NAV2 and SCAF11 have also significantly associated with height in GWAS studies. Pathway and network analysis indicated functional connections between IFT140, NAV2 and SCAF11 and previously associated (tall) stature genes. Knockout of the all-trans retinoic acid responsive gene, neuron navigator 2 NAV2, in Xenopus supports its functional role as a growth promotor. Collectively, our data expand the spectrum of genes with a putative role in tall stature phenotypes and, among other genes, highlight NAV2 as an interesting gene to this phenotype.

Loss-of-function variant in SPIN4 causes an X-linked overgrowth syndrome.

Overgrowth syndromes can be caused by pathogenic genetic variants in epigenetic writers, such as DNA and histone methyltransferases. However, no overgrowth disorder has previously been ascribed to variants in a gene that acts primarily as an epigenetic reader. Here, we studied a male individual with generalized overgrowth of prenatal onset. Exome sequencing identified a hemizygous frameshift variant in Spindlin 4 (SPIN4), with X-linked inheritance. We found evidence that SPIN4 binds specific histone modifications, promotes canonical WNT signaling, and inhibits cell proliferation in vitro and that the identified frameshift variant had lost all of these functions. Ablation of Spin4 in mice recapitulated the human phenotype with generalized overgrowth, including increased longitudinal bone growth. Growth plate analysis revealed increased cell proliferation in the proliferative zone and an increased number of progenitor chondrocytes in the resting zone. We also found evidence of decreased canonical Wnt signaling in growth plate chondrocytes, providing a potential explanation for the increased number of resting zone chondrocytes. Taken together, our findings provide strong evidence that SPIN4 is an epigenetic reader that negatively regulates mammalian body growth and that loss of SPIN4 causes an overgrowth syndrome in humans, expanding our knowledge of the epigenetic regulation of human growth.

A neomorphic variant in SP7 alters sequence specificity and causes a high-turnover bone disorder.

SP7/Osterix is a transcription factor critical for osteoblast maturation and bone formation. Homozygous loss-of-function mutations in SP7 cause osteogenesis imperfecta type XII, but neomorphic (gain-of-new-function) mutations of SP7 have not been reported in humans. Here we describe a de novo dominant neomorphic missense variant (c.926 C > G:p.S309W) in SP7 in a patient with craniosynostosis, cranial hyperostosis, and long bone fragility. Histomorphometry shows increased osteoblasts but decreased bone mineralization. Mice with the corresponding variant also show a complex skeletal phenotype distinct from that of Sp7-null mice. The mutation alters the binding specificity of SP7 from AT-rich motifs to a GC-consensus sequence (typical of other SP family members) and produces an aberrant gene expression profile, including increased expression of Col1a1 and endogenous Sp7, but decreased expression of genes involved in matrix mineralization. Our study identifies a pathogenic mechanism in which a mutation in a transcription factor shifts DNA binding specificity and provides important in vivo evidence that the affinity of SP7 for AT-rich motifs, unique among SP proteins, is critical for normal osteoblast differentiation.

Evidence That Non-Syndromic Familial Tall Stature Has an Oligogenic Origin Including Ciliary Genes.

Human growth is a complex trait. A considerable number of gene defects have been shown to cause short stature, but there are only few examples of genetic causes of non-syndromic tall stature. Besides rare variants with large effects and common risk alleles with small effect size, oligogenic effects may contribute to this phenotype. Exome sequencing was carried out in a tall male (height 3.5 SDS) and his parents. Filtered damaging variants with high CADD scores were validated by Sanger sequencing in the trio and three other affected and one unaffected family members. Network analysis was carried out to assess links between the candidate genes, and the transcriptome of murine growth plate was analyzed by microarray as well as RNA Seq. Heterozygous gene variants in CEP104, CROCC, NEK1, TOM1L2, and TSTD2 predicted as damaging were found to be shared between the four tall family members. Three of the five genes (CEP104, CROCC, and NEK1) belong to the ciliary gene family. All genes are expressed in mouse growth plate. Pathway and network analyses indicated close functional connections. Together, these data expand the spectrum of genes with a role in linear growth and tall stature phenotypes.

Laser Capture Microdissection of Mouse Growth Plate Cartilage.

The ability to identify, isolate, and study pure populations of cells is critical for understanding normal physiology in organs and tissues, which involves spatial regulation of signaling pathways and interactions between cells with different functions, expression profiles, and lineages. Here, we focus on assessing the growth plate cartilage, composed of multiple functionally and histologically distinct zones, to investigate temporally and spatially dependent gene expression differences. In this chapter, we describe the method of laser capture microdissection to isolate chondrocytes from different zones of differentiation in the mouse growth plate cartilage for RNA isolation, and subsequent downstream applications, such as RNA-sequencing and quantitative real-time PCR. We also provide an assessment of different factors contributing to the integrity of the isolated RNA, such as staining methods and procedures in RNA isolation.

DLG2 variants in patients with pubertal disorders.

PURPOSE: Impaired function of gonadotropin-releasing hormone (GnRH) neurons can cause a phenotypic spectrum ranging from delayed puberty to isolated hypogonadotropic hypogonadism (IHH). We sought to identify a new genetic etiology for these conditions. METHODS: Exome sequencing was performed in an extended family with autosomal dominant, markedly delayed puberty. The effects of the variant were studied in a GnRH neuronal cell line. Variants in the same gene were sought in a large cohort of individuals with IHH. RESULTS: We identified a rare missense variant (F900V) in DLG2 (which encodes PSD-93) that cosegregated with the delayed puberty. The variant decreased GnRH expression in vitro. PSD-93 is an anchoring protein of NMDA receptors, a type of glutamate receptor that has been implicated in the control of puberty in laboratory animals. The F900V variant impaired the interaction between PSD-93 and a known binding partner, Fyn, which phosphorylates NMDA receptors. Variants in DLG2 that also decreased GnRH expression were identified in three unrelated families with IHH. CONCLUSION: The findings indicate that variants in DLG2/PSD-93 cause autosomal dominant delayed puberty and may also contribute to IHH. The findings also suggest that the pathogenesis involves impaired NMDA receptor signaling and consequently decreased GnRH secretion.

Home for a rest: stem cell niche of the postnatal growth plate.

The resting zone houses a group of slowly proliferating 'reserve' chondrocytes and has long been speculated to serve as the stem cell niche of the postnatal growth plate. But are these resting chondrocytes bona fide stem cells? Recent technological advances in lineage tracing and next-generation sequencing have finally allowed researchers to answer this question. Several recent studies have also shed light into the signaling pathways and molecular mechanisms involved in the maintenance of resting chondrocytes, thus providing us with important new insights into the role of the resting zone in the paracrine and endocrine regulation of childhood bone growth.

Persistent Sox9 expression in hypertrophic chondrocytes suppresses transdifferentiation into osteoblasts.

Longitudinal bone growth is driven by endochondral ossification, a process in which cartilage tissue is generated by growth plate chondrocytes and then remodeled into bone by osteoblasts. In the postnatal growth plate, as hypertrophic chondrocytes approach the chondro-osseous junction, they may undergo apoptosis, or directly transdifferentiate into osteoblasts. The molecular mechanisms governing this switch in cell lineage are poorly understood. Here we show that the physiological downregulation of Sox9 in hypertrophic chondrocyte is associated with upregulation of osteoblast-associated genes (such as Mmp13, Cola1, Ibsp) in hypertrophic chondrocytes, before they enter the metaphyseal bone. In transgenic mice that continued to express Sox9 in all cells derived from the chondrocytic lineage, upregulation of these osteoblast-associated genes in the hypertrophic zone failed to occur. Furthermore, lineage tracing experiments showed that, in transgenic mice expressing Sox9, the number of chondrocytes transdifferentiating into osteoblasts was markedly reduced. Collectively, our findings suggest that Sox9 downregulation in hypertrophic chondrocytes promotes expression of osteoblast-associated genes in hypertrophic chondrocytes and promotes the subsequent transdifferentiation of these cells into osteoblasts.

Cartilage-Targeted IGF-1 Treatment to Promote Longitudinal Bone Growth.

Recombinant human growth hormone (GH) is commonly used to treat short stature in children. However, GH treatment has limited efficacy, particularly in severe, non-GH-deficient conditions such as chondrodysplasias, and potential off-target effects. Because short stature results from decreased growth plate chondrogenesis, we developed a cartilage-targeting single-chain human antibody fragment (CaAb) aiming to deliver therapeutic molecules to the growth plate, thereby increasing treatment efficacy while minimizing adverse effects on other tissues. To this end, we created fusion proteins of these CaAbs conjugated with insulin-like growth factor 1 (IGF-1), an endocrine and/or paracrine factor that positively regulates chondrogenesis. These CaAb-IGF-1 fusion proteins retained both cartilage binding and IGF-1 biological activity, and they were able to stimulate bone growth in an organ culture system. Using a GH-deficient (lit) mouse model, we found that subcutaneous injections of these CaAb-IGF-1 fusion proteins increased overall growth plate height without increasing proliferation in kidney cortical cells, suggesting on-target efficacy at the growth plate and less off-target effect on the kidney than IGF-1 alone. Alternate-day injections of these fusion proteins, unlike IGF-1 alone, were sufficient to produce a therapeutic effect. Our findings provide proof of principle that targeting therapeutics to growth plate cartilage can potentially improve treatment for childhood growth disorders.

QRICH1 mutations cause a chondrodysplasia with developmental delay.

In many children with short stature, the etiology of the decreased linear growth remains unknown. We sought to identify the underlying genetic etiology in a patient with short stature, irregular growth plates of the proximal phalanges, developmental delay, and mildly dysmorphic facial features. Exome sequencing identified a de novo, heterozygous, nonsense mutation (c.1606C>T:p.R536X) in QRICH1. In vitro studies confirmed that the mutation impaired expression of the QRICH1 protein. SiRNA-mediated knockdown of Qrich1 in primary mouse epiphyseal chondrocytes caused downregulation of gene expression associated with hypertrophic differentiation. We then identified an unrelated individual with another heterozygous de novo nonsense mutation in QRICH1 who had a similar phenotype. A recently published study identified QRICH1 mutations in three patients with developmental delay, one of whom had short stature. Our findings indicate that QRICH1 mutations cause not only developmental delay but also a chondrodysplasia characterized by diminished linear growth and abnormal growth plate morphology due to impaired growth plate chondrocyte hypertrophic differentiation.

Differential aging of growth plate cartilage underlies differences in bone length and thus helps determine skeletal proportions.

Bones at different anatomical locations vary dramatically in size. For example, human femurs are 20-fold longer than the phalanges in the fingers and toes. The mechanisms responsible for these size differences are poorly understood. Bone elongation occurs at the growth plates and advances rapidly in early life but then progressively slows due to a developmental program termed "growth plate senescence." This developmental program includes declines in cell proliferation and hypertrophy, depletion of cells in all growth plate zones, and extensive underlying changes in the expression of growth-regulating genes. Here, we show evidence that these functional, structural, and molecular senescent changes occur earlier in the growth plates of smaller bones (metacarpals, phalanges) than in the growth plates of larger bones (femurs, tibias) and that this differential aging contributes to the disparities in bone length. We also show evidence that the molecular mechanisms that underlie the differential aging between different bones involve modulation of critical paracrine regulatory pathways, including insulin-like growth factor (Igf), bone morphogenetic protein (Bmp), and Wingless and Int-1 (Wnt) signaling. Taken together, the findings reveal that the striking disparities in the lengths of different bones, which characterize normal mammalian skeletal proportions, is achieved in part by modulating the progression of growth plate senescence.

mir-374-5p, mir-379-5p, and mir-503-5p Regulate Proliferation and Hypertrophic Differentiation of Growth Plate Chondrocytes in Male Rats.

Growth plate chondrocytes undergo sequential differentiation to form the resting zone, the proliferative zone (PZ), and the hypertrophic zone (HZ). The important role of microRNAs (miRNAs) in the growth plate was previously revealed by cartilage-specific ablation of Dicer, an enzyme essential for biogenesis of many miRNAs. To identify specific miRNAs that regulate differentiation of PZ chondrocytes to HZ chondrocytes, we microdissected individual growth plate zones from juvenile rats and performed miRNA profiling using a solution hybridization method and miRNA sequencing. Thirty-four miRNAs were differentially expressed between the PZ and the HZ, and we hypothesized that some of the miRNAs that are preferentially expressed in the PZ may promote proliferation and inhibit hypertrophic differentiation. Consistent with this hypothesis, transfection of inhibitors for four of these miRNAs (mir-369-3p, mir-374-5p, mir-379-5p, and mir-503-5p) decreased proliferation in primary epiphyseal chondrocytes. The inhibitors for three of these miRNAs (mir-374-5p, mir-379-5p, and mir-503-5p) also increased expression of multiple genes that are associated with chondrocyte hypertrophic differentiation. We next hypothesized that preferential expression of these miRNAs in the PZ is driven by the parathyroid hormone-related protein (PTHrP) concentration gradient across the growth plate. Consistent with this hypothesis, treatment of primary chondrocytes with a parathyroid hormone (PTH)/PTHrP receptor agonist, PTH1-34, increased expression of mir-374-5p, mir-379-5p, and mir-503-5p. Taken together, our findings suggest that the PTHrP concentration gradient across the growth plate induces differential expression of mir-374-5p, mir-379-5p, and mir-503-5p between the PZ and the HZ. In the PZ, the higher expression levels of these miRNAs promote proliferation and inhibit hypertrophic differentiation. In the HZ, downregulation of these miRNAs inhibits proliferation and promotes hypertrophic differentiation.

Ezh2 Mutations Found in the Weaver Overgrowth Syndrome Cause a Partial Loss of H3K27 Histone Methyltransferase Activity.

Context: Weaver syndrome is characterized by tall stature, advanced bone age, characteristic facies, and variable intellectual disability. It is caused by heterozygous mutations in enhancer of zeste homolog 2 (EZH2), a histone methyltransferase responsible for histone H3 at lysine 27 (H3K27) trimethylation. However, no early truncating mutations have been identified, suggesting that null mutations do not cause Weaver syndrome. Objective: To test alternative hypotheses that EZH2 variants found in Weaver syndrome cause either a gain of function or a partial loss of function. Design: Exome sequencing was performed in a boy with tall stature, advanced bone age, and mild dysmorphic features. Mutant or wild-type EZH2 protein was expressed in mouse growth plate chondrocytes with or without endogenous EZH2, and enzymatic activity was measured. A mouse model was generated, and histone methylation was assessed in heterozygous and homozygous embryos. Results: A de novo missense EZH2 mutation [c.1876G>A (p.Val626Met)] was identified in the proband. When expressed in growth plate chondrocytes, the mutant protein showed decreased histone methyltransferase activity. A mouse model carrying this EZH2 mutation was generated using CRISPR/Cas9. Homozygotes showed perinatal lethality, whereas heterozygotes were viable, fertile, and showed mild overgrowth. Both homozygous and heterozygous embryos showed decreased H3K27 methylation. Conclusion: We generated a mouse model with the same mutation as our patient, found that it recapitulates the Weaver overgrowth phenotype, and demonstrated that EZH2 mutations found in Weaver syndrome cause a partial loss of function.

Spatial regulation of bone morphogenetic proteins (BMPs) in postnatal articular and growth plate cartilage.

Articular and growth plate cartilage both arise from condensations of mesenchymal cells, but ultimately develop important histological and functional differences. Each is composed of three layers-the superficial, mid and deep zones of articular cartilage and the resting, proliferative and hypertrophic zones of growth plate cartilage. The bone morphogenetic protein (BMP) system plays an important role in cartilage development. A gradient in expression of BMP-related genes has been observed across growth plate cartilage, likely playing a role in zonal differentiation. To investigate the presence of a similar expression gradient in articular cartilage, we used laser capture microdissection (LCM) to separate murine growth plate and articular cartilage from the proximal tibia into their six constituent zones, and used a solution hybridization assay with color-coded probes (nCounter) to quantify mRNAs for 30 different BMP-related genes in each zone. In situ hybridization and immunohistochemistry were then used to confirm spatial expression patterns. Expression gradients for Bmp2 and 6 were observed across growth plate cartilage with highest expression in hypertrophic zone. However, intracellular BMP signaling, assessed by phospho-Smad1/5/8 immunohistochemical staining, appeared to be higher in the proliferative zone and prehypertrophic area than in hypertrophic zone, possibly due to high expression of Smad7, an inhibitory Smad, in the hypertrophic zone. We also found BMP expression gradients across the articular cartilage with BMP agonists primarily expressed in the superficial zone and BMP functional antagonists primarily expressed in the deep zone. Phospho-Smad1/5/8 immunohistochemical staining showed a similar gradient. In combination with previous evidence that BMPs regulate chondrocyte proliferation and differentiation, the current findings suggest that BMP signaling gradients exist across both growth plate and articular cartilage and that these gradients may contribute to the spatial differentiation of chondrocytes in the postnatal endochondral skeleton.