Transcriptional enhancers in animal development and evolution

Predicting unrecognized enhancer-mediated genome topology by an ensemble machine learning model

Transcriptional enhancers commonly work over long genomic distances to precisely regulate spatiotemporal gene expression patterns. Dissecting the promoters physically contacted by these distal regulatory elements is essential for understanding developmental processes as well as the role of disease-associated risk variants. Modern proximity-ligation assays, like HiChIP and ChIA-PET, facilitate the accurate identification of long-range contacts between enhancers and promoters. However, these assays are technically challenging, expensive, and time-consuming, making it difficult to investigate enhancer topologies, especially in uncharacterized cell types. To overcome these shortcomings, we therefore designed LoopPredictor, an ensemble machine learning model, to predict genome topology for cell types which lack long-range contact maps. To enrich for functional enhancer-promoter loops over common structural genomic contacts, we trained LoopPredictor with both H3K27ac and YY1 HiChIP data. Moreover, the integration of several related multi-omics features facilitated identifying and annotating the predicted loops. LoopPredictor is able to efficiently identify cell type–specific enhancer-mediated loops, and promoter–promoter interactions, with a modest feature input requirement. Comparable to experimentally generated H3K27ac HiChIP data, we found that LoopPredictor was able to identify functional enhancer loops. Furthermore, to explore the cross-species prediction capability of LoopPredictor, we fed mouse multi-omics features into a model trained on human data and found that the predicted enhancer loops outputs were highly conserved. LoopPredictor enables the dissection of cell type–specific long-range gene regulation and can accelerate the identification of distal disease-associated risk variants.

Distinctive regulatory architectures of germline-active and somatic genes in C. elegans

RNA profiling has provided increasingly detailed knowledge of gene expression patterns, yet the different regulatory architectures that drive them are not well understood. To address this, we profiled and compared transcriptional and regulatory element activities across five tissues of Caenorhabditis elegans, covering ∼90% of cells. We find that the majority of promoters and enhancers have tissue-specific accessibility, and we discover regulatory grammars associated with ubiquitous, germline, and somatic tissue–specific gene expression patterns. In addition, we find that germline-active and soma-specific promoters have distinct features. Germline-active promoters have well-positioned +1 and −1 nucleosomes associated with a periodic 10-bp WW signal (W = A/T). Somatic tissue–specific promoters lack positioned nucleosomes and this signal, have wide nucleosome-depleted regions, and are more enriched for core promoter elements, which largely differ between tissues. We observe the 10-bp periodic WW signal at ubiquitous promoters in other animals, suggesting it is an ancient conserved signal. Our results show fundamental differences in regulatory architectures of germline and somatic tissue–specific genes, uncover regulatory rules for generating diverse gene expression patterns, and provide a tissue-specific resource for future studies.

Quantitative dissection of transcription in development yields evidence for transcription-factor-driven chromatin accessibility

Thermodynamic models of gene regulation can predict transcriptional regulation in bacteria, but in eukaryotes, chromatin accessibility and energy expenditure may call for a different framework. Here, we systematically tested the predictive power of models of DNA accessibility based on the Monod-Wyman-Changeux (MWC) model of allostery, which posits that chromatin fluctuates between accessible and inaccessible states. We dissected the regulatory dynamics of hunchback by the activator Bicoid and the pioneer-like transcription factor Zelda in living Drosophila embryos and showed that no thermodynamic or non-equilibrium MWC model can recapitulate hunchback transcription. Therefore, we explored a model where DNA accessibility is not the result of thermal fluctuations but is catalyzed by Bicoid and Zelda, possibly through histone acetylation, and found that this model can predict hunchback dynamics. Thus, our theory-experiment dialogue uncovered potential molecular mechanisms of transcriptional regulatory dynamics, a key step toward reaching a predictive understanding of developmental decision-making.

Epigenomic landscape of enhancer elements during Hydra head organizer formation

BACKGROUND

Axis patterning during development is accompanied by large-scale gene expression changes. These are brought about by changes in the histone modifications leading to dynamic alterations in chromatin architecture. The cis regulatory DNA elements also play an important role towards modulating gene expression in a context-dependent manner. Hydra belongs to the phylum Cnidaria where the first asymmetry in the body plan was observed and the oral-aboral axis originated. Wnt signaling has been shown to determine the head organizer function in the basal metazoan Hydra.

RESULTS

To gain insights into the evolution of cis regulatory elements and associated chromatin signatures, we ectopically activated the Wnt signaling pathway in Hydra and monitored the genome-wide alterations in key histone modifications. Motif analysis of putative intergenic enhancer elements from Hydra revealed the conservation of bilaterian cis regulatory elements that play critical roles in development. Differentially regulated enhancer elements were identified upon ectopic activation of Wnt signaling and found to regulate many head organizer specific genes. Enhancer activity of many of the identified cis regulatory elements was confirmed by luciferase reporter assay. Quantitative chromatin immunoprecipitation analysis upon activation of Wnt signaling further confirmed the enrichment of H3K27ac on the enhancer elements of Hv_Wnt5a, Hv_Wnt11 and head organizer genes Hv_Bra1, CnGsc and Hv_Pitx1. Additionally, perturbation of the putative H3K27me3 eraser activity using a specific inhibitor affected the ectopic activation of Wnt signaling indicating the importance of the dynamic changes in the H3K27 modifications towards regulation of the genes involved in the head organizer activity.

CONCLUSIONS

The activation-associated histone marks H3K4me3, H3K27ac and H3K9ac mark chromatin in a similar manner as seen in bilaterians. We identified intergenic cis regulatory elements which harbor sites for key transcription factors involved in developmental processes. Differentially regulated enhancers exhibited motifs for many zinc-finger, T-box and ETS related TFs whose homologs have a head specific expression in Hydra and could be a part of the pioneer TF network in the patterning of the head. The ability to differentially modify the H3K27 residue is critical for the patterning of Hydra axis revealing a dynamic acetylation/methylation switch to regulate gene expression and chromatin architecture.

Transcriptional enhancers: from prediction to functional assessment on a genome-wide scale

Enhancers are cis-regulatory sequences located distally to target genes. These sequences consolidate developmental and environmental cues to coordinate gene expression in a tissue-specific manner. Enhancer function and tissue specificity depend on the expressed set of transcription factors, which recognize binding sites and recruit cofactors that regulate local chromatin organization and gene transcription. Unlike other genomic elements, enhancers are challenging to identify because they function independently of orientation, are often distant from their promoters, have poorly defined boundaries, and display no reading frame. In addition, there are no defined genetic or epigenetic features that are unambiguously associated with enhancer activity. Over recent years there have been developments in both empirical assays and computational methods for enhancer prediction. We review genome-wide tools, CRISPR advancements, and high-throughput screening approaches that have improved our ability to both observe and manipulate enhancers in vitro at the level of primary genetic sequences, chromatin states, and spatial interactions. We also highlight contemporary animal models and their importance to enhancer validation. Together, these experimental systems and techniques complement one another and broaden our understanding of enhancer function in development, evolution, and disease.

Single-gene imaging links genome topology, promoter-enhancer communication and transcription control

Transcription activation by distal enhancers is essential for cell-fate specification and maintenance of cellular identities. How long-range gene regulation is physically achieved, especially within complex regulatory landscapes of non-binary enhancer-promoter configurations, remains elusive. Recent nanoscopy advances have quantitatively linked promoter kinetics and ~100- to 200-nm-sized clusters of enhancer-associated regulatory factors (RFs) at important developmental genes. Here, we further dissect mechanisms of RF clustering and transcription activation in mouse embryonic stem cells. RF recruitment into clusters involves specific molecular recognition of cognate DNA and chromatin-binding sites, suggesting underlying cis-element clustering. Strikingly, imaging of tagged genomic loci, with ≤1 kilobase and ~20-nanometer precision, in live cells, reveals distal enhancer clusters over the extended locus in frequent close proximity to target genes-within RF-clustering distances. These high-interaction-frequency enhancer-cluster 'superclusters' create nano-environments wherein clustered RFs activate target genes, providing a structural framework for relating genome organization, focal RF accumulation and transcription activation.

The Drosophila MLR COMPASS-like complex regulates bantam miRNA expression differentially in the context of cell fate

The conserved MLR COMPASS-like complexes are histone modifiers that are recruited by a variety of transcription factors to enhancer regions where they act as necessary epigenetic tools for enhancer establishment and function. A critical in vivo target of the Drosophila MLR complex is the bantam miRNA that regulates cell survival and functions in feedback regulation of cellular signaling pathways during development. We determine that loss of Drosophila MLR complex function in developing wing and eye imaginal discs results in growth and patterning defects that are sensitive to bantam levels. Consistent with an essential regulatory role in modulating bantam transcription, the MLR complex binds to tissue-specific bantam enhancers and contributes to fine-tuning expression levels during larval tissue development. In wing imaginal discs, the MLR complex attenuates bantam enhancer activity by negatively regulating expression; whereas, in differentiating eye discs, the complex exerts either positive or negative regulatory activity on bantam transcription depending on cell fate. Furthermore, while the MLR complex is not required to control bantam levels in undifferentiated eye cells anterior to the morphogenetic furrow, it serves to prepare critical enhancer control of bantam transcription for later regulation upon differentiation. Our investigation into the transcriptional regulation of a single target in a developmental context has provided novel insights as to how the MLR complex contributes to the precise timing of gene expression, and how the complex functions to help orchestrate the regulatory output of conserved signaling pathways during animal development.

Conjugated activation of myocardial-specific transcription of Gja5 by a pair of Nkx2-5-Shox2 co-responsive elements

The sinoatrial node (SAN) is the primary pacemaker in the heart. During cardiogenesis, Shox2 and Nkx2-5 are co-expressed in the junction domain of the SAN and regulate pacemaker cell fate through a Shox2-Nkx2-5 antagonism. Cx40 is a marker of working myocardium and an Nkx2-5 transcriptional output antagonized by Shox2, but the underlying regulatory mechanisms remain elusive. Here we characterized a bona fide myocardial-specific Gja5 (coding gene of Cx40) distal enhancer consisting of a pair of Nkx2-5 and Shox2 co-bound elements in the regulatory region of Gja5. Transgenic reporter assays revealed that neither element alone, but the conjugation of both elements together, drives myocardial-specific transcription. Genetic analyses confirmed that the activation of this enhancer depends on Nkx2-5 but is inhibited by Shox2 in vivo, and its presence is essential for Gja5 expression in the myocardium but not the endothelial cells of the heart. Furthermore, chromatin conformation analysis showed an Nkx2-5-dependent loop formation between these two elements and the Gja5 promoter in vivo, indicating that Nkx2-5 bridges the conjugated activation of this enhancer by pairing the two elements to the Gja5 promoter.

The function of architecture and logic in developmental gene regulatory networks

An important contribution of systems biology is the insight that biological systems depend on the function of molecular interactions and not just on individual molecules. System level mechanisms are particularly important in the development of animals and plants which depends not just on transcription factors and signaling molecules, but also on regulatory circuits and gene regulatory networks (GRNs). However, since GRNs consist of transcription factors, it can be challenging to assess the function of regulatory circuits independently of the function of regulatory factors. The comparison of different GRNs offers a way to do so and leads to several observations. First, similar regulatory circuits operate in various developmental contexts and in different species, and frequently, these circuits are associated with similar developmental functions. Second, given regulatory circuits are often used at particular positions within the GRN hierarchy. Third, in some GRNs, regulatory circuits are organized in a particular order in respect to each other. And fourth, the evolution of GRNs occurs not just by co-option of regulatory genes but also by rewiring of regulatory linkages between conserved regulatory genes, indicating that the organization of interactions is important. Thus, even though in most instances the function of regulatory circuits remains to be discovered, it becomes evident that the architecture and logic of GRNs are functionally important for the control of genome activity and for the specification of the body plan.

Super-enhancers in transcriptional regulation and genome organization

Gene expression is precisely controlled in a stage and cell-type-specific manner, largely through the interaction between cis-regulatory elements and their associated trans-acting factors. Where these components aggregate in promoters and enhancers, they are able to cooperate to modulate chromatin structure and support the engagement in long-range 3D superstructures that shape the dynamics of a cell's genomic architecture. Recently, the term 'super-enhancer' has been introduced to describe a hyper-active regulatory domain comprising a complex array of sequence elements that work together to control the key gene networks involved in cell identity. Here, we survey the unique characteristics of super-enhancers compared to other enhancer types and summarize the recent advances in our understanding of their biological role in gene regulation. In particular, we discuss their capacity to attract the formation of phase-separated condensates, and capacity to generate three-dimensional genome structures that precisely activate their target genes. We also propose a multi-stage transition model to explain the evolutionary pressure driving the development of super-enhancers in complex organisms, and highlight the potential for involvement in tumorigenesis. Finally, we discuss more broadly the role of super-enhancers in human health disorders and related potential in therapeutic interventions.

On the specificity of gene regulatory networks: How does network co-option affect subsequent evolution?

The process of multicellular organismal development hinges upon the specificity of developmental programs: for different parts of the organism to form unique features, processes must exist to specify each part. This specificity is thought to be hardwired into gene regulatory networks, which activate cohorts of genes in particular tissues at particular times during development. However, the evolution of gene regulatory networks sometimes occurs by mechanisms that sacrifice specificity. One such mechanism is network co-option, in which existing gene networks are redeployed in new developmental contexts. While network co-option may offer an efficient mechanism for generating novel phenotypes, losses of tissue specificity at redeployed network genes could restrict the ability of the affected traits to evolve independently. At present, there has not been a detailed discussion regarding how tissue specificity of network genes might be altered due to gene network co-option at its initiation, as well as how trait independence can be retained or restored after network co-option. A lack of clarity about network co-option makes it more difficult to speculate on the long-term evolutionary implications of this mechanism. In this review, we will discuss the possible initial outcomes of network co-option, outline the mechanisms by which networks may retain or subsequently regain specificity after network co-option, and comment on some of the possible evolutionary consequences of network co-option. We place special emphasis on the need to consider selectively-neutral outcomes of network co-option to improve our understanding of the role of this mechanism in trait evolution.

Functional Conservation of Divergent p63-Bound cis-Regulatory Elements

The transcription factor p63 is an essential regulator of vertebrate ectoderm development, including epidermis, limbs, and craniofacial tissues. Here, we have investigated the evolutionary conservation of p63 binding sites (BSs) between zebrafish and human. First, we have analyzed sequence conservation of p63 BSs by comparing ChIP-seq data from human keratinocytes and zebrafish embryos, observing a very poor conservation. Next, we compared the gene regulatory network orchestrated by p63 in both species and found a high overlap between them, suggesting a high degree of functional conservation during evolution despite sequence divergence and the large evolutionary distance. Finally, we used transgenic reporter assays in zebrafish embryos to functionally validate a set of equivalent p63 BSs from zebrafish and human located close to genes involved in epidermal development. Reporter expression was driven by human and zebrafish BSs to many common tissues related to p63 expression domains. Therefore, we conclude that the gene regulatory network controlled by p63 is highly conserved across vertebrates despite the fact that p63-bound regulatory elements show high divergence.

Disparate expression specificities coded by a shared Hox-C enhancer

Can a single regulatory sequence be shared by two genes undergoing functional divergence? Here we describe a single promiscuous enhancer within the Drosophila Antennapedia Complex, EO053, that directs aspects of the expression of two adjacent genes, pb (a Hox2 ortholog) and zen2 (a divergent Hox3 paralog), with disparate spatial and temporal expression patterns. We were unable to separate the pb-like and zen2-like specificities within EO053, and we identify sequences affecting both expression patterns. Importantly, genomic deletion experiments demonstrate that EO053 cooperates with additional pb- and zen2-specific enhancers to regulate the mRNA expression of both genes. We examine sequence conservation of EO053 within the Schizophora, and show that patterns of synteny between the Hox2 and Hox3 orthologs in Arthropods are consistent with a shared regulatory relationship extending prior to the Hox3/zen divergence. Thus, EO053 represents an example of two genes having evolved disparate outputs while utilizing this shared regulatory region.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Functional genomic approaches to elucidate the role of enhancers during development

Successful development depends on the precise tissue-specific regulation of genes by enhancers, genetic elements that act as switches to control when and where genes are expressed. Because enhancers are critical for development, and the majority of disease-associated mutations reside within enhancers, it is essential to understand which sequences within enhancers are important for function. Advances in sequencing technology have enabled the rapid generation of genomic data that predict putative active enhancers, but functionally validating these sequences at scale remains a fundamental challenge. Herein, we discuss the power of genome-wide strategies used to identify candidate enhancers, and also highlight limitations and misconceptions that have arisen from these data. We discuss the use of massively parallel reporter assays to test enhancers for function at scale. We also review recent advances in our ability to study gene regulation during development, including CRISPR-based tools to manipulate genomes and single-cell transcriptomics to finely map gene expression. Finally, we look ahead to a synthesis of complementary genomic approaches that will advance our understanding of enhancer function during development. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Developmental Biology > Developmental Processes in Health and Disease Laboratory Methods and Technologies > Genetic/Genomic Methods.

The notochord gene regulatory network in chordate evolution: Conservation and divergence from Ciona to vertebrates

The notochord is a structure required for support and patterning of all chordate embryos, from sea squirts to humans. An increasing amount of information on notochord development and on the molecular strategies that ensure its proper morphogenesis has been gleaned through studies in the sea squirt Ciona. This invertebrate chordate offers a fortunate combination of experimental advantages, ranging from translucent, fast-developing embryos to a compact genome and impressive biomolecular resources. These assets have enabled the rapid identification of numerous notochord genes and cis-regulatory regions, and provide a rather unique opportunity to reconstruct the gene regulatory network that controls the formation of this developmental and evolutionary chordate landmark. This chapter summarizes the morphogenetic milestones that punctuate notochord formation in Ciona, their molecular effectors, and the current knowledge of the gene regulatory network that ensures the accurate spatial and temporal orchestration of these processes.

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.

Global Analysis of Enhancer Targets Reveals Convergent Enhancer-Driven Regulatory Modules

Single-cell screens enable high-throughput functional assessment of enhancers in their endogenous genomic context. However, the design of current studies limits their application to identifying the primary gene targets of enhancers. Here, we improve the experimental and computational parameters of single-cell enhancer screens to identify the secondary gene targets of enhancers. Our analysis of >500 putative enhancers in K562 cells reveals an interwoven enhancer-driven gene regulatory network. We find that enhancers from distinct genomic loci converge to modulate the expression of common sub-modules, including the α- and β-globin loci, by directly regulating transcription factors. Our analysis suggests that several genetic variants associated with myeloid blood cell traits alter the activity of a distal enhancer of MYB (~140 kb away), with downstream consequences on hemoglobin genes expression and cell state. These data have implications for the understanding of enhancer-associated traits and emphasize the flexibility of controlling transcriptional systems by modifying enhancer activity.

Xie et al. apply improved strategies for single-cell screens to identify an enhancer-driven transcriptional regulatory network in K562 cells. They demonstrate that the same group of genes can be indirectly regulated by enhancers from distinct genomic loci. These data have implications for the understanding of enhancer-associated traits.

Evolution: How Many Phenotypes Do Regulatory Mutations Affect?

Mutations in gene regulatory regions are thought to play an important role in the evolution of morphological structures. This is largely due to their minimal pleiotropic effects, limiting their impact to one particular body part. A recent study finds that one such regulatory mutation may affect two particular morphological structures.

DNA Methylation of Enhancer Elements in Myeloid Neoplasms: Think Outside the Promoters?

Gene regulation through DNA methylation is a well described phenomenon that has a prominent role in physiological and pathological cell-states. This epigenetic modification is usually grouped in regions denominated CpG islands, which frequently co-localize with gene promoters, silencing the transcription of those genes. Recent genome-wide DNA methylation studies have challenged this paradigm, demonstrating that DNA methylation of regulatory regions outside promoters is able to influence cell-type specific gene expression programs under physiologic or pathologic conditions. Coupling genome-wide DNA methylation assays with histone mark annotation has allowed for the identification of specific epigenomic changes that affect enhancer regulatory regions, revealing an additional layer of complexity to the epigenetic regulation of gene expression. In this review, we summarize the novel evidence for the molecular and biological regulation of DNA methylation in enhancer regions and the dynamism of these changes contributing to the fine-tuning of gene expression. We also analyze the contribution of enhancer DNA methylation on the expression of relevant genes in acute myeloid leukemia and chronic myeloproliferative neoplasms. The characterization of the aberrant enhancer DNA methylation provides not only a novel pathogenic mechanism for different tumors but also highlights novel potential therapeutic targets for myeloid derived neoplasms.

QTL mapping of rat blood pressure loci on RNO1 within a homologous region linked to human hypertension on HSA15

Fine-mapping of regions linked to the inheritance of hypertension is accomplished by genetic dissection of blood pressure quantitative trait loci (BP QTLs) in rats. The goal of the current study was to further fine-map two genomic regions on rat chromosome 1 with opposing blood pressure effects (BP QTL1b1 and BP QTL1b1a), the homologous region of which on human chromosome 15 harbors BP QTLs. Two new substrains were constructed and studied from the previously reported BP QTL1b1, one having significantly lower systolic BP by 17 mmHg than that of the salt-sensitive (S) rat (P = 0.007). The new limits of BP QTL1b1 were between 134.09 Mb and 135.40 Mb with a 43% improvement from the previous 2.31 Mb to the current 1.31 Mb interval containing 4 protein-coding genes (Rgma, Chd2, Fam174b, and St8sia2), 2 predicted miRNAs, and 4 lncRNAs. One new substrain was constructed and studied from the previously reported BPQTL1b1a having a significantly higher systolic BP by 22 mmHg (P = 0.006) than that of the S rat. The new limits of BPQTL1b1a were between 133.53 Mb and 134.52 Mb with a 32% improvement from the previous1.45 Mb to the current 990.21 Kb interval containing 1 protein-coding gene, Mctp2, and a lncRNA. The congenic segments of these two BP QTLs overlapped between 134.09 Mb and 134.52 Mb. No exonic variants were detected in any of the genes. These findings reiterate complexity of genetic regulation of BP within QTL regions, where elements beyond protein-coding sequences could be factors in controlling BP.

Transcriptional control by enhancers and enhancer RNAs

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

Reorganization of 3D genome structure may contribute to gene regulatory evolution in primates

A growing body of evidence supports the notion that variation in gene regulation plays a crucial role in both speciation and adaptation. However, a comprehensive functional understanding of the mechanisms underlying regulatory evolution remains elusive. In primates, one of the crucial missing pieces of information towards a better understanding of regulatory evolution is a comparative annotation of interactions between distal regulatory elements and promoters. Chromatin conformation capture technologies have enabled genome-wide quantifications of such distal 3D interactions. However, relatively little comparative research in primates has been done using such technologies. To address this gap, we used Hi-C to characterize 3D chromatin interactions in induced pluripotent stem cells (iPSCs) from humans and chimpanzees. We also used RNA-seq to collect gene expression data from the same lines. We generally observed that lower-order, pairwise 3D genomic interactions are conserved in humans and chimpanzees, but higher order genomic structures, such as topologically associating domains (TADs), are not as conserved. Inter-species differences in 3D genomic interactions are often associated with gene expression differences between the species. To provide additional functional context to our observations, we considered previously published chromatin data from human stem cells. We found that inter-species differences in 3D genomic interactions, which are also associated with gene expression differences between the species, are enriched for both active and repressive marks. Overall, our data demonstrate that, as expected, an understanding of 3D genome reorganization is key to explaining regulatory evolution.

The way in which a genome folds affects the regulation of gene expression. This is often due to loops in the three-dimensional structure that bring linearly distant genes and regulatory elements into close proximity. Most studies examining three-dimensional structure genome-wide are limited to a single species. In this study, we compared three-dimensional structure in the genomes of induced pluripotent stem cells from humans and chimpanzees. We collected gene expression data from the same samples, which allowed us to assess the contribution of three-dimensional chromatin conformation to gene regulatory evolution in primates. Our results demonstrate that gene expression differences between the species may often be mediated by differences in three-dimensional genomic interactions. Our data also suggest that large-scale chromatin structures (i.e. topologically associating domains, TADs) are not well conserved in their placement across species. We hope the analytical paradigms we present here could serve as a basis for future comparative studies of three-dimensional genome organization, elucidating the putative functional regulatory loci driving speciation.

Highly rearranged chromosomes reveal uncoupling between genome topology and gene expression

Chromatin topology is intricately linked to gene expression, yet its functional requirement remains unclear. Here, we comprehensively assessed the interplay between genome topology and gene expression using highly rearranged chromosomes (balancers) spanning ~75% of the Drosophila genome. Using transheterozyte (balancer/wild-type) embryos, we measured allele-specific changes in topology and gene expression in cis, whilst minimizing trans effects. Through genome sequencing, we resolved eight large nested inversions, smaller inversions, duplications, and thousands of deletions. These extensive rearrangements caused many changes to chromatin topology, including long-range loops, TADs and promoter interactions, yet these are not predictive of changes in expression. Gene expression is generally not altered around inversion breakpoints, indicating that mis-appropriate enhancer-promoter activation is a rare event. Similarly, shuffling or fusing TADs, changing intra-TAD connections and disrupting long-range inter-TAD loops, does not alter expression for the majority of genes. Our results suggest that properties other than chromatin topology ensure productive enhancer-promoter interactions.

Dissecting the sharp response of a canonical developmental enhancer reveals multiple sources of cooperativity

Developmental enhancers integrate graded concentrations of transcription factors (TFs) to create sharp gene expression boundaries. Here we examine the hunchback P2 (HbP2) enhancer which drives a sharp expression pattern in the Drosophila blastoderm embryo in response to the transcriptional activator Bicoid (Bcd). We systematically interrogate cis and trans factors that influence the shape and position of expression driven by HbP2, and find that the prevailing model, based on pairwise cooperative binding of Bcd to HbP2 is not adequate. We demonstrate that other proteins, such as pioneer factors, Mediator and histone modifiers influence the shape and position of the HbP2 expression pattern. Comparing our results to theory reveals how higher-order cooperativity and energy expenditure impact boundary location and sharpness. Our results emphasize that the bacterial view of transcription regulation, where pairwise interactions between regulatory proteins dominate, must be reexamined in animals, where multiple molecular mechanisms collaborate to shape the gene regulatory function.

Dedifferentiation by adenovirus E1A due to inactivation of Hippo pathway effectors YAP and TAZ

In this study, Zemke et al. show that E1A inactivates the Hippo pathway-regulated TEAD coactivators YAP and TAZ by causing their sequestration in the cytoplasm. Their findings suggest that YAP/TAZ function in a developmental checkpoint controlled by signaling from the actin cytoskeleton that prevents differentiation of a progenitor cell until it is in the correct cellular and tissue environment.

Adenovirus transformed cells have a dedifferentiated phenotype. Eliminating E1A in transformed human embryonic kidney cells derepressed ∼2600 genes, generating a gene expression profile closely resembling mesenchymal stem cells (MSCs). This was associated with a dramatic change in cell morphology from one with scant cytoplasm and a globular nucleus to one with increased cytoplasm, extensive actin stress fibers, and actomyosin-dependent flattening against the substratum. E1A-induced hypoacetylation at histone H3 Lys27 and Lys18 (H3K27/18) was reversed. Most of the increase in H3K27/18ac was in enhancers near TEAD transcription factors bound by Hippo signaling-regulated coactivators YAP and TAZ. E1A causes YAP/TAZ cytoplasmic sequestration. After eliminating E1A, YAP/TAZ were transported into nuclei, where they associated with poised enhancers with DNA-bound TEAD4 and H3K4me1. This activation of YAP/TAZ required RHO family GTPase signaling and caused histone acetylation by p300/CBP, chromatin remodeling, and cohesin loading to establish MSC-associated enhancers and then superenhancers. Consistent results were also observed in primary rat embryo kidney cells, human fibroblasts, and human respiratory tract epithelial cells. These results together with earlier studies suggest that YAP/TAZ function in a developmental checkpoint controlled by signaling from the actin cytoskeleton that prevents differentiation of a progenitor cell until it is in the correct cellular and tissue environment.

Integrated high-throughput analysis identifies super enhancers associated with chemoresistance in SCLC

BACKGROUND

Chemoresistance is a primary clinical challenge for the management of small cell lung cancer. Additionally, transcriptional regulation by super enhancer (SE) has an important role in tumor evolution. The functions of SEs, a key class of noncoding DNA cis-regulatory elements, have been the subject of many recent studies in the field of cancer research.

METHODS

In this study, using chromatin immunoprecipitation-sequencing and RNA-sequencing (RNA-seq), we aimed to identify SEs associated with chemoresistance from H69AR cells. Through integrated bioinformatics analysis of the MEME chip, we predicted the master transcriptional factors (TFs) binding to SE sites and verified the relationships between TFs of SEs and drug resistance by RNA interference, cell counting kit 8 assays, quantitative real-time reverse transcription polymerase chain reaction.

RESULTS

In total, 108 SEs were screened from H69AR cells. When combining this analysis with RNA-seq data, 45 SEs were suggested to be closely related to drug resistance. Then, 12 master TFs were predicted to localize to regions of those SEs. Subsequently, we selected forkhead box P1 (FOXP1), interferon regulatory factor 1 (IRF1), and specificity protein 1 (SP1) to authenticate the functional relationships of master TFs with chemoresistance via SEs.

CONCLUSIONS

We screened out SEs involved with drug resistance and evaluated the functions of FOXP1, IRF1, and SP1 in chemoresistance. Our findings established a large group of SEs associated with drug resistance in small cell lung cancer, revealed the drug resistance mechanisms of SEs, and provided insights into the clinical applications of SEs.

Identification and Conservation Analysis of Cis-Regulatory Elements in Pig Liver

The liver plays a key role in metabolism and affects pig production. However, the functional annotation of noncoding regions of the pig liver remains poorly understood. We revealed the landscape of cis-regulatory elements and their functional characterization in pig liver. We identified 102,373 cis-regulatory elements in the pig liver, including enhancers, promoters, super-enhancers, and broad H3K4me3 domains, and highlighted 26 core transcription regulatory factors in the pig liver as well. We found similarity of cis-regulatory elements among those of pigs, humans, and cattle. Despite the low proportion of functionally conserved enhancers (~30%) between pig and human liver tissue, ~78% of the pig liver enhancer orthologues sequence could play an enhancer role in other human tissues. Additionally, we observed that the ratio of consistent super-enhancer-associated genes was significantly higher than the ratio of functionally conserved super-enhancers. Approximately 54% of the core regulation factors driven by super-enhancers were consistent across the liver from these three species. Our pig liver annotation and functional characterization studies provide a system and resource for noncoding annotation for future gene regulatory studies in pigs. Furthermore, our study also showed the high level functional conservation of cis-regulatory elements in mammals; it also improved our understanding of regulation function of mammal cis-regulatory elements.

Functional characterization of two enhancers located downstream FOXP2

BACKGROUND

Mutations in the coding region of FOXP2 are known to cause speech and language impairment. However, it is not clear how dysregulation of the gene contributes to language deficit. Interestingly, microdeletions of the region downstream the gene have been associated with cognitive deficits.

METHODS

Here, we investigate changes in FOXP2 expression in the SK-N-MC neuroblastoma human cell line after deletion by CRISPR-Cas9 of two enhancers located downstream of the gene.

RESULTS

Deletion of any of these two functional enhancers downregulates FOXP2, but also upregulates the closest 3' gene MDFIC. Because this effect is not statistically significant in a HEK 293 cell line, derived from the human kidney, both enhancers might confer a tissue specific regulation to both genes. We have also found that the deletion of any of these enhancers downregulates six well-known FOXP2 target genes in the SK-N-MC cell line.

CONCLUSIONS

We expect these findings contribute to a deeper understanding of how FOXP2 and MDFIC are regulated to pace neuronal development supporting cognition, speech and language.

Complex signal processing in synthetic gene circuits using cooperative regulatory assemblies

Eukaryotic genes are regulated by multivalent transcription factor complexes. Through cooperative self-assembly, these complexes perform nonlinear regulatory operations involved in cellular decision-making and signal processing. In this study, we apply this design principle to synthetic networks, testing whether engineered cooperative assemblies can program nonlinear gene circuit behavior in yeast. Using a model-guided approach, we show that specifying the strength and number of assembly subunits enables predictive tuning between linear and nonlinear regulatory responses for single- and multi-input circuits. We demonstrate that assemblies can be adjusted to control circuit dynamics. We harness this capability to engineer circuits that perform dynamic filtering, enabling frequency-dependent decoding in cell populations. Programmable cooperative assembly provides a versatile way to tune the nonlinearity of network connections, markedly expanding the engineerable behaviors available to synthetic circuits.

100 Years of evolving gene-disease complexities and scientific debutants

It's been over 100 years since the word `gene' is around and progressively evolving in several scientific directions. Time-to-time technological advancements have heavily revolutionized the field of genomics, especially when it's about, e.g. triple code development, gene number proposition, genetic mapping, data banks, gene-disease maps, catalogs of human genes and genetic disorders, CRISPR/Cas9, big data and next generation sequencing, etc. In this manuscript, we present the progress of genomics from pea plant genetics to the human genome project and highlight the molecular, technical and computational developments. Studying genome and epigenome led to the fundamentals of development and progression of human diseases, which includes chromosomal, monogenic, multifactorial and mitochondrial diseases. World Health Organization has classified, standardized and maintained all human diseases, when many academic and commercial online systems are sharing information about genes and linking to associated diseases. To efficiently fathom the wealth of this biological data, there is a crucial need to generate appropriate gene annotation repositories and resources. Our focus has been how many gene-disease databases are available worldwide and which sources are authentic, timely updated and recommended for research and clinical purposes. In this manuscript, we have discussed and compared 43 such databases and bioinformatics applications, which enable users to connect, explore and, if possible, download gene-disease data.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

ATAC-seq reveals regional differences in enhancer accessibility during the establishment of spatial coordinates in the Drosophila blastoderm

Establishment of spatial coordinates during Drosophila embryogenesis relies on differential regulatory activity of axis patterning enhancers. Concentration gradients of activator and repressor transcription factors (TFs) provide positional information to each enhancer, which in turn promotes transcription of a target gene in a specific spatial pattern. However, the interplay between an enhancer regulatory activity and its accessibility as determined by local chromatin organization is not well understood. We profiled chromatin accessibility with ATAC-seq in narrow, genetically tagged domains along the antero-posterior axis in the Drosophila blastoderm. We demonstrate that one-quarter of the accessible genome displays significant regional variation in its ATAC-seq signal immediately after zygotic genome activation. Axis patterning enhancers are enriched among the most variable intervals, and their accessibility changes correlate with their regulatory activity. In an embryonic domain where an enhancer receives a net activating TF input and promotes transcription, it displays elevated accessibility in comparison to a domain where it receives a net repressive input. We propose that differential accessibility is a signature of patterning cis-regulatory elements in the Drosophila blastoderm and discuss potential mechanisms by which accessibility of enhancers may be modulated by activator and repressor TFs.

Inference of Developmental Gene Regulatory Networks Beyond Classical Model Systems: New Approaches in the Post-genomic Era

The advent of high-throughput sequencing (HTS) technologies has revolutionized the way we understand the transformation of genetic information into morphological traits. Elucidating the network of interactions between genes that govern cell differentiation through development is one of the core challenges in genome research. These networks are known as developmental gene regulatory networks (dGRNs) and consist largely of the functional linkage between developmental control genes, cis-regulatory modules, and differentiation genes, which generate spatially and temporally refined patterns of gene expression. Over the last 20 years, great advances have been made in determining these gene interactions mainly in classical model systems, including human, mouse, sea urchin, fruit fly, and worm. This has brought about a radical transformation in the fields of developmental biology and evolutionary biology, allowing the generation of high-resolution gene regulatory maps to analyze cell differentiation during animal development. Such maps have enabled the identification of gene regulatory circuits and have led to the development of network inference methods that can recapitulate the differentiation of specific cell-types or developmental stages. In contrast, dGRN research in non-classical model systems has been limited to the identification of developmental control genes via the candidate gene approach and the characterization of their spatiotemporal expression patterns, as well as to the discovery of cis-regulatory modules via patterns of sequence conservation and/or predicted transcription-factor binding sites. However, thanks to the continuous advances in HTS technologies, this scenario is rapidly changing. Here, we give a historical overview on the architecture and elucidation of the dGRNs. Subsequently, we summarize the approaches available to unravel these regulatory networks, highlighting the vast range of possibilities of integrating multiple technical advances and theoretical approaches to expand our understanding on the global gene regulation during animal development in non-classical model systems. Such new knowledge will not only lead to greater insights into the evolution of molecular mechanisms underlying cell identity and animal body plans, but also into the evolution of morphological key innovations in animals.

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.

Redundant and Cryptic Enhancer Activities of the Drosophila yellow Gene

Cis-regulatory sequences known as enhancers play a key role in regulating gene expression. Evolutionary changes in these DNA sequences contribute to phenotypic evolution. The Drosophila yellow gene, which is required for pigmentation, has emerged as a model system for understanding how cis-regulatory sequences evolve, providing some of the most detailed insights available into how activities of orthologous enhancers have diverged between species. Here, we examine the evolution of yellow cis-regulatory sequences on a broader scale, by comparing the distribution and function of yellow enhancer activities throughout the 5' intergenic and intronic sequences of Drosophila melanogaster, D. pseudoobscura, and D. willistoni We find that cis-regulatory sequences driving expression in a particular tissue are not as modular as previously described, but rather have many redundant and cryptic enhancer activities distributed throughout the regions surveyed. Interestingly, cryptic enhancer activities of sequences from one species often drove patterns of expression observed in other species, suggesting that the frequent evolutionary changes in yellow expression observed among Drosophila species may be facilitated by gaining and losing repression of preexisting cis-regulatory sequences.

Genome-wide approaches to unravelling host-virus interactions in Dengue and Zika infections

Genomics approaches are increasingly utilized to probe host-viral interactions and identify mechanisms of viral pathogenesis and host-subversion. Here we review recent studies that utilize Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-Cas9 screens, transcriptomics and epigenomics to gain insight into Dengue and Zika virus infections in humans. We discuss the benefits and limitations of recently utilized techniques that separate virally infected cells from neighboring uninfected cells to identify the mechanisms by which these viruses regulate host responses. We conclude by discussing how these approaches can best advance our understanding of Dengue and Zika virus pathogenesis in humans.

Optimal Decoding of Cellular Identities in a Genetic Network

In developing organisms, spatially prescribed cell identities are thought to be determined by the expression levels of multiple genes. Quantitative tests of this idea, however, require a theoretical framework capable of exposing the rules and precision of cell specification over developmental time. We use the gap gene network in the early fly embryo as an example to show how expression levels of the four gap genes can be jointly decoded into an optimal specification of position with 1% accuracy. The decoder correctly predicts, with no free parameters, the dynamics of pair-rule expression patterns at different developmental time points and in various mutant backgrounds. Precise cellular identities are thus available at the earliest stages of development, contrasting the prevailing view of positional information being slowly refined across successive layers of the patterning network. Our results suggest that developmental enhancers closely approximate a mathematically optimal decoding strategy.

Diverse motif ensembles specify non-redundant DNA binding activities of AP-1 family members in macrophages

Mechanisms by which members of the AP-1 family of transcription factors play non-redundant biological roles despite recognizing the same DNA sequence remain poorly understood. To address this question, here we investigate the molecular functions and genome-wide DNA binding patterns of AP-1 family members in primary and immortalized mouse macrophages. ChIP-sequencing shows overlapping and distinct binding profiles for each factor that were remodeled following TLR4 ligation. Development of a machine learning approach that jointly weighs hundreds of DNA recognition elements yields dozens of motifs predicted to drive factor-specific binding profiles. Machine learning-based predictions are confirmed by analysis of the effects of mutations in genetically diverse mice and by loss of function experiments. These findings provide evidence that non-redundant genomic locations of different AP-1 family members in macrophages largely result from collaborative interactions with diverse, locus-specific ensembles of transcription factors and suggest a general mechanism for encoding functional specificities of their common recognition motif.

Functional dissection of the ARGONAUTE7 promoter

ARGONAUTES are the central effector proteins of RNA silencing which bind target transcripts in a small RNA-guided manner. Arabidopsis thaliana has 10 ARGONAUTE (AGO) genes, with specialized roles in RNA-directed DNA methylation, post-transcriptional gene silencing, and antiviral defense. To better understand specialization among AGO genes at the level of transcriptional regulation we tested a library of 1497 transcription factors for binding to the promoters of AGO1,AGO10, and AGO7 using yeast 1-hybrid assays. A ranked list of candidate DNA-binding TFs revealed binding of the AGO7 promoter by a number of proteins in two families: the miR156-regulated SPL family and the miR319-regulated TCP family, both of which have roles in developmental timing and leaf morphology. Possible functions for SPL and TCP binding are unclear: we showed that these binding sites are not required for the polar expression pattern of AGO7, nor for the function of AGO7 in leaf shape. Normal AGO7 transcription levels and function appear to depend instead on an adjacent 124-bp region. Progress in understanding the structure of this promoter may aid efforts to understand how the conserved AGO7-triggered TAS3 pathway functions in timing and polarity.

RAEdb: a database of enhancers identified by high-throughput reporter assays

High-throughput reporter assays have been recently developed to directly and quantitatively assess enhancer activity for thousands of regulatory elements. However, there is still no database to collect these enhancers. We developed RAEdb, the first database to collect enhancers identified by high-throughput reporter assays. RAEdb includes 538 320 enhancers derived from eight studies, most of which were from six human cell lines. An activity score was assigned to each enhancer based on reporter assays. Based on these enhancers, 7658 epromoters (promoters with enhancer activity) were identified and stored in the database. RAEdb provides two ways of searches: the first is to search studies by species and cell line; the other is to search enhancers or epromoters by position, activity score, sequence and gene. RAEdb also provides a genome browser to query, visualize and compare enhancers. All data in RAEdb is freely available for download.

Buffering and Amplifying Transcriptional Noise During Cell Fate Specification

The molecular processes that drive gene transcription are inherently noisy. This noise often manifests in the form of transcriptional bursts, producing fluctuations in gene activity over time. During cell fate specification, this noise is often buffered to ensure reproducible developmental outcomes. However, sometimes noise is utilized as a “bet-hedging” mechanism to diversify functional roles across a population of cells. Studies of bacteria, yeast, and cultured cells have provided insights into the nature and roles of noise in transcription, yet we are only beginning to understand the mechanisms by which noise influences the development of multicellular organisms. Here we discuss the sources of transcriptional noise and the mechanisms that either buffer noise to drive reproducible fate choices or amplify noise to randomly specify fates.

From pattern to process: studies at the interface of gene regulatory networks, morphogenesis, and evolution

The development of anatomical structures is complex, beginning with patterning of gene expression by multiple gene regulatory networks (GRNs). These networks ultimately regulate the activity of effector molecules, which in turn alter cellular behavior during development. Together these processes biomechanically produce the three-dimensional shape that the anatomical structure adopts over time. However, the interfaces between these processes are often overlooked and also include counter-intuitive feedback mechanisms. In this review, we examine each step in this extraordinarily complex process and explore how evolutionary developmental biology model systems, such as butterfly scales, vertebrate teeth, and the Drosophila dorsal appendage offer a complementary approach to expose the multifactorial integration of genetics and morphogenesis from an alternative perspective.

Enhancer and super-enhancer: Positive regulators in gene transcription

Enhancer is a positive regulator for spatiotemporal development in eukaryotes. As a cluster, super-enhancer is closely related to cell identity- and fate-determined processes. Both of them function tightly depending on their targeted transcription factors, cofactors, and genes through distal genomic interactions. They have been recognized as critical components and played positive roles in transcriptional regulatory network or factory. Recent advances of next-generation sequencing have dramatically expanded our ability and knowledge to interrogate the molecular mechanism of enhancer and super-enhancer for transcription. Here, we review the history, importance, advances and challenges on enhancer and super-enhancer field. This will benefit our understanding of their function mechanism for transcription underlying precise gene expression.

The genome-wide transcriptional regulatory landscape of ecdysone in the silkworm

BACKGROUND

The silkworm, Bombyx mori, a typical representative of metamorphic insects, is of great agricultural and economic importance. The steroid hormone ecdysone (20-hydroxyecdysone, 20E) is the central regulator of insect developmental transitions, and its nuclear receptors are crucial for numerous biological processes, including reproduction, metabolism, and immunity. However, genome-wide DNA regulatory elements and the ecdysone receptor (EcR) that control these programs of gene expression are not well defined.

RESULTS

In this study, we investigated the alterations in three types of histone modification in silkworm embryonic cells treated with 20E by chromatin immunoprecipitation sequencing (ChIP-seq). We identified enhancers using histone modifications and derived genome-wide ecdysone-dependent enhancer activity maps in the silkworm. We found enhancers enriched for monomethylation of histone H3 Lys4 (H3K4me1) that showed dynamic changes in acetylation of histone H3 Lys27 (H3K27ac) after 20E treatment and functioned to regulate the transcription of specific genes. EcR regulated transcription by binding not only to proximal promoters but also to the distal enhancers of target genes. Moreover, only 52.65% EcR peaks contained ecdysone response element (EcRE) motif, suggesting that EcR regulates the expression of target genes not only by binding directly to EcRE, but also by binding with other transcription factor.

CONCLUSIONS

Our findings provide novel insights into the complex regulatory landscape of hormone-responsive cell activity and a basis for understanding the complex transcriptional regulatory processes of ecdysone.

cis-regulatory architecture of a short-range EGFR organizing center in the Drosophila melanogaster leg

We characterized the establishment of an Epidermal Growth Factor Receptor (EGFR) organizing center (EOC) during leg development in Drosophila melanogaster. Initial EGFR activation occurs in the center of leg discs by expression of the EGFR ligand Vn and the EGFR ligand-processing protease Rho, each through single enhancers, vnE and rhoE, that integrate inputs from Wg, Dpp, Dll and Sp1. Deletion of vnE and rhoE eliminates vn and rho expression in the center of the leg imaginal discs, respectively. Animals with deletions of both vnE and rhoE (but not individually) show distal but not medial leg truncations, suggesting that the distal source of EGFR ligands acts at short-range to only specify distal-most fates, and that multiple additional ‘ring’ enhancers are responsible for medial fates. Further, based on the cis-regulatory logic of vnE and rhoE we identified many additional leg enhancers, suggesting that this logic is broadly used by many genes during Drosophila limb development.

The EGFR signaling pathway plays a major role in innumerable developmental processes in all animals and its deregulation leads to different types of cancer, as well as many other developmental diseases in humans. Here we explored the integration of inputs from the Wnt- and TGF-beta signaling pathways and the leg-specifying transcription factors Distal-less and Sp1 at enhancer elements of EGFR ligands. These enhancers trigger a specific EGFR-dependent developmental output in the fly leg that is limited to specifying distal-most fates. Our findings suggest that activation of the EGFR pathway during fly leg development occurs through the activation of multiple EGFR ligand enhancers that are active at different positions along the proximo-distal axis. Similar enhancer elements are likely to control EGFR activation in humans as well. Such DNA elements might be ‘hot spots’ that cause formation of EGFR-dependent tumors if mutations in them occur. Thus, understanding the molecular characteristics of such DNA elements could facilitate the detection and treatment of cancer.

Mutant p53 regulates enhancer-associated H3K4 monomethylation through interactions with the methyltransferase MLL4

Monomethylation of histone H3 lysine 4 (H3K4me1) is enriched at enhancers that are primed for activation and the levels of this histone mark are frequently altered in various human cancers. Yet, how alterations in H3K4me1 are established and the consequences of these epigenetic changes in tumorigenesis are not well understood. Using ChIP-Seq in human colon cancer cells, we demonstrate that mutant p53 depletion results in decreased H3K4me1 levels at active enhancers that reveal a striking colocalization of mutant p53 and the H3K4 monomethyltransferase MLL4 following chronic tumor necrosis factor alpha (TNFα) signaling. We further reveal that mutant p53 forms physiological associations and direct interactions with MLL4 and promotes the enhancer binding of MLL4, which is required for TNFα-inducible H3K4me1 and histone H3 lysine 27 acetylation (H3K27ac) levels, enhancer-derived transcript (eRNA) synthesis, and mutant p53-dependent target gene activation. Complementary in vitro studies with recombinant chromatin and purified proteins demonstrate that binding of the MLL3/4 complex and H3K4me1 deposition is enhanced by mutant p53 and p300-mediated acetylation, which in turn reflects a MLL3/4-dependent enhancement of mutant p53 and p300-dependent transcriptional activation. Collectively, our findings establish a mechanism in which mutant p53 cooperates with MLL4 to regulate aberrant enhancer activity and tumor-promoting gene expression in response to chronic immune signaling.

Functional CRISPR screen identifies AP1-associated enhancer regulating FOXF1 to modulate oncogene-induced senescence

BACKGROUND

Functional characterization of non-coding elements in the human genome is a major genomic challenge and the maturation of genome-editing technologies is revolutionizing our ability to achieve this task. Oncogene-induced senescence, a cellular state of irreversible proliferation arrest that is enforced following excessive oncogenic activity, is a major barrier against cancer transformation; therefore, bypassing oncogene-induced senescence is a critical step in tumorigenesis. Here, we aim at further identification of enhancer elements that are required for the establishment of this state.

RESULTS

We first apply genome-wide profiling of enhancer-RNAs (eRNAs) to systematically identify enhancers that are activated upon oncogenic stress. DNA motif analysis of these enhancers indicates AP-1 as a major regulator of the transcriptional program induced by oncogene-induced senescence. We thus constructed a CRISPR-Cas9 sgRNA library designed to target senescence-induced enhancers that are putatively regulated by AP-1 and used it in a functional screen. We identify a critical enhancer that we name EnhAP1-OIS1 and validate that mutating the AP-1 binding site within this element results in oncogene-induced senescence bypass. Furthermore, we identify FOXF1 as the gene regulated by this enhancer and demonstrate that FOXF1 mediates EnhAP1-OIS1 effect on the senescence phenotype.

CONCLUSIONS

Our study elucidates a novel cascade mediated by AP-1 and FOXF1 that regulates oncogene-induced senescence and further demonstrates the power of CRISPR-based functional genomic screens in deciphering the function of non-coding regulatory elements in the genome.

Functional characterization of enhancer evolution in the primate lineage

BACKGROUND

Enhancers play an important role in morphological evolution and speciation by controlling the spatiotemporal expression of genes. Previous efforts to understand the evolution of enhancers in primates have typically studied many enhancers at low resolution, or single enhancers at high resolution. Although comparative genomic studies reveal large-scale turnover of enhancers, a specific understanding of the molecular steps by which mammalian or primate enhancers evolve remains elusive.

RESULTS

We identified candidate hominoid-specific liver enhancers from H3K27ac ChIP-seq data. After locating orthologs in 11 primates spanning around 40 million years, we synthesized all orthologs as well as computational reconstructions of 9 ancestral sequences for 348 active tiles of 233 putative enhancers. We concurrently tested all sequences for regulatory activity with STARR-seq in HepG2 cells. We observe groups of enhancer tiles with coherent trajectories, most of which can be potentially explained by a single gain or loss-of-activity event per tile. We quantify the correlation between the number of mutations along a branch and the magnitude of change in functional activity. Finally, we identify 84 mutations that correlate with functional changes; these are enriched for cytosine deamination events within CpGs.

CONCLUSIONS

We characterized the evolutionary-functional trajectories of hundreds of liver enhancers throughout the primate phylogeny. We observe subsets of regulatory sequences that appear to have gained or lost activity. We use these data to quantify the relationship between sequence and functional divergence, and to identify CpG deamination as a potentially important force in driving changes in enhancer activity during primate evolution.

Dynamic interplay between enhancer-promoter topology and gene activity

A long-standing question in gene regulation is how remote enhancers communicate with their target promoters, and specifically how chromatin topology dynamically relates to gene activation. Here, we combine genome editing and multi-color live imaging to simultaneously visualize physical enhancer-promoter interaction and transcription at the single-cell level in Drosophila embryos. By examining transcriptional activation of a reporter by the endogenous even-skipped enhancers, which are located 150 kb away, we identify three distinct topological conformation states and measure their transition kinetics. We show that sustained proximity of the enhancer to its target is required for activation. Transcription in turn affects the three-dimensional topology as it enhances the temporal stability of the proximal conformation and is associated with further spatial compaction. Furthermore, the facilitated long-range activation results in transcriptional competition at the locus, causing corresponding developmental defects. Our approach offers quantitative insight into the spatial and temporal determinants of long-range gene regulation and their implications for cellular fates.