Shuffling of cis-regulatory elements is a pervasive feature of the vertebrate lineage

Developmental and evolutionary comparative analysis of a regulatory landscape in mouse and chicken

Modifications in gene regulation are driving forces in the evolution of organisms. Part of these changes involve cis-regulatory elements (CREs), which contact their target genes through higher-order chromatin structures. However, how such architectures and variations in CREs contribute to transcriptional evolvability remains elusive. We use Hoxd genes as a paradigm for the emergence of regulatory innovations, as many relevant enhancers are located in a regulatory landscape highly conserved in amniotes. Here, we analysed their regulation in murine vibrissae and chicken feather primordia, two skin appendages expressing different Hoxd gene subsets, and compared the regulation of these genes in these appendages with that in the elongation of the posterior trunk. In the two former structures, distinct subsets of Hoxd genes are contacted by different lineage-specific enhancers, probably as a result of using an ancestral chromatin topology as an evolutionary playground, whereas the gene regulation that occurs in the mouse and chicken embryonic trunk partially relies on conserved CREs. A high proportion of these non-coding sequences active in the trunk have functionally diverged between species, suggesting that transcriptional robustness is maintained, despite considerable divergence in enhancer sequences.

Summary: Analyses of the relationships between chromatin architecture and regulatory activities at the HoxD locus show that ancestral transcription patterns can be maintained while new regulations evolve.

Conserved noncoding sequences conserve biological networks and influence genome evolution

Comparative genomics approaches have identified numerous conserved cis-regulatory sequences near genes in plant genomes. Despite the identification of these conserved noncoding sequences (CNSs), our knowledge of their functional importance and selection remains limited. Here, we used a combination of DNA methylome analysis, microarray expression analyses, and functional annotation to study these sequences in the model tree Populus trichocarpa. Methylation in CG contexts and non-CG contexts was lower in CNSs, particularly CNSs in the 5'-upstream regions of genes, compared with other sites in the genome. We observed that CNSs are enriched in genes with transcription and binding functions, and this also associated with syntenic genes and those from whole-genome duplications, suggesting that cis-regulatory sequences play a key role in genome evolution. We detected a significant positive correlation between CNS number and protein interactions, suggesting that CNSs may have roles in the evolution and maintenance of biological networks. The divergence of CNSs indicates that duplication-degeneration-complementation drives the subfunctionalization of a proportion of duplicated genes from whole-genome duplication. Furthermore, population genomics confirmed that most CNSs are under strong purifying selection and only a small subset of CNSs shows evidence of adaptive evolution. These findings provide a foundation for future studies exploring these key genomic features in the maintenance of biological networks, local adaptation, and transcription.

Conserved non-coding elements: developmental gene regulation meets genome organization

Comparative genomics has revealed a class of non-protein-coding genomic sequences that display an extraordinary degree of conservation between two or more organisms, regularly exceeding that found within protein-coding exons. These elements, collectively referred to as conserved non-coding elements (CNEs), are non-randomly distributed across chromosomes and tend to cluster in the vicinity of genes with regulatory roles in multicellular development and differentiation. CNEs are organized into functional ensembles called genomic regulatory blocks–dense clusters of elements that collectively coordinate the expression of shared target genes, and whose span in many cases coincides with topologically associated domains. CNEs display sequence properties that set them apart from other sequences under constraint, and have recently been proposed as useful markers for the reconstruction of the evolutionary history of organisms. Disruption of several of these elements is known to contribute to diseases linked with development, and cancer. The emergence, evolutionary dynamics and functions of CNEs still remain poorly understood, and new approaches are required to enable comprehensive CNE identification and characterization. Here, we review current knowledge and identify challenges that need to be tackled to resolve the impasse in understanding extreme non-coding conservation.

Hedgehog signalling acts upstream of Laminin alpha1 transcription in the zebrafish paraxial mesoderm

Laminin-111 (α1β1γ1) is a member of the Laminin family of extra-cellular matrix proteins that comprises 16 members, components of basement membranes. Laminin-111, one of the first Laminin proteins synthesised during embryogenesis, is required for basement membrane deposition and has essential roles in tissue morphogenesis and patterning. Yet, the mechanisms controlling Laminin-111 expression are poorly understood. We generated a zebrafish transgenic reporter line that reproduces faithfully the expression pattern of lama1, the gene encoding Laminin α1, and we used this reporter line to investigate lama1 transcriptional regulation. Our findings established that lama1 expression is controlled by intronic enhancers, including an enhancer directing expression in the paraxial mesoderm, anterior spinal cord and hindbrain, located in intron 1. We show that Hedgehog signalling is necessary and sufficient for lama1 transcription in the paraxial mesoderm and identify putative Gli/Zic binding sites that may mediate this control. These findings uncover a conserved role for Hedgehog signalling in the control of basement membrane assembly via its transcriptional regulation of lama1, and provide a mechanism to coordinate muscle cell fate specification in the zebrafish embryo.

C2H2 zinc finger proteins greatly expand the human regulatory lexicon

Cys2-His2 zinc finger (C2H2-ZF) proteins represent the largest class of putative human transcription factors. However, for most C2H2-ZF proteins it is unknown whether they even bind DNA or, if they do, to which sequences. Here, by combining data from a modified bacterial one-hybrid system with protein-binding microarray and chromatin immunoprecipitation analyses, we show that natural C2H2-ZFs encoded in the human genome bind DNA both in vitro and in vivo, and we infer the DNA recognition code using DNA-binding data for thousands of natural C2H2-ZF domains. In vivo binding data are generally consistent with our recognition code and indicate that C2H2-ZF proteins recognize more motifs than all other human transcription factors combined. We provide direct evidence that most KRAB-containing C2H2-ZF proteins bind specific endogenous retroelements (EREs), ranging from currently active to ancient families. The majority of C2H2-ZF proteins, including KRAB proteins, also show widespread binding to regulatory regions, indicating that the human genome contains an extensive and largely unstudied adaptive C2H2-ZF regulatory network that targets a diverse range of genes and pathways.

Conserved noncoding elements follow power-law-like distributions in several genomes as a result of genome dynamics

Conserved, ultraconserved and other classes of constrained elements (collectively referred as CNEs here), identified by comparative genomics in a wide variety of genomes, are non-randomly distributed across chromosomes. These elements are defined using various degrees of conservation between organisms and several thresholds of minimal length. We here investigate the chromosomal distribution of CNEs by studying the statistical properties of distances between consecutive CNEs. We find widespread power-law-like distributions, i.e. linearity in double logarithmic scale, in the inter-CNE distances, a feature which is connected with fractality and self-similarity. Given that CNEs are often found to be spatially associated with genes, especially with those that regulate developmental processes, we verify by appropriate gene masking that a power-law-like pattern emerges irrespectively of whether elements found close or inside genes are excluded or not. An evolutionary model is put forward for the understanding of these findings that includes segmental or whole genome duplication events and eliminations (loss) of most of the duplicated CNEs. Simulations reproduce the main features of the observed size distributions. Power-law-like patterns in the genomic distributions of CNEs are in accordance with current knowledge about their evolutionary history in several genomes.

Deep conservation of cis-regulatory elements in metazoans

Despite the vast morphological variation observed across phyla, animals share multiple basic developmental processes orchestrated by a common ancestral gene toolkit. These genes interact with each other building complex gene regulatory networks (GRNs), which are encoded in the genome by cis-regulatory elements (CREs) that serve as computational units of the network. Although GRN subcircuits involved in ancient developmental processes are expected to be at least partially conserved, identification of CREs that are conserved across phyla has remained elusive. Here, we review recent studies that revealed such deeply conserved CREs do exist, discuss the difficulties associated with their identification and describe new approaches that will facilitate this search.

Gene transcription in the zebrafish embryo: regulators and networks

The precise spatial and temporal control of gene expression is a key process in the development, maintenance and regeneration of the vertebrate body. A substantial proportion of vertebrate genomes encode genes that control the transcription of the genetic information into mRNA. The zebrafish is particularly well suited to investigate gene regulatory networks underlying the control of gene expression during development due to the external development of its transparent embryos and the increasingly sophisticated tools for genetic manipulation available for this model system. We review here recent data on the analysis of cis-regulatory modules, transcriptional regulators and their integration into gene regulatory networks in the zebrafish, using the developing spinal cord as example.

Highly conserved elements discovered in vertebrates are present in non-syntenic loci of tunicates, act as enhancers and can be transcribed during development

Co-option of cis-regulatory modules has been suggested as a mechanism for the evolution of expression sites during development. However, the extent and mechanisms involved in mobilization of cis-regulatory modules remains elusive. To trace the history of non-coding elements, which may represent candidate ancestral cis-regulatory modules affirmed during chordate evolution, we have searched for conserved elements in tunicate and vertebrate (Olfactores) genomes. We identified, for the first time, 183 non-coding sequences that are highly conserved between the two groups. Our results show that all but one element are conserved in non-syntenic regions between vertebrate and tunicate genomes, while being syntenic among vertebrates. Nevertheless, in all the groups, they are significantly associated with transcription factors showing specific functions fundamental to animal development, such as multicellular organism development and sequence-specific DNA binding. The majority of these regions map onto ultraconserved elements and we demonstrate that they can act as functional enhancers within the organism of origin, as well as in cross-transgenesis experiments, and that they are transcribed in extant species of Olfactores. We refer to the elements as 'Olfactores conserved non-coding elements'.

Transcriptional enhancers in protein-coding exons of vertebrate developmental genes

Many conserved noncoding sequences function as transcriptional enhancers that regulate gene expression. Here, we report that protein-coding DNA also frequently contains enhancers functioning at the transcriptional level. We tested the enhancer activity of 31 protein-coding exons, which we chose based on strong sequence conservation between zebrafish and human, and occurrence in developmental genes, using a Tol2 transposable GFP reporter assay in zebrafish. For each exon we measured GFP expression in hundreds of embryos in 10 anatomies via a novel system that implements the voice-recognition capabilities of a cellular phone. We find that 24/31 (77%) exons drive GFP expression compared to a minimal promoter control, and 14/24 are anatomy-specific (expression in four anatomies or less). GFP expression driven by these coding enhancers frequently overlaps the anatomies where the host gene is expressed (60%), suggesting self-regulation. Highly conserved coding sequences and highly conserved noncoding sequences do not significantly differ in enhancer activity (coding: 24/31 vs. noncoding: 105/147) or tissue-specificity (coding: 14/24 vs. noncoding: 50/105). Furthermore, coding and noncoding enhancers display similar levels of the enhancer-related histone modification H3K4me1 (coding: 9/24 vs noncoding: 34/81). Meanwhile, coding enhancers are over three times as likely to contain an H3K4me1 mark as other exons of the host gene. Our work suggests that developmental transcriptional enhancers do not discriminate between coding and noncoding DNA and reveals widespread dual functions in protein-coding DNA.

Rapid turnover of functional sequence in human and other genomes

The amount of a genome's sequence that is functional has been surprisingly difficult to estimate accurately. This has severely hindered analyses asking whether the amount of functional genomic sequence correlates with organismal complexity. Most studies estimate these amounts by considering nucleotide substitution rates within aligned sequences. These approaches show reduced power to identify sequence that is aligned, functional, and constrained only within narrowly defined phyla. The neutral indel model exploits insertions or deletions (indels) rather than substitutions in predicting functional sequence. Surprisingly, this method indicates that half of all functional sequence is specific to individual eutherian lineages. This review considers the rates at which coding or noncoding and functional or nonfunctional sequence changes among mammalian genomes. In contrast to the slow rate at which protein-coding sequence changes, functional noncoding sequence appears to change or be turned over at rapid rates in mammals.

Fishing for function: zebrafish BAC transgenics for functional genomics

Transgenics using bacterial artificial chromosomes (BACs) offers a great opportunity to look at gene regulation in a developing embryo. The modified BAC containing a reporter inserted just before the translational start site of the gene of interest allows for the visualization of spatio-temporal gene expression. Though this method has been used in the mouse model extensively, its utility in zebrafish studies is relatively new. This review aims to look at the utility of making BAC transgenics in zebrafish and its applications in functional genomics. We look at the various methods to modify the BAC, some limitations and what the future holds.

A catalogue of eukaryotic transcription factor types, their evolutionary origin, and species distribution

Transcription factors (TFs) play key roles in the regulation of gene expression by binding in a sequence-specific manner to genomic DNA. In eukaryotes, DNA binding is achieved by a wide range of structural forms and motifs. TFs are typically classified by their DNA-binding domain (DBD) type. In this chapter, we catalogue and survey 91 different TF DBD types in metazoa, plants, fungi, and protists. We briefly discuss well-characterized TF families representing the major DBD superclasses. We also examine the species distributions and inferred evolutionary histories of the various families, and the potential roles played by TF family expansion and dimerization.

snail gene expression in the medaka, Oryzias latipes

Snail transcription factors have prominent roles during embryonic development of vertebrates. They are often involved in cell migration processes during neural crest development, epithelial-mesenchymal transition and cancer progression. Comparative expression studies of snai gene family members in different vertebrate species are expected to contribute to a better understanding of their roles during development and reflect their evolutionary history. To investigate and to compare the expression patterns of snai genes in a second main fish model we used the medaka fish (Oryzias latipes), a complementary teleost model to zebrafish. We identified three snai gene family members, snai1a, snai1b and snai2. Phylogenetic and synteny analysis show a close relatedness of all family members to other vertebrate snai genes. Surprisingly, no homologue of snai3 could be identified in medaka, although this gene is present in zebrafish and the puffer fishes. Here we demonstrate that while most expression domains of medaka snai genes are comparable to zebrafish, the contribution of the individual paralogs to the overall pattern differs between the two teleosts and indicate lineage specific expression shuffling.

Conserved expression without conserved regulatory sequence: the more things change, the more they stay the same

Regulatory regions with similar transcriptional output often have little overt sequence similarity, both within and between genomes. Although cis- and trans-regulatory changes can contribute to sequence divergence without dramatically altering gene expression outputs, heterologous DNA often functions similarly in organisms that share little regulatory sequence similarities (e.g. human DNA in fish), indicating that trans-regulatory mechanisms tend to diverge more slowly and can accommodate a variety of cis-regulatory configurations. This capacity to 'tinker' with regulatory DNA probably relates to the complexity, robustness and evolvability of regulatory systems, but cause-and-effect relationships among evolutionary processes and properties of regulatory systems remain a topic of debate. The challenge of understanding the concrete mechanisms underlying cis-regulatory evolution - including the conservation of function without the conservation of sequence - relates to the challenge of understanding the function of regulatory systems in general. Currently, we are largely unable to recognize functionally similar regulatory DNA.

Deeply conserved chordate noncoding sequences preserve genome synteny but do not drive gene duplicate retention

Animal genomes possess highly conserved cis-regulatory sequences that are often found near genes that regulate transcription and development. Researchers have proposed that the strong conservation of these sequences may affect the evolution of the surrounding genome, both by repressing rearrangement, and possibly by promoting duplicate gene retention. Conflicting data, however, have made the validity of these propositions unclear. Here, we use a new computational method to identify phylogenetically conserved noncoding elements (PCNEs) in a manner that is not biased by rearrangement and duplication. This method is powerful enough to identify more than a thousand PCNEs that have been conserved between vertebrates and the basal chordate amphioxus. We test 42 of our PCNEs in transgenic zebrafish assays--including examples from vertebrates and amphioxus--and find that the majority are functional enhancers. We find that PCNEs are enriched around genes with ancient synteny conservation, and that this association is strongest for extragenic PCNEs, suggesting that cis-regulatory interdigitation plays a key role in repressing genome rearrangement. Next, we classify mouse and zebrafish genes according to association with PCNEs, synteny conservation, duplication history, and presence in bidirectional promoter pairs, and use these data to cluster gene functions into a series of distinct evolutionary patterns. These results demonstrate that subfunctionalization of conserved cis-regulation has not been the primary determinate of gene duplicate retention in vertebrates. Instead, the data support the gene balance hypothesis, which proposes that duplicate retention has been driven by selection against dosage imbalances in genes with many protein connections.

Conservation of core gene expression in vertebrate tissues

BACKGROUND

Vertebrates share the same general body plan and organs, possess related sets of genes, and rely on similar physiological mechanisms, yet show great diversity in morphology, habitat and behavior. Alteration of gene regulation is thought to be a major mechanism in phenotypic variation and evolution, but relatively little is known about the broad patterns of conservation in gene expression in non-mammalian vertebrates.

RESULTS

We measured expression of all known and predicted genes across twenty tissues in chicken, frog and pufferfish. By combining the results with human and mouse data and considering only ten common tissues, we have found evidence of conserved expression for more than a third of unique orthologous genes. We find that, on average, transcription factor gene expression is neither more nor less conserved than that of other genes. Strikingly, conservation of expression correlates poorly with the amount of conserved nonexonic sequence, even using a sequence alignment technique that accounts for non-collinearity in conserved elements. Many genes show conserved human/fish expression despite having almost no nonexonic conserved primary sequence.

CONCLUSIONS

There are clearly strong evolutionary constraints on tissue-specific gene expression. A major challenge will be to understand the precise mechanisms by which many gene expression patterns remain similar despite extensive cis-regulatory restructuring.

Identification of direct T-box target genes in the developing zebrafish mesoderm

The zebrafish genes spadetail (spt) and no tail (ntl) encode T-box transcription factors that are important for early mesoderm development. Although much has been done to characterize these genes, the identity and location of target regulatory elements remain largely unknown. Here, we survey the genome for downstream target genes of the Spt and Ntl T-box transcription factors. We find evidence for extensive additive interactions towards gene activation and limited evidence for combinatorial and antagonistic interactions between the two factors. Using in vitro binding selection assays to define Spt- and Ntl-binding motifs, we searched for target regulatory sequence via a combination of binding motif searches and comparative genomics. We identified regulatory elements for tbx6 and deltaD, and, using chromatin immunoprecipitation, in vitro DNA binding assays and transgenic methods, we provide evidence that both are directly regulated by T-box transcription factors. We also find that deltaD is directly activated by T-box factors in the tail bud, where it has been implicated in starting the segmentation clock, suggesting that spt and ntl act upstream of this process.

Asymmetrical distribution of non-conserved regulatory sequences at PHOX2B is reflected at the ENCODE loci and illuminates a possible genome-wide trend

BACKGROUND

Transcriptional regulatory elements are central to development and interspecific phenotypic variation. Current regulatory element prediction tools rely heavily upon conservation for prediction of putative elements. Recent in vitro observations from the ENCODE project combined with in vivo analyses at the zebrafish phox2b locus suggests that a significant fraction of regulatory elements may fall below commonly applied metrics of conservation. We propose to explore these observations in vivo at the human PHOX2B locus, and also evaluate the potential evidence for genome-wide applicability of these observations through a novel analysis of extant data.

RESULTS

Transposon-based transgenic analysis utilizing a tiling path proximal to human PHOX2B in zebrafish recapitulates the observations at the zebrafish phox2b locus of both conserved and non-conserved regulatory elements. Analysis of human sequences conserved with previously identified zebrafish phox2b regulatory elements demonstrates that the orthologous sequences exhibit overlapping regulatory control. Additionally, analysis of non-conserved sequences scattered over 135 kb 5' to PHOX2B, provides evidence of non-conserved regulatory elements positively biased with close proximity to the gene. Furthermore, we provide a novel analysis of data from the ENCODE project, finding a non-uniform distribution of regulatory elements consistent with our in vivo observations at PHOX2B. These observations remain largely unchanged when one accounts for the sequence repeat content of the assayed intervals, when the intervals are sub-classified by biological role (developmental versus non-developmental), or by gene density (gene desert versus non-gene desert).

CONCLUSION

While regulatory elements frequently display evidence of evolutionary conservation, a fraction appears to be undetected by current metrics of conservation. In vivo observations at the PHOX2B locus, supported by our analyses of in vitro data from the ENCODE project, suggest that the risk of excluding non-conserved sequences in a search for regulatory elements may decrease as distance from the gene increases. Our data combined with the ENCODE data suggests that this may represent a genome wide trend.

Dispensability of mammalian DNA

In the lab, the cis-regulatory network seems to exhibit great functional redundancy. Many experiments testing enhancer activity of neighboring cis-regulatory elements show largely overlapping expression domains. Of recent interest, mice in which cis-regulatory ultraconserved elements were knocked out showed no obvious phenotype, further suggesting functional redundancy. Here, we present a global evolutionary analysis of mammalian conserved nonexonic elements (CNEs), and find strong evidence to the contrary. Given a set of CNEs conserved between several mammals, we characterize functional dispensability as the propensity for the ancestral element to be lost in mammalian species internal to the spanned species tree. We show that ultraconserved-like elements are over 300-fold less likely than neutral DNA to have been lost during rodent evolution. In fact, many thousands of noncoding loci under purifying selection display near uniform indispensability during mammalian evolution, largely irrespective of nucleotide conservation level. These findings suggest that many genomic noncoding elements possess functions that contribute noticeably to organism fitness in naturally evolving populations.

Evaluation of cis-regulatory function in zebrafish

As increasing numbers of vertebrate genomes are sequenced, comparative genomics offers tremendous promise to unveil mechanisms of transcriptional gene regulation on a large scale. However, the challenge of analysing immense amounts of sequence data and relating primary sequence to function is daunting. Several teleost species occupy crucial niches in the world of comparative genomics, as experimental model organisms of wide utility and living roadmaps of molecular evolution. Extant species have evolved after a teleost-specific genome duplication, and offer the opportunity to examine the evolution of thousands of duplicate gene pairs. Transgenesis in zebrafish is being increasingly employed to functionally examine non-coding sequences, from fish and mammals. Here, we discuss current approaches to the study of gene regulation in teleosts, and the promise of future research.

Evolution of the mammalian transcription factor binding repertoire via transposable elements

Identification of lineage-specific innovations in genomic control elements is critical for understanding transcriptional regulatory networks and phenotypic heterogeneity. We analyzed, from an evolutionary perspective, the binding regions of seven mammalian transcription factors (ESR1, TP53, MYC, RELA, POU5F1, SOX2, and CTCF) identified on a genome-wide scale by different chromatin immunoprecipitation approaches and found that only a minority of sites appear to be conserved at the sequence level. Instead, we uncovered a pervasive association with genomic repeats by showing that a large fraction of the bona fide binding sites for five of the seven transcription factors (ESR1, TP53, POU5F1, SOX2, and CTCF) are embedded in distinctive families of transposable elements. Using the age of the repeats, we established that these repeat-associated binding sites (RABS) have been associated with significant regulatory expansions throughout the mammalian phylogeny. We validated the functional significance of these RABS by showing that they are over-represented in proximity of regulated genes and that the binding motifs within these repeats have undergone evolutionary selection. Our results demonstrate that transcriptional regulatory networks are highly dynamic in eukaryotic genomes and that transposable elements play an important role in expanding the repertoire of binding sites.

Uncovering cis regulatory codes using synthetic promoter shuffling

Revealing the spectrum of combinatorial regulation of transcription at individual promoters is essential for understanding the complex structure of biological networks. However, the computations represented by the integration of various molecular signals at complex promoters are difficult to decipher in the absence of simple cis regulatory codes. Here we synthetically shuffle the regulatory architecture — operator sequences binding activators and repressors — of a canonical bacterial promoter. The resulting library of complex promoters allows for rapid exploration of promoter encoded logic regulation. Among all possible logic functions, NOR and ANDN promoter encoded logics predominate. A simple transcriptional cis regulatory code determines both logics, establishing a straightforward map between promoter structure and logic phenotype. The regulatory code is determined solely by the type of transcriptional regulation combinations: two repressors generate a NOR: NOT (a OR b) whereas a repressor and an activator generate an ANDN: a AND NOT b. Three-input versions of both logics, having an additional repressor as an input, are also present in the library. The resulting complex promoters cover a wide dynamic range of transcriptional strengths. Synthetic promoter shuffling represents a fast and efficient method for exploring the spectrum of complex regulatory functions that can be encoded by complex promoters. From an engineering point of view, synthetic promoter shuffling enables the experimental testing of the functional properties of complex promoters that cannot necessarily be inferred ab initio from the known properties of the individual genetic components. Synthetic promoter shuffling may provide a useful experimental tool for studying naturally occurring promoter shuffling.

The words of the regulatory code are arranged in a variable manner in highly conserved enhancers

The cis-regulatory regions of many developmental regulators and transcription factors are believed to be highly conserved in the genomes of vertebrate species, suggesting specific regulatory mechanisms for these gene classes. We functionally characterized five notochord enhancers, whose sequence is highly conserved, and systematically mutated two of them. Two subregions were identified to be essential for expression in the notochord of the zebrafish embryo. Synthetic enhancers containing the two essential regions in front of a TATA-box drive expression in the notochord while concatemerization of the subregions alone is not sufficient, indicating that the combination of the two sequence elements is required for notochord expression. Both regions are present in the five functionally characterized notochord enhancers. However, the position, the distance and relative orientation of the two sequence motifs can vary substantially within the enhancer sequences. This suggests that the regulatory grammar itself does not dictate the high evolutionary conservation between these orthologous cis-regulatory sequences. Rather, it represents a less well-conserved layer of sequence organization within these sequences.

Functional diversification of sonic hedgehog paralog enhancers identified by phylogenomic reconstruction

BACKGROUND

Cis-regulatory modules of developmental genes are targets of evolutionary changes that underlie the morphologic diversity of animals. Little is known about the 'grammar' of interactions between transcription factors and cis-regulatory modules and therefore about the molecular mechanisms that underlie changes in these modules, particularly after gene and genome duplications. We investigated the ar-C midline enhancer of sonic hedgehog (shh) orthologs and paralogs from distantly related vertebrate lineages, from fish to human, including the basal vertebrate Latimeria menadoensis.

RESULTS

We demonstrate that the sonic hedgehog a (shha) paralogs sonic hedgehog b (tiggy winkle hedgehog; shhb) genes of fishes have a modified ar-C enhancer, which specifies a diverged function at the embryonic midline. We have identified several conserved motifs that are indicative of putative transcription factor binding sites by local alignment of ar-C enhancers of numerous vertebrate sequences. To trace the evolutionary changes among paralog enhancers, phylogenomic reconstruction was carried out and lineage-specific motif changes were identified. The relation between motif composition and observed developmental differences was evaluated through transgenic functional analyses. Altering and exchanging motifs between paralog enhancers resulted in reversal of enhancer specificity in the floor plate and notochord. A model reconstructing enhancer divergence during vertebrate evolution was developed.

CONCLUSION

Our model suggests that the identified motifs of the ar-C enhancer function as binary switches that are responsible for specific activity between midline tissues, and that these motifs are adjusted during functional diversification of paralogs. The unraveled motif changes can also account for the complex interpretation of activator and repressor input signals within a single enhancer.

Organization of conserved elements near key developmental regulators in vertebrate genomes

Sequence conservation has traditionally been used as a means to target functional regions of complex genomes. In addition to its use in identifying coding regions of genes, the recent availability of whole genome data for a number of vertebrates has permitted high-resolution analyses of the noncoding "dark matter" of the genome. This has resulted in the identification of a large number of highly conserved sequence elements that appear to be preserved in all bony vertebrates. Further positional analysis of these conserved noncoding elements (CNEs) in the genome demonstrates that they cluster around genes involved in developmental regulation. This chapter describes the identification and characterization of these elements, with particular reference to their composition and organization.

Conserved non-coding sequences and transcriptional regulation

Genes coding for transcription factors and developmental regulators have a high likelihood to harbour cis-regulatory regions that are structurally conserved among orthologous genes in the vertebrate lineage. These regions can span up to several hundred basepairs with 70 and more percent sequence identity between fish and mammals. Even though this conservation is an efficient tool to discover cis-regulatory regions, we know little about why these specific genes maintain such highly conserved regulatory sequences. Here, we summarise work of the past few years on the regulatory modules of the sonic hedgehog and neurogenin1 genes. We will discuss the high sequence conservation of the regulatory elements in the context of models of enhancer evolution. Our data suggest that conservation of sequence does not necessarily imply a conserved function in other vertebrates.

TFCONES: a database of vertebrate transcription factor-encoding genes and their associated conserved noncoding elements

Background

Transcription factors (TFs) regulate gene transcription and play pivotal roles in various biological processes such as development, cell cycle progression, cell differentiation and tumor suppression. Identifying cis-regulatory elements associated with TF-encoding genes is a crucial step in understanding gene regulatory networks. To this end, we have used a comparative genomics approach to identify putative cis-regulatory elements associated with TF-encoding genes in vertebrates.

Description

We have created a database named TFCONES (Transcription Factor Genes & Associated COnserved Noncoding ElementS) () which contains all human, mouse and fugu TF-encoding genes and conserved noncoding elements (CNEs) associated with them. The CNEs were identified by gene-by-gene alignments of orthologous TF-encoding gene loci using MLAGAN. We also predicted putative transcription factor binding sites within the CNEs. A significant proportion of human-fugu CNEs contain experimentally defined binding sites for transcriptional activators and repressors, indicating that a majority of the CNEs may function as transcriptional regulatory elements. The TF-encoding genes that are involved in nervous system development are generally enriched for human-fugu CNEs. Users can retrieve TF-encoding genes and their associated CNEs by conducting a keyword search or by selecting a family of DNA-binding proteins.

Conclusion

The conserved noncoding elements identified in TFCONES represent a catalog of highly prioritized putative cis-regulatory elements of TF-encoding genes and are candidates for functional assay.

Identification and analysis of a conserved Tcfap2a intronic enhancer element required for expression in facial and limb bud mesenchyme

Tcfap2a, the gene encoding the mouse AP-2alpha transcription factor, is required for normal development of multiple structures during embryogenesis, including the face and limbs. Using comparative sequence analysis and transgenic-mouse experiments we have identified an intronic enhancer within this gene that directs expression to the face and limb mesenchyme. There are two conserved sequence blocks within this intron, and the larger of these directs tissue-specific activity and is found in all vertebrate Tcfap2a genes analyzed. To assess the role of the enhancer in regulating endogenous mouse Tcfap2a expression, we have deleted this cis-regulatory sequence from the genome. Loss of this element severely impairs Tcfap2a expression in the limb bud mesenchyme but generates only a modest reduction in the facial mesenchyme. The reduction in Tcfap2a transcription is accompanied by altered patterning of the forelimb, resulting in postaxial polydactyly. These results indicate that the major role for this enhancer resides within the limb bud, and it serves to maintain a level of Tcfap2a expression that limits the size of the hand plate and the associated number of digit primordia. The potential role of this cis-acting sequence in modeling the size and shape of the face and limbs during evolution is discussed.

CONDOR: a database resource of developmentally associated conserved non-coding elements

Background

Comparative genomics is currently one of the most popular approaches to study the regulatory architecture of vertebrate genomes. Fish-mammal genomic comparisons have proved powerful in identifying conserved non-coding elements likely to be distal cis-regulatory modules such as enhancers, silencers or insulators that control the expression of genes involved in the regulation of early development. The scientific community is showing increasing interest in characterizing the function, evolution and language of these sequences. Despite this, there remains little in the way of user-friendly access to a large dataset of such elements in conjunction with the analysis and the visualization tools needed to study them.

Description

Here we present CONDOR (COnserved Non-coDing Orthologous Regions) available at: . In an interactive and intuitive way the website displays data on > 6800 non-coding elements associated with over 120 early developmental genes and conserved across vertebrates. The database regularly incorporates results of ongoing in vivo zebrafish enhancer assays of the CNEs carried out in-house, which currently number ~100. Included and highlighted within this set are elements derived from duplication events both at the origin of vertebrates and more recently in the teleost lineage, thus providing valuable data for studying the divergence of regulatory roles between paralogs. CONDOR therefore provides a number of tools and facilities to allow scientists to progress in their own studies on the function and evolution of developmental cis-regulation.

Conclusion

By providing access to data with an approachable graphics interface, the CONDOR database presents a rich resource for further studies into the regulation and evolution of genes involved in early development.

Raising the estimate of functional human sequences

While less than 1.5% of the mammalian genome encodes proteins, it is now evident that the vast majority is transcribed, mainly into non-protein-coding RNAs. This raises the question of what fraction of the genome is functional, i.e., composed of sequences that yield functional products, are required for the expression (regulation or processing) of these products, or are required for chromosome replication and maintenance. Many of the observed noncoding transcripts are differentially expressed, and, while most have not yet been studied, increasing numbers are being shown to be functional and/or trafficked to specific subcellular locations, as well as exhibit subtle evidence of selection. On the other hand, analyses of conservation patterns indicate that only approximately 5% (3%-8%) of the human genome is under purifying selection for functions common to mammals. However, these estimates rely on the assumption that reference sequences (usually ancient transposon-derived sequences) have evolved neutrally, which may not be the case, and if so would lead to an underestimate of the fraction of the genome under evolutionary constraint. These analyses also do not detect functional sequences that are evolving rapidly and/or have acquired lineage-specific functions. Indeed, many regulatory sequences and known functional noncoding RNAs, including many microRNAs, are not conserved over significant evolutionary distances, and recent evidence from the ENCODE project suggests that many functional elements show no detectable level of sequence constraint. Thus, it is likely that much more than 5% of the genome encodes functional information, and although the upper bound is unknown, it may be considerably higher than currently thought.

Rodent Aanat: intronic E-box sequences control tissue specificity but not rhythmic expression in the pineal gland

Arylalkylamine N-acetyltransferase (Aanat) is the penultimate enzyme in the serotonin-N-acetylserotonin-melatonin pathway. It is nearly exclusively expressed in the pineal gland and the retina. A marked rhythm of Aanat gene expression in the rat pineal is mediated by cyclic AMP response elements located in the promoter and first intron. Intron 1 also contains E-box elements, which mediate circadian gene expression in other cells. Here we examined whether these elements contribute to rhythmic Aanat expression in the pineal gland. This was done using transgenic rats carrying Aanat transgenes with mutant E-box elements. Circadian expression of Aanat transgenes was not altered by these mutations. However, these mutations enhanced ectopic expression establishing that the intronic Aanat E-box elements contribute to the gene's pineal specific expression. A similar role of the Aanat E-box has been reported in zebrafish, indicating that Aanat E-box mediated silencing is a conserved feature of vertebrate biology.

A new paradigm for developmental biology

It is usually thought that the development of complex organisms is controlled by protein regulatory factors and morphogenetic signals exchanged between cells and differentiating tissues during ontogeny. However, it is now evident that the majority of all animal genomes is transcribed, apparently in a developmentally regulated manner, suggesting that these genomes largely encode RNA machines and that there may be a vast hidden layer of RNA regulatory transactions in the background. I propose that the epigenetic trajectories of differentiation and development are primarily programmed by feed-forward RNA regulatory networks and that most of the information required for multicellular development is embedded in these networks, with cell–cell signalling required to provide important positional information and to correct stochastic errors in the endogenous RNA-directed program.

Multiple non-collinear TF-map alignments of promoter regions

Background

The analysis of the promoter sequence of genes with similar expression patterns is a basic tool to annotate common regulatory elements. Multiple sequence alignments are on the basis of most comparative approaches. The characterization of regulatory regions from co-expressed genes at the sequence level, however, does not yield satisfactory results in many occasions as promoter regions of genes sharing similar expression programs often do not show nucleotide sequence conservation.

Results

In a recent approach to circumvent this limitation, we proposed to align the maps of predicted transcription factors (referred as TF-maps) instead of the nucleotide sequence of two related promoters, taking into account the label of the corresponding factor and the position in the primary sequence. We have now extended the basic algorithm to permit multiple promoter comparisons using the progressive alignment paradigm. In addition, non-collinear conservation blocks might now be identified in the resulting alignments. We have optimized the parameters of the algorithm in a small, but well-characterized collection of human-mouse-chicken-zebrafish orthologous gene promoters.

Conclusion

Results in this dataset indicate that TF-map alignments are able to detect high-level regulatory conservation at the promoter and the 3'UTR gene regions, which cannot be detected by the typical sequence alignments. Three particular examples are introduced here to illustrate the power of the multiple TF-map alignments to characterize conserved regulatory elements in absence of sequence similarity. We consider this kind of approach can be extremely useful in the future to annotate potential transcription factor binding sites on sets of co-regulated genes from high-throughput expression experiments.

Human GLI3 intragenic conserved non-coding sequences are tissue-specific enhancers

The zinc-finger transcription factor GLI3 is a key regulator of development, acting as a primary transducer of Sonic hedgehog (SHH) signaling in a combinatorial context dependent fashion controlling multiple patterning steps in different tissues/organs. A tight temporal and spatial control of gene expression is indispensable, however, cis-acting sequence elements regulating GLI3 expression have not yet been reported. We show that 11 ancient genomic DNA signatures, conserved from the pufferfish Takifugu (Fugu) rubripes to man, are distributed throughout the introns of human GLI3. They map within larger conserved non-coding elements (CNEs) that are found in the tetrapod lineage. Full length CNEs transiently transfected into human cell cultures acted as cell type specific enhancers of gene transcription. The regulatory potential of these elements is conserved and was exploited to direct tissue specific expression of a reporter gene in zebrafish embryos. Assays of deletion constructs revealed that the human-Fugu conserved sequences within the GLI3 intronic CNEs were essential but not sufficient for full-scale transcriptional activation. The enhancer activity of the CNEs is determined by a combinatorial effect of a core sequence conserved between human and teleosts (Fugu) and flanking tetrapod-specific sequences, suggesting that successive clustering of sequences with regulatory potential around an ancient, highly conserved nucleus might be a possible mechanism for the evolution of cis-acting regulatory elements.

Parallel evolution of conserved non-coding elements that target a common set of developmental regulatory genes from worms to humans

BACKGROUND

The human genome contains thousands of non-coding sequences that are often more conserved between vertebrate species than protein-coding exons. These highly conserved non-coding elements (CNEs) are associated with genes that coordinate development, and have been proposed to act as transcriptional enhancers. Despite their extreme sequence conservation in vertebrates, sequences homologous to CNEs have not been identified in invertebrates.

RESULTS

Here we report that nematode genomes contain an alternative set of CNEs that share sequence characteristics, but not identity, with their vertebrate counterparts. CNEs thus represent a very unusual class of sequences that are extremely conserved within specific animal lineages yet are highly divergent between lineages. Nematode CNEs are also associated with developmental regulatory genes, and include well-characterized enhancers and transcription factor binding sites, supporting the proposed function of CNEs as cis-regulatory elements. Most remarkably, 40 of 156 human CNE-associated genes with invertebrate orthologs are also associated with CNEs in both worms and flies.

CONCLUSION

A core set of genes that regulate development is associated with CNEs across three animal groups (worms, flies and vertebrates). We propose that these CNEs reflect the parallel evolution of alternative enhancers for a common set of developmental regulatory genes in different animal groups. This 're-wiring' of gene regulatory networks containing key developmental coordinators was probably a driving force during the evolution of animal body plans. CNEs may, therefore, represent the genomic traces of these 'hard-wired' core gene regulatory networks that specify the development of each alternative animal body plan.

The relationship between non-protein-coding DNA and eukaryotic complexity

There are two intriguing paradoxes in molecular biology--the inconsistent relationship between organismal complexity and (1) cellular DNA content and (2) the number of protein-coding genes--referred to as the C-value and G-value paradoxes, respectively. The C-value paradox may be largely explained by varying ploidy. The G-value paradox is more problematic, as the extent of protein coding sequence remains relatively static over a wide range of developmental complexity. We show by analysis of sequenced genomes that the relative amount of non-protein-coding sequence increases consistently with complexity. We also show that the distribution of introns in complex organisms is non-random. Genes composed of large amounts of intronic sequence are significantly overrepresented amongst genes that are highly expressed in the nervous system, and amongst genes downregulated in embryonic stem cells and cancers. We suggest that the informational paradox in complex organisms may be explained by the expansion of cis-acting regulatory elements and genes specifying trans-acting non-protein-coding RNAs.