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.     

Multiple Promoter Targeting Sequences exist in Abdominal-B to regulate long-range gene activation

In complex genomes, insulators set up chromatin domain boundaries and protect promoters from inappropriate activation by enhancers from neighboring genes. The Drosophila Abdominal-B locus uses insulator elements to organize its large regulatory region into several body segment-specific chromatin domains. This organization leads to a problem in enhancer-promoter communication, that is, how do distal enhancers activate the Abd-B promoter when there are several insulators in between? This issue is partially resolved by the Promoter Targeting Sequence, which can overcome the enhancer blocking effect of an insulator. In this study, we describe a new Promoter Targeting Sequence, PTS-6, from the Abd-B 3' regulatory region. PTS-6, comprised of approximately 200 bp, was found to bypass both homologous Abdominal-B insulators, such as Fab-7 and Fab-8, and a heterologous insulator, suHw. Most importantly, it also overcomes a combination of two insulators such as Fab-7/Fab-8. Thus, PTS-6 could, in principle, target remote enhancers that are separated from the Abd-B promoter by multiple insulators. In addition, PTS-6 selectively targets the distal enhancer to only one transgenic promoter, and it strongly facilitates Abd-B enhancers. These results suggest that promoter targeting is necessary for long-range enhancer-promoter communication in Abd-B, and PTS elements could be a common occurrence in large, complex genetic loci.     

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.     

The role of a retinoic acid response element in establishing the anterior neural expression border of Hoxd4 transgenes

The zebrafish hoxd4a locus was compared to its murine ortholog, Hoxd4. The sequence of regulatory elements, including a DR5 type retinoic acid response element (RARE) required for Hoxd4 neural enhancer activity, are highly conserved. Additionally, zebrafish and mouse neural enhancers function identically in transgenic mouse embryos. We tested whether sequence conservation reflects functional importance by altering the spacing and sequence of the RARE in the Hoxd4 neural enhancer. Stabilizing receptor-DNA interactions did not anteriorize transgene expression. By contrast, conversion of the RARE from a DR5 to a DR2 type element decreased receptor-DNA stability and posteriorized expression. Hence, the setting of the Hox anterior expression border is not a simple function of the affinity of retinoid receptors for their cognate element.     

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.     

Diverse functions of distal regulatory elements at the IFNG locus

Previous studies have identified multiple conserved noncoding sequences (CNS) at the mouse Ifng locus sufficient for enhancer activity in cell-based assays. These studies do not directly address biology of the human IFNG locus in a genomic setting. IFNG enhancers may be functionally redundant or each may be functionally unique. We test the hypothesis that each IFNG enhancer has a unique necessary function using a bacterial artificial chromosome transgenic model. We find that CNS-30, CNS-4, and CNS+20 are required at distinct stages of Th1 differentiation, whereas CNS-16 has a repressive role in Th1 and Th2 cells. CNS+20 is required for IFN-γ expression by memory Th1 cells and NKT cells. CNS-4 is required for IFN-γ expression by effector Th1 cells. In contrast, CNS-16, CNS-4, and CNS+20 are each partially required for human IFN-γ expression by NK cells. Thus, IFNG CNS enhancers have redundant necessary functions in NK cells but unique necessary functions in Th cells. These results also demonstrate that distinct CNSs are required to transcribe IFNG at each stage of the Th1 differentiation pathway.     

A latitudinal cline in the Chinook salmon (Oncorhynchus tshawytscha) Clock gene: evidence for selection on PolyQ length variants

A critical seasonal event for anadromous Chinook salmon (Oncorhynchus tshawytscha) is the time at which adults migrate from the ocean to breed in freshwater. We investigated whether allelic variation at the circadian rhythm genes, OtsClock1a and OtsClock1b, underlies genetic control of migration timing among 42 populations in North America. We identified eight length variants of the functionally important polyglutamine repeat motif (PolyQ) of OtsClock1b while OtsClock1a PolyQ was highly conserved. We found evidence of a latitudinal cline in average allele length and frequency of the two most common OtsClock1b alleles. The shorter 335 bp allele increases in frequency with decreasing latitude while the longer 359 bp allele increases in frequency at higher latitudes. Comparison to 13 microsatellite loci showed that 335 and 359 bp deviate significantly from neutral expectations. Furthermore, a hierarchical gene diversity analysis based on OtsClock1b PolyQ variation revealed that run timing explains 40.9 per cent of the overall genetic variance among populations. By contrast, an analysis based on 13 microsatellite loci showed that run timing explains only 13.2 per cent of the overall genetic variance. Our findings suggest that length polymorphisms in OtsClock1b PolyQ may be maintained by selection and reflect an adaptation to ecological factors correlated with latitude, such as the seasonally changing day length.     

Evolution acts on enhancer organization to fine-tune gradient threshold readouts

The elucidation of principles governing evolution of gene regulatory sequence is critical to the study of metazoan diversification. We are therefore exploring the structure and organizational constraints of regulatory sequences by studying functionally equivalent cis-regulatory modules (CRMs) that have been evolving in parallel across several loci. Such an independent dataset allows a multi-locus study that is not hampered by nonfunctional or constrained homology. The neurogenic ectoderm enhancers (NEEs) of Drosophila melanogaster are one such class of coordinately regulated CRMs. The NEEs share a common organization of binding sites and as a set would be useful to study the relationship between CRM organization and CRM activity across evolving lineages. We used the D. melanogaster transgenic system to screen for functional adaptations in the NEEs from divergent drosophilid species. We show that the individual NEE modules across a genome in any one lineage have independently evolved adaptations to compensate for lineage-specific developmental and/or genomic changes. Specifically, we show that both the site composition and the site organization of NEEs have been finely tuned by distinct, lineage-specific selection pressures in each of the three divergent species that we have examined: D. melanogaster, D. pseudoobscura, and D. virilis. Furthermore, by precisely altering the organization of NEEs with different morphogen gradient threshold readouts, we show that CRM organizational evolution is sufficient for explaining changes in enhancer activity. Thus, evolution can act on CRM organization to fine-tune morphogen gradient threshold readouts over a wide dynamic range. Our study demonstrates that equivalence classes of CRMs are powerful tools for detecting lineage-specific adaptations by gene regulatory sequences.     

Evolution of lineage-specific functions in ancient cis-regulatory modules

Morphological evolution is driven both by coding sequence variation and by changes in regulatory sequences. However, how cis-regulatory modules (CRMs) evolve to generate entirely novel expression domains is largely unknown. Here, we reconstruct the evolutionary history of a lens enhancer located within a CRM that not only predates the lens, a vertebrate innovation, but bilaterian animals in general. Alignments of orthologous sequences from different deuterostomes sub-divide the CRM into a deeply conserved core and a more divergent flanking region. We demonstrate that all deuterostome flanking regions, including invertebrate sequences, activate gene expression in the zebrafish lens through the same ancient cluster of activator sites. However, levels of gene expression vary between species due to the presence of repressor motifs in flanking region and core. These repressor motifs are responsible for the relatively weak enhancer activity of tetrapod flanking regions. Ray-finned fish, however, have gained two additional lineage-specific activator motifs which in combination with the ancient cluster of activators and the core constitute a potent lens enhancer. The exploitation and modification of existing regulatory potential in flanking regions but not in the highly conserved core might represent a more general model for the emergence of novel regulatory functions in complex CRMs.     

Characterization of a cardiac-specific enhancer, which directs {alpha}-cardiac actin gene transcription in the mouse adult heart

Expression of the mouse alpha-cardiac actin gene in skeletal and cardiac muscle is regulated by enhancers lying 5' to the proximal promoter. Here we report the characterization of a cardiac-specific enhancer located within -2.354/-1.36 kbp of the gene, which is active in cardiocytes but not in C2 skeletal muscle cells. In vivo it directs reporter gene expression to the adult heart, where the proximal promoter alone is inactive. An 85-bp region within the enhancer is highly conserved between human and mouse and contains a central AT-rich site, which is essential for enhancer activity. This site binds myocyte enhancer factor (MEF)2 factors, principally MEF2D and MEF2A in cardiocyte nuclear extracts. These results are discussed in the context of MEF2 activity and of the regulation of the alpha-cardiac actin locus.