Bruck syndrome (BS) is a disorder characterized by joint flexion contractures and skeletal dysplasia that shows strong clinical overlap with the brittle bone disease Osteogenesis Imperfecta (OI). BS is caused by bi-allelic mutations in either the FKBP10 or the PLOD2 gene. PLOD2 encodes the lysyl hydroxylase 2 (LH2) enzyme, which is responsible for the hydroxylation of lysine residues in fibrillar collagen telopeptides. This hydroxylation directs cross-linking of collagen fibrils in the extracellular matrix, which is necessary to provide stability and tensile integrity to the collagen fibrils. To further elucidate the function of LH2 in vertebrate skeletal development, we created a zebrafish model harboring a homozygous plod2 nonsense mutation resulting in reduced telopeptide hydroxylation and cross-linking of bone type I collagen. Adult plod2 mutants present with a shortened body axis and severe skeletal abnormalities with evidence of bone fragility and fractures. The vertebral column of plod2 mutants is short and scoliotic with compressed vertebrae that show excessive bone formation at the vertebral end plates, and increased tissue mineral density in the vertebral centra. The muscle fibers of mutant zebrafish have a reduced diameter near the horizontal myoseptum. The endomysium, a layer of connective tissue ensheathing the individual muscle fibers, is enlarged. Transmission electron microscopy of mutant vertebral bone shows type I collagen fibrils that are less organized with loss of the typical plywood-like structure. In conclusion, plod2 mutant zebrafish show molecular and tissue abnormalities in the musculoskeletal system that are concordant with clinical findings in BS patients. Therefore, the plod2 zebrafish mutant is a promising model for the elucidation of the underlying pathogenetic mechanisms leading to BS and the development of novel therapeutic avenues in this syndrome.

Phenomics, which ideally involves in-depth phenotyping at the whole-organism scale, may enhance our functional understanding of genetic variation. Here, we demonstrate methods to profile hundreds of phenotypic measures comprised of morphological and densitometric traits at a large number of sites within the axial skeleton of adult zebrafish. We show the potential for vertebral patterns to confer heightened sensitivity, with similar specificity, in discriminating mutant populations compared to analyzing individual vertebrae in isolation. We identify phenotypes associated with human brittle bone disease and thyroid stimulating hormone receptor hyperactivity. Finally, we develop allometric models and show their potential to aid in the discrimination of mutant phenotypes masked by alterations in growth. Our studies demonstrate virtues of deep phenotyping in a spatially distributed organ system. Analyzing phenotypic patterns may increase productivity in genetic screens, and facilitate the study of genetic variants associated with smaller effect sizes, such as those that underlie complex diseases.

Over the last years the zebrafish imposed itself as a powerful model to study skeletal diseases, but a limit to its use is the poor characterization of collagen type I, the most abundant protein in bone and skin. In tetrapods collagen type I is a trimer mainly composed of two α1 chains and one α2 chain, encoded by COL1A1 and COL1A2 genes, respectively. In contrast, in zebrafish three type I collagen genes exist, col1a1a, col1a1b and col1a2 coding for α1(I), α3(I) and α2(I) chains. During embryonic and larval development the three collagen type I genes showed a similar spatio-temporal expression pattern, indicating their co-regulation and interdependence at these stages. In both embryonic and adult tissues, the presence of the three α(I) chains was demonstrated, although in embryos α1(I) was present in two distinct glycosylated states, suggesting a developmental-specific collagen composition. Even though in adult bone, skin and scales equal amounts of α1(I), α3(I) and α2(I) chains are present, the presented data suggest a tissue-specific stoichiometry and/or post-translational modification status for collagen type I. In conclusion, this data will be useful to properly interpret results and insights gained from zebrafish models of skeletal diseases.

SignificanceType I collagenopathies are a heterogenous group of connective tissue disorders, caused by genetic defects in type I collagen. Inherent to these disorders is a large clinical variability, of which the underlying molecular basis remains undefined. By systematically analyzing skeletal phenotypes in a large set of type I collagen zebrafish mutants, we show that zebrafish models are able to both genocopy and phenocopy different forms of human type I collagenopathies, arguing for a similar pathogenetic basis. This study illustrates the future potential of zebrafish as a tool to further dissect the molecular basis of phenotypic variability in human type I collagenopathies, to improve diagnostic strategies as well as promote the discovery of new targetable pathways for pharmacological intervention of these disorders.