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.

Non-coding deletions identify Maenli lncRNA as a limb-specific En1 regulator.

Long non-coding RNAs (lncRNAs) can be important components in gene-regulatory networks1, but the exact nature and extent of their involvement in human Mendelian disease is largely unknown. Here we show that genetic ablation of a lncRNA locus on human chromosome 2 causes a severe congenital limb malformation. We identified homozygous 27-63-kilobase deletions located 300 kilobases upstream of the engrailed-1 gene (EN1) in patients with a complex limb malformation featuring mesomelic shortening, syndactyly and ventral nails (dorsal dimelia). Re-engineering of the human deletions in mice resulted in a complete loss of En1 expression in the limb and a double dorsal-limb phenotype that recapitulates the human disease phenotype. Genome-wide transcriptome analysis in the developing mouse limb revealed a four-exon-long non-coding transcript within the deleted region, which we named Maenli. Functional dissection of the Maenli locus showed that its transcriptional activity is required for limb-specific En1 activation in cis, thereby fine-tuning the gene-regulatory networks controlling dorso-ventral polarity in the developing limb bud. Its loss results in the En1-related dorsal ventral limb phenotype, a subset of the full En1-associated phenotype. Our findings demonstrate that mutations involving lncRNA loci can result in human Mendelian disease.

Functional dissection of the Sox9-Kcnj2 locus identifies nonessential and instructive roles of TAD architecture.

The genome is organized in three-dimensional units called topologically associating domains (TADs), through a process dependent on the cooperative action of cohesin and the DNA-binding factor CTCF. Genomic rearrangements of TADs have been shown to cause gene misexpression and disease, but genome-wide depletion of CTCF has no drastic effects on transcription. Here, we investigate TAD function in vivo in mouse limb buds at the Sox9-Kcnj2 locus. We show that the removal of all major CTCF sites at the boundary and within the TAD resulted in a fusion of neighboring TADs, without major effects on gene expression. Gene misexpression and disease phenotypes, however, were achieved by redirecting regulatory activity through inversions and/or the repositioning of boundaries. Thus, TAD structures provide robustness and precision but are not essential for developmental gene regulation. Aberrant disease-related gene activation is not induced by a mere loss of insulation but requires CTCF-dependent redirection of enhancer-promoter contacts.

GORAB scaffolds COPI at the trans-Golgi for efficient enzyme recycling and correct protein glycosylation.

COPI is a key mediator of protein trafficking within the secretory pathway. COPI is recruited to the membrane primarily through binding to Arf GTPases, upon which it undergoes assembly to form coated transport intermediates responsible for trafficking numerous proteins, including Golgi-resident enzymes. Here, we identify GORAB, the protein mutated in the skin and bone disorder gerodermia osteodysplastica, as a component of the COPI machinery. GORAB forms stable domains at the trans-Golgi that, via interactions with the COPI-binding protein Scyl1, promote COPI recruitment to these domains. Pathogenic GORAB mutations perturb Scyl1 binding or GORAB assembly into domains, indicating the importance of these interactions. Loss of GORAB causes impairment of COPI-mediated retrieval of trans-Golgi enzymes, resulting in a deficit in glycosylation of secretory cargo proteins. Our results therefore identify GORAB as a COPI scaffolding factor, and support the view that defective protein glycosylation is a major disease mechanism in gerodermia osteodysplastica.

Dynamic 3D chromatin architecture contributes to enhancer specificity and limb morphogenesis.

The regulatory specificity of enhancers and their interaction with gene promoters is thought to be controlled by their sequence and the binding of transcription factors. By studying Pitx1, a regulator of hindlimb development, we show that dynamic changes in chromatin conformation can restrict the activity of enhancers. Inconsistent with its hindlimb-restricted expression, Pitx1 is controlled by an enhancer (Pen) that shows activity in forelimbs and hindlimbs. By Capture Hi-C and three-dimensional modeling of the locus, we demonstrate that forelimbs and hindlimbs have fundamentally different chromatin configurations, whereby Pen and Pitx1 interact in hindlimbs and are physically separated in forelimbs. Structural variants can convert the inactive into the active conformation, thereby inducing Pitx1 misexpression in forelimbs, causing partial arm-to-leg transformation in mice and humans. Thus, tissue-specific three-dimensional chromatin conformation can contribute to enhancer activity and specificity in vivo and its disturbance can result in gene misexpression and disease.

Impaired proteoglycan glycosylation, elevated TGF-ß signaling, and abnormal osteoblast differentiation as the basis for bone fragility in a mouse model for gerodermia osteodysplastica.

Gerodermia osteodysplastica (GO) is characterized by skin laxity and early-onset osteoporosis. GORAB, the responsible disease gene, encodes a small Golgi protein of poorly characterized function. To circumvent neonatal lethality of the GorabNull full knockout, Gorab was conditionally inactivated in mesenchymal progenitor cells (Prx1-cre), pre-osteoblasts (Runx2-cre), and late osteoblasts/osteocytes (Dmp1-cre), respectively. While in all three lines a reduction in trabecular bone density was evident, only GorabPrx1 and GorabRunx2 mutants showed dramatically thinned, porous cortical bone and spontaneous fractures. Collagen fibrils in the skin of GorabNull mutants and in bone of GorabPrx1 mutants were disorganized, which was also seen in a bone biopsy from a GO patient. Measurement of glycosaminoglycan contents revealed a reduction of dermatan sulfate levels in skin and cartilage from GorabNull mutants. In bone from GorabPrx1 mutants total glycosaminoglycan levels and the relative percentage of dermatan sulfate were both strongly diminished. Accordingly, the proteoglycans biglycan and decorin showed reduced glycanation. Also in cultured GORAB-deficient fibroblasts reduced decorin glycanation was evident. The Golgi compartment of these cells showed an accumulation of decorin, but reduced signals for dermatan sulfate. Moreover, we found elevated activation of TGF-ß in GorabPrx1 bone tissue leading to enhanced downstream signalling, which was reproduced in GORAB-deficient fibroblasts. Our data suggest that the loss of Gorab primarily perturbs pre-osteoblasts. GO may be regarded as a congenital disorder of glycosylation affecting proteoglycan synthesis due to delayed transport and impaired posttranslational modification in the Golgi compartment.

Examining tissue composition, whole-bone morphology and mechanical behavior of GorabPrx1 mice tibiae: A mouse model of premature aging.

Gerodermia osteodysplastica (GO) is a segmental progeroid disorder caused by loss-of-function mutations in the GORAB gene, associated with early onset osteoporosis and bone fragility. A conditional mouse model of GO (GorabPrx1) was generated in which the Gorab gene was deleted in long bones. We examined the biomechanical/functional relevance of the GorabPrx1 mutants as a premature aging model by characterizing bone composition, tissue-level strains, and whole-bone morphology and mechanical properties of the tibia. MicroCT imaging showed that GorabPrx1 tibiae had an increased anterior convex curvature and decreased cortical cross-sectional area, cortical thickness and moments of inertia, compared to littermate control (LC) tibiae. Fourier transform infrared (FTIR) imaging indicated a 34% decrease in mineral/matrix ratio and a 27% increase in acid phosphate content in the posterior metaphyseal cortex of the GorabPrx1 tibiae (p < .05), suggesting delayed mineralization. In vivo strain gauge measurement and finite element analysis showed ∼two times higher tissue-level strains within the GorabPrx1 tibiae relative to LC tibiae when subjected to axial compressive loads of the same magnitude. Three-point bending tests suggested that GorabPrx1 tibiae were weaker and more brittle, as indicated by decreasing whole-bone strength (46%), stiffness (55%), work-to-fracture (61%) and post-yield displacement (47%). Many of these morphological and biomechanical characteristics of the GorabPrx1 tibia recapitulated changes in other animal models of skeletal aging. Future studies are necessary to confirm how our observations might guide the way to a better understanding and treatment of GO.

Composition and dosage of a multipartite enhancer cluster control developmental expression of Ihh (Indian hedgehog).

Copy number variations (CNVs) often include noncoding sequences and putative enhancers, but how these rearrangements induce disease is poorly understood. Here we investigate CNVs involving the regulatory landscape of IHH (encoding Indian hedgehog), which cause multiple, highly localized phenotypes including craniosynostosis and synpolydactyly. We show through transgenic reporter and genome-editing studies in mice that Ihh is regulated by a constellation of at least nine enhancers with individual tissue specificities in the digit anlagen, growth plates, skull sutures and fingertips. Consecutive deletions, resulting in growth defects of the skull and long bones, showed that these enhancers function in an additive manner. Duplications, in contrast, caused not only dose-dependent upregulation but also misexpression of Ihh, leading to abnormal phalanges, fusion of sutures and syndactyly. Thus, precise spatiotemporal control of developmental gene expression is achieved by complex multipartite enhancer ensembles. Alterations in the composition of such clusters can result in gene misexpression and disease.

Formation of new chromatin domains determines pathogenicity of genomic duplications.

Chromosome conformation capture methods have identified subchromosomal structures of higher-order chromatin interactions called topologically associated domains (TADs) that are separated from each other by boundary regions. By subdividing the genome into discrete regulatory units, TADs restrict the contacts that enhancers establish with their target genes. However, the mechanisms that underlie partitioning of the genome into TADs remain poorly understood. Here we show by chromosome conformation capture (capture Hi-C and 4C-seq methods) that genomic duplications in patient cells and genetically modified mice can result in the formation of new chromatin domains (neo-TADs) and that this process determines their molecular pathology. Duplications of non-coding DNA within the mouse Sox9 TAD (intra-TAD) that cause female to male sex reversal in humans, showed increased contact of the duplicated regions within the TAD, but no change in the overall TAD structure. In contrast, overlapping duplications that extended over the next boundary into the neighbouring TAD (inter-TAD), resulted in the formation of a new chromatin domain (neo-TAD) that was isolated from the rest of the genome. As a consequence of this insulation, inter-TAD duplications had no phenotypic effect. However, incorporation of the next flanking gene, Kcnj2, in the neo-TAD resulted in ectopic contacts of Kcnj2 with the duplicated part of the Sox9 regulatory region, consecutive misexpression of Kcnj2, and a limb malformation phenotype. Our findings provide evidence that TADs are genomic regulatory units with a high degree of internal stability that can be sculptured by structural genomic variations. This process is important for the interpretation of copy number variations, as these variations are routinely detected in diagnostic tests for genetic disease and cancer. This finding also has relevance in an evolutionary setting because copy-number differences are thought to have a crucial role in the evolution of genome complexity.

Exome sequencing and CRISPR/Cas genome editing identify mutations of ZAK as a cause of limb defects in humans and mice.

The CRISPR/Cas technology enables targeted genome editing and the rapid generation of transgenic animal models for the study of human genetic disorders. Here we describe an autosomal recessive human disease in two unrelated families characterized by a split-foot defect, nail abnormalities of the hands, and hearing loss, due to mutations disrupting the SAM domain of the protein kinase ZAK. ZAK is a member of the MAPKKK family with no known role in limb development. We show that Zak is expressed in the developing limbs and that a CRISPR/Cas-mediated knockout of the two Zak isoforms is embryonically lethal in mice. In contrast, a deletion of the SAM domain induces a complex hindlimb defect associated with down-regulation of Trp63, a known split-hand/split-foot malformation disease gene. Our results identify ZAK as a key player in mammalian limb patterning and demonstrate the rapid utility of CRISPR/Cas genome editing to assign causality to human mutations in the mouse in <10 wk.

GORAB Missense Mutations Disrupt RAB6 and ARF5 Binding and Golgi Targeting.

Gerodermia osteodysplastica is a hereditary segmental progeroid disorder affecting skin, connective tissues, and bone that is caused by loss-of-function mutations in GORAB. The golgin, RAB6-interacting (GORAB) protein localizes to the Golgi apparatus and interacts with the small GTPase RAB6. In this study, we used different approaches to shed more light on the recruitment of GORAB to this compartment. We show that GORAB best colocalizes with trans-Golgi markers and is rapidly displaced upon Brefeldin A exposition, indicating a loose association with Golgi membranes. A yeast two-hybrid screening revealed a specific interaction with the small GTPase ADP-ribosylation factor (ARF5) in its active, GTP-bound form. ARF5 and RAB6 bind to GORAB via the same internal Golgi-targeting RAB6 and ARF5 binding (IGRAB) domain. Two GORAB missense mutations identified in gerodermia osteodysplastica patients fall within this IGRAB domain. GORAB carrying the mutation p.Ala220Pro had a cytoplasmic distribution and failed to interact with both RAB6 and ARF5. In contrast, the p.Ser175Phe mutation displaced GORAB from the Golgi compartment to vesicular structures and selectively impaired ARF5 binding. Our findings indicate that the IGRAB domain is crucial for the Golgi localization of GORAB and that loss of this localization impairs its physiological function.

Deletions, Inversions, Duplications: Engineering of Structural Variants using CRISPR/Cas in Mice.

Structural variations (SVs) contribute to the variability of our genome and are often associated with disease. Their study in model systems was hampered until now by labor-intensive genetic targeting procedures and multiple mouse crossing steps. Here we present the use of CRISPR/Cas for the fast (10 weeks) and efficient generation of SVs in mice. We specifically produced deletions, inversions, and also duplications at six different genomic loci ranging from 1.1 kb to 1.6 Mb with efficiencies up to 42%. After PCR-based selection, clones were successfully used to create mice via aggregation. To test the practicability of the method, we reproduced a human 500 kb disease-associated deletion and were able to recapitulate the human phenotype in mice. Furthermore, we evaluated the regulatory potential of a large genomic interval by deleting a 1.5 Mb fragment. The method presented permits rapid in vivo modeling of genomic rearrangements.

Transport activity and presence of ClC-7/Ostm1 complex account for different cellular functions.

Loss of the lysosomal ClC-7/Ostm1 2Cl(-)/H(+) exchanger causes lysosomal storage disease and osteopetrosis in humans and additionally changes fur colour in mice. Its conversion into a Cl(-) conductance in Clcn7(unc/unc) mice entails similarly severe lysosomal storage, but less severe osteopetrosis and no change in fur colour. To elucidate the basis for these phenotypical differences, we generated Clcn7(td/td) mice expressing an ion transport-deficient mutant. Their osteopetrosis was as severe as in Clcn7(-/-) mice, suggesting that the electric shunt provided by ClC-7(unc) can partially rescue osteoclast function. The normal coat colour of Clcn7(td/td) mice and their less severe neurodegeneration suggested that the ClC-7 protein, even when lacking measurable ion transport activity, is sufficient for hair pigmentation and that the conductance of ClC-7(unc) is harmful for neurons. Our in vivo structure-function analysis of ClC-7 reveals that both protein-protein interactions and ion transport must be considered in the pathogenesis of ClC-7-related diseases.

Properties of beta-propeller phytase expressed in transgenic tobacco.

Phytases are enzymes that liberate inorganic phosphates from phytate. In a previous study, a beta-propeller phytase (168phyA) from Bacillus subtilis was introduced into transgenic tobacco, which resulted in certain phenotypic changes. In the study described herein, the recombinant phytase (t168phyA) was purified from transgenic tobacco to near homogeneity by a three-step purification scheme. The biochemical properties and kinetic parameters of t168phyA were compared with those of its counterpart from B. subtilis. t168phyA was glycosylated, and it showed a 4 kDa increase in molecular size in SDS-PAGE (44 kDa vs. 40 kDa). Although its thermostability remained unchanged, its temperature optimum shifted from 60 degrees C to 45-50 degrees C and its pH optimum shifted from pH 5.5 to 6.0. Kinetic data showed that the t168phyA had a lower Kcat, but a higher Km than the native enzyme. Despite these changes, t168phyA remained catalytically active and has a specific activity of 2.3 U/mg protein. These results verify the activity of recombinant Bacillus phytase that is expressed in plants.