Search for enhancers: teleost models in comparative genomic and transgenic analysis of cis regulatory elements

Zebrafish Embryos and Larvae as Alternative Animal Models for Toxicity Testing

Prerequisite to any biological laboratory assay employing living animals is consideration about its necessity, feasibility, ethics and the potential harm caused during an experiment. The imperative of these thoughts has led to the formulation of the 3R-principle, which today is a pivotal scientific standard of animal experimentation worldwide. The rising amount of laboratory investigations utilizing living animals throughout the last decades, either for regulatory concerns or for basic science, demands the development of alternative methods in accordance with 3R to help reduce experiments in mammals. This demand has resulted in investigation of additional vertebrate species displaying favourable biological properties. One prominent species among these is the zebrafish (Danio rerio), as these small laboratory ray-finned fish are well established in science today and feature outstanding biological characteristics. In this review, we highlight the advantages and general prerequisites of zebrafish embryos and larvae before free-feeding stages for toxicological testing, with a particular focus on cardio-, neuro, hepato- and nephrotoxicity. Furthermore, we discuss toxicokinetics, current advances in utilizing zebrafish for organ toxicity testing and highlight how advanced laboratory methods (such as automation, advanced imaging and genetic techniques) can refine future toxicological studies in this species.

On the importance of evolutionary constraint for regulatory sequence identification

Regulation of gene expression relies on the activity of specialized genomic elements, enhancers or silencers, distributed over sometimes large distance from their target gene promoters. A significant part of vertebrate genomes consists in such regulatory elements, but their identification and that of their target genes remains challenging, due to the lack of clear signature at the nucleotide level. For many years the main hallmark used for identifying functional elements has been their sequence conservation between genomes of distant species, indicative of purifying selection. More recently, genome-wide biochemical assays have opened new avenues for detecting regulatory regions, shifting attention away from evolutionary constraints. Here, we review the respective contributions of comparative genomics and biochemical assays for the definition of regulatory elements and their targets and advocate that both sequence conservation and preserved synteny, taken as signature of functional constraint, remain essential tools in this task.

A Perspective on the Value of Aquatic Models in Biomedical Research

For at least 150 years, biological scientists have congregated at marine laboratories, located at the edge of the sea, to explore aquatic life. The purpose of this minireview is to offer a brief perspective on the relevance of this activity to our knowledge of human physiology and disease, drawing heavily on the experience of the authors and without attempting to offer a comprehensive history of the many contributions of marine models to biomedical research.

Initiation of cyp26a1 Expression in the Zebrafish Anterior Neural Plate by a Novel Cis-Acting Element

Early patterning of the vertebrate neural plate involves a complex hierarchy of inductive interactions orchestrated by signalling molecules and their antagonists. The morphogen retinoic acid, together with the Cyp26 enzymes which degrade it, play a central role in this process. The cyp26a1 gene expressed in the anterior neural plate thus contributes to the fine modulation of the rostrocaudal retinoic acid gradient. Despite this important role of cyp26a1 in early brain formation, the mechanisms that control its expression in the anterior neural plate are totally unknown. Here, we present the isolation of a 310-base-pair DNA element adjacent to cyp26a1 promoter, displaying enhancer activity restricted to the anterior neural plate of the zebrafish gastrula. We show that unlike that of cyp26a1, expression driven by this cyp26a1 anterior neural plate element (cANE) is independent of retinoic acid. Through deletion analysis, we identify a 12-nucleotide motif essential for cANE activity. A consensus bipartite binding site for SoxB:Oct transcription factors overlaps with this motif. Mutational analysis suggests that SoxB binding is essential for its activity. We discuss the contribution of this study to the elucidation of the regulatory hierarchy involved in early neural plate patterning.

Advancing biology through a deeper understanding of zebrafish ecology and evolution

Over the last two decades, the zebrafish has joined the ranks of premier model organisms for biomedical research, with a full suite of tools and genomic resources. Yet we still know comparatively little about its natural history. Here I review what is known about the natural history of the zebrafish, where significant gaps in our knowledge remain, and how a fuller appreciation of this organism's ecology and behavior, population genetics, and phylogeny can inform a variety of research endeavors.

Making reporter gene constructs to analyze cis-regulatory elements

Cis-regulatory sequences control when, where, and how much genes are transcribed. A better understanding on these elements is a fundamental keystone to better understand development, cell differentiation, and morphogenesis. Several methods based on in silico analysis or ChIP-seq experiments have been developed to detect cis-acting sequences. Here, we describe a protocol to isolate such sequences from genomic DNA and to clone them into expression vectors for functional assays using the Gateway cloning technology.

Characterization of zebrafish intestinal smooth muscle development using a novel sm22α-b promoter

Smooth muscle cells provide structural support for many tissues and control essential physiological processes, such as blood pressure and gastrointestinal motility. Relatively little is known about the early stages of intestinal smooth muscle development and its relationship to the development of the enteric nervous system, which regulates intestinal motility. Here, we report an evolutionarily conserved 523 base pair regulatory element within the promoter of the zebrafish sm22α-b (transgelin1) gene that directs transgene expression in smooth muscle cells of the intestine and other tissues. Comparative genomic analysis identified a conserved motif within this element consisting of two Serum Response Factor binding sites that is also present in the promoters of many mammalian smooth muscle genes. We established a stable line expressing GFP in smooth muscle cell and used this line to describe lineage relationships among cells within different intestinal smooth muscle layers and their co-development with the enteric nervous system (ENS).

Transcriptional enhancers in animal development and evolution

Regulatory DNAs serve as templates to bring weakly interacting transcription factors into close proximity so they can work synergistically to switch genes on and off in time and space. Most of these regulatory DNAs are enhancers that can work over long distances--a million base pairs or more in mammals--to control gene expression. Critical enhancers are sometimes even found within the introns of neighboring genes. This review summarizes well-defined examples of enhancers controlling key processes in animal development. Potential mechanisms of transcriptional synergy are discussed with regard to enhancer structure and contemporary ChIP-sequencing assays, whereby just a small fraction of the observed binding sites represent bona fide regulatory DNAs. Finally, there is a discussion of how enhancer evolution can produce novelty in animal morphology and of the prospects for reconstructing transitions in animal evolution by introducing derived enhancers in basal ancestors.

Conservation of shh cis-regulatory architecture of the coelacanth is consistent with its ancestral phylogenetic position

Background

The modern coelacanth (Latimeria) is the extant taxon of a basal sarcopterygian lineage and sister group to tetrapods. Apart from certain apomorphic traits, its morphology is characterized by a high degree of retention of ancestral vertebrate structures and little morphological change. An insight into the molecular evolution that may explain the unchanged character of Latimeria morphology requires the analysis of the expression patterns of developmental regulator genes and their cis-regulatory modules (CRMs).

Results

We describe the comparative and functional analysis of the sonic hedgehog (shh) genomic region of Latimeria menadoensis. Several putative enhancers in the Latimeria shh locus have been identified by comparisons to sarcopterygian and actinopterygian extant species. Specific sequence conservation with all known actinopterygian enhancer elements has been detected. However, these elements are selectively missing in more recently diverged actinopterygian and sarcopterygian species. The functionality of the putative Latimeria enhancers was confirmed by reporter gene expression analysis in transient transgenic zebrafish and chick embryos.

Conclusions

Latimeria shh CRMs represent the ancestral set of enhancers that have emerged before the split of lobe-finned and ray-finned fishes. In contrast to lineage-specific losses and differentiations in more derived lineages, Latimeria shh enhancers reveal low levels of sequence diversification. High overall sequence conservation of shh conserved noncoding elements (CNE) is consistent with the general trend of high levels of conservation of noncoding DNA in the slowly evolving Latimeria genome.

Comparative genomics in teleost species: Knowledge transfer by linking the genomes of model and non-model fish species

Comparative genomics is a powerful tool to transfer knowledge coming from model fish species to non-model fish species of economic or/and evolutionary interest. Such transfer is of importance as functional studies either are difficult to perform with most non-model species. The first comparative map constructed using the human and the chimpanzee genome allowed the identification of putative orthologues. Although comparative mapping in teleosts is still in its infancy, five model teleost genomes from different orders have been fully sequenced to date and the sequencing of several commercially important species are also underway or near completion. The accessibility of these whole genome sequences and rapid developments in genomics of fish species are paving the way towards new and valuable research in comparative genetics and genomics. With the accumulation of information in model species, the genetic and genomic characterization of non-model, but economically, physiologically or evolutionary important species is now feasible. Furthermore, comparison of low coverage gene maps of non-model fish species against fully sequenced fish species will enhance the efficiency of candidate gene identification projected for quantitative trait loci (QTL) scans for traits of special interest.

A cis-regulatory signature in ascidians and flies, independent of transcription factor binding sites

BACKGROUND

Transcription initiation is controlled by cis-regulatory modules. Although these modules are usually made of clusters of short transcription factor binding sites, a small minority of such clusters in the genome have cis-regulatory activity. This paradox is currently unsolved.

RESULTS

To identify what discriminates active from inactive clusters, we focused our attention on short topologically unconstrained clusters of two ETS and two GATA binding sites, similar to the early neural enhancer of Ciona intestinalis Otx. We first computationally identified 55 such clusters, conserved between the two Ciona genomes. In vivo assay of the activity of 19 hits identified three novel early neural enhancers, all located next to genes coexpressed with Otx. Optimization of ETS and GATA binding sites was not always sufficient to confer activity to inactive clusters. Rather, a dinucleotide sequence code associated to nucleosome depletion showed a robust correlation with enhancer potential. Identification of a large collection of Ciona regulatory regions revealed that predicted nucleosome depletion constitutes a general signature of Ciona enhancers, which is conserved between orthologous loci in the two Ciona genomes and which partitions conserved noncoding sequences into a major nucleosome-bound fraction and a minor nucleosome-free fraction with higher cis-regulatory potential. We also found this signature in a large fraction of short Drosophila cis-regulatory modules.

CONCLUSION

This study indicates that a sequence-based dinucleotide signature, previously associated with nucleosome depletion and independent of transcription factor binding sites, contributes to the definition of a local cis-regulatory potential in two metazoa, Ciona intestinalis and Drosophila melanogaster.

A systems biology approach to understanding cis-regulatory module function

The genomic instructions used to regulate development are encoded within a set of functional DNA elements called cis-regulatory modules (CRMs). These elements determine the precise patterns of temporal and spatial gene expression. Here we summarize recent progress made towards cataloguing and characterizing the complete repertoire of CRMs. We describe CRMs as genomic information processing devices containing clusters of transcription factor binding sites and we position CRMs as nodes within large gene regulatory networks. We define CRM architecture and describe how these genomic elements process the information they encode to their target genes. Furthermore, we present an overview describing high-throughput techniques to identify CRMs genome wide and experimental methodologies to validate their function on a large scale. This review emphasizes the advantages and power of a systems biology approach which integrates computational and experimental technologies to further our understanding of CRM function.

Functional and comparative genomics analyses of pmp22 in medaka fish

Background

Pmp22, a member of the junction protein family Claudin/EMP/PMP22, plays an important role in myelin formation. Increase of pmp22 transcription causes peripheral neuropathy, Charcot-Marie-Tooth disease type1A (CMT1A). The pathophysiological phenotype of CMT1A is aberrant axonal myelination which induces a reduction in nerve conduction velocity (NCV). Several CMT1A model rodents have been established by overexpressing pmp22. Thus, it is thought that pmp22 expression must be tightly regulated for correct myelin formation in mammals. Interestingly, the myelin sheath is also present in other jawed vertebrates. The purpose of this study is to analyze the evolutionary conservation of the association between pmp22 transcription level and vertebrate myelin formation, and to find the conserved non-coding sequences for pmp22 regulation by comparative genomics analyses between jawed fishes and mammals.

Results

A transgenic pmp22 over-expression medaka fish line was established. The transgenic fish had approximately one fifth the peripheral NCV values of controls, and aberrant myelination of transgenic fish in the peripheral nerve system (PNS) was observed. We successfully confirmed that medaka fish pmp22 has the same exon-intron structure as mammals, and identified some known conserved regulatory motifs. Furthermore, we found novel conserved sequences in the first intron and 3'UTR.

Conclusion

Medaka fish undergo abnormalities in the PNS when pmp22 transcription increases. This result indicates that an adequate pmp22 transcription level is necessary for correct myelination of jawed vertebrates. Comparison of pmp22 orthologs between distantly related species identifies evolutionary conserved sequences that contribute to precise regulation of pmp22 expression.

Use of a novel induced spawning technique for the first reported captive spawning of Tetraodon nigroviridis

The spotted green pufferfish (Tetraodon nigroviridis) is an important genetics model animal due to its small, well-mapped genome. However, only wild-caught juveniles and adults are available to researchers. A lack of gametes, fertilized eggs, developing embryos, and other early life stages hampers development of the full potential of T. nigroviridis as a model research species. We report on successful spawning trials using a novel induced spawning technique, ovarian lavage. Chorulon(®) (human chorionic gonadotropin, hCG) was injected into a catheter inserted into the oviduct at a rate of 3 μl/g body weight. In one trial, a female paired with a male spawned in an aquarium at about 72 h post-treatment. In other trials, females were hand-stripped of eggs at 36 h post-treatment. There were 3680 eggs/g of eggs and females produced up to 24% of their body weight in eggs. Hatch resulted from all trials on the 4th day post-fertilization. Ovarian lavage is a simple method for administering spawning hormones, uses a catheter technique similar to that frequently performed to determine egg maturity in broodstock, and eliminates the need for injection.

Zebrafish embryos as models for embryotoxic and teratological effects of chemicals

The experimental virtues of the zebrafish embryo such as small size, development outside of the mother, cheap maintenance of the adult made the zebrafish an excellent model for phenotypic genetic and more recently also chemical screens. The availability of a genome sequence and several thousand mutants and transgenic lines together with gene arrays and a broad spectrum of techniques to manipulate gene functions add further to the experimental strength of this model. Pioneering studies suggest that chemicals can have in many cases very similar toxicological and teratological effects in zebrafish embryos and humans. In certain areas such as cardiotoxicity, the zebrafish appears to outplay the traditional rodent models of toxicity testing. Several pilot projects used zebrafish embryos to identify new chemical entities with specific biological functions. In combination with the establishment of transgenic sensor lines and the further development of existing and new automated imaging systems, the zebrafish embryos could therefore be used as cost-effective and ethically acceptable animal models for drug screening as well as toxicity testing.

Axonal projections originating from raphe serotonergic neurons in the developing and adult zebrafish, Danio rerio, using transgenics to visualize raphe-specific pet1 expression

Serotonin is a major central nervous modulator of physiology and behavior and plays fundamental roles during development and plasticity of the vertebrate central nervous system (CNS). Understanding the developmental control and functions of serotonergic neurons is therefore an important task. In all vertebrates, prominent serotonergic neurons are found in the superior and inferior raphe nuclei in the hindbrain innervating most CNS regions. In addition, all vertebrates except for mammals harbor other serotonergic centers, including several populations in the diencephalon. This, in combination with the intricate and wide distribution of serotonergic fibers, makes it difficult to sort out serotonergic innervation originating from the raphe from that of other serotonergic cell populations. To resolve this issue, we isolated the regulatory elements of the zebrafish raphe-specific gene pet1 and used them to drive expression of an eGFP transgene in the raphe population of serotonergic neurons. With this approach together with retrograde tracing we 1) describe in detail the development, anatomical organization, and projection pattern of zebrafish pet1-positive neurons compared with their mammalian counterparts, 2) identify a new serotonergic population in the ventrolateral zebrafish hindbrain, and 3) reveal some extent of functional subdivisions within the zebrafish superior raphe complex. Together, our results reveal for the first time the specific innervation pattern of the zebrafish raphe and, thus, provide a new model and various tools to investigate further the role of raphe serotonergic neurons in vertebrates.

Characterization of a nervous system-specific promoter for growth-associated protein 43 gene in Medaka (Oryzias latipes)

Genes expressed by neurons are controlled by specific, interacting cis-regulatory elements and trans-acting factors within their promoters. In the present study, we asked whether the transcriptional machinery regulating neuron-specific gene expression was conserved in evolution. We identified a GAP-43 homolog in Medaka (Oryzias latipes), and analyzed its expression during various stages of development. Compared with the amino acid sequences of GAP-43 homologs in other vertebrates, the amino-terminus of GAP-43 was highly conserved evolutionarily, but the carboxy-terminus exhibited significant variability. Expression of GAP-43 predominantly occurred in cells of the central and peripheral nervous systems as determined by in situ hybridization and by RT-PCR. Expression of GAP-43 increased throughout development and significant levels continued to be expressed into adulthood. We also showed that a proximal approximately 2.0 kbp fragment in the 5'-flanking region had promoter activity as determined by in vivo reporter assays. Furthermore, based upon computational analysis of transcription factor binding sites and an in vivo reporter analysis using sequentially deleted promoters, we demonstrated that cis-regulatory elements for neuronal expression were widely distributed in this region. In mammals, a TATA-box, E-box and neuronal repressive elements have been thought to contribute to neuronal expression. However, these features were not found in the orthologous region of the Medaka GAP-43 promoter. Our results suggest that the arrangement of cis-regulatory elements of the GAP-43 ortholog in Medaka is different from that in mammals, yet maintains neuron-specific regulation.

Functional and structural characterization of the zebrafish Na+-sulfate cotransporter 1 (NaS1) cDNA and gene (slc13a1)

Sulfate plays an essential role during growth, development and cellular metabolism. In this study, we characterized the function and structure of the zebrafish ( Danio rerio) Na+-sulfate cotransporter 1 (NaS1) cDNA and gene ( slc13a1). Zebrafish NaS1 encodes a protein of 583 amino acids with 13 putative transmembrane domains. Expression of zebrafish NaS1 protein in Xenopus oocytes led to Na+-sulfate cotransport, which was significantly inhibited by thiosulfate, selenate, molybdate, and tungstate. Zebrafish NaS1 transport kinetics were: Vmax = 1,731.670 ± 92.853 pmol sulfate/oocyte·hour and Km = 1.414 ± 0.275 mM for sulfate and Vmax = 307.016 ± 32.992 pmol sulfate/oocyte·hour, Km = 24.582 ± 4.547 mM and n (Hill coefficient) = 1.624 ± 0.354 for sodium. Zebrafish NaS1 mRNA is developmentally expressed in embryos from day 1 postfertilization and in the intestine, kidney, brain, and eye of adult zebrafish. The zebrafish NaS1 gene slc13a1 contains 15 exons spanning 8,716 bp. Characterization of the zebrafish NaS1 contributes to a greater understanding of sulfate transporters in a well-defined genetic model and will allow the elucidation of evolutionary and functional relationships among vertebrate sulfate transporters.

Retracted: Gene duplication and functional evolution of Hox genes in fishes

With their power to shape animal morphology, few genes have captured the imagination of biologists as much as the evolutionarily conserved members of the Hox clusters. Hox genes encode transcription factors that play a key role in specifying the body plan in metazoans and are therefore essential in explaining patterns of evolutionary diversity. While each Hox cluster contains the same genes among the different mammalian species, this does not happen in ray-finned fish, in which both the number and organization of Hox genes and even Hox clusters are variable. Teleost fishes provide the first unambiguous support for ancient whole-genome duplication (third round) in an animal lineage. The number of genes differs in each cluster as a result of increased freedom to mutate after duplication. This has also allowed them to diverge and to adopt novel developmental roles. In this review, the authors have firstly focused on broadly outlining the duplication of Hoxgenes in fishes and discussing how comparative genomics is elucidating the molecular changes associated with the evolution of Hox genes expression and developmental function in the teleost fishes.Additional related research aspects, such as imaging of roles of microRNAs, chromatin regulation and evolutionary findings are also discussed.

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.

Integrated approaches to uncovering transcription regulatory networks in mammalian cells

Integrative systems biology has emerged as an exciting research approach in molecular biology and functional genomics that involves the integration of genomics, proteomics, and metabolomics datasets. These endeavors establish a systematic paradigm by which to interrogate, model, and iteratively refine our knowledge of the regulatory events within a cell. Here we review the latest technologies available to collect high-throughput measurements of a cellular state as well as the most successful methods for the integration and interrogation of these measurements. In particular we will focus on methods available to infer transcription regulatory networks in mammals.

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.

Transposons as tools for enhancer trap screens in vertebrates

DNA transposons are efficient tools in transgenesis and have therefore become popular in the analysis of the regulatory genome in vertebrates via enhancer trap screens. Here, I discuss recent progress in this field of research, with a focus on the application of one of these transposons, namely the medaka fish derived Tol2, to enhancer trapping in zebrafish, and how this approach compares with others that have a similar objective.

Convergence of a head-field selector Otx2 and Notch signaling: a mechanism for lens specification

Xenopus is ideal for systematic decoding of cis-regulatory networks because its evolutionary position among vertebrates allows one to combine comparative genomics with efficient transgenic technology in one system. Here, we have identified and analyzed the major enhancer of FoxE3 (Lens1), a gene essential for lens formation that is activated in the presumptive lens ectoderm (PLE) when commitment to the lens fate occurs. Deletion and mutation analyses of the enhancer based on comparison of Xenopus and mammalian sequences and in vitro and in vivo binding assays identified two essential transcriptional regulators: Otx2, a homeodomain protein expressed broadly in head ectoderm including the PLE, and Su(H), a nuclear signal transducer of Notch signaling. A Notch ligand, Delta2, is expressed in the optic vesicle adjacent to the PLE, and inhibition of its activity led to loss, or severe reduction, of FoxE3 expression followed by failure of placode formation. Ectopic activation of Notch signaling induced FoxE3 expression within head ectoderm expressing Otx2, and additional misexpression of Otx2 in trunk ectoderm extended the Notch-induced FoxE3 expression posteriorly. These data provide the first direct evidence of the involvement of Notch signaling in lens induction. The obligate integration of inputs of a field-selector (Otx2) and localized signaling (Notch) within target cis-regulatory elements might be a general mechanism of organ-field specification in vertebrates (as it is in Drosophila). This concept is also consistent with classical embryological studies of many organ systems involving a ;multiple-step induction'.

A highly conserved Wnt-dependent TCF4 binding site within the proximal enhancer of the anti-myogenic Msx1 gene supports expression within Pax3-expressing limb bud muscle precursor cells

The product of the Msx1 gene is a potent inhibitor of muscle differentiation. Msx1 is expressed in muscle precursor cells of the limb bud that also express Pax3. It is thought that Msx1 may facilitate distal migration by delaying myogenesis in these cells. Despite the role played by Msx1 in inhibiting muscle differentiation, nothing is known of the mechanisms that support the expression of the Msx1 gene within limb bud muscle precursor cells. In the present study we have used a combination of comparative genomics, mouse transgenic analysis, in situ hybridisation and immunohistochemistry to identify a highly conserved and tissue-specific regulatory sub-domain within the previously characterised Msx1 gene proximal enhancer element that supports the expression of the Msx1 gene in Pax3-expressing mouse limb pre-muscle masses. Furthermore, using a combination of in situ hybridisation, in vivo ChIP assay and transgenic explant culture analysis we provide evidence that Msx1 expression in limb bud muscle precursor cells is dependent on the canonical Wnt/TCF signalling pathway that is important in muscle shape formation. The results of these studies provide evidence of a mechanistic link between the Wnt/TCF and the Msx1/Pax3/MyoD pathways within limb bud muscle precursor cells.

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

Alignment of orthologous vertebrate loci reveals that a significant proportion of conserved cis-regulatory elements have undergone shuffling during evolution.

Background

All vertebrates share a remarkable degree of similarity in their development as well as in the basic functions of their cells. Despite this, attempts at unearthing genome-wide regulatory elements conserved throughout the vertebrate lineage using BLAST-like approaches have thus far detected noncoding conservation in only a few hundred genes, mostly associated with regulation of transcription and development.

Results

We used a unique combination of tools to obtain regional global-local alignments of orthologous loci. This approach takes into account shuffling of regulatory regions that are likely to occur over evolutionary distances greater than those separating mammalian genomes. This approach revealed one order of magnitude more vertebrate conserved elements than was previously reported in over 2,000 genes, including a high number of genes found in the membrane and extracellular regions. Our analysis revealed that 72% of the elements identified have undergone shuffling. We tested the ability of the elements identified to enhance transcription in zebrafish embryos and compared their activity with a set of control fragments. We found that more than 80% of the elements tested were able to enhance transcription significantly, prevalently in a tissue-restricted manner corresponding to the expression domain of the neighboring gene.

Conclusion

Our work elucidates the importance of shuffling in the detection of cis-regulatory elements. It also elucidates how similarities across the vertebrate lineage, which go well beyond development, can be explained not only within the realm of coding genes but also in that of the sequences that ultimately govern their expression.

SCL-GFP transgenic zebrafish: in vivo imaging of blood and endothelial development and identification of the initial site of definitive hematopoiesis

The bHLH transcription factor SCL plays a central role in the generation of hematopoietic cells in vertebrates. We modified a PAC containing the whole zebrafish scl locus, inserting GFP into the first coding exon of scl. In germline-transgenic zebrafish generated using this construct, GFP expression completely recapitulates the endogenous expression of scl in blood, endothelium and CNS. We performed in vivo timelapse imaging of blood and endothelial precursor migration at the single-cell level and show that these cells migrate from the posterior lateral plate mesoderm to their site of differentiation in the intermediate cell mass. The anterior lateral plate domain of GFP expression gives rise to primitive macrophages and the blood vessels of the head. In later embryos, GFP expression identifies clusters of hematopoietic cells that develop between the dorsal aorta and posterior cardinal veins after primitive erythrocytes have entered circulation. Two treatments that block definitive hematopoiesis (treatment with dioxin (TCDD), and injection of an antisense morpholino oligonucleotide targeted to runx1) ablate these hematopoietic clusters. This indicates that these clusters represent the first site of definitive hematopoiesis in zebrafish. This site is anatomically homologous to the proposed source of hematopoietic stem cells in amniotes, the aorta-gonad-mesonephros (AGM) region. A second transgenic line, containing the promoter of scl driving GFP, lacks expression in the definitive clusters.

Cooperation of sonic hedgehog enhancers in midline expression

In zebrafish, as in other vertebrates, the secreted signalling molecule Sonic hedgehog (Shh) is expressed in organiser regions such as the embryonic midline and the zona limitans intrathalamica (zli). To investigate the regulatory mechanisms underlying the pattern of shh expression, we carried out a systematic analysis of the intronic regulatory sequences of zebrafish shh using stable transgenesis. Deletion analysis identified the modules responsible for expression in the embryonic shield, the hypothalamus and the zli and confirmed the activities of previously identified notochord and floor plate enhancers. We detected a strong synergism between regulatory regions. The degree of synergy varied over time in the hypothalamus suggesting different mechanisms for initiation and maintenance of expression. Our data show that the pattern of shh expression in the embryonic central nervous system involves an intricate crosstalk of at least 4 different regulatory regions. When compared to the enhancer activities of the mouse Shh gene, we observed a remarkable divergence of function of structurally conserved enhancer sequences. The activating region ar-C (61% identical to SFPE2 in mouse Shh), for example, mediates floor plate expression in the mouse embryo while it directs expression in the forebrain and the notochord and only weakly in the floor plate in the zebrafish embryo. This raises doubts on the predictive power of phylogenetic footprinting and indicates a stunning divergence of function of structurally conserved regulatory modules during vertebrate evolution.

Sequence analyses to study the evolutionary history and cis-regulatory elements of Hedgehog genes

Sequence analysis and comparative genomics are powerful tools to gain knowledge on multiple aspects of gene and protein regulation and function. These have been widely used to understand the evolutionary history and the biochemistry of Hedgehog (Hh) proteins, and the molecular control of Hedgehog gene expression. Here, we report on some of the methods available to retrieve protein and genomic sequences. We describe how protein sequence comparison can produce information on the evolutionary history of Hh proteins. Moreover, we describe the use of genomic sequence analysis including phylogenetic footprinting and transcription factor-binding site search tools, techniques that allow for the characterization of cis-regulatory elements of developmental genes such as the Hedgehog genes.

Comprehensive characterization of the cis-regulatory code responsible for the spatio-temporal expression of olSix3.2 in the developing medaka forebrain

BACKGROUND

Embryonic development is coordinated by sets of cis-regulatory elements that are collectively responsible for the precise spatio-temporal organization of regulatory gene networks. There is little information on how these elements, which are often associated with highly conserved noncoding sequences, are combined to generate precise gene expression patterns in vertebrates. To address this issue, we have focused on Six3, an important regulator of vertebrate forebrain development.

RESULTS

Using computational analysis and exploiting the diversity of teleost genomes, we identified a cluster of highly conserved noncoding sequences surrounding the Six3 gene. Transgenesis in medaka fish demonstrates that these sequences have enhancer, silencer, and silencer blocker activities that are differentially combined to control the entire distribution of Six3.

CONCLUSION

This report provides the first example of the precise regulatory code necessary for the expression of a vertebrate gene, and offers a unique framework for defining the interplay of trans-acting factors that control the evolutionary conserved use of Six3.

Cracking the genome's second code: enhancer detection by combined phylogenetic footprinting and transgenic fish and frog embryos

Genes involved in vertebrate development are unusually enriched for highly conserved non-coding sequence elements. These regions are readily detected in silico, by genome-wide sequence comparisons between different vertebrates, from mammals to fish (phylogenetic footprinting). It follows that sequence conservation must be the result of positive selection for an essential physiological role. An obvious possibility is that these conserved sequences possess regulatory or structural functions important for gene expression and, thus, an in vivo assay becomes necessary. We have developed a rapid testing system using zebrafish and Xenopus laevis embryos that allows us to assign transcriptional regulatory functions to conserved non-coding sequence elements. The sequences are cloned into a vector containing a minimal promoter and the GFP reporter, and are assayed for their putative cis-regulatory activity in zebrafish or Xenopus transgenic experiments. Vectors used include plasmid DNA and the Tol2 transposon system in fish and X. laevis. We have followed this logic to detect and analyze conserved elements in an intergenic region present in the Iroquois (Irx) gene clusters of zebrafish, Xenopus tropicalis, Fugu rubripes and mouse. We have assayed approximately 50 of these conserved elements and shown that the majority behave as modular positive regulatory elements (enhancers) that contribute to specific temporal and spatial domains that are part of the endogenous gene expression pattern. Moreover, comparison of the activity of cognate Irx enhancers from different organisms demonstrates that conservation of sequence is accompanied by in vivo functional conservation across species. Finally, for some of the most conserved elements, we have been able to identify a critical core sequence, essential for correct enhancer function.

New technologies, new findings, and new concepts in the study of vertebrate cis-regulatory sequences

All vertebrates share a similar early embryonic body plan and use the same regulatory genes for their development. The availability of numerous sequenced vertebrate genomes and significant advances in bioinformatics have resulted in the finding that the genomic regions of many of these developmental regulatory genes also contain highly conserved noncoding sequence. In silico discovery of conserved noncoding regions and of transcription factor binding sites as well as the development of methods for high throughput transgenesis in Xenopus and zebrafish are dramatically increasing the speed with which regulatory elements can be discovered, characterized, and tested in the context of whole live embryos. We review here some of the recent technological developments that will likely lead to a surge in research on how vertebrate genomes encode regulation of transcriptional activity, how regulatory sequences constrain genomic architecture, and ultimately how vertebrate form has evolved.

Conserved co-regulation and promoter sharing of hoxb3a and hoxb4a in zebrafish

The expression of zebrafish hoxb3a and hoxb4a has been found to be mediated through five transcripts, hoxb3a transcripts I-III and hoxb4a transcripts I-II, driven by four promoters. A "master" promoter, located about 2 kb downstream of hoxb5a, controls transcription of a pre-mRNA comprising exon sequences of both genes. This unique gene structure is proposed to provide a novel mechanism to ensure overlapping, tissue-specific expression of both genes in the posterior hindbrain and spinal cord. Transgenic approaches were used to analyze the functions of zebrafish hoxb3a/hoxb4a promoters and enhancer sequences containing regions of homology that were previously identified by comparative genomics. Two neural enhancers were shown to establish specific anterior expression borders within the hindbrain and mediate expression in defined neuronal populations derived from hindbrain rhombomeres (r) 5 to 8, suggesting a late role of the genes in neuronal cell lineage specification. Species comparison showed that the zebrafish hoxb3a r5 and r6 enhancer corresponded to a sequence within the mouse HoxA cluster controlling activity of Hoxa3 in r5 and r6, whereas a homologous region within the HoxB cluster activated Hoxb3 expression but limited to r5. We conclude that the similarity of hoxb3a/Hoxa3 regulatory mechanisms reflect the shared descent of both genes from a single ancestral paralog group 3 gene.

High-affinity peptide transporter PEPT2 (SLC15A2) of the zebrafish Danio rerio: functional properties, genomic organization, and expression analysis

Solute carrier 15 (SLC15) membrane proteins PEPT1 (SLC15A1) and PEPT2 (SLC15A2) have been described in great detail in mammals. In contrast, information in lower vertebrates is limited. We characterized the functional properties of a novel zebrafish peptide transporter orthologous to mammalian and avian PEPT2, described its gene (pept2) structure, and determined mRNA tissue distribution. An expressed sequence tag (EST) cDNA (Integrated Molecular Analysis of Gene Expression; IMAGE) corresponding to zebrafish pept2 was completed by inserting a stretch of 75 missing nucleotides in the coding sequence to obtain a 3,238-bp functional clone. The complete open reading frame (ORF) was 2,160 bp and encoded a 719-amino acid protein. Electrophysiological analysis after cRNA injection in Xenopus laevis oocytes suggested that zebrafish PEPT2 is a high-affinity/low-capacity transporter (K(0.5) for glycyl-L-glutamine approximately 18 microM at -120 mV and pH 7.5). Zebrafish pept2 gene was 19,435 kb, thus being the shortest vertebrate pept2 fully characterized so far. Also, zebrafish pept2 exhibited 23 exons and 22 introns, whereas human and rodent pept2 genes contain 22 exons and 21 introns only. Zebrafish pept2 mRNA was mainly detected in brain, kidney, gut, and, interestingly, otic vesicle, the embryonic structure that develops into the auditory/vestibular organ, homolog to the higher vertebrate inner ear, of the adult fish. Characterization of zebrafish pept2 will contribute to the investigation of peptide transporters using a well-established genetic model and will allow the elucidation of the evolutionary and functional relationships among vertebrate peptide transporters. Moreover, it can represent a useful marker to screen mutations that affect choroid plexus and inner ear development.

A functional survey of the enhancer activity of conserved non-coding sequences from vertebrate Iroquois cluster gene deserts

Recent studies of the genome architecture of vertebrates have uncovered two unforeseen aspects of its organization. First, large regions of the genome, called gene deserts, are devoid of protein-coding sequences and have no obvious biological role. Second, comparative genomics has highlighted the existence of an array of highly conserved non-coding regions (HCNRs) in all vertebrates. Most surprisingly, these structural features are strongly associated with genes that have essential functions during development. Among these, the vertebrate Iroquois (Irx) genes stand out on both fronts. Mammalian Irx genes are organized in two clusters (IrxA and IrxB) that span >1 Mb each with no other genes interspersed. Additionally, a large number of HCNRs exist within Irx clusters. We have systematically examined the enhancer activity of HCNRs from the IrxB cluster using transgenic Xenopus and zebrafish embryos. Most of these HCNRs are active in subdomains of endogenous Irx expression, and some are candidates to contain shared enhancers of neighboring genes, which could explain the evolutionary conservation of Irx clusters. Furthermore, HCNRs present in tetrapod IrxB but not in fish may be responsible for novel Irx expression domains that appeared after their divergence. Finally, we have performed a more detailed analysis on two IrxB ultraconserved non-coding regions (UCRs) duplicated in IrxA clusters in similar relative positions. These four regions share a core region highly conserved among all of them and drive expression in similar domains. However, inter-species conserved sequences surrounding the core, specific for each of these UCRs, are able to modulate their expression.

Identification of a conserved 125 base-pair Hb9 enhancer that specifies gene expression to spinal motor neurons

The homeobox gene Hb9 is expressed selectively by motor neurons (MNs) in the developing CNS. Previous studies have identified a 9-kb 5' fragment of the mouse Hb9 gene that is sufficient to direct gene expression to spinal MNs in vivo. Here, we sought to identify more discrete MN-specifying elements, using homology searches between genomic sequences of evolutionarily distant species. Based on homology screening of the mouse and human Hb9 promoters, we identified a 3.6-kb Hb9 enhancer that proved sufficient to drive MN-specific lacZ expression. We then compared mouse, human, and pufferfish (Fugu rubripes) genomic sequences, and identified a conserved 438-bp sequence, consisting of noncontiguous 313-bp and 125-bp fragments, residing within the 3.6-kb Hb9 enhancer. The zebrafish (Danio rerio) Hb9 genomic region was then found to have two identical copies of the 125-bp sequence, but no counterpart for the 313-bp sequence. Transgenic analysis showed that the 125-bp alone was both necessary and sufficient to direct spinal MN-specific lacZ expression, whereas the 313-bp sequence had no such enhancer activity. Moreover, the 125-bp Hb9 enhancer was found to harbor two Hox/Pbx consensus-binding sequences, mutations of which completely disrupted thoracolumbar Hb9 expression. These data suggest that Hox/Pbx plays a critical role in the segmental specification of spinal MNs. Together, these results indicate that the molecular pathways regulating Hb9 expression are evolutionarily conserved, and that MN-specific gene expression may be directed and achieved using a small 125-bp 5' enhancer.

Prediction of functional modules based on comparative genome analysis and Gene Ontology application

We present a computational method for the prediction of functional modules encoded in microbial genomes. In this work, we have also developed a formal measure to quantify the degree of consistency between the predicted and the known modules, and have carried out statistical significance analysis of consistency measures. We first evaluate the functional relationship between two genes from three different perspectives--phylogenetic profile analysis, gene neighborhood analysis and Gene Ontology assignments. We then combine the three different sources of information in the framework of Bayesian inference, and we use the combined information to measure the strength of gene functional relationship. Finally, we apply a threshold-based method to predict functional modules. By applying this method to Escherichia coli K12, we have predicted 185 functional modules. Our predictions are highly consistent with the previously known functional modules in E.coli. The application results have demonstrated that our approach is highly promising for the prediction of functional modules encoded in a microbial genome.

Phylogenetic footprinting and genome scanning identify vertebrate BMP response elements and new target genes

The complex gene regulatory networks governed by growth factor signaling are still poorly understood. In order to accelerate the rate of progress in uncovering these networks, we explored the usefulness of interspecies sequence comparison (phylogenetic footprinting) to identify conserved growth factor response elements. The promoter regions of two direct target genes of Bone Morphogenetic Protein (BMP) signaling in Xenopus, Xvent2 and XId3, were compared with the corresponding human and/or mouse counterparts to identify conserved sequences. A comparison between the Xenopus and human Vent2 promoter sequences revealed a highly conserved 21 bp sequence that overlaps the previously reported Xvent2 BMP response element (BRE). Reporter gene assays using Xenopus animal pole ectodermal explants (animal caps) revealed that this conserved 21 bp BRE is both necessary and sufficient for BMP responsiveness. We combine the same phylogenetic footprinting approach with luciferase assays to identify a highly conserved 49 bp BMP responsive region in the Xenopus Id3 promoter. GFP reporters containing multimers of either the Xvent2 or XId3 BREs appear to recapitulate endogenous BMP signaling activity in transgenic Xenopus embryos. Comparison of the Xvent2 and the XId3 BRE revealed core sequence features that are both necessary and sufficient for BMP responsiveness: a Smad binding element (SBE) and a GC-rich element resembling an OAZ binding site. Based on these findings, we have implemented genome scanning to identify over 100 additional putative target genes containing 2 or more BRE-like sequences which are conserved between human and mouse. RT-PCR and in situ analyses revealed that this in silico approach can effectively be used to identify potential BMP target genes.

bHLH-dependent and -independent modes of Ath5 gene regulation during retinal development

In a wide range of vertebrate species, the bHLH transcription factor Ath5 is tightly associated with both the initiation of neurogenesis in the retina and the genesis of retinal ganglion cells. Here, we describe at least two modes of regulating the expression of Ath5 during retinal development. We have found that a proximal cis-regulatory region of the Xenopus Ath5 gene (Xath5) is highly conserved across vertebrate species and is sufficient to drive retinal-specific reporter gene expression in transgenic Xenopus embryos. Xath5 proximal transgene expression depended upon two highly conserved bHLH factor binding sites (E-boxes) as well as bHLH factor activity in vivo. However, we found that bHLH activity was not required for expression of a longer Xath5 transgene, suggesting that additional mechanisms contribute to Xath5 expression in vivo. Consistent with this, we showed that a more distal fragment that does not include the conserved proximal region is sufficient to promote transgene expression in the developing retina. In mouse, we found that a longer fragment of the cis-regulatory region of either the mouse or Xenopus Ath5 gene was necessary for transgene expression, and that expression of a mouse Math5 (Atoh7) transgene was not dependent upon autoregulation. Thus, despite extensive conservation in the proximal region, the importance of these elements may be species dependent.

Comparative functional genomics revealed conservation and diversification of three enhancers of the isl1 gene for motor and sensory neuron-specific expression

Islet-1 (Isl1) is a member of the Isl1 family of LIM-homeodomain transcription factors (LIM-HD) that is expressed in a defined subset of motor and sensory neurons during vertebrate embryogenesis. To investigate how this specific expression of isl1 is regulated, we searched for enhancers of the isl1 gene that are conserved in vertebrate evolution. Initially, two enhancer elements, CREST1 and CREST2, were identified downstream of the isl1 locus in the genomes of fugu, chick, mouse, and human by BLAST searching for highly similar elements to those originally identified as motor and sensory neuron-specific enhancers in the zebrafish genome. The combined action of these elements is sufficient for completely recapitulating the subtype-specific expression of the isl1 gene in motor neurons of the mouse spinal cord. Furthermore, by direct comparison of the upstream flanking regions of the zebrafish and human isl1 genes, we identified another highly conserved noncoding element, CREST3, and subsequently C3R, a similar element to CREST3 with two CDP CR1 recognition motifs, in the upstream regions of all other isl1 family members. In mouse and human, CRESTs are located as far as more than 300 kb away from the isl1 locus, while they are much closer to the isl1 locus in zebrafish. Although all of zebrafish CREST2, CREST3, and C3R activate gene expression in the sensory neurons of zebrafish, CREST2 of mouse and human does not have the sequence necessary for sensory neuron-specific expression. Our results revealed both a remarkable conservation of the regulatory elements regulating subtype-specific gene expression in motor and sensory neurons and the dynamic process of reorganization of these elements whereby each element increases the level of cell-type specificity by losing redundant functions with the other elements during vertebrate evolution.

Comparative genomics using fugu: a tool for the identification of conserved vertebrate cis-regulatory elements

With the imminent completion of the whole genome sequence of humans, increasing attention is being focused on the annotation of cis-regulatory elements in the human genome. Comparative genomics approaches based on evolutionary conservation have proved useful in the detection of conserved cis-regulatory elements. The pufferfish, Fugu rubripes, is an attractive vertebrate model for comparative genomics, by virtue of its compact genome and maximal phylogenetic distance from mammals. Fugu has lost a large proportion of nonessential DNA, and retained single orthologs for many duplicate genes that arose in the fish lineage. Non-coding sequences conserved between fugu and mammals have been shown to be functional cis-regulatory elements. Thus, fugu is a model fish genome of choice for discovering evolutionarily conserved regulatory elements in the human genome. Such evolutionarily conserved elements are likely to be shared by all vertebrates, and related to regulatory interactions fundamental to all vertebrates. The functions of these conserved vertebrate elements can be rapidly assayed in mammalian cell lines or in transgenic systems such as zebrafish/medaka and Xenopus, followed by validation of crucial elements in transgenic rodents.

Monorail/Foxa2 regulates floorplate differentiation and specification of oligodendrocytes, serotonergic raphé neurones and cranial motoneurones

In this study, we elucidate the roles of the winged-helix transcription factor Foxa2 in ventral CNS development in zebrafish. Through cloning of monorail (mol), which we find encodes the transcription factor Foxa2, and phenotypic analysis of mol-/- embryos, we show that floorplate is induced in the absence of Foxa2 function but fails to further differentiate. In mol-/- mutants, expression of Foxa and Hh family genes is not maintained in floorplate cells and lateral expansion of the floorplate fails to occur. Our results suggest that this is due to defects both in the regulation of Hh activity in medial floorplate cells as well as cell-autonomous requirements for Foxa2 in the prospective laterally positioned floorplate cells themselves. Foxa2 is also required for induction and/or patterning of several distinct cell types in the ventral CNS. Serotonergic neurones of the raphenucleus and the trochlear motor nucleus are absent in mol-/- embryos, and oculomotor and facial motoneurones ectopically occupy ventral CNS midline positions in the midbrain and hindbrain. There is also a severe reduction of prospective oligodendrocytes in the midbrain and hindbrain. Finally, in the absence of Foxa2, at least two likely Hh pathway target genes are ectopically expressed in more dorsal regions of the midbrain and hindbrain ventricular neuroepithelium, raising the possibility that Foxa2 activity may normally be required to limit the range of action of secreted Hh proteins.

Phylogenetic footprinting reveals extensive conservation of Sonic Hedgehog (SHH) regulatory elements

Sonic Hedgehog (SHH) plays a fundamental role in numerous developmental processes including morphogenesis of limbs, nervous system, and teeth. Using a Bayesian alignment algorithm for phylogenetic footprinting we analyzed approximately 28 kb of noncoding DNA in the SHH locus of human and mouse. This showed that the length of conserved noncoding sequences (4196 nt) shared by these species was approximately 3 times larger than the SHH coding sequence (1386 nt). Most segments were located in introns (53%) or within 2-kb regions upstream (16%) or downstream (20%) of the first and last SHH codon. Even though regions more than 2 kb upstream or downstream of the first and last SHH codon represented 57% (16 kb) of the sequence compared, they accounted for only 11% (494 nt) of the total length of conserved noncoding segments. One region of 650 nt downstream of SHH was identified as a putative scaffold/matrix attachment region (SMAR). Human-mouse analysis was complemented by sequencing in apes, monkeys, rodents, and bats, thus further confirming the evolutionary conservation of some segments. Gel-shift assays indicated that conserved segments are targeted by nuclear proteins and showed differences between two cell types that expressed different levels of SHH, namely human endothelial cells and breast cancer cells. The relevance of these findings with respect to regulation of SHH expression during normal and pathologic development is discussed.

Conserved and acquired features of neurogenin1 regulation

The telencephalon shows vast morphological variations among different vertebrate groups. The transcription factor neurogenin1(ngn1) controls neurogenesis in the mouse pallium and is also expressed in the dorsal telencephalon of the evolutionary distant zebrafish. The upstream regions of the zebrafish and mammalian ngn1 loci harbour several stretches of conserved sequences. Here, we show that the upstream region of zebrafish ngn1 is capable of faithfully recapitulating endogenous expression in the zebrafish and mouse telencephalon. A single conserved regulatory region is essential for dorsal telencephalic expression in the zebrafish, and for expression in the dorsal pallium of the mouse. However, a second conserved region that is inactive in the fish telencephalon is necessary for expression in the lateral pallium of mouse embryos. This regulatory region, which drives expression in the zebrafish diencephalon and hindbrain, is dependent on Pax6 activity and binds recombinant Pax6 in vitro. Thus, the regulatory elements of ngn1 appear to be conserved among vertebrates, with certain differences being incorporated in the utilisation of these enhancers, for the acquisition of more advanced features in amniotes. Our data provide evidence for the co-option of regulatory regions as a mechanism of evolutionary diversification of expression patterns, and suggest that an alteration in Pax6expression was crucial in neocortex evolution.

Surveying phylogenetic footprints in large gene clusters: applications to Hox cluster duplications

Evolutionarily conserved non-coding genomic sequences represent a potentially rich source for the discovery of gene regulatory regions. Since these elements are subject to stabilizing selection they evolve much more slowly than adjacent non-functional DNA. These so-called phylogenetic footprints can be detected by comparison of the sequences surrounding orthologous genes in different species. Therefore the loss of phylogenetic footprints as well as the acquisition of conserved non-coding sequences in some lineages, but not in others, can provide evidence for the evolutionary modification of cis-regulatory elements. We introduce here a statistical model of footprint evolution that allows us to estimate the loss of sequence conservation that can be attributed to gene loss and other structural reasons. This approach to studying the pattern of cis-regulatory element evolution, however, requires the comparison of relatively long sequences from many species. We have therefore developed an efficient software tool for the identification of corresponding footprints in long sequences from multiple species. We apply this novel method to the published sequences of HoxA clusters of shark, human, and the duplicated zebrafish and Takifugu clusters as well as the published HoxB cluster sequences. We find that there is a massive loss of sequence conservation in the intergenic region of the HoxA clusters, consistent with the finding in [Chiu et al., PNAS 99 (2002) 5492]. The loss of conservation after cluster duplication is more extensive than expected from structural reasons. This suggests that binding site turnover and/or adaptive modification may also contribute to the loss of sequence conservation.

Phylogenetic analysis of 5'-noncoding regions from the ABA-responsive rab16/17 gene family of sorghum, maize and rice provides insight into the composition, organization and function of cis-regulatory modules

Phylogenetic analysis of sequences from gene families and homologous genes from species of varying divergence can be used to identify conserved noncoding regulatory elements. In this study, phylogenetic analysis of 5'-noncoding sequences was optimized using rab17, a well-characterized ABA-responsive gene from maize, and five additional rab16/17 homologs from sorghum and rice. Conserved 5'-noncoding sequences among the maize, sorghum, and rice rab16/17 homologs were identified with the aid of the software program FootPrinter and by screening for known transcription-factor-binding sites. Searches for 7 of 8 (7/8)bp sequence matches within aligned 5'-noncoding segments of the rab genes identified many of the cis-elements previously characterized by biochemical analysis in maize rab17 plus several additional putative regulatory elements. Differences in the composition of conserved noncoding sequences among rab16/17 genes were related to variation in rab gene mRNA levels in different tissues and to response to ABA treatment using qRT-PCR. Absence of a GRA-like element in the promoter of sorghum dhn2 relative to maize rab17 was correlated with an approximately 85-fold reduction of dhn2 RNA in sorghum shoots. Overall, we conclude that phylogenetic analysis of gene families among rice, sorghum, and maize will help identify regulatory sequences in the noncoding regions of genes and contribute to our understanding of grass gene regulatory networks.

Regulation of Otx2 expression and its functions in mouse forebrain and midbrain

Otx2 expression in the forebrain and midbrain was found to be regulated by two distinct enhancers (FM and FM2) located at 75 kb 5′upstream and 115 kb 3′ downstream. The activities of these two enhancers were absent in anterior neuroectoderm earlier than E8.0; however, at E9.5 their regions of activity spanned the entire mesencephalon and diencephalon with their caudal limits at the boundary with the metencephalon or isthmus. In telencephalon, activities were found only in the dorsomedial aspect. Potential binding sites of OTX and TCF were essential to FM activity, and TCF sites were also essential to FM2 activity. The FM2 enhancer appears to be unique to rodent; however, the FM enhancer region is deeply conserved in gnathostomes. Studies of mutants lacking FM or FM2 enhancer demonstrated that these enhancers indeed regulate Otx2 expression in forebrain and midbrain. Development of mesencephalic and diencephalic regions was differentially regulated in a dose-dependent manner by the cooperation between Otx1and Otx2 under FM and FM2 enhancers: the more caudal the structure the higher the OTX dose requirement. At E10.5 Otx1–/–Otx2ΔFM/ΔFMmutants, in which Otx2 expression under the FM2 enhancer remained,exhibited almost complete loss of the entire diencephalon and mesencephalon;the telencephalon did, however, develop.