Opposing functions of the ETS factor family define Shh spatial expression in limb buds and underlie polydactyly

Topology of mammalian developmental enhancers and their regulatory landscapes

How a complex animal can arise from a fertilized egg is one of the oldest and most fascinating questions of biology, the answer to which is encoded in the genome. Body shape and organ development, and their integration into a functional organism all depend on the precise expression of genes in space and time. The orchestration of transcription relies mostly on surrounding control sequences such as enhancers, millions of which form complex regulatory landscapes in the non-coding genome. Recent research shows that high-order chromosome structures make an important contribution to enhancer functionality by triggering their physical interactions with target genes.

Massively parallel decoding of mammalian regulatory sequences supports a flexible organizational model

Despite continual progress in the cataloging of vertebrate regulatory elements, little is known about their organization and regulatory architecture. Here we describe a massively parallel experiment to systematically test the impact of copy number, spacing, combination and order of transcription factor binding sites on gene expression. A complex library of ∼5,000 synthetic regulatory elements containing patterns from 12 liver-specific transcription factor binding sites was assayed in mice and in HepG2 cells. We find that certain transcription factors act as direct drivers of gene expression in homotypic clusters of binding sites, independent of spacing between sites, whereas others function only synergistically. Heterotypic enhancers are stronger than their homotypic analogs and favor specific transcription factor binding site combinations, mimicking putative native enhancers. Exhaustive testing of binding site permutations suggests that there is flexibility in binding site order. Our findings provide quantitative support for a flexible model of regulatory element activity and suggest a framework for the design of synthetic tissue-specific enhancers.

SP8 regulates signaling centers during craniofacial development

Much of the bone, cartilage and smooth muscle of the vertebrate face is derived from neural crest (NC) cells. During craniofacial development, the anterior neural ridge (ANR) and olfactory pit (OP) signaling centers are responsible for driving the outgrowth, survival, and differentiation of NC populated facial prominences, primarily via FGF. While much is known about the functional importance of signaling centers, relatively little is understood of how these signaling centers are made and maintained. In this report we describe a dramatic craniofacial malformation in mice mutant for the zinc finger transcription factor gene Sp8. At E14.5 they show facial prominences that are reduced in size and underdeveloped, giving an almost faceless phenotype. At later times they show severe midline defects, excencephaly, hyperterlorism, cleft palate, and a striking loss of many NC and paraxial mesoderm derived cranial bones. Sp8 expression was primarily restricted to the ANR and OP regions during craniofacial development. Analysis of an extensive series of conditional Sp8 mutants confirmed the critical role of Sp8 in signaling centers, and not directly in the NC and paraxial mesoderm cells. The NC cells of the Sp8 mutants showed increased levels of apoptosis and decreased cell proliferation, thereby explaining the reduced sizes of the facial prominences. Perturbed gene expression in the Sp8 mutants was examined by laser capture microdissection coupled with microarrays, as well as in situ hybridization and immunostaining. The most dramatic differences included striking reductions in Fgf8 and Fgf17 expression in the ANR and OP signaling centers. We were also able to achieve genetic and pharmaceutical partial rescue of the Sp8 mutant phenotype by reducing Sonic Hedgehog (SHH) signaling. These results show that Sp8 primarily functions to promote Fgf expression in the ANR and OP signaling centers that drive the survival, proliferation, and differentiation of the NC and paraxial mesoderm that make the face.

The spatial organization of the human genome

In vivo, the human genome functions as a complex, folded, three-dimensional chromatin polymer. Understanding how the human genome is spatially organized and folded inside the cell nucleus is therefore central to understanding how genes are regulated in normal development and dysregulated in disease. Established light microscopy-based approaches and more recent molecular chromosome conformation capture methods are now combining to give us unprecedented insight into this fascinating aspect of human genomics.

Functional anatomy of distant-acting mammalian enhancers

Transcriptional enhancers are a major class of functional element embedded in the vast non-coding portion of the human genome. Acting over large genomic distances, enhancers play critical roles in the tissue and cell type-specific regulation of genes, and there is mounting evidence that they contribute to the aetiology of many human diseases. Methods for genome-wide mapping of enhancer regions are now available, but the functional architecture contained within human enhancer elements remains unclear. Here, we review recent approaches aimed at understanding the functional anatomy of individual enhancer elements, using systematic qualitative and quantitative assessments of mammalian enhancer variants in cultured cells and in vivo. These studies provide direct insight into common architectural characteristics of enhancers including the presence of multiple transcription factor-binding sites and the mixture of both transcriptionally activating and repressing domains within the same enhancer. Despite such progress in understanding the functional composition of enhancers, the inherent complexities of enhancer anatomy continue to limit our ability to predict the impact of sequence changes on in vivo enhancer function. While providing an initial glimpse into the mutability of mammalian enhancers, these observations highlight the continued need for experimental enhancer assessment as genome sequencing becomes routine in the clinic.

Alterations to the remote control of Shh gene expression cause congenital abnormalities

Multi-species conserved non-coding elements occur in the vertebrate genome and are clustered in the vicinity of developmentally regulated genes. Many are known to act as cis-regulators of transcription and may reside at long distances from the genes they regulate. However, the relationship of conserved sequence to encoded regulatory information and indeed, the mechanism by which these contribute to long-range transcriptional regulation is not well understood. The ZRS, a highly conserved cis-regulator, is a paradigm for such long-range gene regulation. The ZRS acts over approximately 1 Mb to control spatio-temporal expression of Shh in the limb bud and mutations within it result in a number of limb abnormalities, including polydactyly, tibial hypoplasia and syndactyly. We describe the activity of this developmental regulator and discuss a number of mechanisms by which regulatory mutations in this enhancer function to cause congenital abnormalities.

Structural variations, the regulatory landscape of the genome and their alteration in human disease

High-throughput genomic technologies are revolutionizing human genetics. So far the focus has been on the 1.5% of the genome, which is coding, in spite of the fact that the great majority of genomic variants fall outside the coding regions. Recent efforts to annotate the non-coding sequence show that over 80% of the genome is biochemically active. The genome is divided into regulatory domains consisting of sequence regions that enhance and/or silence the expression of nearby genes and are, in some cases, separated by boundaries with insulator activity. In this paper, we review the recent advances in the identification of variations that influence gene regulation and their consequences for human disease. We hypothesize that structural variations outside of the coding genome can interfere with normal gene regulation by disrupting the regulatory landscape. Therefore, the regulatory landscape of the genome has also to be taken into consideration when investigating the pathology of human disease.

Phenotypic impact of genomic structural variation: insights from and for human disease

Genomic structural variants have long been implicated in phenotypic diversity and human disease, but dissecting the mechanisms by which they exert their functional impact has proven elusive. Recently however, developments in high-throughput DNA sequencing and chromosomal engineering technology have facilitated the analysis of structural variants in human populations and model systems in unprecedented detail. In this Review, we describe how structural variants can affect molecular and cellular processes, leading to complex organismal phenotypes, including human disease. We further present advances in delineating disease-causing elements that are affected by structural variants, and we discuss future directions for research on the functional consequences of structural variants.

Complex effects of nucleotide variants in a mammalian cis-regulatory element

Cis-regulatory elements (CREs) control gene expression by recruiting transcription factors (TFs) and other DNA binding proteins. We aim to understand how individual nucleotides contribute to the function of CREs. Here we introduce CRE analysis by sequencing (CRE-seq), a high-throughput method for producing and testing large numbers of reporter genes in mammalian cells. We used CRE-seq to assay >1,000 single and double nucleotide mutations in a 52-bp CRE in the Rhodopsin promoter that drives strong and specific expression in mammalian photoreceptors. We find that this particular CRE is remarkably complex. The majority (86%) of single nucleotide substitutions in this sequence exert significant effects on regulatory activity. Although changes in the affinity of known TF binding sites explain some of these expression changes, we present evidence for complex phenomena, including binding site turnover and TF competition. Analysis of double mutants revealed complex, nucleotide-specific interactions between residues in different TF binding sites. We conclude that some mammalian CREs are finely tuned by evolution and function through complex, nonadditive interactions between bound TFs. CRE-seq will be an important tool to uncover the rules that govern these interactions.

Driving developmental and evolutionary change: A systems biology view

Embryonic development is underpinned by ∼50 core processes that drive morphogenesis, growth, patterning and differentiation, and each is the functional output of a complex molecular network. Processes are thus the natural and parsimonious link between genotype and phenotype and the obvious focus for any discussion of biological change. Here, the implications of this approach are explored. One is that many features of developmental change can be modeled as mathematical graphs, or sets of connected triplets of the general form <noun><verb><noun>. In these, the verbs (edges) are the outputs of the processes that drive change and the nouns (nodes) are the time-dependent states of biological entities (from molecules to tissues). Such graphs help unpick the multi-level complexity of developmental phenomena and may help suggest new experiments. Another comes from analyzing the effect of mutation that lead to tinkering with the dynamic properties of these processes and to congenital abnormalities; if these changes are both inherited and advantageous, they become evolutionary modifications. In this context, protein networks often represents what classical evolutionary genetics sees as genes, and the realization that traits reflect the output processes of complex networks, particularly for growth, patterning and pigmentation, rather than anything simpler clarifies some problems that the evolutionary synthesis of the 1950s has found hard to solve. In the wider context, most processes are used many times in development and cooperate to produce tissue modules (bones, branching duct systems, muscles etc.). Their underlying generative networks can thus be thought of as genomic modules or subroutines.

Anterior-posterior differences in HoxD chromatin topology in limb development

A late phase of HoxD activation is crucial for the patterning and growth of distal structures across the anterior-posterior (A-P) limb axis of mammals. Polycomb complexes and chromatin compaction have been shown to regulate Hox loci along the main body axis in embryonic development, but the extent to which they have a role in limb-specific HoxD expression, an evolutionary adaptation defined by the activity of distal enhancer elements that drive expression of 5' Hoxd genes, has yet to be fully elucidated. We reveal two levels of chromatin topology that differentiate distal limb A-P HoxD activity. Using both immortalised cell lines derived from posterior and anterior regions of distal E10.5 mouse limb buds, and analysis in E10.5 dissected limb buds themselves, we show that there is a loss of polycomb-catalysed H3K27me3 histone modification and a chromatin decompaction over HoxD in the distal posterior limb compared with anterior. Moreover, we show that the global control region (GCR) long-range enhancer spatially colocalises with the 5' HoxD genomic region specifically in the distal posterior limb. This is consistent with the formation of a chromatin loop between 5' HoxD and the GCR regulatory module at the time and place of distal limb bud development when the GCR participates in initiating Hoxd gene quantitative collinearity and Hoxd13 expression. This is the first example of A-P differences in chromatin compaction and chromatin looping in the development of the mammalian secondary body axis (limb).

Conserved cis-regulatory regions in a large genomic landscape control SHH and BMP-regulated Gremlin1 expression in mouse limb buds

Background

Mouse limb bud is a prime model to study the regulatory interactions that control vertebrate organogenesis. Major aspects of limb bud development are controlled by feedback loops that define a self-regulatory signalling system. The SHH/GREM1/AER-FGF feedback loop forms the core of this signalling system that operates between the posterior mesenchymal organiser and the ectodermal signalling centre. The BMP antagonist Gremlin1 (GREM1) is a critical node in this system, whose dynamic expression is controlled by BMP, SHH, and FGF signalling and key to normal progression of limb bud development. Previous analysis identified a distant cis-regulatory landscape within the neighbouring Formin1 (Fmn1) locus that is required for Grem1 expression, reminiscent of the genomic landscapes controlling HoxD and Shh expression in limb buds.

Results

Three highly conserved regions (HMCO1-3) were identified within the previously defined critical genomic region and tested for their ability to regulate Grem1 expression in mouse limb buds. Using a combination of BAC and conventional transgenic approaches, a 9 kb region located ~70 kb downstream of the Grem1 transcription unit was identified. This region, termed Grem1 Regulatory Sequence 1 (GRS1), is able to recapitulate major aspects of Grem1 expression, as it drives expression of a LacZ reporter into the posterior and, to a lesser extent, in the distal-anterior mesenchyme. Crossing the GRS1 transgene into embryos with alterations in the SHH and BMP pathways established that GRS1 depends on SHH and is modulated by BMP signalling, i.e. integrates inputs from these pathways. Chromatin immunoprecipitation revealed interaction of endogenous GLI3 proteins with the core cis-regulatory elements in the GRS1 region. As GLI3 is a mediator of SHH signal transduction, these results indicated that SHH directly controls Grem1 expression through the GRS1 region. Finally, all cis-regulatory regions within the Grem1 genomic landscape locate to the DNAse I hypersensitive sites identified in this genomic region by the ENCODE consortium.

Conclusions

This study establishes that distant cis-regulatory regions scattered through a larger genomic landscape control the highly dynamic expression of Grem1, which is key to normal progression of mouse limb bud development.

Transcription factors: from enhancer binding to developmental control

Developmental progression is driven by specific spatiotemporal domains of gene expression, which give rise to stereotypically patterned embryos even in the presence of environmental and genetic variation. Views of how transcription factors regulate gene expression are changing owing to recent genome-wide studies of transcription factor binding and RNA expression. Such studies reveal patterns that, at first glance, seem to contrast with the robustness of the developmental processes they encode. Here, we review our current knowledge of transcription factor function from genomic and genetic studies and discuss how different strategies, including extensive cooperative regulation (both direct and indirect), progressive priming of regulatory elements, and the integration of activities from multiple enhancers, confer specificity and robustness to transcriptional regulation during development.

Integration of the transcriptional networks regulating limb morphogenesis

The developing limb is one of the best described vertebrate systems for understanding how coordinated gene expression during embryogenesis leads to the structures present in the mature organism. This knowledge, derived from decades of research, is largely based upon gain- and loss-of-function experiments. These studies have provided limited information about how the key signaling pathways interact with each other and the downstream effectors of these pathways. We summarize our current understanding of known genetic interactions in the context of three temporally defined gene regulatory networks. These networks crystallize our current knowledge, depicting a dynamic process involving multiple feedback loops between the ectoderm and mesoderm. At the same time, they highlight the fact that many essential processes are still largely undescribed. Much of the dynamic transcriptional activity occurring during development is regulated by distal cis-regulatory elements. Modern genomic tools have provided new approaches for studying the function of cis-regulatory elements and we discuss the results of these studies in regard to understanding limb development. Ultimately, these genomic techniques will allow scientists to understand how multiple signaling pathways are integrated in space and time to drive gene expression and regulate the formation of the limb.

The Sonic hedgehog gradient in the developing limb

A gradient of Sonic hedgehog (Shh) plays a major role in specifying the antero-posterior pattern of structures that develop in the distal part of the vertebrate limb, in particular, the antero-posterior pattern of the digits. Classical embryological experiments identified the polarizing region (or zone of polarizing activity, ZPA), a signaling region at the posterior margin of the early chick wing bud and, consistent with a model in which production of a diffusible morphogen specifies antero-posterior positional information, polarizing region signaling was shown to be dose dependent and long range. It is now well established that the vertebrate hedgehog gene, Sonic hedgehog (Shh), which encodes a secreted protein, is expressed in the polarizing region of the chick wing and that Shh signaling has the same characteristics as polarizing region signaling. Shh expression at the posterior of the early limb bud and the mechanism of Shh signal transduction are conserved among vertebrates including mammals. However, it is unlikely that a simple Shh gradient is responsible for digit pattern formation in mammalian limbs and there is still little understanding of how positional information specified by Shh signaling is encoded and translated into digit anatomy. Alterations in Shh signaling underlie some congenital limb abnormalities and also changes in timing and extent of Shh signaling appear to be related to the evolution of morphological diversity of vertebrate limbs.