Systematic human/zebrafish comparative identification of cis-regulatory activity around vertebrate developmental transcription factor genes

Genomic insights into the seawater adaptation in Cyprinidae

BACKGROUND

Cyprinidae, the largest fish family, encompasses approximately 367 genera and 3006 species. While they exhibit remarkable adaptability to diverse aquatic environments, it is exceptionally rare to find them in seawater, with the Far Eastern daces being of few exceptions. Therefore, the Far Eastern daces serve as a valuable model for studying the genetic mechanisms underlying seawater adaptation in Cyprinidae.

RESULTS

Here, we sequenced the chromosome-level genomes of two Far Eastern daces (Pseudaspius brandtii and P. hakonensis), the two known cyprinid fishes found in seawater, and performed comparative genomic analyses to investigate their genetic mechanism of seawater adaptation. Demographic history reconstruction of the two species reveals that their population dynamics are correlated with the glacial-interglacial cycles and sea level changes. Genomic analyses identified Pseudaspius-specific genetic innovations related to seawater adaptation, including positively selected genes, rapidly evolving genes, and conserved non-coding elements (CNEs). Functional assays of Pseudaspius-specific variants of the prolactin (prl) gene showed enhanced cell adaptation to greater osmolarity. Functional assays of Pseudaspius specific CNEs near atg7 and usp45 genes suggest that they exhibit higher promoter activity and significantly induced at high osmolarity.

CONCLUSIONS

Our results reveal the genome-wide evidence for the evolutionary adaptation of cyprinid fishes to seawater, offering valuable insights into the molecular mechanisms supporting the survival of migratory fish in marine environments. These findings are significant as they contribute to our understanding of how cyprinid fishes navigate and thrive in diverse aquatic habitats, providing useful implications for the conservation and management of marine ecosystems.

The Density of Regulatory Information Is a Major Determinant of Evolutionary Constraint on Noncoding DNA in Drosophila

Evolutionary analyses have estimated that ∼60% of nucleotides in intergenic regions of the Drosophila melanogaster genome are functionally relevant, suggesting that regulatory information may be encoded more densely in intergenic regions than has been revealed by most functional dissections of regulatory DNA. Here, we approached this issue through a functional dissection of the regulatory region of the gene shavenbaby (svb). Most of the ∼90 kb of this large regulatory region is highly conserved in the genus Drosophila, though characterized enhancers occupy a small fraction of this region. By analyzing the regulation of svb in different contexts of Drosophila development, we found that the regulatory information that drives svb expression in the abdominal pupal epidermis is organized in a different way than the elements that drive svb expression in the embryonic epidermis. While in the embryonic epidermis svb is activated by compact enhancers separated by large inactive DNA regions, svb expression in the pupal epidermis is driven by regulatory information distributed over broader regions of svb cis-regulatory DNA. In the same vein, we observed that other developmental genes also display a dense distribution of putative regulatory elements in their regulatory regions. Furthermore, we found that a large percentage of conserved noncoding DNA of the Drosophila genome is contained within regions of open chromatin. These results suggest that part of the evolutionary constraint on noncoding DNA of Drosophila is explained by the density of regulatory information, which may be greater than previously appreciated.

Genetic predisposition to neuroblastoma results from a regulatory polymorphism that promotes the adrenergic cell state

Childhood neuroblastomas exhibit plasticity between an undifferentiated neural crest-like mesenchymal cell state and a more differentiated sympathetic adrenergic cell state. These cell states are governed by autoregulatory transcriptional loops called core regulatory circuitries (CRCs), which drive the early development of sympathetic neuronal progenitors from migratory neural crest cells during embryogenesis. The adrenergic cell identity of neuroblastoma requires LMO1 as a transcriptional cofactor. Both LMO1 expression levels and the risk of developing neuroblastoma in children are associated with a single nucleotide polymorphism, G/T, that affects a GATA motif in the first intron of LMO1. Here, we showed that WT zebrafish with the GATA genotype developed adrenergic neuroblastoma, while knock-in of the protective TATA allele at this locus reduced the penetrance of MYCN-driven tumors, which were restricted to the mesenchymal cell state. Whole genome sequencing of childhood neuroblastomas demonstrated that TATA/TATA tumors also exhibited a mesenchymal cell state and were low risk at diagnosis. Thus, conversion of the regulatory GATA to a TATA allele in the first intron of LMO1 reduced the neuroblastoma-initiation rate by preventing formation of the adrenergic cell state. This mechanism was conserved over 400 million years of evolution, separating zebrafish and humans.

Identification of clade-wide putative cis-regulatory elements from conserved non-coding sequences in Cucurbitaceae genomes

Cis-regulatory elements regulate gene expression and play an essential role in the development and physiology of organisms. Many conserved non-coding sequences (CNSs) function as cis-regulatory elements. They control the development of various lineages. However, predicting clade-wide cis-regulatory elements across several closely related species remains challenging. Based on the relationship between CNSs and cis-regulatory elements, we present a computational approach that predicts the clade-wide putative cis-regulatory elements in 12 Cucurbitaceae genomes. Using 12-way whole-genome alignment, we first obtained 632 112 CNSs in Cucurbitaceae. Next, we identified 16 552 Cucurbitaceae-wide cis-regulatory elements based on collinearity among all 12 Cucurbitaceae plants. Furthermore, we predicted 3 271 potential regulatory pairs in the cucumber genome, of which 98 were verified using integrative RNA sequencing and ChIP sequencing datasets from samples collected during various fruit development stages. The CNSs, Cucurbitaceae-wide cis-regulatory elements, and their target genes are accessible at http://cmb.bnu.edu.cn/cisRCNEs_cucurbit/. These elements are valuable resources for functionally annotating CNSs and their regulatory roles in Cucurbitaceae genomes.

Genetic Predisposition to Neuroblastoma Results from a Regulatory Polymorphism that Promotes the Adrenergic Cell State

Childhood neuroblastomas exhibit plasticity between an undifferentiated neural crest-like "mesenchymal" cell state and a more differentiated sympathetic "adrenergic" cell state. These cell states are governed by autoregulatory transcriptional loops called core regulatory circuitries (CRCs), which drive the early development of sympathetic neuronal progenitors from migratory neural crest cells during embryogenesis. The adrenergic cell identity of neuroblastoma requires LMO1 as a transcriptional co-factor. Both LMO1 expression levels and the risk of developing neuroblastoma in children are associated with a single nucleotide polymorphism G/T that affects a G ATA motif in the first intron of LMO1. Here we show that wild-type zebrafish with the G ATA genotype develop adrenergic neuroblastoma, while knock-in of the protective T ATA allele at this locus reduces the penetrance of MYCN-driven tumors, which are restricted to the mesenchymal cell state. Whole genome sequencing of childhood neuroblastomas demonstrates that T ATA/ T ATA tumors also exhibit a mesenchymal cell state and are low risk at diagnosis. Thus, conversion of the regulatory G ATA to a T ATA allele in the first intron of LMO1 reduces the neuroblastoma initiation rate by preventing formation of the adrenergic cell state, a mechanism that is conserved over 400 million years of evolution separating zebrafish and humans.

Chromosome-level genome assembly of the Muscovy duck provides insight into fatty liver susceptibility

The Muscovy duck (Cairina moschata) is an economically important poultry species, which is susceptible to fatty liver. Thus, the Muscovy duck may serve as an excellent candidate animal model of non-alcoholic fatty liver disease. However, the mechanisms underlying fatty liver development in this species are poorly understood. In this study, we report a chromosome-level genome assembly of the Muscovy duck, with a contig N50 of 11.8 Mb and scaffold N50 of 83.16 Mb. The susceptibility of Muscovy duck to fatty liver was mainly attributed to weak lipid catabolism capabilities (fatty acid β-oxidation and lipolysis). Furthermore, conserved noncoding elements (CNEs) showing accelerated evolution contributed to fatty liver formation by down-regulating the expression of genes involved in hepatic lipid catabolism. We propose that the susceptibility of Muscovy duck to fatty liver is an evolutionary by-product. In conclusion, this study revealed the potential mechanisms underlying the susceptibility of Muscovy duck to fatty liver.

The tissue-specificity associated region and motif of an emx2 downstream enhancer CNE2.04 in zebrafish

BACKGROUND

Expression level of EMX2 plays an important role in the development of nervous system and cancers. CNE2.04, a conserved enhancer downstream of emx2, drives fluorescent protein expression in the similar pattern of emx2.

METHODS

CNE2.04 truncated or motif-mutated transgenic reporter plasmids were constructed and injected into the zebrafish fertilized egg with Tol2 mRNA at the unicellular stage of zebrafish eggs. The green fluorescence expression patterns were observed at 24, 48, and 72 hpf, and the fluorescence rates of different tissues were counted at 48 hpf.

RESULTS

Compared to CNE2.04, CNE2.04-R400 had comparable enhancer activity, while the tissue specificity of CNE2.04-L400 was obviously changed. Motif CCCCTC mutation obviously changed the enhancer activity, while motif CCGCTC mutations also changed it.

CONCLUSION

Due to their correlation with tissue specificity, CNE2.04-R400 is associated with the tissue-specificity of CNE2.04, and motif CCCCTC plays an important role in the enhancer activity of CNE2.04.

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.

Chromosome-level genome assembly of grass carp (Ctenopharyngodon idella) provides insights into its genome evolution

Background

The grass carp has great economic value and occupies an important evolutionary position. Genomic information regarding this species could help better understand its rapid growth rate as well as its unique body plan and environmental adaptation.

Results

We assembled the chromosome-level grass carp genome using the PacBio sequencing and chromosome structure capture technique. The final genome assembly has a total length of 893.2 Mb with a contig N50 of 19.3 Mb and a scaffold N50 of 35.7 Mb. About 99.85% of the assembled contigs were anchored into 24 chromosomes. Based on the prediction, this genome contained 30,342 protein-coding genes and 43.26% repetitive sequences. Furthermore, we determined that the large genome size can be attributed to the DNA-mediated transposable elements which accounted for 58.9% of the repetitive sequences in grass carp. We identified that the grass carp has only 24 pairs of chromosomes due to the fusion of two ancestral chromosomes. Enrichment analyses of significantly expanded and positively selected genes reflected evolutionary adaptation of grass carp to the feeding habits. We also detected the loss of conserved non-coding regulatory elements associated with the development of the immune system, nervous system, and digestive system, which may be critical for grass carp herbivorous traits.

Conclusions

The high-quality reference genome reported here provides a valuable resource for the genetic improvement and molecular-guided breeding of the grass carp.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08503-x.

Perfect and imperfect views of ultraconserved sequences

Across the human genome, there are nearly 500 'ultraconserved' elements: regions of at least 200 contiguous nucleotides that are perfectly conserved in both the mouse and rat genomes. Remarkably, the majority of these sequences are non-coding, and many can function as enhancers that activate tissue-specific gene expression during embryonic development. From their first description more than 15 years ago, their extreme conservation has both fascinated and perplexed researchers in genomics and evolutionary biology. The intrigue around ultraconserved elements only grew with the observation that they are dispensable for viability. Here, we review recent progress towards understanding the general importance and the specific functions of ultraconserved sequences in mammalian development and human disease and discuss possible explanations for their extreme conservation.

The Genomes of Two Billfishes Provide Insights into the Evolution of Endothermy in Teleosts

Endothermy is a typical convergent phenomenon which has evolved independently at least eight times in vertebrates, and is of significant advantage to organisms in extending their niches. However, how vertebrates other than mammals or birds, especially teleosts, achieve endothermy has not previously been fully understood. In this study, we sequenced the genomes of two billfishes (swordfish and sailfish), members of a representative lineage of endothermic teleosts. Convergent amino acid replacements were observed in proteins related to heat production and the visual system in two endothermic teleost lineages, billfishes and tunas. The billfish-specific genetic innovations were found to be associated with heat exchange, thermoregulation, and the specialized morphology, including elongated bill, enlarged dorsal fin in sailfish and loss of the pelvic fin in swordfish.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

Chromosome-level assembly of southern catfish (silurus meridionalis) provides insights into visual adaptation to nocturnal and benthic lifestyles

The Southern catfish (Silurus meridionalis) is a nocturnal and benthic freshwater fish endemic to the Yangtze River and its tributaries. In this study, we constructed a chromosome-level draft genome of S. meridionalis using 69.7-Gb Nanopore long reads and 49.5-Gb Illumina short reads. The genome assembly was 741.2 Mb in size with a contig N50 of 13.19 Mb. An additional 116.4 Gb of Bionano and 77.4 Gb of Hi-C data were applied to assemble contigs into scaffolds and further into 29 chromosomes, resulting in a 738.9-Mb genome with a scaffold N50 of 28.04 Mb. A total of 22,965 protein-coding genes were predicted from the genome with 22,519 (98.06%) genes functionally annotated. Comparative genomic and transcriptomic analyses revealed a rod-dominated visual system which was responsible for scotopic vision. The absence of cone opsins SWS1 and SWS2 resulted in the lack of ultraviolet and blue violet sensitivity. Mutations at key amino acid sites of RH1.1, RH1.2 and RH2 resulted in spectral tuning good for dim light vision and narrow colour vision. A higher expression level of rod phototransduction genes than that of cone genes and higher rod-to-cone ratio led to higher optical sensitivity under dim light conditions. In addition, analysis of the genes involved in eye morphogenesis and development revealed the loss of some conserved noncoding elements, which might be associated with the small eyes in catfish. Together, our study provides important clues for the adaptation of the catfish visual system to the nocturnal and benthic lifestyles. The draft genome of S. meridionalis represents a valuable resource for studies of the molecular mechanisms of ecological adaptation.

Hemodynamic Forces Regulate Cardiac Regeneration-Responsive Enhancer Activity during Ventricle Regeneration

Cardiac regenerative capacity varies widely among vertebrates. Zebrafish can robustly regenerate injured hearts and are excellent models to study the mechanisms of heart regeneration. Recent studies have shown that enhancers are able to respond to injury and regulate the regeneration process. However, the mechanisms to activate these regeneration-responsive enhancers (RREs) remain poorly understood. Here, we utilized transient and transgenic analysis combined with a larval zebrafish ventricle ablation model to explore the activation and regulation of a representative RRE. lepb-linked enhancer sequence (LEN) directed enhanced green fluorescent protein (EGFP) expression in response to larval ventricle regeneration and such activation was attenuated by hemodynamic force alteration and mechanosensation pathway modulation. Further analysis revealed that Notch signaling influenced the endocardial LEN activity as well as endogenous lepb expression. Altogether, our work has established zebrafish models for rapid characterization of cardiac RREs in vivo and provides novel insights on the regulation of LEN by hemodynamic forces and other signaling pathways during heart regeneration.

Evaluating Enhancer Function and Transcription

Cell-type- and condition-specific profiles of gene expression require coordination between protein-coding gene promoters and cis-regulatory sequences called enhancers. Enhancers can stimulate gene activity at great genomic distances from their targets, raising questions about how enhancers communicate with specific gene promoters and what molecular mechanisms underlie enhancer function. Characterization of enhancer loci has identified the molecular features of active enhancers that accompany the binding of transcription factors and local opening of chromatin. These characteristics include coactivator recruitment, histone modifications, and noncoding RNA transcription. However, it remains unclear which of these features functionally contribute to enhancer activity. Here, we discuss what is known about how enhancers regulate their target genes and how enhancers and promoters communicate. Further, we describe recent data demonstrating many similarities between enhancers and the gene promoters they control, and we highlight unanswered questions in the field, such as the potential roles of transcription at enhancers.

A novel measure of non-coding genome conservation identifies genomic regulatory blocks within primates

MOTIVATION

Clusters of extremely conserved non-coding elements (CNEs) mark genomic regions devoted to cis-regulation of key developmental genes in Metazoa. We have recently shown that their span coincides with that of topologically associating domains (TADs), making them useful for estimating conserved TAD boundaries in the absence of Hi-C data. The standard approach-detecting CNEs in genome alignments and then establishing the boundaries of their clusters-requires tuning of several parameters and breaks down when comparing closely related genomes.

RESULTS

We present a novel, kurtosis-based measure of pairwise non-coding conservation that requires no pre-set thresholds for conservation level and length of CNEs. We show that it performs robustly across a large span of evolutionary distances, including across the closely related genomes of primates for which standard approaches fail. The method is straightforward to implement and enables detection and comparison of clusters of CNEs and estimation of underlying TADs across a vastly increased range of Metazoan genomes.

AVAILABILITY AND IMPLEMENTATION

The data generated for this study, and the scripts used to generate the data, can be found at https://github.com/alexander-nash/kurtosis_conservation.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Bone morphogenetic protein signaling regulates Id1-mediated neural stem cell quiescence in the adult zebrafish brain via a phylogenetically conserved enhancer module

In the telencephalon of adult zebrafish, the inhibitor of DNA binding 1 (id1) gene is expressed in radial glial cells (RGCs), behaving as neural stem cells (NSCs), during constitutive and regenerative neurogenesis. Id1 controls the balance between resting and proliferating states of RGCs by promoting quiescence. Here, we identified a phylogenetically conserved cis-regulatory module (CRM) mediating the specific expression of id1 in RGCs. Systematic deletion mapping and mutation of conserved transcription factor binding sites in stable transgenic zebrafish lines reveal that this CRM operates via conserved smad1/5 and 4 binding motifs under both homeostatic and regenerative conditions. Transcriptome analysis of injured and uninjured telencephala as well as pharmacological inhibition experiments identify a crucial role of bone morphogenetic protein (BMP) signaling for the function of the CRM. Our data highlight that BMP signals control id1 expression and thus NSC proliferation during constitutive and induced neurogenesis.

Studies of mice deleted for Sox3 and uc482: relevance to X-linked hypoparathyroidism

Hypoparathyroidism is genetically heterogeneous and characterized by low plasma calcium and parathyroid hormone (PTH) concentrations. X-linked hypoparathyroidism (XLHPT) in two American families, is associated with interstitial deletion-insertions involving deletions of chromosome Xq27.1 downstream of SOX3 and insertions of predominantly non-coding DNA from chromosome 2p25.3. These could result in loss, gain, or movement of regulatory elements, which include ultraconserved element uc482, that could alter SOX3 expression,. To investigate this, we analysed SOX3 expression in EBV-transformed lymphoblastoid cells from 3 affected males, 3 unaffected males, and 4 carrier females from one XLHPT family. SOX3 expression was similar in all individuals, indicating that the spatiotemporal effect of the interstitial deletion-insertion on SOX3 expression postulated to occur in developing parathyroids did not manifest in lymphoblastoids. Expression of SNTG2, which is duplicated and inserted into the X chromosome, and ATP11C, which is moved telomerically, were also similarly expressed in all individuals. Investigation of male hemizygous (Sox3-/Y and uc482-/Y) and female heterozygous (Sox3+/- and uc482+/-) knock-out mice, together with wild-type littermates (male Sox3+/Y and uc482+/Y, and female Sox3+/+ and uc482+/+), revealed Sox3-/Y, Sox3+/-, uc482-/Y, and uc482+/- mice to have normal plasma biochemistry, compared to their respective wild-type littermates. When challenged with a low calcium diet, all mice had hypocalcaemia, and elevated plasma PTH concentrations and alkaline phosphatase activities, and Sox3-/Y, Sox3+/-, uc482-/Y, and uc482+/- mice had similar plasma biochemistry, compared to wild-type littermates. Thus, these results indicate that absence of Sox3 or uc482 does not cause hypoparathyroidism, and that XLHPT likely reflects a more complex mechanism.

Understanding the genetics of neuropsychiatric disorders: the potential role of genomic regulatory blocks

Recent genome-wide association studies have identified numerous loci associated with neuropsychiatric disorders. The majority of these are in non-coding regions, and are commonly assigned to the nearest gene along the genome. However, this approach neglects the three-dimensional organisation of the genome, and the fact that the genome contains arrays of extremely conserved non-coding elements termed genomic regulatory blocks (GRBs), which can be utilized to detect genes under long-range developmental regulation. Here we review a GRB-based approach to assign loci in non-coding regions to potential target genes, and apply it to reanalyse the results of one of the largest schizophrenia GWAS (SWG PGC, 2014). We further apply this approach to GWAS data from two related neuropsychiatric disorders-autism spectrum disorder and bipolar disorder-to show that it is applicable to developmental disorders in general. We find that disease-associated SNPs are overrepresented in GRBs and that the GRB model is a powerful tool for linking these SNPs to their correct target genes under long-range regulation. Our analysis identifies novel genes not previously implicated in schizophrenia and corroborates a number of predicted targets from the original study. The results are available as an online resource in which the genomic context and the strength of enhancer-promoter associations can be browsed for each schizophrenia-associated SNP.

A screen for deeply conserved non-coding GWAS SNPs uncovers a MIR-9-2 functional mutation associated to retinal vasculature defects in human

Thousands of human disease-associated single nucleotide polymorphisms (SNPs) lie in the non-coding genome, but only a handful have been demonstrated to affect gene expression and human biology. We computationally identified risk-associated SNPs in deeply conserved non-exonic elements (CNEs) potentially contributing to 45 human diseases. We further demonstrated that human CNE1/rs17421627 associated with retinal vasculature defects showed transcriptional activity in the zebrafish retina, while introducing the risk-associated allele completely abolished CNE1 enhancer activity. Furthermore, deletion of CNE1 led to retinal vasculature defects and to a specific downregulation of microRNA-9, rather than MEF2C as predicted by the original genome-wide association studies. Consistent with these results, miR-9 depletion affects retinal vasculature formation, demonstrating MIR-9-2 as a critical gene underpinning the associated trait. Importantly, we validated that other CNEs act as transcriptional enhancers that can be disrupted by conserved non-coding SNPs. This study uncovers disease-associated non-coding mutations that are deeply conserved, providing a path for in vivo testing to reveal their cis-regulated genes and biological roles.

The Divergent Genomes of Teleosts

Boasting nearly 30,000 species, teleosts account for half of all extant vertebrates and approximately 98% of all ray-finned fish species (Actinopterygii). Teleosts are also the largest and most diverse group of vertebrates, exhibiting an astonishing level of morphological, physiological, and behavioral diversity. Previous studies had indicated that the teleost lineage has experienced an additional whole-genome duplication event. Recent comparative genomic analyses of teleosts and other bony vertebrates using spotted gar (a nonteleost ray-finned fish) and elephant shark (a cartilaginous fish) as outgroups have revealed several divergent features of teleost genomes. These include an accelerated evolutionary rate of protein-coding and nucleotide sequences, a higher rate of intron turnover, loss of many potential cis-regulatory elements and shorter conserved syntenic blocks. A combination of these divergent genomic features might have contributed to the evolution of the amazing phenotypic diversity and morphological innovations of teleosts.

The HMG box transcription factors Sox1a and b specify a new class of glycinergic interneurons in the spinal cord of zebrafish embryos

Specification of neurons in the spinal cord relies on extrinsic and intrinsic signals, which in turn are interpreted by expression of transcription factors. V2 interneurons develop from the ventral aspects of the spinal cord. We report here a novel neuronal V2 subtype, named V2s, in zebrafish embryos. Formation of these neurons depends on the transcription factors sox1a and sox1b. They develop from common gata2a/gata3 dependent precursors co-expressing markers of V2b and V2s interneurons. Chemical blockage of Notch signaling causes a decrease of V2s and an increase of V2b cells. Our results are consistent with the existence of at least two types of precursors arranged in a hierarchical manner in the V2 domain. V2s neurons grow long ipsilateral descending axonal projections with a short branch at the ventral midline. They acquire a glycinergic neurotransmitter type during the second day of development. Unilateral ablation of V2s interneurons causes a delay in touch-provoked escape behavior suggesting that V2s interneurons are involved in fast motor responses.

Identification of novel cis-regulatory elements of Eya1 in Xenopus laevis using BAC recombineering

The multifunctional Eya1 protein plays important roles during the development of cranial sensory organs and ganglia, kidneys, hypaxial muscles and several other organs in vertebrates. Eya1 is encoded by a complex locus with candidate cis-regulatory elements distributed over a 329 kbp wide genomic region in Xenopus. Consequently, very little is currently known about how expression of Eya1 is controlled by upstream regulators. Here we use a library of Xenopus tropicalis genomic sequences in bacterial artificial chromosomes (BAC) to analyze the genomic region surrounding the Eya1 locus for enhancer activity. We used BAC recombineering to first create GFP reporter constructs, which were analysed for enhancer activity by injection into Xenopus laevis embryos. We then used a second round of BAC recombineering to create deletion constructs of these BAC reporters to localize enhancer activity more precisely. This double recombineering approach allowed us to probe a large genomic region for enhancer activity without assumptions on sequence conservation. Using this approach we were able to identify two novel cis-regulatory regions, which direct Eya1 expression to the somites, pharyngeal pouches, the preplacodal ectoderm (the common precursor region of many cranial sensory organs and ganglia), and other ectodermal domains.

Topologically associating domains are ancient features that coincide with Metazoan clusters of extreme noncoding conservation

Developmental genes in metazoan genomes are surrounded by dense clusters of conserved noncoding elements (CNEs). CNEs exhibit unexplained extreme levels of sequence conservation, with many acting as developmental long-range enhancers. Clusters of CNEs define the span of regulatory inputs for many important developmental regulators and have been described previously as genomic regulatory blocks (GRBs). Their function and distribution around important regulatory genes raises the question of how they relate to 3D conformation of these loci. Here, we show that clusters of CNEs strongly coincide with topological organisation, predicting the boundaries of hundreds of topologically associating domains (TADs) in human and Drosophila. The set of TADs that are associated with high levels of noncoding conservation exhibit distinct properties compared to TADs devoid of extreme noncoding conservation. The close correspondence between extreme noncoding conservation and TADs suggests that these TADs are ancient, revealing a regulatory architecture conserved over hundreds of millions of years.

Metazoan genomes contain many clusters of conserved noncoding elements. Here, the authors provide evidence that these clusters coincide with distinct topologically associating domains in humans and Drosophila, revealing a conserved regulatory genomic architecture.

The seahorse genome and the evolution of its specialized morphology

Seahorses have a specialized morphology that includes a toothless tubular mouth, a body covered with bony plates, a male brood pouch, and the absence of caudal and pelvic fins. Here we report the sequencing and de novo assembly of the genome of the tiger tail seahorse, Hippocampus comes. Comparative genomic analysis identifies higher protein and nucleotide evolutionary rates in H. comes compared with other teleost fish genomes. We identified an astacin metalloprotease gene family that has undergone expansion and is highly expressed in the male brood pouch. We also find that the H. comes genome lacks enamel matrix protein-coding proline/glutamine-rich secretory calcium-binding phosphoprotein genes, which might have led to the loss of mineralized teeth. tbx4, a regulator of hindlimb development, is also not found in H. comes genome. Knockout of tbx4 in zebrafish showed a ‘pelvic fin-loss’ phenotype similar to that of seahorses.

Here, the genome sequence of the tiger tail seahorse is reported and comparative genomic analyses with other ray-finned fishes are used to explore the genetic basis of the unique morphology and reproductive system of the seahorse.

Seahorses are prime examples of the exuberance of evolution and are unique among bony fish on several counts, including their equine body shape and male brood pouch. An international collaboration reporting in this issue of Nature has determined the genome sequence of a seahorse (Hippocampus comes, the tiger tail seahorse). They find it to be the most rapidly evolving fish genome studied so far. H. comes is among the most commonly traded seahorse species—dried for traditional medicines and live for the aquarium trade—and is on the IUCN Red List as a 'vulnerable' species. Analysis of the genomic sequence provides insights into the evolution of its unique morphology. Of note is the absence of a master control gene, tbx4, which functions in the development of hindlimbs and pelvic fins. Pelvic fins are missing in seahorses, and tbx4-knockout mutant zebrafish also lack pelvic fins.

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.

Functions of noncoding sequences in mammalian genomes

Most of the mammalian genome consists of nucleotide sequences not coding for proteins. Exons of genes make up only 3% of the human genome, while the significance of most other sequences remains unknown. Recent genome studies with high-throughput methods demonstrate that the so-called noncoding part of the genome may perform important functions. This hypothesis is supported by three groups of experimental data: 1) approximately 10% of the sequences, most of which are located in noncoding parts of the genome, is evolutionarily conserved and thus can be of functional importance; 2) up to 99% of the mammalian genome is being transcribed forming short and long noncoding RNAs in addition to common mRNA; and 3) mutations in noncoding parts of the genome can be accompanied by progression of pathological states of the organism. In the light of these data, in the review we consider the functional role of numerous known sequences of noncoding parts of the genome including introns, DNA methylation regions, enhancers and locus control regions, insulators, S/MAR sequences, pseudogenes, and genes of noncoding RNAs, as well as transposons and simple repeats of centromeric and telomeric regions of chromosomes. The assumption is made that the intergenic noncoding sequences without definite/clear functions can be involved in spatial organization of genetic loci in interphase nuclei.

DENdb: database of integrated human enhancers

Enhancers are cis-acting DNA regulatory regions that play a key role in distal control of transcriptional activities. Identification of enhancers, coupled with a comprehensive functional analysis of their properties, could improve our understanding of complex gene transcription mechanisms and gene regulation processes in general. We developed DENdb, a centralized on-line repository of predicted enhancers derived from multiple human cell-lines. DENdb integrates enhancers predicted by five different methods generating an enriched catalogue of putative enhancers for each of the analysed cell-lines. DENdb provides information about the overlap of enhancers with DNase I hypersensitive regions, ChIP-seq regions of a number of transcription factors and transcription factor binding motifs, means to explore enhancer interactions with DNA using several chromatin interaction assays and enhancer neighbouring genes. DENdb is designed as a relational database that facilitates fast and efficient searching, browsing and visualization of information. Database URL: http://www.cbrc.kaust.edu.sa/dendb/.

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.

Phylum-Level Conservation of Regulatory Information in Nematodes despite Extensive Non-coding Sequence Divergence

Gene regulatory information guides development and shapes the course of evolution. To test conservation of gene regulation within the phylum Nematoda, we compared the functions of putative cis-regulatory sequences of four sets of orthologs (unc-47, unc-25, mec-3 and elt-2) from distantly-related nematode species. These species, Caenorhabditis elegans, its congeneric C. briggsae, and three parasitic species Meloidogyne hapla, Brugia malayi, and Trichinella spiralis, represent four of the five major clades in the phylum Nematoda. Despite the great phylogenetic distances sampled and the extensive sequence divergence of nematode genomes, all but one of the regulatory elements we tested are able to drive at least a subset of the expected gene expression patterns. We show that functionally conserved cis-regulatory elements have no more extended sequence similarity to their C. elegans orthologs than would be expected by chance, but they do harbor motifs that are important for proper expression of the C. elegans genes. These motifs are too short to be distinguished from the background level of sequence similarity, and while identical in sequence they are not conserved in orientation or position. Functional tests reveal that some of these motifs contribute to proper expression. Our results suggest that conserved regulatory circuitry can persist despite considerable turnover within cis elements.

To explore the phylogenetic limits of conservation of cis-regulatory elements, we used transgenesis to test the functions of enhancers of four genes from several species spanning the phylum Nematoda. While we found a striking degree of functional conservation among the examined cis elements, their DNA sequences lacked apparent conservation with the C. elegans orthologs. In fact, sequence similarity between C. elegans and the distantly related nematodes was no greater than would be expected by chance. Short motifs, similar to known regulatory sequences in C. elegans, can be detected in most of the cis elements. When tested, some of these sites appear to mediate regulatory function. However, they seem to have originated through motif turnover, rather than to have been preserved from a common ancestor. Our results suggest that gene regulatory networks are broadly conserved in the phylum Nematoda, but this conservation persists despite substantial reorganization of regulatory elements and could not be detected using naïve comparisons of sequence similarity.

A PAX1 enhancer locus is associated with susceptibility to idiopathic scoliosis in females

Idiopathic scoliosis (IS) is a common paediatric musculoskeletal disease that displays a strong female bias. By performing a genome-wide association study (GWAS) of 3,102 individuals, we identify significant associations with 20p11.22 SNPs for females (P=6.89 × 10⁻⁹) but not males (P=0.71). This association with IS is also found in independent female cohorts from the United States of America and Japan (overall P=2.15 × 10⁻¹⁰, OR=1.30 (rs6137473)). Unexpectedly, the 20p11.22 IS risk alleles were previously associated with protection from early-onset alopecia, another sexually dimorphic condition. The 174-kb associated locus is distal to PAX1, which encodes paired box 1, a transcription factor involved in spine development. We identify a sequence in the associated locus with enhancer activity in zebrafish somitic muscle and spinal cord, an activity that is abolished by IS-associated SNPs. We thus identify a sexually dimorphic IS susceptibility locus, and propose the first functionally defined candidate mutations in an enhancer that may regulate expression in specific spinal cells.

Girls are tenfold more likely than boys to require surgical treatment for idiopathic scoliosis, a common paediatric skeletal disorder. Here, Sharma et al. identify the first sexually dimorphic idiopathic scoliosis risk locus, and demonstrate that it may play a role in the regulation of spinal cells.

Identification and computational analysis of gene regulatory elements

Over the last two decades, advances in experimental and computational technologies have greatly facilitated genomic research. Next-generation sequencing technologies have made de novo sequencing of large genomes affordable, and powerful computational approaches have enabled accurate annotations of genomic DNA sequences. Charting functional regions in genomes must account for not only the coding sequences, but also noncoding RNAs, repetitive elements, chromatin states, epigenetic modifications, and gene regulatory elements. A mix of comparative genomics, high-throughput biological experiments, and machine learning approaches has played a major role in this truly global effort. Here we describe some of these approaches and provide an account of our current understanding of the complex landscape of the human genome. We also present overviews of different publicly available, large-scale experimental data sets and computational tools, which we hope will prove beneficial for researchers working with large and complex genomes.

Immediate and long-term consequences of vascular toxicity during zebrafish development

Proper formation of the vascular system is necessary for embryogenesis, and chemical disruption of vascular development may be a key event driving developmental toxicity. In order to test the effect of environmental chemicals on this critical process, we evaluated a quantitative assay in transgenic zebrafish using angiogenesis inhibitors that target VEGFR2 (PTK787) or EGFR (AG1478). Both PTK787 and AG1478 exposure impaired intersegmental vessel (ISV) sprouting, while AG1478 also produced caudal and pectoral fin defects at concentrations below those necessary to blunt ISV morphogenesis. The functional consequences of vessel toxicity during early development included decreased body length and survival in juvenile cohorts developmentally exposed to inhibitor concentrations sufficient to completely block ISV sprouting angiogenesis. These data show that concentration-dependent disruption of the presumed targets for these inhibitors produce adverse outcomes at advanced life stages.

Conserved noncoding elements follow power-law-like distributions in several genomes as a result of genome dynamics

Conserved, ultraconserved and other classes of constrained elements (collectively referred as CNEs here), identified by comparative genomics in a wide variety of genomes, are non-randomly distributed across chromosomes. These elements are defined using various degrees of conservation between organisms and several thresholds of minimal length. We here investigate the chromosomal distribution of CNEs by studying the statistical properties of distances between consecutive CNEs. We find widespread power-law-like distributions, i.e. linearity in double logarithmic scale, in the inter-CNE distances, a feature which is connected with fractality and self-similarity. Given that CNEs are often found to be spatially associated with genes, especially with those that regulate developmental processes, we verify by appropriate gene masking that a power-law-like pattern emerges irrespectively of whether elements found close or inside genes are excluded or not. An evolutionary model is put forward for the understanding of these findings that includes segmental or whole genome duplication events and eliminations (loss) of most of the duplicated CNEs. Simulations reproduce the main features of the observed size distributions. Power-law-like patterns in the genomic distributions of CNEs are in accordance with current knowledge about their evolutionary history in several genomes.

Elucidation of distinctive genomic DNA structures in patients with 46,XX testicular disorders of sex development using genome wide analyses

PURPOSE

Although several genes, including the SRY gene, are involved in testicular differentiation, the entire mechanism of this differentiation remains unclear. We performed genome wide analysis in patients with 46,XX testicular disorders of sex development to comprehensively elucidate the mechanisms of testicular differentiation.

MATERIALS AND METHODS

Whole genomic DNA was extracted from the peripheral blood of 4 patients with 46,XX testicular disorders of sex development who were SRY negative. Genomic DNA was hybridized to a GeneChip® human mapping 250K array set. Compared to normal female data, we detected common loss of heterozygosity and copy number variation regions in 4 patients using Genotyping Console™ software.

RESULTS

Loss of heterozygosity was detected in 19 regions of 11 chromosomes. A total of 27 genes or nearby genomic areas were included in the applicable regions. Copy number loss was recognized in 13 regions of 10 chromosomes, and these regions included 55 genes. Copy number gain was detected in 6 regions of 4 chromosomes, which included the upstream region of the SOX3 gene.

CONCLUSIONS

The regions with loss of heterozygosity did not contain genes associated with testicular differentiation. However, the upstream area of the SOX3 gene, which is located in Xq27.1, was included in the region of copy number gain. These results suggest that high expression of the SOX3 gene led to testicular differentiation despite SRY gene loss. As this applicable area is not within a coding region, genome wide analyses were valuable for detecting the novel regions associated with testicular differentiation.

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.

A reporter assay in lamprey embryos reveals both functional conservation and elaboration of vertebrate enhancers

The sea lamprey is an important model organism for investigating the evolutionary origins of vertebrates. As more vertebrate genome sequences are obtained, evolutionary developmental biologists are becoming increasingly able to identify putative gene regulatory elements across the breadth of the vertebrate taxa. The identification of these regions makes it possible to address how changes at the genomic level have led to changes in developmental gene regulatory networks and ultimately to the evolution of morphological diversity. Comparative genomics approaches using sea lamprey have already predicted a number of such regulatory elements in the lamprey genome. Functional characterisation of these sequences and other similar elements requires efficient reporter assays in lamprey. In this report, we describe the development of a transient transgenesis method for lamprey embryos. Focusing on conserved non-coding elements (CNEs), we use this method to investigate their functional conservation across the vertebrate subphylum. We find instances of both functional conservation and lineage-specific functional evolution of CNEs across vertebrates, emphasising the utility of functionally testing homologous CNEs in their host species.

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.

Computational methods to detect conserved non-genic elements in phylogenetically isolated genomes: application to zebrafish

Many important model organisms for biomedical and evolutionary research have sequenced genomes, but occupy a phylogenetically isolated position, evolutionarily distant from other sequenced genomes. This phylogenetic isolation is exemplified for zebrafish, a vertebrate model for cis-regulation, development and human disease, whose evolutionary distance to all other currently sequenced fish exceeds the distance between human and chicken. Such large distances make it difficult to align genomes and use them for comparative analysis beyond gene-focused questions. In particular, detecting conserved non-genic elements (CNEs) as promising cis-regulatory elements with biological importance is challenging. Here, we develop a general comparative genomics framework to align isolated genomes and to comprehensively detect CNEs. Our approach integrates highly sensitive and quality-controlled local alignments and uses alignment transitivity and ancestral reconstruction to bridge large evolutionary distances. We apply our framework to zebrafish and demonstrate substantially improved CNE detection and quality compared with previous sets. Our zebrafish CNE set comprises 54 533 CNEs, of which 11 792 (22%) are conserved to human or mouse. Our zebrafish CNEs (http://zebrafish.stanford.edu) are highly enriched in known enhancers and extend existing experimental (ChIP-Seq) sets. The same framework can now be applied to the isolated genomes of frog, amphioxus, Caenorhabditis elegans and many others.

Conserved non-coding elements and cis regulation: actions speak louder than words

It is a truth (almost) universally acknowledged that conserved non-coding genomic sequences function in the cis regulation of neighbouring genes. But is this a misconception? The literature is strewn with examples of conserved non-coding sequences being able to drive reporter expression, but the extent to which such sequences are actually used endogenously in vivo is only now being rigorously explored using unbiased genome-scale approaches. Here, we review the emerging picture, examining the extent to which conserved non-coding sequences equivalently regulate gene expression in different species, or at different developmental stages, and how genomics approaches are revealing the relationship between sequence conservation and functional use of cis-regulatory elements.

The African coelacanth genome provides insights into tetrapod evolution

It was a zoological sensation when a living specimen of the coelacanth was first discovered in 1938, as this lineage of lobe-finned fish was thought to have gone extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features . Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain, and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues demonstrate the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

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.

Co-orthology of Pax4 and Pax6 to the fly eyeless gene: molecular phylogenetic, comparative genomic, and embryological analyses

The functional equivalence of Pax6/eyeless genes across distantly related animal phyla has been one of central findings on which evo-devo studies is based. In this study, we show that Pax4, in addition to Pax6, is a vertebrate ortholog of the fly eyeless gene (and its duplicate, twin of eyeless [toy] gene, unique to Insecta). Molecular phylogenetic trees published to date placed the Pax4 gene outside the Pax6/eyeless subgroup as if the Pax4 gene originated from a gene duplication before the origin of bilaterians. However, Pax4 genes had only been reported for mammals. Our molecular phylogenetic analysis, including previously unidentified teleost fish pax4 genes, equally supported two scenarios: one with the Pax4-Pax6 duplication early in vertebrate evolution and the other with this duplication before the bilaterian radiation. We then investigated gene compositions in the genomic regions containing Pax4 and Pax6, and identified (1) conserved synteny between these two regions, suggesting that the Pax4-Pax6 split was caused by a large-scale duplication and (2) its timing within early vertebrate evolution based on the duplication timing of the members of neighboring gene families. Our results are consistent with the so-called two-round genome duplications in early vertebrates. Overall, the Pax6/eyeless ortholog is merely part of a 2:2 orthology relationship between vertebrates (with Pax4 and Pax6) and the fly (with eyeless and toy). In this context, evolution of transcriptional regulation associated with the Pax4-Pax6 split is also discussed in light of the zebrafish pax4 expression pattern that is analyzed here for the first time.

Sequencing of Pax6 loci from the elephant shark reveals a family of Pax6 genes in vertebrate genomes, forged by ancient duplications and divergences

Pax6 is a developmental control gene essential for eye development throughout the animal kingdom. In addition, Pax6 plays key roles in other parts of the CNS, olfactory system, and pancreas. In mammals a single Pax6 gene encoding multiple isoforms delivers these pleiotropic functions. Here we provide evidence that the genomes of many other vertebrate species contain multiple Pax6 loci. We sequenced Pax6-containing BACs from the cartilaginous elephant shark (Callorhinchus milii) and found two distinct Pax6 loci. Pax6.1 is highly similar to mammalian Pax6, while Pax6.2 encodes a paired-less Pax6. Using synteny relationships, we identify homologs of this novel paired-less Pax6.2 gene in lizard and in frog, as well as in zebrafish and in other teleosts. In zebrafish two full-length Pax6 duplicates were known previously, originating from the fish-specific genome duplication (FSGD) and expressed in divergent patterns due to paralog-specific loss of cis-elements. We show that teleosts other than zebrafish also maintain duplicate full-length Pax6 loci, but differences in gene and regulatory domain structure suggest that these Pax6 paralogs originate from a more ancient duplication event and are hence renamed as Pax6.3. Sequence comparisons between mammalian and elephant shark Pax6.1 loci highlight the presence of short- and long-range conserved noncoding elements (CNEs). Functional analysis demonstrates the ancient role of long-range enhancers for Pax6 transcription. We show that the paired-less Pax6.2 ortholog in zebrafish is expressed specifically in the developing retina. Transgenic analysis of elephant shark and zebrafish Pax6.2 CNEs with homology to the mouse NRE/Pα internal promoter revealed highly specific retinal expression. Finally, morpholino depletion of zebrafish Pax6.2 resulted in a "small eye" phenotype, supporting a role in retinal development. In summary, our study reveals that the pleiotropic functions of Pax6 in vertebrates are served by a divergent family of Pax6 genes, forged by ancient duplication events and by independent, lineage-specific gene losses.

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.

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.

Pax4 is not essential for beta-cell differentiation in zebrafish embryos but modulates alpha-cell generation by repressing arx gene expression

BACKGROUND

Genetic studies in mouse have demonstrated the crucial function of PAX4 in pancreatic cell differentiation. This transcription factor specifies β- and δ-cell fate at the expense of α-cell identity by repressing Arx gene expression and ectopic expression of PAX4 in α-cells is sufficient to convert them into β-cells. Surprisingly, no Pax4 orthologous gene can be found in chicken and Xenopus tropicalis raising the question of the function of pax4 gene in lower vertebrates such as in fish. In the present study, we have analyzed the expression and the function of the orthologous pax4 gene in zebrafish.

RESULTS

pax4 gene is transiently expressed in the pancreas of zebrafish embryos and is mostly restricted to endocrine precursors as well as to some differentiating δ- and ε-cells but was not detected in differentiating β-cells. pax4 knock-down in zebrafish embryos caused a significant increase in α-cells number while having no apparent effect on β- and δ-cell differentiation. This rise of α-cells is due to an up-regulation of the Arx transcription factor. Conversely, knock-down of arx caused to a complete loss of α-cells and a concomitant increase of pax4 expression but had no effect on the number of β- and δ-cells. In addition to the mutual repression between Arx and Pax4, these two transcription factors negatively regulate the transcription of their own gene. Interestingly, disruption of pax4 RNA splicing or of arx RNA splicing by morpholinos targeting exon-intron junction sites caused a blockage of the altered transcripts in cell nuclei allowing an easy characterization of the arx- and pax4-deficient cells. Such analyses demonstrated that arx knock-down in zebrafish does not lead to a switch of cell fate, as reported in mouse, but rather blocks the cells in their differentiation process towards α-cells.

CONCLUSIONS

In zebrafish, pax4 is not required for the generation of the first β- and δ-cells deriving from the dorsal pancreatic bud, unlike its crucial role in the differentiation of these cell types in mouse. On the other hand, the mutual repression between Arx and Pax4 is observed in both mouse and zebrafish. These data suggests that the main original function of Pax4 during vertebrate evolution was to modulate the number of pancreatic α-cells and its role in β-cells differentiation appeared later in vertebrate evolution.