Efficient identification of regulatory sequences in the chicken genome by a powerful combination of embryo electroporation and genome comparison

Recently expanded knowledge of gene regulation clearly indicates that the regulatory sequences of a gene, usually identified as enhancers, are widely distributed in the gene locus, revising the classical view that they are clustered in the vicinity of genes. To identify regulatory sequences for Sox2 expression governing early neurogenesis, we scanned the 50-kb region of the chicken Sox2 locus for enhancer activity utilizing embryo electroporation, resulting in identification of a number of enhancers scattered throughout the analyzed genomic span. The 'pan-neural' Sox2 expression in early embryos is actually brought about by the composite activities of five separate enhancers with distinct spatio-temporal specificities. These and other functionally defined enhancers exactly correspond to extragenic sequence blocks that are conspicuously conserved between the chicken and mammalian genomes and that are embedded in sequences with a wide range of sequence conservation between humans and mice. The sequences conserved between amniotes and teleosts correspond to subregions of the enhancer subsets which presumably represent core motifs of the enhancers, and the limited conservation partly reflects divergent expression patterns of the gene. The phylogenic distance between the chicken and mammals appears optimal for identifying a battery of genetic regulatory elements as conserved sequence blocks, and chicken embryo electroporation facilitates functional characterization of these elements.     

A survey of ancient conserved non-coding elements in the PAX6 locus reveals a landscape of interdigitated cis-regulatory archipelagos

Biological differences between cell types and developmental processes are characterised by differences in gene expression profiles. Gene-distal enhancers are key components of the regulatory networks that specify the tissue-specific expression patterns driving embryonic development and cell fate decisions, and variations in their sequences are a major contributor to genetic disease and disease susceptibility. Despite advances in the methods for discovery of putative cis-regulatory sequences, characterisation of their spatio-temporal enhancer activities in a mammalian model system remains a major bottle-neck. We employed a strategy that combines gnathostome sequence conservation with transgenic mouse and zebrafish reporter assays to survey the genomic locus of the developmental control gene PAX6 for the presence of novel cis-regulatory elements. Sequence comparison between human and the cartilaginous elephant shark (Callorhinchus milii) revealed several ancient gnathostome conserved non-coding elements (agCNEs) dispersed widely throughout the PAX6 locus, extending the range of the known PAX6 cis-regulatory landscape to contain the full upstream PAX6-RCN1 intergenic region. Our data indicates that ancient conserved regulatory sequences can be tested effectively in transgenic zebrafish even when not conserved in zebrafish themselves. The strategy also allows efficient dissection of compound regulatory regions previously assessed in transgenic mice. Remarkable overlap in expression patterns driven by sets of agCNEs indicates that PAX6 resides in a landscape of multiple tissue-specific regulatory archipelagos.     

Identification and experimental validation of splicing regulatory elements in Drosophila melanogaster reveals functionally conserved splicing enhancers in metazoans

RNA sequence elements involved in the regulation of pre-mRNA splicing have previously been identified in vertebrate genomes by computational methods. Here, we apply such approaches to predict splicing regulatory elements in Drosophila melanogaster and compare them with elements previously found in the human, mouse, and pufferfish genomes. We identified 99 putative exonic splicing enhancers (ESEs) and 231 putative intronic splicing enhancers (ISEs) enriched near weak 5' and 3' splice sites of constitutively spliced introns, distinguishing between those found near short and long introns. We found that a significant proportion (58%) of fly enhancer sequences were previously reported in at least one of the vertebrates. Furthermore, 20% of putative fly ESEs were previously identified as ESEs in human, mouse, and pufferfish; while only two fly ISEs, CTCTCT and TTATAA, were identified as ISEs in all three vertebrate species. Several putative enhancer sequences are similar to characterized binding-site motifs for Drosophila and mammalian splicing regulators. To provide additional evidence for the function of putative ISEs, we separately identified 298 intronic hexamers significantly enriched within sequences phylogenetically conserved among 15 insect species. We found that 73 putative ISEs were among those enriched in conserved regions of the D. melanogaster genome. The functions of nine enhancer sequences were verified in a heterologous splicing reporter, demonstrating that these sequences are sufficient to enhance splicing in vivo. Taken together, these data identify a set of predicted positive-acting splicing regulatory motifs in the Drosophila genome and reveal regulatory sequences that are present in distant metazoan genomes.     

Conservation and diversity in the cis-regulatory networks that integrate information controlling expression of Hoxa2 in hindbrain and cranial neural crest cells in vertebrates

The Hoxa2 and Hoxb2 genes are members of paralogy group II and display segmental patterns of expression in the developing vertebrate hindbrain and cranial neural crest cells. Functional analyses have demonstrated that these genes play critical roles in regulating morphogenetic pathways that direct the regional identity and anteroposterior character of hindbrain rhombomeres and neural crest-derived structures. Transgenic regulatory studies have also begun to characterize enhancers and cis-elements for those mouse and chicken genes that direct restricted patterns of expression in the hindbrain and neural crest. In light of the conserved role of Hoxa2 in neural crest patterning in vertebrates and the similarities between paralogs, it is important to understand the extent to which common regulatory networks and elements have been preserved between species and between paralogs. To investigate this problem, we have cloned and sequenced the intergenic region between Hoxa2 and Hoxa3 in the chick HoxA complex and used it for making comparative analyses with the respective human, mouse, and horn shark regions. We have also used transgenic assays in mouse and chick embryos to test the functional activity of Hoxa2 enhancers in heterologous species. Our analysis reveals that three of the critical individual components of the Hoxa2 enhancer region from mouse necessary for hindbrain expression (Krox20, BoxA, and TCT motifs) have been partially conserved. However, their number and organization are highly varied for the same gene in different species and between paralogs within a species. Other essential mouse elements appear to have diverged or are absent in chick and shark. We find the mouse r3/r5 enhancer fails to work in chick embryos and the chick enhancer works poorly in mice. This implies that new motifs have been recruited or utilized to mediate restricted activity of the enhancer in other species. With respect to neural crest regulation, cis-components are embedded among the hindbrain control elements and are highly diverged between species. Hence, there has been no widespread conservation of sequence identity over the entire enhancer domain from shark to humans, despite the common function of these genes in head patterning. This provides insight into how apparently equivalent regulatory regions from the same gene in different species have evolved different components to potentiate their activity in combination with a selection of core components.     

Enhancer evolution in chordates: Lessons from functional analyses of cephalochordate cis-regulatory modules

Chordates comprise three major groups, cephalochordates (amphioxus), tunicates (urochordates), and vertebrates. Since cephalochordates were the early branching group, comparisons between amphioxus and other chordates help us to speculate about ancestral chordates. Here, I summarize accumulating data from functional studies analyzing amphioxus cis-regulatory modules (CRMs) in model systems of other chordate groups, such as mice, chickens, clawed frogs, fish, and ascidians. Conservatism and variability of CRM functions illustrate how gene regulatory networks have evolved in chordates. Amphioxus CRMs, which correspond to CRMs deeply conserved among animal phyla, govern reporter gene expression in conserved expression domains of the putative target gene in host animals. In addition, some CRMs located in similar genomic regions (intron, upstream, or downstream) also possess conserved activity, even though their sequences are divergent. These conservative CRM functions imply ancestral genomic structures and gene regulatory networks in chordates. However, interestingly, if expression patterns of amphioxus genes do not correspond to those of orthologs of experimental models, some amphioxus CRMs recapitulate expression patterns of amphioxus genes, but not those of endogenous genes, suggesting that these amphioxus CRMs are close to the ancestral states of chordate CRMs, while vertebrates/tunicates innovated new CRMs to reconstruct gene regulatory networks subsequent to the divergence of the cephalochordates. Alternatively, amphioxus CRMs may have secondarily lost ancestral CRM activity and evolved independently. These data help to solve fundamental questions of chordate evolution, such as neural crest cells, placodes, a forebrain/midbrain, and genome duplication. Experimental validation is crucial to verify CRM functions and evolution.     

Dissection of cis-regulatory elements in the C. elegans Hox gene egl-5 promoter

Hox genes are highly conserved segmental identity genes well known for their complex expression patterns and divergent targets. Here we present an analysis of cis-regulatory elements in the Caenorhabditis elegans Hox gene egl-5, which is expressed in multiple tissues in the posterior region of the nematode. We have utilized phylogenetic footprinting to efficiently identify cis-regulatory elements and have characterized these with gfp reporters and tissue-specific rescue experiments. We have found that the complex expression pattern of egl-5 is the cumulative result of the activities of multiple tissue or local region-specific activator sequences that are conserved both in sequence and near-perfect order in the related nematode Caenorhabditis briggsae. Two conserved regulatory blocks analyzed in detail contain multiple sites for both positively and negatively acting factors. One of these regions may promote activation of egl-5 in certain cells via the Wnt pathway. Positively acting regions are repressed in inappropriate tissues by additional negative pathways acting at other sites within the promoter. Our analysis has allowed us to implicate several new regulatory factors significant to the control of egl-5 expression.     

Validation of skeletal muscle cis-regulatory module predictions reveals nucleotide composition bias in functional enhancers

We performed a genome-wide scan for muscle-specific cis-regulatory modules (CRMs) using three computational prediction programs. Based on the predictions, 339 candidate CRMs were tested in cell culture with NIH3T3 fibroblasts and C2C12 myoblasts for capacity to direct selective reporter gene expression to differentiated C2C12 myotubes. A subset of 19 CRMs validated as functional in the assay. The rate of predictive success reveals striking limitations of computational regulatory sequence analysis methods for CRM discovery. Motif-based methods performed no better than predictions based only on sequence conservation. Analysis of the properties of the functional sequences relative to inactive sequences identifies nucleotide sequence composition can be an important characteristic to incorporate in future methods for improved predictive specificity. Muscle-related TFBSs predicted within the functional sequences display greater sequence conservation than non-TFBS flanking regions. Comparison with recent MyoD and histone modification ChIP-Seq data supports the validity of the functional regions.