Characterization of two ENU-induced mutations affecting mouse skeletal morphology

Interspecies transcriptomics identify genes that underlie disproportionate foot growth in jerboas

Despite the great diversity of vertebrate limb proportion and our deep understanding of the genetic mechanisms that drive skeletal elongation, little is known about how individual bones reach different lengths in any species. Here, we directly compare the transcriptomes of homologous growth cartilages of the mouse (Mus musculus) and bipedal jerboa (Jaculus jaculus), the latter of which has "mouse-like" arms but extremely long metatarsals of the feet. Intersecting gene-expression differences in metatarsals and forearms of the two species revealed that about 10% of orthologous genes are associated with the disproportionately rapid elongation of neonatal jerboa feet. These include genes and enriched pathways not previously associated with endochondral elongation as well as those that might diversify skeletal proportion in addition to their known requirements for bone growth throughout the skeleton. We also identified transcription regulators that might act as "nodes" for sweeping differences in genome expression between species. Among these, Shox2, which is necessary for proximal limb elongation, has gained expression in jerboa metatarsals where it has not been detected in other vertebrates. We show that Shox2 is sufficient to increase mouse distal limb length, and a nearby putative cis-regulatory region is preferentially accessible in jerboa metatarsals. In addition to mechanisms that might directly promote growth, we found evidence that jerboa foot elongation may occur in part by de-repressing latent growth potential. The genes and pathways that we identified here provide a framework to understand the modular genetic control of skeletal growth and the remarkable malleability of vertebrate limb proportion.

Genetic animal models of scoliosis: A systematical review

Scoliosis is a complex disease with undetermined pathogenesis and has a strong relationship with genetics. Models of scoliosis in animals have been established for better comprehending its pathogenesis and treatment. In this review, we searched all the genetic animal models with body curvature in databases, and reviewed the related genes and scoliosis types. Meanwhile, we also summarized the pathogenesis of scoliosis reported so far. Summarizing the positive phenotypic animal models contributes to a better understanding on the pathogenesis of scoliosis and facilitates the selection of experimental models when a possible pathogenic factor is concerned.

Prioritization of Osteoporosis-Associated Genome-wide Association Study (GWAS) Single-Nucleotide Polymorphisms (SNPs) Using Epigenomics and Transcriptomics

Genetic risk factors for osteoporosis, a prevalent disease associated with aging, have been examined in many genome-wide association studies (GWASs). A major challenge is to prioritize transcription-regulatory GWAS-derived variants that are likely to be functional. Given the critical role of epigenetics in gene regulation, we have used an unusual epigenetics-based and transcription-based approach to identify some of the credible regulatory single-nucleotide polymorphisms (SNPs) relevant to osteoporosis from 38 reported bone mineral density (BMD) GWASs. Using Roadmap databases, we prioritized SNPs based upon their overlap with strong enhancer or promoter chromatin preferentially in osteoblasts relative to 12 heterologous cell culture types. We also required that these SNPs overlap open chromatin (Deoxyribonuclease I [DNaseI]-hypersensitive sites) and DNA sequences predicted to bind to osteoblast-relevant transcription factors in an allele-specific manner. From >50,000 GWAS-derived SNPs, we identified 14 novel and credible regulatory SNPs (Tier-1 SNPs) for osteoporosis risk. Their associated genes, BICC1, LGR4, DAAM2, NPR3, or HMGA2, are involved in osteoblastogenesis or bone homeostasis and regulate cell signaling or enhancer function. Four of these genes are preferentially expressed in osteoblasts. BICC1, LGR4, and DAAM2 play important roles in canonical Wnt signaling, a pathway critical for bone formation and repair. The transcription factors predicted to bind to the Tier-1 SNP-containing DNA sequences also have bone-related functions. We present evidence that some of the Tier-1 SNPs exert their effects on BMD risk indirectly through little-studied long noncoding RNA (lncRNA) genes, which may, in turn, control the nearby bone-related protein-encoding gene. Our study illustrates a method to identify novel BMD-related causal regulatory SNPs for future study and to prioritize candidate regulatory GWAS-derived SNPs, in general. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Bi-allelic Loss-of-Function Mutations in the NPR-C Receptor Result in Enhanced Growth and Connective Tissue Abnormalities

The natriuretic peptide signaling pathway has been implicated in many cellular processes, including endochondral ossification and bone growth. More precisely, different mutations in the NPR-B receptor and the CNP ligand have been identified in individuals with either short or tall stature. In this study we show that the NPR-C receptor (encoded by NPR3) is also important for the regulation of linear bone growth. We report four individuals, originating from three different families, with a phenotype characterized by tall stature, long digits, and extra epiphyses in the hands and feet. In addition, aortic dilatation was observed in two of these families. In each affected individual, we identified a bi-allelic loss-of-function mutation in NPR3. The missense mutations (c.442T>C [p.Ser148Pro] and c.1088A>T [p.Asp363Val]) resulted in intracellular retention of the NPR-C receptor and absent localization on the plasma membrane, whereas the nonsense mutation (c.1524delC [p.Tyr508^∗]) resulted in nonsense-mediated mRNA decay. Biochemical analysis of plasma from two affected and unrelated individuals revealed a reduced NTproNP/NP ratio for all ligands and also high cGMP levels. These data strongly suggest a reduced clearance of natriuretic peptides by the defective NPR-C receptor and consequently increased activity of the NPR-A/B receptors. In conclusion, this study demonstrates that loss-of-function mutations in NPR3 result in increased NPR-A/B signaling activity and cause a phenotype marked by enhanced bone growth and cardiovascular abnormalities.

Osteocrin, a peptide secreted from the heart and other tissues, contributes to cranial osteogenesis and chondrogenesis in zebrafish

The heart is an endocrine organ, as cardiomyocytes (CMs) secrete natriuretic peptide (NP) hormones. Since the discovery of NPs, no other peptide hormones that affect remote organs have been identified from the heart. We identified osteocrin (Ostn) as an osteogenesis/chondrogenesis regulatory hormone secreted from CMs in zebrafish. ostn mutant larvae exhibit impaired membranous and chondral bone formation. The impaired bones were recovered by CM-specific overexpression of OSTN. We analyzed the parasphenoid (ps) as a representative of membranous bones. In the shortened ps of ostn morphants, nuclear Yap1/Wwtr1-dependent transcription was increased, suggesting that Ostn might induce the nuclear export of Yap1/Wwtr1 in osteoblasts. Although OSTN is proposed to bind to NPR3 (clearance receptor for NPs) to enhance the binding of NPs to NPR1 or NPR2, OSTN enhanced C-type NP (CNP)-dependent nuclear export of YAP1/WWTR1 of cultured mouse osteoblasts stimulated with saturable CNP. OSTN might therefore activate unidentified receptors that augment protein kinase G signaling mediated by a CNP-NPR2 signaling axis. These data demonstrate that Ostn secreted from the heart contributes to bone formation as an endocrine hormone.

Mice with an N-Ethyl-N-Nitrosourea (ENU) Induced Tyr209Asn Mutation in Natriuretic Peptide Receptor 3 (NPR3) Provide a Model for Kyphosis Associated with Activation of the MAPK Signaling Pathway

Non-syndromic kyphosis is a common disorder that is associated with significant morbidity and has a strong genetic involvement; however, the causative genes remain to be identified, as such studies are hampered by genetic heterogeneity, small families and various modes of inheritance. To overcome these limitations, we investigated 12 week old progeny of mice treated with the chemical mutagen N-ethyl-N-nitrosourea (ENU) using phenotypic assessments including dysmorphology, radiography, and dual-energy X-ray absorptiometry. This identified a mouse with autosomal recessive kyphosis (KYLB). KYLB mice, when compared to unaffected littermates, had: thoraco-lumbar kyphosis, larger vertebrae, and increased body length and increased bone area. In addition, female KYLB mice had increases in bone mineral content and plasma alkaline phosphatase activity. Recombination mapping localized the Kylb locus to a 5.5Mb region on chromosome 15A1, which contained 51 genes, including the natriuretic peptide receptor 3 (Npr3) gene. DNA sequence analysis of Npr3 identified a missense mutation, Tyr209Asn, which introduced an N-linked glycosylation consensus sequence. Expression of wild-type NPR3 and the KYLB-associated Tyr209Asn NPR3 mutant in COS-7 cells demonstrated the mutant to be associated with abnormal N-linked glycosylation and retention in the endoplasmic reticulum that resulted in its absence from the plasma membrane. NPR3 is a decoy receptor for C-type natriuretic peptide (CNP), which also binds to NPR2 and stimulates mitogen-activated protein kinase (MAPK) signaling, thereby increasing the number and size of hypertrophic chondrocytes. Histomorphometric analysis of KYLB vertebrae and tibiae showed delayed endochondral ossification and expansion of the hypertrophic zones of the growth plates, and immunohistochemistry revealed increased p38 MAPK phosphorylation throughout the growth plates of KYLB vertebrae. Thus, we established a model of kyphosis due to a novel NPR3 mutation, in which loss of plasma membrane NPR3 expression results in increased MAPK pathway activation, causing elongation of the vertebrae and resulting in kyphosis.

Micro CT Analysis of Spine Architecture in a Mouse Model of Scoliosis

OBJECTIVE

Mice homozygous for targeted deletion of the gene encoding fibroblast growth factor receptor 3 (FGFR3(-/-)) develop kyphoscoliosis by 2 months of age. The first objective of this study was to use high resolution X-ray to characterize curve progression in vivo and micro CT to quantify spine architecture ex vivo in FGFR3(-/-) mice. The second objective was to determine if slow release of the bone anabolic peptide parathyroid hormone related protein (PTHrP-1-34) from a pellet placed adjacent to the thoracic spine could inhibit progressive kyphoscoliosis.

MATERIALS AND METHODS

Pellets loaded with placebo or PTHrP-1-34 were implanted adjacent to the thoracic spine of 1-month-old FGFR3(-/-) mice obtained from in house breeding. X rays were captured at monthly intervals up to 4 months to quantify curve progression using the Cobb method. High resolution post-mortem scans of FGFR3(-/-) and FGFR3(+/+) spines, from C5/6 to L4/5, were captured to evaluate the 3D structure, rotation, and micro-architecture of the affected vertebrae. Un-decalcified and decalcified histology were performed on the apical and adjacent vertebrae of FGFR3(-/-) spines, and the corresponding vertebrae from FGFR3(+/+) spines.

RESULTS

The mean Cobb angle was significantly greater at all ages in FGFR3(-/-) mice compared with wild type mice and appeared to stabilize around skeletal maturity at 4 months. 3D reconstructions of the thoracic spine of 4-month-old FGFR3(-/-) mice treated with PTHrP-1-34 revealed correction of left/right asymmetry, vertebral rotation, and lateral displacement compared with mice treated with placebo. Histologic analysis of the apical vertebrae confirmed correction of the asymmetry in PTHrP-1-34 treated mice, in the absence of any change in bone volume, and a significant reduction in the wedging of intervertebral disks (IVD) seen in placebo treated mice.

CONCLUSION

Local treatment of the thoracic spine of juvenile FGFR3(-/-) mice with a bone anabolic agent inhibited progression of scoliosis, but with little impact on kyphosis. The significant improvement in IVD integrity suggests PTHrP-1-34 might also be considered as a therapeutic agent for degenerative disk disorders.

High-throughput screening of mouse gene knockouts identifies established and novel skeletal phenotypes

Screening gene function in vivo is a powerful approach to discover novel drug targets. We present high-throughput screening (HTS) data for 3 762 distinct global gene knockout (KO) mouse lines with viable adult homozygous mice generated using either gene-trap or homologous recombination technologies. Bone mass was determined from DEXA scans of male and female mice at 14 weeks of age and by microCT analyses of bones from male mice at 16 weeks of age. Wild-type (WT) cagemates/littermates were examined for each gene KO. Lethality was observed in an additional 850 KO lines. Since primary HTS are susceptible to false positive findings, additional cohorts of mice from KO lines with intriguing HTS bone data were examined. Aging, ovariectomy, histomorphometry and bone strength studies were performed and possible non-skeletal phenotypes were explored. Together, these screens identified multiple genes affecting bone mass: 23 previously reported genes (Calcr, Cebpb, Crtap, Dcstamp, Dkk1, Duoxa2, Enpp1, Fgf23, Kiss1/Kiss1r, Kl (Klotho), Lrp5, Mstn, Neo1, Npr2, Ostm1, Postn, Sfrp4, Slc30a5, Slc39a13, Sost, Sumf1, Src, Wnt10b), five novel genes extensively characterized (Cldn18, Fam20c, Lrrk1, Sgpl1, Wnt16), five novel genes with preliminary characterization (Agpat2, Rassf5, Slc10a7, Slc26a7, Slc30a10) and three novel undisclosed genes coding for potential osteoporosis drug targets.