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

CRISPR-Cas9 editing of non-coding genomic loci as a means of controlling gene expression in the sea urchin

We seek to manipulate gene function here through CRISPR-Cas9 editing of cis-regulatory sequences, rather than the more typical mutation of coding regions. This approach would minimize secondary effects of cellular responses to nonsense mediated decay pathways or to mutant protein products by premature stops. This strategy also allows for reducing gene activity in cases where a complete gene knockout would result in lethality, and it can be applied to the rapid identification of key regulatory sites essential for gene expression. We tested this strategy here with genes of known function as a proof of concept, and then applied it to examine the upstream genomic region of the germline gene Nanos2 in the sea urchin, Strongylocentrotus purpuratus. We first used CRISPR-Cas9 to target established genomic cis-regulatory regions of the skeletogenic cell transcription factor, Alx1, and the TGF-β signaling ligand, Nodal, which produce obvious developmental defects when altered in sea urchin embryos. Importantly, mutation of cis-activator sites (Alx1) and cis-repressor sites (Nodal) result in the predicted decreased and increased transcriptional output, respectively. Upon identification of efficient gRNAs by genomic mutations, we then used the same validated gRNAs to target a deadCas9-VP64 transcriptional activator to increase Nodal transcription directly. Finally, we paired these new methodologies with a more traditional, GFP reporter construct approach to further our understanding of the transcriptional regulation of Nanos2, a key gene required for germ cell identity in S. purpuratus. With a series of reporter assays, upstream Cas9-promoter targeted mutagenesis, coupled with qPCR and in situ RNA hybridization, we concluded that the promoter of Nanos2 drives strong mRNA expression in the sea urchin embryo, indicating that its primordial germ cell (PGC)-specific restriction may rely instead on post-transcriptional regulation. Overall, we present a proof-of-principle tool-kit of Cas9-mediated manipulations of promoter regions that should be applicable in most cells and embryos for which CRISPR-Cas9 is employed.

ATAC-Seq for Assaying Chromatin Accessibility Protocol Using Echinoderm Embryos

Cis-regulatory elements (CREs) and transcription factors (TFs) associated with them determine temporal and spatial domains of gene expression. Therefore, identification of these CREs and TFs is crucial to elucidating transcriptional programs across taxa. With chromatin accessibility facilitating transcription factor access to DNA, the identification of regions of open chromatin sheds light both on the function of the regulatory elements and their evolution, thus allowing the recognition of potential CREs. Buenrostro and colleagues have developed a novel method for exploring chromatin accessibility: assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), which can be used for the purpose of identifying putative CREs. This method was shown to have considerable advantages when compared to traditional methods such as sequence conservation analyses or functional assays. Here we present the adaptation of the ATAC-seq method to echinoderm species and discuss how it can be used for CRE discovery.

Evaluation of Whole-Exome Enrichment Solutions: Lessons from the High-End of the Short-Read Sequencing Scale

Whole-exome sequencing has become a popular technique in research and clinical settings, assisting in disease diagnosis and increasing the understanding of disease pathogenesis. In this study, we aimed to compare common enrichment capture solutions available in the market. Peripheral blood-purified DNA samples were enriched with SureSelect^QXT V6 (Agilent) and various Illumina solutions: TruSeq DNA Nano, TruSeq DNA Exome, Nextera DNA Exome, and Illumina DNA Prep with Enrichment, and sequenced on a HiSeq 4000. We found that their percentage of duplicate reads was as much as 2 times higher than previously reported values for the previous HiSeq series. SureSelect^QXT and Illumina DNA Prep with Enrichment showed the best average on-target coverage, which improved when off-target regions were included. At high coverage levels and in shared bases, these two solutions and TruSeq DNA Exome provided three of the best performances. With respect to the number of small variants detected, SureSelect^QXT presented the lowest number of detected variants in target regions. When off-target regions were considered, its ability equalized to other solutions. Our results show SureSelect^QXT and Illumina DNA Prep with Enrichment to be the best enrichment capture solutions.

Functional interrogation of non-coding DNA through CRISPR genome editing

Methodologies to interrogate non-coding regions have lagged behind coding regions despite comprising the vast majority of the genome. However, the rapid evolution of clustered regularly interspaced short palindromic repeats (CRISPR)-based genome editing has provided a multitude of novel techniques for laboratory investigation including significant contributions to the toolbox for studying non-coding DNA. CRISPR-mediated loss-of-function strategies rely on direct disruption of the underlying sequence or repression of transcription without modifying the targeted DNA sequence. CRISPR-mediated gain-of-function approaches similarly benefit from methods to alter the targeted sequence through integration of customized sequence into the genome as well as methods to activate transcription. Here we review CRISPR-based loss- and gain-of-function techniques for the interrogation of non-coding DNA.

Enhancer Turnover Is Associated with a Divergent Transcriptional Response to Glucocorticoid in Mouse and Human Macrophages

Phenotypic differences between individuals and species are controlled in part through differences in expression of a relatively conserved set of genes. Genes expressed in the immune system are subject to especially powerful selection. We have investigated the evolution of both gene expression and candidate enhancers in human and mouse macrophages exposed to glucocorticoid (GC), a regulator of innate immunity and an important therapeutic agent. Our analyses revealed a very limited overlap in the repertoire of genes responsive to GC in human and mouse macrophages. Peaks of inducible binding of the GC receptor (GR) detected by chromatin immunoprecipitation-Seq correlated with induction, but not repression, of target genes in both species, occurred at distal regulatory sites not promoters, and were strongly enriched for the consensus GR-binding motif. Turnover of GR binding between mice and humans was associated with gain and loss of the motif. There was no detectable signal of positive selection at species-specific GR binding sites, but clear evidence of purifying selection at the small number of conserved sites. We conclude that enhancer divergence underlies the difference in transcriptional activation after GC treatment between mouse and human macrophages. Only the shared inducible loci show evidence of selection, and therefore these loci may be important for the subset of responses to GC that is shared between species.

Analysis of 5' gene regions reveals extraordinary conservation of novel non-coding sequences in a wide range of animals

Background

Phylogenetic footprinting is a comparative method based on the principle that functional sequence elements will acquire fewer mutations over time than non-functional sequences. Successful comparisons of distantly related species will thus yield highly important sequence elements likely to serve fundamental biological roles. RNA regulatory elements are less well understood than those in DNA. In this study we use the emerging model organism Nasonia vitripennis, a parasitic wasp, in a comparative analysis against 12 insect genomes to identify deeply conserved non-coding elements (CNEs) conserved in large groups of insects, with a focus on 5’ UTRs and promoter sequences.

Results

We report the identification of 322 CNEs conserved across a broad range of insect orders. The identified regions are associated with regulatory and developmental genes, and contain short footprints revealing aspects of their likely function in translational regulation. The most ancient regions identified in our analysis were all found to overlap transcribed regions of genes, reflecting stronger conservation of translational regulatory elements than transcriptional elements. Further expanding sequence analyses to non-insect species we also report the discovery of, to our knowledge, the two oldest and most ubiquitous CNE’s yet described in the animal kingdom (700 MYA). These ancient conserved non-coding elements are associated with the two ribosomal stalk genes, RPLP1 and RPLP2, and were very likely functional in some of the earliest animals.

Conclusions

We report the identification of the most deeply conserved CNE’s found to date, and several other deeply conserved elements which are without exception, part of 5’ untranslated regions of transcripts, and occur in a number of key translational regulatory genes, highlighting translational regulation of translational regulators as a conserved feature of insect genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0499-6) contains supplementary material, which is available to authorized users.

Evidence for deep regulatory similarities in early developmental programs across highly diverged insects

Many genes familiar from Drosophila development, such as the so-called gap, pair-rule, and segment polarity genes, play important roles in the development of other insects and in many cases appear to be deployed in a similar fashion, despite the fact that Drosophila-like "long germband" development is highly derived and confined to a subset of insect families. Whether or not these similarities extend to the regulatory level is unknown. Identification of regulatory regions beyond the well-studied Drosophila has been challenging as even within the Diptera (flies, including mosquitoes) regulatory sequences have diverged past the point of recognition by standard alignment methods. Here, we demonstrate that methods we previously developed for computational cis-regulatory module (CRM) discovery in Drosophila can be used effectively in highly diverged (250-350 Myr) insect species including Anopheles gambiae, Tribolium castaneum, Apis mellifera, and Nasonia vitripennis. In Drosophila, we have successfully used small sets of known CRMs as "training data" to guide the search for other CRMs with related function. We show here that although species-specific CRM training data do not exist, training sets from Drosophila can facilitate CRM discovery in diverged insects. We validate in vivo over a dozen new CRMs, roughly doubling the number of known CRMs in the four non-Drosophila species. Given the growing wealth of Drosophila CRM annotation, these results suggest that extensive regulatory sequence annotation will be possible in newly sequenced insects without recourse to costly and labor-intensive genome-scale experiments. We develop a new method, Regulus, which computes a probabilistic score of similarity based on binding site composition (despite the absence of nucleotide-level sequence alignment), and demonstrate similarity between functionally related CRMs from orthologous loci. Our work represents an important step toward being able to trace the evolutionary history of gene regulatory networks and defining the mechanisms underlying insect evolution.

Identifying transcriptional cis-regulatory modules in animal genomes

Gene expression is regulated through the activity of transcription factors (TFs) and chromatin-modifying proteins acting on specific DNA sequences, referred to as cis-regulatory elements. These include promoters, located at the transcription initiation sites of genes, and a variety of distal cis-regulatory modules (CRMs), the most common of which are transcriptional enhancers. Because regulated gene expression is fundamental to cell differentiation and acquisition of new cell fates, identifying, characterizing, and understanding the mechanisms of action of CRMs is critical for understanding development. CRM discovery has historically been challenging, as CRMs can be located far from the genes they regulate, have few readily identifiable sequence characteristics, and for many years were not amenable to high-throughput discovery methods. However, the recent availability of complete genome sequences and the development of next-generation sequencing methods have led to an explosion of both computational and empirical methods for CRM discovery in model and nonmodel organisms alike. Experimentally, CRMs can be identified through chromatin immunoprecipitation directed against TFs or histone post-translational modifications, identification of nucleosome-depleted 'open' chromatin regions, or sequencing-based high-throughput functional screening. Computational methods include comparative genomics, clustering of known or predicted TF-binding sites, and supervised machine-learning approaches trained on known CRMs. All of these methods have proven effective for CRM discovery, but each has its own considerations and limitations, and each is subject to a greater or lesser number of false-positive identifications. Experimental confirmation of predictions is essential, although shortcomings in current methods suggest that additional means of validation need to be developed. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

Using giant scarlet runner bean embryos to uncover regulatory networks controlling suspensor gene activity

One of the major unsolved issues in plant development is understanding the regulatory networks that control the differential gene activity that is required for the specification and development of the two major embryonic regions, the embryo proper and suspensor. Historically, the giant embryo of scarlet runner bean (SRB), Phaseolus coccineus, has been used as a model system to investigate the physiological events that occur early in embryogenesis-focusing on the question of what role the suspensor region plays. A major feature distinguishing SRB embryos from those of other plants is a highly enlarged suspensor containing at least 200 cells that synthesize growth regulators required for subsequent embryonic development. Recent studies have exploited the giant size of the SRB embryo to micro-dissect the embryo proper and suspensor regions in order to use genomics-based approaches to identify regulatory genes that may be involved in controlling suspensor and embryo proper differentiation, as well as the cellular processes that may be unique to each embryonic region. Here we review the current genomics resources that make SRB embryos a compelling model system for studying the early events required to program embryo development.

Shared transcription factors contribute to distinct cell fates

Genome-wide transcription factor (TF) binding profiles differ dramatically between cell types. However, not much is known about the relationship between cell-type-specific binding patterns and gene expression. A recent study demonstrated how the same TFs can have functional roles when binding to largely non-overlapping genomic regions in hematopoietic progenitor and mast cells. Cell-type specific binding profiles of shared TFs are therefore not merely the consequence of opportunistic and functionally irrelevant binding to accessible chromatin, but instead have the potential to make meaningful contributions to cell-type specific transcriptional programs.

Chromatin accessibility mapping identifies mediators of basal transcription and retinoid-induced repression of OTX2 in medulloblastoma

Despite an emerging understanding of the genetic alterations giving rise to various tumors, the mechanisms whereby most oncogenes are overexpressed remain unclear. Here we have utilized an integrated approach of genomewide regulatory element mapping via DNase-seq followed by conventional reporter assays and transcription factor binding site discovery to characterize the transcriptional regulation of the medulloblastoma oncogene Orthodenticle Homeobox 2 (OTX2). Through these studies we have revealed that OTX2 is differentially regulated in medulloblastoma at the level of chromatin accessibility, which is in part mediated by DNA methylation. In cell lines exhibiting chromatin accessibility of OTX2 regulatory regions, we found that autoregulation maintains OTX2 expression. Comparison of medulloblastoma regulatory elements with those of the developing brain reveals that these tumors engage a developmental regulatory program to drive OTX2 transcription. Finally, we have identified a transcriptional regulatory element mediating retinoid-induced OTX2 repression in these tumors. This work characterizes for the first time the mechanisms of OTX2 overexpression in medulloblastoma. Furthermore, this study establishes proof of principle for applying ENCODE datasets towards the characterization of upstream trans-acting factors mediating expression of individual genes.

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.

Targeted CNS Delivery Using Human MiniPromoters and Demonstrated Compatibility with Adeno-Associated Viral Vectors

Critical for human gene therapy is the availability of small promoter tools to drive gene expression in a highly specific and reproducible manner. We tackled this challenge by developing human DNA MiniPromoters using computational biology and phylogenetic conservation. MiniPromoters were tested in mouse as single-copy knock-ins at the Hprt locus on the X Chromosome, and evaluated for lacZ reporter expression in CNS and non-CNS tissue. Eighteen novel MiniPromoters driving expression in mouse brain were identified, two MiniPromoters for driving pan-neuronal expression, and 17 MiniPromoters for the mouse eye. Key areas of therapeutic interest were represented in this set: the cerebral cortex, embryonic hypothalamus, spinal cord, bipolar and ganglion cells of the retina, and skeletal muscle. We also demonstrated that three retinal ganglion cell MiniPromoters exhibit similar cell-type specificity when delivered via adeno-associated virus (AAV) vectors intravitreally. We conclude that our methodology and characterization has resulted in desirable expression characteristics that are intrinsic to the MiniPromoter, not dictated by copy number effects or genomic location, and results in constructs predisposed to success in AAV. These MiniPromoters are immediately applicable for pre-clinical studies towards gene therapy in humans, and are publicly available to facilitate basic and clinical research, and human gene therapy.

Evolution of transcriptional enhancers and animal diversity

Deciphering the genetic bases that drive animal diversity is one of the major challenges of modern biology. Although four decades ago it was proposed that animal evolution was mainly driven by changes in cis-regulatory DNA elements controlling gene expression rather than in protein-coding sequences, only now are powerful bioinformatics and experimental approaches available to accelerate studies into how the evolution of transcriptional enhancers contributes to novel forms and functions. In the introduction to this Theme Issue, we start by defining the general properties of transcriptional enhancers, such as modularity and the coexistence of tight sequence conservation with transcription factor-binding site shuffling as different mechanisms that maintain the enhancer grammar over evolutionary time. We discuss past and current methods used to identify cell-type-specific enhancers and provide examples of how enhancers originate de novo, change and are lost in particular lineages. We then focus in the central part of this Theme Issue on analysing examples of how the molecular evolution of enhancers may change form and function. Throughout this introduction, we present the main findings of the articles, reviews and perspectives contributed to this Theme Issue that together illustrate some of the great advances and current frontiers in the field.

Study of Cis-regulatory Elements in the Ascidian Ciona intestinalis

The ascidian (sea squirt) C. intestinalis has become an important model organism for the study of cis-regulation. This is largely due to the technology that has been developed for assessing cis-regulatory activity through the use of transient reporter transgenes introduced into fertilized eggs. This technique allows the rapid and inexpensive testing of endogenous or altered DNA for regulatory activity in vivo. This review examines evidence that C. intestinaliscis-regulatory elements are located more closely to coding regions than in other model organisms. I go on to compare the organization of cis-regulatory elements and conserved non-coding sequences in Ciona, mammals, and other deuterostomes for three representative C.intestinalis genes, Pax6, FoxAa, and the DlxA-B cluster, along with homologs in the other species. These comparisons point out some of the similarities and differences between cis-regulatory elements and their study in the various model organisms. Finally, I provide illustrations of how C. intestinalis lends itself to detailed study of the structure of cis-regulatory elements, which have led, and promise to continue to lead, to important insights into the fundamentals of transcriptional regulation.

Finding of IFNγ gene enhancers and their core sequences

DNA segmentation methods were used to study which fragments of the human IFNγ gene possess enhancer activity. The human IFNγ gene was divided into 240-bp fragments, which were inserted between the GFP gene and the Alu tandem sequence to determine whether the inserted sequences eliminate the inhibition induced by the Alu tandem sequence. We found that five different 240-bp fragments (FUIFN3F3R, IFN4F4R, IFN6F6R, IFN21F21R, and IFN22F22R) and two 60-bp core sequences (IFN6-2F2R and IFN21-3-4F3-4R) derived from the IFNγ gene contain enhancers that can activate the GFP reporter gene. These enhancers may be targets of IFNγ gene expression regulation.

Intronic cis-regulatory modules mediate tissue-specific and microbial control of angptl4/fiaf transcription

The intestinal microbiota enhances dietary energy harvest leading to increased fat storage in adipose tissues. This effect is caused in part by the microbial suppression of intestinal epithelial expression of a circulating inhibitor of lipoprotein lipase called Angiopoietin-like 4 (Angptl4/Fiaf). To define the cis-regulatory mechanisms underlying intestine-specific and microbial control of Angptl4 transcription, we utilized the zebrafish system in which host regulatory DNA can be rapidly analyzed in a live, transparent, and gnotobiotic vertebrate. We found that zebrafish angptl4 is transcribed in multiple tissues including the liver, pancreatic islet, and intestinal epithelium, which is similar to its mammalian homologs. Zebrafish angptl4 is also specifically suppressed in the intestinal epithelium upon colonization with a microbiota. In vivo transgenic reporter assays identified discrete tissue-specific regulatory modules within angptl4 intron 3 sufficient to drive expression in the liver, pancreatic islet β-cells, or intestinal enterocytes. Comparative sequence analyses and heterologous functional assays of angptl4 intron 3 sequences from 12 teleost fish species revealed differential evolution of the islet and intestinal regulatory modules. High-resolution functional mapping and site-directed mutagenesis defined the minimal set of regulatory sequences required for intestinal activity. Strikingly, the microbiota suppressed the transcriptional activity of the intestine-specific regulatory module similar to the endogenous angptl4 gene. These results suggest that the microbiota might regulate host intestinal Angptl4 protein expression and peripheral fat storage by suppressing the activity of an intestine-specific transcriptional enhancer. This study provides a useful paradigm for understanding how microbial signals interact with tissue-specific regulatory networks to control the activity and evolution of host gene transcription.

Transcriptional enhancers in ascidian development

The study of cis-regulatory DNAs that control developmental gene expression is integral to the modeling of comprehensive genomic regulatory networks for embryogenesis. Ascidian embryos provide a unique opportunity for the analysis of cis-regulatory DNAs with cellular resolution in the context of a simple but typical chordate body plan. Here, we review landmark studies that have laid the foundations for the study of transcriptional enhancers, among other cis-regulatory DNAs, and their roles in ascidian development. The studies using ascidians of the Ciona genus have capitalized on a unique electroporation technique that permits the simultaneous transfection of hundreds of fertilized eggs, which develop rapidly and express transgenes with little mosaicism. Current studies using the ascidian embryo benefit from extensively annotated genomic resources to characterize transcript models in silico. The search for functional noncoding sequences can be guided by bioinformatic analyses combining evolutionary conservation, gene coexpression, and combinations of overrepresented short-sequence motifs. The power of the transient transfection assays has allowed thorough dissection of numerous cis-regulatory modules, which provided insights into the functional constraints that shape enhancer architecture and diversification. Future studies will benefit from pioneering stable transgenic lines and the analysis of chromatin states. Whole genome expression, functional and DNA binding data are being integrated into comprehensive genomic regulatory network models of early ascidian cell specification with a single-cell resolution that is unique among chordate model systems.

Computational strategies for the genome-wide identification of cis-regulatory elements and transcriptional targets

Transcription factors (TFs) are key proteins that decode the information in our genome to express a precise and unique set of proteins and RNA molecules in each cell type in our body. These factors play a pivotal role in all biological processes, including the determination of a cell's fate during development and the maintenance of a cell's physiological function. To achieve this, a TF binds to specific DNA sequences in the noncoding part of the genome, recruits chromatin modifiers and cofactors, and directs the transcription initiation rate of its "target genes." Therefore, a key challenge in deciphering a transcriptional switch is to identify the direct target genes of the master regulators that control the switch, the cis-regulatory elements implementing (auto-)regulatory loops, and the target genes of all the TFs in the downstream regulatory network. A better knowledge of a TF's targetome during specification and differentiation of a particular cell type will generate mechanistic insight into its developmental program. Here, I review computational strategies and methods to predict transcriptional targets by genome-wide searches for TF binding sites using position weight matrices, motif clusters, phylogenetic footprinting, chromatin binding and accessibility data, enhancer classification, motif enrichment, and gene expression signatures.

Erroneous attribution of relevant transcription factor binding sites despite successful prediction of cis-regulatory modules

Background

Cis-regulatory modules are bound by transcription factors to regulate gene expression. Characterizing these DNA sequences is central to understanding gene regulatory networks and gaining insight into mechanisms of transcriptional regulation, but genome-scale regulatory module discovery remains a challenge. One popular approach is to scan the genome for clusters of transcription factor binding sites, especially those conserved in related species. When such approaches are successful, it is typically assumed that the activity of the modules is mediated by the identified binding sites and their cognate transcription factors. However, the validity of this assumption is often not assessed.

Results

We successfully predicted five new cis-regulatory modules by combining binding site identification with sequence conservation and compared these to unsuccessful predictions from a related approach not utilizing sequence conservation. Despite greatly improved predictive success, the positive set had similar degrees of sequence and binding site conservation as the negative set. We explored the reasons for this by mutagenizing putative binding sites in three cis-regulatory modules. A large proportion of the tested sites had little or no demonstrable role in mediating regulatory element activity. Examination of loss-of-function mutants also showed that some transcription factors supposedly binding to the modules are not required for their function.

Conclusions

Our results raise important questions about interpreting regulatory module predictions obtained by finding clusters of conserved binding sites. Attribution of function to these sites and their cognate transcription factors may be incorrect even when modules are successfully identified. Our study underscores the importance of empirical validation of computational results even when these results are in line with expectation.

High-resolution mapping of open chromatin in the rice genome

Gene expression is controlled by the complex interaction of transcription factors binding to promoters and other regulatory DNA elements. One common characteristic of the genomic regions associated with regulatory proteins is a pronounced sensitivity to DNase I digestion. We generated genome-wide high-resolution maps of DNase I hypersensitive (DH) sites from both seedling and callus tissues of rice (Oryza sativa). Approximately 25% of the DH sites from both tissues were found in putative promoters, indicating that the vast majority of the gene regulatory elements in rice are not located in promoter regions. We found 58% more DH sites in the callus than in the seedling. For DH sites detected in both the seedling and callus, 31% displayed significantly different levels of DNase I sensitivity within the two tissues. Genes that are differentially expressed in the seedling and callus were frequently associated with DH sites in both tissues. The DNA sequences contained within the DH sites were hypomethylated, consistent with what is known about active gene regulatory elements. Interestingly, tissue-specific DH sites located in the promoters showed a higher level of DNA methylation than the average DNA methylation level of all the DH sites located in the promoters. A distinct elevation of H3K27me3 was associated with intergenic DH sites. These results suggest that epigenetic modifications play a role in the dynamic changes of the numbers and DNase I sensitivity of DH sites during development.

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

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