Lens development depends on a pair of highly conserved Sox21 regulatory elements

Highly conserved non-coding elements (CNEs) linked to genes involved in embryonic development have been hypothesised to correspond to cis-regulatory modules due to their ability to induce tissue-specific expression patterns. However, attempts to prove their requirement for normal development or for the correct expression of the genes they are associated with have yielded conflicting results. Here, we show that CNEs at the vertebrate Sox21 locus are crucial for Sox21 expression in the embryonic lens and that loss of Sox21 function interferes with normal lens development. Using different expression assays in zebrafish we find that two CNEs linked to Sox21 in all vertebrates contain lens enhancers and that their removal from a reporter BAC abolishes lens expression. Furthermore inhibition of Sox21 function after the injection of a sox21b morpholino into zebrafish leads to defects in lens development. These findings identify a direct link between sequence conservation and genomic function of regulatory sequences. In addition to this we provide evidence that putative Sox binding sites in one of the CNEs are essential for induction of lens expression as well as enhancer function in the CNS. Our results show that CNEs identified in pufferfish-mammal whole-genome comparisons are crucial developmental enhancers and hence essential components of gene regulatory networks underlying vertebrate embryogenesis.     

Early Evolution of Conserved Regulatory Sequences Associated with Development in Vertebrates

Comparisons between diverse vertebrate genomes have uncovered thousands of highly conserved non-coding sequences, an increasing number of which have been shown to function as enhancers during early development. Despite their extreme conservation over 500 million years from humans to cartilaginous fish, these elements appear to be largely absent in invertebrates, and, to date, there has been little understanding of their mode of action or the evolutionary processes that have modelled them. We have now exploited emerging genomic sequence data for the sea lamprey, Petromyzon marinus, to explore the depth of conservation of this type of element in the earliest diverging extant vertebrate lineage, the jawless fish (agnathans). We searched for conserved non-coding elements (CNEs) at 13 human gene loci and identified lamprey elements associated with all but two of these gene regions. Although markedly shorter and less well conserved than within jawed vertebrates, identified lamprey CNEs are able to drive specific patterns of expression in zebrafish embryos, which are almost identical to those driven by the equivalent human elements. These CNEs are therefore a unique and defining characteristic of all vertebrates. Furthermore, alignment of lamprey and other vertebrate CNEs should permit the identification of persistent sequence signatures that are responsible for common patterns of expression and contribute to the elucidation of the regulatory language in CNEs. Identifying the core regulatory code for development, common to all vertebrates, provides a foundation upon which regulatory networks can be constructed and might also illuminate how large conserved regulatory sequence blocks evolve and become fixed in genomic DNA.     

Functional analysis of conserved non-coding regions around the short stature hox gene (shox) in whole zebrafish embryos

BackgroundMutations in the SHOX gene are responsible for Leri-Weill Dyschondrosteosis, a disorder characterised by mesomelic limb shortening. Recent investigations into regulatory elements surrounding SHOX have shown that deletions of conserved non-coding elements (CNEs) downstream of the SHOX gene produce a phenotype indistinguishable from Leri-Weill Dyschondrosteosis. As this gene is not found in rodents, we used zebrafish as a model to characterise the expression pattern of the shox gene across the whole embryo and characterise the enhancer domains of different CNEs associated with this gene.Conclusion/SignificanceOur results using whole zebrafish embryos have provided a more comprehensive picture of the expression pattern of the shox gene, and a better understanding of its regulation via deeply conserved noncoding elements. In particular, we identify additional tissues under the regulatory control of previously identified SHOX CNEs. We also demonstrate the importance of these CNEs in evolution by identifying duplicated shox CNEs and more deeply conserved sub-sequences within already identified CNEs.     

Ancient vertebrate conserved noncoding elements have been evolving rapidly in teleost fishes

Vertebrate genomes contain thousands of conserved noncoding elements (CNEs) that often function as tissue-specific enhancers. In this study, we have identified CNEs in human, dog, chicken, Xenopus, and four teleost fishes (zebrafish, stickleback, medaka, and fugu) using elephant shark, a cartilaginous vertebrate, as the base genome and investigated the evolution of these ancient vertebrate CNEs (aCNEs) in bony vertebrate lineages. Our analysis shows that aCNEs have been evolving at different rates in different bony vertebrate lineages. Although 78-83% of CNEs have diverged beyond recognition ("lost") in different teleost fishes, only 24% and 40% have been lost in the chicken and mammalian lineages, respectively. Relative rate tests of substitution rates in CNEs revealed that the teleost fish CNEs have been evolving at a significantly higher rate than those in other bony vertebrates. In the ray-finned fish lineage, 68% of aCNEs were lost before the divergence of the four teleosts. This implicates the "fish-specific" whole-genome duplication in the accelerated evolution and the loss of a large number of both copies of duplicated CNEs in teleost fishes. The aCNEs are rich in tissue-specific enhancers and thus many of them are likely to be evolutionarily constrained cis-regulatory elements. The rapid evolution of aCNEs might have affected the expression patterns driven by them. Transgenic zebrafish assay of some human CNE enhancers that have been lost in teleosts has indicated instances of conservation or changes in trans-acting factors between mammals and fishes.     

Highly conserved regulatory elements around the SHH gene may contribute to the maintenance of conserved synteny across human chromosome 7q36.3

Comparative genomic analysis reveals an exceptionally large section of conserved shared synteny between the human 7q36 chromosomal region and the pufferfish (Fugu rubripes) genome. Remarkably, this conservation extends not only to gene order across 16 genes, but also to the position and orientation of a number of prominent conserved noncoding elements (CNEs). A functional assay using zebrafish has shown that most of the CNEs have reproducible and specific enhancer activity. This enhancer activity is often detected in a subset of tissues which reflect the endogenous expression pattern of a proximal gene, though some CNEs may act over a long range. We propose that the distribution of CNEs, and their probable association with a number of genes throughout the region, imposes a critical constraint on genome architecture, resulting in the maintenance of such a large section of conserved synteny across the vertebrate lineage.     

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.     

Function and Regulation of AUTS2, a Gene Implicated in Autism and Human Evolution

Nucleotide changes in the AUTS2 locus, some of which affect only noncoding regions, are associated with autism and other neurological disorders, including attention deficit hyperactivity disorder, epilepsy, dyslexia, motor delay, language delay, visual impairment, microcephaly, and alcohol consumption. In addition, AUTS2 contains the most significantly accelerated genomic region differentiating humans from Neanderthals, which is primarily composed of noncoding variants. However, the function and regulation of this gene remain largely unknown. To characterize auts2 function, we knocked it down in zebrafish, leading to a smaller head size, neuronal reduction, and decreased mobility. To characterize AUTS2 regulatory elements, we tested sequences for enhancer activity in zebrafish and mice. We identified 23 functional zebrafish enhancers, 10 of which were active in the brain. Our mouse enhancer assays characterized three mouse brain enhancers that overlap an ASD–associated deletion and four mouse enhancers that reside in regions implicated in human evolution, two of which are active in the brain. Combined, our results show that AUTS2 is important for neurodevelopment and expose candidate enhancer sequences in which nucleotide variation could lead to neurological disease and human-specific traits.     

Cis-regulatory characterization of sequence conservation surrounding the Hox4 genes

Hox genes are key regulators of anterior-posterior axis patterning and have a major role in hindbrain development. The zebrafish Hox4 paralogs have strong overlapping activities in hindbrain rhombomeres 7 and 8, in the spinal cord and in the pharyngeal arches. With the aim to predict enhancers that act on the hoxa4a, hoxb4a, hoxc4a and hoxd4a genes, we used sequence conservation around the Hox4 genes to analyze all fish:human conserved non-coding sequences by reporter assays in stable zebrafish transgenesis. Thirty-four elements were functionally tested in GFP reporter gene constructs and more than 100 F1 lines were analyzed to establish a correlation between sequence conservation and cis-regulatory function, constituting a catalog of Hox4 CNEs. Sixteen tissue-specific enhancers could be identified. Multiple alignments of the CNEs revealed paralogous cis-regulatory sequences, however, the CNE sequence similarities were found not to correlate with tissue specificity. To identify ancestral enhancers that direct Hox4 gene activity, genome sequence alignments of mammals, teleosts, horn shark and the cephalochordate amphioxus, which is the most basal extant chordate possessing a single prototypical Hox cluster, were performed. Three elements were identified and two of them exhibited regulatory activity in transgenic zebrafish, however revealing no specificity. Our data show that the approach to identify cis-regulatory sequences by genome sequence alignments and subsequent testing in zebrafish transgenesis can be used to define enhancers within the Hox clusters and that these have significantly diverged in their function during evolution.     

Otx2 expression in anterior neuroectoderm and forebrain/midbrain is directed by more than six enhancers

Otx2 plays essential roles in each site at each step of head development. We previously identified the AN1 enhancer at 91 kb 5' upstream for the Otx2 expressions in anterior neuroectoderm (AN) at neural plate stage before E8.5, and the FM1 enhancer at 75 kb 5' upstream and the FM2 enhancer at 122 kb 3' downstream for the expression in forebrain/midbrain (FM) at brain vesicle stage after E8.5. The present study identified a second AN enhancer (AN2) at 88 kb 5' upstream; the AN2 enhancer also recapitulates the endogenous Otx2 expression in choroid plexus, cortical hem and choroidal roof. However, the enhancer mutants indicated the presence of another AN enhancer. The study also identified a third FM enhancer (FM3) at 153 kb 5' upstream. Thus, the Otx2 expressions in anterior neuroectoderm and forebrain/midbrain are regulated by more than six enhancers located far from the coding region. The enhancers identified are differentially conserved among vertebrates; none of the AN enhancers has activities in caudal forebrain and midbrain at brain vesicle stage after E8.5, nor do any of the FM enhancers in anterior neuroectoderm at neural plate stage before E8.5.     

Three enhancer regions regulate gbx2 gene expression in the isthmic region during zebrafish development

In vertebrate embryos, positioning of the boundary between the midbrain and hindbrain (MHB) and subsequent isthmus formation are dependent upon the interaction between the Otx2 and Gbx genes. In zebrafish, sequential expression of gbx1 and gbx2 in the anterior hindbrain contributes to this process, whereas in mouse embryos, a single Gbx gene (Gbx2) is responsible for MHB development. In the present study, to investigate the regulatory mechanism of gbx2 in the MHB/isthmic region of zebrafish embryos, we cloned the gene and showed that its organization is conserved among different vertebrates. Promoter analyses revealed three enhancers that direct reporter gene expression after the end of epiboly in the anterior-most hindbrain, which is a feature of the zebrafish gbx2 gene. One of the enhancers is located upstream of gbx2 (AMH1), while the other two enhancers are located downstream of gbx2 (AMH2 and AMH3). Detailed analysis of the AMH1 enhancer showed that it directs expression in the rhombomere 1 (r1) region and the dorsal thalamus, as has been shown for gbx2, whereas no expression was induced by the AMH1 enhancer in other embryonic regions in which gbx2 is expressed. The AMH1 enhancer is composed of multiple regulatory subregions that share the same spatial specificity. The most active of the regulatory subregions is a 291-bp region that contains at least two Pax2-binding sites, both of which are necessary for the function of the main component (PB1-A region) of the AMH1 enhancer. In accordance with these results, enhancer activity in the PB1-A region, as well as gbx2 expression in r1, was missing in no isthmus mutant embryos that lacked functional pax2a. In addition, we identified an upstream conserved sequence of 227bp that suppresses the enhancer activity of AMH1. Taken together, these findings suggest that gbx2 expression during the somitogenesis stage in zebrafish is regulated by a complex mechanism involving Pax2 as well as activators and suppressors in the regions flanking the gene.     

Cloning and functional analysis of the proximal promoter region of the three GnRH genes from the silver sea bream (Sparus sarba)

Gonadotropin-releasing hormone (GnRH) is a neuropeptide that plays a major role in releasing pituitary gonadotropin and controlling vertebrate reproduction. In this study, three GnRH cDNAs, GnRH-I (sbGnRH; 348 bp), GnRH-II (cGnRH-II; 557 bp), and GnRH-III (sGnRH; 483 bp), were cloned from the brain of the silver sea bream (Sparus sarba). In order to understand how the expression of the GnRH isoforms was regulated in the brain, the promoter of each gene was cloned and analyzed. We found regulatory motifs in the promoters that were conserved in the GnRH promoters of tilapia and zebrafish, suggesting that these motifs play a critical role in GnRH regulation. We performed functional analyses and examined tissue-specific expression for each GnRH promoter using EGFP reporter fusions in zebrafish. The GnRH-I promoter was active in the forebrain area, including the olfactory bulb-terminal nerve area and peripheral preoptic areas; the GnRH-II promoter was active in the midbrain; and the GnRH-III promoter was active in the olfactory bulb. These results show that the GnRH promoters of the silver sea bream GnRH genes exhibit tissue-specific activity.     

Functional analysis of limb enhancers in the developing fin

Despite diverging ～365 million years ago, tetrapod limbs and pectoral fins express similar genes that could be regulated by shared regulatory elements. In this study, we set out to analyze the ability of enhancers to maintain tissue specificity in these two divergent structures. We tested 22 human sequences that were previously reported as mouse limb enhancers for their enhancer activity in zebrafish (Danio rerio). Using a zebrafish enhancer assay, we found that 10/22 (45 %) were positive for pectoral fin activity. Analysis of the various criteria that correlated with positive fin activity found that both spatial limb activity and evolutionary conservation are not good predictors of fin enhancer activity. These results suggest that zebrafish enhancer assays may be limited in detecting human limb enhancers, and this limitation does not improve by the use of limb spatial expression or evolutionary conservation.