Targeting TGF-β signaling, oxidative stress, and cellular senescence rescues osteoporosis in gerodermia osteodysplastica

GORAB is a key regulator of Golgi vesicle transport and protein glycanation. Loss of GORAB function in gerodermia osteodysplastica (GO) causes shortening of glycosaminoglycan chains, leading to extracellular matrix disorganization that results in wrinkled skin, osteoporosis and elevated TGF-β signaling. In this study, we investigated the role of TGF-β-signaling, oxidative stress, and resulting cellular senescence in the osteoporosis phenotype of GO. Treatment of GorabPrx1 conditional knockouts with the TGF-β neutralizing antibody 1D11 rescued the trabecular bone loss, indicating that TGF-β overactivation causes osteoporosis in GO. Using an inducible knockout system, we demonstrated that TGF-β dysregulation was not a cell-intrinsic effect of GORAB inactivation, but a consequence of a disorganized extracellular matrix. Enhanced TGF-β signaling caused elevated Nox4 expression in GorabPrx1 mutants and in GO patients' fibroblasts, resulting in overproduction of mitochondrial superoxide. The resulting oxidative stress was detected in GORAB null cells and also in wildtype bystander cells. The same effect was observed in zebrafish after TALEN-mediated gorab inactivation, indicating that the pathway is evolutionarily conserved. Treating GorabPrx1 mutants with the antioxidant N-acetylcysteine ameliorated the osteoporosis phenotype. TGF-β induced oxidative stress coincided with accumulation of DNA damage and elevated expression of senescence markers. Inactivation of Cdkn2a in the GorabPrx1 rescued the osteoporosis phenotype. Reduced colony formation and altered subpopulations of bone marrow stromal cells were normalized upon inactivation of Cdkn2a, thus further demonstrating the relevance of cellular senescence in the pathogenesis. Our results shed light on the causative role of a TGF-β-Nox4-senescence axis and therapeutic strategies for GO.

Next-generation phenotyping integrated in a national framework for patients with ultrarare disorders improves genetic diagnostics and yields new molecular findings

Individuals with ultrarare disorders pose a structural challenge for healthcare systems since expert clinical knowledge is required to establish diagnoses. In TRANSLATE NAMSE, a 3-year prospective study, we evaluated a novel diagnostic concept based on multidisciplinary expertise in Germany. Here we present the systematic investigation of the phenotypic and molecular genetic data of 1,577 patients who had undergone exome sequencing and were partially analyzed with next-generation phenotyping approaches. Molecular genetic diagnoses were established in 32% of the patients totaling 370 distinct molecular genetic causes, most with prevalence below 1:50,000. During the diagnostic process, 34 novel and 23 candidate genotype-phenotype associations were identified, mainly in individuals with neurodevelopmental disorders. Sequencing data of the subcohort that consented to computer-assisted analysis of their facial images with GestaltMatcher could be prioritized more efficiently compared with approaches based solely on clinical features and molecular scores. Our study demonstrates the synergy of using next-generation sequencing and phenotyping for diagnosing ultrarare diseases in routine healthcare and discovering novel etiologies by multidisciplinary teams.

Multi-gene panel sequencing in highly consanguineous families and patients with congenital forms of skeletal dysplasias

Skeletal dysplasias (SKDs) are a heterogeneous group of more than 750 genetic disorders characterized by abnormal development, growth, and maintenance of bones or cartilage in the human skeleton. SKDs are often caused by variants in early patterning genes and in many cases part of multiple malformation syndromes and occur in combination with non-skeletal phenotypes. The aim of this study was to investigate the underlying genetic cause of congenital SKDs in highly consanguineous Pakistani families, as well as in sporadic and familial SKD cases from India using multigene panel sequencing analysis. Therefore, we performed panel sequencing of 386 bone-related genes in 7 highly consanguineous families from Pakistan and 27 cases from India affected with SKDs. In the highly consanguineous families, we were able to identify the underlying genetic cause in five out of seven families, resulting in a diagnostic yield of 71%. Whereas, in the sporadic and familial SKD cases, we identified 12 causative variants, corresponding to a diagnostic yield of 44%. The genetic heterogeneity in our cohorts was very high and we were able to detect various types of variants, including missense, nonsense, and frameshift variants, across multiple genes known to cause different types of SKDs. In conclusion, panel sequencing proved to be a highly effective way to decipher the genetic basis of SKDs in highly consanguineous families as well as sporadic and or familial cases from South Asia. Furthermore, our findings expand the allelic spectrum of skeletal dysplasias.

Loss-of-function variants affecting the STAGA complex component SUPT7L cause a developmental disorder with generalized lipodystrophy

Generalized lipodystrophy is a feature of various hereditary disorders, often leading to a progeroid appearance. In the present study we identified a missense and a frameshift variant in a compound heterozygous state in SUPT7L in a boy with intrauterine growth retardation, generalized lipodystrophy, and additional progeroid features. SUPT7L encodes a component of the transcriptional coactivator complex STAGA. By transcriptome sequencing, we showed the predicted missense variant to cause aberrant splicing, leading to exon truncation and thereby to a complete absence of SUPT7L in dermal fibroblasts. In addition, we found altered expression of genes encoding DNA repair pathway components. This pathway was further investigated and an increased rate of DNA damage was detected in proband-derived fibroblasts and genome-edited HeLa cells. Finally, we performed transient overexpression of wildtype SUPT7L in both cellular systems, which normalizes the number of DNA damage events. Our findings suggest SUPT7L as a novel disease gene and underline the link between genome instability and progeroid phenotypes.

Zebrafish as a model to investigate a biallelic gain-of-function variant in MSGN1, associated with a novel skeletal dysplasia syndrome

BACKGROUND/OBJECTIVES

Rare genetic disorders causing specific congenital developmental abnormalities often manifest in single families. Investigation of disease-causing molecular features are most times lacking, although these investigations may open novel therapeutic options for patients. In this study, we aimed to identify the genetic cause in an Iranian patient with severe skeletal dysplasia and to model its molecular function in zebrafish embryos.

RESULTS

The proband displays short stature and multiple skeletal abnormalities, including mesomelic dysplasia of the arms with complete humero-radio-ulna synostosis, arched clavicles, pelvic dysplasia, short and thin fibulae, proportionally short vertebrae, hyperlordosis and mild kyphosis. Exome sequencing of the patient revealed a novel homozygous c.374G > T, p.(Arg125Leu) missense variant in MSGN1 (NM_001105569). MSGN1, a basic-Helix-Loop-Helix transcription factor, plays a crucial role in formation of presomitic mesoderm progenitor cells/mesodermal stem cells during early developmental processes in vertebrates. Initial in vitro experiments show protein stability and correct intracellular localization of the novel variant in the nucleus and imply retained transcription factor function. To test the pathogenicity of the detected variant, we overexpressed wild-type and mutant msgn1 mRNA in zebrafish embryos and analyzed tbxta (T/brachyury/ntl). Overexpression of wild-type or mutant msgn1 mRNA significantly reduces tbxta expression in the tailbud compared to control embryos. Mutant msgn1 mRNA injected embryos depict a more severe effect, implying a gain-of-function mechanism. In vivo analysis on embryonic development was performed by clonal msgn1 overexpression in zebrafish embryos further demonstrated altered cell compartments in the presomitic mesoderm, notochord and pectoral fin buds. Detection of ectopic tbx6 and bmp2 expression in these embryos hint to affected downstream signals due to Msgn1 gain-of-function.

CONCLUSION

In contrast to loss-of-function effects described in animal knockdown models, gain-of-function of MSGN1 explains the only mildly affected axial skeleton of the proband and rather normal vertebrae. In this context we observed notochord bending and potentially disruption of pectoral fin buds/upper extremity after overexpression of msgn1 in zebrafish embryos. The latter might result from Msgn1 function on mesenchymal stem cells or on chondrogenesis in these regions. In addition, we detected ectopic tbx6 and bmp2a expression after gain of Msgn1 function in zebrafish, which are interconnected to short stature, congenital scoliosis, limb shortening and prominent skeletal malformations in patients. Our findings highlight a rare, so far undescribed skeletal dysplasia syndrome associated with a gain-of-function mutation in MSGN1 and hint to its molecular downstream effectors.

Single-cell, whole-embryo phenotyping of mammalian developmental disorders

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

AXIN1 bi-allelic variants disrupting the C-terminal DIX domain cause craniometadiaphyseal osteosclerosis with hip dysplasia

Sclerosing skeletal dysplasias result from an imbalance between bone formation and resorption. We identified three homozygous, C-terminally truncating AXIN1 variants in seven individuals from four families affected by macrocephaly, cranial hyperostosis, and vertebral endplate sclerosis. Other frequent findings included hip dysplasia, heart malformations, variable developmental delay, and hematological anomalies. In line with AXIN1 being a central component of the β-catenin destruction complex, analyses of primary and genome-edited cells harboring the truncating variants revealed enhanced basal canonical Wnt pathway activity. All three AXIN1-truncating variants resulted in reduced protein levels and impaired AXIN1 polymerization mediated by its C-terminal DIX domain but partially retained Wnt-inhibitory function upon overexpression. Addition of a tankyrase inhibitor attenuated Wnt overactivity in the AXIN1-mutant model systems. Our data suggest that AXIN1 coordinates the action of osteoblasts and osteoclasts and that tankyrase inhibitors can attenuate the effects of AXIN1 hypomorphic variants.

The enhancer landscape predetermines the skeletal regeneration capacity of stromal cells

Multipotent stromal cells are considered attractive sources for cell therapy and tissue engineering. Despite numerous experimental and clinical studies, broad application of stromal cell therapeutics is not yet emerging. A major challenge is the functional diversity of available cell sources. Here, we investigated the regenerative potential of clinically relevant human stromal cells from bone marrow (BMSCs), white adipose tissue, and umbilical cord compared with mature chondrocytes and skin fibroblasts in vitro and in vivo. Although all stromal cell types could express transcription factors related to endochondral ossification, only BMSCs formed cartilage discs in vitro that fully regenerated critical-size femoral defects after transplantation into mice. We identified cell type-specific epigenetic landscapes as the underlying molecular mechanism controlling transcriptional stromal differentiation networks. Binding sites of commonly expressed transcription factors in the enhancer and promoter regions of ossification-related genes, including Runt and bZIP families, were accessible only in BMSCs but not in extraskeletal stromal cells. This suggests an epigenetically predetermined differentiation potential depending on cell origin that allows common transcription factors to trigger distinct organ-specific transcriptional programs, facilitating forward selection of regeneration-competent cell sources. Last, we demonstrate that viable human BMSCs initiated defect healing through the secretion of osteopontin and contributed to transient mineralized bone hard callus formation after transplantation into immunodeficient mice, which was eventually replaced by murine recipient bone during final tissue remodeling.

Dissecting the influence of cellular senescence on cell mechanics and extracellular matrix formation in vitro

Tissue formation and healing both require cell proliferation and migration, but also extracellular matrix production and tensioning. In addition to restricting proliferation of damaged cells, increasing evidence suggests that cellular senescence also has distinct modulatory effects during wound healing and fibrosis. Yet, a direct role of senescent cells during tissue formation beyond paracrine signaling remains unknown. We here report how individual modules of the senescence program differentially influence cell mechanics and ECM expression with relevance for tissue formation. We compared DNA damage-mediated and DNA damage-independent senescence which was achieved through over-expression of either p16^Ink4a or p21^Cip1 cyclin-dependent kinase inhibitors in primary human skin fibroblasts. Cellular senescence modulated focal adhesion size and composition. All senescent cells exhibited increased single cell forces which led to an increase in tissue stiffness and contraction in an in vitro 3D tissue formation model selectively for p16 and p21-overexpressing cells. The mechanical component was complemented by an altered expression profile of ECM-related genes including collagens, lysyl oxidases, and MMPs. We found that particularly the lack of collagen and lysyl oxidase expression in the case of DNA damage-mediated senescence foiled their intrinsic mechanical potential. These observations highlight the active mechanical role of cellular senescence during tissue formation as well as the need to synthesize a functional ECM network capable of transferring and storing cellular forces.

CLCN7, a gene shared by autosomal recessive and autosomal dominant osteopetrosis

After the discovery of abundant v-ATPase complexes in the osteoclast ruffled membrane it was obvious that in parallel a negative counter-ion needs to be transported across this membrane to allow for efficient transport of protons into the resorption lacuna. While different candidate proteins were discussed the osteopetrosis phenotype of Clcn7 knockout mice suggested that the chloride/proton-exchanger ClC-7 might be responsible for transporting the negative charge. In the following, individuals with autosomal recessive osteopetrosis (ARO) were found to carry biallelic CLCN7 pathogenic variants. Shortly thereafter, heterozygous pathogenic variants were identified as the exclusive cause of autosomal dominant osteopetrosis type 2 (ADO2). Since in most cell types other than osteoclasts ClC-7 resides in late endosomes and lysosomes, it took some time until the electrophysiological properties of ClC-7 were elucidated. Whereas most missense variants lead to reduced chloride currents, several variants with accelerated kinetics have been identified. Evidence for folding problems is also known for several missense variants. Paradoxically, a heterozygous activating variant in ClC-7 was described to cause lysosomal alteration, pigmentation defects, and intellectual disability without osteopetrosis. The counter-intuitive 2 Cl^-/H⁺ exchange function of ClC-7 was shown to be physiologically important for intravesicular ion homeostasis. The lysosomal function of ClC-7 is also the reason why individuals with CLCN7-ARO can develop a storage disorder and neurodegeneration, a feature that is variable and difficult to predict. Furthermore, the low penetrance of heterozygous pathogenic CLCN7 variants and the clinical variability of ADO2 are incompletely understood. We aim to give an overview not only of the current knowledge about ClC-7 and its related pathologies, but also of the scientists and clinicians that paved the way for these discoveries.

Biallelic variants in ADAMTS15 cause a novel form of distal arthrogryposis

PURPOSE

We aimed to identify the underlying genetic cause for a novel form of distal arthrogryposis.

METHODS

Rare variant family-based genomics, exome sequencing, and disease-specific panel sequencing were used to detect ADAMTS15 variants in affected individuals. Adamts15 expression was analyzed at the single-cell level during murine embryogenesis. Expression patterns were characterized using in situ hybridization and RNAscope.

RESULTS

We identified homozygous rare variant alleles of ADAMTS15 in 5 affected individuals from 4 unrelated consanguineous families presenting with congenital flexion contractures of the interphalangeal joints and hypoplastic or absent palmar creases. Radiographic investigations showed physiological interphalangeal joint morphology. Additional features included knee, Achilles tendon, and toe contractures, spinal stiffness, scoliosis, and orthodontic abnormalities. Analysis of mouse whole-embryo single-cell sequencing data revealed a tightly regulated Adamts15 expression in the limb mesenchyme between embryonic stages E11.5 and E15.0. A perimuscular and peritendinous expression was evident in in situ hybridization in the developing mouse limb. In accordance, RNAscope analysis detected a significant coexpression with Osr1, but not with markers for skeletal muscle or joint formation.

CONCLUSION

In aggregate, our findings provide evidence that rare biallelic recessive trait variants in ADAMTS15 cause a novel autosomal recessive connective tissue disorder, resulting in a distal arthrogryposis syndrome.

Early-Onset Osteoporosis: Rare Monogenic Forms Elucidate the Complexity of Disease Pathogenesis Beyond Type I Collagen

Early-onset osteoporosis (EOOP), characterized by low bone mineral density (BMD) and fractures, affects children, premenopausal women and men aged <50 years. EOOP may be secondary to a chronic illness, long-term medication, nutritional deficiencies, etc. If no such cause is identified, EOOP is regarded primary and may then be related to rare variants in genes playing a pivotal role in bone homeostasis. If the cause remains unknown, EOOP is considered idiopathic. The scope of this review is to guide through clinical and genetic diagnostics of EOOP, summarize the present knowledge on rare monogenic forms of EOOP, and describe how analysis of bone biopsy samples can lead to a better understanding of the disease pathogenesis. The diagnostic pathway of EOOP is often complicated and extensive assessments may be needed to reliably exclude secondary causes. Due to the genetic heterogeneity and overlapping features in the various genetic forms of EOOP and other bone fragility disorders, the genetic diagnosis usually requires the use of next-generation sequencing to investigate several genes simultaneously. Recent discoveries have elucidated the complexity of disease pathogenesis both regarding genetic architecture and bone tissue-level pathology. Two rare monogenic forms of EOOP are due to defects in genes partaking in the canonical WNT pathway: LRP5 and WNT1. Variants in the genes encoding plastin-3 (PLS3) and sphingomyelin synthase 2 (SGMS2) have also been found in children and young adults with skeletal fragility. The molecular mechanisms leading from gene defects to clinical manifestations are often not fully understood. Detailed analysis of patient-derived transiliac bone biopsies gives valuable information to understand disease pathogenesis, distinguishes EOOP from other bone fragility disorders, and guides in patient management, but is not widely available in clinical settings. Despite the great advances in this field, EOOP remains an insufficiently explored entity and further research is needed to optimize diagnostic and therapeutic approaches. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

ClearCNV: CNV calling from NGS panel data in the presence of ambiguity and noise

MOTIVATION

While the identification of small variants in panel sequencing data can be considered a solved problem, the identification of larger, multi-exon copy number variants (CNVs) still poses a considerable challenge. Thus, CNV calling has not been established in all laboratories performing panel sequencing. At the same time, such laboratories have accumulated large datasets and thus have the need to identify CNVs on their data to close the diagnostic gap.

RESULTS

In this article, we present our method clearCNV that addresses this need in two ways. First, it helps laboratories to properly assign datasets to enrichment kits. Based on homogeneous subsets of data, clearCNV identifies CNVs affecting the targeted regions. Using real-world datasets and validation, we show that our method is highly competitive with previous methods and preferable in terms of specificity.

AVAILABILITY AND IMPLEMENTATION

The software is available for free under a permissible license at https://github.com/bihealth/clear-cnv.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Clinical Spectrum of Hereditary Hypophosphatemic Rickets With Hypercalciuria (HHRH)

Hereditary hypophosphatemic rickets with hypercalciuria (HHRH) represents an FGF23-independent disease caused by biallelic variants in the solute carrier family 34-member 3 (SLC34A3) gene. HHRH is characterized by chronic hypophosphatemia and an increased risk for nephrocalcinosis and rickets/osteomalacia, muscular weakness, and secondary limb deformity. Biochemical changes, but no relevant skeletal changes, have been reported for heterozygous SLC34A3 carriers. Therefore, we assessed the characteristics of individuals with biallelic and monoallelic SLC34A3 variants. In 8 index patients and 5 family members, genetic analysis was performed using a custom gene panel. The skeletal assessment comprised biochemical parameters, areal bone mineral density (aBMD), and bone microarchitecture. Pathogenic SLC34A3 variants were revealed in 7 of 13 individuals (2 homozygous, 5 heterozygous), whereas 3 of 13 carried monoallelic variants of unknown significance. Whereas both homozygous individuals had nephrocalcinosis, only one displayed a skeletal phenotype consistent with HHRH. Reduced to low-normal phosphate levels, decreased tubular reabsorption of phosphate (TRP), and high-normal to elevated values of 1,25-OH₂ -D₃ accompanied by normal cFGF23 levels were revealed independently of mutational status. Interestingly, individuals with nephrocalcinosis showed significantly increased calcium excretion and 1,25-OH₂ -D₃ levels but normal phosphate reabsorption. Furthermore, aBMD Z-score <-2.0 was revealed in 4 of 8 heterozygous carriers, and HR-pQCT analysis showed a moderate decrease in structural parameters. Our findings highlight the clinical relevance also of monoallelic SLC34A3 variants, including their potential skeletal impairment. Calcium excretion and 1,25-OH₂ -D₃ levels, but not TRP, were associated with nephrocalcinosis. Future studies should investigate the effects of distinct SLC34A3 variants and optimize treatment and monitoring regimens to prevent nephrocalcinosis and skeletal deterioration. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Genetic Diagnostics in Routine Osteological Assessment of Adult Low Bone Mass Disorders

CONTEXT

Many different inherited and acquired conditions can result in premature bone fragility/low bone mass disorders (LBMDs).

OBJECTIVE

We aimed to elucidate the impact of genetic testing on differential diagnosis of adult LBMDs and at defining clinical criteria for predicting monogenic forms.

METHODS

Four clinical centers broadly recruited a cohort of 394 unrelated adult women before menopause and men younger than 55 years with a bone mineral density (BMD) Z-score < -2.0 and/or pathological fractures. After exclusion of secondary causes or unequivocal clinical/biochemical hallmarks of monogenic LBMDs, all participants were genotyped by targeted next-generation sequencing.

RESULTS

In total, 20.8% of the participants carried rare disease-causing variants (DCVs) in genes known to cause osteogenesis imperfecta (COL1A1, COL1A2), hypophosphatasia (ALPL), and early-onset osteoporosis (LRP5, PLS3, and WNT1). In addition, we identified rare DCVs in ENPP1, LMNA, NOTCH2, and ZNF469. Three individuals had autosomal recessive, 75 autosomal dominant, and 4 X-linked disorders. A total of 9.7% of the participants harbored variants of unknown significance. A regression analysis revealed that the likelihood of detecting a DCV correlated with a positive family history of osteoporosis, peripheral fractures (> 2), and a high normal body mass index (BMI). In contrast, mutation frequencies did not correlate with age, prevalent vertebral fractures, BMD, or biochemical parameters. In individuals without monogenic disease-causing rare variants, common variants predisposing for low BMD (eg, in LRP5) were overrepresented.

CONCLUSION

The overlapping spectra of monogenic adult LBMD can be easily disentangled by genetic testing and the proposed clinical criteria can help to maximize the diagnostic yield.

Alternative splicing of BUD13 determines the severity of a developmental disorder with lipodystrophy and progeroid features

PURPOSE

In this study we aimed to identify the molecular genetic cause of a progressive multisystem disease with prominent lipodystrophy.

METHODS

In total, 5 affected individuals were investigated using exome sequencing. Dermal fibroblasts were characterized using RNA sequencing, proteomics, immunoblotting, immunostaining, and electron microscopy. Subcellular localization and rescue studies were performed.

RESULTS

We identified a lipodystrophy phenotype with a typical facial appearance, corneal clouding, achalasia, progressive hearing loss, and variable severity. Although 3 individuals showed stunted growth, intellectual disability, and died within the first decade of life (A1, A2, and A3), 2 are adults with normal intellectual development (A4 and A5). All individuals harbored an identical homozygous nonsense variant affecting the retention and splicing complex component BUD13. The nucleotide substitution caused alternative splicing of BUD13 leading to a stable truncated protein whose expression positively correlated with disease expression and life expectancy. In dermal fibroblasts, we found elevated intron retention, a global reduction of spliceosomal proteins, and nuclei with multiple invaginations, which were more pronounced in A1, A2, and A3. Overexpression of both BUD13 isoforms normalized the nuclear morphology.

CONCLUSION

Our results define a hitherto unknown syndrome and show that the alternative splice product converts a loss-of-function into a hypomorphic allele, thereby probably determining the severity of the disease and the survival of affected individuals.

Autosomal recessive hypophosphatemic rickets type 2 (ARHR2) due to ENPP1-deficiency

Awareness for hypophosphatemic rickets has increased in the last years, based on the availability of specific medical treatments. Autosomal recessive hypophosphatemic rickets type 2 (ARHR2) is a rare form of hypophosphatemic rickets, which is known to develop in survivors of generalized arterial calcification of infancy (GACI). Both disorders are based on a deficiency of ectonucleotide pyrophosphatase/phosphodiesterase 1 (ENPP1) and present with a high clinical variability and a lack of a phenotype-genotype association. ARHR2 is characterized by phosphate wasting due to elevated fibroblast growth factor 23 (FGF23) levels and might represent a response of the organism to minimize ectopic calcification in individuals with ENPP1-deficiency. This report reviews the recent clinical and preclinical data on this ultra-rare disease in childhood.

Efficient generation of osteoclasts from human induced pluripotent stem cells and functional investigations of lethal CLCN7-related osteopetrosis

Human induced pluripotent stem cells (hiPSCs) hold great potential for modeling human diseases and the development of innovative therapeutic approaches. Here, we report on a novel, simplified differentiation method for forming functional osteoclasts from hiPSCs. The three-step protocol starts with embryoid body formation, followed by hematopoietic specification, and finally osteoclast differentiation. We observed continuous production of monocyte-like cells over a period of up to 9 weeks, generating sufficient material for several osteoclast differentiations. The analysis of stage-specific gene and surface marker expression proved mesodermal priming, the presence of monocyte-like cells, and of terminally differentiated multinucleated osteoclasts, able to form resorption pits and trenches on bone and dentine in vitro. In comparison to peripheral blood mononuclear cell (PBMC)-derived osteoclasts hiPSC-derived osteoclasts were larger and contained a higher number of nuclei. Detailed functional studies on the resorption behavior of hiPSC-osteoclasts indicated a trend towards forming more trenches than pits and an increase in pseudoresorption. We used hiPSCs from an autosomal recessive osteopetrosis (ARO) patient (BIHi002-A, ARO hiPSCs) with compound heterozygous missense mutations p.(G292E) and p.(R403Q) in CLCN7, coding for the Cl^- /H⁺ -exchanger ClC-7, for functional investigations. The patient's leading clinical feature was a brain malformation due to defective neuronal migration. Mutant ClC-7 displayed residual expression and retained lysosomal co-localization with OSTM1, the gene coding for the osteopetrosis-associated transmembrane protein 1, but only ClC-7 harboring the mutation p.(R403Q) gave strongly reduced ion currents. An increased autophagic flux in spite of unchanged lysosomal pH was evident in undifferentiated ARO hiPSCs. ARO hiPSC-derived osteoclasts showed an increased size compared to hiPSCs of healthy donors. They were not able to resorb bone, underlining a loss-of-function effect of the mutations. In summary, we developed a highly reproducible, straightforward hiPSC-osteoclast differentiation protocol. We demonstrated that osteoclasts differentiated from ARO hiPSCs can be used as a disease model for ARO and potentially also other osteoclast-related diseases. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Expanding the clinical and molecular spectrum of ATP6V1A related metabolic cutis laxa

Several inborn errors of metabolism show cutis laxa as a highly recognizable feature. One group of these metabolic cutis laxa conditions is autosomal recessive cutis laxa type 2 caused by defects in v-ATPase components or the mitochondrial proline cycle. Besides cutis laxa, muscular hypotonia and cardiac abnormalities are hallmarks of autosomal recessive cutis laxa type 2D (ARCL2D) due to pathogenic variants in ATP6V1A encoding subunit A of the v-ATPase. Here, we report on three affected individuals from two families with ARCL2D in whom we performed whole exome and Sanger sequencing. We performed functional studies in fibroblasts from one individual, summarized all known probands' clinical, molecular, and biochemical features and compared them, also to other metabolic forms of cutis laxa. We identified novel missense and the first nonsense variant strongly affecting ATP6V1A expression. All six ARCL2D affected individuals show equally severe cutis laxa and dysmorphism at birth. While for one no information was available, two died in infancy and three are now adolescents with mild or absent intellectual disability. Muscular weakness, ptosis, contractures, and elevated muscle enzymes indicated a persistent myopathy. In cellular studies, a fragmented Golgi compartment, a delayed Brefeldin A-induced retrograde transport and glycosylation abnormalities were present in fibroblasts from two individuals. This is the second and confirmatory report on pathogenic variants in ATP6V1A as the cause of this extremely rare condition and the first to describe a nonsense allele. Our data highlight the tremendous clinical variability of ATP6V1A related phenotypes even within the same family.

Biallelic truncating variants in ATP9A cause a novel neurodevelopmental disorder involving postnatal microcephaly and failure to thrive

Background

Genes implicated in the Golgi and endosomal trafficking machinery are crucial for brain development, and mutations in them are particularly associated with postnatal microcephaly (POM).

Methods

Exome sequencing was performed in three affected individuals from two unrelated consanguineous families presenting with delayed neurodevelopment, intellectual disability of variable degree, POM and failure to thrive. Patient-derived fibroblasts were tested for functional effects of the variants.

Results

We detected homozygous truncating variants in ATP9A. While the variant in family A is predicted to result in an early premature termination codon, the variant in family B affects a canonical splice site. Both variants lead to a substantial reduction of ATP9A mRNA expression. It has been shown previously that ATP9A localises to early and recycling endosomes, whereas its depletion leads to altered gene expression of components from this compartment. Consistent with previous findings, we also observed overexpression of ARPC3 and SNX3, genes strongly interacting with ATP9A.

Conclusion

In aggregate, our findings show that pathogenic variants in ATP9A cause a novel autosomal recessive neurodevelopmental disorder with POM. While the physiological function of endogenous ATP9A is still largely elusive, our results underline a crucial role of this gene in endosomal transport in brain tissue.

Compound Heterozygous Frameshift Mutations in MESD Cause a Lethal Syndrome Suggestive of Osteogenesis Imperfecta Type XX

Multiple genes are known to be associated with osteogenesis imperfecta (OI), a phenotypically and genetically heterogenous bone disorder, marked predominantly by low bone mineral density and increased risk of fractures. Recently, mutations affecting MESD, which encodes for a chaperone required for trafficking of the low-density lipoprotein receptors LRP5 and LRP6 in the endoplasmic reticulum, were described to cause autosomal-recessive OI XX in homozygous children. In the present study, whole-exome sequencing of three stillbirths in one family was performed to evaluate the presence of a hereditary disorder. To further characterize the skeletal phenotype, fetal autopsy, bone histology, and quantitative backscattered electron imaging (qBEI) were performed, and the results were compared with those from an age-matched control with regular skeletal phenotype. In each of the affected individuals, compound heterozygous mutations in MESD exon 2 and exon 3 were detected. Based on the skeletal phenotype, which was characterized by multiple intrauterine fractures and severe skeletal deformity, OI XX was diagnosed in these individuals. Histological evaluation of MESD specimens revealed an impaired osseous development with an altered osteocyte morphology and reduced canalicular connectivity. Moreover, analysis of bone mineral density distribution by qBEI indicated an impaired and more heterogeneous matrix mineralization in individuals with MESD mutations than in controls. In contrast to the previously reported phenotypes of individuals with OI XX, the more severe phenotype in the present study is likely explained by a mutation in exon 2, located within the chaperone domain of MESD, that leads to a complete loss of function, which indicates the relevance of MESD in early skeletal development. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR)..

Relevant genetic variants are common in women with pregnancy and lactation-associated osteoporosis (PLO) and predispose to more severe clinical manifestations

Pregnancy and lactation-associated osteoporosis (PLO) is a rare skeletal disorder characterized by early-onset osteoporosis typically manifestating with vertebral compression fractures or transient osteoporosis of the hip. We hypothesized that genetic variants may play a role in the development of PLO. This study aimed to analyze the presence of genetic variants and a potential association with the clinical presentation in PLO. 42 women with PLO were included from 2013 to 2019 in a multicenter study in Germany. All cases underwent comprehensive genetic analysis based on a custom-designed gene panel including genes relevant for skeletal disorders. The skeletal status was assessed using dual-energy X-ray absorptiometry (DXA). Subgroups were further analyzed by serum bone turnover markers (n = 31) and high-resolution peripheral computed tomography (HR-pQCT; n = 23). We detected relevant genetic variants in 21 women (50%), with LRP5, WNT1 and COL1A1/A2 being the most commonly involved genes. The mean number of vertebral compression fractures was 3.3 ± 3.4 per case with a significantly higher occurrence in the subgroup with genetic variants (4.8 ± 3.7 vs. 1.8 ± 2.3, p = 0.02). Among the total cohort, DXA Z-scores were significantly lower at the lumbar spine compared to the femoral neck (p = 0.002). HR-pQCT revealed a pronounced reduction of trabecular and cortical thickness, while trabecular number was within the reference range. Eighteen women (43%) received a bone-specific therapy (primarily teriparatide). Overall, a steep increase in bone mass (+37.7%) was observed after 3 years. In conclusion, pregnancy and lactation represent skeletal risk factors, which may unmask hereditary bone disorders leading to PLO. These cases were affected more severely. Nevertheless, a timely diagnosis and adequate treatment can ensure a substantial recovery potential even without specific therapy. Patients with genetically induced low bone turnover (e.g.; LRP5, WNT1) may especially benefit from osteo-anabolic medication.