Cis-regulation and chromosomal rearrangement of the fgf8 locus after the teleost/tetrapod split

The importance of considering regulatory domains in genome-wide analyses - the nearest gene is often wrong!

The expression of a large number of genes is regulated by regulatory elements that are located far away from their promoters. Identifying which gene is the target of a specific regulatory element or is affected by a non-coding mutation is often accomplished by assigning these regions to the nearest gene in the genome. However, this heuristic ignores key features of genome organisation and gene regulation; in that the genome is partitioned into regulatory domains, which at some loci directly coincide with the span of topologically associated domains (TADs), and that genes are regulated by enhancers located throughout these regions, even across intervening genes. In this review, we examine the results from genome-wide studies using chromosome conformation capture technologies and from those dissecting individual gene regulatory domains, to highlight that the phenomenon of enhancer skipping is pervasive and affects multiple types of genes. We discuss how simply assigning a genomic region of interest to its nearest gene is problematic and often leads to incorrect predictions and highlight that where possible information on both the conservation and topological organisation of the genome should be used to generate better hypotheses.

The article has an associated Future Leader to Watch interview.

Summary: Identifying which gene is the target of an enhancer is often accomplished by assigning it to the nearest gene, here we discuss how this heuristic can lead to incorrect predictions.

The bowfin genome illuminates the developmental evolution of ray-finned fishes

The bowfin (Amia calva) is a ray-finned fish that possesses a unique suite of ancestral and derived phenotypes, which are key to understanding vertebrate evolution. The phylogenetic position of bowfin as a representative of neopterygian fishes, its archetypical body plan and its unduplicated and slowly evolving genome make bowfin a central species for the genomic exploration of ray-finned fishes. Here we present a chromosome-level genome assembly for bowfin that enables gene-order analyses, settling long-debated neopterygian phylogenetic relationships. We examine chromatin accessibility and gene expression through bowfin development to investigate the evolution of immune, scale, respiratory and fin skeletal systems and identify hundreds of gene-regulatory loci conserved across vertebrates. These resources connect developmental evolution among bony fishes, further highlighting the bowfin's importance for illuminating vertebrate biology and diversity in the genomic era.

Dissection of the Fgf8 regulatory landscape by in vivo CRISPR-editing reveals extensive intra- and inter-enhancer redundancy

Developmental genes are often regulated by multiple elements with overlapping activity. Yet, in most cases, the relative function of those elements and their contribution to endogenous gene expression remain poorly characterized. An example of this phenomenon is that distinct sets of enhancers have been proposed to direct Fgf8 in the limb apical ectodermal ridge and the midbrain-hindbrain boundary. Using in vivo CRISPR/Cas9 genome engineering, we functionally dissect this complex regulatory ensemble and demonstrate two distinct regulatory logics. In the apical ectodermal ridge, the control of Fgf8 expression appears distributed between different enhancers. In contrast, we find that in the midbrain-hindbrain boundary, one of the three active enhancers is essential while the other two are dispensable. We further dissect the essential midbrain-hindbrain boundary enhancer to reveal that it is also composed by a mixture of essential and dispensable modules. Cross-species transgenic analysis of this enhancer suggests that its composition may have changed in the vertebrate lineage.

Comparative Analysis of the Organization of the Catecholaminergic Systems in the Brain of Holostean Fishes (Actinopterygii/Neopterygii)

Living holosteans, comprising 8 species of bowfins and gars, form a small monophyletic group of actinopterygian fishes, which are currently considered as the sister group to the enormously numerous teleosts and have largely been neglected in neuroanatomical studies. We have studied the catecholaminergic (CAergic) systems by means of antibodies against tyrosine hydroxylase (TH) and dopamine (DA) in the brain of representative species of the 3 genera included in the 2 orders of holostean fishes: Amia calva (Amiiformes) and Lepisosteus platyrhincus, Lepisosteus oculatus, and Atractosteus spatula (Lepisosteiformes). Different groups of TH/DA-immunoreactive (ir) cells were observed in the olfactory bulb, subpallium, and preoptic area of the telencephalon. Hypothalamic groups were labeled in the suprachiasmatic nucleus, tuberal (only in A. calva), retrotuberal, and retromamillary areas; specifically, the paraventricular organ showed only DA immunoreactivity. In the diencephalon, TH/DA-ir groups were detected in the prethalamus, posterior tubercle, and pretectum. In the caudal hindbrain, the solitary tract nucleus and area postrema presented TH/DA-ir cell groups, and also the spinal cord and the retina. Only in A. calva, particular CAergic cell groups were observed in the habenula, the mesencephalic tegmentum, and in the locus coeruleus. Following a neuromeric analysis, the comparison of these results with those obtained in other classes of fishes and tetrapods shows many common traits of CAergic systems shared by most vertebrates and in addition highlights unique features of actinopterygian fishes.

Common genetic causes of holoprosencephaly are limited to a small set of evolutionarily conserved driver genes of midline development coordinated by TGF-β, hedgehog, and FGF signaling

Here, we applied targeted capture to examine 153 genes representative of all the major vertebrate developmental pathways among 333 probands to rank their relative significance as causes for holoprosencephaly (HPE). We now show that comparisons of variant transmission versus nontransmission among 136 HPE Trios indicates some reported genes now lack confirmation, while novel genes are implicated. Furthermore, we demonstrate that variation of modest intrinsic effect can synergize with these driver mutations as gene modifiers.

Robust Identification of Developmentally Active Endothelial Enhancers in Zebrafish Using FANS-Assisted ATAC-Seq

Identification of tissue-specific and developmentally active enhancers provides insights into mechanisms that control gene expression during embryogenesis. However, robust detection of these regulatory elements remains challenging, especially in vertebrate genomes. Here, we apply fluorescent-activated nuclei sorting (FANS) followed by Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) to identify developmentally active endothelial enhancers in the zebrafish genome. ATAC-seq of nuclei from Tg(fli1a:egfp)^y1 transgenic embryos revealed expected patterns of nucleosomal positioning at transcriptional start sites throughout the genome and association with active histone modifications. Comparison of ATAC-seq from GFP-positive and -negative nuclei identified more than 5,000 open elements specific to endothelial cells. These elements flanked genes functionally important for vascular development and that displayed endothelial-specific gene expression. Importantly, a majority of tested elements drove endothelial gene expression in zebrafish embryos. Thus, FANS-assisted ATAC-seq using transgenic zebrafish embryos provides a robust approach for genome-wide identification of active tissue-specific enhancer elements.

Copy-number variants and candidate gene mutations in isolated split hand/foot malformation

Split hand/foot malformation (SHFM) is a congenital limb deficiency with missing or shortened central digits. Some SHFM genes have been identified but the cause of many SHFM cases is unknown. We used single-nucleotide polymorphism (SNP) microarray analysis to detect copy-number variants (CNVs) in 25 SHFM cases without other birth defects from New York State (NYS), prioritized CNVs absent from population CNV databases, and validated these CNVs using quantitative real-time polymerase chain reaction (qPCR). We tested for the validated CNVs in seven cases from Iowa using qPCR, and also sequenced 36 SHFM candidate genes in all the subjects. Seven NYS cases had a potentially deleterious variant: two had a p.R225H or p.R225L mutation in TP63, one had a 17q25 microdeletion, one had a 10q24 microduplication and three had a 17p13.3 microduplication. In addition, one Iowa case had a de novo 10q24 microduplication. The 17q25 microdeletion has not been reported previously in SHFM and included two SHFM candidate genes (SUMO2 and GRB2), while the 10q24 and 17p13.3 CNVs had breakpoints within genomic regions that contained putative regulatory elements and a limb development gene. In SHFM pathogenesis, the microdeletion may cause haploinsufficiency of SHFM genes and/or deletion of their regulatory regions, and the microduplications could disrupt regulatory elements that control transcription of limb development genes.

Regulatory evolution of Tbx5 and the origin of paired appendages

The diversification of paired appendages has been a major factor in the evolutionary radiation of vertebrates. Despite its importance, an understanding of the origin of paired appendages has remained elusive. To address this problem, we focused on T-box transcription factor 5 (Tbx5), a gene indispensable for pectoral appendage initiation and development. Comparison of gene expression in jawless and jawed vertebrates reveals that the Tbx5 expression in jawed vertebrates is derived in having an expression domain that extends caudal to the heart and gills. Chromatin profiling, phylogenetic footprinting, and functional assays enabled the identification of a Tbx5 fin enhancer associated with this apomorphic pattern of expression. Comparative functional analysis of reporter constructs reveals that this enhancer activity is evolutionarily conserved among jawed vertebrates and is able to rescue the finless phenotype of tbx5a mutant zebrafish. Taking paleontological evidence of early vertebrates into account, our results suggest that the gain of apomorphic patterns of Tbx5 expression and regulation likely contributed to the morphological transition from a finless to finned condition at the base of the vertebrate lineage.

Genetic basis of hindlimb loss in a naturally occurring vertebrate model

Here we genetically characterise pelvic finless, a naturally occurring model of hindlimb loss in zebrafish that lacks pelvic fin structures, which are homologous to tetrapod hindlimbs, but displays no other abnormalities. Using a hybrid positional cloning and next generation sequencing approach, we identified mutations in the nuclear localisation signal (NLS) of T-box transcription factor 4 (Tbx4) that impair nuclear localisation of the protein, resulting in altered gene expression patterns during pelvic fin development and the failure of pelvic fin development. Using a TALEN-induced tbx4 knockout allele we confirm that mutations within the Tbx4 NLS (A78V; G79A) are sufficient to disrupt pelvic fin development. By combining histological, genetic, and cellular approaches we show that the hindlimb initiation gene tbx4 has an evolutionarily conserved, essential role in pelvic fin development. In addition, our novel viable model of hindlimb deficiency is likely to facilitate the elucidation of the detailed molecular mechanisms through which Tbx4 functions during pelvic fin and hindlimb development.

Summary: Here we genetically characterise mutations in tbx4 which underlie pelvic finless, a naturally occurring model of hindlimb loss in zebrafish that lacks pelvic fin structures.

Cis-regulatory programs in the development and evolution of vertebrate paired appendages

Differential gene expression is the core of development, mediating the genetic changes necessary for determining cell identity. The regulation of gene activity by cis-acting elements (e.g., enhancers) is a crucial mechanism for determining differential gene activity by precise control of gene expression in embryonic space and time. Modifications to regulatory regions can have profound impacts on phenotype, and therefore developmental and evolutionary biologists have increasingly focused on elucidating the transcriptional control of genes that build and pattern body plans. Here, we trace the evolutionary history of transcriptional control of three loci key to vertebrate appendage development (Fgf8, Shh, and HoxD/A). Within and across these regulatory modules, we find both complex and flexible regulation in contrast with more fixed enhancers that appear unchanged over vast timescales of vertebrate evolution. The transcriptional control of vertebrate appendage development was likely already incredibly complex in the common ancestor of fish, implying that subtle changes to regulatory networks were more likely responsible for alterations in phenotype rather than the de novo addition of whole regulatory domains. Finally, we discuss the dangers of relying on inter-species transgenesis when testing enhancer function, and call for more controlled regulatory swap experiments when inferring the evolutionary history of enhancer elements.

Diversity in cell motility reveals the dynamic nature of the formation of zebrafish taste sensory organs

Taste buds are sensory organs in jawed vertebrates, composed of distinct cell types that detect and transduce specific taste qualities. Taste bud cells differentiate from oropharyngeal epithelial progenitors localized mainly in proximity of the forming organs. Despite recent progress in elucidating the molecular interactions required for taste bud cell development and function, the cell behaviour underlying the organ assembly is poorly defined. Here, we used time-lapse imaging to observe the formation of taste buds in live zebrafish larvae. We found that tg(fgf8a.dr17) expressing cells form taste buds and get rearranged within the forming organs. In addition, differentiating cells move from the epithelium to the forming organs and can be displaced between developing organs. During organ formation, taste bud tg(fgf8a.dr17) and Type-II cells are displaced in random, directed or confined mode relative to the taste bud they join or are maintained. Finally, ascl1a activity in the 5-HT/Type-III cell is required to direct and maintain tg(fgf8a.dr17) expressing cells into the taste bud. We propose diversity in displacement modes of differentiating cells as a key mechanism for the highly dynamic process of taste bud assembly.

Identification of Critical Region Responsible for Split Hand/Foot Malformation Type 3 (SHFM3) Phenotype through Systematic Review of Literature and Mapping of Breakpoints Using Microarray Data

Split hand/foot malformation (SHFM) is a limb malformation with underdeveloped or absent central digital rays, clefts of hands and feet, and variable syndactyly of the remaining digits. There are six types of SHFM. Here, we report a boy with SHFM type 3 having normal 4th and 5th digits, absent 2nd and 3rd digits, and a 4th finger flexion deformity, as well as absent 2nd, 3rd and 4th toes bilaterally. His father, two paternal uncles, and two paternal first cousins have similar phenotype. Chromosome analysis showed a normal male karyotype. A 514 kb gain at 10q24.31-q24.32 (chr10:102,962,134-103,476,346, hg19) was identified using 6.0 Single nucleotide polymorphism (SNP) microarray, resulting in the duplication of nine genes, including BTRC and FBXW4. A detailed systematic review of literature and mapping of breakpoints using microarray data from all reported cases in PubMed and DECIPHER were conducted, and exon 1 of BTRC gene was identified as the critical region responsible for the SHFM3 phenotype. The potential mechanism and future studies of this critical region causing the SHFM3 phenotype are discussed.

Copy number variants in patients with intellectual disability affect the regulation of ARX transcription factor gene

Protein-coding mutations in the transcription factor-encoding gene ARX cause various forms of intellectual disability (ID) and epilepsy. In contrast, variations in surrounding non-coding sequences are correlated with milder forms of non-syndromic ID and autism and had suggested the importance of ARX gene regulation in the etiology of these disorders. We compile data on several novel and some already identified patients with or without ID that carry duplications of ARX genomic region and consider likely genetic mechanisms underlying the neurodevelopmental defects. We establish the long-range regulatory domain of ARX and identify its brain region-specific autoregulation. We conclude that neurodevelopmental disturbances in the patients may not simply arise from increased dosage due to ARX duplication. This is further exemplified by a small duplication involving a non-functional ARX copy, but with duplicated enhancers. ARX enhancers are located within a 504-kb region and regulate expression specifically in the forebrain in developing and adult zebrafish. Transgenic enhancer-reporter lines were used as in vivo tools to delineate a brain region-specific negative and positive autoregulation of ARX. We find autorepression of ARX in the telencephalon and autoactivation in the ventral thalamus. Fluorescently labeled brain regions in the transgenic lines facilitated the identification of neuronal outgrowth and pathfinding disturbances in the ventral thalamus and telencephalon that occur when arxa dosage is diminished. In summary, we have established a model for how breakpoints in long-range gene regulation alter the expression levels of a target gene brain region-specifically, and how this can cause subtle neuronal phenotypes relating to the etiology of associated neuropsychiatric disease.

The architecture of gene expression: integrating dispersed cis-regulatory modules into coherent regulatory domains

Specificity and precision of expression are essential for the genes that regulate developmental processes. The specialized cis-acting modules, such as enhancers, that define gene expression patterns can be distributed across large regions, raising questions about the nature of the mechanisms that underline their action. Recent data has exposed the structural 3D context in which these long-range enhancers are operating. Here, we present how these studies shed new light on principles driving long-distance regulatory relationships. We discuss the molecular mechanisms that enable and accompany the action of long-range acting elements and the integration of multiple distributed regulatory inputs into the coherent and specific regulatory programs that are key to embryonic development.

Leucophores are similar to xanthophores in their specification and differentiation processes in medaka

Animal body color is generated primarily by neural crest-derived pigment cells in the skin. Mammals and birds have only melanocytes on the surface of their bodies; however, fish have a variety of pigment cell types or chromatophores, including melanophores, xanthophores, and iridophores. The medaka has a unique chromatophore type called the leucophore. The genetic basis of chromatophore diversity remains poorly understood. Here, we report that three loci in medaka, namely, leucophore free (lf), lf-2, and white leucophore (wl), which affect leucophore and xanthophore differentiation, encode solute carrier family 2, member 15b (slc2a15b), paired box gene 7a (pax7a), and solute carrier family 2 facilitated glucose transporter, member 11b (slc2a11b), respectively. Because lf-2, a loss-of-function mutant for pax7a, causes defects in the formation of xanthophore and leucophore precursor cells, pax7a is critical for the development of the chromatophores. This genetic evidence implies that leucophores are similar to xanthophores, although it was previously thought that leucophores were related to iridophores, as these chromatophores have purine-dependent light reflection. Our identification of slc2a15b and slc2a11b as genes critical for the differentiation of leucophores and xanthophores in medaka led to a further finding that the existence of these two genes in the genome coincides with the presence of xanthophores in nonmammalian vertebrates: birds have yellow-pigmented irises with xanthophore-like intracellular organelles. Our findings provide clues for revealing diverse evolutionary mechanisms of pigment cell formation in animals.

Foxi transcription factors promote pharyngeal arch development by regulating formation of FGF signaling centers

The bones of the vertebrate face develop from transient embryonic branchial arches that are populated by cranial neural crest cells. We have characterized a mouse mutant for the Forkhead family transcription factor Foxi3, which is expressed in branchial ectoderm and endoderm. Foxi3 mutant mice are not viable and display severe branchial arch-derived facial skeleton defects, including absence of all but the most distal tip of the mandible and complete absence of the inner, middle and external ear structures. Although cranial neural crest cells of Foxi3 mutants are able to migrate, populate the branchial arches, and display some elements of correct proximo-distal patterning, they succumb to apoptosis from embryonic day 9.75 onwards. We show this cell death correlates with a delay in expression of Fgf8 in branchial arch ectoderm and a failure of neural crest cells in the arches to express FGF-responsive genes. Zebrafish foxi1 is also expressed in branchial arch ectoderm and endoderm, and morpholino knock-down of foxi1 also causes apoptosis of neural crest in the branchial arches. We show that heat shock induction of fgf3 in zebrafish arch tissue can rescue cell death in foxi1 morphants. Our results suggest that Foxi3 may play a role in the establishment of signaling centers in the branchial arches that are required for neural crest survival, patterning and the subsequent development of branchial arch derivatives.

Deep conservation of cis-regulatory elements in metazoans

Despite the vast morphological variation observed across phyla, animals share multiple basic developmental processes orchestrated by a common ancestral gene toolkit. These genes interact with each other building complex gene regulatory networks (GRNs), which are encoded in the genome by cis-regulatory elements (CREs) that serve as computational units of the network. Although GRN subcircuits involved in ancient developmental processes are expected to be at least partially conserved, identification of CREs that are conserved across phyla has remained elusive. Here, we review recent studies that revealed such deeply conserved CREs do exist, discuss the difficulties associated with their identification and describe new approaches that will facilitate this search.

The mystery of extreme non-coding conservation

Regions of several dozen to several hundred base pairs of extreme conservation have been found in non-coding regions in all metazoan genomes. The distribution of these elements within and across genomes has suggested that many have roles as transcriptional regulatory elements in multi-cellular organization, differentiation and development. Currently, there is no known mechanism or function that would account for this level of conservation at the observed evolutionary distances. Previous studies have found that, while these regions are under strong purifying selection, and not mutational coldspots, deletion of entire regions in mice does not necessarily lead to identifiable changes in phenotype during development. These opposing findings lead to several questions regarding their functional importance and why they are under strong selection in the first place. In this perspective, we discuss the methods and techniques used in identifying and dissecting these regions, their observed patterns of conservation, and review the current hypotheses on their functional significance.

Using zebrafish transgenesis to test human genomic sequences for specific enhancer activity

We detail an approach for the identification of human tissue-specific transcriptional enhancers involving three steps: delineation of search space around a locus or target gene, in silico identification and size definition of putative candidate sequences, and testing through several independent genomic insertions in a transgenic zebrafish reporter assay. Candidate sequences are defined through evolutionary conservation, transcription factor binding and chromatin marks (e.g. ENCODE data) and are amplified from genomic DNA, cloned into basal promoter:fluorescent protein reporter vectors based on the Tol2 transposon system and are microinjected into fertilized zebrafish eggs. After raising injected founders to sexual maturity, fluorescent screening identifies positive founder fish whose offspring undergo a detailed expression analysis to determine tissue specificity and reproducibility of specific enhancers.

An integrated holo-enhancer unit defines tissue and gene specificity of the Fgf8 regulatory landscape

Fgf8 encodes a key signaling factor, and its precise regulation is essential for embryo patterning. Here, we identified the regulatory modules that control Fgf8 expression during mammalian embryogenesis. These enhancers are interspersed with unrelated genes along a large region of 220 kb; yet they act on Fgf8 only. Intriguingly, this region also contains additional genuine enhancer activities that are not transformed into gene expression. Using genomic engineering strategies, we showed that these multiple and distinct regulatory modules act as a coherent unit and influence genes depending on their position rather than on their promoter sequence. These findings highlight how the structure of a locus regulates the autonomous intrinsic activities of the regulatory elements it contains and contributes to their tissue and target specificities. We discuss the implications of such regulatory systems regarding the evolution of gene expression and the impact of human genomic structural variations.

Differences in enhancer activity in mouse and zebrafish reporter assays are often associated with changes in gene expression

Background

Phenotypic evolution in animals is thought to be driven in large part by differences in gene expression patterns, which can result from sequence changes in cis-regulatory elements (cis-changes) or from changes in the expression pattern or function of transcription factors (trans-changes). While isolated examples of trans-changes have been identified, the scale of their overall contribution to regulatory and phenotypic evolution remains unclear.

Results

Here, we attempt to examine the prevalence of trans-effects and their potential impact on gene expression patterns in vertebrate evolution by comparing the function of identical human tissue-specific enhancer sequences in two highly divergent vertebrate model systems, mouse and zebrafish. Among 47 human conserved non-coding elements (CNEs) tested in transgenic mouse embryos and in stable zebrafish lines, at least one species-specific expression domain was observed in the majority (83%) of cases, and 36% presented dramatically different expression patterns between the two species. Although some of these discrepancies may be due to the use of different transgenesis systems in mouse and zebrafish, in some instances we found an association between differences in enhancer activity and changes in the endogenous gene expression patterns between mouse and zebrafish, suggesting a potential role for trans-changes in the evolution of gene expression.

Conclusions

In total, our results: (i) serve as a cautionary tale for studies investigating the role of human enhancers in different model organisms, and (ii) suggest that changes in the trans environment may play a significant role in the evolution of gene expression in vertebrates.

Zebrafish as a genomics model for human neurological and polygenic disorders

Whole exome sequencing and, to a lesser extent, genome-wide association studies, have provided unprecedented advances in identifying genes and candidate genomic regions involved in the development of human disease. Further progress will come from sequencing the entire genome of multiple patients and normal controls to evaluate overall mutational burden and disease risk. A major challenge will be the interpretation of the resulting data and distinguishing true pathogenic mutations from rare benign variants.While in model organisms such as the zebrafish,mutants are sought that disrupt the function of individual genes, human mutations that cause, or are associated with, the development of disease, are often not acting in a Mendelian fashion, are frequently of small effect size, are late onset, and may reside in noncoding parts of the genome. The zebrafish model is uniquely poised for understanding human coding- and noncoding variants because of its sequenced genome, a large body of knowledge on gene expression and function, rapid generation time, and easy access to embryos. A critical advantage is the ease of zebrafish transgenesis, both for the testing of human regulatory DNA driving expression of fluorescent reporter proteins, and the expression of mutated disease-associated human proteins in specific neurons to rapidly model aspects of neurological disorders. The zebrafish affords progress both through its model genome and it is rapidly developing transparent model vertebrate embryo.

A SINE-derived element constitutes a unique modular enhancer for mammalian diencephalic Fgf8

Transposable elements, including short interspersed repetitive elements (SINEs), comprise nearly half the mammalian genome. Moreover, they are a major source of conserved non-coding elements (CNEs), which play important functional roles in regulating development-related genes, such as enhancing and silencing, serving for the diversification of morphological and physiological features among species. We previously reported a novel SINE family, AmnSINE1, as part of mammalian-specific CNEs. One AmnSINE1 locus, named AS071, showed an enhancer property in the developing mouse diencephalon. Indeed, AS071 appears to recapitulate the expression of diencephalic fibroblast growth factor 8 (Fgf8). Here we established three independent lines of AS071-transgenic mice and performed detailed expression profiling of AS071-enhanced lacZ in comparison with that of Fgf8 across embryonic stages. We demonstrate that AS071 is a distal enhancer that directs Fgf8 expression in the developing diencephalon. Furthermore, enhancer assays with constructs encoding partially deleted AS071 sequence revealed a unique modular organization in which AS071 contains at least three functionally distinct sub-elements that cooperatively direct the enhancer activity in three diencephalic domains, namely the dorsal midline and the lateral wall of the diencephalon, and the ventral midline of the hypothalamus. Interestingly, the AmnSINE1-derived sub-element was found to specify the enhancer activity to the ventral midline of the hypothalamus. To our knowledge, this is the first discovery of an enhancer element that could be separated into respective sub-elements that determine regional specificity and/or the core enhancing activity. These results potentiate our understanding of the evolution of retroposon-derived cis-regulatory elements as well as the basis for future studies of the molecular mechanism underlying the determination of domain-specificity of an enhancer.

An ancient genomic regulatory block conserved across bilaterians and its dismantling in tetrapods by retrogene replacement

Developmental genes are regulated by complex, distantly located cis-regulatory modules (CRMs), often forming genomic regulatory blocks (GRBs) that are conserved among vertebrates and among insects. We have investigated GRBs associated with Iroquois homeobox genes in 39 metazoans. Despite 600 million years of independent evolution, Iroquois genes are linked to ankyrin-repeat-containing Sowah genes in nearly all studied bilaterians. We show that Iroquois-specific CRMs populate the Sowah locus, suggesting that regulatory constraints underlie the maintenance of the Iroquois-Sowah syntenic block. Surprisingly, tetrapod Sowah orthologs are intronless and not associated with Iroquois; however, teleost and elephant shark data demonstrate that this is a derived feature, and that many Iroquois-CRMs were ancestrally located within Sowah introns. Retroposition, gene, and genome duplication have allowed selective elimination of Sowah exons from the Iroquois regulatory landscape while keeping associated CRMs, resulting in large associated gene deserts. These results highlight the importance of CRMs in imposing constraints to genome architecture, even across large phylogenetic distances, and of gene duplication-mediated genetic redundancy to disentangle these constraints, increasing genomic plasticity.

Analysis of genome-wide association study data using the protein knowledge base

BACKGROUND

Genome-wide association studies (GWAS) aim to identify causal variants and genes for complex disease by independently testing a large number of SNP markers for disease association. Although genes have been implicated in these studies, few utilise the multiple-hit model of complex disease to identify causal candidates. A major benefit of multi-locus comparison is that it compensates for some shortcomings of current statistical analyses that test the frequency of each SNP in isolation for the phenotype population versus control.

RESULTS

Here we developed and benchmarked several protocols for GWAS data analysis using different in-silico gene prediction and prioritisation methodologies. We adopted a high sensitivity approach to the data, using less conservative statistical SNP associations. Multiple gene search spaces, either of fixed-widths or proximity-based, were generated around each SNP marker. We used the candidate disease gene prediction system Gentrepid to identify candidates based on shared biomolecular pathways or domain-based protein homology. Predictions were made either with phenotype-specific known disease genes as input; or without a priori knowledge, by exhaustive comparison of genes in distinct loci. Because Gentrepid uses biomolecular data to find interactions and common features between genes in distinct loci of the search spaces, it takes advantage of the multi-locus aspect of the data.

CONCLUSIONS

Results suggest testing multiple SNP-to-gene search spaces compensates for differences in phenotypes, populations and SNP platforms. Surprisingly, domain-based homology information was more informative when benchmarked against gene candidates reported by GWA studies compared to previously determined disease genes, possibly suggesting a larger contribution of gene homologs to complex diseases than Mendelian diseases.

Dissecting the transcriptional regulatory properties of human chromosome 16 highly conserved non-coding regions

Non-coding DNA conservation across species has been often used as a predictor for transcriptional enhancer activity. However, only a few systematic analyses of the function of these highly conserved non-coding regions (HCNRs) have been performed. Here we use zebrafish transgenic assays to perform a systematic study of 113 HCNRs from human chromosome 16. By comparing transient and stable transgenesis, we show that the first method is highly inefficient, leading to 40% of false positives and 20% of false negatives. When analyzed in stable transgenic lines, a great majority of HCNRs were active in the central nervous system, although some of them drove expression in other organs such as the eye and the excretory system. Finally, by testing a fraction of the HCNRs lacking enhancer activity for in vivo insulator activity, we find that 20% of them may contain enhancer-blocking function. Altogether our data indicate that HCNRs may contain different types of cis-regulatory activity, including enhancer, insulators as well as other not yet discovered functions.

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.

Identification of cis regulatory features in the embryonic zebrafish genome through large-scale profiling of H3K4me1 and H3K4me3 binding sites

An organism's genome sequence serves as a blueprint for the proteins and regulatory RNAs essential for cellular function. The genome also harbors cis-acting non-coding sequences that control gene expression and are essential to coordinate regulatory programs during embryonic development. However, the genome sequence is largely identical between cell types within a multi-cellular organism indicating that factors such as DNA accessibility and chromatin structure play a crucial role in governing cell-specific gene expression. Recent studies have identified particular chromatin modifications that define functionally distinct cis regulatory elements. Among these are forms of histone 3 that are mono- or tri-methylated at lysine 4 (H3K4me1 or H3K4me3, respectively), which bind preferentially to promoter and enhancer elements in the mammalian genome. In this work, we investigated whether these modified histones could similarly identify cis regulatory elements within the zebrafish genome. By applying chromatin immunoprecipitation followed by deep sequencing, we find that H3K4me1 and H3K4me3 are enriched at transcriptional start sites in the genome of the developing zebrafish embryo and that this association correlates with gene expression. We further find that these modifications associate with distal non-coding conserved elements, including known active enhancers. Finally, we demonstrate that it is possible to utilize H3K4me1 and H3K4me3 binding profiles in combination with available expression data to computationally identify relevant cis regulatory sequences flanking syn-expressed genes in the developing embryo. Taken together, our results indicate that H3K4me1 and H3K4me3 generally mark cis regulatory elements within the zebrafish genome and indicate that further characterization of the zebrafish using this approach will prove valuable in defining transcriptional networks in this model system.

Zebrafish: an integrative system for neurogenomics and neurosciences

Rapid technological advances over the past decade have moved us closer to a high throughput molecular approach to neurobiology, where we see the merging of neurogenetics, genomics, physiology, imaging and pharmacology. This is the case more in zebrafish than in any other model organism commonly used. Recent improvements in the generation of transgenic zebrafish now allow genetic manipulation and live imaging of neuronal development and function in early embryonic, larval, and adult animals. The sequenced zebrafish genome and comparative genomics give unprecedented insights into genome evolution and its relation to genome structure and function. There is now information on embryonic and larval expression of over 12,000 genes and just under 1000 mutant phenotypes. We review the remarkable similarity of the zebrafish genetic blueprint for the nervous system to that of mammals and assess recent technological advances that make the zebrafish a model of choice for elucidating the development and function of neuronal circuitry, transgene-based neuroanatomy, and small molecule neuropharmacology.

Midline signaling and evolution of the forebrain in chordates: a focus on the lamprey Hedgehog case

Lampreys are agnathans (vertebrates without jaws). They occupy a key phylogenetic position in the emergence of novelties and in the diversification of morphology at the dawn of vertebrates. We have used lampreys to investigate the possibility that embryonic midline signaling systems have been a driving force for the evolution of the forebrain in vertebrates. We have focused on Sonic Hedgehog/Hedgehog (Shh/Hh) signaling. In this article, we first review and summarize our recent work on the comparative analysis of embryonic expression patterns for Shh/Hh, together with Fgf8 (fibroblast growth factor 8) and Wnt (wingless-Int) pathway components, in the embryonic lamprey forebrain. Comparison with nonvertebrate chordates on one hand, and jawed vertebrates on the other hand, shows that these morphogens/growth factors acquired new expression domains in the most rostral part of the neural tube in lampreys compared to nonvertebrate chordates, and in jawed vertebrates compared to lampreys. These data are consistent with the idea that changes in Shh, Fgf8 or Wnt signaling in the course of evolution have been instrumental for the emergence and diversification of the telencephalon, a part of the forebrain that is unique to vertebrates. We have then used comparative genomics on Shh/Hh loci to identify commonalities and differences in noncoding regulatory sequences across species and phyla. Conserved noncoding elements (CNEs) can be detected in lamprey Hh introns, even though they display unique structural features and need adjustments of parameters used for in silico alignments to be detected, because of lamprey-specific properties of the genome. The data also show conservation of a ventral midline enhancer located in Shh/Hh intron 2 of all chordates, the very species which possess a notochord and a floor plate, but not in earlier emerged deuterostomes or protostomes. These findings exemplify how the Shh/Hh locus is one of the best loci to study genome evolution with regards to developmental events.

A Hypomorphic Allele in the FGF8 Gene Contributes to Holoprosencephaly and Is Allelic to Gonadotropin-Releasing Hormone Deficiency in Humans

Holoprosencephaly (HPE), the most common malformation of the human forebrain, may arise due to interacting genetic and environmental factors. To date, at least 12 contributory genes have been identified. Fibroblast growth factor 8 (Fgf8) belongs to the FGF family of genes expressed in several developmental signaling centers, including the anterior neural ridge, which is implicated in midline anomalies in mice. In humans, FGF8 mutations have been previously reported in facial clefting and in hypogonadotropic hypogonadism, but have not been reported in patients with HPE. We screened 360 probands with HPE for sequence variations in FGF8 using High Resolution DNA Melting (HRM) and sequenced all identified variations. Here we describe a total of 8 sequence variations in HPE patients, including a putative loss-of-function mutation in 3 members of a family with variable forms of classic HPE, and relate these findings to the phenotypes seen in other conditions.

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.