Repression and 3D-restructuring resolves regulatory conflicts in evolutionarily rearranged genomes.

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1's intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42's unresponsiveness. Rather, Zfp42's promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes.

Disruptions of topological chromatin domains cause pathogenic rewiring of gene-enhancer interactions.

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.

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.

Combinatorial microRNA activity is essential for the transition of pluripotent cells from proliferation into dormancy.

Dormancy is a key feature of stem cell function in adult tissues as well as in embryonic cells in the context of diapause. The establishment of dormancy is an active process that involves extensive transcriptional, epigenetic, and metabolic rewiring. How these processes are coordinated to successfully transition cells to the resting dormant state remains unclear. Here we show that microRNA activity, which is otherwise dispensable for preimplantation development, is essential for the adaptation of early mouse embryos to the dormant state of diapause. In particular, the pluripotent epiblast depends on miRNA activity, the absence of which results in the loss of pluripotent cells. Through the integration of high-sensitivity small RNA expression profiling of individual embryos and protein expression of miRNA targets with public data of protein-protein interactions, we constructed the miRNA-mediated regulatory network of mouse early embryos specific to diapause. We find that individual miRNAs contribute to the combinatorial regulation by the network, and the perturbation of the network compromises embryo survival in diapause. We further identified the nutrient-sensitive transcription factor TFE3 as an upstream regulator of diapause-specific miRNAs, linking cytoplasmic MTOR activity to nuclear miRNA biogenesis. Our results place miRNAs as a critical regulatory layer for the molecular rewiring of early embryos to establish dormancy.

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.

BMP4 triggers regulatory circuits specifying the cardiac mesoderm lineage.

Cardiac lineage specification in the mouse is controlled by TGFß and WNT signaling. From fly to fish, BMP has been identified as an indispensable heart inducer. A detailed analysis of the role of Bmp4 and its effectors Smad1/5, however, was still missing. We show that Bmp4 induces cardiac mesoderm formation in murine embryonic stem cells in vitro. Bmp4 first activates Wnt3 and upregulates Nodal. pSmad1/5 and the WNT effector Tcf3 form a complex, and together with pSmad2/3 activate mesoderm enhancers and Eomes. They then cooperate with Eomes to consolidate the expression of many mesoderm factors, including T. Eomes and T form a positive- feedback loop and open additional enhancers regulating early mesoderm genes, including the transcription factor Mesp1, establishing the cardiac mesoderm lineage. In parallel, the neural fate is suppressed. Our data confirm the pivotal role of Bmp4 in cardiac mesoderm formation in the mouse. We describe in detail the consecutive and cooperative actions of three signaling pathways, BMP, WNT and Nodal, and their effector transcription factors, during cardiac mesoderm specification.

Genome-wide identification of notochord enhancers comprising the regulatory landscape of the brachyury locus in mouse.

The node and notochord are important signaling centers organizing the dorso-ventral patterning of cells arising from neuro-mesodermal progenitors forming the embryonic body anlage. Owing to the scarcity of notochord progenitors and notochord cells, a comprehensive identification of regulatory elements driving notochord-specific gene expression has been lacking. Here, we have used ATAC-seq analysis of FACS-purified notochord cells from Theiler stage 12-13 mouse embryos to identify 8921 putative notochord enhancers. In addition, we established a new model for generating notochord-like cells in culture, and found 3728 of these enhancers occupied by the essential notochord control factors brachyury (T) and/or Foxa2. We describe the regulatory landscape of the T locus, comprising ten putative enhancers occupied by these factors, and confirmed the regulatory activity of three of these elements. Moreover, we characterized seven new elements by knockout analysis in embryos and identified one new notochord enhancer, termed TNE2. TNE2 cooperates with TNE in the trunk notochord, and is essential for notochord differentiation in the tail. Our data reveal an essential role of Foxa2 in directing T-expressing cells towards the notochord lineage.

Induction of kidney-related gene programs through co-option of SALL1 in mole ovotestes.

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.

Epigenetic regulator function through mouse gastrulation.

During ontogeny, proliferating cells become restricted in their fate through the combined action of cell-type-specific transcription factors and ubiquitous epigenetic machinery, which recognizes universally available histone residues or nucleotides in a context-dependent manner1,2. The molecular functions of these regulators are generally well understood, but assigning direct developmental roles to them is hampered by complex mutant phenotypes that often emerge after gastrulation3,4. Single-cell RNA sequencing and analytical approaches have explored this highly conserved, dynamic period across numerous model organisms5-8, including mouse9-18. Here we advance these strategies using a combined zygotic perturbation and single-cell RNA-sequencing platform in which many mutant mouse embryos can be assayed simultaneously, recovering robust  morphological and transcriptional information across a panel of ten essential regulators. Deeper analysis of central Polycomb repressive complex (PRC) 1 and 2 components indicates substantial cooperativity, but distinguishes a dominant role for PRC2 in restricting the germline. Moreover, PRC mutant phenotypes emerge after gross epigenetic and transcriptional changes within the initial conceptus prior to gastrulation. Our experimental framework may eventually lead to a fully quantitative view of how cellular diversity emerges using an identical genetic template and from a single totipotent cell.

Polycomb-mediated genome architecture enables long-range spreading of H3K27 methylation.

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.

A 37†kb region upstream of brachyury comprising a notochord enhancer is essential for notochord and tail development.

The node-streak border region comprising notochord progenitor cells (NPCs) at the posterior node and neuro-mesodermal progenitor cells (NMPs) in the adjacent epiblast is the prime organizing center for axial elongation in mouse embryos. The T-box transcription factor brachyury (T) is essential for both formation of the notochord and maintenance of NMPs, and thus is a key regulator of trunk and tail development. The T promoter controlling T expression in NMPs and nascent mesoderm has been characterized in detail; however, control elements for T expression in the notochord have not been identified yet. We have generated a series of deletion alleles by CRISPR/Cas9 genome editing in mESCs, and analyzed their effects in mutant mouse embryos. We identified a 37 kb region upstream of T that is essential for notochord function and tailbud outgrowth. Within that region, we discovered a T-binding enhancer required for notochord cell specification and differentiation. Our data reveal a complex regulatory landscape controlling cell type-specific expression and function of T in NMP/nascent mesoderm and node/notochord, allowing proper trunk and tail development.

A CRISPR-Cas9-engineered mouse model for GPI-anchor deficiency mirrors human phenotypes and exhibits hippocampal synaptic dysfunctions.

Pathogenic germline mutations in PIGV lead to glycosylphosphatidylinositol biosynthesis deficiency (GPIBD). Individuals with pathogenic biallelic mutations in genes of the glycosylphosphatidylinositol (GPI)-anchor pathway exhibit cognitive impairments, motor delay, and often epilepsy. Thus far, the pathophysiology underlying the disease remains unclear, and suitable rodent models that mirror all symptoms observed in human patients have not been available. Therefore, we used CRISPR-Cas9 to introduce the most prevalent hypomorphic missense mutation in European patients, Pigv:c.1022C > A (p.A341E), at a site that is conserved in mice. Mirroring the human pathology, mutant Pigv341E mice exhibited deficits in motor coordination, cognitive impairments, and alterations in sociability and sleep patterns, as well as increased seizure susceptibility. Furthermore, immunohistochemistry revealed reduced synaptophysin immunoreactivity in Pigv341E mice, and electrophysiology recordings showed decreased hippocampal synaptic transmission that could underlie impaired memory formation. In single-cell RNA sequencing, Pigv341E-hippocampal cells exhibited changes in gene expression, most prominently in a subtype of microglia and subicular neurons. A significant reduction in Abl1 transcript levels in several cell clusters suggested a link to the signaling pathway of GPI-anchored ephrins. We also observed elevated levels of Hdc transcripts, which might affect histamine metabolism with consequences for circadian rhythm. This mouse model will not only open the doors to further investigation into the pathophysiology of GPIBD, but will also deepen our understanding of the role of GPI-anchor-related pathways in brain development.

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.

Preformed chromatin topology assists transcriptional robustness of Shh during limb development.

Long-range gene regulation involves physical proximity between enhancers and promoters to generate precise patterns of gene expression in space and time. However, in some cases, proximity coincides with gene activation, whereas, in others, preformed topologies already exist before activation. In this study, we investigate the preformed configuration underlying the regulation of the Shh gene by its unique limb enhancer, the ZRS, in vivo during mouse development. Abrogating the constitutive transcription covering the ZRS region led to a shift within the Shh-ZRS contacts and a moderate reduction in Shh transcription. Deletion of the CTCF binding sites around the ZRS resulted in the loss of the Shh-ZRS preformed interaction and a 50% decrease in Shh expression but no phenotype, suggesting an additional, CTCF-independent mechanism of promoter-enhancer communication. This residual activity, however, was diminished by combining the loss of CTCF binding with a hypomorphic ZRS allele, resulting in severe Shh loss of function and digit agenesis. Our results indicate that the preformed chromatin structure of the Shh locus is sustained by multiple components and acts to reinforce enhancer-promoter communication for robust transcription.

BRACHYURY directs histone acetylation to target loci during mesoderm development.

T-box transcription factors play essential roles in multiple aspects of vertebrate development. Here, we show that cooperative function of BRACHYURY (T) with histone-modifying enzymes is essential for mouse embryogenesis. A single point mutation (TY88A) results in decreased histone 3 lysine 27 acetylation (H3K27ac) at T target sites, including the T locus, suggesting that T autoregulates the maintenance of its expression and functions by recruiting permissive chromatin modifications to putative enhancers during mesoderm specification. Our data indicate that T mediates H3K27ac recruitment through a physical interaction with p300. In addition, we determine that T plays a prominent role in the specification of hematopoietic and endothelial cell types. Hematopoietic and endothelial gene expression programs are disrupted in TY88A mutant embryos, leading to a defect in the differentiation of hematopoietic progenitors. We show that this role of T is mediated, at least in part, through activation of a distal Lmo2 enhancer.

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.

Noncoding copy-number variations are associated with congenital limb malformation.

PurposeCopy-number variants (CNVs) are generally interpreted by linking the effects of gene dosage with phenotypes. The clinical interpretation of noncoding CNVs remains challenging. We investigated the percentage of disease-associated CNVs in patients with congenital limb malformations that affect noncoding cis-regulatory sequences versus genes sensitive to gene dosage effects.MethodsWe applied high-resolution copy-number analysis to 340 unrelated individuals with isolated limb malformation. To investigate novel candidate CNVs, we re-engineered human CNVs in mice using clustered regularly interspaced short palindromic repeats (CRISPR)-based genome editing.ResultsOf the individuals studied, 10% harbored CNVs segregating with the phenotype in the affected families. We identified 31 CNVs previously associated with congenital limb malformations and four novel candidate CNVs. Most of the disease-associated CNVs (57%) affected the noncoding cis-regulatory genome, while only 43% included a known disease gene and were likely to result from gene dosage effects. In transgenic mice harboring four novel candidate CNVs, we observed altered gene expression in all cases, indicating that the CNVs had a regulatory effect either by changing the enhancer dosage or altering the topological associating domain architecture of the genome.ConclusionOur findings suggest that CNVs affecting noncoding regulatory elements are a major cause of congenital limb malformations.

The tissue-specific transcriptomic landscape of the mid-gestational mouse embryo.

Differential gene expression is a prerequisite for the formation of multiple cell types from the fertilized egg during embryogenesis. Understanding the gene regulatory networks controlling cellular differentiation requires the identification of crucial differentially expressed control genes and, ideally, the determination of the complete transcriptomes of each individual cell type. Here, we have analyzed the transcriptomes of six major tissues dissected from mid-gestational (TS12) mouse embryos. Approximately one billion reads derived by RNA-seq analysis provided extended transcript lengths, novel first exons and alternative transcripts of known genes. We have identified 1375 genes showing tissue-specific expression, providing gene signatures for each of the six tissues. In addition, we have identified 1403 novel putative long noncoding RNA gene loci, 439 of which show differential expression. Our analysis provides the first complete transcriptome data for the mouse embryo. It offers a rich data source for the analysis of individual genes and gene regulatory networks controlling mid-gestational development.

Different Concentrations of FGF Ligands, FGF2 or FGF8 Determine Distinct States of WNT-Induced Presomitic Mesoderm.

Presomitic mesoderm (PSM) cells are the precursors of the somites, which flank both sides of the neural tube and give rise to the musculo-skeletal system shaping the vertebrate body. WNT and FGF signaling control the formation of both the PSM and the somites and show a graded distribution with highest levels in the posterior PSM. We have used reporters for the mesoderm/PSM control genes T, Tbx6, and Msgn1 to investigate the differentiation of mouse ESCs from the naïve state via EpiSCs to PSM cells. Here we show that the activation of WNT signaling by CHIR99021 (CH) in combination with FGF ligand induces embryo-like PSM at high efficiency. By varying the FGF ligand concentration, the state of PSM cells formed can be altered. High FGF concentration supports posterior PSM formation, whereas low FGF generates anterior/differentiating PSM, in line with in vivo data. Furthermore, the level of Msgn1 expression depends on the FGF ligand concentration. We also show that Activin/Nodal signaling inhibits CH-mediated PSM induction in EpiSCs, without affecting T-expression. Inversely, Activin/Nodal inhibition enhances PSM induction by WNT/high FGF signaling. The ability to generate PSM cells of either posterior or anterior PSM identity with high efficiency in vitro will promote the investigation of the gene regulatory networks controlling the formation of nascent PSM cells and their switch to differentiating/somitic paraxial mesoderm. Stem Cells 2016;34:1790-1800.

Genome-wide association study and mouse expression data identify a highly conserved 32 kb intergenic region between WNT3 and WNT9b as possible susceptibility locus for isolated classic exstrophy of the bladder.

Bladder exstrophy-epispadias complex (BEEC), the severe end of the urorectal malformation spectrum, has a profound impact on continence as well as sexual and renal functions. It is widely accepted that for the majority of cases the genetic basis appears to be multifactorial. Here, we report the first study which utilizes genome-wide association methods to analyze a cohort comprising patients presenting the most common BEEC form, classic bladder exstrophy (CBE), to identify common variation associated with risk for isolated CBE. We employed discovery and follow-up samples comprising 218 cases/865 controls and 78 trios in total, all of European descent. Our discovery sample identified a marker near SALL1, showing genome-wide significant association with CBE. However, analyses performed on follow-up samples did not add further support to these findings. We were also able to identify an association with CBE across our study samples (discovery: P = 8.88 × 10(-5); follow-up: P = 0.0025; combined: 1.09 × 10(-6)) in a highly conserved 32 kb intergenic region containing regulatory elements between WNT3 and WNT9B. Subsequent analyses in mice revealed expression for both genes in the genital region during stages relevant to the development of CBE in humans. Unfortunately, we were not able to replicate the suggestive signal for WNT3 and WNT9B in a sample that was enriched for non-CBE BEEC cases (P = 0.51). Our suggestive findings support the hypothesis that larger samples are warranted to identify association of common variation with CBE.