Osteocyte transcriptome mapping identifies a molecular landscape controlling skeletal homeostasis and susceptibility to skeletal disease.

Osteocytes are master regulators of the skeleton. We mapped the transcriptome of osteocytes from different skeletal sites, across age and sexes in mice to reveal genes and molecular programs that control this complex cellular-network. We define an osteocyte transcriptome signature of 1239 genes that distinguishes osteocytes from other cells. 77% have no previously known role in the skeleton and are enriched for genes regulating neuronal network formation, suggesting this programme is important in osteocyte communication. We evaluated 19 skeletal parameters in 733 knockout mouse lines and reveal 26 osteocyte transcriptome signature genes that control bone structure and function. We showed osteocyte transcriptome signature genes are enriched for human orthologs that cause monogenic skeletal disorders (P = 2.4 × 10-22) and are associated with the polygenic diseases osteoporosis (P = 1.8 × 10-13) and osteoarthritis (P = 1.6 × 10-7). Thus, we reveal the molecular landscape that regulates osteocyte network formation and function and establish the importance of osteocytes in human skeletal disease.

Osteoclasts recycle via osteomorphs during RANKL-stimulated bone resorption.

Osteoclasts are large multinucleated bone-resorbing cells formed by the fusion of monocyte/macrophage-derived precursors that are thought to undergo apoptosis once resorption is complete. Here, by intravital imaging, we reveal that RANKL-stimulated osteoclasts have an alternative cell fate in which they fission into daughter cells called osteomorphs. Inhibiting RANKL blocked this cellular recycling and resulted in osteomorph accumulation. Single-cell RNA sequencing showed that osteomorphs are transcriptionally distinct from osteoclasts and macrophages and express a number of non-canonical osteoclast genes that are associated with structural and functional bone phenotypes when deleted in mice. Furthermore, genetic variation in human orthologs of osteomorph genes causes monogenic skeletal disorders and associates with bone mineral density, a polygenetic skeletal trait. Thus, osteoclasts recycle via osteomorphs, a cell type involved in the regulation of bone resorption that may be targeted for the treatment of skeletal diseases.

A molecular quantitative trait locus map for osteoarthritis.

Osteoarthritis causes pain and functional disability for over 500 million people worldwide. To develop disease-stratifying tools and modifying therapies, we need a better understanding of the molecular basis of the disease in relevant tissue and cell types. Here, we study primary cartilage and synovium from 115 patients with osteoarthritis to construct a deep molecular signature map of the disease. By integrating genetics with transcriptomics and proteomics, we discover molecular trait loci in each tissue type and omics level, identify likely effector genes for osteoarthritis-associated genetic signals and highlight high-value targets for drug development and repurposing. These findings provide insights into disease aetiopathology, and offer translational opportunities in response to the global clinical challenge of osteoarthritis.

Accelerating functional gene discovery in osteoarthritis.

Osteoarthritis causes debilitating pain and disability, resulting in a considerable socioeconomic burden, yet no drugs are available that prevent disease onset or progression. Here, we develop, validate and use rapid-throughput imaging techniques to identify abnormal joint phenotypes in randomly selected mutant mice generated by the International Knockout Mouse Consortium. We identify 14 genes with functional involvement in osteoarthritis pathogenesis, including the homeobox gene Pitx1, and functionally characterize 6 candidate human osteoarthritis genes in mouse models. We demonstrate sensitivity of the methods by identifying age-related degenerative joint damage in wild-type mice. Finally, we phenotype previously generated mutant mice with an osteoarthritis-associated polymorphism in the Dio2 gene by CRISPR/Cas9 genome editing and demonstrate a protective role in disease onset with public health implications. We hope this expanding resource of mutant mice will accelerate functional gene discovery in osteoarthritis and offer drug discovery opportunities for this common, incapacitating chronic disease.

IGSF1 Deficiency Results in Human and Murine Somatotrope Neurosecretory Hyperfunction.

CONTEXT: The X-linked immunoglobulin superfamily, member 1 (IGSF1), gene is highly expressed in the hypothalamus and in pituitary cells of the POU1F1 lineage. Human loss-of-function mutations in IGSF1 cause central hypothyroidism, hypoprolactinemia, and macroorchidism. Additionally, most affected adults exhibit higher than average IGF-1 levels and anecdotal reports describe acromegaloid features in older subjects. However, somatotrope function has not yet been formally evaluated in this condition. OBJECTIVE: We aimed to evaluate the role of IGSF1 in human and murine somatotrope function. PATIENTS, DESIGN, AND SETTING: We evaluated 21 adult males harboring hemizygous IGSF1 loss-of-function mutations for features of GH excess, in an academic clinical setting. MAIN OUTCOME MEASURES: We compared biochemical and tissue markers of GH excess in patients and controls, including 24-hour GH profile studies in 7 patients. Parallel studies were undertaken in male Igsf1-deficient mice and wild-type littermates. RESULTS: IGSF1-deficient adult male patients demonstrated acromegaloid facial features with increased head circumference as well as increased finger soft-tissue thickness. Median serum IGF-1 concentrations were elevated, and 24-hour GH profile studies confirmed 2- to 3-fold increased median basal, pulsatile, and total GH secretion. Male Igsf1-deficient mice also demonstrated features of GH excess with increased lean mass, organ size, and skeletal dimensions and elevated mean circulating IGF-1 and pituitary GH levels. CONCLUSIONS: We demonstrate somatotrope neurosecretory hyperfunction in IGSF1-deficient humans and mice. These observations define a hitherto uncharacterized role for IGSF1 in somatotropes and indicate that patients with IGSF1 mutations should be evaluated for long-term consequences of increased GH exposure.

PYY is a negative regulator of bone mass and strength.

OBJECTIVE: Bone loss in anorexia nervosa and following bariatric surgery is associated with an elevated circulating concentration of the gastrointestinal, anorexigenic hormone, peptide YY (PYY). Selective deletion of the PYY receptor Y1R in osteoblasts or Y2R in the hypothalamus results in high bone mass, but deletion of PYY in mice has resulted in conflicting skeletal phenotypes leading to uncertainty regarding its role in the regulation of bone mass. As PYY analogs are under development for treatment of obesity, we aimed to clarify the relationship between PYY and bone mass. METHODS: The skeletal phenotype of Pyy knockout (KO) mice was investigated during growth (postnatal day P14) and adulthood (P70 and P186) using X-ray microradiography, micro-CT, back-scattered electron scanning electron microscopy (BSE-SEM), histomorphometry and biomechanical testing. RESULTS: Bones from juvenile and Pyy KO mice were longer (P < 0.001), with decreased bone mineral content (P < 0.001). Whereas, bones from adult Pyy KO mice had increased bone mineral content (P < 0.05) with increased mineralisation of both cortical (P < 0.001) and trabecular (P < 0.001) compartments. Long bones from adult Pyy KO mice were stronger (maximum load P < 0.001), with increased stiffness (P < 0.01) and toughness (P < 0.05) compared to wild-type (WT) control mice despite increased cortical vascularity and porosity (P < 0.001). The increased bone mass and strength in Pyy KO mice resulted from increases in trabecular (P < 0.01) and cortical bone formation (P < 0.05). CONCLUSIONS: These findings demonstrate that PYY acts as a negative regulator of osteoblastic bone formation, implicating increased PYY levels in the pathogenesis of bone loss during anorexia or following bariatric surgery.

Author Correction: An atlas of genetic influences on osteoporosis in humans and mice.

In the version of this article initially published, in Fig. 5a, the data in the right column of 'DAAM2 gRNA1' were incorrectly plotted as circles indicating 'untreated' rather than as squares indicating 'treated'. The error has been corrected in the HTML and PDF versions of the article.

Slc20a2, Encoding the Phosphate Transporter PiT2, Is an Important Genetic Determinant of Bone Quality and Strength.

Osteoporosis is characterized by low bone mineral density (BMD) and fragility fracture and affects over 200 million people worldwide. Bone quality describes the material properties that contribute to strength independently of BMD, and its quantitative analysis is a major priority in osteoporosis research. Tissue mineralization is a fundamental process requiring calcium and phosphate transporters. Here we identify impaired bone quality and strength in Slc20a2-/- mice lacking the phosphate transporter SLC20A2. Juveniles had abnormal endochondral and intramembranous ossification, decreased mineral accrual, and short stature. Adults exhibited only small reductions in bone mass and mineralization but a profound impairment of bone strength. Bone quality was severely impaired in Slc20a2-/- mice: yield load (-2.3 SD), maximum load (-1.7 SD), and stiffness (-2.7 SD) were all below values predicted from their bone mineral content as determined in a cohort of 320 wild-type controls. These studies identify Slc20a2 as a physiological regulator of tissue mineralization and highlight its critical role in the determination of bone quality and strength. © 2019 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals Inc.

Genetic and Pharmacological Targeting of Transcriptional Repression in Resistance to Thyroid Hormone Alpha.

Background: Thyroid hormones act in bone and cartilage via thyroid hormone receptor alpha (TRα). In the absence of triiodothyronine (T3), TRα interacts with co-repressors, including nuclear receptor co-repressor-1 (NCoR1), which recruit histone deacetylases (HDACs) and mediate transcriptional repression. Dominant-negative mutations of TRα cause resistance to thyroid hormone alpha (RTHα; OMIM 614450), characterized by excessive repression of T3 target genes leading to delayed skeletal development, growth retardation, and bone dysplasia. Treatment with thyroxine has been of limited benefit, even in mildly affected individuals, and there is a need for new therapeutic strategies. It was hypothesized that (i) the skeletal manifestations of RTHα are mediated by the persistent TRα/NCoR1/HDAC repressor complex containing mutant TRα, and (ii) treatment with the HDAC inhibitor suberoylanilide hydroxamic acid (SAHA) would ameliorate these manifestations. Methods: The skeletal phenotypes of (i) Thra1PV/+ mice, a well characterized model of RTHα; (ii) Ncor1ΔID/ΔID mice, which express an NCoR1 mutant that fails to interact with TRα; and (iii) Thra1PV/+Ncor1ΔID/ΔID double-mutant adult mice were determined. Wild-type, Thra1PV/+, Ncor1ΔID/ΔID, and Thra1PV/+Ncor1ΔID/ΔID double-mutant mice were also treated with SAHA to determine whether HDAC inhibition results in amelioration of skeletal abnormalities. Results:Thra1PV/+ mice had a severe skeletal dysplasia, characterized by short stature, abnormal bone morphology, and increased bone mineral content. Despite normal bone length, Ncor1ΔID/ΔID mice displayed increased cortical bone mass, mineralization, and strength. Thra1PV/+Ncor1ΔID/ΔID double-mutant mice displayed only a small improvement of skeletal abnormalities compared to Thra1PV/+ mice. Treatment with SAHA to inhibit histone deacetylation had no beneficial or detrimental effects on bone structure, mineralization, or strength in wild-type or mutant mice. Conclusions: These studies indicate treatment with SAHA is unlikely to improve the skeletal manifestations of RTHα. Nevertheless, the findings (i) confirm that TRα1 has a critical role in the regulation of skeletal development and adult bone mass, (ii) suggest a physiological role for alternative co-repressors that interact with TR in skeletal cells, and (iii) demonstrate a novel role for NCoR1 in the regulation of adult bone mass and strength.

Quantitative X-Ray Imaging of Mouse Bone by Faxitron.

This chapter describes the use of point projection digital microradiography for rapid imaging and quantitation of bone mineral content in mice.

An atlas of genetic influences on osteoporosis in humans and mice.

Osteoporosis is a common aging-related disease diagnosed primarily using bone mineral density (BMD). We assessed genetic determinants of BMD as estimated by heel quantitative ultrasound in 426,824 individuals, identifying 518 genome-wide significant loci (301 novel), explaining 20% of its variance. We identified 13 bone fracture loci, all associated with estimated BMD (eBMD), in ~1.2 million individuals. We then identified target genes enriched for genes known to influence bone density and strength (maximum odds ratio (OR) = 58, P = 1 × 10-75) from cell-specific features, including chromatin conformation and accessible chromatin sites. We next performed rapid-throughput skeletal phenotyping of 126 knockout mice with disruptions in predicted target genes and found an increased abnormal skeletal phenotype frequency compared to 526 unselected lines (P < 0.0001). In-depth analysis of one gene, DAAM2, showed a disproportionate decrease in bone strength relative to mineralization. This genetic atlas provides evidence linking associated SNPs to causal genes, offers new insight into osteoporosis pathophysiology, and highlights opportunities for drug development.

Type 2 deiodinase polymorphism causes ER stress and hypothyroidism in the brain.

Levothyroxine (LT4) is a form of thyroid hormone used to treat hypothyroidism. In the brain, T4 is converted to the active form T3 by type 2 deiodinase (D2). Thus, it is intriguing that carriers of the Thr92Ala polymorphism in the D2 gene (DIO2) exhibit clinical improvement when liothyronine (LT3) is added to LT4 therapy. Here, we report that D2 is a cargo protein in ER Golgi intermediary compartment (ERGIC) vesicles, recycling between ER and Golgi. The Thr92-to-Ala substitution (Ala92-D2) caused ER stress and activated the unfolded protein response (UPR). Ala92-D2 accumulated in the trans-Golgi and generated less T3, which was restored by eliminating ER stress with the chemical chaperone 4-phenyl butyric acid (4-PBA). An Ala92-Dio2 polymorphism-carrying mouse exhibited UPR and hypothyroidism in distinct brain areas. The mouse refrained from physical activity, slept more, and required additional time to memorize objects. Enhancing T3 signaling in the brain with LT3 improved cognition, whereas restoring proteostasis with 4-PBA eliminated the Ala92-Dio2 phenotype. In contrast, primary hypothyroidism intensified the Ala92-Dio2 phenotype, with only partial response to LT4 therapy. Disruption of cellular proteostasis and reduced Ala92-D2 activity may explain the failure of LT4 therapy in carriers of Thr92Ala-DIO2.

Identification of 153 new loci associated with heel bone mineral density and functional involvement of GPC6 in osteoporosis.

Osteoporosis is a common disease diagnosed primarily by measurement of bone mineral density (BMD). We undertook a genome-wide association study (GWAS) in 142,487 individuals from the UK Biobank to identify loci associated with BMD as estimated by quantitative ultrasound of the heel. We identified 307 conditionally independent single-nucleotide polymorphisms (SNPs) that attained genome-wide significance at 203 loci, explaining approximately 12% of the phenotypic variance. These included 153 previously unreported loci, and several rare variants with large effect sizes. To investigate the underlying mechanisms, we undertook (1) bioinformatic, functional genomic annotation and human osteoblast expression studies; (2) gene-function prediction; (3) skeletal phenotyping of 120 knockout mice with deletions of genes adjacent to lead independent SNPs; and (4) analysis of gene expression in mouse osteoblasts, osteocytes and osteoclasts. The results implicate GPC6 as a novel determinant of BMD, and also identify abnormal skeletal phenotypes in knockout mice associated with a further 100 prioritized genes.

Pitx1 determines characteristic hindlimb morphologies in cartilage micromass culture.

The shapes of homologous skeletal elements in the vertebrate forelimb and hindlimb are distinct, with each element exquisitely adapted to their divergent functions. Many of the signals and signalling pathways responsible for patterning the developing limb bud are common to both forelimb and hindlimb. How disparate morphologies are generated from common signalling inputs during limb development remains poorly understood. We show that, similar to what has been shown in the chick, characteristic differences in mouse forelimb and hindlimb cartilage morphology are maintained when chondrogenesis proceeds in vitro away from the endogenous limb bud environment. Chondrogenic nodules that form in high-density micromass cultures derived from forelimb and hindlimb buds are consistently different in size and shape. We described analytical tools we have developed to quantify these differences in nodule morphology and demonstrate that characteristic hindlimb nodule morphology is lost in the absence of the hindlimb-restricted limb modifier gene Pitx1. Furthermore, we show that ectopic expression of Pitx1 in the forelimb is sufficient to generate nodule patterns characteristic of the hindlimb. We also demonstrate that hindlimb cells are less adhesive to the tissue culture substrate and, within the limb environment, to the extracellular matrix and to each other. These results reveal autonomously programmed differences in forelimb and hindlimb cartilage precursors of the limb skeleton are controlled, at least in part, by Pitx1 and suggest this has an important role in generating distinct limb-type morphologies. Our results demonstrate that the micromass culture system is ideally suited to study cues governing morphogenesis of limb skeletal elements in a simple and experimentally tractable in vitro system that reflects in vivo potential.

Tbx5 Buffers Inherent Left/Right Asymmetry Ensuring Symmetric Forelimb Formation.

The forelimbs and hindlimbs of vertebrates are bilaterally symmetric. The mechanisms that ensure symmetric limb formation are unknown but they can be disrupted in disease. In Holt-Oram Syndrome (HOS), caused by mutations in TBX5, affected individuals have left-biased upper/forelimb defects. We demonstrate a role for the transcription factor Tbx5 in ensuring the symmetric formation of the left and right forelimb. In our mouse model, bilateral hypomorphic levels of Tbx5 produces asymmetric forelimb defects that are consistently more severe in the left limb than the right, phenocopying the left-biased limb defects seen in HOS patients. In Tbx hypomorphic mutants maintained on an INV mutant background, with situs inversus, the laterality of defects is reversed. Our data demonstrate an early, inherent asymmetry in the left and right limb-forming regions and that threshold levels of Tbx5 are required to overcome this asymmetry to ensure symmetric forelimb formation.

Mutations in mouse Ift144 model the craniofacial, limb and rib defects in skeletal ciliopathies.

Mutations in components of the intraflagellar transport (IFT) machinery required for assembly and function of the primary cilium cause a subset of human ciliopathies characterized primarily by skeletal dysplasia. Recently, mutations in the IFT-A gene IFT144 have been described in patients with Sensenbrenner and Jeune syndromes, which are associated with short ribs and limbs, polydactyly and craniofacial defects. Here, we describe an N-ethyl-N-nitrosourea-derived mouse mutant with a hypomorphic missense mutation in the Ift144 gene. The mutant twinkle-toes (Ift144(twt)) phenocopies a number of the skeletal and craniofacial anomalies seen in patients with human skeletal ciliopathies. Like other IFT-A mouse mutants, Ift144 mutant embryos display a generalized ligand-independent expansion of hedgehog (Hh) signalling, in spite of defective ciliogenesis and an attenuation of the ability of mutant cells to respond to upstream stimulation of the pathway. This enhanced Hh signalling is consistent with cleft palate and polydactyly phenotypes in the Ift144(twt) mutant, although extensive rib branching, fusion and truncation phenotypes correlate with defects in early somite patterning and may reflect contributions from multiple signalling pathways. Analysis of embryos harbouring a second allele of Ift144 which represents a functional null, revealed a dose-dependent effect on limb outgrowth consistent with the short-limb phenotypes characteristic of these ciliopathies. This allelic series of mouse mutants provides a unique opportunity to uncover the underlying mechanistic basis of this intriguing subset of ciliopathies.

The molecular regulation of vertebrate limb patterning.

The limb has long been considered a paradigm for organogenesis because of its simplicity and ease of manipulation. However, it has become increasingly clear that the processes required to produce a perfectly formed limb involve complex molecular interactions across all three axes of limb development. Old models have evolved with acquisition of molecular knowledge, and in more recent times mathematical modeling approaches have been invoked to explain the precise spatio-temporal regulation of gene networks that coordinate limb patterning and outgrowth. This review focuses on recent advances in our understanding of vertebrate limb development, highlighting the signaling interactions required to lay down the pattern on which the processes of differentiation will act to ultimately produce the final limb.