Ancora: a web resource for exploring highly conserved noncoding elements and their association with developmental regulatory genes

Metazoan genomes contain arrays of highly conserved noncoding elements (HCNEs) that span developmental regulatory genes and define regulatory domains. We describe Ancora , a web resource that provides data and tools for exploring genomic organization of HCNEs for multiple genomes. Ancora includes a genome browser that shows HCNE locations and features novel HCNE density plots as a powerful tool to discover developmental regulatory genes and distinguish their regulatory elements and domains.     

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.     

Highly Conserved Non-Coding Sequences and the 18q Critical Region for Short Stature: A Common Mechanism of Disease?

Background

Isolated growth hormone deficiency (IGHD) and multiple pituitary hormone deficiency (MPHD) are heterogeneous disorders with several different etiologies and they are responsible for most cases of short stature. Mutations in different genes have been identified but still many patients did not present mutations in any of the known genes. Chromosomal rearrangements may also be involved in short stature and, among others, deletions of 18q23 defined a critical region for the disorder. No gene was yet identified.

Methodology/Principal Findings

We now report a balanced translocation X;18 in a patient presenting a breakpoint in 18q23 that was surprisingly mapped about 500 Kb distal from the short stature critical region. It separated from the flanking SALL3 gene a region enriched in highly conserved non-coding elements (HCNE) that appeared to be regulatory sequences, active as enhancers or silencers during embryonic development.

Conclusion

We propose that, during pituitary development, the 18q rearrangement may alter expression of 18q genes or of X chromosome genes that are translocated next to the HCNEs. Alteration of expression of developmentally regulated genes by translocation of HCNEs may represent a common mechanism for disorders associated to isolated chromosomal rearrangements.     

A systematic approach to identify functional motifs within vertebrate developmental enhancers

Uncovering the cis-regulatory logic of developmental enhancers is critical to understanding the role of non-coding DNA in development. However, it is cumbersome to identify functional motifs within enhancers, and thus few vertebrate enhancers have their core functional motifs revealed. Here we report a combined experimental and computational approach for discovering regulatory motifs in developmental enhancers. Making use of the zebrafish gene expression database, we computationally identified conserved non-coding elements (CNEs) likely to have a desired tissue-specificity based on the expression of nearby genes. Through a high throughput and robust enhancer assay, we tested the activity of ∼ 100 such CNEs and efficiently uncovered developmental enhancers with desired spatial and temporal expression patterns in the zebrafish brain. Application of de novo motif prediction algorithms on a group of forebrain enhancers identified five top-ranked motifs, all of which were experimentally validated as critical for forebrain enhancer activity. These results demonstrate a systematic approach to discover important regulatory motifs in vertebrate developmental enhancers. Moreover, this dataset provides a useful resource for further dissection of vertebrate brain development and function.     

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.     

Shadow enhancers flanking the HoxB cluster direct dynamic Hox expression in early heart and endoderm development

The products of Hox genes function in assigning positional identity along the anterior–posterior body axis during animal development. In mouse embryos, Hox genes located at the 3′ end of HoxA and HoxB complexes are expressed in nested patterns in the progenitors of the secondary heart field during early cardiogenesis and the combined activities of both of these clusters are required for proper looping of the heart. Using Hox bacterial artificial chromosomes (BACs), transposon reporters, and transgenic analyses in mice, we present the identification of several novel enhancers flanking the HoxB complex which can work over a long range to mediate dynamic reporter expression in the endoderm and embryonic heart during development. These enhancers respond to exogenously added retinoic acid and we have identified two retinoic acid response elements (RAREs) within these control modules that play a role in potentiating their regulatory activity. Deletion analysis in HoxB BAC reporters reveals that these control modules, spread throughout the flanking intergenic region, have regulatory activities that overlap with other local enhancers. This suggests that they function as shadow enhancers to modulate the expression of genes from the HoxB complex during cardiac development. Regulatory analysis of the HoxA complex reveals that it also has enhancers in the 3′ flanking region which contain RAREs and have the potential to modulate expression in endoderm and heart tissues. Together, the similarities in their location, enhancer output, and dependence on retinoid signaling suggest that a conserved cis-regulatory cassette located in the 3′ proximal regions adjacent to the HoxA and HoxB complexes evolved to modulate Hox gene expression during mammalian cardiac and endoderm development. This suggests a common regulatory mechanism, whereby the conserved control modules act over a long range on multiple Hox genes to generate nested patterns of HoxA and HoxB expression during cardiogenesis.     

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.     

Identification of a dopaminergic enhancer indicates complexity in vertebrate dopamine neuron phenotype specification

The dopaminergic neurons of the basal ganglia play critical roles in CNS function and human disease, but specification of dopamine neuron phenotype is poorly understood in vertebrates. We performed an in vivo screen in zebrafish to identify dopaminergic neuron enhancers, in order to facilitate studies on the specification of neuronal identity, connectivity, and function in the basal ganglia. Based primarily on identification of conserved non-coding elements, we tested 54 DNA elements from four species (zebrafish, pufferfish, mouse, and rat), that included 21 genes with known or putative roles in dopaminergic neuron specification or function. Most elements failed to drive CNS expression or did not express specifically in dopaminergic neurons. However, we did isolate a discrete enhancer from the otpb gene that drove specific expression in diencephalic dopaminergic neurons, although it did not share sequence conservation with regulatory regions of otpa or other dopamine-specific genes. For the otpb enhancer, regulation of expression in dopamine neurons requires multiple elements spread across a large genomic area. In addition, we compared our in vivo testing with in silico analysis of genomic regions for genes involved in dopamine neuron function, but failed to find conserved regions that functioned as enhancers. We conclude that regulation of dopaminergic neuron phenotype in vertebrates is regulated by dispersed regulatory elements.     

cis-Regulatory Complexity within a Large Non-Coding Region in the Drosophila Genome

Analysis of cis-regulatory enhancers has revealed that they consist of clustered blocks of highly conserved sequences. Although most characterized enhancers reside near their target genes, a growing number of studies have shown that enhancers located over 50 kb from their minimal promoter(s) are required for appropriate gene expression and many of these ‘long-range’ enhancers are found in genomic regions that are devoid of identified exons. To gain insight into the complexity of Drosophila cis-regulatory sequences within exon-poor regions, we have undertaken an evolutionary analysis of 39 of these regions located throughout the genome. This survey revealed that within these genomic expanses, clusters of conserved sequence blocks (CSBs) are positioned once every 1.1 kb, on average, and that a typical cluster contains multiple (5 to 30 or more) CSBs that have been maintained for at least 190 My of evolutionary divergence. As an initial step toward assessing the cis-regulatory activity of conserved clusters within gene-free genomic expanses, we have tested the in-vivo enhancer activity of 19 consecutive CSB clusters located in the middle of a 115 kb gene-poor region on the 3^rd chromosome. Our studies revealed that each cluster functions independently as a specific spatial/temporal enhancer. In total, the enhancers possess a diversity of regulatory functions, including dynamically activating expression in defined patterns within subsets of cells in discrete regions of the embryo, larvae and/or adult. We also observed that many of the enhancers are multifunctional–that is, they activate expression during multiple developmental stages. By extending these results to the rest of the Drosophila genome, which contains over 70,000 non-coding CSB clusters, we suggest that most function as enhancers.