Exonic remnants of whole-genome duplication reveal cis-regulatory function of coding exons

Cis-regulatory landscapes in the evolution and development of the mammalian skull

Extensive morphological variation found in mammals reflects the wide spectrum of their ecological adaptations. The highest morphological diversity is present in the craniofacial region, where geometry is mainly dictated by the bony skull. Mammalian craniofacial development represents complex multistep processes governed by numerous conserved genes that require precise spatio-temporal control. A central question in contemporary evolutionary biology is how a defined set of conserved genes can orchestrate formation of fundamentally different structures, and therefore how morphological variability arises. In principle, differential gene expression patterns during development are the source of morphological variation. With the emergence of multicellular organisms, precise regulation of gene expression in time and space is attributed to cis-regulatory elements. These elements contribute to higher-order chromatin structure and together with trans-acting factors control transcriptional landscapes that underlie intricate morphogenetic processes. Consequently, divergence in cis-regulation is believed to rewire existing gene regulatory networks and form the core of morphological evolution. This review outlines the fundamental principles of the genetic code and genomic regulation interplay during development. Recent work that deepened our comprehension of cis-regulatory element origin, divergence and function is presented here to illustrate the state-of-the-art research that uncovered the principles of morphological novelty. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Characterization of Two Transposable Elements and an Ultra-Conserved Element Isolated in the Genome of Zootoca vivipara (Squamata, Lacertidae)

Transposable elements (TEs) constitute a considerable fraction of eukaryote genomes representing a major source of genetic variability. We describe two DNA sequences isolated in the lizard Zootoca vivipara, here named Zv516 and Zv817. Both sequences are single-copy nuclear sequences, including a truncation of two transposable elements (TEs), SINE Squam1 in Zv516 and a Tc1/Mariner-like DNA transposon in Zv817. FISH analyses with Zv516 showed the occurrence of interspersed signals of the SINE Squam1 sequence on all chromosomes of Z. vivipara and quantitative dot blot indicated that this TE is present with about 4700 copies in the Z. vivipara genome. FISH and dot blot with Zv817 did not produce clear hybridization signals. Bioinformatic analysis showed the presence of active SINE Squam 1 copies in the genome of different lacertids, in different mRNAs, and intronic and coding regions of various genes. The Tc1/Mariner-like DNA transposon occurs in all reptiles, excluding Sphenodon and Archosauria. Zv817 includes a trait of 284 bp, representing an amniote ultra-conserved element (UCE). Using amniote UCE homologous sequences from available whole genome sequences of major amniote taxonomic groups, we performed a phylogenetic analysis which retrieved Prototheria as the sister group of Metatheria and Eutheria. Within diapsids, Testudines are the sister group to Aves + Crocodylia (Archosauria), and Sphenodon is the sister group to Squamata. Furthermore, large trait regions flanking the UCE are conserved at family level.

A map of cis-regulatory modules and constituent transcription factor binding sites in 80% of the mouse genome

BACKGROUND

Mouse is probably the most important model organism to study mammal biology and human diseases. A better understanding of the mouse genome will help understand the human genome, biology and diseases. However, despite the recent progress, the characterization of the regulatory sequences in the mouse genome is still far from complete, limiting its use to understand the regulatory sequences in the human genome.

RESULTS

Here, by integrating binding peaks in ~ 9,000 transcription factor (TF) ChIP-seq datasets that cover 79.9% of the mouse mappable genome using an efficient pipeline, we were able to partition these binding peak-covered genome regions into a cis-regulatory module (CRM) candidate (CRMC) set and a non-CRMC set. The CRMCs contain 912,197 putative CRMs and 38,554,729 TF binding sites (TFBSs) islands, covering 55.5% and 24.4% of the mappable genome, respectively. The CRMCs tend to be under strong evolutionary constraints, indicating that they are likely cis-regulatory; while the non-CRMCs are largely selectively neutral, indicating that they are unlikely cis-regulatory. Based on evolutionary profiles of the genome positions, we further estimated that 63.8% and 27.4% of the mouse genome might code for CRMs and TFBSs, respectively.

CONCLUSIONS

Validation using experimental data suggests that at least most of the CRMCs are authentic. Thus, this unprecedentedly comprehensive map of CRMs and TFBSs can be a good resource to guide experimental studies of regulatory genomes in mice and humans.

Towards unveiling the nature of short SERPINA1 transcripts: Avoiding the main ORF control to translate alpha1-antitrypsin C-terminal peptides

Alternative ORFs in-frame with the known genes are challenging to reveal. Yet they may contribute significantly to proteome diversity. Here we focused on the individual expression of the SERPINA1 gene exon 5 leading to direct translation of alpha1-antitrypsin (AAT) C-terminal peptides. The discovery of alternative ways for their production may expand the current understanding of the serpin gene's functioning. We detected short transcripts expressed primarily in hepatocytes. We identified four variants of hepatocyte-specific SERPINA1 short transcripts and individually probed their potential to be translated in living cells. The long mRNA gave the full-length AAT-eGFP fusion, while in case of short transcripts we deduced four active SERPINA1 in-frame alternative ORFs encoding 10, 21, 153 and 169 amino acids AAT C-terminal oligo- and polypeptides. Unlike secretory AAT-eGFP fusion exhibiting classical AAT behavior, truncated AAT-fusions differ by intracellular retention and nuclear enrichment. Immunofluorescence on the endogenous AAT C-terminal epitope showed its accumulation in the cell nucleoli, indicating that short transcripts may be translated in vivo. FANTOM5 CAGE data on SERPINA1 suggest that short transcripts originate from the post-transcriptional cleavage of the spliced mRNA, initiated mainly from the hepatocyte-specific promoter. CONCLUSION: Short SERPINA1 transcripts may represent a source for the direct synthesis of AAT C-terminal peptides with properties uncommon to AAT.

Using a Machine Learning Approach to Identify Low-Frequency and Rare FLG Alleles Associated with Remission of Atopic Dermatitis

Atopic dermatitis (AD) is a common relapsing inflammatory skin disease. FLG is the gene most consistently associated with AD. Loss-of-function variants in FLG have been previously associated with AD. Low-frequency and rare alleles (minor allele frequency < 5%) in this gene have been given less attention than loss-of-function variants. We fine sequenced the FLG gene in a cohort of individuals with AD. We developed a machine learning‒based algorithm to associate low-frequency and rare alleles with the disease. We then applied this algorithm to the FLG data, searching for associations between groups of low-frequency and rare FLG alleles and AD remission. A group of 46 rare and low-frequency FLG alleles was associated with increased AD remission (P = 2.76e-11). Overall, 16 of these 46 FLG variants were identified in an independent cohort and were associated with decreased AD incidence (P = 0.0007). This study presents an application of statistical methods in AD genetics and suggests that low-frequency and rare alleles may play a larger role in AD pathogenesis than previously appreciated.

Structural and functional analysis of somatic coding and UTR indels in breast and lung cancer genomes

Insertions and deletions (Indels) represent one of the major variation types in the human genome and have been implicated in diseases including cancer. To study the features of somatic indels in different cancer genomes, we investigated the indels from two large samples of cancer types: invasive breast carcinoma (BRCA) and lung adenocarcinoma (LUAD). Besides mapping somatic indels in both coding and untranslated regions (UTRs) from the cancer whole exome sequences, we investigated the overlap between these indels and transcription factor binding sites (TFBSs), the key elements for regulation of gene expression that have been found in both coding and non-coding sequences. Compared to the germline indels in healthy genomes, somatic indels contain more coding indels with higher than expected frame-shift (FS) indels in cancer genomes. LUAD has a higher ratio of deletions and higher coding and FS indel rates than BRCA. More importantly, these somatic indels in cancer genomes tend to locate in sequences with important functions, which can affect the core secondary structures of proteins and have a bigger overlap with predicted TFBSs in coding regions than the germline indels. The somatic CDS indels are also enriched in highly conserved nucleotides when compared with germline CDS indels.

Mechanisms of enhancer action: the known and the unknown

Differential gene expression mechanisms ensure cellular differentiation and plasticity to shape ontogenetic and phylogenetic diversity of cell types. A key regulator of differential gene expression programs are the enhancers, the gene-distal cis-regulatory sequences that govern spatiotemporal and quantitative expression dynamics of target genes. Enhancers are widely believed to physically contact the target promoters to effect transcriptional activation. However, our understanding of the full complement of regulatory proteins and the definitive mechanics of enhancer action is incomplete. Here, we review recent findings to present some emerging concepts on enhancer action and also outline a set of outstanding questions.

Fish-Ing for Enhancers in the Heart

Precise control of gene expression is crucial to ensure proper development and biological functioning of an organism. Enhancers are non-coding DNA elements which play an essential role in regulating gene expression. They contain specific sequence motifs serving as binding sites for transcription factors which interact with the basal transcription machinery at their target genes. Heart development is regulated by intricate gene regulatory network ensuring precise spatiotemporal gene expression program. Mutations affecting enhancers have been shown to result in devastating forms of congenital heart defect. Therefore, identifying enhancers implicated in heart biology and understanding their mechanism is key to improve diagnosis and therapeutic options. Despite their crucial role, enhancers are poorly studied, mainly due to a lack of reliable way to identify them and determine their function. Nevertheless, recent technological advances have allowed rapid progress in enhancer discovery. Model organisms such as the zebrafish have contributed significant insights into the genetics of heart development through enabling functional analyses of genes and their regulatory elements in vivo. Here, we summarize the current state of knowledge on heart enhancers gained through studies in model organisms, discuss various approaches to discover and study their function, and finally suggest methods that could further advance research in this field.

Transcriptional control by enhancers and enhancer RNAs

The regulation of gene expression is a fundamental cellular process and its misregulation is a key component of disease. Enhancers are one of the most salient regulatory elements in the genome and help orchestrate proper spatiotemporal gene expression during development, in homeostasis, and in response to signaling. Notably, molecular aberrations at enhancers, such as translocations and single nucleotide polymorphisms, are emerging as an important source of human variation and susceptibility to disease. Herein we discuss emerging paradigms addressing how genes are regulated by enhancers, common features of active enhancers, and how non-coding enhancer RNAs (eRNAs) can direct gene expression programs that underlie cellular phenotypes. We survey the current evidence, which suggests that eRNAs can bind to transcription factors, mediate enhancer-promoter interactions, influence RNA Pol II elongation, and act as decoys for repressive cofactors. Furthermore, we discuss current methodologies for the identification of eRNAs and novel approaches to elucidate their functions.

Identification of a 42-bp heart-specific enhancer of the notch1b gene in zebrafish embryos

BACKGROUND

NOTCH1 plays a key role in the differentiation of ventricles, and mutations are strongly associated with both sporadic and familial bicuspid aortic valves. However, few heart-specific enhancers have been identified to date.

RESULTS

In this study, we investigated evolutionary conserved regions (ECRs) that might act as potential enhancers within the region approximately 150-kb upstream and downstream of the NOTCH1 gene. Functional validation revealed that one 127-bp ECR located ~85-kb downstream of the NOTCH1 gene drives green fluorescent protein expression in the zebrafish embryo heart. Transcription factor (TF) prediction and core TF distribution analyses were performed to identify the core region. Dissection of ECR3 was performed to identify the 42-bp sequence, which is sufficient for heart-specific expression. In situ hybridization experiments showed that notch1b is expressed in the heart. Overexpression experiments in cells indicated that NKX2-5 is critical for enhancer activity. Mutation of the NKX-5 binding site significantly decreased reporter gene expression. Next, compared with the commonly used myocardium-labeled zebrafish transgenic strain Tg(cmlc2: mCherry), this 42-bp enhancer-labeled stable line mediated a similar expression pattern but with a smaller core region.

CONCLUSION

This study identified a 42-bp heart-specific enhancer near the NOTCH1 gene and further verified its functional targeting by NKX2-5.

Ancient duons may underpin spatial patterning of gene expression in C4 leaves

If the highly efficient C4 photosynthesis pathway could be transferred to crops with the C3 pathway there could be yield gains of up to 50%. It has been proposed that the multiple metabolic and developmental modifications associated with C4 photosynthesis are underpinned by relatively few master regulators that have allowed the evolution of C4 photosynthesis more than 60 times in flowering plants. Here we identify a component of one such regulator that consists of a pair of cis-elements located in coding sequence of multiple genes that are preferentially expressed in bundle sheath cells of C4 leaves. These motifs represent duons as they play a dual role in coding for amino acids as well as controlling the spatial patterning of gene expression associated with the C4 leaf. They act to repress transcription of C4 photosynthesis genes in mesophyll cells. These duons are also present in the C3 model Arabidopsis thaliana, and, in fact, are conserved in all land plants and even some algae that use C3 photosynthesis. C4 photosynthesis therefore appears to have coopted an ancient regulatory code to generate the spatial patterning of gene expression that is a hallmark of C4 photosynthesis. This intragenic transcriptional regulatory sequence could be exploited in the engineering of efficient photosynthesis of crops.

Conserved non-coding elements: developmental gene regulation meets genome organization

Comparative genomics has revealed a class of non-protein-coding genomic sequences that display an extraordinary degree of conservation between two or more organisms, regularly exceeding that found within protein-coding exons. These elements, collectively referred to as conserved non-coding elements (CNEs), are non-randomly distributed across chromosomes and tend to cluster in the vicinity of genes with regulatory roles in multicellular development and differentiation. CNEs are organized into functional ensembles called genomic regulatory blocks–dense clusters of elements that collectively coordinate the expression of shared target genes, and whose span in many cases coincides with topologically associated domains. CNEs display sequence properties that set them apart from other sequences under constraint, and have recently been proposed as useful markers for the reconstruction of the evolutionary history of organisms. Disruption of several of these elements is known to contribute to diseases linked with development, and cancer. The emergence, evolutionary dynamics and functions of CNEs still remain poorly understood, and new approaches are required to enable comprehensive CNE identification and characterization. Here, we review current knowledge and identify challenges that need to be tackled to resolve the impasse in understanding extreme non-coding conservation.

Population Genomics of Paramecium Species

Population-genomic analyses are essential to understanding factors shaping genomic variation and lineage-specific sequence constraints. The dearth of such analyses for unicellular eukaryotes prompted us to assess genomic variation in Paramecium, one of the most well-studied ciliate genera. The Paramecium aurelia complex consists of ∼15 morphologically indistinguishable species that diverged subsequent to two rounds of whole-genome duplications (WGDs, as long as 320 MYA) and possess extremely streamlined genomes. We examine patterns of both nuclear and mitochondrial polymorphism, by sequencing whole genomes of 10-13 worldwide isolates of each of three species belonging to the P. aurelia complex: P. tetraurelia, P. biaurelia, P. sexaurelia, as well as two outgroup species that do not share the WGDs: P. caudatum and P. multimicronucleatum. An apparent absence of global geographic population structure suggests continuous or recent dispersal of Paramecium over long distances. Intergenic regions are highly constrained relative to coding sequences, especially in P. caudatum and P. multimicronucleatum that have shorter intergenic distances. Sequence diversity and divergence are reduced up to ∼100-150 bp both upstream and downstream of genes, suggesting strong constraints imposed by the presence of densely packed regulatory modules. In addition, comparison of sequence variation at non-synonymous and synonymous sites suggests similar recent selective pressures on paralogs within and orthologs across the deeply diverging species. This study presents the first genome-wide population-genomic analysis in ciliates and provides a valuable resource for future studies in evolutionary and functional genetics in Paramecium.

Association of HOMER1 rs2290639 with suicide attempts in Hong Kong Chinese and the potentially functional role of this polymorphism

OBJECTIVE

Animal evidence and genetic studies suggest that HOMER1 (homer homolog 1) is involved in the etiology of suicidal behavior and major depression disorder (MDD). However, most of genetic studies were performed in Caucasians and the potentially functional role of associated polymorphisms in HOMER1 was seldom reported. The purpose of this study was to investigate the association of a HOMER1 polymorphism rs2290639 with suicide attempts (SA) and MDD in Hong Kong Chinese, and then briefly elucidate the potentially functional role of the associated polymorphism.

METHODS

NEO personality inventory, impulsiveness and depression rating scales were completed by the subjects. The association studies of HOMER1 rs2290639 with SA or MDD were performed by case-control association studies. The bioinformatics analyses were adapted to predict potential transcription factors binding sites for the associated polymorphism.

RESULTS

The association studies and meta-analysis suggested that the HOMER1 rs2290639 was significantly associated with susceptibility to SA but seemed not to be associated with MDD in Hong Kong Chinese. This polymorphism might affect the transcription of the HOMER1 gene through interacting with a reliable transcription factor as found by three of four bioinformatics tools. In addition, close correlations between impulsiveness and NEO personality five factors were found in SA and MDD patients, which provide a possible way to assess the impulsiveness of patients through subjects' personality profiles for Hong Kong Chinese.

CONCLUSIONS

The HOMER1 rs2290639 polymorphism was significantly associated with susceptibility to SA in Hong Kong Chinese affected by psychiatric disorders, which might be explained by the potentially functional role of this polymorphism.

Exonic enhancers: proceed with caution in exome and genome sequencing studies

Exonic enhancers (eExons) are coding exons that also function as enhancers of the gene in which they reside or (a) nearby gene(s). Mutations that affect the enhancer activity of these eExons have been associated with human disease. Therefore, eExon mutations should be taken into account in exome and genome sequencing projects, not only because of the ability of these mutations to modify the encoded proteins but also because of their effects on enhancer activity.

Favorable genomic environments for cis-regulatory evolution: A novel theoretical framework

Cis-regulatory changes are arguably the primary evolutionary source of animal morphological diversity. With the recent explosion of genome-wide comparisons of the cis-regulatory content in different animal species is now possible to infer general principles underlying enhancer evolution. However, these studies have also revealed numerous discrepancies and paradoxes, suggesting that the mechanistic causes and modes of cis-regulatory evolution are still not well understood and are probably much more complex than generally appreciated. Here, we argue that the mutational mechanisms and genomic regions generating new regulatory activities must comply with the constraints imposed by the molecular properties of cis-regulatory elements (CREs) and the organizational features of long-range chromatin interactions. Accordingly, we propose a new integrative evolutionary framework for cis-regulatory evolution based on two major premises for the origin of novel enhancer activity: (i) an accessible chromatin environment and (ii) compatibility with the 3D structure and interactions of pre-existing CREs. Mechanisms and DNA sequences not fulfilling these premises, will be less likely to have a measurable impact on gene expression and as such, will have a minor contribution to the evolution of gene regulation. Finally, we discuss current comparative cis-regulatory data under the light of this new evolutionary model, and propose that the two most prominent mechanisms for the evolution of cis-regulatory changes are the overprinting of ancestral CREs and the exaptation of transposable elements.

Individual genome sequencing identified a novel enhancer element in exon 7 of the CSFR1 gene by shift of expressed allele ratios

The sequencing of individual genetic information may provide a powerful tool for elucidating the mechanism by which individual SNPs affect promoter function. Here, we assessed the genome of a Russian male that was previously sequenced. The RNA-Seq data from blood cells revealed 234 candidate transcripts with shifts of greater than 1.5-fold from equal biallelic transcription. Of these genes, the CSF1R gene had variations in genic regions that affected the association of RORalpha with its target binding site in vivo. The results of a reporter assay confirmed that a single nucleotide substitution, rs2228422, within the RORalpha recognition motif altered the ability of the enhancer to regulate CSF1R gene transcription. Notably, 31% of Europeans and only 3% of Asians are homozygous for a RORalpha responsive "A" allele, but no association with diseases of rs2228422 has been found thus far.

Systematic dissection of coding exons at single nucleotide resolution supports an additional role in cell-specific transcriptional regulation

In addition to their protein coding function, exons can also serve as transcriptional enhancers. Mutations in these exonic-enhancers (eExons) could alter both protein function and transcription. However, the functional consequence of eExon mutations is not well known. Here, using massively parallel reporter assays, we dissect the enhancer activity of three liver eExons (SORL1 exon 17, TRAF3IP2 exon 2, PPARG exon 6) at single nucleotide resolution in the mouse liver. We find that both synonymous and non-synonymous mutations have similar effects on enhancer activity and many of the deleterious mutation clusters overlap known liver-associated transcription factor binding sites. Carrying a similar massively parallel reporter assay in HeLa cells with these three eExons found differences in their mutation profiles compared to the liver, suggesting that enhancers could have distinct operating profiles in different tissues. Our results demonstrate that eExon mutations could lead to multiple phenotypes by disrupting both the protein sequence and enhancer activity and that enhancers can have distinct mutation profiles in different cell types.

Exons that code for protein can also have additional functions, such as regulating gene transcription through enhancer activity. Here, we changed every nucleotide in three different exons that also function as enhancers, and examined their enhancer activity to test whether nucleotide changes in these exons can affect both the protein sequence and enhancer function. We found that mutations with a significant effect on enhancer function can reside both in regions that change the protein sequence (non-synonymous) and regions that do not change it (synonymous). When we conducted a similar analysis in a different cell type, we observed a difference in the nucleotide changes that cause a significant effect on enhancer activity, suggesting that the enhancer functional units can differ between tissues.

Computational characterization of modes of transcriptional regulation of nuclear receptor genes

Background

Nuclear receptors are a large structural class of transcription factors that act with their co-regulators and repressors to maintain a variety of biological and physiological processes such as metabolism, development and reproduction. They are activated through the binding of small ligands, which can be replaced by drug molecules, making nuclear receptors promising drug targets. Transcriptional regulation of the genes that encode them is central to gaining a deeper understanding of the diversity of their biochemical and biophysical roles and their role in disease and therapy. Even though they share evolutionary history, nuclear receptor genes have fundamentally different expression patterns, ranging from ubiquitously expressed to tissue-specific and spatiotemporally complex. However, current understanding of regulation in nuclear receptor gene family is still nascent.

Methodology/Principal Findings

In this study, we investigate the relationship between long-range regulation of nuclear receptor family and their known functionality. Towards this goal, we identify the nuclear receptor genes that are potential targets based on counts of highly conserved non-coding elements. We validate our results using publicly available expression (RNA-seq) and histone modification (ChIP-seq) data from the ENCODE project. We find that nuclear receptor genes involved in developmental roles show strong evidence of long-range mechanism of transcription regulation with distinct cis-regulatory content they feature clusters of highly conserved non-coding elements distributed in regions spanning several Megabases, long and multiple CpG islands, bivalent promoter marks and statistically significant higher enrichment of enhancer mark around their gene loci. On the other hand nuclear receptor genes that are involved in tissue-specific roles lack these features, having simple transcriptional controls and a greater variety of mechanisms for producing paralogs. We further examine the combinatorial patterns of histone maps associated with dynamic functional elements in order to explore the regulatory landscape of the gene family. The results show that our proposed classification capturing long-range regulation is strongly indicative of the functional roles of the nuclear receptors compared to existing classifications.

Conclusions/Significanc

We present a new classification for nuclear receptor gene family capturing whether a nuclear receptor is a possible target of long-range regulation or not. We compare our classification to existing structural (mechanism of action) and homology-based classifications. Our results show that understanding long-range regulation of nuclear receptors can provide key insight into their functional roles as well as evolutionary history; and this strongly merits further study.

Comparison of ultra-conserved elements in drosophilids and vertebrates

Metazoan genomes contain many ultra-conserved elements (UCEs), long sequences identical between distant species. In this study we identified UCEs in drosophilid and vertebrate species with a similar level of phylogenetic divergence measured at protein-coding regions, and demonstrated that both the length and number of UCEs are larger in vertebrates. The proportion of non-exonic UCEs declines in distant drosophilids whilst an opposite trend was observed in vertebrates. We generated a set of 2,126 Sophophora UCEs by merging elements identified in several drosophila species and compared these to the eutherian UCEs identified in placental mammals. In contrast to vertebrates, the Sophophora UCEs are depleted around transcription start sites. Analysis of 52,954 P-element, piggyBac and Minos insertions in the D. melanogaster genome revealed depletion of the P-element and piggyBac insertions in and around the Sophophora UCEs. We examined eleven fly strains with transposon insertions into the intergenic UCEs and identified associated phenotypes in five strains. Four insertions behave as recessive lethals, and in one case we observed a suppression of the marker gene within the transgene, presumably by silenced chromatin around the integration site. To confirm the lethality is caused by integration of transposons we performed a phenotype rescue experiment for two stocks and demonstrated that the excision of the transposons from the intergenic UCEs restores viability. Sequencing of DNA after the transposon excision in one fly strain with the restored viability revealed a 47 bp insertion at the original transposon integration site suggesting that the nature of the mutation is important for the appearance of the phenotype. Our results suggest that the UCEs in flies and vertebrates have both common and distinct features, and demonstrate that a significant proportion of intergenic drosophila UCEs are sensitive to disruption.

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.

Chromatin and epigenetic features of long-range gene regulation

The precise regulation of gene transcription during metazoan development is controlled by a complex system of interactions between transcription factors, histone modifications and modifying enzymes and chromatin conformation. Developments in chromosome conformation capture technologies have revealed that interactions between regions of chromatin are pervasive and highly cell-type specific. The movement of enhancers and promoters in and out of higher-order chromatin structures within the nucleus are associated with changes in expression and histone modifications. However, the factors responsible for mediating these changes and determining enhancer:promoter specificity are still not completely known. In this review, we summarize what is known about the patterns of epigenetic and chromatin features characteristic of elements involved in long-range interactions. In addition, we review the insights into both local and global patterns of chromatin interactions that have been revealed by the latest experimental and computational methods.

Highly conserved elements discovered in vertebrates are present in non-syntenic loci of tunicates, act as enhancers and can be transcribed during development

Co-option of cis-regulatory modules has been suggested as a mechanism for the evolution of expression sites during development. However, the extent and mechanisms involved in mobilization of cis-regulatory modules remains elusive. To trace the history of non-coding elements, which may represent candidate ancestral cis-regulatory modules affirmed during chordate evolution, we have searched for conserved elements in tunicate and vertebrate (Olfactores) genomes. We identified, for the first time, 183 non-coding sequences that are highly conserved between the two groups. Our results show that all but one element are conserved in non-syntenic regions between vertebrate and tunicate genomes, while being syntenic among vertebrates. Nevertheless, in all the groups, they are significantly associated with transcription factors showing specific functions fundamental to animal development, such as multicellular organism development and sequence-specific DNA binding. The majority of these regions map onto ultraconserved elements and we demonstrate that they can act as functional enhancers within the organism of origin, as well as in cross-transgenesis experiments, and that they are transcribed in extant species of Olfactores. We refer to the elements as 'Olfactores conserved non-coding elements'.

Hundreds of conserved non-coding genomic regions are independently lost in mammals

Conserved non-protein-coding DNA elements (CNEs) often encode cis-regulatory elements and are rarely lost during evolution. However, CNE losses that do occur can be associated with phenotypic changes, exemplified by pelvic spine loss in sticklebacks. Using a computational strategy to detect complete loss of CNEs in mammalian genomes while strictly controlling for artifacts, we find >600 CNEs that are independently lost in at least two mammalian lineages, including a spinal cord enhancer near GDF11. We observed several genomic regions where multiple independent CNE loss events happened; the most extreme is the DIAPH2 locus. We show that CNE losses often involve deletions and that CNE loss frequencies are non-uniform. Similar to less pleiotropic enhancers, we find that independently lost CNEs are shorter, slightly less constrained and evolutionarily younger than CNEs without detected losses. This suggests that independently lost CNEs are less pleiotropic and that pleiotropic constraints contribute to non-uniform CNE loss frequencies. We also detected 35 CNEs that are independently lost in the human lineage and in other mammals. Our study uncovers an interesting aspect of the evolution of functional DNA in mammalian genomes. Experiments are necessary to test if these independently lost CNEs are associated with parallel phenotype changes in mammals.

Discovery and characterization of human exonic transcriptional regulatory elements

We sought exonic transcriptional regulatory elements by shotgun cloning human cDNA fragments into luciferase reporter vectors and measuring the resulting expression levels in liver cells. We uncovered seven regulatory elements within coding regions and three within 3' untranslated regions (UTRs). Two of the putative regulatory elements were enhancers and eight were silencers. The regulatory elements were generally but not consistently evolutionarily conserved and also showed a trend toward decreased population diversity. Furthermore, the exonic regulatory elements were enriched in known transcription factor binding sites (TFBSs) and were associated with several histone modifications and transcriptionally relevant chromatin. Evidence was obtained for bidirectional cis-regulation of a coding region element within a tubulin gene, TUBA1B, by the transcription factors PPARA and RORA. We estimate that hundreds of exonic transcriptional regulatory elements exist, an unexpected finding that highlights a surprising multi-functionality of sequences in the human genome.

Polyploidy and the evolution of complex traits

We explore how whole-genome duplications (WGDs) may have given rise to complex innovations in cellular networks, innovations that could not have evolved through sequential single-gene duplications. We focus on two classical WGD events, one in bakers' yeast and the other at the base of vertebrates (i.e., two rounds of whole-genome duplication: 2R-WGD). Two complex adaptations are discussed in detail: aerobic ethanol fermentation in yeast and the rewiring of the vertebrate developmental regulatory network through the 2R-WGD. These two examples, derived from diverged branches on the eukaryotic tree, boldly underline the evolutionary potential of WGD in facilitating major evolutionary transitions. We close by arguing that the evolutionary importance of WGD may require updating certain aspects of modern evolutionary theory, perhaps helping to synthesize a new evolutionary systems biology.

Coding exons function as tissue-specific enhancers of nearby genes

Enhancers are essential gene regulatory elements whose alteration can lead to morphological differences between species, developmental abnormalities, and human disease. Current strategies to identify enhancers focus primarily on noncoding sequences and tend to exclude protein coding sequences. Here, we analyzed 25 available ChIP-seq data sets that identify enhancers in an unbiased manner (H3K4me1, H3K27ac, and EP300) for peaks that overlap exons. We find that, on average, 7% of all ChIP-seq peaks overlap coding exons (after excluding for peaks that overlap with first exons). By using mouse and zebrafish enhancer assays, we demonstrate that several of these exonic enhancer (eExons) candidates can function as enhancers of their neighboring genes and that the exonic sequence is necessary for enhancer activity. Using ChIP, 3C, and DNA FISH, we further show that one of these exonic limb enhancers, Dync1i1 exon 15, has active enhancer marks and physically interacts with Dlx5/6 promoter regions 900 kb away. In addition, its removal by chromosomal abnormalities in humans could cause split hand and foot malformation 1 (SHFM1), a disorder associated with DLX5/6. These results demonstrate that DNA sequences can have a dual function, operating as coding exons in one tissue and enhancers of nearby gene(s) in another tissue, suggesting that phenotypes resulting from coding mutations could be caused not only by protein alteration but also by disrupting the regulation of another gene.

Transcriptional enhancers in protein-coding exons of vertebrate developmental genes

Many conserved noncoding sequences function as transcriptional enhancers that regulate gene expression. Here, we report that protein-coding DNA also frequently contains enhancers functioning at the transcriptional level. We tested the enhancer activity of 31 protein-coding exons, which we chose based on strong sequence conservation between zebrafish and human, and occurrence in developmental genes, using a Tol2 transposable GFP reporter assay in zebrafish. For each exon we measured GFP expression in hundreds of embryos in 10 anatomies via a novel system that implements the voice-recognition capabilities of a cellular phone. We find that 24/31 (77%) exons drive GFP expression compared to a minimal promoter control, and 14/24 are anatomy-specific (expression in four anatomies or less). GFP expression driven by these coding enhancers frequently overlaps the anatomies where the host gene is expressed (60%), suggesting self-regulation. Highly conserved coding sequences and highly conserved noncoding sequences do not significantly differ in enhancer activity (coding: 24/31 vs. noncoding: 105/147) or tissue-specificity (coding: 14/24 vs. noncoding: 50/105). Furthermore, coding and noncoding enhancers display similar levels of the enhancer-related histone modification H3K4me1 (coding: 9/24 vs noncoding: 34/81). Meanwhile, coding enhancers are over three times as likely to contain an H3K4me1 mark as other exons of the host gene. Our work suggests that developmental transcriptional enhancers do not discriminate between coding and noncoding DNA and reveals widespread dual functions in protein-coding DNA.

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.

Making enhancers from spare parts of the genome

New studies show that novel long-range enhancers of developmental genes can emerge by exaptation of protein-coding sequences with no previous regulatory function.

De novo genesis of enhancers in vertebrates

Whole genome duplication in teleost fish reveals that a few changes in non-regulatory genomic sequences are a source for generating new enhancers.

Evolutionary innovation relies partially on changes in gene regulation. While a growing body of evidence demonstrates that such innovation is generated by functional changes or translocation of regulatory elements via mobile genetic elements, the de novo generation of enhancers from non-regulatory/non-mobile sequences has, to our knowledge, not previously been demonstrated. Here we show evidence for the de novo genesis of enhancers in vertebrates. For this, we took advantage of the massive gene loss following the last whole genome duplication in teleosts to systematically identify regions that have lost their coding capacity but retain sequence conservation with mammals. We found that these regions show enhancer activity while the orthologous coding regions have no regulatory activity. These results demonstrate that these enhancers have been de novo generated in fish. By revealing that minor changes in non-regulatory sequences are sufficient to generate new enhancers, our study highlights an important playground for creating new regulatory variability and evolutionary innovation.

The genome of each living organism contains thousands of genes, and the precise control of the timing and location of expression of these genes is key for normal development and homeostasis of each individual. Despite the oftentimes high genetic similarity between organisms, the source of phenotypic differences, for example between human and mouse, is thought to originate mainly from changes in how and when genes are expressed. This is partially determined by enhancers, that contribute to the control of gene expression. For decades, duplication of existing genomic enhancers, mobile elements, and changes in the sequence of existing enhancers were believed to be the major ways of increasing the number and modifying the activity of enhancers. In this study, we show that enhancers don't have to be derived from pre-existing ones but can also appear de novo in regions of the genome that were previously not regulating gene expression. We analyzed teleost fish genomes and found three regions for which a limited number of changes in the DNA sequence was sufficient to generate new enhancers. We predict that such a process is frequent in vertebrate genomes, making de novo generation of enhancers an important mechanism for creating variation in gene expression.

Locating protein-coding sequences under selection for additional, overlapping functions in 29 mammalian genomes

The degeneracy of the genetic code allows protein-coding DNA and RNA sequences to simultaneously encode additional, overlapping functional elements. A sequence in which both protein-coding and additional overlapping functions have evolved under purifying selection should show increased evolutionary conservation compared to typical protein-coding genes--especially at synonymous sites. In this study, we use genome alignments of 29 placental mammals to systematically locate short regions within human ORFs that show conspicuously low estimated rates of synonymous substitution across these species. The 29-species alignment provides statistical power to locate more than 10,000 such regions with resolution down to nine-codon windows, which are found within more than a quarter of all human protein-coding genes and contain ∼2% of their synonymous sites. We collect numerous lines of evidence that the observed synonymous constraint in these regions reflects selection on overlapping functional elements including splicing regulatory elements, dual-coding genes, RNA secondary structures, microRNA target sites, and developmental enhancers. Our results show that overlapping functional elements are common in mammalian genes, despite the vast genomic landscape.

Identification of cis-regulatory sequence variations in individual genome sequences

Functional contributions of cis-regulatory sequence variations to human genetic disease are numerous. For instance, disrupting variations in a HNF4A transcription factor binding site upstream of the Factor IX gene contributes causally to hemophilia B Leyden. Although clinical genome sequence analysis currently focuses on the identification of protein-altering variation, the impact of cis-regulatory mutations can be similarly strong. New technologies are now enabling genome sequencing beyond exomes, revealing variation across the non-coding 98% of the genome responsible for developmental and physiological patterns of gene activity. The capacity to identify causal regulatory mutations is improving, but predicting functional changes in regulatory DNA sequences remains a great challenge. Here we explore the existing methods and software for prediction of functional variation situated in the cis-regulatory sequences governing gene transcription and RNA processing.

Genomic selection identifies vertebrate transcription factor Fezf2 binding sites and target genes

Identification of transcription factor targets is critical to understanding gene regulatory networks. Here, we uncover transcription factor binding sites and target genes employing systematic evolution of ligands by exponential enrichment (SELEX). Instead of selecting randomly synthesized DNA oligonucleotides as in most SELEX studies, we utilized zebrafish genomic DNA to isolate fragments bound by Fezf2, an evolutionarily conserved gene critical for vertebrate forebrain development. This is, to our knowledge, the first time that SELEX is applied to a vertebrate genome. Computational analysis of bound genomic fragments predicted a core consensus binding site, which identified response elements that mediated Fezf2-dependent transcription both in vitro and in vivo. Fezf2-bound fragments were enriched for conserved sequences. Surprisingly, ∼20% of these fragments overlapped well annotated protein-coding exons. Through loss of function, gain of function, and chromatin immunoprecipitation, we further identified and validated eomesa/tbr2 and lhx2b as biologically relevant target genes of Fezf2. Mutations in eomesa/tbr2 cause microcephaly in humans, whereas lhx2b is a critical regulator of cell fate and axonal targeting in the developing forebrain. These results demonstrate the feasibility of employing genomic SELEX to identify vertebrate transcription factor binding sites and target genes and reveal Fezf2 as a transcription activator and a candidate for evaluation in human microcephaly.

DNA shape, genetic codes, and evolution

Although the three-letter genetic code that maps nucleotide sequence to protein sequence is well known, there must exist other codes that are embedded in the human genome. Recent work points to sequence-dependent variation in DNA shape as one mechanism by which regulatory and other information could be encoded in DNA. Recent advances include the discovery of shape-dependent recognition of DNA that depends on minor groove width and electrostatics, the existence of overlapping codes in protein-coding regions of the genome, and evolutionary selection for compensatory changes in nucleotide composition that facilitate nucleosome occupancy. It is becoming clear that DNA shape is important to biological function, and therefore will be subject to evolutionary constraint.

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.

The zebrafish: scalable in vivo modeling for systems biology

The zebrafish offers a scalable vertebrate model for many areas of biologic investigation. There is substantial conservation of genetic and genomic features and, at a higher order, conservation of intermolecular networks, as well as physiologic systems and phenotypes. We highlight recent work demonstrating the extent of this homology, and efforts to develop high-throughput phenotyping strategies suited to genetic or chemical screening on a scale compatible with in vivo validation for systems biology. We discuss the implications of these approaches for functional annotation of the genome, elucidation of multicellular processes in vivo, and mechanistic exploration of hypotheses generated by a broad range of 'unbiased' 'omic technologies such as expression profiling and genome-wide association. Finally, we outline potential strategies for the application of the zebrafish to the systematic study of phenotypic architecture, disease heterogeneity and drug responses.