A regulatory module embedded in the coding region of Hoxa2 controls expression in rhombomere 2

A multiple super-enhancer region establishes inter-TAD interactions and controls Hoxa function in cranial neural crest

Enhancer-promoter interactions preferentially occur within boundary-insulated topologically associating domains (TADs), limiting inter-TAD interactions. Enhancer clusters in linear proximity, termed super-enhancers (SEs), ensure high target gene expression levels. Little is known about SE topological regulatory impact during craniofacial development. Here, we identify 2232 genome-wide putative SEs in mouse cranial neural crest cells (CNCCs), 147 of which target genes establishing CNCC positional identity during face formation. In second pharyngeal arch (PA2) CNCCs, a multiple SE-containing region, partitioned into Hoxa Inter-TAD Regulatory Element 1 and 2 (HIRE1 and HIRE2), establishes long-range inter-TAD interactions selectively with Hoxa2, that is required for external and middle ear structures. HIRE2 deletion in a Hoxa2 haploinsufficient background results in microtia. HIRE1 deletion phenocopies the full homeotic Hoxa2 knockout phenotype and induces PA3 and PA4 CNCC abnormalities correlating with Hoxa2 and Hoxa3 transcriptional downregulation. Thus, SEs can overcome TAD insulation and regulate anterior Hoxa gene collinear expression in a CNCC subpopulation-specific manner during craniofacial development.

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.

Transcriptional Regulation and Implications for Controlling Hox Gene Expression

Hox genes play key roles in axial patterning and regulating the regional identity of cells and tissues in a wide variety of animals from invertebrates to vertebrates. Nested domains of Hox expression generate a combinatorial code that provides a molecular framework for specifying the properties of tissues along the A–P axis. Hence, it is important to understand the regulatory mechanisms that coordinately control the precise patterns of the transcription of clustered Hox genes required for their roles in development. New insights are emerging about the dynamics and molecular mechanisms governing transcriptional regulation, and there is interest in understanding how these may play a role in contributing to the regulation of the expression of the clustered Hox genes. In this review, we summarize some of the recent findings, ideas and emerging mechanisms underlying the regulation of transcription in general and consider how they may be relevant to understanding the transcriptional regulation of Hox genes.

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.

HOX gene cluster (de)regulation in brain: from neurodevelopment to malignant glial tumours

HOX genes encode a family of evolutionarily conserved homeodomain transcription factors that are crucial both during development and adult life. In humans, 39 HOX genes are arranged in four clusters (HOXA, B, C, and D) in chromosomes 7, 17, 12, and 2, respectively. During embryonic development, particular epigenetic states accompany their expression along the anterior-posterior body axis. This tightly regulated temporal-spatial expression pattern reflects their relative chromosomal localization, and is critical for normal embryonic brain development when HOX genes are mainly expressed in the hindbrain and mostly absent in the forebrain region. Epigenetic marks, mostly polycomb-associated, are dynamically regulated at HOX loci and regulatory regions to ensure the finely tuned HOX activation and repression, highlighting a crucial epigenetic plasticity necessary for homeostatic development. HOX genes are essentially absent in healthy adult brain, whereas they are detected in malignant brain tumours, namely gliomas, where HOX genes display critical roles by regulating several hallmarks of cancer. Here, we review the major mechanisms involved in HOX genes (de)regulation in the brain, from embryonic to adult stages, in physiological and oncologic conditions. We focus particularly on the emerging causes of HOX gene deregulation in glioma, as well as on their functional and clinical implications.

Oxaliplatin-induced neuropathic pain involves HOXA6 via a TET1-dependent demethylation of the SOX10 promoter

Chemotherapy-induced neuropathic pain is a common dose-limiting side effect of cancer treatment but the underlying mechanisms are largely unknown. Here, we used a whole-genome expression microarray and gene ontology analysis to identify the upregulation of a sequence-specific DNA-binding protein, HOXA6, in the spinal dorsal horn on Day 10 after injection of rats with oxaliplatin. Genetic disruption of HOXA6 with siRNAs alleviated mechanical allodynia after oxaliplatin administration. Reduced representation bisulfite sequencing assays indicated that oxaliplatin decreased the methylation levels of the SOX10 promoter but not of HOXA6. TET1 was also upregulated by oxaliplatin. Genetic disruption of TET1 with siRNA blocked the promoter demethylation of SOX10 and the upregulation of HOXA6 and SOX10. Importantly, inhibition of SOX10 by intrathecal application of SOX10 siRNA ameliorated the mechanical allodynia induced by oxaliplatin and downregulated the expression of HOXA6. Consistently, overexpression of SOX10 through intraspinal injection of AAV-SOX10-EGFP produced mechanical allodynia and upregulated the expression of spinal dorsal horn HOXA6. Moreover, chromatin immunoprecipitation assays demonstrated that oxaliplatin increased the binding of SOX10 to the promoter region of HOXA6. Taken together, our data suggest that HOXA6 upregulation through the TET1-mediated promoter demethylation of SOX10 may contribute to oxaliplatin-induced neuropathic pain.

Hindbrain induction and patterning during early vertebrate development

The hindbrain is a key relay hub of the central nervous system (CNS), linking the bilaterally symmetric half-sides of lower and upper CNS centers via an extensive network of neural pathways. Dedicated neural assemblies within the hindbrain control many physiological processes, including respiration, blood pressure, motor coordination and different sensations. During early development, the hindbrain forms metameric segmented units known as rhombomeres along the antero-posterior (AP) axis of the nervous system. These compartmentalized units are highly conserved during vertebrate evolution and act as the template for adult brainstem structure and function. TALE and HOX homeodomain family transcription factors play a key role in the initial induction of the hindbrain and its specification into rhombomeric cell fate identities along the AP axis. Signaling pathways, such as canonical-Wnt, FGF and retinoic acid, play multiple roles to initially induce the hindbrain and regulate Hox gene-family expression to control rhombomeric identity. Additional transcription factors including Krox20, Kreisler and others act both upstream and downstream to Hox genes, modulating their expression and protein activity. In this review, we will examine the earliest embryonic signaling pathways that induce the hindbrain and subsequent rhombomeric segmentation via Hox and other gene expression. We will examine how these signaling pathways and transcription factors interact to activate downstream targets that organize the segmented AP pattern of the embryonic vertebrate hindbrain.

Presynaptic Inhibition Selectively Gates Auditory Transmission to the Brainstem Startle Circuit

Filtering mechanisms prevent a continuous stream of sensory information from swamping perception, leading to diminished focal attention and cognitive processing. Mechanisms for sensory gating are commonly studied using prepulse inhibition, a paradigm that measures the regulated transmission of auditory information to the startle circuit; however, the underlying neuronal pathways are unresolved. Using large-scale calcium imaging, optogenetics, and laser ablations, we reveal a cluster of 30 morphologically identified neurons in zebrafish that suppress the transmission of auditory signals during prepulse inhibition. These neurons project to a key sensorimotor interface in the startle circuit-the termination zone of auditory afferents on the dendrite of a startle command neuron. Direct measurement of auditory nerve neurotransmitter release revealed selective presynaptic inhibition of sensory transmission to the startle circuit, sparing signaling to other brain regions. Our results provide the first cellular resolution circuit for prepulse inhibition in a vertebrate, revealing a central role for presynaptic gating of sensory information to a brainstem motor circuit.

Coupling the roles of Hox genes to regulatory networks patterning cranial neural crest

The neural crest is a transient population of cells that forms within the developing central nervous system and migrates away to generate a wide range of derivatives throughout the body during vertebrate embryogenesis. These cells are of evolutionary and clinical interest, constituting a key defining trait in the evolution of vertebrates and alterations in their development are implicated in a high proportion of birth defects and craniofacial abnormalities. In the hindbrain and the adjacent cranial neural crest cells (cNCCs), nested domains of Hox gene expression provide a combinatorial'Hox-code' for specifying regional properties in the developing head. Hox genes have been shown to play important roles at multiple stages in cNCC development, including specification, migration, and differentiation. However, relatively little is known about the underlying gene-regulatory mechanisms involved, both upstream and downstream of Hox genes. Furthermore, it is still an open question as to how the genes of the neural crest GRN are linked to Hox-dependent pathways. In this review, we describe Hox gene expression, function and regulation in cNCCs with a view to integrating these genes within the emerging gene regulatory network for cNCC development. We highlight early roles for Hox1 genes in cNCC specification, proposing that this may be achieved, in part, by regulation of the balance between pluripotency and differentiation in precursor cells within the neuro-epithelium. We then describe what is known about the regulation of Hox gene expression in cNCCs and discuss this from the perspective of early vertebrate evolution.

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.

Segmental arithmetic: summing up the Hox gene regulatory network for hindbrain development in chordates

Organization and development of the early vertebrate hindbrain are controlled by a cascade of regulatory interactions that govern the process of segmentation and patterning along the anterior-posterior axis via Hox genes. These interactions can be assembled into a gene regulatory network that provides a framework to interpret experimental data, generate hypotheses, and identify gaps in our understanding of the progressive process of hindbrain segmentation. The network can be broadly separated into a series of interconnected programs that govern early signaling, segmental subdivision, secondary signaling, segmentation, and ultimately specification of segmental identity. Hox genes play crucial roles in multiple programs within this network. Furthermore, the network reveals properties and principles that are likely to be general to other complex developmental systems. Data from vertebrate and invertebrate chordate models are shedding light on the origin and diversification of the network. Comprehensive cis-regulatory analyses of vertebrate Hox gene regulation have enabled powerful cross-species gene regulatory comparisons. Such an approach in the sea lamprey has revealed that the network mediating segmental Hox expression was present in ancestral vertebrates and has been maintained across diverse vertebrate lineages. Invertebrate chordates lack hindbrain segmentation but exhibit conservation of some aspects of the network, such as a role for retinoic acid in establishing nested Hox expression domains. These comparisons lead to a model in which early vertebrates underwent an elaboration of the network between anterior-posterior patterning and Hox gene expression, leading to the gene-regulatory programs for segmental subdivision and rhombomeric segmentation. WIREs Dev Biol 2017, 6:e286. doi: 10.1002/wdev.286 For further resources related to this article, please visit the WIREs website.

A Conserved Structural Signature of the Homeobox Coding DNA in HOX genes

The homeobox encodes a DNA-binding domain found in transcription factors regulating key developmental processes. The most notable examples of homeobox containing genes are the Hox genes, arranged on chromosomes in the same order as their expression domains along the body axis. The mechanisms responsible for the synchronous regulation of Hox genes and the molecular function of their colinearity remain unknown. Here we report the discovery of a conserved structural signature of the 180-base pair DNA fragment comprising the homeobox. We demonstrate that the homeobox DNA has a characteristic 3-base-pair periodicity in the hydroxyl radical cleavage pattern. This periodic pattern is significant in most of the 39 mammalian Hox genes and in other homeobox-containing transcription factors. The signature is present in segmented bilaterian animals as evolutionarily distant as humans and flies. It remains conserved despite the fact that it would be disrupted by synonymous mutations, which raises the possibility of evolutionary selective pressure acting on the structure of the coding DNA. The homeobox coding DNA may therefore have a secondary function, possibly as a regulatory element. The existence of such element may have important consequences for understanding how these genes are regulated.

Functional and Comparative Genomics of Hoxa2 Gene cis-Regulatory Elements: Evidence for Evolutionary Modification of Ancestral Core Element Activity

Hoxa2 is an evolutionarily conserved developmental regulatory gene that functions to specify rhombomere (r) and pharyngeal arch (PA) identities throughout the Osteichthyes. Japanese medaka (Oryzias latipes) hoxa2a, like orthologous Hoxa2 genes from other osteichthyans, is expressed during embryogenesis in r2–7 and PA2-7, whereas the paralogous medaka pseudogene, ψhoxa2b, is expressed in noncanonical Hoxa2 domains, including the pectoral fin buds. To understand the evolution of cis-regulatory element (CRE) control of gene expression, we conducted eGFP reporter gene expression studies with extensive functional mapping of several conserved CREs upstream of medaka hoxa2a and ψhoxa2b in transient and stable-line transgenic medaka embryos. The CREs tested were previously shown to contribute to directing mouse Hoxa2 gene expression in r3, r5, and PA2-4. Our results reveal the presence of sequence elements embedded in the medaka hoxa2a and ψhoxa2b upstream enhancer regions (UERs) that mediate expression in r4 and the PAs (hoxa2a r4/CNCC element) or in r3–7 and the PAs ψhoxa2b r3–7/CNCC element), respectively. Further, these elements were shown to be highly conserved among osteichthyans, which suggests that the r4 specifying element embedded in the UER of Hoxa2 is a deeply rooted rhombomere specifying element and the activity of this element has been modified by the evolution of flanking sequences that redirect its activity to alternative developmental compartments.

Uncovering drug-responsive regulatory elements

Nucleotide changes in gene regulatory elements can have a major effect on interindividual differences in drug response. For example, by reviewing all published pharmacogenomic genome-wide association studies, we show here that 96.4% of the associated single nucleotide polymorphisms reside in noncoding regions. We discuss how sequencing technologies are improving our ability to identify drug response-associated regulatory elements genome-wide and to annotate nucleotide variants within them. We highlight specific examples of how nucleotide changes in these elements can affect drug response and illustrate the techniques used to find them and functionally characterize them. Finally, we also discuss challenges in the field of drug-responsive regulatory elements that need to be considered in order to translate these findings into the clinic.

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.

The gene regulatory networks underlying formation of the auditory hindbrain

Development and evolution of auditory hindbrain nuclei are two major unsolved issues in hearing research. Recent characterization of transgenic mice identified the rhombomeric origins of mammalian auditory nuclei and unraveled genes involved in their formation. Here, we provide an overview on these data by assembling them into rhombomere-specific gene regulatory networks (GRNs), as they underlie developmental and evolutionary processes. To explore evolutionary mechanisms, we compare the GRNs operating in the mammalian auditory hindbrain with data available from the inner ear and other vertebrate groups. Finally, we propose that the availability of genomic sequences from all major vertebrate taxa and novel genetic techniques for non-model organisms provide an unprecedented opportunity to investigate development and evolution of the auditory hindbrain by comparative molecular approaches. The dissection of the molecular mechanisms leading to auditory structures will also provide an important framework for auditory processing disorders, a clinical problem difficult to tackle so far. These data will, therefore, foster basic and clinical hearing research alike.

Enhancer variants: evaluating functions in common disease

Gene enhancer elements are noncoding segments of DNA that play a central role in regulating transcriptional programs that control development, cell identity, and evolutionary processes. Recent studies have shown that noncoding single nucleotide polymorphisms (SNPs) that have been associated with risk for numerous common diseases through genome-wide association studies frequently lie in cell-type-specific enhancer elements. These enhancer variants probably influence transcriptional output, thereby offering a mechanistic basis to explain their association with risk for many common diseases. This review focuses on the identification and interpretation of disease-susceptibility variants that influence enhancer function. We discuss strategies for prioritizing the study of functional enhancer SNPs over those likely to be benign, review experimental and computational approaches to identifying the gene targets of enhancer variants, and highlight efforts to quantify the impact of enhancer variants on target transcript levels and cellular phenotypes. These studies are beginning to provide insights into the mechanistic basis of many common diseases, as well as into how we might translate this knowledge for improved disease diagnosis, prevention and treatments. Finally, we highlight five major challenges often associated with interpreting enhancer variants, and discuss recent technical advances that may help to surmount these challenges.

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.

A Hox regulatory network of hindbrain segmentation is conserved to the base of vertebrates

A defining feature governing head patterning of jawed vertebrates is a highly conserved gene regulatory network that integrates hindbrain segmentation with segmentally restricted domains of Hox gene expression. Although non-vertebrate chordates display nested domains of axial Hox expression, they lack hindbrain segmentation. The sea lamprey, a jawless fish, can provide unique insights into vertebrate origins owing to its phylogenetic position at the base of the vertebrate tree. It has been suggested that lamprey may represent an intermediate state where nested Hox expression has not been coupled to the process of hindbrain segmentation. However, little is known about the regulatory network underlying Hox expression in lamprey or its relationship to hindbrain segmentation. Here, using a novel tool that allows cross-species comparisons of regulatory elements between jawed and jawless vertebrates, we report deep conservation of both upstream regulators and segmental activity of enhancer elements across these distant species. Regulatory regions from diverse gnathostomes drive segmental reporter expression in the lamprey hindbrain and require the same transcriptional inputs (for example, Kreisler (also known as Mafba), Krox20 (also known as Egr2a)) in both lamprey and zebrafish. We find that lamprey hox genes display dynamic segmentally restricted domains of expression; we also isolated a conserved exonic hox2 enhancer from lamprey that drives segmental expression in rhombomeres 2 and 4. Our results show that coupling of Hox gene expression to segmentation of the hindbrain is an ancient trait with origin at the base of vertebrates that probably led to the formation of rhombomeric compartments with an underlying Hox code.

Sox2 acts as a rheostat of epithelial to mesenchymal transition during neural crest development

Precise control of self-renewal and differentiation of progenitor cells into the cranial neural crest (CNC) pool ensures proper head development, guided by signaling pathways such as BMPs, FGFs, Shh and Notch. Here, we show that murine Sox2 plays an essential role in controlling progenitor cell behavior during craniofacial development. A “Conditional by Inversion” Sox2 allele (Sox2^COIN) has been employed to generate an epiblast ablation of Sox2 function (Sox2^EpINV). Sox2^EpINV/+(H) haploinsufficient and conditional (Sox2^EpINV/mosaic) mutant embryos proceed beyond gastrulation and die around E11. These mutant embryos exhibit severe anterior malformations, with hydrocephaly and frontonasal truncations, which could be attributed to the deregulation of CNC progenitor cells during their epithelial to mesenchymal transition. This irregularity results in an exacerbated and aberrant migration of Sox10⁺ NCC in the branchial arches and frontonasal process of the Sox2 mutant embryos. These results suggest a novel role for Sox2 as a regulator of the epithelial to mesenchymal transitions (EMT) that are important for the cell flow in the developing head.

Synonymous constraint elements show a tendency to encode intrinsically disordered protein segments

Synonymous constraint elements (SCEs) are protein-coding genomic regions with very low synonymous mutation rates believed to carry additional, overlapping functions. Thousands of such potentially multi-functional elements were recently discovered by analyzing the levels and patterns of evolutionary conservation in human coding exons. These elements provide a good opportunity to improve our understanding of how the redundant nature of the genetic code is exploited in the cell. Our premise is that the protein segments encoded by such elements might better comply with the increased functional demands if they are structurally less constrained (i.e. intrinsically disordered). To test this idea, we investigated the protein segments encoded by SCEs with computational tools to describe the underlying structural properties. In addition to SCEs, we examined the level of disorder, secondary structure, and sequence complexity of protein regions overlapping with experimentally validated splice regulatory sites. We show that multi-functional gene regions translate into protein segments that are significantly enriched in structural disorder and compositional bias, while they are depleted in secondary structure and domain annotations compared to reference segments of similar lengths. This tendency suggests that relaxed protein structural constraints provide an advantage when accommodating multiple overlapping functions in coding regions.

Gene transcription in the zebrafish embryo: regulators and networks

The precise spatial and temporal control of gene expression is a key process in the development, maintenance and regeneration of the vertebrate body. A substantial proportion of vertebrate genomes encode genes that control the transcription of the genetic information into mRNA. The zebrafish is particularly well suited to investigate gene regulatory networks underlying the control of gene expression during development due to the external development of its transparent embryos and the increasingly sophisticated tools for genetic manipulation available for this model system. We review here recent data on the analysis of cis-regulatory modules, transcriptional regulators and their integration into gene regulatory networks in the zebrafish, using the developing spinal cord as example.

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.

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.

A high-resolution map of human evolutionary constraint using 29 mammals

Comparison of related genomes has emerged as a powerful lens for genome interpretation. Here, we report the sequencing and comparative analysis of 29 eutherian genomes. We confirm that at least 5.5% of the human genome has undergone purifying selection, and report constrained elements covering ~4.2% of the genome. We use evolutionary signatures and comparison with experimental datasets to suggest candidate functions for ~60% of constrained bases. These elements reveal a small number of new coding exons, candidate stop codon readthrough events, and over 10,000 regions of overlapping synonymous constraint within protein-coding exons. We find 220 candidate RNA structural families, and nearly a million elements overlapping potential promoter, enhancer and insulator regions. We report specific amino acid residues that have undergone positive selection, 280,000 non-coding elements exapted from mobile elements, and ~1,000 primate- and human-accelerated elements. Overlap with disease-associated variants suggests our findings will be relevant for studies of human biology and health.

Identification of two novel mutations in Shh long-range regulator associated with familial pre-axial polydactyly

Pre-axial polydactyly type II (PPDII, MIM #174500), Werner mesomelic syndrome (MIM %188770) and Haas polysyndactyly (MIM #186200) are a group of closely related conditions caused by mutations in a long-range Sonic hedgehog (SHH, MIM *600725) regulator called ZRS. To date, 19 point mutations, 10 duplications and 1 triplication of the ZRS associated with those pre-axial polydactylies have been reported in humans, mice, cats and chickens. Some of these have been shown to cause ectopic expression of Shh in the limb bud in mice, leading to a polydactylous phenotype, but the precise mechanism by which ZRS mutations generate this phenotype remains unknown. We present two PPDII families with fully penetrant point mutations in ultra-conserved predicted binding sites for transcription factors SOX9 and PAX3, two possible candidates for regulating SHH expression. Screening for point mutations or copy-number variation of the ZRS, high-resolution array-CGH, and screening of other conserved non-coding sequences (CNS) surrounding SHH in a third family are negative. This is the sixth PPDII pedigree with possible linkage to 7q36 that presents with no detectable ZRS mutation. We hypothesize that another nearby regulatory sequence, or an undetected position effect between ZRS and SHH, could be responsible for negative familial cases linked to 7q36.

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.

Regionalization of the shark hindbrain: a survey of an ancestral organization

Cartilaginous fishes (chondrichthyans) represent an ancient radiation of vertebrates currently considered the sister group of the group of gnathostomes with a bony skeleton that gave rise to land vertebrates. This out-group position makes chondrichthyans essential in assessing the ancestral organization of the brain of jawed vertebrates. To gain knowledge about hindbrain evolution we have studied its development in a shark, the lesser spotted dogfish Scyliorhinus canicula by analyzing the expression of some developmental genes and the origin and distribution of specific neuronal populations, which may help to identify hindbrain subdivisions and boundaries and the topology of specific cell groups. We have characterized three developmental periods that will serve as a framework to compare the development of different neuronal systems and may represent a suitable tool for comparing the absolute chronology of development among vertebrates. The expression patterns of Pax6, Wnt8, and HoxA2 genes in early embryos of S. canicula showed close correspondence to what has been described in other vertebrates and helped to identify the anterior rhombomeres. Also in these early embryos, the combination of Pax6 with protein markers of migrating neuroblasts (DCX) and early differentiating neurons (general: HuC/D; neuron type specific: GAD, the GABA synthesizing enzyme) revealed the organization of S. canicula hindbrain in both transverse segmental units corresponding to visible rhombomeres and longitudinal columns. Later in development, when the interrhombomeric boundaries fade away, accurate information about S. canicula hindbrain subdivisions was achieved by comparing the expression patterns of Pax6 and GAD, serotonin (serotoninergic neurons), tyrosine hydroxylase (catecholaminergic neurons), choline acetyltransferase (cholinergic neurons), and calretinin (a calcium-binding protein). The patterns observed revealed many topological correspondences with other vertebrates and led to reconsideration of the current view of the elasmobranch hindbrain segmentation as peculiar among vertebrates.

cis-regulatory mutations are a genetic cause of human limb malformations

The underlying mutations that cause human limb malformations are often difficult to determine, particularly for limb malformations that occur as isolated traits. Evidence from a variety of studies shows that cis-regulatory mutations, specifically in enhancers, can lead to some of these isolated limb malformations. Here, we provide a review of human limb malformations that have been shown to be caused by enhancer mutations and propose that cis-regulatory mutations will continue to be identified as the cause of additional human malformations as our understanding of regulatory sequences improves.

When needles look like hay: how to find tissue-specific enhancers in model organism genomes

A major prerequisite for the investigation of tissue-specific processes is the identification of cis-regulatory elements. No generally applicable technique is available to distinguish them from any other type of genomic non-coding sequence. Therefore, researchers often have to identify these elements by elaborate in vivo screens, testing individual regions until the right one is found. Here, based on many examples from the literature, we summarize how functional enhancers have been isolated from other elements in the genome and how they have been characterized in transgenic animals. Covering computational and experimental studies, we provide an overview of the global properties of cis-regulatory elements, like their specific interactions with promoters and target gene distances. We describe conserved non-coding elements (CNEs) and their internal structure, nucleotide composition, binding site clustering and overlap, with a special focus on developmental enhancers. Conflicting data and unresolved questions on the nature of these elements are highlighted. Our comprehensive overview of the experimental shortcuts that have been found in the different model organism communities and the new field of high-throughput assays should help during the preparation phase of a screen for enhancers. The review is accompanied by a list of general guidelines for such a project.

ChIPing the cistrome of PXR in mouse liver

The pregnane X receptor (PXR) is a key regulator of xenobiotic metabolism and disposition in liver. However, little is known about the PXR DNA-binding signatures in vivo, or how PXR regulates novel direct targets on a genome-wide scale. Therefore, we generated a roadmap of hepatic PXR bindings in the entire mouse genome [chromatin immunoprecipitation (ChIP)-Seq]. The most frequent PXR DNA-binding motif is the AGTTCA-like direct repeat with a 4 bp spacer [direct repeat (DR)-4)]. Surprisingly, there are also high motif occurrences with spacers of a periodicity of 5 bp, forming a novel DR-(5 n+4) pattern for PXR binding. PXR-binding overlaps with the epigenetic mark for gene activation (histone-H3K4-di-methylation), but not with epigenetic marks for gene suppression (DNA methylation or histone-H3K27-tri-methylation) (ChIP-on-chip). After administering a PXR agonist, changes in mRNA of most PXR-direct target genes correlate with increased PXR binding. Specifically, increased PXR binding triggers the trans-activation of critical drug-metabolizing enzymes and transporters. The mRNA induction of these genes is absent in PXR-null mice. The current work provides the first in vivo evidence of PXR DNA-binding signatures in the mouse genome, paving the path for predicting and further understanding the multifaceted roles of PXR in liver.

Proteome analysis and genome-wide regulatory motif prediction identify novel potentially sex-hormone regulated proteins in rat efferent ducts

The efferent ducts are a series of tubules that conduct sperm from the rete testis to the epididymis. They absorb most fluid and proteins originating from the rete testis during concentration of spermatozoa prior to their entry into the epididymis. Proteome analysis of micro-dissected efferent duct samples from adult rats was combined with genome-wide computational prediction of conserved hormone response elements to identify factors likely regulated by oestrogens and androgens. We identified 165 proteins and found subsets of the promoters controlling their corresponding genes to contain androgen- and oestrogen response elements (ARE/EREs) at similar frequencies. Moreover, EREs were significantly enriched among the loci identified compared with their genome-wide occurrence. The expression and localization of Anxa6, Ckb, Krt19, Park7, Pdzk1 and Tpt1 in the efferent ducts and other related hormone controlled tissues was further validated at the RNA or protein level. This study identifies many novel proteins predicted to play roles in sperm maturation and male fertility and provides significant computational evidence that the efferent ducts express genes transcriptionally controlled by sex hormones.

Molecular mechanisms of cranial neural crest cell migration and patterning in craniofacial development

During vertebrate craniofacial development, neural crest cells (NCCs) contribute much of the cartilage, bone and connective tissue that make up the developing head. Although the initial patterns of NCC segmentation and migration are conserved between species, the variety of vertebrate facial morphologies that exist indicates that a complex interplay occurs between intrinsic genetic NCC programs and extrinsic environmental signals during morphogenesis. Here, we review recent work that has begun to shed light on the molecular mechanisms that govern the spatiotemporal patterning of NCC-derived skeletal structures – advances that are central to understanding craniofacial development and its evolution.

Mouse and zebrafish Hoxa3 orthologues have nonequivalent in vivo protein function

Hox genes play evolutionarily conserved roles in specifying axial position during embryogenesis. A prevailing paradigm is that changes in Hox gene expression drive evolution of metazoan body plans. Conservation of Hox function across species, and among paralogous Hox genes within a species, supports a model of functional equivalence. In this report, we demonstrate that zebrafish hoxa3a (zfhoxa3a) expressed from the mouse Hoxa3 locus can substitute for mouse Hoxa3 in some tissues, but has distinct or null phenotypes in others. We further show, by using an allele encoding a chimeric protein, that this difference maps primarily to the zfhoxa3a C-terminal domain. Our data imply that the mouse and zebrafish proteins have diverged considerably since their last common ancestor, and that the major difference between them resides in the C-terminal domain. Our data further show that Hox protein function can evolve independently in different cell types or for specific functions. The inability of zfhoxa3a to perform all of the normal roles of mouse Hoxa3 illustrates that Hox orthologues are not always functionally interchangeable.

Conserved expression of mouse Six1 in the pre-placodal region (PPR) and identification of an enhancer for the rostral PPR

All cranial sensory organs and sensory neurons of vertebrates develop from cranial placodes. In chick, amphibians and zebrafish, all placodes originate from a common precursor domain, the pre-placodal region (PPR), marked by the expression of Six1/4 and Eya1/2. However, the PPR has never been described in mammals and the mechanism involved in the formation of PPR is poorly defined. Here, we report the expression of Six1 in the horseshoe-shaped mouse ectoderm surrounding the anterior neural plate in a pattern broadly similar to that of non-mammalian vertebrates. To elucidate the identity of Six1-positive mouse ectoderm, we searched for enhancers responsible for Six1 expression by in vivo enhancer assays. One conserved non-coding sequence, Six1-14, showed specific enhancer activity in the rostral PPR of chick and Xenopus and in the mouse ectoderm. These results strongly suggest the presence of PPR in mouse and that it is conserved in vertebrates. Moreover, we show the importance of the homeodomain protein-binding sites of Six1-14, the Six1 rostral PPR enhancer, for enhancer activity, and that Dlx5, Msx1 and Pax7 are candidate binding factors that regulate the level and area of Six1 expression, and thereby the location of the PPR. Our findings provide critical information and tools to elucidate the molecular mechanism of early sensory development and have implications for the development of sensory precursor/stem cells.

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

Using a comparative genomics approach to reconstruct the fate of genomic regulatory blocks (GRBs) and identify exonic remnants that have survived the disappearance of their host genes after whole-genome duplication (WGD) in teleosts, we discover a set of 38 candidate cis-regulatory coding exons (RCEs) with predicted target genes. These elements demonstrate evolutionary separation of overlapping protein-coding and regulatory information after WGD in teleosts. We present evidence that the corresponding mammalian exons are still under both coding and non-coding selection pressure, are more conserved than other protein coding exons in the host gene and several control sets, and share key characteristics with highly conserved non-coding elements in the same regions. Their dual function is corroborated by existing experimental data. Additionally, we show examples of human exon remnants stemming from the vertebrate 2R WGD. Our findings suggest that long-range cis-regulatory inputs for developmental genes are not limited to non-coding regions, but can also overlap the coding sequence of unrelated genes. Thus, exonic regulatory elements in GRBs might be functionally equivalent to those in non-coding regions, calling for a re-evaluation of the sequence space in which to look for long-range regulatory elements and experimentally test their activity.

A systems biology approach to understanding cis-regulatory module function

The genomic instructions used to regulate development are encoded within a set of functional DNA elements called cis-regulatory modules (CRMs). These elements determine the precise patterns of temporal and spatial gene expression. Here we summarize recent progress made towards cataloguing and characterizing the complete repertoire of CRMs. We describe CRMs as genomic information processing devices containing clusters of transcription factor binding sites and we position CRMs as nodes within large gene regulatory networks. We define CRM architecture and describe how these genomic elements process the information they encode to their target genes. Furthermore, we present an overview describing high-throughput techniques to identify CRMs genome wide and experimental methodologies to validate their function on a large scale. This review emphasizes the advantages and power of a systems biology approach which integrates computational and experimental technologies to further our understanding of CRM function.

Transcriptional coupling of synaptic transmission and energy metabolism: role of nuclear respiratory factor 1 in co-regulating neuronal nitric oxide synthase and cytochrome c oxidase genes in neurons

Neuronal activity is highly dependent on energy metabolism; yet, the two processes have traditionally been regarded as independently regulated at the transcriptional level. Recently, we found that the same transcription factor, nuclear respiratory factor 1 (NRF-1) co-regulates an important energy-generating enzyme, cytochrome c oxidase, as well as critical subunits of glutamatergic receptors. The present study tests our hypothesis that the co-regulation extends to the next level of glutamatergic synapses, namely, neuronal nitric oxide synthase, which generates nitric oxide as a downstream signaling molecule. Using in silico analysis, electrophoretic mobility shift assay, chromatin immunoprecipitation, promoter mutations, and NRF-1 silencing, we documented that NRF-1 functionally bound to Nos1, but not Nos2 (inducible) and Nos3 (endothelial) gene promoters. Both COX and Nos1 transcripts were up-regulated by depolarizing KCl treatment and down-regulated by TTX-mediated impulse blockade in neurons. However, NRF-1 silencing blocked the up-regulation of both Nos1 and COX induced by KCl depolarization, and over-expression of NRF-1 rescued both Nos1 and COX transcripts down-regulated by TTX. These findings are consistent with our hypothesis that synaptic neuronal transmission and energy metabolism are tightly coupled at the molecular level.

The G allele of SNP E1/A118G at the µ-opioid receptor gene locus shows genomic evidence of recent positive selection

Opioid drug response and pain perception differs greatly amongst different individuals. The µ-opioid receptor (MOR) is the main receptor target for important opioid analgesics. As SNPs may contribute to interindividual differences in drug response, in silico signatures of recent positive selection (RPS) were utilized to seek out potentially functional SNPs in the MOR gene in order to facilitate the prioritization of SNPs for evaluation in genetic association studies. Out of over 1000 SNPs at the MOR locus, 184 high-frequency SNPs were interrogated for signatures of RPS. A total of five SNPs (four noncoding and one nonsynonymous coding) demonstrated in silico evidence of RPS. Significantly, the nonsynonymous E1/A118G SNP, which was previously reported to be functionally important, showed in silico evidence of RPS. This reaffirms the feasibility of utilizing in silico signatures of RPS to identify potentially functionally significant SNPs for association studies. Interestingly, the positively selected G allele of this RPS SNP was also predicted to create a novel exon splice enhancer as well as p53 binding sites.

Conserved elements within open reading frames of mammalian Hox genes

A recent study in BMC Evolutionary Biology shows that many of the open reading frames in mammalian Hox genes are more conserved than expected on the basis of their protein sequence. The presence of highly conserved DNA elements is thus not confined to the noncoding DNA in neighboring regions but clearly overlaps with coding sequences. These findings support an emerging view that gene regulatory and coding sequences are likely to be more intermingled than once believed.