Regulation of Myc transcription by an enhancer cluster dedicated to pluripotency and early embryonic expression.

MYC plays various roles in pluripotent stem cells, including the promotion of somatic cell reprogramming to pluripotency, the regulation of cell competition and the control of embryonic diapause. However, how Myc expression is regulated in this context remains unknown. The Myc gene lies within a ~ 3-megabase gene desert with multiple cis-regulatory elements. Here we use genomic rearrangements, transgenesis and targeted mutation to analyse Myc regulation in early mouse embryos and pluripotent stem cells. We identify a topologically-associated region that homes enhancers dedicated to Myc transcriptional regulation in stem cells of the pre-implantation and early post-implantation embryo. Within this region, we identify elements exclusively dedicated to Myc regulation in pluripotent cells, with distinct enhancers that sequentially activate during naive and formative pluripotency. Deletion of pluripotency-specific enhancers dampens embryonic stem cell competitive ability. These results identify a topologically defined enhancer cluster dedicated to early embryonic expression and uncover a modular mechanism for the regulation of Myc expression in different states of pluripotency.

HYENA detects oncogenes activated by distal enhancers in cancer.

Somatic structural variations (SVs) in cancer can shuffle DNA content in the genome, relocate regulatory elements, and alter genome organization. Enhancer hijacking occurs when SVs relocate distal enhancers to activate proto-oncogenes. However, most enhancer hijacking studies have only focused on protein-coding genes. Here, we develop a computational algorithm "HYENA" to identify candidate oncogenes (both protein-coding and non-coding) activated by enhancer hijacking based on tumor whole-genome and transcriptome sequencing data. HYENA detects genes whose elevated expression is associated with somatic SVs by using a rank-based regression model. We systematically analyze 1,146 tumors across 25 types of adult tumors and identify a total of 108 candidate oncogenes including many non-coding genes. A long non-coding RNA TOB1-AS1 is activated by various types of SVs in 10% of pancreatic cancers through altered 3-dimensional genome structure. We find that high expression of TOB1-AS1 can promote cell invasion and metastasis. Our study highlights the contribution of genetic alterations in non-coding regions to tumorigenesis and tumor progression.

Shifting landscapes: the role of 3D genomic organizations in gene regulatory strategies.

3D genome folding enables the physical storage of chromosomes into the compact volume of a cell's nucleus, allows for the accurate segregation of chromatin to daughter cells, and has been shown to be tightly coupled to the way in which genetic information is converted into transcriptional programs [1-3]. Importantly, this link between chromatin architecture and gene regulation is a selectable feature in which modifications to chromatin organization accompany, or perhaps even drive the establishment of new regulatory strategies with enduring impacts on animal body plan complexity. Here, we discuss the nature of different 3D genome folding systems found across the tree of life, with particular emphasis on metazoans, and the relative influence of these systems on gene regulation. We suggest how the properties of these folding systems have influenced regulatory strategies employed by different lineages and may have catalyzed the partitioning and specialization of genetic programs that enabled multicellularity and organ-grade body plan complexity.

TATTOO-seq delineates spatial and cell type-specific regulatory programs in the developing limb.

The coordinated differentiation of progenitor cells into specialized cell types and their spatial organization into distinct domains is central to embryogenesis. Here, we developed and applied an unbiased spatially resolved single-cell transcriptomics method to identify the genetic programs underlying the emergence of specialized cell types during mouse limb development and their spatial integration. We identify multiple transcription factors whose expression patterns are predominantly associated with cell type specification or spatial position, suggesting two parallel yet highly interconnected regulatory systems. We demonstrate that the embryonic limb undergoes a complex multiscale reorganization upon perturbation of one of its spatial organizing centers, including the loss of specific cell populations, alterations of preexisting cell states' molecular identities, and changes in their relative spatial distribution. Our study shows how multidimensional single-cell, spatially resolved molecular atlases can allow the deconvolution of spatial identity and cell fate and reveal the interconnected genetic networks that regulate organogenesis and its reorganization upon genetic alterations.

Genetic dissection identifies Necdin as a driver gene in a mouse model of paternal 15q duplications.

Maternally inherited duplication of chromosome 15q11-q13 (Dup15q) is a pathogenic copy number variation (CNV) associated with autism spectrum disorder (ASD). Recently, paternally derived duplication has also been shown to contribute to the development of ASD. The molecular mechanism underlying paternal Dup15q remains unclear. Here, we conduct genetic and overexpression-based screening and identify Necdin (Ndn) as a driver gene for paternal Dup15q resulting in the development of ASD-like phenotypes in mice. An excess amount of Ndn results in enhanced spine formation and density as well as hyperexcitability of cortical pyramidal neurons. We generate 15q dupΔNdn mice with a normalized copy number of Ndn by excising its one copy from Dup15q mice using a CRISPR-Cas9 system. 15q dupΔNdn mice do not show ASD-like phenotypes and show dendritic spine dynamics and cortical excitatory-inhibitory balance similar to wild type animals. Our study provides an insight into the role of Ndn in paternal 15q duplication and a mouse model of paternal Dup15q syndrome.

Local retinoic acid signaling directs emergence of the extraocular muscle functional unit.

Coordinated development of muscles, tendons, and their attachment sites ensures emergence of functional musculoskeletal units that are adapted to diverse anatomical demands among different species. How these different tissues are patterned and functionally assembled during embryogenesis is poorly understood. Here, we investigated the morphogenesis of extraocular muscles (EOMs), an evolutionary conserved cranial muscle group that is crucial for the coordinated movement of the eyeballs and for visual acuity. By means of lineage analysis, we redefined the cellular origins of periocular connective tissues interacting with the EOMs, which do not arise exclusively from neural crest mesenchyme as previously thought. Using 3D imaging approaches, we established an integrative blueprint for the EOM functional unit. By doing so, we identified a developmental time window in which individual EOMs emerge from a unique muscle anlage and establish insertions in the sclera, which sets these muscles apart from classical muscle-to-bone type of insertions. Further, we demonstrate that the eyeballs are a source of diffusible all-trans retinoic acid (ATRA) that allow their targeting by the EOMs in a temporal and dose-dependent manner. Using genetically modified mice and inhibitor treatments, we find that endogenous local variations in the concentration of retinoids contribute to the establishment of tendon condensations and attachment sites that precede the initiation of muscle patterning. Collectively, our results highlight how global and site-specific programs are deployed for the assembly of muscle functional units with precise definition of muscle shapes and topographical wiring of their tendon attachments.

Signals from the brain and olfactory epithelium control shaping of the mammalian nasal capsule cartilage.

Facial shape is the basis for facial recognition and categorization. Facial features reflect the underlying geometry of the skeletal structures. Here, we reveal that cartilaginous nasal capsule (corresponding to upper jaw and face) is shaped by signals generated by neural structures: brain and olfactory epithelium. Brain-derived Sonic Hedgehog (SHH) enables the induction of nasal septum and posterior nasal capsule, whereas the formation of a capsule roof is controlled by signals from the olfactory epithelium. Unexpectedly, the cartilage of the nasal capsule turned out to be important for shaping membranous facial bones during development. This suggests that conserved neurosensory structures could benefit from protection and have evolved signals inducing cranial cartilages encasing them. Experiments with mutant mice revealed that the genomic regulatory regions controlling production of SHH in the nervous system contribute to facial cartilage morphogenesis, which might be a mechanism responsible for the adaptive evolution of animal faces and snouts.

Cnidarian Cell Type Diversity and Regulation Revealed by Whole-Organism Single-Cell RNA-Seq.

The emergence and diversification of cell types is a leading factor in animal evolution. So far, systematic characterization of the gene regulatory programs associated with cell type specificity was limited to few cell types and few species. Here, we perform whole-organism single-cell transcriptomics to map adult and larval cell types in the cnidarian Nematostella vectensis, a non-bilaterian animal with complex tissue-level body-plan organization. We uncover eight broad cell classes in Nematostella, including neurons, cnidocytes, and digestive cells. Each class comprises different subtypes defined by the expression of multiple specific markers. In particular, we characterize a surprisingly diverse repertoire of neurons, which comparative analysis suggests are the result of lineage-specific diversification. By integrating transcription factor expression, chromatin profiling, and sequence motif analysis, we identify the regulatory codes that underlie Nematostella cell-specific expression. Our study reveals cnidarian cell type complexity and provides insights into the evolution of animal cell-specific genomic regulation.

Author Correction: A Myc enhancer cluster regulates normal and leukaemic haematopoietic stem cell hierarchies.

In the originally published version of this Letter, ref. 43 was erroneously provided twice. In the 'Estimation of relative cell-type-specific composition of AML samples' section in the Methods, the citation to ref. 43 after the GEO dataset GSE24759 is correct. However, in the 'Mice' section of the Methods, the citation to ref. 43 after 'TAMERE' should have been associated with a new reference1. The original Letter has been corrected online (with the new reference included as ref. 49).

A Myc enhancer cluster regulates normal and leukaemic haematopoietic stem cell hierarchies.

The transcription factor Myc is essential for the regulation of haematopoietic stem cells and progenitors and has a critical function in haematopoietic malignancies. Here we show that an evolutionarily conserved region located 1.7 megabases downstream of the Myc gene that has previously been labelled as a 'super-enhancer' is essential for the regulation of Myc expression levels in both normal haematopoietic and leukaemic stem cell hierarchies in mice and humans. Deletion of this region in mice leads to a complete loss of Myc expression in haematopoietic stem cells and progenitors. This caused an accumulation of differentiation-arrested multipotent progenitors and loss of myeloid and B cells, mimicking the phenotype caused by Mx1-Cre-mediated conditional deletion of the Myc gene in haematopoietic stem cells. This super-enhancer comprises multiple enhancer modules with selective activity that recruits a compendium of transcription factors, including GFI1b, RUNX1 and MYB. Analysis of mice carrying deletions of individual enhancer modules suggests that specific Myc expression levels throughout most of the haematopoietic hierarchy are controlled by the combinatorial and additive activity of individual enhancer modules, which collectively function as a 'blood enhancer cluster' (BENC). We show that BENC is also essential for the maintenance of MLL-AF9-driven leukaemia in mice. Furthermore, a BENC module, which controls Myc expression in mouse haematopoietic stem cells and progenitors, shows increased chromatin accessibility in human acute myeloid leukaemia stem cells compared to blasts. This difference correlates with MYC expression and patient outcome. We propose that clusters of enhancers, such as BENC, form highly combinatorial systems that allow precise control of gene expression across normal cellular hierarchies and which also can be hijacked in malignancies.

Two independent modes of chromatin organization revealed by cohesin removal.

Imaging and chromosome conformation capture studies have revealed several layers of chromosome organization, including segregation into megabase-sized active and inactive compartments, and partitioning into sub-megabase domains (TADs). It remains unclear, however, how these layers of organization form, interact with one another and influence genome function. Here we show that deletion of the cohesin-loading factor Nipbl in mouse liver leads to a marked reorganization of chromosomal folding. TADs and associated Hi-C peaks vanish globally, even in the absence of transcriptional changes. By contrast, compartmental segregation is preserved and even reinforced. Strikingly, the disappearance of TADs unmasks a finer compartment structure that accurately reflects the underlying epigenetic landscape. These observations demonstrate that the three-dimensional organization of the genome results from the interplay of two independent mechanisms: cohesin-independent segregation of the genome into fine-scale compartments, defined by chromatin state; and cohesin-dependent formation of TADs, possibly by loop extrusion, which helps to guide distant enhancers to their target genes.

The Shh Topological Domain Facilitates the Action of Remote Enhancers by Reducing the Effects of Genomic Distances.

Gene expression often requires interaction between promoters and distant enhancers, which occur within the context of highly organized topologically associating domains (TADs). Using a series of engineered chromosomal rearrangements at the Shh locus, we carried out an extensive fine-scale characterization of the factors that govern the long-range regulatory interactions controlling Shh expression. We show that Shh enhancers act pervasively, yet not uniformly, throughout the TAD. Importantly, changing intra-TAD distances had no impact on Shh expression. In contrast, inversions disrupting the TAD altered global folding of the region and prevented regulatory contacts in a distance-dependent manner. Our data indicate that the Shh TAD promotes distance-independent contacts between distant regions that would otherwise interact only sporadically, enabling functional communication between them. In large genomes where genomic distances per se can limit regulatory interactions, this function of TADs could be as essential for gene expression as the formation of insulated neighborhoods.

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.

Gene regulation at a distance: From remote enhancers to 3D regulatory ensembles.

Large-scale identification of elements associated with gene expression revealed that many of them are located extremely far from gene transcriptional start sites. We review here the growing evidence that show that distal cis-acting elements provide key instructions to genes, as genetic variation affecting them is growingly identified as an importance source of phenotypic diversity and disease. We discuss the different mechanisms that allow these elements to exert their regulatory functions, in a robust and specific manner, despite the large genomic distances separating them from their target genes. We particularly focus on the role of the structural organization of the genome in guiding such regulatory interactions.

Cis-regulatory architecture of a brain signaling center predates the origin of chordates.

Genomic approaches have predicted hundreds of thousands of tissue-specific cis-regulatory sequences, but the determinants critical to their function and evolutionary history are mostly unknown. Here we systematically decode a set of brain enhancers active in the zona limitans intrathalamica (zli), a signaling center essential for vertebrate forebrain development via the secreted morphogen Sonic hedgehog (Shh). We apply a de novo motif analysis tool to identify six position-independent sequence motifs together with their cognate transcription factors that are essential for zli enhancer activity and Shh expression in the mouse embryo. Using knowledge of this regulatory lexicon, we discover new Shh zli enhancers in mice and a functionally equivalent element in hemichordates, indicating an ancient origin of the Shh zli regulatory network that predates the chordate phylum. These findings support a strategy for delineating functionally conserved enhancers in the absence of overt sequence homologies and over extensive evolutionary distances.

Gene regulation during development in the light of topologically associating domains.

During embryonic development, complex transcriptional programs govern the precision of gene expression. Many key developmental genes are regulated via cis-regulatory elements that are located far away in the linear genome. How sequences located hundreds of kilobases away from a promoter can influence its activity has been the subject of numerous speculations, which all underline the importance of the 3D-organization of the genome. The recent advent of chromosome conformation capture techniques has put into focus the subdivision of the genome into topologically associating domains (TADs). TADs may influence regulatory activities on multiple levels. The relative invariance of TAD limits across cell types suggests that they may form fixed structural domains that could facilitate and/or confine long-range regulatory interactions. However, most recent studies suggest that interactions within TADs are more variable and dynamic than initially described. Hence, different models are emerging regarding how TADs shape the complex 3D conformations, and thereafter influence the networks of cis-interactions that govern gene expression during development. For further resources related to this article, please visit the WIREs website.

Transcriptome profiling of white adipose tissue in a mouse model for 15q duplication syndrome.

Obesity is not only associated with unhealthy lifestyles, but also linked to genetic predisposition. Previously, we generated an autism mouse model (patDp/+) that carries a 6.3 Mb paternal duplication homologous to the human 15q11-q13 locus. Chromosomal abnormalities in this region are known to cause autism spectrum disorder, Prader-Willi syndrome, and Angelman syndrome in humans. We found that, in addition to autistic-like behaviors, patDp/+ mice display late-onset obesity and hypersensitivity to a high-fat diet. These phenotypes are likely to be the results of genetic perturbations since the energy expenditures and food intakes of patDp/+ mice do not significantly differ from those of wild-type mice. Intriguingly, we found that an enlargement of adipose cells precedes the onset of obesity in patDp/+ mice. To understand the underlying molecular networks responsible for this pre-obese phenotype, we performed transcriptome profiling of white adipose tissue from patDp/+ and wild-type mice using microarray. We identified 230 genes as differentially expressed genes. Sfrp5 - a gene whose expression is positively correlated with adipocyte size, was found to be up-regulated, and Fndc5, a potent inducer of brown adipogenesis was identified to be the top down-regulated gene. Subsequent pathway analysis highlighted a set of 35 molecules involved in energy production, lipid metabolism, and small molecule biochemistry as the top candidate biological network responsible for the pre-obese phenotype of patDp/+. The microarray data were deposited in NCBI Gene Expression Omnibus database with accession number GSE58191. Ultimately, our dataset provides novel insights into the molecular mechanism of obesity and demonstrated that patDp/+ is a valuable mouse model for obesity research.

Model mice for 15q11-13 duplication syndrome exhibit late-onset obesity and altered lipid metabolism.

Copy number variations on human chromosome 15q11-q13 have been implicated in several neurodevelopmental disorders. A paternal loss or duplication of the Prader-Willi syndrome/Angelman syndrome (PWS/AS) region confers a risk of obesity, although the mechanism remains a mystery due to a lack of an animal model that accurately recreates the obesity phenotype. We performed detailed analyses of mice with duplication of PWS/AS locus (6 Mb) generated by chromosome engineering and found that animals with a paternal duplication of this region (patDp/+) show late-onset obesity, high sensitivity for high-fat diet, high levels of blood leptin and insulin without an increase in food intake. We show that prior to becoming obese, young patDp/+ mice already had enlarged white adipocytes. Transcriptome analysis of adipose tissue revealed an up-regulation of Secreted frizzled-related protein 5 (Sfrp5), known to promote adipogenesis. We additionally generated a new mouse model of paternal duplication focusing on a 3 Mb region (3 Mb patDp/+) within the PWS/AS locus. These mice recapitulate the obese phenotypes including expansion of visceral adipose tissue. Our results suggest paternally expressed genes in PWS/AS locus play a significant role for the obesity and identify new potential targets for future research and treatment of obesity.