Metrics of sequence constraint overlook regulatory sequences in an exhaustive analysis at phox2b

A functional schizophrenia-associated genetic variant near the TSNARE1 and ADGRB1 genes

Recent collaborative genome-wide association studies (GWAS) have identified >200 independent loci contributing to risk for schizophrenia (SCZ). The genes closest to these loci have diverse functions, supporting the potential involvement of multiple relevant biological processes, yet there is no direct evidence that individual variants are functional or directly linked to specific genes. Nevertheless, overlap with certain epigenetic marks suggest that most GWAS-implicated variants are regulatory. Based on the strength of association with SCZ and the presence of regulatory epigenetic marks, we chose one such variant near TSNARE1 and ADGRB1, rs4129585, to test for functional potential and assay differences that may drive the pathogenicity of the risk allele. We observed that the variant-containing sequence drives reporter expression in relevant neuronal populations in zebrafish. Next, we introduced each allele into human induced pluripotent cells and differentiated four isogenic clones homozygous for the risk allele and five clones homozygous for the non-risk allele into neural progenitor cells. Employing RNA sequencing, we found that the two alleles yield significant transcriptional differences in the expression of 109 genes at a false discovery rate (FDR) of <0.05 and 259 genes at a FDR of <0.1. We demonstrate that these genes are highly interconnected in pathways enriched for synaptic proteins, axon guidance, and regulation of synapse assembly. Exploration of genes near rs4129585 suggests that this variant does not regulate TSNARE1 transcripts, as previously thought, but may regulate the neighboring ADGRB1, a regulator of synaptogenesis. Our results suggest that rs4129585 is a functional common variant that functions in specific pathways likely involved in SCZ risk.

A Functional Schizophrenia-associated genetic variant near the TSNARE1 and ADGRB1 genes

Recent collaborative genome wide association studies (GWAS) have identified >200 independent loci contributing to risk for schizophrenia (SCZ). The genes closest to these loci have diverse functions, supporting the potential involvement of multiple relevant biological processes; yet there is no direct evidence that individual variants are functional or directly linked to specific genes. Nevertheless, overlap with certain epigenetic marks suggest that most GWAS-implicated variants are regulatory. Based on the strength of association with SCZ and the presence of regulatory epigenetic marks, we chose one such variant near TSNARE1 and ADGRB1, rs4129585, to test for functional potential and assay differences that may drive the pathogenicity of the risk allele. We observed that the variant-containing sequence drives reporter expression in relevant neuronal populations in zebrafish. Next, we introduced each allele into human induced pluripotent cells and differentiated 4 isogenic clones homozygous for the risk allele and 5 clones homozygous for the non-risk allele into neural precursor cells. Employing RNA-seq, we found that the two alleles yield significant transcriptional differences in the expression of 109 genes at FDR <0.05 and 259 genes at FDR <0.1. We demonstrate that these genes are highly interconnected in pathways enriched for synaptic proteins, axon guidance, and regulation of synapse assembly. Exploration of genes near rs4129585 suggests that this variant does not regulate TSNARE1 transcripts, as previously thought, but may regulate the neighboring ADGRB1, a regulator of synaptogenesis. Our results suggest that rs4129585 is a functional common variant that functions in specific pathways likely involved in SCZ risk.

Conservation of peripheral nervous system formation mechanisms in divergent ascidian embryos

Ascidians with very similar embryos but highly divergent genomes are thought to have undergone extensive developmental system drift. We compared, in four species (Ciona and Phallusia for Phlebobranchia, Molgula and Halocynthia for Stolidobranchia), gene expression and gene regulation for a network of six transcription factors regulating peripheral nervous system (PNS) formation in Ciona. All genes, but one in Molgula, were expressed in the PNS with some differences correlating with phylogenetic distance. Cross-species transgenesis indicated strong levels of conservation, except in Molgula, in gene regulation despite lack of sequence conservation of the enhancers. Developmental system drift in ascidians is thus higher for gene regulation than for gene expression and is impacted not only by phylogenetic distance, but also in a clade-specific manner and unevenly within a network. Finally, considering that Molgula is divergent in our analyses, this suggests deep conservation of developmental mechanisms in ascidians after 390 My of separate evolution.

Hereditary hearing loss; about the known and the unknown

Hereditary hearing loss is both clinically and genetically very heterogeneous. Despite the large number of genes that have been associated with the condition, many cases remain unexplained. Novel gene associations with hearing loss are to be expected but also are defects of regulatory regions of the genome which are currently not routinely addressed in molecular genetic testing and research. Inheritance patterns other than monogenic might be more common than assumed in isolated cases and diagnoses might have been missed because of misinterpretation of identified DNA variants. This review summarizes current insights in the genetics of hearing loss, the next steps that are being taken in research, and their challenges. Furthermore, genotype-phenotype correlations and modifying factors are discussed as these are instrumental in counselling hearing impaired individuals and/or their family members.

Individual long non-coding RNAs have no overt functions in zebrafish embryogenesis, viability and fertility

Hundreds of long non-coding RNAs (lncRNAs) have been identified as potential regulators of gene expression, but their functions remain largely unknown. To study the role of lncRNAs during vertebrate development, we selected 25 zebrafish lncRNAs based on their conservation, expression profile or proximity to developmental regulators, and used CRISPR-Cas9 to generate 32 deletion alleles. We observed altered transcription of neighboring genes in some mutants, but none of the lncRNAs were required for embryogenesis, viability or fertility. Even RNAs with previously proposed non-coding functions (cyrano and squint) and other conserved lncRNAs (gas5 and lnc-setd1ba) were dispensable. In one case (lnc-phox2bb), absence of putative DNA regulatory-elements, but not of the lncRNA transcript itself, resulted in abnormal development. LncRNAs might have redundant, subtle, or context-dependent roles, but extrapolation from our results suggests that the majority of individual zebrafish lncRNAs have no overt roles in embryogenesis, viability and fertility.

Gut feelings: Studying enteric nervous system development, function, and disease in the zebrafish model system

The enteric nervous system (ENS) is the largest part of the peripheral nervous system and is entirely neural crest-derived. It provides the intrinsic innervation of the gut, controlling different aspects of gut function, such as motility. In this review, we will discuss key points of Zebrafish ENS development, genes, and signaling pathways regulating ENS development, as well as contributions of the Zebrafish model system to better understand ENS disorders. During their migration, enteric progenitor cells (EPCs) display a gradient of developmental states based on their proliferative and migratory characteristics, and show spatiotemporal heterogeneity based on gene expression patterns. Many genes and signaling pathways that regulate the migration and proliferation of EPCs have been identified, but later stages of ENS development, especially steps of neuronal and glial differentiation, remain poorly understood. In recent years, Zebrafish have become increasingly important to test candidate genes for ENS disorders (e.g., from genome-wide association studies), to identify environmental influences on ENS development (e.g., through large-scale drug screens), and to investigate the role the gut microbiota play in ENS development and disease. With its unique advantages as a model organism, Zebrafish will continue to contribute to a better understanding of ENS development, function, and disease. Developmental Dynamics 247:268-278, 2018. © 2017 Wiley Periodicals, Inc.

Genomic dissection of conserved transcriptional regulation in intestinal epithelial cells

The intestinal epithelium serves critical physiologic functions that are shared among all vertebrates. However, it is unknown how the transcriptional regulatory mechanisms underlying these functions have changed over the course of vertebrate evolution. We generated genome-wide mRNA and accessible chromatin data from adult intestinal epithelial cells (IECs) in zebrafish, stickleback, mouse, and human species to determine if conserved IEC functions are achieved through common transcriptional regulation. We found evidence for substantial common regulation and conservation of gene expression regionally along the length of the intestine from fish to mammals and identified a core set of genes comprising a vertebrate IEC signature. We also identified transcriptional start sites and other putative regulatory regions that are differentially accessible in IECs in all 4 species. Although these sites rarely showed sequence conservation from fish to mammals, surprisingly, they drove highly conserved IEC expression in a zebrafish reporter assay. Common putative transcription factor binding sites (TFBS) found at these sites in multiple species indicate that sequence conservation alone is insufficient to identify much of the functionally conserved IEC regulatory information. Among the rare, highly sequence-conserved, IEC-specific regulatory regions, we discovered an ancient enhancer upstream from her6/HES1 that is active in a distinct population of Notch-positive cells in the intestinal epithelium. Together, these results show how combining accessible chromatin and mRNA datasets with TFBS prediction and in vivo reporter assays can reveal tissue-specific regulatory information conserved across 420 million years of vertebrate evolution. We define an IEC transcriptional regulatory network that is shared between fish and mammals and establish an experimental platform for studying how evolutionarily distilled regulatory information commonly controls IEC development and physiology.

The epithelium lining the intestine is an ancient animal tissue that serves as a primary site of nutrient absorption and interaction with microbiota. Its formation and function require complex patterns of gene transcription that vary along the intestine and in specialized intestinal epithelial cell (IEC) subtypes. However, it is unknown how the underlying transcriptional regulatory mechanisms have changed over the course of vertebrate evolution. Here, we used genome-wide profiling of mRNA levels and chromatin accessibility to identify conserved IEC genes and regulatory regions in 4 vertebrate species (zebrafish, stickleback, mouse, and human) separated from a common ancestor by 420 million years. We identified substantial similarities in genes expressed along the vertebrate intestine. These data disclosed putative conserved transcription factor binding sites (TFBS) enriched in accessible chromatin near IEC genes and in regulatory sites with accessibility restricted to IECs. Fluorescent reporter assays in transparent zebrafish showed that these regions, which frequently lacked sequence conservation, were still capable of driving conserved expression patterns. We also found a highly conserved region near mammalian and fish hes1 sufficient to drive expression in a specific population of IECs with active Notch signaling. These results establish a platform to define the conserved transcriptional networks underlying vertebrate IEC physiology.

Transcriptional Dynamics at Brain Enhancers: from Functional Specialization to Neurodegeneration

Over the last decade, the noncoding part of the genome has been shown to harbour thousands of cis-regulatory elements, such as enhancers, that activate well-defined gene expression programs. Driven by the development of numerous techniques, many of these elements are now identified in multiple tissues and cell types, and their characteristics as well as importance in development and disease are becoming increasingly clear. Here, we provide an overview of the insights that were gained from the analysis of noncoding gene regulatory elements in the brain and describe their potential contribution to cell type specialization, brain function and neurodegenerative disease.

A sequence-based method to predict the impact of regulatory variants using random forest

Background

Most disease-associated variants identified by genome-wide association studies (GWAS) exist in noncoding regions. In spite of the common agreement that such variants may disrupt biological functions of their hosting regulatory elements, it remains a great challenge to characterize the risk of a genetic variant within the implicated genome sequence. Therefore, it is essential to develop an effective computational model that is not only capable of predicting the potential risk of a genetic variant but also valid in interpreting how the function of the genome is affected with the occurrence of the variant.

Results

We developed a method named kmerForest that used a random forest classifier with k-mer counts to predict accessible chromatin regions purely based on DNA sequences. We demonstrated that our method outperforms existing methods in distinguishing known accessible chromatin regions from random genomic sequences. Furthermore, the performance of our method can further be improved with the incorporation of sequence conservation features. Based on this model, we assessed importance of the k-mer features by a series of permutation experiments, and we characterized the risk of a single nucleotide polymorphism (SNP) on the function of the genome using the difference between the importance of the k-mer features affected by the occurrence of the SNP. We conducted a series of experiments and showed that our model can well discriminate between pathogenic and normal SNPs. Particularly, our model correctly prioritized SNPs that are proved to be enriched for the binding sites of FOXA1 in breast cancer cell lines from previous studies.

Conclusions

We presented a novel method to interpret functional genetic variants purely base on DNA sequences. The proposed k-mer based score offers an effective means of measuring the impact of SNPs on the function of the genome, and thus shedding light on the identification of genetic risk factors underlying complex traits and diseases.

Modeling intestinal disorders using zebrafish

Although the zebrafish was initially developed as a model system to study embryonic development, it has gained increasing attention as an advantageous system to investigate human diseases, including intestinal disorders. Zebrafish embryos develop rapidly, and their digestive system is fully functional and visible by 5days post fertilization. There is a large degree of homology between the intestine of zebrafish and higher vertebrate organisms in terms of its cellular composition and function as both a digestive and immune organ. Furthermore, molecular pathways regulating injury and immune responses are highly conserved. In this chapter, we provide an overview of studies addressing developmental and physiological processes relevant to human intestinal disease. These studies include those related to congenital disorders, host-microbiota interactions, inflammatory diseases, motility disorders, and intestinal cancer. We also highlight the utility of zebrafish to functionally validate candidate genes identified through mutational analyses and genome-wide association studies, and discuss methodologies to investigate the intestinal biology that are unique to zebrafish.

The Promoter and Multiple Enhancers of the pou4f3 Gene Regulate Expression in Inner Ear Hair Cells

Few enhancers that target gene expression to inner ear hair cells (HCs) have been identified. Using transgenic analysis of enhanced green fluorescent protein (eGFP) reporter constructs and bioinformatics, we evaluated the control of pou4f3 gene expression, since it is expressed only in HCs within the inner ear and continues to be expressed throughout life. An 8.5-kb genomic DNA fragment 5' to the start codon, containing three regions of high cross-species homology, drove expression in all embryonic and neonatal HCs, and adult vestibular and inner HCs, but not adult outer HCs. Transgenes with 0.4, 0.8, 2.5, or 6.5 kb of 5' DNA did not produce HC expression. However, addition of the region from 6.5 to 7.2 kb produced expression in vestibular HCs and neonatal basal turn outer HCs, which also implicated the region from 7.2 to 8.5 kb in inner and apical outer HC expression. Deletion of the region from 0.4 to 5.5 kb 5' from the 8.5-kb construct did not affect HC expression, further indicating lack of HC regulatory elements. When the region from 1 to 0.4 kb was replaced with the minimal promoter of the Ela1 gene, HC expression was maintained but at a drastically reduced level. Bioinformatics identified regions of highly conserved sequence outside of the 8.5 kb, which contained POU4F3-, GFI1-, and LHX3-binding sites. These regions may be involved in maintaining POU4F3 expression in adult outer HCs. Our results identify separate enhancers at various locations that direct expression to different HC types at different ages and determine that 0.4 kb of upstream sequence determines expression level. These data will assist in the identification of mutations in noncoding, regulatory regions of this deafness gene.

Enteric nervous system development in avian and zebrafish models

Our current understanding of the developmental biology of the enteric nervous system (ENS) and the genesis of ENS diseases is founded almost entirely on studies using model systems. Although genetic studies in the mouse have been at the forefront of this field over the last 20 years or so, historically it was the easy accessibility of the chick embryo for experimental manipulations that allowed the first descriptions of the neural crest origins of the ENS in the 1950s. More recently, studies in the chick and other non-mammalian model systems, notably zebrafish, have continued to advance our understanding of the basic biology of ENS development, with each animal model providing unique experimental advantages. Here we review the basic biology of ENS development in chick and zebrafish, highlighting conserved and unique features, and emphasising novel contributions to our general understanding of ENS development due to technical or biological features.

Functionally conserved enhancers with divergent sequences in distant vertebrates

BACKGROUND

To examine the contributions of sequence and function conservation in the evolution of enhancers, we systematically identified enhancers whose sequences are not conserved among distant groups of vertebrate species, but have homologous function and are likely to be derived from a common ancestral sequence. Our approach combined comparative genomics and epigenomics to identify potential enhancer sequences in the genomes of three groups of distantly related vertebrate species.

RESULTS

We searched for sequences that were conserved within groups of closely related species but not between groups of more distant species, and were associated with an epigenetic mark of enhancer activity. To facilitate inferring orthology between non-conserved sequences, we limited our search to introns whose orthology could be unambiguously established by mapping the bracketing exons. We show that a subset of these non-conserved but syntenic sequences from the mouse and zebrafish genomes have homologous functions in a zebrafish transgenic enhancer assay. The conserved expression patterns driven by these enhancers are probably associated with short transcription factor-binding motifs present in the divergent sequences.

CONCLUSIONS

We have identified numerous potential enhancers with divergent sequences but a conserved function. These results indicate that selection on function, rather than sequence, may be a common mode of enhancer evolution; evidence for selection at the sequence level is not a necessary criterion to define a gene regulatory element.

Variation in vertebrate cis-regulatory elements in evolution and disease

Much of the genetic information that drives animal diversity lies within the vast non-coding regions of the genome. Multi-species sequence conservation in non-coding regions of the genome flags important regulatory elements and more recently, techniques that look for functional signatures predicted for regulatory sequences have added to the identification of thousands more. For some time, biologists have argued that changes in cis-regulatory sequences creates the basic genetic framework for evolutionary change. Recent advances support this notion and show that there is extensive genomic variability in non-coding regulatory elements associated with trait variation, speciation and disease.

A single nucleotide polymorphism associated with isolated cleft lip and palate, thyroid cancer and hypothyroidism alters the activity of an oral epithelium and thyroid enhancer near FOXE1

Three common diseases, isolated cleft lip and cleft palate (CLP), hypothyroidism and thyroid cancer all map to the FOXE1 locus, but causative variants have yet to be identified. In patients with CLP, the frequency of coding mutations in FOXE1 fails to account for the risk attributable to this locus, suggesting that the common risk alleles reside in nearby regulatory elements. Using a combination of zebrafish and mouse transgenesis, we screened 15 conserved non-coding sequences for enhancer activity, identifying three that regulate expression in a tissue specific pattern consistent with endogenous foxe1 expression. These three, located -82.4, -67.7 and +22.6 kb from the FOXE1 start codon, are all active in the oral epithelium or branchial arches. The -67.7 and +22.6 kb elements are also active in the developing heart, and the -67.7 kb element uniquely directs expression in the developing thyroid. Within the -67.7 kb element is the SNP rs7850258 that is associated with all three diseases. Quantitative reporter assays in oral epithelial and thyroid cell lines show that the rs7850258 allele (G) associated with CLP and hypothyroidism has significantly greater enhancer activity than the allele associated with thyroid cancer (A). Moreover, consistent with predicted transcription factor binding differences, the -67.7 kb element containing rs7850258 allele G is significantly more responsive to both MYC and ARNT than allele A. By demonstrating that this common non-coding variant alters FOXE1 expression, we have identified at least in part the functional basis for the genetic risk of these seemingly disparate disorders.

A Phox2b BAC Transgenic Rat Line Useful for Understanding Respiratory Rhythm Generator Neural Circuitry

The key role of the respiratory neural center is respiratory rhythm generation to maintain homeostasis through the control of arterial blood pCO₂/pH and pO₂ levels. The neuronal network responsible for respiratory rhythm generation in neonatal rat resides in the ventral side of the medulla and is composed of two groups; the parafacial respiratory group (pFRG) and the pre-Bötzinger complex group (preBötC). The pFRG partially overlaps in the retrotrapezoid nucleus (RTN), which was originally identified in adult cats and rats. Part of the pre-inspiratory (Pre-I) neurons in the RTN/pFRG serves as central chemoreceptor neurons and the CO₂ sensitive Pre-I neurons express homeobox gene Phox2b. Phox2b encodes a transcription factor and is essential for the development of the sensory-motor visceral circuits. Mutations in human PHOX2B cause congenital hypoventilation syndrome, which is characterized by blunted ventilatory response to hypercapnia. Here we describe the generation of a novel transgenic (Tg) rat harboring fluorescently labeled Pre-I neurons in the RTN/pFRG. In addition, the Tg rat showed fluorescent signals in autonomic enteric neurons and carotid bodies. Because the Tg rat expresses inducible Cre recombinase in PHOX2B-positive cells during development, it is a potentially powerful tool for dissecting the entire picture of the respiratory neural network during development and for identifying the CO₂/O₂ sensor molecules in the adult central and peripheral nervous systems.

A large genomic deletion leads to enhancer adoption by the lamin B1 gene: a second path to autosomal dominant adult-onset demyelinating leukodystrophy (ADLD)

Chromosomal rearrangements with duplication of the lamin B1 (LMNB1) gene underlie autosomal dominant adult-onset demyelinating leukodystrophy (ADLD), a rare neurological disorder in which overexpression of LMNB1 causes progressive central nervous system demyelination. However, we previously reported an ADLD family (ADLD-1-TO) without evidence of duplication or other mutation in LMNB1 despite linkage to the LMNB1 locus and lamin B1 overexpression. By custom array-CGH, we further investigated this family and report here that patients carry a large (∼660 kb) heterozygous deletion that begins 66 kb upstream of the LMNB1 promoter. Lamin B1 overexpression was confirmed in further ADLD-1-TO tissues and in a postmortem brain sample, where lamin B1 was increased in the frontal lobe. Through parallel studies, we investigated both loss of genetic material and chromosomal rearrangement as possible causes of LMNB1 overexpression, and found that ADLD-1-TO plausibly results from an enhancer adoption mechanism. The deletion eliminates a genome topological domain boundary, allowing normally forbidden interactions between at least three forebrain-directed enhancers and the LMNB1 promoter, in line with the observed mainly cerebral localization of lamin B1 overexpression and myelin degeneration. This second route to LMNB1 overexpression and ADLD is a new example of the relevance of regulatory landscape modifications in determining Mendelian phenotypes.

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

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

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

Identification and functional studies of regulatory variants responsible for the association of NRG3 with a delusion phenotype in schizophrenia

We previously reported genetic linkage for Schizophrenia (SZ) (NPL of 4.7) at 10q22 in the Ashkenazi Jewish (AJ) population. In follow up fine mapping we found strong evidence of association between three intronic single nucleotide variants (SNVs) in the 5' end of Neuregulin 3 (NRG3) and the delusion factor score of our phenotypic principal component analysis. Two independent groups replicated these findings, indicating that variants in NRG3 confer risk for a delusion-rich SZ subtype. To identify the causative variants, we sequenced the 162 kb linkage disequilibrium (LD) block covering the NRG3 5' end in 47 AJ SZ patients at the extremes of the delusion factor quantitative trait distribution. Among the identified variants we found 5 noncoding SNVs present on the high delusion factor haplotype and significantly overrepresented in high delusion factor subjects. We tested these for regulatory effects and found that risk alleles of rs10883866 and rs60827755 decreased and increased, respectively, the expression of a reporter gene as compared to the reference allele. In post-mortem brain RNA quantification experiments we found the same variants also perturb relative expression of alternative NRG3 isoforms. In summary, we have identified regulatory SNVs contributing to the association of NRG3 with delusion symptoms in SZ.

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

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

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

In vivo visualization of the development of the enteric nervous system using a Tg(-8.3bphox2b:Kaede) transgenic zebrafish

The phox2b gene encodes a transcription factor that is expressed in the developing enteric nervous system (ENS). An enhancer element has been identified in the zebrafish phox2b locus that can drive tissue specific expression of reporter genes in enteric neuron precursor cells. We have generated a transgenic zebrafish line in which the Kaede fluorescent protein is under the control of this phox2b enhancer. This line has stable expression of the Kaede protein in enteric neuron precursor cells over three generations. To demonstrate the utility of this line we compared the migration and division rates of enteric neuron precursor cells in wild type and the zebrafish ENS mutant lessen.

Enhancer interaction networks as a means for singular olfactory receptor expression

The transcriptional activation of one out of ?2800 olfactory receptor (OR) alleles is a poorly understood process. Here, we identify a plethora of putative OR enhancers and study their in vivo activity in olfactory neurons. Distinguished by an unusual epigenetic signature, candidate OR enhancers are characterized by extensive interchromosomal interactions associated with OR transcription and share a similar pattern of transcription factor footprints. In particular, we establish the role of the transcription factor Bptf as a facilitator of both enhancer interactions and OR transcription. Our observations agree with the model whereby OR transcription occurs in the context of multiple interacting enhancers. Disruption of these interchromosomal interactions results in weak and multigenic OR expression, suggesting that the rare coincidence of numerous enhancers over a stochastically chosen OR may account for the singularity and robustness in OR transcription.

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.

Pervasive divergence of transcriptional gene regulation in Caenorhabditis nematodes

Because there is considerable variation in gene expression even between closely related species, it is clear that gene regulatory mechanisms evolve relatively rapidly. Because primary sequence conservation is an unreliable proxy for functional conservation of cis-regulatory elements, their assessment must be carried out in vivo. We conducted a survey of cis-regulatory conservation between C. elegans and closely related species C. briggsae, C. remanei, C. brenneri, and C. japonica. We tested enhancers of eight genes from these species by introducing them into C. elegans and analyzing the expression patterns they drove. Our results support several notable conclusions. Most exogenous cis elements direct expression in the same cells as their C. elegans orthologs, confirming gross conservation of regulatory mechanisms. However, the majority of exogenous elements, when placed in C. elegans, also directed expression in cells outside endogenous patterns, suggesting functional divergence. Recurrent ectopic expression of different promoters in the same C. elegans cells may reflect biases in the directions in which expression patterns can evolve due to shared regulatory logic of coexpressed genes. The fact that, despite differences between individual genes, several patterns repeatedly emerged from our survey, encourages us to think that general rules governing regulatory evolution may exist and be discoverable.

Integrating diverse datasets improves developmental enhancer prediction

Gene-regulatory enhancers have been identified using various approaches, including evolutionary conservation, regulatory protein binding, chromatin modifications, and DNA sequence motifs. To integrate these different approaches, we developed EnhancerFinder, a two-step method for distinguishing developmental enhancers from the genomic background and then predicting their tissue specificity. EnhancerFinder uses a multiple kernel learning approach to integrate DNA sequence motifs, evolutionary patterns, and diverse functional genomics datasets from a variety of cell types. In contrast with prediction approaches that define enhancers based on histone marks or p300 sites from a single cell line, we trained EnhancerFinder on hundreds of experimentally verified human developmental enhancers from the VISTA Enhancer Browser. We comprehensively evaluated EnhancerFinder using cross validation and found that our integrative method improves the identification of enhancers over approaches that consider a single type of data, such as sequence motifs, evolutionary conservation, or the binding of enhancer-associated proteins. We find that VISTA enhancers active in embryonic heart are easier to identify than enhancers active in several other embryonic tissues, likely due to their uniquely high GC content. We applied EnhancerFinder to the entire human genome and predicted 84,301 developmental enhancers and their tissue specificity. These predictions provide specific functional annotations for large amounts of human non-coding DNA, and are significantly enriched near genes with annotated roles in their predicted tissues and lead SNPs from genome-wide association studies. We demonstrate the utility of EnhancerFinder predictions through in vivo validation of novel embryonic gene regulatory enhancers from three developmental transcription factor loci. Our genome-wide developmental enhancer predictions are freely available as a UCSC Genome Browser track, which we hope will enable researchers to further investigate questions in developmental biology.

The human genome contains an immense amount of non-protein-coding DNA with unknown function. Some of this DNA regulates when, where, and at what levels genes are active during development. Enhancers, one type of regulatory element, are short stretches of DNA that can act as “switches” to turn a gene on or off at specific times in specific cells or tissues. Understanding where in the genome enhancers are located can provide insight into the genetic basis of development and disease. Enhancers are hard to identify, but clues about their locations are found in different types of data including DNA sequence, evolutionary history, and where proteins bind to DNA. Here, we introduce a new tool, called EnhancerFinder, which combines these data to predict the location and activity of enhancers during embryonic development. We trained EnhancerFinder on a large set of functionally validated human enhancers, and it proved to be very accurate. We used EnhancerFinder to predict tens of thousands of enhancers in the human genome and validated several of the predictions near three important developmental genes in mouse or zebrafish. EnhancerFinder's predictions will be useful in understanding functional regions hidden in the vast amounts of human non-coding DNA.

No gene in the genome makes sense except in the light of evolution

Evolutionary conservation has been an accurate predictor of functional elements across the first decade of metazoan genomics. More recently, there has been a move to define functional elements instead from biochemical annotations. Evolutionary methods are, however, more comprehensive than biochemical approaches can be and can assess quantitatively, especially for subtle effects, how biologically important--how injurious after mutation--different types of elements are. Evolutionary methods are thus critical for understanding the large fraction (up to 10%) of the human genome that does not encode proteins and yet might convey function. These methods can also capture the ephemeral nature of much noncoding functional sequence, with large numbers of functional elements having been gained and lost rapidly along each mammalian lineage. Here, we review how different strengths of purifying selection have impacted on protein-coding and non-protein-coding loci and on transcription factor binding sites in mammalian and fruit fly genomes.

Defining functional DNA elements in the human genome

With the completion of the human genome sequence, attention turned to identifying and annotating its functional DNA elements. As a complement to genetic and comparative genomics approaches, the Encyclopedia of DNA Elements Project was launched to contribute maps of RNA transcripts, transcriptional regulator binding sites, and chromatin states in many cell types. The resulting genome-wide data reveal sites of biochemical activity with high positional resolution and cell type specificity that facilitate studies of gene regulation and interpretation of noncoding variants associated with human disease. However, the biochemically active regions cover a much larger fraction of the genome than do evolutionarily conserved regions, raising the question of whether nonconserved but biochemically active regions are truly functional. Here, we review the strengths and limitations of biochemical, evolutionary, and genetic approaches for defining functional DNA segments, potential sources for the observed differences in estimated genomic coverage, and the biological implications of these discrepancies. We also analyze the relationship between signal intensity, genomic coverage, and evolutionary conservation. Our results reinforce the principle that each approach provides complementary information and that we need to use combinations of all three to elucidate genome function in human biology and disease.

Evolution of transcription factor binding in metazoans - mechanisms and functional implications

Differences in transcription factor binding can contribute to organismal evolution by altering downstream gene expression programmes. Genome-wide studies in Drosophila melanogaster and mammals have revealed common quantitative and combinatorial properties of in vivo DNA binding, as well as marked differences in the rate and mechanisms of evolution of transcription factor binding in metazoans. Here, we review the recently discovered rapid 're-wiring' of in vivo transcription factor binding between related metazoan species and summarize general principles underlying the observed patterns of evolution. We then consider what might explain the differences in genome evolution between metazoan phyla and outline the conceptual and technological challenges facing this research field.

A survey of ancient conserved non-coding elements in the PAX6 locus reveals a landscape of interdigitated cis-regulatory archipelagos

Biological differences between cell types and developmental processes are characterised by differences in gene expression profiles. Gene-distal enhancers are key components of the regulatory networks that specify the tissue-specific expression patterns driving embryonic development and cell fate decisions, and variations in their sequences are a major contributor to genetic disease and disease susceptibility. Despite advances in the methods for discovery of putative cis-regulatory sequences, characterisation of their spatio-temporal enhancer activities in a mammalian model system remains a major bottle-neck. We employed a strategy that combines gnathostome sequence conservation with transgenic mouse and zebrafish reporter assays to survey the genomic locus of the developmental control gene PAX6 for the presence of novel cis-regulatory elements. Sequence comparison between human and the cartilaginous elephant shark (Callorhinchus milii) revealed several ancient gnathostome conserved non-coding elements (agCNEs) dispersed widely throughout the PAX6 locus, extending the range of the known PAX6 cis-regulatory landscape to contain the full upstream PAX6-RCN1 intergenic region. Our data indicates that ancient conserved regulatory sequences can be tested effectively in transgenic zebrafish even when not conserved in zebrafish themselves. The strategy also allows efficient dissection of compound regulatory regions previously assessed in transgenic mice. Remarkable overlap in expression patterns driven by sets of agCNEs indicates that PAX6 resides in a landscape of multiple tissue-specific regulatory archipelagos.

Enhancer turnover and conserved regulatory function in vertebrate evolution

Mutations in regulatory regions including enhancers are an important source of variation and innovation during evolution. Enhancers can evolve by changes in the sequence, arrangement and repertoire of transcription factor binding sites, but whole enhancers can also be lost or gained in certain lineages in a process of turnover. The proopiomelanocortin gene (Pomc), which encodes a prohormone, is expressed in the pituitary and hypothalamus of all jawed vertebrates. We have previously described that hypothalamic Pomc expression in mammals is controlled by two enhancers-nPE1 and nPE2-that are derived from transposable elements and that presumably replaced the ancestral neuronal Pomc regulatory regions. Here, we show that nPE1 and nPE2, even though they are mammalian novelties with no homologous counterpart in other vertebrates, nevertheless can drive gene expression specifically to POMC neurons in the hypothalamus of larval and adult transgenic zebrafish. This indicates that when neuronal Pomc enhancers originated de novo during early mammalian evolution, the newly created cis- and trans-codes were similar to the ancestral ones. We also identify the neuronal regulatory region of zebrafish pomca and confirm that it is not homologous to the mammalian enhancers. Our work sheds light on the process of gene regulatory evolution by showing how a locus can undergo enhancer turnover and nevertheless maintain the ancestral transcriptional output.

Functional analysis of limb enhancers in the developing fin

Despite diverging ∼365 million years ago, tetrapod limbs and pectoral fins express similar genes that could be regulated by shared regulatory elements. In this study, we set out to analyze the ability of enhancers to maintain tissue specificity in these two divergent structures. We tested 22 human sequences that were previously reported as mouse limb enhancers for their enhancer activity in zebrafish (Danio rerio). Using a zebrafish enhancer assay, we found that 10/22 (45 %) were positive for pectoral fin activity. Analysis of the various criteria that correlated with positive fin activity found that both spatial limb activity and evolutionary conservation are not good predictors of fin enhancer activity. These results suggest that zebrafish enhancer assays may be limited in detecting human limb enhancers, and this limitation does not improve by the use of limb spatial expression or evolutionary conservation.

Regulation of MEIS1 by distal enhancer elements in acute leukemia

Aberrant activation of the three-amino-acid-loop extension (TALE) homeobox gene MEIS1 shortens the latency and accelerates the onset and progression of acute leukemia, yet the molecular mechanism underlying persistent activation of the MEIS1 gene in leukemia remains poorly understood. Here we used a combined comparative genomics analysis and an in vivo transgenic zebrafish assay to identify 6 regulatory DNA elements that are able to direct GFP expression in a spatiotemporal manner during zebrafish embryonic hematopoiesis. Analysis of chromatin characteristics and regulatory signatures suggest that many of these predicted elements are potential enhancers in mammalian hematopoiesis. Strikingly, one of the enhancer elements (E9) is a frequent integration site in retroviral induced mouse acute leukemia. The genomic region corresponding to enhancer E9 is differentially marked by H3K4 mono-methylation and H3K27 acetylation, hallmarks of active enhancers, in multiple leukemia cell lines. Decreased enrichment of these histone marks is associated with downregulation of MEIS1 expression during hematopoietic differentiation. Furthermore, MEIS1/HOXA9 transactivate this enhancer via a conserved binding motif in vitro, and participate in an autoregulatory loop that modulates MEIS1 expression in vivo. Our results suggest that an intronic enhancer regulates the expression of MEIS1 in hematopoiesis and contributes to its aberrant expression in acute leukemia.

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

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

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.

Computational identification of active enhancers in model organisms

As a class of cis-regulatory elements, enhancers were first identified as the genomic regions that are able to markedly increase the transcription of genes nearly 30years ago. Enhancers can regulate gene expression in a cell-type specific and developmental stage specific manner. Although experimental technologies have been developed to identify enhancers genome-wide, the design principle of the regulatory elements and the way they rewire the transcriptional regulatory network tempo-spatially are far from clear. At present, developing predictive methods for enhancers, particularly for the cell-type specific activity of enhancers, is central to computational biology. In this review, we survey the current computational approaches for active enhancer prediction and discuss future directions.

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

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

Genome-Level Analysis of Selective Constraint without Apparent Sequence Conservation

Conservation of function can be accompanied by obvious similarity of homologous sequences which may persist for billions of years (Iyer LM, Leipe DD, Koonin EV, Aravind L. 2004. Evolutionary history and higher order classification of AAA+ ATPases. J Struct Biol. 146:11–31.). However, presumably homologous segments of noncoding DNA can also retain their ancestral function even after their sequences diverge beyond recognition (Fisher S, Grice EA, Vinton RM, Bessling SL, McCallion AS. 2006. Conservation of RET regulatory function from human to zebrafish without sequence similarity. Science 312:276–279.). To investigate this phenomenon at the genomic scale, we studied homologous introns in a quartet of insect species, and in a quartet of vertebrate species. Each quartet consisted of two pairs of moderately distant genomes, with a much larger evolutionary distance between the pairs. In both quartets, we found that introns that carry a regulatory segment or a conserved segment in the first pair tend to carry a conserved segment in the second pair, even though no similarity of these segments could be detected between the two pairs. Furthermore, introns from one pair that are preserved in the other pair tend to carry a conserved segment within the first pair, and be longer in the first pair, compared with the introns that were lost between pairs, even though no similarity between pairs could be detected in such preserved introns. These results indicate that selective constraint, presumably caused by conservation of the ancestral function, often persists even after the homologous DNA segments become unalignable.

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

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

What does our genome encode?

In its first production phase, The ENCODE Project Consortium (ENCODE) has generated thousands of genome-scale data sets, resulting in a genomic “parts list” that encompasses transcripts, sites of transcription factor binding, and other functional features that now number in the millions of distinct elements. These data are reshaping many long-held beliefs concerning the information content of the human and other complex genomes, including the very definition of the gene. Here I discuss and place in context many of the leading findings of ENCODE, as well as trends that are shaping the generation and interpretation of ENCODE data. Finally, I consider prospects for the future, including maximizing the accuracy, completeness, and utility of ENCODE data for the community.

Function and regulation of AUTS2, a gene implicated in autism and human evolution

Nucleotide changes in the AUTS2 locus, some of which affect only noncoding regions, are associated with autism and other neurological disorders, including attention deficit hyperactivity disorder, epilepsy, dyslexia, motor delay, language delay, visual impairment, microcephaly, and alcohol consumption. In addition, AUTS2 contains the most significantly accelerated genomic region differentiating humans from Neanderthals, which is primarily composed of noncoding variants. However, the function and regulation of this gene remain largely unknown. To characterize auts2 function, we knocked it down in zebrafish, leading to a smaller head size, neuronal reduction, and decreased mobility. To characterize AUTS2 regulatory elements, we tested sequences for enhancer activity in zebrafish and mice. We identified 23 functional zebrafish enhancers, 10 of which were active in the brain. Our mouse enhancer assays characterized three mouse brain enhancers that overlap an ASD-associated deletion and four mouse enhancers that reside in regions implicated in human evolution, two of which are active in the brain. Combined, our results show that AUTS2 is important for neurodevelopment and expose candidate enhancer sequences in which nucleotide variation could lead to neurological disease and human-specific traits.

Conserved regulatory architecture underlies parallel genetic changes and convergent phenotypic evolution

Similar morphological, physiological, and behavioral features have evolved independently in different species, a pattern known as convergence. It is known that morphological convergence can occur through changes in orthologous genes. In some cases of convergence, cis-regulatory changes generate parallel modifications in the expression patterns of orthologous genes. Our understanding of how changes in cis-regulatory regions contribute to convergence is hampered, usually, by a limited understanding of the global cis-regulatory structure of the evolving genes. Here we examine the genetic causes of a case of precise phenotypic convergence between Drosophila sechellia and Drosophila ezoana, species that diverged ~40 Mya. Previous studies revealed that changes in multiple transcriptional enhancers of shavenbaby (svb, a transcript of the ovo locus) caused phenotypic evolution in the D. sechellia lineage. It has also been shown that the convergent phenotype of D. ezoana was likely caused by cis-regulatory evolution of svb. Here we show that the large-scale cis-regulatory architecture of svb is conserved between these Drosophila species. Furthermore, we show that the D. ezoana orthologs of the evolved D. sechellia enhancers have also evolved expression patterns that correlate precisely with the changes in the phenotype. Our results suggest that phenotypic convergence resulted from multiple noncoding changes that occurred in parallel in the D. sechellia and D. ezoana lineages.

Long distance relationships: enhancer-promoter communication and dynamic gene transcription

The three-dimensional regulation of gene transcription involves loop formation between enhancer and promoter elements, controlling spatiotemporal gene expression in multicellular organisms. Enhancers are usually located in non-coding DNA and can activate gene transcription by recruiting transcription factors, chromatin remodeling factors and RNA Polymerase II. Research over the last few years has revealed that enhancers have tell-tale characteristics that facilitate their detection by several approaches, although the hallmarks of enhancers are not always uniform. Enhancers likely play an important role in the activation of genes by functioning as a primary point of contact for transcriptional activators, and by making physical contact with gene promoters often by means of a chromatin loop. Although numerous transcriptional regulators participate in the formation of chromatin loops that bring enhancers into proximity with promoters, the mechanism(s) of enhancer-promoter connectivity remain enigmatic. Here we discuss enhancer function, review some of the many proteins shown to be involved in establishing enhancer-promoter loops, and describe the dynamics of enhancer-promoter contacts during development, differentiation and in specific cell types.

Identifying and characterizing regulatory sequences in the human genome with chromatin accessibility assays

After finishing a human genome reference sequence in 2002, the genomics community has turned to the task of interpreting it. A primary focus is to identify and characterize not only protein-coding genes, but all functional elements in the genome. The effort includes both individual investigators and large-scale projects like the Encyclopedia of DNA Elements (ENCODE) project. As part of the ENCODE project, several groups have identified millions of regulatory elements in hundreds of human cell-types using DNase-seq and FAIRE-seq experiments that detect regions of nucleosome-free open chromatin. ChIP-seq experiments have also been used to discover transcription factor binding sites and map histone modifications. Nearly all identified elements are found in non-coding DNA, hypothesizing a function for previously unannotated sequence. In this review, we provide an overview of the ENCODE effort to define regulatory elements, summarize the main results, and discuss implications of the millions of regulatory elements distributed throughout the genome.

Integration of ChIP-seq and machine learning reveals enhancers and a predictive regulatory sequence vocabulary in melanocytes

We take a comprehensive approach to the study of regulatory control of gene expression in melanocytes that proceeds from large-scale enhancer discovery facilitated by ChIP-seq; to rigorous validation in silico, in vitro, and in vivo; and finally to the use of machine learning to elucidate a regulatory vocabulary with genome-wide predictive power. We identify 2489 putative melanocyte enhancer loci in the mouse genome by ChIP-seq for EP300 and H3K4me1. We demonstrate that these putative enhancers are evolutionarily constrained, enriched for sequence motifs predicted to bind key melanocyte transcription factors, located near genes relevant to melanocyte biology, and capable of driving reporter gene expression in melanocytes in culture (86%; 43/50) and in transgenic zebrafish (70%; 7/10). Next, using the sequences of these putative enhancers as a training set for a supervised machine learning algorithm, we develop a vocabulary of 6-mers predictive of melanocyte enhancer function. Lastly, we demonstrate that this vocabulary has genome-wide predictive power in both the mouse and human genomes. This study provides deep insight into the regulation of gene expression in melanocytes and demonstrates a powerful approach to the investigation of regulatory sequences that can be applied to other cell types.

The chromatin fingerprint of gene enhancer elements

Different cell types within a single organism are generally distinguished by strikingly different patterns of gene expression, which are dynamic throughout development and adult life. Distal enhancer elements are key drivers of spatiotemporal specificity in gene regulation. Often located tens of kilobases from their target promoters and functioning in an orientation-independent manner, the identification of bona fide enhancers has proved a formidable challenge. With the development of ChIP-seq, global cataloging of putative enhancers has become feasible. Here, we review the current understanding of the chromatin landscape at enhancers and how these chromatin features enable robust identification of tissue-specific enhancers.

Interferon regulatory factor 6 promotes differentiation of the periderm by activating expression of Grainyhead-like 3

Interferon Regulatory Factor 6 (IRF6) is a transcription factor that, in mammals, is required for the differentiation of skin, breast epithelium, and oral epithelium. However, the transcriptional targets that mediate these effects are currently unknown. In zebrafish and frog embryos Irf6 is necessary for differentiation of the embryonic superficial epithelium, or periderm. Here we use microarrays to identify genes that are expressed in the zebrafish periderm and whose expression is inhibited by a dominant-negative variant of Irf6 (dnIrf6). These methods identify Grhl3, an ancient regulator of the epidermal permeability barrier, as acting downstream of Irf6. In human keratinocytes, IRF6 binds conserved elements near the GHRL3 promoter. We show that one of these elements has enhancer activity in human keratinocytes and zebrafish periderm, suggesting that Irf6 directly stimulates Grhl3 expression in these tissues. Simultaneous inhibition of grhl1 and grhl3 disrupts periderm differentiation in zebrafish, and, intriguingly, forced grhl3 expression restores periderm markers in both zebrafish injected with dnIrf6 and frog embryos depleted of Irf6. Finally, in Irf6 deficient mouse embryos, Grhl3 expression in the periderm and oral epithelium is virtually absent. These results indicate that Grhl3 is a key effector of Irf6 in periderm differentiation.

Functional characterization of tissue-specific enhancers in the DLX5/6 locus

Disruption of distaless homeobox 5 and 6 (Dlx5/6) in mice results in brain, craniofacial, genital, ear and limb defects. In humans, chromosomal aberrations in the DLX5/6 region, some of which do not encompass DLX5/6, are associated with split hand/foot malformation 1 (SHFM1) as well as intellectual disability, craniofacial anomalies and hearing loss, suggesting that the disruption of DLX5/6 regulatory elements could lead to these abnormalities. Here, we characterized enhancers in the DLX5/6 locus whose tissue-specific expression and genomic location along with previously characterized enhancers correlate with phenotypes observed in individuals with chromosomal abnormalities. By analyzing chromosomal aberrations at 7q21, we refined the minimal SHFM1 critical region and used comparative genomics to select 26 evolutionary conserved non-coding sequences in this critical region for zebrafish enhancer assays. Eight of these sequences were shown to function as brain, olfactory bulb, branchial arch, otic vesicle and fin enhancers, recapitulating dlx5a/6a expression. Using a mouse enhancer assay, several of these zebrafish enhancers showed comparable expression patterns in the branchial arch, otic vesicle, forebrain and/or limb at embryonic day 11.5. Examination of the coordinates of various chromosomal rearrangements in conjunction with the genomic location of these tissue-specific enhancers showed a correlation with the observed clinical abnormalities. Our findings suggest that chromosomal abnormalities that disrupt the function of these tissue-specific enhancers could be the cause of SHFM1 and its associated phenotypes. In addition, they highlight specific enhancers in which mutations could lead to non-syndromic hearing loss, craniofacial defects or limb malformations.

Systematic elucidation and in vivo validation of sequences enriched in hindbrain transcriptional control

Illuminating the primary sequence encryption of enhancers is central to understanding the regulatory architecture of genomes. We have developed a machine learning approach to decipher motif patterns of hindbrain enhancers and identify 40,000 sequences in the human genome that we predict display regulatory control that includes the hindbrain. Consistent with their roles in hindbrain patterning, MEIS1, NKX6-1, as well as HOX and POU family binding motifs contributed strongly to this enhancer model. Predicted hindbrain enhancers are overrepresented at genes expressed in hindbrain and associated with nervous system development, and primarily reside in the areas of open chromatin. In addition, 77 (0.2%) of these predictions are identified as hindbrain enhancers on the VISTA Enhancer Browser, and 26,000 (60%) overlap enhancer marks (H3K4me1 or H3K27ac). To validate these putative hindbrain enhancers, we selected 55 elements distributed throughout our predictions and six low scoring controls for evaluation in a zebrafish transgenic assay. When assayed in mosaic transgenic embryos, 51/55 elements directed expression in the central nervous system. Furthermore, 30/34 (88%) predicted enhancers analyzed in stable zebrafish transgenic lines directed expression in the larval zebrafish hindbrain. Subsequent analysis of sequence fragments selected based upon motif clustering further confirmed the critical role of the motifs contributing to the classifier. Our results demonstrate the existence of a primary sequence code characteristic to hindbrain enhancers. This code can be accurately extracted using machine-learning approaches and applied successfully for de novo identification of hindbrain enhancers. This study represents a critical step toward the dissection of regulatory control in specific neuronal subtypes.