A perspective on muscle phenotyping in musculoskeletal research.

Musculoskeletal research should synergistically investigate bone and muscle to inform approaches for maintaining mobility and to avoid bone fractures. The relationship between sarcopenia and osteoporosis, integrated in the term 'osteosarcopenia', is underscored by the close association shown between these two conditions in many studies, whereby one entity emerges as a predictor of the other. In a recent workshop of Working Group (WG) 2 of the EU Cooperation in Science and Technology (COST) Action 'Genomics of MusculoSkeletal traits Translational Network' (GEMSTONE) consortium (CA18139), muscle characterization was highlighted as being important, but currently under-recognized in the musculoskeletal field. Here, we summarize the opinions of the Consortium and research questions around translational and clinical musculoskeletal research, discussing muscle phenotyping in human experimental research and in two animal models: zebrafish and mouse.

Osteomodulin downregulation is associated with osteoarthritis development.

Abnormal subchondral bone remodeling leading to sclerosis is a main feature of osteoarthritis (OA), and osteomodulin (OMD), a proteoglycan involved in extracellular matrix mineralization, is associated with the sclerotic phenotype. However, the functions of OMD remain poorly understood, specifically in vivo. We used Omd knockout and overexpressing male mice and mutant zebrafish to study its roles in bone and cartilage metabolism and in the development of OA. The expression of Omd is deeply correlated with bone and cartilage microarchitectures affecting the bone volume and the onset of subchondral bone sclerosis and spontaneous cartilage lesions. Mechanistically, OMD binds to RANKL and inhibits osteoclastogenesis, thus controlling the balance of bone remodeling. In conclusion, OMD is a key factor in subchondral bone sclerosis associated with OA. It participates in bone and cartilage homeostasis by acting on the regulation of osteoclastogenesis. Targeting OMD may be a promising new and personalized approach for OA.

Finding the genes for fragile bones.

Combining transcriptomic data with the analysis of large genome-wide association studies helps identify genes that are likely important for regulating bone mineral density.

The genetic overlap between osteoporosis and craniosynostosis.

Osteoporosis is the most prevalent bone condition in the ageing population. This systemic disease is characterized by microarchitectural deterioration of bone, leading to increased fracture risk. In the past 15 years, genome-wide association studies (GWAS), have pinpointed hundreds of loci associated with bone mineral density (BMD), helping elucidate the underlying molecular mechanisms and genetic architecture of fracture risk. However, the challenge remains in pinpointing causative genes driving GWAS signals as a pivotal step to drawing the translational therapeutic roadmap. Recently, a skull BMD-GWAS uncovered an intriguing intersection with craniosynostosis, a congenital anomaly due to premature suture fusion in the skull. Here, we recapitulate the genetic contribution to both osteoporosis and craniosynostosis, describing the biological underpinnings of this overlap and using zebrafish models to leverage the functional investigation of genes associated with skull development and systemic skeletal homeostasis.

Functional Validation of Osteoporosis Genetic Findings Using Small Fish Models.

The advancement of human genomics has revolutionized our understanding of the genetic architecture of many skeletal diseases, including osteoporosis. However, interpreting results from human association studies remains a challenge, since index variants often reside in non-coding regions of the genome and do not possess an obvious regulatory function. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary, such as the one offered by animal models. These models enable us to identify causal mechanisms, clarify the underlying biology, and apply interventions. Over the past several decades, small teleost fishes, mostly zebrafish and medaka, have emerged as powerful systems for modeling the genetics of human diseases. Due to their amenability to genetic intervention and the highly conserved genetic and physiological features, fish have become indispensable for skeletal genomic studies. The goal of this review is to summarize the evidence supporting the utility of Zebrafish (Danio rerio) for accelerating our understanding of human skeletal genomics and outlining the remaining gaps in knowledge. We provide an overview of zebrafish skeletal morphophysiology and gene homology, shedding light on the advantages of human skeletal genomic exploration and validation. Knowledge of the biology underlying osteoporosis through animal models will lead to the translation into new, better and more effective therapeutic approaches.

Dominating lengthscales of zebrafish collective behaviour.

Collective behaviour in living systems is observed across many scales, from bacteria to insects, to fish shoals. Zebrafish have emerged as a model system amenable to laboratory study. Here we report a three-dimensional study of the collective dynamics of fifty zebrafish. We observed the emergence of collective behaviour changing between ordered to randomised, upon adaptation to new environmental conditions. We quantify the spatial and temporal correlation functions of the fish and identify two length scales, the persistence length and the nearest neighbour distance, that capture the essence of the behavioural changes. The ratio of the two length scales correlates robustly with the polarisation of collective motion that we explain with a reductionist model of self-propelled particles with alignment interactions.

Perspective of the GEMSTONE Consortium on Current and Future Approaches to Functional Validation for Skeletal Genetic Disease Using Cellular, Molecular and Animal-Modeling Techniques.

The availability of large human datasets for genome-wide association studies (GWAS) and the advancement of sequencing technologies have boosted the identification of genetic variants in complex and rare diseases in the skeletal field. Yet, interpreting results from human association studies remains a challenge. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary. Multiple unknowns exist for putative causal genes, including cellular localization of the molecular function. Intermediate traits ("endophenotypes"), e.g. molecular quantitative trait loci (molQTLs), are needed to identify mechanisms of underlying associations. Furthermore, index variants often reside in non-coding regions of the genome, therefore challenging for interpretation. Knowledge of non-coding variance (e.g. ncRNAs), repetitive sequences, and regulatory interactions between enhancers and their target genes is central for understanding causal genes in skeletal conditions. Animal models with deep skeletal phenotyping and cell culture models have already facilitated fine mapping of some association signals, elucidated gene mechanisms, and revealed disease-relevant biology. However, to accelerate research towards bridging the current gap between association and causality in skeletal diseases, alternative in vivo platforms need to be used and developed in parallel with the current -omics and traditional in vivo resources. Therefore, we argue that as a field we need to establish resource-sharing standards to collectively address complex research questions. These standards will promote data integration from various -omics technologies and functional dissection of human complex traits. In this mission statement, we review the current available resources and as a group propose a consensus to facilitate resource sharing using existing and future resources. Such coordination efforts will maximize the acquisition of knowledge from different approaches and thus reduce redundancy and duplication of resources. These measures will help to understand the pathogenesis of osteoporosis and other skeletal diseases towards defining new and more efficient therapeutic targets.

Bone Phenotyping Approaches in Human, Mice and Zebrafish - Expert Overview of the EU Cost Action GEMSTONE ("GEnomics of MusculoSkeletal traits TranslatiOnal NEtwork").

A synoptic overview of scientific methods applied in bone and associated research fields across species has yet to be published. Experts from the EU Cost Action GEMSTONE ("GEnomics of MusculoSkeletal Traits translational Network") Working Group 2 present an overview of the routine techniques as well as clinical and research approaches employed to characterize bone phenotypes in humans and selected animal models (mice and zebrafish) of health and disease. The goal is consolidation of knowledge and a map for future research. This expert paper provides a comprehensive overview of state-of-the-art technologies to investigate bone properties in humans and animals - including their strengths and weaknesses. New research methodologies are outlined and future strategies are discussed to combine phenotypic with rapidly developing -omics data in order to advance musculoskeletal research and move towards "personalised medicine".

3D assessment of intervertebral disc degeneration in zebrafish identifies changes in bone density that prime disc disease.

Back pain is a common condition with a high social impact and represents a global health burden. Intervertebral disc disease (IVDD) is one of the major causes of back pain; no therapeutics are currently available to reverse this disease. The impact of bone mineral density (BMD) on IVDD has been controversial, with some studies suggesting osteoporosis as causative for IVDD and others suggesting it as protective for IVDD. Functional studies to evaluate the influence of genetic components of BMD in IVDD could highlight opportunities for drug development and repurposing. By taking a holistic 3D approach, we established an aging zebrafish model for spontaneous IVDD. Increased BMD in aging, detected by automated computational analysis, is caused by bone deformities at the endplates. However, aged zebrafish spines showed changes in bone morphology, microstructure, mineral heterogeneity, and increased fragility that resembled osteoporosis. Elements of the discs recapitulated IVDD symptoms found in humans: the intervertebral ligament (equivalent to the annulus fibrosus) showed disorganized collagen fibers and herniation, while the disc center (nucleus pulposus equivalent) showed dehydration and cellular abnormalities. We manipulated BMD in young zebrafish by mutating sp7 and cathepsin K, leading to low and high BMD, respectively. Remarkably, we detected IVDD in both groups, demonstrating that low BMD does not protect against IVDD, and we found a strong correlation between high BMD and IVDD. Deep learning was applied to high-resolution synchrotron µCT image data to analyze osteocyte 3D lacunar distribution and morphology, revealing a role of sp7 in controlling the osteocyte lacunar 3D profile. Our findings suggest potential avenues through which bone quality can be targeted to identify beneficial therapeutics for IVDD.

Welfare Assessment of Adult Laboratory Zebrafish: A Practical Guide.

Teleost fish such as Danio rerio (zebrafish) have been successfully used in biomedical research since decades. Genetically altered fish lines obtained by state-of-the-art genetic technologies are serving as well-known model organisms. In Europe, following Directive 2010/63/EU, generation, breeding, and husbandry of new genetically altered lines of laboratory animals require governmental state approval in case pain, suffering, distress, or long-lasting harm to the offspring derived by breeding of these lines cannot be excluded. The identification and assessment of pain, distress, or harm, according to a severity classification of mild, moderate, severe, or humane endpoint, became a new challenging task for all scientists, animal technicians, and veterinarians for daily work with laboratory zebrafish. In this study, we describe the performance of the assessment of welfare parameters of selected pathologic phenotypes and abnormalities frequently found in laboratory fish facilities based on veterinary, biological, and physiological aspects by using a dedicated score sheet. In a colony of zebrafish, we evaluated the frequency of genotype-independent abnormalities observed within 3 years. We give examples for severity classification and measures once an abnormality has been identified according to the 3Rs (Replacement, Reduction and Refinement).

Wnt16 Elicits a Protective Effect Against Fractures and Supports Bone Repair in Zebrafish.

Bone homeostasis is a dynamic, multicellular process that is required throughout life to maintain bone integrity, prevent fracture, and respond to skeletal damage. WNT16 has been linked to bone fragility and osteoporosis in human genome wide-association studies, as well as the functional hematopoiesis of leukocytes in vivo. However, the mechanisms by which WNT16 promotes bone health and repair are not fully understood. In this study, CRISPR-Cas9 was used to generate mutant zebrafish lacking Wnt16 (wnt16 -/- ) to study its effect on bone dynamically. The wnt16 mutants displayed variable tissue mineral density (TMD) and were susceptible to spontaneous fractures and the accumulation of bone calluses at an early age. Fractures were induced in the lepidotrichia of the caudal fins of wnt16 -/- and WT zebrafish; this model was used to probe the mechanisms by which Wnt16 regulates skeletal and immune cell dynamics in vivo. In WT fins, wnt16 expression increased significantly during the early stages for bone repair. Mineralization of bone during fracture repair was significantly delayed in wnt16 mutants compared with WT zebrafish. Surprisingly, there was no evidence that the recruitment of innate immune cells to fractures or soft callus formation was altered in wnt16 mutants. However, osteoblast recruitment was significantly delayed in wnt16 mutants postfracture, coinciding with precocious activation of the canonical Wnt signaling pathway. In situ hybridization suggests that canonical Wnt-responsive cells within fractures are osteoblast progenitors, and that osteoblast differentiation during bone repair is coordinated by the dynamic expression of runx2a and wnt16. This study highlights zebrafish as an emerging model for functionally validating osteoporosis-associated genes and investigating fracture repair dynamically in vivo. Using this model, it was found that Wnt16 protects against fracture and supports bone repair, likely by modulating canonical Wnt activity via runx2a to facilitate osteoblast differentiation and bone matrix deposition. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC. on behalf of American Society for Bone and Mineral Research.

Transformed notochordal cells trigger chronic wounds in zebrafish, destabilizing the vertebral column and bone homeostasis.

Notochordal cells play a pivotal role in vertebral column patterning, contributing to the formation of the inner architecture of intervertebral discs (IVDs). Their disappearance during development has been associated with reduced repair capacity and IVD degeneration. Notochord cells can give rise to chordomas, a highly invasive bone cancer associated with late diagnosis. Understanding the impact of neoplastic cells during development and on the surrounding vertebral column could open avenues for earlier intervention and therapeutics. We investigated the impact of transformed notochord cells in the zebrafish skeleton using a line expressing RAS in the notochord under the control of the kita promoter, with the advantage of adulthood endurance. Transformed cells caused damage in the notochord and destabilised the sheath layer, triggering a wound repair mechanism, with enrolment of sheath cells (col9a2+) and expression of wt1b, similar to induced notochord wounds. Moreover, increased recruitment of neutrophils and macrophages, displaying abnormal behaviour in proximity to the notochord sheath and transformed cells, supported parallels between chordomas, wound and inflammation. Cancerous notochordal cells interfere with differentiation of sheath cells to form chordacentra domains, leading to fusions and vertebral clefts during development. Adults displayed IVD irregularities reminiscent of degeneration, including reduced bone mineral density and increased osteoclast activity, along with disorganised osteoblasts and collagen, indicating impaired bone homeostasis. By depleting inflammatory cells, we abrogated chordoma development and rescued the skeletal features of the vertebral column. Therefore, we showed that transformed notochord cells alter the skeleton during life, causing a wound-like phenotype and activating chronic wound response, suggesting parallels between chordoma, wound, IVD degeneration and inflammation, highlighting inflammation as a promising target for future therapeutics. This article has an associated First Person interview with the first author of the paper.

Genome-wide meta-analysis of muscle weakness identifies 15 susceptibility loci in older men and women.

Low muscle strength is an important heritable indicator of poor health linked to morbidity and mortality in older people. In a genome-wide association study meta-analysis of 256,523 Europeans aged 60 years and over from 22 cohorts we identify 15 loci associated with muscle weakness (European Working Group on Sarcopenia in Older People definition: n = 48,596 cases, 18.9% of total), including 12 loci not implicated in previous analyses of continuous measures of grip strength. Loci include genes reportedly involved in autoimmune disease (HLA-DQA1 p = 4 × 10-17), arthritis (GDF5 p = 4 × 10-13), cell cycle control and cancer protection, regulation of transcription, and others involved in the development and maintenance of the musculoskeletal system. Using Mendelian randomization we report possible overlapping causal pathways, including diabetes susceptibility, haematological parameters, and the immune system. We conclude that muscle weakness in older adults has distinct mechanisms from continuous strength, including several pathways considered to be hallmarks of ageing.

NCOA3 identified as a new candidate to explain autosomal dominant progressive hearing loss.

Hearing loss is a frequent sensory impairment in humans and genetic factors account for an elevated fraction of the cases. We have investigated a large family of five generations, with 15 reported individuals presenting non-syndromic, sensorineural, bilateral and progressive hearing loss, segregating as an autosomal dominant condition. Linkage analysis, using SNP-array and selected microsatellites, identified a region of near 13 cM in chromosome 20 as the best candidate to harbour the causative mutation. After exome sequencing and filtering of variants, only one predicted deleterious variant in the NCOA3 gene (NM_181659, c.2810C > G; p.Ser937Cys) fit in with our linkage data. RT-PCR, immunostaining and in situ hybridization showed expression of ncoa3 in the inner ear of mice and zebrafish. We generated a stable homozygous zebrafish mutant line using the CRISPR/Cas9 system. ncoa3-/- did not display any major morphological abnormalities in the ear, however, anterior macular hair cells showed altered orientation. Surprisingly, chondrocytes forming the ear cartilage showed abnormal behaviour in ncoa3-/-, detaching from their location, invading the ear canal and blocking the cristae. Adult mutants displayed accumulation of denser material wrapping the otoliths of ncoa3-/- and increased bone mineral density. Altered zebrafish swimming behaviour corroborates a potential role of ncoa3 in hearing loss. In conclusion, we identified a potential candidate gene to explain hereditary hearing loss, and our functional analyses suggest subtle and abnormal skeletal behaviour as mechanisms involved in the pathogenesis of progressive sensory function impairment.

Opportunities and Challenges in Functional Genomics Research in Osteoporosis: Report From a Workshop Held by the Causes Working Group of the Osteoporosis and Bone Research Academy of the Royal Osteoporosis Society on October 5th 2020.

The discovery that sclerostin is the defective protein underlying the rare heritable bone mass disorder, sclerosteosis, ultimately led to development of anti-sclerostin antibodies as a new treatment for osteoporosis. In the era of large scale GWAS, many additional genetic signals associated with bone mass and related traits have since been reported. However, how best to interrogate these signals in order to identify the underlying gene responsible for these genetic associations, a prerequisite for identifying drug targets for further treatments, remains a challenge. The resources available for supporting functional genomics research continues to expand, exemplified by "multi-omics" database resources, with improved availability of datasets derived from bone tissues. These databases provide information about potential molecular mediators such as mRNA expression, protein expression, and DNA methylation levels, which can be interrogated to map genetic signals to specific genes based on identification of causal pathways between the genetic signal and the phenotype being studied. Functional evaluation of potential causative genes has been facilitated by characterization of the "osteocyte signature", by broad phenotyping of knockout mice with deletions of over 7,000 genes, in which more detailed skeletal phenotyping is currently being undertaken, and by development of zebrafish as a highly efficient additional in vivo model for functional studies of the skeleton. Looking to the future, this expanding repertoire of tools offers the hope of accurately defining the major genetic signals which contribute to osteoporosis. This may in turn lead to the identification of additional therapeutic targets, and ultimately new treatments for osteoporosis.

Finite element and deformation analyses predict pattern of bone failure in loaded zebrafish spines.

The spine is the central skeletal support structure in vertebrates consisting of repeated units of bone, the vertebrae, separated by intervertebral discs (IVDs) that enable the movement of the spine. Spinal pathologies such as idiopathic back pain, vertebral compression fractures and IVD failure affect millions of people worldwide. Animal models can help us to understand the disease process, and zebrafish are increasingly used as they are highly genetically tractable, their spines are axially loaded like humans, and they show similar pathologies to humans during ageing. However, biomechanical models for the zebrafish are largely lacking. Here, we describe the results of loading intact zebrafish spinal motion segments on a material testing stage within a micro-computed tomography machine. We show that vertebrae and their arches show predictable patterns of deformation prior to their ultimate failure, in a pattern dependent on their position within the segment. We further show using geometric morphometrics which regions of the vertebra deform the most during loading, and that finite-element models of the trunk subjected reflect the real patterns of deformation and strain seen during loading and can therefore be used as a predictive model for biomechanical performance.

Scleraxis genes are required for normal musculoskeletal development and for rib growth and mineralization in zebrafish.

Tendons are an essential part of the musculoskeletal system, connecting muscle and skeletal elements to enable force generation. The transcription factor scleraxis marks vertebrate tendons from early specification. Scleraxis-null mice are viable and have a range of tendon and bone defects in the trunk and limbs but no described cranial phenotype. We report the expression of zebrafish scleraxis orthologs: scleraxis homolog (scx)-a and scxb in cranial and intramuscular tendons and in other skeletal elements. Single mutants for either scxa or scxb, generated by clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9), are viable and fertile as adult fish. Although scxb mutants show no obvious phenotype, scxa mutant embryos have defects in cranial tendon maturation and muscle misalignment. Mutation of both scleraxis genes results in more severe defects in cranial tendon differentiation, muscle and cartilage dysmorphogenesis and paralysis, and lethality by 2-5 wk, which indicates an essential function of scleraxis for craniofacial development. At juvenile and adult stages, ribs in scxa mutants fail to mineralize and/or are small and heavily fractured. Scxa mutants also have smaller muscle volume, abnormal swim movement, and defects in bone growth and composition. Scleraxis function is therefore essential for normal craniofacial form and function and vital for fish development.-Kague, E., Hughes, S. M., Lawrence, E. A., Cross, S., Martin-Silverstone, E., Hammond, C. L., Hinits, Y. Scleraxis genes are required for normal musculoskeletal development and for rib growth and mineralization in zebrafish.

Craniofrontonasal Syndrome Caused by Introduction of a Novel uATG in the 5'UTR of EFNB1.

Craniofrontonasal syndrome (CFNS) is an X-linked disorder caused by EFNB1 mutations in which females are more severely affected than males. Severe male phenotypes are associated with mosaicism, supporting cellular interference for sex bias in this disease. Although many variants have been found in the coding region of EFNB1, only 2 pathogenic variants have been identified in the same nucleotide in 5'UTR, disrupting the stop codon of an upstream open reading frame (uORF). uORFs are known to be part of a wide range of post-transcriptional regulation processes, and just recently, their association with human diseases has come to light. In the present study, we analyzed EFNB1 in a female patient with typical features of CFNS. We identified a variant, located at c.-411, creating a new upstream ATG (uATG) in the 5'UTR of EFNB1, which is predicted to alter an existing uORF. Dual-luciferase reporter assays showed significant reduction in protein translation, but no difference in the mRNA levels. Our study demonstrates, for the first time, the regulatory impact of uATG formation on EFNB1 levels and suggests that this should be the target region in molecular diagnosis of CFNS cases without pathogenic variants in the coding and splice sites regions of EFNB1.