Insulators and promoters: closer than we think

Beyond modular enhancers: new questions in cis-regulatory evolution

Our understanding of how cis-regulatory elements work has advanced rapidly, outpacing our evolutionary models. In this review, we consider the implications of new mechanistic findings for evolutionary developmental biology. We focus on three different debates: whether evolutionary innovation occurs more often via the modification of old cis-regulatory elements or the emergence of new ones; the extent to which individual elements are specific and autonomous or multifunctional and interdependent; and how the robustness of cis-regulatory architectures influences the rate of trait evolution. These discussions lead us to propose new questions for the evo-devo of cis-regulation.

Transcriptional silencing in Saccharomyces cerevisiae: known unknowns

Transcriptional silencing in Saccharomyces cerevisiae is a persistent and highly stable form of gene repression. It involves DNA silencers and repressor proteins that bind nucleosomes. The silenced state is influenced by numerous factors including the concentration of repressors, nature of activators, architecture of regulatory elements, modifying enzymes and the dynamics of chromatin.Silencers function to increase the residence time of repressor Sir proteins at silenced domains while clustering of silenced domains enables increased concentrations of repressors and helps facilitate long-range interactions. The presence of an accessible NDR at the regulatory regions of silenced genes, the cycling of chromatin configurations at regulatory sites, the mobility of Sir proteins, and the non-uniform distribution of the Sir proteins across the silenced domain, all result in silenced chromatin that only stably silences weak promoters and enhancers via changes in transcription burst duration and frequency.These data collectively suggest that silencing is probabilistic and the robustness of silencing is achieved through sub-optimization of many different nodes of action such that a stable expression state is generated and maintained even though individual constituents are in constant flux.

Yeast heterochromatin stably silences only weak regulatory elements by altering burst duration

Transcriptional silencing in Saccharomyces cerevisiae involves the generation of a chromatin state that stably represses transcription. Using multiple reporter assays, a diverse set of upstream activating sequence enhancers and core promoters were investigated for their susceptibility to silencing. We show that heterochromatin stably silences only weak and stress-induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements, and the partial repression of these elements did not result in bistable expression states. Permutation analysis of enhancers and promoters indicates that both elements are targets of repression. Chromatin remodelers help specific regulatory elements to resist repression, most probably by altering nucleosome mobility and changing transcription burst duration. The strong enhancers/promoters can be repressed if silencer-bound Sir1 is increased. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating-type gene regulatory elements but not strong housekeeping gene regulatory sequences.

Yeast Heterochromatin Only Stably Silences Weak Regulatory Elements by Altering Burst Duration

The interplay between nucleosomes and transcription factors leads to programs of gene expression. Transcriptional silencing involves the generation of a chromatin state that represses transcription and is faithfully propagated through DNA replication and cell division. Using multiple reporter assays, including directly visualizing transcription in single cells, we investigated a diverse set of UAS enhancers and core promoters for their susceptibility to heterochromatic gene silencing. These results show that heterochromatin only stably silences weak and stress induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements and the partial repression did not result in bistable expression states. Permutation analysis of different UAS enhancers and core promoters indicate that both elements function together to determine the susceptibility of regulatory sequences to repression. Specific histone modifiers and chromatin remodellers function in an enhancer specific manner to aid these elements to resist repression suggesting that Sir proteins likely function in part by reducing nucleosome mobility. We also show that the strong housekeeping regulatory elements can be repressed if silencer bound Sir1 is increased, suggesting that Sir1 is a limiting component in silencing. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating type gene regulatory elements but not strong housekeeping gene regulatory sequences which could help explain why these genes are often found at the boundaries of silenced domains.

Cell factory engineering: Challenges and opportunities for synthetic biology applications

The production of high-quality recombinant proteins is critical to maintaining a continuous supply of biopharmaceuticals, such as therapeutic antibodies. Engineering mammalian cell factories presents a number of limitations typically associated with the proteotoxic stress induced upon aberrant accumulation of off-pathway protein folding intermediates, which eventually culminate in the induction of apoptosis. In this review, we will discuss advances in cell engineering and their applications at different hierarchical levels of control of the expression of recombinant proteins, from transcription and translational to posttranslational modifications and subcellular trafficking. We also highlight challenges and unique opportunities to apply modern synthetic biology tools to the design of programmable cell factories for improved biomanufacturing of therapeutic proteins.

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

The sound of silence: Transgene silencing in mammalian cell engineering

To elucidate principles operating in native biological systems and to develop novel biotechnologies, synthetic biology aims to build and integrate synthetic gene circuits within native transcriptional networks. The utility of synthetic gene circuits for cell engineering relies on the ability to control the expression of all constituent transgene components. Transgene silencing, defined as the loss of expression over time, persists as an obstacle for engineering primary cells and stem cells with transgenic cargos. In this review, we highlight the challenge that transgene silencing poses to the robust engineering of mammalian cells, outline potential molecular mechanisms of silencing, and present approaches for preventing transgene silencing. We conclude with a perspective identifying future research directions for improving the performance of synthetic gene circuits.

CTCF DNA binding domain undergoes dynamic and selective protein–protein interactions

CTCF is a predominant insulator protein required for three-dimensional chromatin organization. However, the roles of its insulation of enhancers in a 3D nuclear organization have not been fully explained. Here, we found that the CTCF DNA-binding domain (DBD) forms dynamic self-interacting clusters. Strikingly, CTCF DBD clusters were found to incorporate other insulator proteins but are not coenriched with transcriptional activators in the nucleus. This property is not observed in other domains of CTCF or the DBDs of other transcription factors. Moreover, endogenous CTCF shows a phenotype consistent with the DBD by forming small protein clusters and interacting with CTCF motif arrays that have fewer transcriptional activators bound. Our results reveal an interesting phenomenon in which CTCF DBD interacts with insulator proteins and selectively localizes to nuclear positions with lower concentrations of transcriptional activators, providing insights into the insulation function of CTCF.

•

The CTCF DNA-binding domain forms protein clusters in vivo and in vitro

•

CTCF DBD clusters colocalize with insulator proteins but not with activators

•

Arginine residues of CTCF DBD are frequently mutated in cancers

•

Multiple transcription factor DBDs form protein clusters

Dynamic spreading of chromatin-mediated gene silencing and reactivation between neighboring genes in single cells

In mammalian cells genes that are in close proximity can be transcriptionally coupled: silencing or activating one gene can affect its neighbors. Understanding these dynamics is important for natural processes, such as heterochromatin spreading during development and aging, and when designing synthetic gene regulation circuits. Here, we systematically dissect this process in single cells by recruiting and releasing repressive chromatin regulators at dual-gene synthetic reporters, and measuring how fast gene silencing and reactivation spread as a function of intergenic distance and configuration of insulator elements. We find that silencing by KRAB, associated with histone methylation, spreads between two genes within hours, with a time delay that increases with distance. This fast KRAB-mediated spreading is not blocked by the classical cHS4 insulators. Silencing by histone deacetylase HDAC4 of the upstream gene can also facilitate background silencing of the downstream gene by PRC2, but with a days-long delay that does not change with distance. This slower silencing can sometimes be stopped by insulators. Gene reactivation of neighboring genes is also coupled, with strong promoters and insulators determining the order of reactivation. Our data can be described by a model of multi-gene regulation that builds upon previous knowledge of heterochromatin spreading, where both gene silencing and gene reactivation can act at a distance, allowing for coordinated dynamics via chromatin regulator recruitment.

Enhancer RNA: What we know and what we can achieve

Enhancers are important cis-acting elements that can regulate gene transcription and cell fate alongside promoters. In fact, many human cancers and diseases are associated with the malfunction of enhancers. Recent studies have shown that enhancers can produce enhancer RNAs (eRNAs) by RNA polymerase II. In this review, we discuss eRNA production, characteristics, functions and mechanics. eRNAs can determine chromatin accessibility, histone modification and gene expression by constructing a 'chromatin loop', thereby bringing enhancers to their target gene. eRNA can also be involved in the phase separation with enhancers and other proteins. eRNAs are abundant, and importantly, tissue-specific in tumours, various diseases and stem cells; thus, eRNAs can be a potential target for disease diagnosis and treatment. As eRNA is produced from the active transcription of enhancers and is involved in the regulation of cell fate, its manipulation will influence cell function, and therefore, it can be a new target for biological therapy.

Insulators in Plants: Progress and Open Questions

The genomes of higher eukaryotes are partitioned into topologically associated domains or TADs, and insulators (also known as boundary elements) are the key elements responsible for their formation and maintenance. Insulators were first identified and extensively studied in Drosophila as well as mammalian genomes, and have also been described in yeast and plants. In addition, many insulator proteins are known in Drosophila, and some have been investigated in mammals. However, much less is known about this important class of non-coding DNA elements in plant genomes. In this review, we take a detailed look at known plant insulators across different species and provide an overview of potential determinants of plant insulator functions, including cis-elements and boundary proteins. We also discuss methods previously used in attempts to identify plant insulators, provide a perspective on their importance for research and biotechnology, and discuss areas of potential future research.

B2 SINE Copies Serve as a Transposable Boundary of DNA Methylation and Histone Modifications in the Mouse

More than one million copies of short interspersed elements (SINEs), a class of retrotransposons, are present in the mammalian genomes, particularly within gene-rich genomic regions. Evidence has accumulated that ancient SINE sequences have acquired new binding sites for transcription factors (TFs) through multiple mutations following retrotransposition, and as a result have rewired the host regulatory network during the course of evolution. However, it remains unclear whether currently active SINEs contribute to the expansion of TF binding sites. To study the mobility, expression, and function of SINE copies, we first identified about 2,000 insertional polymorphisms of SINE B1 and B2 families within Mus musculus. Using a novel RNA sequencing method designated as melRNA-seq, we detected the expression of SINEs in male germ cells at both the subfamily and genomic copy levels: the vast majority of B1 RNAs originated from evolutionarily young subfamilies, whereas B2 RNAs originated from both young and old subfamilies. DNA methylation and chromatin immunoprecipitation-sequencing (ChIP-seq) analyses in liver revealed that polymorphic B2 insertions served as a boundary element inhibiting the expansion of DNA hypomethylated and histone hyperacetylated regions, and decreased the expression of neighboring genes. Moreover, genomic B2 copies were enriched at the boundary of various histone modifications, and chromatin insulator protein, CCCTC-binding factor, a well-known chromatin boundary protein, bound to >100 polymorphic and >10,000 non-polymorphic B2 insertions. These results suggest that the currently active B2 copies are mobile boundary elements that can modulate chromatin modifications and gene expression, and are likely involved in epigenomic and phenotypic diversification of the mouse species.

Integrating distal and proximal information to predict gene expression via a densely connected convolutional neural network

MOTIVATION

Interactions among cis-regulatory elements such as enhancers and promoters are main driving forces shaping context-specific chromatin structure and gene expression. Although there have been computational methods for predicting gene expression from genomic and epigenomic information, most of them neglect long-range enhancer-promoter interactions, due to the difficulty in precisely linking regulatory enhancers to target genes. Recently, HiChIP, a novel high-throughput experimental approach, has generated comprehensive data on high-resolution interactions between promoters and distal enhancers. Moreover, plenty of studies suggest that deep learning achieves state-of-the-art performance in epigenomic signal prediction, and thus promoting the understanding of regulatory elements. In consideration of these two factors, we integrate proximal promoter sequences and HiChIP distal enhancer-promoter interactions to accurately predict gene expression.

RESULTS

We propose DeepExpression, a densely connected convolutional neural network, to predict gene expression using both promoter sequences and enhancer-promoter interactions. We demonstrate that our model consistently outperforms baseline methods, not only in the classification of binary gene expression status but also in regression of continuous gene expression levels, in both cross-validation experiments and cross-cell line predictions. We show that the sequential promoter information is more informative than the experimental enhancer information; meanwhile, the enhancer-promoter interactions within ±100 kbp around the TSS of a gene are most beneficial. We finally visualize motifs in both promoter and enhancer regions and show the match of identified sequence signatures with known motifs. We expect to see a wide spectrum of applications using HiChIP data in deciphering the mechanism of gene regulation.

AVAILABILITY AND IMPLEMENTATION

DeepExpression is freely available at https://github.com/wanwenzeng/DeepExpression.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Overlapping but Distinct Sequences Play Roles in the Insulator and Promoter Activities of the Drosophila BEAF-Dependent scs' Insulator

Chromatin domain insulators are thought to help partition the genome into genetic units called topologically associating domains (TADs). In Drosophila, TADs are often separated by inter-TAD regions containing active housekeeping genes and associated insulator binding proteins. This raises the question of whether insulator binding proteins are involved primarily in chromosomal TAD architecture or gene activation, or if these two activities are linked. The Boundary Element-Associated Factor of 32 kDa (BEAF-32, or BEAF for short) is usually found in inter-TADs. BEAF was discovered based on binding to the scs' insulator, and is important for the insulator activity of scs' and other BEAF binding sites. There are divergent promoters in scs' with a BEAF binding site by each. Here, we dissect the scs' insulator to identify DNA sequences important for insulator and promoter activity, focusing on the half of scs' with a high affinity BEAF binding site. We find that the BEAF binding site is important for both insulator and promoter activity, as is another sequence we refer to as LS4. Aside from that, different sequences play roles in insulator and promoter activity. So while there is overlap and BEAF is important for both, insulator and promoter activity can be separated.

Towards a comprehensive catalogue of validated and target-linked human enhancers

The human gene catalogue is essentially complete, but we lack an equivalently vetted inventory of bona fide human enhancers. Hundreds of thousands of candidate enhancers have been nominated via biochemical annotations; however, only a handful of these have been validated and confidently linked to their target genes. Here we review emerging technologies for discovering, characterizing and validating human enhancers at scale. We furthermore propose a new framework for operationally defining enhancers that accommodates the heterogeneous and complementary results that are emerging from reporter assays, biochemical measurements and CRISPR screens.

Characterization of the cHS4 insulator in mouse embryonic stem cells

Synthetic biology circuits are often constructed with multiple gene expression units assembled in close proximity, and they can be used to perform complex functions in embryonic stem cells (ESCs). However, mutual interference between transcriptional units has not been well studied in mouse ESCs. To assess the efficiency of insulators at suppressing promoter interference in mouse ESCs, we used an evaluation scheme in which a tunable tetracycline response element promoter is connected to a constant Nanog promoter. The chicken hypersensitive site 4 (cHS4) insulator, widely used both for enhancer blocking and for barrier insulation in vitro and in vivo, was positioned between the two expression units for assessment. By inserting the cassette into various loci of the mouse ESC genome with PiggyBac transposon, we were able to quantitatively examine the protective effect of cHS4 by gradually increasing the transcriptional activity of the tetracycline response element promoter with doxycycline and then measuring the transcriptional activity of the Nanog promoter. Our results indicate that the cHS4 insulator has minimal insulating effects on promoter interference in mouse ESCs. Further studies show that the cHS4 insulation effect may be promoter specific and related to interaction with CCCTC‐binding factor‐mediated loop formation. In addition, we also compared DNA transposition and transgene expression with or without the cHS4 insulator using well‐established ESC reporters. The results indicate that cHS4 has no apparent effects on DNA transposition and transgene expression levels, but exerts modest protective effects on long‐term transgene silencing.

Silencers, Enhancers, and the Multifunctional Regulatory Genome

Negative regulation of gene expression by transcriptional silencers has been difficult to study due to limited defined examples. A new study by Gisselbrecht et al. has dramatically increased the number of identified silencers and reveals that they are bifunctional regulatory sequences that also act as gene expression-promoting enhancers.

A CRISPR-Cas9 screen identifies essential CTCF anchor sites for estrogen receptor-driven breast cancer cell proliferation

Estrogen receptor α (ERα) is an enhancer activating transcription factor, a key driver of breast cancer and a main target for cancer therapy. ERα-mediated gene regulation requires proper chromatin-conformation to facilitate interactions between ERα-bound enhancers and their target promoters. A major determinant of chromatin structure is the CCCTC-binding factor (CTCF), that dimerizes and together with cohesin stabilizes chromatin loops and forms the boundaries of topologically associated domains. However, whether CTCF-binding elements (CBEs) are essential for ERα-driven cell proliferation is unknown. To address this question in a global manner, we implemented a CRISPR-based functional genetic screen targeting CBEs located in the vicinity of ERα-bound enhancers. We identified four functional CBEs and demonstrated the role of one of them in inducing chromatin conformation changes in favor of activation of PREX1, a key ERα target gene in breast cancer. Indeed, high PREX1 expression is a bona-fide marker of ERα-dependency in cell lines, and is associated with good outcome after anti-hormonal treatment. Altogether, our data show that distinct CTCF-mediated chromatin structures are required for ERα- driven breast cancer cell proliferation.

AIKYATAN: mapping distal regulatory elements using convolutional learning on GPU

Background

The data deluge can leverage sophisticated ML techniques for functionally annotating the regulatory non-coding genome. The challenge lies in selecting the appropriate classifier for the specific functional annotation problem, within the bounds of the hardware constraints and the model’s complexity. In our system Aikyatan, we annotate distal epigenomic regulatory sites, e.g., enhancers. Specifically, we develop a binary classifier that classifies genome sequences as distal regulatory regions or not, given their histone modifications’ combinatorial signatures. This problem is challenging because the regulatory regions are distal to the genes, with diverse signatures across classes (e.g., enhancers and insulators) and even within each class (e.g., different enhancer sub-classes).

Results

We develop a suite of ML models, under the banner Aikyatan, including SVM models, random forest variants, and deep learning architectures, for distal regulatory element (DRE) detection. We demonstrate, with strong empirical evidence, deep learning approaches have a computational advantage. Plus, convolutional neural networks (CNN) provide the best-in-class accuracy, superior to the vanilla variant. With the human embryonic cell line H1, CNN achieves an accuracy of 97.9% and an order of magnitude lower runtime than the kernel SVM. Running on a GPU, the training time is sped up 21x and 30x (over CPU) for DNN and CNN, respectively. Finally, our CNN model enjoys superior prediction performance vis-‘a-vis the competition. Specifically, Aikyatan-CNN achieved 40% higher validation rate versus CSIANN and the same accuracy as RFECS.

Conclusions

Our exhaustive experiments using an array of ML tools validate the need for a model that is not only expressive but can scale with increasing data volumes and diversity. In addition, a subset of these datasets have image-like properties and benefit from spatial pooling of features. Our Aikyatan suite leverages diverse epigenomic datasets that can then be modeled using CNNs with optimized activation and pooling functions. The goal is to capture the salient features of the integrated epigenomic datasets for deciphering the distal (non-coding) regulatory elements, which have been found to be associated with functional variants. Our source code will be made publicly available at: https://bitbucket.org/cellsandmachines/aikyatan.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-3049-1) contains supplementary material, which is available to authorized users.

The Synthetic Biology Toolkit for Photosynthetic Microorganisms

Photosynthetic microorganisms offer novel characteristics as synthetic biology chassis, and the toolbox of components and techniques for cyanobacteria and algae is rapidly increasing.

Gene Therapy Tools for Brain Diseases

Neurological disorders affecting the central nervous system (CNS) are still incompletely understood. Many of these disorders lack a cure and are seeking more specific and effective treatments. In fact, in spite of advancements in knowledge of the CNS function, the treatment of neurological disorders with modern medical and surgical approaches remains difficult for many reasons, such as the complexity of the CNS, the limited regenerative capacity of the tissue, and the difficulty in conveying conventional drugs to the organ due to the blood-brain barrier. Gene therapy, allowing the delivery of genetic materials that encodes potential therapeutic molecules, represents an attractive option. Gene therapy can result in a stable or inducible expression of transgene(s), and can allow a nearly specific expression in target cells. In this review, we will discuss the most commonly used tools for the delivery of genetic material in the CNS, including viral and non-viral vectors; their main applications; their advantages and disadvantages. We will discuss mechanisms of genetic regulation through cell-specific and inducible promoters, which allow to express gene products only in specific cells and to control their transcriptional activation. In addition, we will describe the applications to CNS diseases of post-transcriptional regulation systems (RNA interference); of systems allowing spatial or temporal control of expression [optogenetics and Designer Receptors Exclusively Activated by Designer Drugs (DREADDs)]; and of gene editing technologies (CRISPR/Cas9, Zinc finger proteins). Particular attention will be reserved to viral vectors derived from herpes simplex type 1, a potential tool for the delivery and expression of multiple transgene cassettes simultaneously.

tRNA Genes Affect Chromosome Structure and Function via Local Effects

The genome is packaged and organized in an ordered, nonrandom manner, and specific chromatin segments contact nuclear substructures to mediate this organization. tRNA genes (tDNAs) are binding sites for transcription factors and architectural proteins and are thought to play an important role in the organization of the genome. In this study, we investigate the roles of tDNAs in genomic organization and chromosome function by editing a chromosome so that it lacked any tDNAs. Surprisingly our analyses of this tDNA-less chromosome show that loss of tDNAs does not grossly affect chromatin architecture or chromosome tethering and mobility. However, loss of tDNAs affects local nucleosome positioning and the binding of SMC proteins at these loci. The absence of tDNAs also leads to changes in centromere clustering and a reduction in the frequency of long-range HML-HMR heterochromatin clustering with concomitant effects on gene silencing. We propose that the tDNAs primarily affect local chromatin structure, which results in effects on long-range chromosome architecture.

Sequence Characteristics Distinguish Transcribed Enhancers from Promoters and Predict Their Breadth of Activity

Enhancers and promoters both regulate gene expression by recruiting transcription factors (TFs); however, the degree to which enhancer vs. promoter activity is due to differences in their sequences or to genomic context is the subject of ongoing debate. We examined this question by analyzing the sequences of thousands of transcribed enhancers and promoters from hundreds of cellular contexts previously identified by cap analysis of gene expression. Support vector machine classifiers trained on counts of all possible 6-bp-long sequences (6-mers) were able to accurately distinguish promoters from enhancers and distinguish their breadth of activity across tissues. Classifiers trained to predict enhancer activity also performed well when applied to promoter prediction tasks, but promoter-trained classifiers performed poorly on enhancers. This suggests that the learned sequence patterns predictive of enhancer activity generalize to promoters, but not vice versa. Our classifiers also indicate that there are functionally relevant differences in enhancer and promoter GC content beyond the influence of CpG islands. Furthermore, sequences characteristic of broad promoter or broad enhancer activity matched different TFs, with predicted ETS- and RFX-binding sites indicative of promoters, and AP-1 sites indicative of enhancers. Finally, we evaluated the ability of our models to distinguish enhancers and promoters defined by histone modifications. Separating these classes was substantially more difficult, and this difference may contribute to ongoing debates about the similarity of enhancers and promoters. In summary, our results suggest that high-confidence transcribed enhancers and promoters can largely be distinguished based on biologically relevant sequence properties.

Identifying individual risk rare variants using protein structure guided local tests (POINT)

Rare variants are of increasing interest to genetic association studies because of their etiological contributions to human complex diseases. Due to the rarity of the mutant events, rare variants are routinely analyzed on an aggregate level. While aggregation analyses improve the detection of global-level signal, they are not able to pinpoint causal variants within a variant set. To perform inference on a localized level, additional information, e.g., biological annotation, is often needed to boost the information content of a rare variant. Following the observation that important variants are likely to cluster together on functional domains, we propose a protein structure guided local test (POINT) to provide variant-specific association information using structure-guided aggregation of signal. Constructed under a kernel machine framework, POINT performs local association testing by borrowing information from neighboring variants in the 3-dimensional protein space in a data-adaptive fashion. Besides merely providing a list of promising variants, POINT assigns each variant a p-value to permit variant ranking and prioritization. We assess the selection performance of POINT using simulations and illustrate how it can be used to prioritize individual rare variants in PCSK9, ANGPTL4 and CETP in the Action to Control Cardiovascular Risk in Diabetes (ACCORD) clinical trial data.

Genome organization and DNA accessibility control antigenic variation in trypanosomes

Many evolutionarily distant pathogenic organisms have evolved similar survival strategies to evade the immune responses of their hosts. These include antigenic variation, through which an infecting organism prevents clearance by periodically altering the identity of proteins that are visible to the immune system of the host1. Antigenic variation requires large reservoirs of immunologically diverse antigen genes, which are often generated through homologous recombination, as well as mechanisms to ensure the expression of one or very few antigens at any given time. Both homologous recombination and gene expression are affected by three-dimensional genome architecture and local DNA accessibility2,3. Factors that link three-dimensional genome architecture, local chromatin conformation and antigenic variation have, to our knowledge, not yet been identified in any organism. One of the major obstacles to studying the role of genome architecture in antigenic variation has been the highly repetitive nature and heterozygosity of antigen-gene arrays, which has precluded complete genome assembly in many pathogens. Here we report the de novo haplotype-specific assembly and scaffolding of the long antigen-gene arrays of the model protozoan parasite Trypanosoma brucei, using long-read sequencing technology and conserved features of chromosome folding4. Genome-wide chromosome conformation capture (Hi-C) reveals a distinct partitioning of the genome, with antigen-encoding subtelomeric regions that are folded into distinct, highly compact compartments. In addition, we performed a range of analyses-Hi-C, fluorescence in situ hybridization, assays for transposase-accessible chromatin using sequencing and single-cell RNA sequencing-that showed that deletion of the histone variants H3.V and H4.V increases antigen-gene clustering, DNA accessibility across sites of antigen expression and switching of the expressed antigen isoform, via homologous recombination. Our analyses identify histone variants as a molecular link between global genome architecture, local chromatin conformation and antigenic variation.

Spermiogenesis and Male Fertility Require the Function of Suppressor of Hairy-Wing in Somatic Cyst Cells of Drosophila

Drosophila Suppressor of Hairy-wing [Su(Hw)] protein is an example of a multivalent transcription factor. Although best known for its role in establishing the chromatin insulator of the gypsy retrotransposon, Su(Hw) functions as an activator and repressor at non-gypsy genomic sites. It remains unclear how the different regulatory activities of Su(Hw) are utilized during development. Motivated from observations of spatially restricted expression of Su(Hw) in the testis, we investigated the role of Su(Hw) in spermatogenesis to advance an understanding of its developmental contributions as an insulator, repressor, and activator protein. We discovered that Su(Hw) is required for sustained male fertility. Although dynamics of Su(Hw) expression coincide with changes in nuclear architecture and activation of coregulated testis-specific gene clusters, we show that loss of Su(Hw) does not disrupt meiotic chromosome pairing or transcription of testis-specific genes, suggesting that Su(Hw) has minor architectural or insulator functions in the testis. Instead, Su(Hw) has a prominent role as a repressor of neuronal genes, consistent with suggestions that Su(Hw) is a functional homolog of mammalian REST, a repressor of neuronal genes in non-neuronal tissues. We show that Su(Hw) regulates transcription in both germline and somatic cells. Surprisingly, the essential spermatogenesis function of Su(Hw) resides in somatic cyst cells, implying context-specific consequences due to loss of this transcription factor. Together, our studies highlight that Su(Hw) has a major developmental function as a transcriptional repressor, with the effect of its loss dependent upon the cell-specific factors.

A Zoom-Focus algorithm (ZFA) to locate the optimal testing region for rare variant association tests

Motivation

Increasing amounts of whole exome or genome sequencing data present the challenge of analysing rare variants with extremely small minor allele frequencies. Various statistical tests have been proposed, which are specifically configured to increase power for rare variants by conducting the test within a certain bin, such as a gene or a pathway. However, a gene may contain from several to thousands of markers, and not all of them are related to the phenotype. Combining functional and non-functional variants in an arbitrary genomic region could impair the testing power.

Results

We propose a Zoom-Focus algorithm (ZFA) to locate the optimal testing region within a given genomic region. It can be applied as a wrapper function in existing rare variant association tests to increase testing power. The algorithm consists of two steps. In the first step, Zooming, a given genomic region is partitioned by an order of two, and the best partition is located. In the second step, Focusing, the boundaries of the zoomed region are refined. Simulation studies showed that ZFA substantially increased the statistical power of rare variants' tests, including the SKAT, SKAT-O, burden test and the W-test. The algorithm was applied on real exome sequencing data of hypertensive disorder, and identified biologically relevant genetic markers to metabolic disorders that were undetectable by a gene-based method. The proposed algorithm is an efficient and powerful tool to enhance the power of association study for whole exome or genome sequencing data.

Availability and Implementation

The ZFA software is available at: http://www2.ccrb.cuhk.edu.hk/statgene/software.html.

Contact

maggiew@cuhk.edu.hk or bzee@cuhk.edu.hk.

Supplementary information

Supplementary data are available at Bioinformatics online.

The untold stories of the speech gene, the FOXP2 cancer gene

FOXP2 encodes a transcription factor involved in speech and language acquisition. Growing evidence now suggests that dysregulated FOXP2 activity may also be instrumental in human oncogenesis, along the lines of other cardinal developmental transcription factors such as DLX5 and DLX6 [1-4]. Several FOXP familymembers are directly involved during cancer initiation, maintenance and progression in the adult [5-8]. This may comprise either a pro-oncogenic activity or a deficient tumor-suppressor role, depending upon cell types and associated signaling pathways. While FOXP2 is expressed in numerous cell types, its expression has been found to be down-regulated in breast cancer [9], hepatocellular carcinoma [8] and gastric cancer biopsies [10]. Conversely, overexpressed FOXP2 has been reported in multiple myelomas, MGUS (Monoclonal Gammopathy of Undetermined Significance), several subtypes of lymphomas [5,11], as well as in neuroblastomas [12] and ERG fusion-negative prostate cancers [13]. According to functional evidences reported in breast cancer [9] and survey of recent transcriptomic and proteomic analyses of different tumor biopsies, we postulate that FOXP2 dysregulation may play a main role throughout cancer initiation and progression. In some cancer conditions, FOXP2 levels are now considered as a critical diagnostic marker of neoplastic cells, and in many situations, they even bear strong prognostic value [5]. Whether FOXP2 may further become a therapeutic target is an actively explored lead. Knowledge reviewed here may help improve our understanding of FOXP2 roles during oncogenesis and provide cues for diagnostic, prognostic and therapeutic analyses.

Opbp is a new architectural/insulator protein required for ribosomal gene expression

A special class of poorly characterized architectural proteins is required for chromatin topology and enhancer–promoter interactions. Here, we identify Opbp as a new Drosophila architectural protein, interacting with CP190 both in vivo and in vitro. Opbp binds to a very restrictive set of genomic regions, through a rare sequence specific motif. These sites are co-bound by CP190 in vivo, and generally located at bidirectional promoters of ribosomal protein genes. We show that Opbp is essential for viability, and loss of opbp function, or destruction of its motif, leads to reduced ribosomal protein gene expression, indicating a functional role in promoter activation. As characteristic of architectural/insulator proteins, the Opbp motif is sufficient for distance-dependent reporter gene activation and enhancer-blocking activity, suggesting an Opbp-mediated enhancer–promoter interaction. Rather than having a constitutive role, Opbp represents a new type of architectural protein with a very restricted, yet essential, function in regulation of housekeeping gene expression.

Epigenetic regulation of noncoding RNA transcription by mammalian RNA polymerase III

RNA polymerase III (Pol III) synthesizes a range of medium-sized noncoding RNAs (collectively 'Pol III genes') whose early established biological roles were so essential that they were considered 'housekeeping genes'. Besides these fundamental functions, diverse unconventional roles of mammalian Pol III genes have recently been recognized and their expression must be exquisitely controlled. In this review, we summarize the epigenetic regulation of Pol III genes by chromatin structure, histone modification and CpG DNA methylation. We also recapitulate the association between dysregulation of Pol III genes and diseases such as cancer and neurological disorders. Additionally, we will discuss why in-depth molecular studies of Pol III genes have not been attempted and how nc886, a Pol III gene, may resolve this issue.

Condensin II is anchored by TFIIIC and H3K4me3 in the mammalian genome and supports the expression of active dense gene clusters

Structural maintenance of chromosome complexes, such as cohesin, have been implicated in a wide variety of chromatin-dependent functions such as genome organization, replication, and gene expression. How these complexes find their sites of association and affect local chromosomal processes is not well understood. We report that condensin II, a complex distinct from cohesin, physically interacts with TFIIIC, and they both colocalize at active gene promoters in the mouse and human genomes, facilitated by interaction between NCAPD3 and the epigenetic mark H3K4me3. Condensin II is important for maintaining high levels of expression of the histone gene clusters as well as the interaction between these clusters in the mouse genome. Our findings suggest that condensin II is anchored to the mammalian genome by a combination of H3K4me3 and the sequence-specific binding of TFIIIC, and that condensin supports the expression of active gene-dense regions found at the boundaries of topological domains. Together, our results support a working model in which condensin II contributes to topological domain boundary-associated gene activity in the mammalian genome.

On the association analysis of genome-sequencing data: A spatial clustering approach for partitioning the entire genome into nonoverlapping windows

For the association analysis of whole-genome sequencing (WGS) studies, we propose an efficient and fast spatial-clustering algorithm. Compared to existing analysis approaches for WGS data, that define the tested regions either by sliding or consecutive windows of fixed sizes along variants, a meaningful grouping of nearby variants into consecutive regions has the advantage that, compared to sliding window approaches, the number of tested regions is likely to be smaller. In comparison to consecutive, fixed-window approaches, our approach is likely to group nearby variants together. Given existing biological evidence that disease-associated mutations tend to physically cluster in specific regions along the chromosome, the identification of meaningful groups of nearby located variants could thus lead to a potential power gain for association analysis. Our algorithm defines consecutive genomic regions based on the physical positions of the variants, assuming an inhomogeneous Poisson process and groups together nearby variants. As parameters are estimated locally, the algorithm takes the differing variant density along the chromosome into account and provides locally optimal partitioning of variants into consecutive regions. An R-implementation of the algorithm is provided. We discuss the theoretical advances of our algorithm compared to existing, window-based approaches and show the performance and advantage of our introduced algorithm in a simulation study and by an application to Alzheimer's disease WGS data. Our analysis identifies a region in the ITGB3 gene that potentially harbors disease susceptibility loci for Alzheimer's disease. The region-based association signal of ITGB3 replicates in an independent data set and achieves formally genome-wide significance. Software Implementation: An implementation of the algorithm in R is available at: https://github.com/heidefier/cluster_wgs_data.

Safe and Effective Gene Therapy for Murine Wiskott-Aldrich Syndrome Using an Insulated Lentiviral Vector

Wiskott-Aldrich syndrome (WAS) is a life-threatening immunodeficiency caused by mutations within the WAS gene. Viral gene therapy to restore WAS protein (WASp) expression in hematopoietic cells of patients with WAS has the potential to improve outcomes relative to the current standard of care, allogeneic bone marrow transplantation. However, the development of viral vectors that are both safe and effective has been problematic. While use of viral transcriptional promoters may increase the risk of insertional mutagenesis, cellular promoters may not achieve WASp expression levels necessary for optimal therapeutic effect. Here we evaluate a self-inactivating (SIN) lentiviral vector combining a chromatin insulator upstream of a viral MND (MPSV LTR, NCR deleted, dl587 PBS) promoter driving WASp expression. Used as a gene therapeutic in Was^−/− mice, this vector resulted in stable WASp⁺ cells in all hematopoietic lineages and rescue of T and B cell defects with a low number of viral integrations per cell, without evidence of insertional mutagenesis in serial bone marrow transplants. In a gene transfer experiment in non-human primates, the insulated MND promoter (driving GFP expression) demonstrated long-term polyclonal engraftment of GFP⁺ cells. These observations demonstrate that the insulated MND promoter safely and efficiently reconstitutes clinically effective WASp expression and should be considered for future WAS therapy.

The Sea Urchin sns5 Chromatin Insulator Shapes the Chromatin Architecture of a Lentivirus Vector Integrated in the Mammalian Genome

Lentivirus vectors are presently the favorite vehicles for therapeutic gene transfer in hematopoietic cells. Nonetheless, these vectors integrate randomly throughout the genome, exhibiting variegation of transgene expression due to the spreading of heterochromatin into the vector sequences. Moreover, the cis-regulatory elements harbored by the vector could disturb the proper transcription of resident genes neighboring the integration site. The incorporation of chromatin insulators in flanking position to the transferred unit can alleviate both the above-mentioned dangerous effects, due to the insulator-specific barrier and enhancer-blocking activities. In this study, we report the valuable properties of the sea urchin-derived sns5 insulator in improving the expression efficiency of a lentivirus vector integrated in the mammalian erythroid genome. We show that these results neither reflect an intrinsic sns5 enhancer activity nor rely on the recruitment of the erythroid-specific GATA-1 factor to sns5. Furthermore, by using the Chromosome Conformation Capture technology, we report that a single copy of the sns5-insulated vector is specifically organized into an independent chromatin loop at the provirus locus. Our results not only provide new clues concerning the molecular mechanism of sns5 function in the erythroid genome but also reassure the use of sns5 to improve the performance of gene therapy vectors.

Architectural proteins Pita, Zw5,and ZIPIC contain homodimerization domain and support specific long-range interactions in Drosophila

According to recent models, as yet poorly studied architectural proteins appear to be required for local regulation of enhancer-promoter interactions, as well as for global chromosome organization. Transcription factors ZIPIC, Pita and Zw5 belong to the class of chromatin insulator proteins and preferentially bind to promoters near the TSS and extensively colocalize with cohesin and condensin complexes. ZIPIC, Pita and Zw5 are structurally similar in containing the N-terminal zinc finger-associated domain (ZAD) and different numbers of C2H2-type zinc fingers at the C-terminus. Here we have shown that the ZAD domains of ZIPIC, Pita and Zw5 form homodimers. In Drosophila transgenic lines, these proteins are able to support long-distance interaction between GAL4 activator and the reporter gene promoter. However, no functional interaction between binding sites for different proteins has been revealed, suggesting that such interactions are highly specific. ZIPIC facilitates long-distance stimulation of the reporter gene by GAL4 activator in yeast model system. Many of the genomic binding sites of ZIPIC, Pita and Zw5 are located at the boundaries of topologically associated domains (TADs). Thus, ZAD-containing zinc-finger proteins can be attributed to the class of architectural proteins.

EAST affects the activity of Su(Hw) insulators by two different mechanisms in Drosophila melanogaster

Recent data suggest that insulators organize chromatin architecture in the nucleus. The best characterized Drosophila insulator, found in the gypsy retrotransposon, contains 12 binding sites for the Su(Hw) protein. Enhancer blocking, along with Su(Hw), requires BTB/POZ domain proteins, Mod(mdg4)-67.2 and CP190. Inactivation of Mod(mdg4)-67.2 leads to a direct repression of the yellow gene promoter by the gypsy insulator. Here, we have shown that such repression is regulated by the level of the EAST protein, which is an essential component of the interchromatin compartment. Deletion of the EAST C-terminal domain suppresses Su(Hw)-mediated repression. Partial inactivation of EAST by mutations in the east gene suppresses the enhancer-blocking activity of the gypsy insulator. The binding of insulator proteins to chromatin is highly sensitive to the level of EAST expression. These results suggest that EAST, one of the main components of the interchromatin compartment, can regulate the activity of chromatin insulators.

Simple Elastic Network Models for Exhaustive Analysis of Long Double-Stranded DNA Dynamics with Sequence Geometry Dependence

Simple elastic network models of DNA were developed to reveal the structure-dynamics relationships for several nucleotide sequences. First, we propose a simple all-atom elastic network model of DNA that can explain the profiles of temperature factors for several crystal structures of DNA. Second, we propose a coarse-grained elastic network model of DNA, where each nucleotide is described only by one node. This model could effectively reproduce the detailed dynamics obtained with the all-atom elastic network model according to the sequence-dependent geometry. Through normal-mode analysis for the coarse-grained elastic network model, we exhaustively analyzed the dynamic features of a large number of long DNA sequences, approximately ∼150 bp in length. These analyses revealed positive correlations between the nucleosome-forming abilities and the inter-strand fluctuation strength of double-stranded DNA for several DNA sequences.

Synthetic promoters in planta

This paper reviews the importance, prospective and development of synthetic promoters reported in planta. A review of the synthetic promoters developed in planta would help researchers utilize the available resources and design new promoters to benefit fundamental research and agricultural applications. The demand for promoters for the improvement and application of transgenic techniques in research and agricultural production is increasing. Native/naturally occurring promoters have some limitations in terms of their induction conditions, transcription efficiency and size. The strength and specificity of native promoter can be tailored by manipulating its 'cis-architecture' by the use of several recombinant DNA technologies. Newly derived chimeric promoters with specific attributes are emerging as an efficient tool for plant molecular biology. In the last three decades, synthetic promoters have been used to regulate plant gene expression. To better understand synthetic promoters, in this article, we reviewed promoter structure, the scope of cis-engineering, strategies for their development, their importance in plant biology and the total number of such promoters (188) developed in planta to date; we then categorized them under different functional regimes as biotic stress-inducible, abiotic stress-inducible, light-responsive, chemical-inducible, hormone-inducible, constitutive and tissue-specific. Furthermore, we identified a set of 36 synthetic promoters that control multiple types of expression in planta. Additionally, we illustrated the differences between native and synthetic promoters and among different synthetic promoter in each group, especially in terms of efficiency and induction conditions. As a prospective of this review, the use of ideal synthetic promoters is one of the prime requirements for generating transgenic plants suitable for promoting sustainable agriculture and plant molecular farming.

EAST Organizes Drosophila Insulator Proteins in the Interchromosomal Nuclear Compartment and Modulates CP190 Binding to Chromatin

Recent data suggest that insulators organize chromatin architecture in the nucleus. The best studied Drosophila insulator proteins, dCTCF (a homolog of the vertebrate insulator protein CTCF) and Su(Hw), are DNA-binding zinc finger proteins. Different isoforms of the BTB-containing protein Mod(mdg4) interact with Su(Hw) and dCTCF. The CP190 protein is a cofactor for the dCTCF and Su(Hw) insulators. CP190 is required for the functional activity of insulator proteins and is involved in the aggregation of the insulator proteins into specific structures named nuclear speckles. Here, we have shown that the nuclear distribution of CP190 is dependent on the level of EAST protein, an essential component of the interchromatin compartment. EAST interacts with CP190 and Mod(mdg4)-67.2 proteins in vitro and in vivo. Over-expression of EAST in S2 cells leads to an extrusion of the CP190 from the insulator bodies containing Su(Hw), Mod(mdg4)-67.2, and dCTCF. In consistent with the role of the insulator bodies in assembly of protein complexes, EAST over-expression led to a striking decrease of the CP190 binding with the dCTCF and Su(Hw) dependent insulators and promoters. These results suggest that EAST is involved in the regulation of CP190 nuclear localization.

The identification of cis-regulatory elements: A review from a machine learning perspective

The majority of the human genome consists of non-coding regions that have been called junk DNA. However, recent studies have unveiled that these regions contain cis-regulatory elements, such as promoters, enhancers, silencers, insulators, etc. These regulatory elements can play crucial roles in controlling gene expressions in specific cell types, conditions, and developmental stages. Disruption to these regions could contribute to phenotype changes. Precisely identifying regulatory elements is key to deciphering the mechanisms underlying transcriptional regulation. Cis-regulatory events are complex processes that involve chromatin accessibility, transcription factor binding, DNA methylation, histone modifications, and the interactions between them. The development of next-generation sequencing techniques has allowed us to capture these genomic features in depth. Applied analysis of genome sequences for clinical genetics has increased the urgency for detecting these regions. However, the complexity of cis-regulatory events and the deluge of sequencing data require accurate and efficient computational approaches, in particular, machine learning techniques. In this review, we describe machine learning approaches for predicting transcription factor binding sites, enhancers, and promoters, primarily driven by next-generation sequencing data. Data sources are provided in order to facilitate testing of novel methods. The purpose of this review is to attract computational experts and data scientists to advance this field.

Boundary Associated Long Noncoding RNA Mediates Long-Range Chromosomal Interactions

CCCTC binding factor (CTCF) is involved in organizing chromosomes into mega base-sized, topologically associated domains (TADs) along with other factors that define sub-TAD organization. CTCF-Cohesin interactions have been shown to be critical for transcription insulation activity as it stabilizes long-range interactions to promote proper gene expression. Previous studies suggest that heterochromatin boundary activity of CTCF may be independent of Cohesin, and there may be additional mechanisms for defining topological domains. Here, we show that a boundary site we previously identified known as CTCF binding site 5 (CBS5) from the homeotic gene cluster A (HOXA) locus exhibits robust promoter activity. This promoter activity from the CBS5 boundary element generates a long noncoding RNA that we designate as boundary associated long noncoding RNA-1 (blncRNA1). Functional characterization of this RNA suggests that the transcript stabilizes long-range interactions at the HOXA locus and promotes proper expression of HOXA genes. Additionally, our functional analysis also shows that this RNA is not needed in the stabilization of CTCF-Cohesin interactions however CTCF-Cohesin interactions are critical in the transcription of blncRNA1. Thus, the CTCF-associated boundary element, CBS5, employs both Cohesin and noncoding RNA to establish and maintain topologically associated domains at the HOXA locus.

Functional role of dimerization and CP190 interacting domains of CTCF protein in Drosophila melanogaster

Background

Insulators play a central role in gene regulation, chromosomal architecture and genome function in higher eukaryotes. To learn more about how insulators carry out their diverse functions, we have begun an analysis of the Drosophila CTCF (dCTCF). CTCF is one of the few insulator proteins known to be conserved from flies to man.

Results

In the studies reported here we have focused on the identification and characterization of two dCTCF protein interaction modules. The first mediates dCTCF multimerization, while the second mediates dCTCF–CP190 interactions. The multimerization domain maps in the N-terminus of the dCTCF protein and likely mediates the formation of tetrameric complexes. The CP190 interaction module encompasses a sequence ~200 amino acids long that spans the C-terminal and mediates interactions with the N-terminal BTB domain of the CP190 protein. Transgene rescue experiments showed that a dCTCF protein lacking sequences critical for CP190 interactions was almost as effective as wild type in rescuing the phenotypic effects of a dCTCF null allele. The mutation did, however, affect CP190 recruitment to specific Drosophila insulator elements and had a modest effect on dCTCF chromatin association. A protein lacking the N-terminal dCTCF multimerization domain incompletely rescued the zygotic and maternal effect lethality of the null and did not rescue the defects in Abd-B regulation evident in surviving adult dCTCF mutant flies. Finally, we show that elimination of maternally contributed dCTCF at the onset of embryogenesis has quite different effects on development and Abd-B regulation than is observed when the homozygous mutant animals develop in the presence of maternally derived dCTCF activity.

Conclusions

Our results indicate that dCTCF–CP190 interactions are less critical for the in vivo functions of the dCTCF protein than the N-terminal dCTCF–dCTCF interaction domain. We also show that the phenotypic consequences of dCTCF mutations differ depending upon when and how dCTCF activity is lost.

Electronic supplementary material

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

Why the activity of a gene depends on its neighbors

Sixty years ago, the position of a gene on a chromosome was seen to be a major determinant of gene activity; however, position effects are rarely central to current discussions of gene expression. We describe a comprehensive and simplifying view of how position in 1D sequence and 3D nuclear space underlies expression. We suggest that apparently-different regulatory motifs including enhancers, silencers, insulators, barriers, and boundaries act similarly - they are active promoters that tether target genes close to, or distant from, appropriate transcription sites or 'factories'. We also suggest that any active transcription unit regulates the firing of its neighbors - and thus can be categorized as one or other type of motif; this is consistent with expression quantitative trait loci (eQTLs) being widely dispersed.

A unified architecture of transcriptional regulatory elements

Gene expression is precisely controlled in time and space through the integration of signals that act at gene promoters and gene-distal enhancers. Classically, promoters and enhancers are considered separate classes of regulatory elements, often distinguished by histone modifications. However, recent studies have revealed broad similarities between enhancers and promoters, blurring the distinction: active enhancers often initiate transcription, and some gene promoters have the potential to enhance transcriptional output of other promoters. Here, we propose a model in which promoters and enhancers are considered a single class of functional element, with a unified architecture for transcription initiation. The context of interacting regulatory elements and the surrounding sequences determine local transcriptional output as well as the enhancer and promoter activities of individual elements.

Promoter or enhancer, what's the difference? Deconstruction of established distinctions and presentation of a unifying model

Gene transcription is strictly controlled by the interplay of regulatory events at gene promoters and gene-distal regulatory elements called enhancers. Despite extensive studies of enhancers, we still have a very limited understanding of their mechanisms of action and their restricted spatio-temporal activities. A better understanding would ultimately lead to fundamental insights into the control of gene transcription and the action of regulatory genetic variants involved in disease. Here, I review and discuss pros and cons of state-of-the-art genomics methods to localize and infer the activity of enhancers. Among the different approaches, profiling of enhancer RNAs yields the highest specificity and may be superior in detecting in vivo activity. I discuss their apparent similarities to promoters, which challenge the established view of enhancers and promoters as distinct entities, and present a unifying model of regulatory elements in transcriptional regulation, in which activity, transcriptional output and regulatory function is context specific.