Loss of Type I Collagen Telopeptide Lysyl Hydroxylation Causes Musculoskeletal Abnormalities in a Zebrafish Model of Bruck Syndrome

Disruption of the moonlighting function of CTF18 in a patient with T-lymphopenia

Introduction

Newborn screening for immunodeficiency has led to the identification of numerous cases for which the causal etiology is unknown.

Methods

Here we report the diagnosis of T lymphopenia of unknown etiology in a male proband. Whole exome sequencing (WES) was employed to nominate candidate variants, which were then analyzed functionally in zebrafish and in mice bearing orthologous mutations.

Results

WES revealed missense mutations in CHTF18 that were inherited in an autosomal recessive manner. CTF18, encoded by the CHTF18 gene, is a component of a secondary clamp loader, which is primarily thought to function by promoting DNA replication. We determined that the patient's variants in CHTF18 (CTF18 R751W and E851Q) were damaging to function and severely attenuated the capacity of CTF18 to support hematopoiesis and lymphoid development, strongly suggesting that they were responsible for his T lymphopenia; however, the function of CTF18 appeared to be unrelated to its role as a clamp loader. DNA-damage, expected when replication is impaired, was not evident by expression profiling in murine Chtf18 mutant hematopoietic stem and progenitor cells (HSPC), nor was development of Ctf18-deficient progenitors rescued by p53 loss. Instead, we observed an expression signature suggesting disruption of HSPC positioning and migration. Indeed, the positioning of HSPC in ctf18 morphant zebrafish embryos was perturbed, suggesting that HSPC function was impaired through disrupted positioning in hematopoietic organs.

Discussion

Accordingly, we propose that T lymphopenia in our patient resulted from disturbed cell-cell contacts and migration of HSPC, caused by a non-canonical function of CHTF18 in regulating gene expression.

MicroRNAs as Biomarker in Rheumatoid Arthritis: Pathogenesis to Clinical Relevance

MicroRNAs (miRNAs) have emerged as intricate players in rheumatoid arthritis (RA), holding promise as discerning biomarkers for diagnostic and prognostic purposes. The lack of sensitivity and specificity in current diagnostic techniques, such as rheumatoid factor (RF) and anti-citrullinated protein antibodies (ACPA), causes diagnosis delays in RA. The miR-146a and miR-155 act in inflammatory cascades and reduce joint deterioration, and miR-223 is paradoxical, acting differently in different illness scenarios. The microenvironment of RA is shaped by the complex modulation of gene expression and cytokine dynamics by miR-126 and miR-24. miRNAs serve as a promising candidate for precision medicine in the management of RA. There are obstacles encountered in validation, delivery optimization, and off-target effect mitigation before miRNA-based biomarkers may be applied in clinical settings. Machine learning (ML) and artificial intelligence (AI) have been used to integrate miRNA expression patterns with clinical data to greatly advance the treatment of RA. Because of the disease's inherent complexity and variability, these state-of-the-art models provide accurate predictions regarding the onset, development, and response to treatment of RA. By using clinical information and miRNA expression data, ML algorithms are revolutionizing the treatment of RA by predicting the onset and course of the disease with remarkably high accuracy. The development of therapeutic modalities and miRNA profiling has great potential to transform the diagnosis, prognosis, and treatment of RA, providing fresh hope for better patient outcomes.

Bmpr1aa modulates the severity of the skeletal phenotype in an fkbp10-deficient Bruck syndrome zebrafish model

Rare monogenic disorders often exhibit significant phenotypic variability among individuals sharing identical genetic mutations. Bruck syndrome (BS), a prime example, is characterized by bone fragility and congenital contractures, although with a pronounced variability among family members. BS arises from recessive biallelic mutations in FKBP10 or PLOD2. FKBP65, the protein encoded by FKBP10, collaborates with the LH2 enzyme (PLOD2) in type I collagen telopeptide lysine hydroxylation, crucial for collagen cross-linking. To identify potential modifier genes and to investigate the mechanistic role of FKBP10 in BS pathogenesis, we established an fkbp10a knockout zebrafish model. Mass-spectrometry analysis in fkbp10a-/- mutants revealed a generally decreased type I collagen lysyl hydroxylation, paralleled by a wide skeletal variability similar to human patients. Ultrastructural examination of the skeleton in severely affected mutants showed enlarged type I collagen fibrils and disturbed elastin layers. Whole-exome sequencing of 7 mildly and 7 severely affected mutant zebrafish siblings, followed by single nucleotide polymorphism-based linkage analysis, indicated a linked region on chromosome 13, which segregates with phenotypic severity. Transcriptome analysis identified 6 differentially expressed genes (DEGs) between mildly and severely affected mutants. The convergence of genes within the linked region and DEGs highlighted bmpr1aa as a potential modifier gene, as its reduced expression correlates with increased skeletal severity. In summary, our study provides deeper insights into the role of FKBP10 in BS pathogenesis. Additionally, we identified a pivotal gene that influences phenotypic severity in a zebrafish model of BS. These findings hold promise for novel treatments in the field of bone diseases.

Cellular and Molecular Effects of the Bruck Syndrome-Associated Mutation in the PLOD2 Gene

Bruck syndrome is a rare autosomal recessive disorder characterized by increased bone fragility and joint contractures similar to those in arthrogryposis and is known to be associated with mutations in the FKBP10 (FKBP prolyl isomerase 10) and PLOD2 (Procollagen-Lysine,2-Oxoglutarate 5-Dioxygenase 2) genes. These genes encode endoplasmic reticulum proteins that play an important role in the biosynthesis of type I collagen, which in turn affects the structure and strength of connective tissues and bones in the body. Mutations are associated with disturbances in both the primary collagen chain and its post-translational formation, but the mechanism by which mutations lead to Bruck syndrome phenotypes has not been determined, not only because of the small number of patients who come to the attention of researchers but also because of the lack of disease models. In our work, we investigated the cellular effects of two forms of the wild-type PLOD2 gene, as well as the PLOD2 gene with homozygous mutation c.1885A>G (p.Thr629Ala). The synthesized genetic constructs were transfected into HEK293 cell line and human skin fibroblasts (DF2 line). The localization of PLOD2 protein in cells and the effects caused by the expression of different isoforms-long, short, and long with mutation-were analyzed. In addition, the results of the transcriptome analysis of a patient with Bruck syndrome, in whom this mutation was detected, are presented.

Zebrafish Models for Skeletal and Extraskeletal Osteogenesis Imperfecta Features: Unveiling Pathophysiology and Paving the Way for Drug Discovery

In the last decades, the easy genetic manipulation, the external fertilization, the high percentage of homology with human genes and the reduced husbandry costs compared to rodents, made zebrafish a valid model for studying human diseases and for developing new therapeutical strategies. Since zebrafish shares with mammals the same bone cells and ossification types, it became widely used to dissect mechanisms and possible new therapeutic approaches in the field of common and rare bone diseases, such as osteoporosis and osteogenesis imperfecta (OI), respectively. OI is a heritable skeletal disorder caused by defects in gene encoding collagen I or proteins/enzymes necessary for collagen I synthesis and secretion. Nevertheless, OI patients can be also characterized by extraskeletal manifestations such as dentinogenesis imperfecta, muscle weakness, cardiac valve and pulmonary abnormalities and skin laxity. In this review, we provide an overview of the available zebrafish models for both dominant and recessive forms of OI. An updated description of all the main similarities and differences between zebrafish and mammal skeleton, muscle, heart and skin, will be also discussed. Finally, a list of high- and low-throughput techniques available to exploit both larvae and adult OI zebrafish models as unique tools for the discovery of new therapeutic approaches will be presented.

Generation of bone-specific lysyl hydroxylase 2 knockout mice and their phenotypes

Lysyl hydroxylase 2 (LH2) catalyzes the hydroxylation of lysine residues in the telopeptides of type I collagen. This modification is critical for the formation of stable hydroxylysine-aldehyde derived collagen cross-links, thus, for the stability of collagen fibrils. Though dysfunction of LH2 causes Bruck syndrome, recessive osteogenesis imperfecta with joint contracture, the molecular mechanisms by which LH2 affects bone formation are still not well understood. Since the Plod2 knockout mice are embryonically lethal, we generated bone-specific LH2 conditional knockout mice (bsLH2-cKO) using the osteocalcin-Cre/loxP system, and evaluated phenotypes of femurs. LH2 mRNA and protein levels assessed by qPCR, immunohistochemistry and Data Independent Acquisition proteomics were all markedly low in bsLH2-cKO femurs when compared to controls. Lysine hydroxylation of both carboxy- and amino-terminal telopeptides of an α1(I) chain were significantly diminished resulting in reduction of the hydroxylysine-aldehyde derived cross-links. The collagen fibrils in bsLH2-cKO appeared to be thicker, often fused and irregular when compared to controls. In addition, bone mineral density and mechanical properties of bsLH2-cKO femurs were significantly impaired. Taken together, these data demonstrate that LH2-catalyzed modification and consequent cross-linking of collagen are critical for proper bone formation and mechanical strength.

ColPTMScape: An open access knowledge base for tissue-specific collagen PTM maps

Collagen is a key component of the extracellular matrix (ECM). In the remodeling of ECM, a remarkable variation in collagen post-translational modifications (PTMs) occurs. This makes collagen a potential target for understanding extracellular matrix remodeling during pathological conditions. Over the years, scientists have gathered a huge amount of data about collagen PTM during extracellular matrix remodeling. To make such information easily accessible in a consolidated space, we have developed ColPTMScape (https://colptmscape.iitmandi.ac.in/), a dedicated knowledge base for collagen PTMs. The identified site-specific PTMs, quantitated PTM sites, and PTM maps of collagen chains are deliverables to the scientific community, especially to matrix biologists. Through this knowledge base, users can easily gain information related to the difference in the collagen PTMs across different tissues in different organisms.

Zebrafish as a Model for Osteoporosis: Functional Validations of Genome-Wide Association Studies

PURPOSE OF REVIEW

GWAS, as a largely correlational analysis, requires in vitro or in vivo validation. Zebrafish (Danio rerio) have many advantages for studying the genetics of human diseases. Since gene editing in zebrafish has been highly valuable for studying embryonic skeletal developmental processes that are prenatally or perinatally lethal in mammalian models, we are reviewing pros and cons of this model.

RECENT FINDINGS

The true power for the use of zebrafish is the ease by which the genome can be edited, especially using the CRISPR/Cas9 system. Gene editing, followed by phenotyping, for complex traits such as BMD, is beneficial, but the major physiological differences between the fish and mammals must be considered. Like mammals, zebrafish do have main bone cells; thus, both in vivo stem cell analyses and in vivo imaging are doable. Yet, the "long" bones of fish are peculiar, and their bone cavities do not contain bone marrow. Partial duplication of the zebrafish genome should be taken into account. Overall, small fish toolkit can provide unmatched opportunities for genetic modifications and morphological investigation as a follow-up to human-first discovery.

Significance of Premature Vertebral Mineralization in Zebrafish Models in Mechanistic and Pharmaceutical Research on Hereditary Multisystem Diseases

Zebrafish are increasingly becoming an important model organism for studying the pathophysiological mechanisms of human diseases and investigating how these mechanisms can be effectively targeted using compounds that may open avenues to novel treatments for patients. The zebrafish skeleton has been particularly instrumental in modeling bone diseases as-contrary to other model organisms-the lower load on the skeleton of an aquatic animal enables mutants to survive to early adulthood. In this respect, the axial skeletons of zebrafish have been a good read-out for congenital spinal deformities such as scoliosis and degenerative disorders such as osteoporosis and osteoarthritis, in which aberrant mineralization in humans is reflected in the respective zebrafish models. Interestingly, there have been several reports of hereditary multisystemic diseases that do not affect the vertebral column in human patients, while the corresponding zebrafish models systematically show anomalies in mineralization and morphology of the spine as their leading or, in some cases, only phenotype. In this review, we describe such examples, highlighting the underlying mechanisms, the already-used or potential power of these models to help us understand and amend the mineralization process, and the outstanding questions on how and why this specific axial type of aberrant mineralization occurs in these disease models.

Features of Congenital Arthrogryposis Due to Abnormalities in Collagen Homeostasis, a Scoping Review

Congenital arthrogryposis (CA) refers to the presence of multiple contractures at birth. It is a feature of several inherited syndromes, notable amongst them are disorders of collagen formation. This review aims to characterize disorders that directly or indirectly impact collagen structure and function leading to CA in search for common phenotypic or pathophysiological features, possible genotype-phenotype correlation, and potential novel treatment approaches based on a better understanding of the underlying pathomechanism. Nine genes, corresponding to five clinical phenotypes, were identified after a literature search. The most notable trend was the extreme phenotype variability. Clinical features across all syndromes ranged from subtle with minimal congenital contractures, to severe with multiple congenital contractures and extra-articular features including skin, respiratory, or other manifestations. Five of the identified genes were involved in the function of the Lysyl Hydroxylase 2 or 3 enzymes, which enable the hydroxylation and/or glycosylation of lysyl residues to allow the formation of the collagen superstructure. Whilst current treatment approaches are post-natal surgical correction, there are also potential in-utero therapies being developed. Cyclosporin A showed promise in treating collagen VI disorders although there is an associated risk of immunosuppression. The treatments that could be in the clinical trials soon are the splice correction therapies in collagen VI-related disorders.

Why Animal Experiments Are Still Indispensable in Bone Research: A Statement by the European Calcified Tissue Society

Major achievements in bone research have always relied on animal models and in vitro systems derived from patient and animal material. However, the use of animals in research has drawn intense ethical debate and the complete abolition of animal experimentation is demanded by fractions of the population. This phenomenon is enhanced by the reproducibility crisis in science and the advance of in vitro and in silico techniques. 3D culture, organ-on-a-chip, and computer models have improved enormously over the last few years. Nevertheless, the overall complexity of bone tissue cross-talk and the systemic and local regulation of bone physiology can often only be addressed in entire vertebrates. Powerful genetic methods such as conditional mutagenesis, lineage tracing, and modeling of the diseases enhanced the understanding of the entire skeletal system. In this review endorsed by the European Calcified Tissue Society (ECTS), a working group of investigators from Europe and the US provides an overview of the strengths and limitations of experimental animal models, including rodents, fish, and large animals, as well the potential and shortcomings of in vitro and in silico technologies in skeletal research. We propose that the proper combination of the right animal model for a specific hypothesis and state-of-the-art in vitro and/or in silico technology is essential to solving remaining important questions in bone research. This is crucial for executing most efficiently the 3R principles to reduce, refine, and replace animal experimentation, for enhancing our knowledge of skeletal biology, and for the treatment of bone diseases that affect a large part of society. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Examining craniofacial variation among crispant and mutant zebrafish models of human skeletal diseases

Genetic diseases affecting the skeletal system present with a wide range of symptoms that make diagnosis and treatment difficult. Genome-wide association and sequencing studies have identified genes linked to human skeletal diseases. Gene editing of zebrafish models allows researchers to further examine the link between genotype and phenotype, with the long-term goal of improving diagnosis and treatment. While current automated tools enable rapid and in-depth phenotyping of the axial skeleton, characterizing the effects of mutations on the craniofacial skeleton has been more challenging. The objective of this study was to evaluate a semi-automated screening tool can be used to quantify craniofacial variations in zebrafish models using four genes that have been associated with human skeletal diseases (meox1, plod2, sost, and wnt16) as test cases. We used traditional landmarks to ground truth our dataset and pseudolandmarks to quantify variation across the 3D cranial skeleton between the groups (somatic crispant, germline mutant, and control fish). The proposed pipeline identified variation between the crispant or mutant fish and control fish for four genes. Variation in phenotypes parallel human craniofacial symptoms for two of the four genes tested. This study demonstrates the potential as well as the limitations of our pipeline as a screening tool to examine multi-dimensional phenotypes associated with the zebrafish craniofacial skeleton.

Expanding the phenotype of Bruck syndrome: Severe limb deformity, arthrogryposis, congenital cardiac disease and pulmonary hemorrhage

Bruck syndrome is a rare collagen disorder with autosomal recessive inheritance caused by pathogenic variants in either FKBP10 or PLOD2 genes. It is characterized by bone fragility and fractures similar in severity and variability to osteogenesis-imperfecta as well as congenital joint contractures. This article describes an infant with a homozygous (partial) gene deletion of PLOD2 that includes the start codon and would be expected to lead to nonfunctional protein product. The infant had a severe phenotype of Bruck syndrome and is the only reported case of Bruck syndrome with congenital cardiac disease (triscuspid valve dysplasia with severe regurgitation, mitral valve prolapses with moderate regurgitation, and pulmonary hypertension) and pulmonary hemorrhage. We hypothesize that the additional feature of congenital cardiac disease in this case was due to the underlying defect in type I collagen, and that the pulmonary hemorrhage was multifactorial, with underlying vessel fragility, rib fractures, and high pulmonary pressures likely to be major contributing factors. Management was largely supportive with the use of bisphosphonates to assist in pain management. Care was complicated by comorbid cardiopulmonary compromise, limited evidence-base guiding care, and difficulties in discussing end-of-life care.

An Update on Animal Models of Osteogenesis Imperfecta

Osteogenesis imperfecta (OI) is a heterogeneous disorder characterized by bone fragility, multiple fractures, bone deformity, and short stature. In recent years, the application of next generation sequencing has triggered the discovery of many new genetic causes for OI. Until now, more than 25 genetic causes of OI and closely related disorders have been identified. However, the mechanisms of many genes on skeletal fragility in OI are not entirely clear. Animal models of OI could help to understand the cellular, signaling, and metabolic mechanisms contributing to the disease, and how targeting these pathways can provide therapeutic targets. To date, a lot of animal models, mainly mice and zebrafish, have been described with defects in 19 OI-associated genes. In this review, we summarize the known genetic causes and animal models that recapitulate OI with a main focus on engineered mouse and zebrafish models. Additionally, we briefly discuss domestic animals with naturally occurring OI phenotypes. Knowledge of the specific molecular basis of OI will advance clinical diagnosis and potentially stimulate targeted therapeutic approaches.

Lysyl hydroxylase 2 mediated collagen post-translational modifications and functional outcomes

Lysyl hydroxylase 2 (LH2) is a member of LH family that catalyzes the hydroxylation of lysine (Lys) residues on collagen, and this particular isozyme has been implicated in various diseases. While its function as a telopeptidyl LH is generally accepted, several fundamental questions remain unanswered: 1. Does LH2 catalyze the hydroxylation of all telopeptidyl Lys residues of collagen? 2. Is LH2 involved in the helical Lys hydroxylation? 3. What are the functional consequences when LH2 is completely absent? To answer these questions, we generated LH2-null MC3T3 cells (LH2KO), and extensively characterized the type I collagen phenotypes in comparison with controls. Cross-link analysis demonstrated that the hydroxylysine-aldehyde (Hyl^ald)-derived cross-links were completely absent from LH2KO collagen with concomitant increases in the Lys^ald-derived cross-links. Mass spectrometric analysis revealed that, in LH2KO type I collagen, telopeptidyl Lys hydroxylation was completely abolished at all sites while helical Lys hydroxylation was slightly diminished in a site-specific manner. Moreover, di-glycosylated Hyl was diminished at the expense of mono-glycosylated Hyl. LH2KO collagen was highly soluble and digestible, fibril diameters were diminished, and mineralization impaired when compared to controls. Together, these data underscore the critical role of LH2-catalyzed collagen modifications in collagen stability, organization and mineralization in MC3T3 cells.

Bruck syndrome in 13 new patients: Identification of five novel FKBP10 and PLOD2 variants and further expansion of the phenotypic spectrum

Bruck Syndrome (BS) is a very rare disorder characterized by osteogenesis imperfecta (OI) associated with congenital contractures and is caused by mutations in FKBP10 or PLOD2 genes. Herein, we describe 13 patients from 9 unrelated Egyptian families with BS. All patients had white sclerae, recurrent fractures, kyphoscoliosis and osteoporosis with variable degrees of severity. Large joint contractures were seen in 11 patients, one patient had contractures of small interphalangeal joints, and one patient had no contractures. Unusual findings noted in individual patients included microcephaly, dental malocclusion, enamel hypoplasia, unilateral congenital dislocation of knee joint, prominent tailbone, and myopathy. Nine different variants were identified in FKBP10 and PLOD2 including five novel ones. FKBP10 variants were found in six families (67%) while PLOD2 variants were identified in three families (33%). The four families, with two affected sibs each, showed inter- and intrafamilial phenotypic variability. In conclusion, we report five novel variants in FKBP10 and PLOD2 thus, expanding the mutational spectrum of BS. In addition, our results expand the phenotypic spectrum, describe newly associated orodental findings, and further illustrate the phenotypic overlap between OI and Bruck syndrome supporting the suggestion of considering BS as a variant of OI rather than a separate entity.

Transcriptome analysis of molecular response to UVC irradiation in zebrafish embryos

Ultraviolet (UV) rays can be both harmful and beneficial to humans. This study aimed to investigate the toxicity and safety of ultraviolet C (UVC) exposure in living organisms and the corresponding biodefense molecular mechanisms. Zebrafish embryos, at an early developmental stage (5-6 h post-fertilization), were irradiated with increasing UVC dosages using high-efficiency deep-ultraviolet light-emitting diodes (278 nm). Morphological phenotypes including survival rate, hatching rate, heart rate, and malformation rate were evaluated. Compared to un-irradiated controls, all zebrafish embryos exposed to 4.5 mJ/cm² UVC survived and showed no significant difference in hatching and heart rate. However, 7.5 mJ/cm² of UVC irradiation caused a significantly decreased survival rate (37.5%) and an increased malformation rate (81.8%). Therefore, 4.5 mJ/cm² was chosen as the limit dosage that the internal biodefense system of zebrafish embryos can protect against UVC radiation. Transcriptome analysis (RNA sequencing) performed on 3 min and 3 days post-irradiation embryos (4.5 mJ/cm²) revealed the molecular mechanisms underlying the response of zebrafish embryos to irradiation. The embryos quickly responded to UVC-induced stress by activating the p53 signaling pathway. In addition, after 3 days of recuperation, the embryos showed activation of signal transducer and activator of transcription (STAT) signaling pathway. To our knowledge, this is the first study to evaluate the toxicological effects and the molecular mechanism of biodefense in zebrafish embryos upon 278 nm UVC irradiation.

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.

Genetic Analysis and Functional Study of a Pedigree With Bruck Syndrome Caused by PLOD2 Variant

BACKGROUND

Bruck syndrome (BS) is a rare autosomal recessive inherited osteogenesis imperfecta disease characterized by increased bone fragility and joint contracture. The pathogenic gene of type I BS is FKBPl0, whereas that of type II BS is PLOD2. No significant difference has been found in the clinical phenotype between the two types of BS. In this study, we performed genetic analysis of a BS pedigree caused by PLOD2 variant and studied the corresponding cellular function.

METHODS

Serum biochemistry, parathyroid hormone (PTH), 25-hydroxyvitamin D [25-(OH) D], osteocalcin, and 24-h urinary calcium levels of a family member with BS was assessed. The genes of the proband were analyzed by second-generation sequencing and exon capture techniques. Sanger sequencing was also performed for the suspected responsible variant of the family member. Wild- and variant-type lentivirus plasmids were constructed by gene cloning and transfected into HEK293T cells. Cell function was verified by real-time quantitative polymerase chain reaction, western blotting, and immunofluorescence detection.

RESULTS

In this pedigree, the proband was found to have a homozygous variant c.1856G > A (p.Arg619His) in exon 17 of PLOD2 (NM_182943.3). His consanguineous parents and sisters were p.Arg619His heterozygous carriers. The mRNA expression of PLOD2 in the constructed p.Arg619His variant cells was significantly upregulated, while the expression of PLOD2 and collagen I protein in the cell lysate was significantly downregulated. Immunofluorescence revealed that the wild-type PLOD2 was mainly located in the cytoplasm, and the expression of the PLOD2 protein after c.1856G > A variant was significantly downregulated, with almost no expression, aligning with the western blot results. The serum sodium, potassium, calcium, phosphorus, magnesium, alkaline phosphatase, PTH, 25-(OH) D, osteocalcin, and 24 h urinary calcium levels of the proband, his parents, and sisters were normal.

CONCLUSION

Through gene and cell function analyses, PLOD2 Arg619His missense variant was preliminarily confirmed to cause BS by reducing protein expression.

Decrease of lysyl hydroxylase 2 activity causes abnormal collagen molecular phenotypes, defective mineralization and compromised mechanical properties of bone

Lysyl hydroxylase 2 (LH2) is an enzyme that catalyzes the hydroxylation of lysine (Lys) residues in fibrillar collagen telopeptides, a critical post-translational modification for the stability of intermolecular cross-links. Though abnormal LH2 activities have been implicated in various diseases including Bruck syndrome, the molecular basis of the pathologies is still not well understood. Since LH2 null mice die at early embryonic stage, we generated LH2 heterozygous (LH2^+/-) mice in which LH2 level is significantly diminished, and characterized collagen and bone phenotypes using femurs. Compared to the wild-type (WT), LH2^+/- collagen showed a significant decrease in the ratio of hydroxylysine (Hyl)- to the Lys-aldehyde-derived collagen cross-links without affecting the total number of aldehydes involved in cross-links. Mass spectrometric analysis revealed that, in LH2^+/- type I collagen, the extent of hydroxylation of all telopeptidyl Lys residues was significantly decreased. In the helical domain, Lys hydroxylation at the cross-linking sites was either unaffected or slightly lower, but other sites were significantly diminished compared to WT. In LH2^+/- femurs, mineral densities of cortical and cancellous bones were significantly decreased and the mechanical properties of cortical bones evaluated by nanoindentation analysis were compromised. When cultured, LH2^+/- osteoblasts poorly produced mineralized nodules compared to WT osteoblasts. These data provide insight into the functionality of LH2 in collagen molecular phenotype and its critical role in bone matrix mineralization and mechanical properties.

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".

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.

Collagen Scaffolds Treated by Hydrogen Peroxide for Cell Cultivation

Collagen in the body is exposed to a range of influences, including free radicals, which can lead to a significant change in its structure. Modeling such an effect on collagen fibrils will allow one to get a native structure in vitro, which is important for modern tissue engineering. The aim of this work is to study the effect of free radicals on a solution of hydrogen peroxide with a concentration of 0.006–0.15% on the structure of collagen fibrils in vitro, and the response of cells to such treatment. SEM measurements show a decrease in the diameter of the collagen fibrils with an increase in the concentration of hydrogen peroxide. Such treatment also leads to an increase in the wetting angle of the collagen surface. Fourier transform infrared spectroscopy demonstrates a decrease in the signal with wave number 1084 cm⁻¹ due to the detachment of glucose and galactose linked to hydroxylysine, connected to the collagen molecule through the -C-O-C- group. During the first day of cultivating ASCs, MG-63, and A-431 cells, an increase in cell adhesion on collagen fibrils treated with H₂O₂ (0.015, 0.03%) was observed. Thus the effect of H₂O₂ on biologically relevant extracellular matrices for the formation of collagen scaffolds was shown.

Clinical and Molecular Delineation of Cutis Laxa Syndromes: Paradigms for Homeostasis

Cutis laxa (CL) syndromes are a large and heterogeneous group of rare connective tissue disorders that share loose redundant skin as a hallmark clinical feature, which reflects dermal elastic fiber fragmentation. Both acquired and congenital-Mendelian- forms exist. Acquired forms are progressive and often preceded by inflammatory triggers in the skin, but may show systemic elastolysis. Mendelian forms are often pleiotropic in nature and classified upon systemic manifestations and mode of inheritance. Though impaired elastogenesis is a common denominator in all Mendelian forms of CL, the underlying gene defects are diverse and affect structural components of the elastic fiber or impair metabolic pathways interfering with cellular trafficking, proline synthesis, or mitochondrial functioning. In this chapter we provide a detailed overview of the clinical and molecular characteristics of the different cutis laxa types and review the latest insights on elastic fiber assembly and homeostasis from both human and animal studies.

Exome sequencing identifies procollagen-lysine 2-oxoglutarate 5-dioxygenase 2 mutations in primary congenital and juvenile glaucoma

PURPOSE

To report the association of procollagen-lysine 2-oxoglutarate 5-dioxygenase 2 (PLOD2) mutations with bilateral primary congenital glaucoma (PCG) in monozygotic twins and with nondominant juvenile-onset primary open-angle glaucoma (JOAG).

METHODS

We utilized family-based whole-exome sequencing to detect disease-causing mutations in a pair of monozygotic twins with de-novo PCG and compared its existence in 50 nonfamilial cases of JOAG and 30 healthy controls. To validate the identified mutations, direct Sanger sequencing was performed. For further evaluation of gene expression in the ocular tissues, we performed whole-mount in situ hybridization in zebrafish embryos.

RESULTS

We identified a novel missense mutation (c.1925A>G, p.Tyr642Cys) in the PLOD2 gene in the monozygotic twin pair with PCG and another missense mutation (c.1880G>A, p.Arg627Gln) in one JOAG patient. Both mutations identified were heterozygous. Neither the parents of the twins nor the parents of the JOAG patient harbored the mutation and it was probably a de-novo change. The zebrafish in situ hybridization revealed expression of the PLOD2 gene during embryogenesis of the eye.

CONCLUSION

We observed an association of PLOD2 mutations with PCG and with nonfamilial JOAG. This new gene needs to be further investigated for its role in pathways associated with glaucoma pathogenesis.

Zebrafish Models for Human Skeletal Disorders

In 2019, the Nosology Committee of the International Skeletal Dysplasia Society provided an updated version of the Nosology and Classification of Genetic Skeletal Disorders. This is a reference list of recognized diseases in humans and their causal genes published to help clinician diagnosis and scientific research advances. Complementary to mammalian models, zebrafish has emerged as an interesting species to evaluate chemical treatments against these human skeletal disorders. Due to its versatility and the low cost of experiments, more than 80 models are currently available. In this article, we review the state-of-art of this “aquarium to bedside” approach describing the models according to the list provided by the Nosology Committee. With this, we intend to stimulate research in the appropriate direction to efficiently meet the actual needs of clinicians under the scope of the Nosology Committee.

Loss of zebrafish atp6v1e1b, encoding a subunit of vacuolar ATPase, recapitulates human ARCL type 2C syndrome and identifies multiple pathobiological signatures

The inability to maintain a strictly regulated endo(lyso)somal acidic pH through the proton-pumping action of the vacuolar-ATPases (v-ATPases) has been associated with various human diseases including heritable connective tissue disorders. Autosomal recessive (AR) cutis laxa (CL) type 2C syndrome is associated with genetic defects in the ATP6V1E1 gene and is characterized by skin wrinkles or loose redundant skin folds with pleiotropic systemic manifestations. The underlying pathological mechanisms leading to the clinical presentations remain largely unknown. Here, we show that loss of atp6v1e1b in zebrafish leads to early mortality, associated with craniofacial dysmorphisms, vascular anomalies, cardiac dysfunction, N-glycosylation defects, hypotonia, and epidermal structural defects. These features are reminiscent of the phenotypic manifestations in ARCL type 2C patients. Our data demonstrates that loss of atp6v1e1b alters endo(lyso)somal protein levels, and interferes with non-canonical v-ATPase pathways in vivo. In order to gain further insights into the processes affected by loss of atp6v1e1b, we performed an untargeted analysis of the transcriptome, metabolome, and lipidome in early atp6v1e1b-deficient larvae. We report multiple affected pathways including but not limited to oxidative phosphorylation, sphingolipid, fatty acid, and energy metabolism together with profound defects on mitochondrial respiration. Taken together, our results identify complex pathobiological effects due to loss of atp6v1e1b in vivo.

Cutis laxa syndromes are pleiotropic disorders of the connective tissue, characterized by skin redundancy and variable systemic manifestations. Cutis laxa syndromes are caused by pathogenic variants in genes encoding structural and regulatory components of the extracellular matrix or in genes encoding components of cellular trafficking, metabolism, and mitochondrial function. Pathogenic variants in genes coding for vacuolar-ATPases, a multisubunit complex responsible for the acidification of multiple intracellular vesicles, cause type 2 cutis laxa syndromes, a group of cutis laxa subtypes further characterized by neurological, skeletal, and rarely cardiopulmonary manifestations. To investigate the pathomechanisms of vacuolar-ATPase dysfunction, we generated zebrafish models that lack a crucial subunit of the vacuolar-ATPases. The mutant zebrafish models show morphological and functional features reminiscent of the phenotypic manifestations in cutis laxa patients carrying pathogenic variants in ATP6V1E1. In-depth analysis at multiple -omic levels identified biological signatures that indicate impairment of signaling pathways, lipid metabolism, and mitochondrial respiration. We anticipate that these data will contribute to a better understanding of the pathogenesis of cutis laxa syndromes and other disorders involving defective v-ATPase function, which may eventually improve patient treatment and management.

Bi-allelic premature truncating variants in LTBP1 cause cutis laxa syndrome

Latent transforming growth factor β (TGFβ)-binding proteins (LTBPs) are microfibril-associated proteins essential for anchoring TGFβ in the extracellular matrix (ECM) as well as for correct assembly of ECM components. Variants in LTBP2, LTBP3, and LTBP4 have been identified in several autosomal recessive Mendelian disorders with skeletal abnormalities with or without impaired development of elastin-rich tissues. Thus far, the human phenotype associated with LTBP1 deficiency has remained enigmatic. In this study, we report homozygous premature truncating LTBP1 variants in eight affected individuals from four unrelated consanguineous families. Affected individuals present with connective tissue features (cutis laxa and inguinal hernia), craniofacial dysmorphology, variable heart defects, and prominent skeletal features (craniosynostosis, short stature, brachydactyly, and syndactyly). In vitro studies on proband-derived dermal fibroblasts indicate distinct molecular mechanisms depending on the position of the variant in LTBP1. C-terminal variants lead to an altered LTBP1 loosely anchored in the microfibrillar network and cause increased ECM deposition in cultured fibroblasts associated with excessive TGFβ growth factor activation and signaling. In contrast, N-terminal truncation results in a loss of LTBP1 that does not alter TGFβ levels or ECM assembly. In vivo validation with two independent zebrafish lines carrying mutations in ltbp1 induce abnormal collagen fibrillogenesis in skin and intervertebral ligaments and ectopic bone formation on the vertebrae. In addition, one of the mutant zebrafish lines shows voluminous and hypo-mineralized vertebrae. Overall, our findings in humans and zebrafish show that LTBP1 function is crucial for skin and bone ECM assembly and homeostasis.

Transcriptomic meta-analysis of disuse muscle atrophy vs. resistance exercise-induced hypertrophy in young and older humans

BACKGROUND

Skeletal muscle atrophy manifests across numerous diseases; however, the extent of similarities/differences in causal mechanisms between atrophying conditions in unclear. Ageing and disuse represent two of the most prevalent and costly atrophic conditions, with resistance exercise training (RET) being the most effective lifestyle countermeasure. We employed gene-level and network-level meta-analyses to contrast transcriptomic signatures of disuse and RET, plus young and older RET to establish a consensus on the molecular features of, and therapeutic targets against, muscle atrophy in conditions of high socio-economic relevance.

METHODS

Integrated gene-level and network-level meta-analysis was performed on publicly available microarray data sets generated from young (18-35 years) m. vastus lateralis muscle subjected to disuse (unilateral limb immobilization or bed rest) lasting ≥7 days or RET lasting ≥3 weeks, and resistance-trained older (≥60 years) muscle.

RESULTS

Disuse and RET displayed predominantly separate transcriptional responses, and transcripts altered across conditions were mostly unidirectional. However, disuse and RET induced directly inverted expression profiles for mitochondrial function and translation regulation genes, with COX4I1, ENDOG, GOT2, MRPL12, and NDUFV2, the central hub components of altered mitochondrial networks, and ZMYND11, a hub gene of altered translation regulation. A substantial number of genes (n = 140) up-regulated post-RET in younger muscle were not similarly up-regulated in older muscle, with young muscle displaying a more pronounced extracellular matrix (ECM) and immune/inflammatory gene expression response. Both young and older muscle exhibited similar RET-induced ubiquitination/RNA processing gene signatures with associated PWP1, PSMB1, and RAF1 hub genes.

CONCLUSIONS

Despite limited opposing gene profiles, transcriptional signatures of disuse are not simply the converse of RET. Thus, the mechanisms of unloading cannot be derived from studying muscle loading alone and provides a molecular basis for understanding why RET fails to target all transcriptional features of disuse. Loss of RET-induced ECM mechanotransduction and inflammatory profiles might also contribute to suboptimal ageing muscle adaptations to RET. Disuse and age-dependent molecular candidates further establish a framework for understanding and treating disuse/ageing atrophy.

Plasticity of the skeleton and skeletal deformities in zebrafish (Danio rerio) linked to rearing density

The teleost zebrafish (Danio rerio), an established model for human skeletal diseases, is reared under controlled conditions with defined parameters for temperature and photoperiod. Studies aimed at defining the proper rearing density have been performed with regard to behavioural and physiological stress response, sex ratio and reproduction. Studies concerning the effect of rearing density on the skeletal phenotype are lacking. This study analyses the response of the skeleton to different rearing densities and describes the skeletal deformities. Wild-type zebrafish were reared up to 30 dpf (days post-fertilization) in a common environment. From 30 to 90 dpf, animals were reared at three different densities: high density (HD), 32 fish l^-1 ; medium density (MD), 8 fish l^-1 and low density (LD), 2 fish l^-1 . Animals at 30 and 90 dpf were collected and whole-mount stained with Alizarin red S to visualize mineralized tissues. The entire skeleton was analysed for meristic counts and 172 types of deformities. The results showed that the rearing density significantly influenced the specimens' average standard length, which decreased with the increase in the rearing density. Differences in meristic counts among the three groups were not observed. Rearing density-independent malformations affected the ribs, neural arches and the spines of the abdominal region, as well as vertebrae of the caudal complex. The HD group showed the highest number of deformities per specimen, the highest number of observed types of deformities and, together with the MD group, the highest frequency of specimens affected by severe deformities. In particular, the HD group showed deformities affecting arches, spines and vertebral centra in the caudal region of the vertebral column. This study provides evidence of an effect of the rearing density on the development of different skeletal phenotypes.

Skeletal Biology and Disease Modeling in Zebrafish

Zebrafish are teleosts (bony fish) that share with mammals a common ancestor belonging to the phylum Osteichthyes, from which their endoskeletal systems have been inherited. Indeed, teleosts and mammals have numerous genetically conserved features in terms of skeletal elements, ossification mechanisms, and bone matrix components in common. Yet differences related to bone morphology and function need to be considered when investigating zebrafish in skeletal research. In this review, we focus on zebrafish skeletal architecture with emphasis on the morphology of the vertebral column and associated anatomical structures. We provide an overview of the different ossification types and osseous cells in zebrafish and describe bone matrix composition at the microscopic tissue level with a focus on assessing mineralization. Processes of bone formation also strongly depend on loading in zebrafish, as we elaborate here. Furthermore, we illustrate the high regenerative capacity of zebrafish bones and present some of the technological advantages of using zebrafish as a model. We highlight zebrafish axial and fin skeleton patterning mechanisms, metabolic bone disease such as after immunosuppressive glucocorticoid treatment, as well as osteogenesis imperfecta (OI) and osteopetrosis research in zebrafish. We conclude with a view of why larval zebrafish xenografts are a powerful tool to study bone metastasis. © 2021 American Society for Bone and Mineral Research (ASBMR).

Abnormal Bone Collagen Cross-Linking in Osteogenesis Imperfecta/Bruck Syndrome Caused by Compound Heterozygous PLOD2 Mutations

Bruck syndrome (BS) is a congenital disorder characterized by joint flexion contractures, skeletal dysplasia, and increased bone fragility, which overlaps clinically with osteogenesis imperfecta (OI). On a genetic level, BS is caused by biallelic mutations in either FKBP10 or PLOD2. PLOD2 encodes the lysyl hydroxylase 2 (LH2) enzyme, which is responsible for the hydroxylation of cross‐linking lysine residues in fibrillar collagen telopeptide domains. This modification enables collagen to form chemically stable (permanent) intermolecular cross‐links in the extracellular matrix. Normal bone collagen develops a unique mix of such stable and labile lysyl‐oxidase–mediated cross‐links, which contribute to bone strength, resistance to microdamage, and crack propagation, as well as the ordered deposition of mineral nanocrystals within the fibrillar collagen matrix. Bone from patients with BS caused by biallelic FKBP10 mutations has been shown to have abnormal collagen cross‐linking; however, to date, no direct studies of human bone from BS caused by PLOD2 mutations have been reported. Here the results from a study of a 4‐year‐old boy with BS caused by compound heterozygous mutations in PLOD2 are discussed. Diminished hydroxylation of type I collagen telopeptide lysines but normal hydroxylation at triple‐helical sites was found. Consequently, stable trivalent cross‐links were essentially absent. Instead, allysine aldol dimeric cross‐links dominated as in normal skin collagen. Furthermore, in contrast to the patient's bone collagen, telopeptide lysines in cartilage type II collagen cross‐linked peptides from the patient's urine were normally hydroxylated. These findings shed light on the complex mechanisms that control the unique posttranslational chemistry and cross‐linking of bone collagen, and how, when defective, they can cause brittle bones and related connective tissue problems. © 2020 The Authors. JBMR Plus published by Wiley Periodicals LLC. on behalf of American Society for Bone and Mineral Research.

b3galt6 Knock-Out Zebrafish Recapitulate β3GalT6-Deficiency Disorders in Human and Reveal a Trisaccharide Proteoglycan Linkage Region

Proteoglycans are structurally and functionally diverse biomacromolecules found abundantly on cell membranes and in the extracellular matrix. They consist of a core protein linked to glycosaminoglycan chains via a tetrasaccharide linkage region. Here, we show that CRISPR/Cas9-mediated b3galt6 knock-out zebrafish, lacking galactosyltransferase II, which adds the third sugar in the linkage region, largely recapitulate the phenotypic abnormalities seen in human β3GalT6-deficiency disorders. These comprise craniofacial dysmorphism, generalized skeletal dysplasia, skin involvement and indications for muscle hypotonia. In-depth TEM analysis revealed disturbed collagen fibril organization as the most consistent ultrastructural characteristic throughout different affected tissues. Strikingly, despite a strong reduction in glycosaminoglycan content, as demonstrated by anion-exchange HPLC, subsequent LC-MS/MS analysis revealed a small amount of proteoglycans containing a unique linkage region consisting of only three sugars. This implies that formation of glycosaminoglycans with an immature linkage region is possible in a pathogenic context. Our study, therefore unveils a novel rescue mechanism for proteoglycan production in the absence of galactosyltransferase II, hereby opening new avenues for therapeutic intervention.

IFT20 is critical for collagen biosynthesis in craniofacial bone formation

Intraflagellar transport (IFT) is essential for assembling primary cilia required for bone formation. Disruption of IFT frequently leads to bone defects in humans. While it has been well studied about the function of IFT in osteogenic cell proliferation and differentiation, little is known about its role in collagen biosynthesis during bone formation. Here we show that IFT20, the smallest IFT protein in the IFT-B complex, is important for collagen biosynthesis in mice. Deletion of Ift20 in craniofacial osteoblasts displayed bone defects in the face. While collagen protein levels are unaffected by loss of Ift20, collagen cross-linking was significantly altered. In both Ift20:Wnt1-Cre and Ift20:Ocn-Cre mice the bones exhibit increased hydroxylysine-aldehyde deived cross-linking, and decreased lysine-aldehyde derived cross-linking. To obtain insight into the molecular mechanisms, we examined the expression levels of telopeptidyl lysyl hydroxylase 2 (LH2), and associated chaperone complexes. The results demonstrated that, while LH2 levels were unaffected by loss of Ift20, its chaperone, FKBP65, was significantly increased in Ift20:Wnt1-Cre and Ift20:Ocn-Cre mouse calvaria as well as femurs. These results suggest that IFT20 plays a pivotal role in collagen biosynthesis by regulating, in part, telopeptidyl lysine hydroxylation and cross-linking in bone. To the best of our knowledge, this is the first to demonstrate that the IFT components control collagen post-translational modifications. This provides a novel insight into the craniofacial bone defects associated with craniofacial skeletal ciliopathies.

Zebrafish: An Emerging Model for Orthopedic Research

Advances in next-generation sequencing have transformed our ability to identify genetic variants associated with clinical disorders of the musculoskeletal system. However, the means to functionally validate and analyze the physiological repercussions of genetic variation have lagged behind the rate of genetic discovery. The zebrafish provides an efficient model to leverage genetic analysis in an in vivo context. Its utility for orthopedic research is becoming evident in regard to both candidate gene validation as well as therapeutic discovery in tissues such as bone, tendon, muscle, and cartilage. With the development of new genetic and analytical tools to better assay aspects of skeletal tissue morphology, mineralization, composition, and biomechanics, researchers are emboldened to systematically approach how the skeleton develops and to identify the root causes, and potential treatments, of skeletal disease. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 38:925-936, 2020.

Phenomics-Based Quantification of CRISPR-Induced Mosaicism in Zebrafish

Genetic mosaicism can manifest as spatially variable phenotypes that vary from site to site within an organism. Here, we use imaging-based phenomics to quantitate phenotypes at many sites within the axial skeleton of CRISPR-edited G0 zebrafish. Through characterization of loss-of-function cell clusters in the developing skeleton, we identify a distinctive size distribution shown to arise from clonal fragmentation and merger events. We quantitate the phenotypic mosaicism produced by somatic mutations of two genes, plod2 and bmp1a, implicated in human osteogenesis imperfecta. Comparison of somatic, CRISPR-generated G0 mutants to homozygous germline mutants reveals phenotypic convergence, suggesting that CRISPR screens of G0 animals can faithfully recapitulate the biology of inbred disease models. We describe statistical frameworks for phenomic analysis of spatial phenotypic variation present in somatic G0 mutants. In sum, this study defines an approach for decoding spatially variable phenotypes generated during CRISPR-based screens.

Global Gene Expression Analysis Identifies Age-Related Differences in Knee Joint Transcriptome during the Development of Post-Traumatic Osteoarthritis in Mice

Aging and injury are two major risk factors for osteoarthritis (OA). Yet, very little is known about how aging and injury interact and contribute to OA pathogenesis. In the present study, we examined age- and injury-related molecular changes in mouse knee joints that could contribute to OA. Using RNA-seq, first we profiled the knee joint transcriptome of 10-week-old, 62-week-old, and 95-week-old mice and found that the expression of several inflammatory-response related genes increased as a result of aging, whereas the expression of several genes involved in cartilage metabolism decreased with age. To determine how aging impacts post-traumatic arthritis (PTOA) development, the right knee joints of 10-week-old and 62-week-old mice were injured using a non-invasive tibial compression injury model and injury-induced structural and molecular changes were assessed. At six-week post-injury, 62-week-old mice displayed significantly more cartilage degeneration and osteophyte formation compared with young mice. Although both age groups elicited similar transcriptional responses to injury, 62-week-old mice had higher activation of inflammatory cytokines than 10-week-old mice, whereas cartilage/bone metabolism genes had higher expression in 10-week-old mice, suggesting that the differential expression of these genes might contribute to the differences in PTOA severity observed between these age groups.

Zebrafish: A Resourceful Vertebrate Model to Investigate Skeletal Disorders

Animal models are essential tools for addressing fundamental scientific questions about skeletal diseases and for the development of new therapeutic approaches. Traditionally, mice have been the most common model organism in biomedical research, but their use is hampered by several limitations including complex generation, demanding investigation of early developmental stages, regulatory restrictions on breeding, and high maintenance cost. The zebrafish has been used as an efficient alternative vertebrate model for the study of human skeletal diseases, thanks to its easy genetic manipulation, high fecundity, external fertilization, transparency of rapidly developing embryos, and low maintenance cost. Furthermore, zebrafish share similar skeletal cells and ossification types with mammals. In the last decades, the use of both forward and new reverse genetics techniques has resulted in the generation of many mutant lines carrying skeletal phenotypes associated with human diseases. In addition, transgenic lines expressing fluorescent proteins under bone cell- or pathway- specific promoters enable in vivo imaging of differentiation and signaling at the cellular level. Despite the small size of the zebrafish, many traditional techniques for skeletal phenotyping, such as x-ray and microCT imaging and histological approaches, can be applied using the appropriate equipment and custom protocols. The ability of adult zebrafish to remodel skeletal tissues can be exploited as a unique tool to investigate bone formation and repair. Finally, the permeability of embryos to chemicals dissolved in water, together with the availability of large numbers of small-sized animals makes zebrafish a perfect model for high-throughput bone anabolic drug screening. This review aims to discuss the techniques that make zebrafish a powerful model to investigate the molecular and physiological basis of skeletal disorders.

Rmrp Mutation Disrupts Chondrogenesis and Bone Ossification in Zebrafish Model of Cartilage-Hair Hypoplasia via Enhanced Wnt/β-Catenin Signaling

Cartilage-hair hypoplasia (CHH) is an autosomal recessive metaphyseal chondrodysplasia characterized by bone dysplasia and many other highly variable features. The gene responsible for CHH is the RNA component of the mitochondrial RNA-processing endoribonuclease (RMRP) gene. Currently, the pathogenesis of osteochondrodysplasia and extraskeletal manifestations in CHH patients remains incompletely understood; in addition, there are no viable animal models for CHH. We generated an rmrp KO zebrafish model to study the developmental mechanisms of CHH. We found that rmrp is required for the patterning and shaping of pharyngeal arches. Rmrp mutation inhibits the intramembranous ossification of skull bones and promotes vertebrae ossification. The abnormalities of endochondral bone ossification are variable, depending on the degree of dysregulated chondrogenesis. Moreover, rmrp mutation inhibits cell proliferation and promotes apoptosis through dysregulating the expressions of cell-cycle- and apoptosis-related genes. We also demonstrate that rmrp mutation upregulates canonical Wnt/β-catenin signaling; the pharmacological inhibition of Wnt/β-catenin could partially alleviate the chondrodysplasia and increased vertebrae mineralization in rmrp mutants. Our study, by establishing a novel zebrafish model for CHH, partially reveals the underlying mechanism of CHH, hence deepening our understanding of the role of rmrp in skeleton development.

Using zebrafish to study skeletal genomics

While genome-wide association studies (GWAS) have revolutionized our understanding of the genetic architecture of skeletal diseases, animal models are required to identify causal mechanisms and to translate underlying biology into new therapies. Despite large-scale knockout mouse phenotyping efforts, the skeletal functions of most genes residing at GWAS-identified loci remain unknown, highlighting a need for complementary model systems to accelerate gene discovery. Over the past several decades, zebrafish (Danio rerio) has emerged as a powerful system for modeling the genetics of human diseases. In this review, our goal is to outline evidence supporting the utility of zebrafish for accelerating our understanding of human skeletal genomics, as well as gaps in knowledge that need to be filled for this purpose. We do this by providing a basic foundation of the zebrafish skeletal morphophysiology and phenotypes, and surveying evidence of skeletal gene homology and the use of zebrafish for post-GWAS analysis in other tissues and organs. We also outline challenges in translating zebrafish mutant phenotypes. Finally, we conclude with recommendations of future directions and how to leverage the large body of tools and knowledge of skeletal genetics in zebrafish for the needs of human skeletal genomic exploration. Due to their amenability to rapid genetic approaches, as well as the large number of conserved genetic and phenotypic features, there is a strong rationale supporting the use of zebrafish for human skeletal genomic studies.

Zebrafish and medaka as models for biomedical research of bone diseases

The identification of disease-causing mutations has in recent years progressed immensely due to whole genome sequencing approaches using patient material. The task accordingly is shifting from gene identification to functional analysis of putative disease-causing genes, preferably in an in vivo setting which also allows testing of drug candidates or biotherapeutics in whole animal disease models. In this review, we highlight the advances made in the field of bone diseases using small laboratory fish, focusing on zebrafish and medaka. We particularly highlight those human conditions where teleost models are available.

Zebrafish type I collagen mutants faithfully recapitulate human type I collagenopathies

The type I collagenopathies are a group of heterogeneous connective tissue disorders, that are caused by mutations in the genes encoding type I collagen and include specific forms of osteogenesis imperfecta (OI) and the Ehlers-Danlos syndrome (EDS). These disorders present with a broad disease spectrum and large clinical variability of which the underlying genetic basis is still poorly understood. In this study, we systematically analyzed skeletal phenotypes in a large set of zebrafish, with diverse mutations in the genes encoding type I collagen, representing different genetic forms of human OI, and a zebrafish model resembling human EDS, which harbors a number of soft connective tissues defects, typical of EDS. Furthermore, we provide insight into how zebrafish and human type I collagen are compositionally and functionally related, which is relevant in the interpretation of human type I collagen-related disease models. Our studies reveal a high degree of intergenotype variability in phenotypic expressivity that closely correlates with associated OI severity. Furthermore, we demonstrate the potential for select mutations to give rise to phenotypic variability, mirroring the clinical variability associated with human disease pathology. Therefore, our work suggests the future potential for zebrafish to aid in identifying unknown genetic modifiers and mechanisms underlying the phenotypic variability in OI and related disorders. This will improve diagnostic strategies and enable the discovery of new targetable pathways for pharmacological intervention.

Severely Impaired Bone Material Quality in Chihuahua Zebrafish Resembles Classical Dominant Human Osteogenesis Imperfecta

Excessive skeletal deformations and brittle fractures in the vast majority of patients suffering from osteogenesis imperfecta (OI) are a result of substantially reduced bone quality. Because the mechanical competence of bone is dependent on the tissue characteristics at small length scales, it is of crucial importance to assess how OI manifests at the micro- and nanoscale of bone. In this context, the Chihuahua (Chi/+) zebrafish, carrying a heterozygous glycine substitution in the α1 chain of collagen type I, has recently been proposed as a suitable animal model of classical dominant OI, showing skeletal deformities, altered mineralization patterns, and a smaller body size. This study assessed the bone quality properties of Chi/+ at multiple length scales using micro-computed tomography (micro-CT), histomorphometry, quantitative back-scattered electron imaging, Fourier-transform infrared spectroscopy, nanoindentation, and X-ray microscopy. At the skeletal level, the Chi/+ displays smaller body size, deformities, and fracture calli in the ribs. Morphological changes at the whole bone level showed that the vertebrae in Chi/+ had a smaller size, smaller thickness, and distorted shape. At the tissue level, Chi/+ displayed a higher degree of mineralization, lower collagen maturity, lower mineral maturity, altered osteoblast morphology, and lower osteocyte lacunar density compared to wild-type zebrafish. The alterations in the cellular, compositional, and structural properties of Chi/+ bones bear an explanation for the impaired local mechanical properties, which promote an increase in overall bone fragility in Chi/+. The quantitative assessment of bone quality in Chi/+ thus further validates this mutant as an important model reflecting osseous characteristics associated with human classical dominant OI. © 2018 American Society for Bone and Mineral Research.

Roles of PLODs in Collagen Synthesis and Cancer Progression

Collagen is the major component of extracellular matrix. Collagen cross-link and deposition depend on lysyl hydroxylation, which is catalyzed by procollagen-lysine, 2-oxoglutarate 5-dioxygenase (PLOD). Aberrant lysyl hydroxylation and collagen cross-link contributes to the progression of many collagen-related diseases, such as fibrosis and cancer. Three lysyl hydroxylases (LH1, LH2, and LH3) are identified, encoded by PLOD1, PLOD2, and PLOD3 genes. Expression of PLODs is regulated by multiple cytokines, transcription factors and microRNAs. Dysregulation of PLODs promotes cancer progression and metastasis, suggesting that targeting PLODs is potential strategy for cancer treatment. Here, we summarize the recent progress in the investigation of function and regulation of PLODs in normal tissue development and disease progression, especially in cancer.

Expanding the Clinical Spectrum of Phenotypes Caused by Pathogenic Variants in PLOD2

Osteogenesis imperfecta (OI) is a strikingly heterogeneous group of disorders with a broad range of phenotypic variations. It is also one of the differential diagnoses in bent bone dysplasias along with campomelic dysplasia and thanatophoric dysplasia and can usually be distinguished by decreased bone mineralization and bone fractures. Bent bone dysplasias also include syndromes such as kyphomelic dysplasia (MIM:211350) and mesomelic dysplasia Kozlowski-Reardon (MIM249710), both of which have been under debate regarding whether or not they are a real entity or simply a phenotypic manifestation of another dysplasia including OI. Bruck syndrome type 2 (BRKS2; MIM:609220) is a rare form of autosomal recessive OI caused by biallelic PLOD2 variants and is associated with congenital joint contractures with pterygia. In this report, we present six patients from four families with novel PLOD2 variants. All cases had multiple fractures. Other features ranged from prenatal lethal severe angulation of the long bones as in kyphomelic dysplasia and mesomelic dysplasia Kozlowski-Reardon through classical Bruck syndrome to moderate OI with normal joints. Two siblings with a kyphomelic dysplasia-like phenotype who were stillborn had compound heterozygous variants in PLOD2 (p.Asp585Val and p.Ser166*). One infant who succumbed at age 4 months had a bent bone phenotype phenotypically like skeletal dysplasia Kozlowski-Reardon (with mesomelic shortening, camptodactyly, retrognathia, cleft palate, skin dimples, but also with fractures). He was homozygous for the nonsense variant (p.Trp561*). Two siblings had various degrees of Bruck syndrome caused by the homozygous missense variant, p.His687Arg. Furthermore a boy with a clinical presentation of moderate OI had a possibly pathogenic homozygous variant p.Trp588Cys. Our experience of six patients with biallelic pathogenic variants in PLOD2 expands the phenotypic spectrum in the PLOD2-related phenotypes. © 2017 American Society for Bone and Mineral Research.

Novel Mutations in PLOD2 Cause Rare Bruck Syndrome

Bruck syndrome is a rare autosomal recessive form of osteogenesis imperfecta (OI), which is mainly characterized by joint contractures and recurrent fragility fractures. Mutations in FKBP10 and PLOD2 were identified as the underlying genetic defects of Bruck syndrome. Here we investigated the phenotypes and the pathogenic mutations of three unrelated Chinese patients with Bruck syndrome. Clinical fractures, bone mineral density (BMD), bone turnover biomarkers, and skeletal images were evaluated in detail. The pathogenic mutations were identified by targeted next-generation sequencing and subsequently confirmed by Sanger sequencing and cosegregation analysis. We also evaluated the effects of zoledronic acid on bone fracture incidence and BMD of the patients. Three patients had congenital joint contractures, recurrent fragility fractures, camptodactyly, clubfoot, scoliosis, but without dentinogenesis imperfecta and hearing loss. Five novel heterozygous mutations were detected in PLOD2, including three heterozygous missense mutations (c.1138C>T, p.Arg380Cys; c.1153T>C, p.Cys385Arg; and c.1982G>A, p.Gly661Asp), one heterozygous nonsense mutation (c.2038C>T, p.Arg680X), and one heterozygous splice-site mutation (c.503-2A>G). Their parents were all heterozygous carriers of these mutations in PLOD2. No clear genotype-phenotype correlations were found in these patients with PLOD2 mutations. Z-score of BMD was significantly increased, but scoliosis progressed and new bone fractures occurred during the treatment of zoledronic acid. Our findings expanded the spectrum of gene mutations of Bruck syndrome.

The chaperone activity of 4PBA ameliorates the skeletal phenotype of Chihuahua, a zebrafish model for dominant osteogenesis imperfecta

Classical osteogenesis imperfecta (OI) is a bone disease caused by type I collagen mutations and characterized by bone fragility, frequent fractures in absence of trauma and growth deficiency. No definitive cure is available for OI and to develop novel drug therapies, taking advantage of a repositioning strategy, the small teleost zebrafish (Danio rerio) is a particularly appealing model. Its small size, high proliferative rate, embryo transparency and small amount of drug required make zebrafish the model of choice for drug screening studies, when a valid disease model is available. We performed a deep characterization of the zebrafish mutant Chihuahua, that carries a G574D (p.G736D) substitution in the α1 chain of type I collagen. We successfully validated it as a model for classical OI. Growth of mutants was delayed compared with WT. X-ray, µCT, alizarin red/alcian blue and calcein staining revealed severe skeletal deformity, presence of fractures and delayed mineralization. Type I collagen extracted from different tissues showed abnormal electrophoretic migration and low melting temperature. The presence of endoplasmic reticulum (ER) enlargement due to mutant collagen retention in osteoblasts and fibroblasts of mutant fish was shown by electron and confocal microscopy. Two chemical chaperones, 4PBA and TUDCA, were used to ameliorate the cellular stress and indeed 4PBA ameliorated bone mineralization in larvae and skeletal deformities in adult, mainly acting on reducing ER cisternae size and favoring collagen secretion. In conclusion, our data demonstrated that ER stress is a novel target to ameliorate OI phenotype; chemical chaperones such as 4PBA may be, alone or in combination, a new class of molecules to be further investigated for OI treatment.

Heritable Skeletal Disorders Arising from Defects in Processing and Transport of Type I Procollagen from the ER: Perspectives on Possible Therapeutic Approaches

Rare bone disorders are a heterogeneous group of diseases, initially associated with mutations in type I procollagen (PC) genes. Recent developments from dissection at the molecular and cellular level have expanded the list of disease-causing proteins, revealing that disruption of the machinery that handles protein secretion can lead to failure in PC secretion and in several cases result in skeletal dysplasia. In parallel, cell-based in vitro studies of PC trafficking pathways offer clues to the identification of new disease candidate genes. Together, this raises the prospect of heritable bone disorders as a paradigm for biosynthetic protein traffic-related diseases, and an avenue through which therapeutic strategies can be explored.Here, we focus on human syndromes linked to defects in type I PC secretion with respect to the landscape of biosynthetic and protein transport steps within the early secretory pathway. We provide a perspective on possible therapeutic interventions for associated heritable craniofacial and skeletal disorders, considering different orders of complexity, from the cellular level by manipulation of proteostasis pathways to higher levels involving cell-based therapies for bone repair and regeneration.