Search

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.

Changes in gene expression represent an important source of phenotypic innovation. Yet, how such changes emerge and impact the evolution of traits remains elusive. Here, we explore the molecular mechanisms associated with the development of masculinizing ovotestes in female moles. By performing integrative analyses of epigenetic and transcriptional data in mole and mouse, we identified the co-option of SALL1 expression in mole ovotestes formation. Chromosome conformation capture analyses highlight a striking conservation of the 3D organization at the SALL1 locus, but an evolutionary divergence of enhancer activity. Interspecies reporter assays support the capability of mole-specific enhancers to activate transcription in urogenital tissues. Through overexpression experiments in transgenic mice, we further demonstrate the capability of SALL1 to induce kidney-related gene programs, which are a signature of mole ovotestes. Our results highlight the co-option of gene expression, through changes in enhancer activity, as a plausible mechanism for the evolution of traits.

Polycomb-group proteins play critical roles in gene silencing through the deposition of histone H3 lysine 27 trimethylation (H3K27me3) and chromatin compaction. This process is essential for embryonic stem cell (ESC) pluripotency, differentiation, and development. Polycomb repressive complex 2 (PRC2) can both read and write H3K27me3, enabling progressive spreading of H3K27me3 on the linear genome. Long-range Polycomb-associated DNA contacts have also been described, but their regulation and role in gene silencing remain unclear. Here, we apply H3K27me3 HiChIP, a protein-directed chromosome conformation method, and optical reconstruction of chromatin architecture to profile long-range Polycomb-associated DNA loops that span tens to hundreds of megabases across multiple topological associated domains in mouse ESCs and human induced pluripotent stem cells. We find that H3K27me3 loop anchors are enriched for Polycomb nucleation points and coincide with key developmental genes. Genetic deletion of H3K27me3 loop anchors results in disruption of spatial contact between distant loci and altered H3K27me3 in cis, both locally and megabases away on the same chromosome. In mouse embryos, loop anchor deletion leads to ectopic activation of the partner gene, suggesting that Polycomb-associated loops control gene silencing during development. Further, we find that alterations in PRC2 occupancy resulting from an RNA binding–deficient EZH2 mutant are accompanied by loss of Polycomb-associated DNA looping. Together, these results suggest PRC2 uses RNA binding to enhance long-range chromosome folding and H3K27me3 spreading. Developmental gene loci have unique roles in Polycomb spreading, emerging as important architectural elements of the epigenome.

Linking genomic variation to phenotypical traits remains a major challenge in evolutionary genetics. In this study, we use phylogenomic strategies to investigate a distinctive trait among mammals: the development of masculinizing ovotestes in female moles. By combining a chromosome-scale genome assembly of the Iberian mole, Talpa occidentalis, with transcriptomic, epigenetic, and chromatin interaction datasets, we identify rearrangements altering the regulatory landscape of genes with distinct gonadal expression patterns. These include a tandem triplication involving CYP17A1, a gene controlling androgen synthesis, and an intrachromosomal inversion involving the pro-testicular growth factor gene FGF9, which is heterochronically expressed in mole ovotestes. Transgenic mice with a knock-in mole CYP17A1 enhancer or overexpressing FGF9 showed phenotypes recapitulating mole sexual features. Our results highlight how integrative genomic approaches can reveal the phenotypic impact of noncoding sequence changes.

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.

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.

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.

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.

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

Mammalian genomes are organized into megabase-scale topologically associated domains (TADs). We demonstrate that disruption of TADs can rewire long-range regulatory architecture and result in pathogenic phenotypes. We show that distinct human limb malformations are caused by deletions, inversions, or duplications altering the structure of the TAD-spanning WNT6/IHH/EPHA4/PAX3 locus. Using CRISPR/Cas genome editing, we generated mice with corresponding rearrangements. Both in mouse limb tissue and patient-derived fibroblasts, disease-relevant structural changes cause ectopic interactions between promoters and non-coding DNA, and a cluster of limb enhancers normally associated with Epha4 is misplaced relative to TAD boundaries and drives ectopic limb expression of another gene in the locus. This rewiring occurred only if the variant disrupted a CTCF-associated boundary domain. Our results demonstrate the functional importance of TADs for orchestrating gene expression via genome architecture and indicate criteria for predicting the pathogenicity of human structural variants, particularly in non-coding regions of the human genome.