Drosophila Enhancer Research Articles

Cell-type-directed design of synthetic enhancers.

Transcriptional enhancers act as docking stations for combinations of transcription factors and thereby regulate spatiotemporal activation of their target genes1. It has been a long-standing goal in the field to decode the regulatory logic of an enhancer and to understand the details of how spatiotemporal gene expression is encoded in an enhancer sequence. Here we show that deep learning models2-6, can be used to efficiently design synthetic, cell-type-specific enhancers, starting from random sequences, and that this optimization process allows detailed tracing of enhancer features at single-nucleotide resolution. We evaluate the function of fully synthetic enhancers to specifically target Kenyon cells or glial cells in the fruit fly brain using transgenic animals. We further exploit enhancer design to create 'dual-code' enhancers that target two cell types and minimal enhancers smaller than 50 base pairs that are fully functional. By examining the state space searches towards local optima, we characterize enhancer codes through the strength, combination and arrangement of transcription factor activator and transcription factor repressor motifs. Finally, we apply the same strategies to successfully design human enhancers, which adhere to enhancer rules similar to those of Drosophila enhancers. Enhancer design guided by deep learning leads to better understanding of how enhancers work and shows that their code can be exploited to manipulate cell states.

An explainable artificial intelligence approach for decoding the enhancer histone modifications code and identification of novel enhancers in Drosophila

BackgroundEnhancers are non-coding regions of the genome that control the activity of target genes. Recent efforts to identify active enhancers experimentally and in silico have proven effective. While these tools can predict the locations of enhancers with a high degree of accuracy, the mechanisms underpinning the activity of enhancers are often unclear.ResultsUsing machine learning (ML) and a rule-based explainable artificial intelligence (XAI) model, we demonstrate that we can predict the location of known enhancers in Drosophila with a high degree of accuracy. Most importantly, we use the rules of the XAI model to provide insight into the underlying combinatorial histone modifications code of enhancers. In addition, we identified a large set of putative enhancers that display the same epigenetic signature as enhancers identified experimentally. These putative enhancers are enriched in nascent transcription, divergent transcription and have 3D contacts with promoters of transcribed genes. However, they display only intermediary enrichment of mediator and cohesin complexes compared to previously characterised active enhancers. We also found that 10–15% of the predicted enhancers display similar characteristics to super enhancers observed in other species.ConclusionsHere, we applied an explainable AI model to predict enhancers with high accuracy. Most importantly, we identified that different combinations of epigenetic marks characterise different groups of enhancers. Finally, we discovered a large set of putative enhancers which display similar characteristics with previously characterised active enhancers.

Deciphering enhancer sequence using thermodynamics-based models and convolutional neural networks.

Deciphering the sequence-function relationship encoded in enhancers holds the key to interpreting non-coding variants and understanding mechanisms of transcriptomic variation. Several quantitative models exist for predicting enhancer function and underlying mechanisms; however, there has been no systematic comparison of these models characterizing their relative strengths and shortcomings. Here, we interrogated a rich data set of neuroectodermal enhancers in Drosophila, representing cis- and trans- sources of expression variation, with a suite of biophysical and machine learning models. We performed rigorous comparisons of thermodynamics-based models implementing different mechanisms of activation, repression and cooperativity. Moreover, we developed a convolutional neural network (CNN) model, called CoNSEPT, that learns enhancer ‘grammar’ in an unbiased manner. CoNSEPT is the first general-purpose CNN tool for predicting enhancer function in varying conditions, such as different cell types and experimental conditions, and we show that such complex models can suggest interpretable mechanisms. We found model-based evidence for mechanisms previously established for the studied system, including cooperative activation and short-range repression. The data also favored one hypothesized activation mechanism over another and suggested an intriguing role for a direct, distance-independent repression mechanism. Our modeling shows that while fundamentally different models can yield similar fits to data, they vary in their utility for mechanistic inference. CoNSEPT is freely available at: https://github.com/PayamDiba/CoNSEPT.

Evolution of hes gene family in vertebrates: the hes5 cluster genes have specifically increased in frogs

Backgroundhes genes are chordate homologs of Drosophila genes, hairy and enhancer of split, which encode a basic helix-loop-helix (bHLH) transcriptional repressor with a WRPW motif. Various developmental functions of hes genes, including early embryogenesis and neurogenesis, have been elucidated in vertebrates. However, their orthologous relationships remain unclear partly because of less conservation of relatively short amino acid sequences, the fact that the genome was not analyzed as it is today, and species-specific genome duplication. This results in complicated gene names in vertebrates, which are not consistent in orthologs. We previously revealed that Xenopus frogs have two clusters of hes5, named “the hes5.1 cluster” and “the hes5.3 cluster”, but the origin and the conservation have not yet been revealed.ResultsHere, we elucidated the orthologous and paralogous relationships of all hes genes of human, mouse, chicken, gecko, zebrafish, medaka, coelacanth, spotted gar, elephant shark and three species of frogs, Xenopus tropicalis (X. tropicalis), X. laevis, Nanorana parkeri, by phylogenetic and synteny analyses. Any duplicated hes5 were not found in mammals, whereas hes5 clusters in teleost were conserved although not as many genes as the three frog species. In addition, hes5 cluster-like structure was found in the elephant shark genome, but not found in cyclostomata.ConclusionThese data suggest that the hes5 cluster existed in the gnathostome ancestor but became a single gene in mammals. The number of hes5 cluster genes were specifically large in frogs.

A small UTX stabilization domain of Trr is conserved within mammalian MLL3-4/COMPASS and is sufficient to rescue loss of viability in null animals.

Catalytic-inactivating mutations within the Drosophila enhancer H3K4 mono-methyltransferase Trr and its mammalian homologs, MLL3/4, cause only minor changes in gene expression compared with whole-gene deletions for these COMPASS members. To identify essential histone methyltransferase-independent functions of Trr, we screened to identify a minimal Trr domain sufficient to rescue Trr-null lethality and demonstrate that this domain binds and stabilizes Utx in vivo. Using the homologous MLL3/MLL4 human sequences, we mapped a short ∼80-amino-acid UTX stabilization domain (USD) that promotes UTX stability in the absence of the rest of MLL3/4. Nuclear UTX stability is enhanced when the USD is fused with the MLL4 HMG-box. Thus, COMPASS-dependent UTX stabilization is an essential noncatalytic function of Trr/MLL3/MLL4, suggesting that stabilizing UTX could be a therapeutic strategy for cancers with MLL3/4 loss-of-function mutations.

Memory of Divisional History Directs the Continuous Process of Primitive Hematopoietic Lineage Commitment.

SummaryHematopoietic stem cells (HSCs) exist in a dormant state and progressively lose regenerative potency as they undergo successive divisions. Why this functional decline occurs and how this information is encoded is unclear. To better understand how this information is stored, we performed RNA sequencing on HSC populations differing only in their divisional history. Comparative analysis revealed that genes upregulated with divisions are enriched for lineage genes and regulated by cell-cycle-associated transcription factors, suggesting that proliferation itself drives lineage priming. Downregulated genes are, however, associated with an HSC signature and targeted by the Polycomb Repressive Complex 2 (PRC2). The PRC2 catalytic subunits Ezh1 and Ezh2 promote and suppress the HSC state, respectively, and successive divisions cause a switch from Ezh1 to Ezh2 dominance. We propose that cell divisions drive lineage priming and Ezh2 accumulation, which represses HSC signature genes to consolidate information on divisional history into memory.

Differential efficacy of genetically swapping GAL4

Several large or mid-scale collections of Drosophila enhancer traps have been recently created to allow for genetic swapping of GAL4 coding sequences to versatile transcription activators or suppressors such as LexA, QF, split-GAL4 (GAL4-AD and GAL4-DBD), GAL80 and QS. Yet a systematic analysis of the feasibility and reproducibility of these tools is lacking. Here we focused on InSITE GAL4 drivers that specifically label different subpopulations of olfactory neurons, particularly local interneurons (LNs), and genetically swapped the GAL4 domain for LexA, GAL80 or QF at the same locus. We found that the major utility-limiting factor for these genetic swaps is that many do not fully reproduce the original GAL4 expression patterns. Different donors exhibit distinct efficacies for reproducing original GAL4 expression patterns. The successfully swapped lines reported here will serve as valuable reagents and expand the genetic toolkits of Drosophila olfactory circuit research.

Memory enhancement by ferulic acid ester across species.

Cognitive impairments can be devastating for quality of life, and thus, preventing or counteracting them is of great value. To this end, the present study exploits the potential of the plant Rhodiola rosea and identifies the constituent ferulic acid eicosyl ester [icosyl-(2E)-3-(4-hydroxy-3-methoxyphenyl)-prop-2-enoate (FAE-20)] as a memory enhancer. We show that food supplementation with dried root material from R. rosea dose-dependently improves odor-taste reward associative memory scores in larval Drosophila and prevents the age-related decline of this appetitive memory in adult flies. Task-relevant sensorimotor faculties remain unaltered. From a parallel approach, a list of candidate compounds has been derived, including R. rosea-derived FAE-20. Here, we show that both R. rosea-derived FAE-20 and synthetic FAE-20 are effective as memory enhancers in larval Drosophila. Synthetic FAE-20 also partially compensates for age-related memory decline in adult flies, as well as genetically induced early-onset loss of memory function in young flies. Furthermore, it increases excitability in mouse hippocampal CA1 neurons, leads to more stable context-shock aversive associative memory in young adult (3-month-old) mice, and increases memory scores in old (>2-year-old) mice. Given these effects, and given the utility of R. rosea-the plant from which we discovered FAE-20-as a memory enhancer, these results may hold potential for clinical applications.

Enhancer transcription: what, where, when, and why?

Following the discovery of widespread enhancer transcription, enhancers and promoters have been found to be far more similar than previously thought. In this issue of Genes & Development, two studies (Henriques and colleagues [pp. 26-41] and Mikhaylichenko and colleagues [pp. 42-57]) shine new light on the transcriptional nature of promoters and enhancers in Drosophila Together, these studies support recent work in mammalian cells that indicates that most active enhancers drive local transcription using factors and mechanisms similar to those of promoters. Intriguingly, enhancer transcription is shown to be coordinated by SPT5- and P-TEFb-mediated pause-release, but the pause half-life is shorter, and termination is more rapid at enhancers than at promoters. Moreover, bidirectional transcription from promoters is associated with enhancer activity, lending further credence to models in which regulatory elements exist along a spectrum of promoter-ness and enhancer-ness. We propose a general unified model to explain possible functions of transcription at enhancers.

Widespread transcriptional pausing and elongation control at enhancers.

Regulation by gene-distal enhancers is critical for cell type-specific and condition-specific patterns of gene expression. Thus, to understand the basis of gene activity in a given cell type or tissue, we must identify the precise locations of enhancers and functionally characterize their behaviors. Here, we demonstrate that transcription is a nearly universal feature of enhancers in Drosophila and mammalian cells and that nascent RNA sequencing strategies are optimal for identification of both enhancers and superenhancers. We dissect the mechanisms governing enhancer transcription and discover remarkable similarities to transcription at protein-coding genes. We show that RNA polymerase II (RNAPII) undergoes regulated pausing and release at enhancers. However, as compared with mRNA genes, RNAPII at enhancers is less stable and more prone to early termination. Furthermore, we found that the level of histone H3 Lys4 (H3K4) methylation at enhancers corresponds to transcriptional activity such that highly active enhancers display H3K4 trimethylation rather than the H3K4 monomethylation considered a hallmark of enhancers. Finally, our work provides insights into the unique characteristics of superenhancers, which stimulate high-level gene expression through rapid pause release; interestingly, this property renders associated genes resistant to the loss of factors that stabilize paused RNAPII.

Short DNA sequence patterns accurately identify broadly active human enhancers

BackgroundEnhancers are DNA regulatory elements that influence gene expression. There is substantial diversity in enhancers’ activity patterns: some enhancers drive expression in a single cellular context, while others are active across many. Sequence characteristics, such as transcription factor (TF) binding motifs, influence the activity patterns of regulatory sequences; however, the regulatory logic through which specific sequences drive enhancer activity patterns is poorly understood. Recent analysis of Drosophila enhancers suggested that short dinucleotide repeat motifs (DRMs) are general enhancer sequence features that drive broad regulatory activity. However, it is not known whether the regulatory role of DRMs is conserved across species.ResultsWe performed a comprehensive analysis of the relationship between short DNA sequence patterns, including DRMs, and human enhancer activity in 38,538 enhancers across 411 different contexts. In a machine-learning framework, the occurrence patterns of short sequence motifs accurately predicted broadly active human enhancers. However, DRMs alone were weakly predictive of broad enhancer activity in humans and showed different enrichment patterns than in Drosophila. In general, GC-rich sequence motifs were significantly associated with broad enhancer activity, and consistent with this enrichment, broadly active human TFs recognize GC-rich motifs.ConclusionsOur results reveal the importance of specific sequence motifs in broadly active human enhancers, demonstrate the lack of evolutionary conservation of the role of DRMs, and provide a computational framework for investigating the logic of enhancer sequences.

MiR-203 Interplays with Polycomb Repressive Complexes to Regulate the Proliferation of Neural Stem/Progenitor Cells.

SummaryThe polycomb repressive complexes 1 (PRC1) and 2 (PRC2) are two distinct polycomb group (PcG) proteins that maintain the stable silencing of specific sets of genes through chromatin modifications. Although the PRC2 component EZH2 has been known as an epigenetic regulator in promoting the proliferation of neural stem/progenitor cells (NSPCs), the regulatory network that controls this process remains largely unknown. Here we show that miR-203 is repressed by EZH2 in both embryonic and adult NSPCs. MiR-203 negatively regulates the proliferation of NSPCs. One of PRC1 components, Bmi1, is a downstream target of miR-203 in NSPCs. Conditional knockout of Ezh2 results in decreased proliferation ability of both embryonic and adult NSPCs. Meanwhile, ectopic overexpression of BMI1 rescues the proliferation defects exhibited by miR-203 overexpression or EZH2 deficiency in NSPCs. Therefore, this study provides evidence for coordinated function of the EZH2-miR-203-BMI1 regulatory axis that regulates the proliferation of NSPCs.

Drosophila Enhancer of Rudimentary Homolog, ERH, Is a Binding Partner of RPS3, RPL19, and DDIT4, Suggesting a Mechanism for the Nuclear Localization of ERH.

The protein enhancer of rudimentary homolog, ERH, is a small, highly conserved protein that has been found in animals, plants, and protists. Genetic and biochemical interactions have implicated ERH in the regulation of pyrimidine biosynthesis, DNA replication, transcription, mRNA splicing, cellular proliferation, tumorigenesis, and the Notch signaling pathway. In vertebrates and insects, ERH is nuclearly localized; however, an examination of the ERH amino-acid sequence does not reveal any nuclear localization signals. In this paper we show that the first 24 amino acids contain sequences necessary and sufficient for nuclear localization. Through yeast two-hybrid screens, three new binding partners of ERH, RPS3, RPL19, and DDIT4, were identified. RPS3 was isolated from both human and Drosophila screens. These interactions suggest functions of ERH in cell growth, cancer, and DNA repair. The ERH sequences necessary for the interactions between ERH and RPS3 and RPL19 are mapped onto the same 24-amino-acid region in ERH which are necessary for nuclear localization, suggesting that ERH is localizing to the nucleus through binding to one of its DNA-binding partners, such as RPS3 or RPL19.

The Capacity to Act in Trans Varies Among Drosophila Enhancers.

The interphase nucleus is organized such that genomic segments interact in cis, on the same chromosome, and in trans, between different chromosomes. In Drosophila and other Dipterans, extensive interactions are observed between homologous chromosomes, which can permit enhancers and promoters to communicate in trans Enhancer action in trans has been observed for a handful of genes in Drosophila, but it is as yet unclear whether this is a general property of all enhancers or specific to a few. Here, we test a collection of well-characterized enhancers for the capacity to act in trans Specifically, we tested 18 enhancers that are active in either the eye or wing disc of third instar Drosophila larvae and, using two different assays, found evidence that each enhancer can act in trans However, the degree to which trans-action was supported varied greatly between enhancers. Quantitative analysis of enhancer activity supports a model wherein an enhancer's strength of transcriptional activation is a major determinant of its ability to act in trans, but that additional factors may also contribute to an enhancer's trans-activity. In sum, our data suggest that a capacity to activate a promoter on a paired chromosome is common among Drosophila enhancers.

The appeasement of Doug: a synthetic approach to enhancer biology.

Genetic approaches have been instrumental in dissecting developmental enhancers by characterizing their transcription factor binding sites. Though some enhancers have been well-studied in this regard, we cannot currently build developmental enhancers from scratch. Reconstitution experiments can provide important complementary tests of our understanding of enhancer function, but these experiments are exceedingly rare in the literature, possibly due to the difficulty of publishing negative results. In this perspective, we argue that the time is right for a synthetic approach to enhancer biology. Focusing primarily on Drosophila enhancers as examples, we review classic and modern methods for dissecting enhancer function as well as computational tools for enhancer design. We include our own negative results from attempts to reconstitute the stripe 2 enhancer from the even-skipped locus and discuss possible ways forward. We believe that with a communal effort in open data sharing, we can make substantial progress toward a complete understanding of enhancer function.

The effect of transcription on enhancer activity in Drosophila melanogaster

In higher eukaryotes, the level of gene transcription is under the control of DNA regulatory elements, such as promoter, from which transcription is initiated with the participation of RNA polymerase II and general transcription factors, as well as the enhancer, which increase the rate of transcription with the involvement of activator proteins and cofactors. It was demonstrated that enhancers are often located in the transcribed regions of the genome. We showed earlier that transcription negatively affected the activity of enhancers in Drosophila in model transgenic systems. In this study, we tested the effect of the distance between the leading promoter, enhancer, and target promoter on the inhibitory effect of transcriptions of different strengths. It was demonstrated that the negative effect of transcription remained, but weakened with increased distance between the leading promoter and enhancer and with decreased distance between the enhancer and target promoter. Thus, transcription can modulate the activity of enhancers by controlling its maximum level.

Genome-wide identification of hypoxia-induced enhancer regions.

Here we present a genome-wide method for de novo identification of enhancer regions. This approach enables massively parallel empirical investigation of DNA sequences that mediate transcriptional activation and provides a platform for discovery of regulatory modules capable of driving context-specific gene expression. The method links fragmented genomic DNA to the transcription of randomer molecule identifiers and measures the functional enhancer activity of the library by massively parallel sequencing. We transfected a Drosophila melanogaster library into S2 cells in normoxia and hypoxia, and assayed 4,599,881 genomic DNA fragments in parallel. The locations of the enhancer regions strongly correlate with genes up-regulated after hypoxia and previously described enhancers. Novel enhancer regions were identified and integrated with RNAseq data and transcription factor motifs to describe the hypoxic response on a genome-wide basis as a complex regulatory network involving multiple stress-response pathways. This work provides a novel method for high-throughput assay of enhancer activity and the genome-scale identification of 31 hypoxia-activated enhancers in Drosophila.

A Unified Model Describing The Architecture And Creation Of Promoters And Enhancers

Promoters and enhancers have been conventionally regarded as structurally and functionally distinct genomic entities. Nonetheless, in mammals, they share many features, including divergent transcription and common transcription factors. Here, we examine the architecture of transcription initiation through comprehensive mapping of transcription start sites (TSSs) in human cell lines. Using a nuclear run‐on protocol called GRO‐cap, which captures TSSs for both stable and unstable transcripts, we conduct detailed comparisons of thousands of promoters and enhancers. These analyses reveal a common architecture of initiation, including tightly spaced (110 bp) divergent initiation, similar frequencies of core‐promoter sequence elements, highly positioned flanking nucleosomes, and two modes of transcription factor binding. Post‐initiation transcript stability provides a more fundamental distinction between promoters and enhancers than patterns of histone modifications, transcription factors or co‐activators. These results support a unified model of transcription initiation and architecture at promoters and enhancers. We have now investigated how these patterns of divergent transcription are generated on enhancers in both Drosophila and mammalian cells. These new results will be presented and the implications of these findings on genome function and evolution will be discussed.

The Fusarium graminearum Histone H3 K27 Methyltransferase KMT6 Regulates Development and Expression of Secondary Metabolite Gene Clusters

The cereal pathogen Fusarium graminearum produces secondary metabolites toxic to humans and animals, yet coordinated transcriptional regulation of gene clusters remains largely a mystery. By chromatin immunoprecipitation and high-throughput DNA sequencing (ChIP-seq) we found that regions with secondary metabolite clusters are enriched for trimethylated histone H3 lysine 27 (H3K27me3), a histone modification associated with gene silencing. H3K27me3 was found predominantly in regions that lack synteny with other Fusarium species, generally subtelomeric regions. Di- or trimethylated H3K4 (H3K4me2/3), two modifications associated with gene activity, and H3K27me3 are predominantly found in mutually exclusive regions of the genome. To find functions for H3K27me3, we deleted the gene for the putative H3K27 methyltransferase, KMT6, a homolog of Drosophila Enhancer of zeste, E(z). The kmt6 mutant lacks H3K27me3, as shown by western blot and ChIP-seq, displays growth defects, is sterile, and constitutively expresses genes for mycotoxins, pigments and other secondary metabolites. Transcriptome analyses showed that 75% of 4,449 silent genes are enriched for H3K27me3. A subset of genes that were enriched for H3K27me3 in WT gained H3K4me2/3 in kmt6. A largely overlapping set of genes showed increased expression in kmt6. Almost 95% of the remaining 2,720 annotated silent genes showed no enrichment for either H3K27me3 or H3K4me2/3 in kmt6. In these cases mere absence of H3K27me3 was insufficient for expression, which suggests that additional changes are required to activate genes. Taken together, we show that absence of H3K27me3 allowed expression of an additional 14% of the genome, resulting in derepression of genes predominantly involved in secondary metabolite pathways and other species-specific functions, including putative secreted pathogenicity factors. Results from this study provide the framework for novel targeted strategies to control the “cryptic genome”, specifically secondary metabolite expression.

The DrosophilaEnhancer of splitGene Complex: Architecture and Coordinate Regulation by Notch, Cohesin, and Polycomb Group Proteins

The cohesin protein complex functionally interacts with Polycomb group (PcG) silencing proteins to control expression of several key developmental genes, such as the Drosophila Enhancer of split gene complex [E(spl)-C]. The E(spl)-C contains 12 genes that inhibit neural development. In a cell line derived from the central nervous system, cohesin and the PRC1 PcG protein complex bind and repress E (spl)-C transcription, but the repression mechanisms are unknown. The genes in the E(spl)-C are directly activated by the Notch receptor. Here we show that depletion of cohesin or PRC1 increases binding of the Notch intracellular fragment to genes in the E(spl)-C, correlating with increased transcription. The increased transcription likely reflects both direct effects of cohesin and PRC1 on RNA polymerase activity at the E(spl)-C, and increased expression of Notch ligands. By chromosome conformation capture we find that the E(spl)-C is organized into a self-interactive architectural domain that is co-extensive with the region that binds cohesin and PcG complexes. The self-interactive architecture is formed independently of cohesin or PcG proteins. We posit that the E(spl)-C architecture dictates where cohesin and PcG complexes bind and act when they are recruited by as yet unidentified factors, thereby controlling the E(spl)-C as a coordinated domain.