Transcriptional Silencers in Drosophila Serve a Dual Role as Transcriptional Enhancers in Alternate Cellular Contexts

CREATE: cell-type-specific cis-regulatory element identification via discrete embedding

Cis-regulatory elements (CREs), including enhancers, silencers, promoters and insulators, play pivotal roles in orchestrating gene regulatory mechanisms that drive complex biological traits. However, current approaches for CRE identification are predominantly sequence-based and typically focus on individual CRE types, limiting insights into their cell-type-specific functions and regulatory dynamics. Here, we present CREATE, a multimodal deep learning framework based on Vector Quantized Variational AutoEncoder, tailored for comprehensive CRE identification and characterization. CREATE integrates genomic sequences, chromatin accessibility, and chromatin interaction data to generate discrete CRE embeddings, enabling accurate multi-class classification and robust characterization of CREs. CREATE excels in identifying cell-type-specific CREs, and provides quantitative and interpretable insights into CRE-specific features, uncovering the underlying regulatory codes. By facilitating large-scale prediction of CREs in specific cell types, CREATE enhances the recognition of disease- or phenotype-associated biological variabilities of CREs, thus advancing our understanding of gene regulatory landscapes and their roles in health and disease.

Heritable maintenance of chromatin modifications confers transcriptional memory of interferon-γ signaling

Interferon-γ (IFNγ) transiently activates genes related to inflammation and innate immunity. A subset of targets retain a mitotically heritable memory of prior IFNγ exposure, resulting in hyperactivation upon re-exposure through poorly understood mechanisms. Here, we discover that the transcriptionally permissive chromatin marks H3K4me1, H3K14ac and H4K16ac are established during IFNγ priming and are selectively maintained on a cluster of guanylate-binding protein (GBP) genes in dividing human cells in the absence of transcription. The histone acetyltransferase KAT7 is required for H3K14ac deposition at GBP genes and for accelerated GBP reactivation upon re-exposure to IFNγ. In naive cells, the GBP cluster is maintained in a low-level repressive chromatin state, marked by H3K27me3, limiting priming through a PRC2-dependent mechanism. Unexpectedly, IFNγ priming results in transient accumulation of this repressive mark despite active gene expression. However, during the memory phase, H3K27 methylation is selectively depleted from primed GBP genes, facilitating hyperactivation. Furthermore, we identified a cis-regulatory element that forms transient, long-range contacts across the GBP cluster and acts as a repressor, curbing hyperactivation of previously IFNγ-primed cells. Our results provide insight into the chromatin basis for the long-term transcriptional memory of IFNγ signaling, which might contribute to enhanced innate immunity.

Defining the regulatory logic of breast cancer using single-cell epigenetic and transcriptome profiling

Annotation of cis-regulatory elements that drive transcriptional dysregulation in cancer cells is critical to understanding tumor biology. Herein, we present matched chromatin accessibility (single-cell assay for transposase-accessible chromatin by sequencing [scATAC-seq]) and transcriptome (single-cell RNA sequencing [scRNA-seq]) profiles at single-cell resolution from human breast tumors and healthy mammary tissues processed immediately following surgical resection. We identify the most likely cell of origin for subtype-specific breast tumors and implement linear mixed-effects modeling to quantify associations between regulatory elements and gene expression in malignant versus normal cells. These data unveil cancer-specific regulatory elements and putative silencer-to-enhancer switching events in cells that lead to the upregulation of clinically relevant oncogenes. In addition, we generate matched scATAC-seq and scRNA-seq profiles for breast cancer cell lines, revealing a conserved oncogenic gene expression program between in vitro and in vivo cells. This work highlights the importance of non-coding regulatory mechanisms that underlie oncogenic processes and the ability of single-cell multi-omics to define the regulatory logic of cancer cells.

Abundant repressor binding sites in human enhancers are associated with the fine-tuning of gene regulation

The regulation of gene expression relies on the coordinated action of transcription factors (TFs) at enhancers, including both activator and repressor TFs. We employed deep learning (DL) to dissect HepG2 enhancers into positive (PAR), negative (NAR), and neutral activity regions. Sharpr-MPRA and STARR-seq highlight the dichotomy impact of NARs and PARs on modulating and catalyzing the activity of enhancers, respectively. Approximately 22% of HepG2 enhancers, termed "repressive impact enhancers" (RIEs), are predominantly populated by NARs and transcriptional repression motifs. Genes flanking RIEs exhibit a stage-specific decline in expression during late development, suggesting RIEs' role in trimming enhancer activities. About 16.7% of human NARs emerge from neutral rhesus macaque DNA. This gain of repressor binding sites in RIEs is associated with a 30% decrease in the average expression of flanking genes in humans compared to rhesus macaque. Our work reveals modulated enhancer activity and adaptable gene regulation through the evolutionary dynamics of TF binding sites.

Active learning of enhancers and silencers in the developing neural retina

Deep learning is a promising strategy for modeling cis-regulatory elements. However, models trained on genomic sequences often fail to explain why the same transcription factor can activate or repress transcription in different contexts. To address this limitation, we developed an active learning approach to train models that distinguish between enhancers and silencers composed of binding sites for the photoreceptor transcription factor cone-rod homeobox (CRX). After training the model on nearly all bound CRX sites from the genome, we coupled synthetic biology with uncertainty sampling to generate additional rounds of informative training data. This allowed us to iteratively train models on data from multiple rounds of massively parallel reporter assays. The ability of the resulting models to discriminate between CRX sites with identical sequence but opposite functions establishes active learning as an effective strategy to train models of regulatory DNA. A record of this paper's transparent peer review process is included in the supplemental information.

Silent silencers: Leave no trace

The regulatory mechanisms of silencers have remained poorly understood. In this issue, Hofbauer et al.1 conduct a genome-wide screen in Drosophila melanogaster and reveal three silencer types that appear to work alone-without the need for combinatorial action, traditional chromatin marks, or open chromatin regions.

A genome-wide screen identifies silencers with distinct chromatin properties and mechanisms of repression

Differential gene transcription enables development and homeostasis in all animals and is regulated by two major classes of distal cis-regulatory DNA elements (CREs): enhancers and silencers. Although enhancers have been thoroughly characterized, the properties and mechanisms of silencers remain largely unknown. By an unbiased genome-wide functional screen in Drosophila melanogaster S2 cells, we discover a class of silencers that bind one of three transcription factors (TFs) and are generally not included in chromatin-defined CRE catalogs as they mostly lack detectable DNA accessibility. The silencer-binding TF CG11247, which we term Saft, safeguards cell fate decisions in vivo and functions via a highly conserved domain we term zinc-finger-associated C-terminal (ZAC) and the corepressor G9a, independently of G9a's H3K9-methyltransferase activity. Overall, our identification of silencers with unexpected properties and mechanisms has important implications for the understanding and future study of repressive CREs, as well as the functional annotation of animal genomes.

Beyond modular enhancers: new questions in cis-regulatory evolution

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

Chromatin structure and 3D architecture define the differential functions of PU.1 regulatory elements in blood cell lineages

The precise spatiotemporal expression of the hematopoietic ETS transcription factor PU.1, a key determinant of hematopoietic cell fates, is tightly regulated at the chromatin level. However, how chromatin signatures are linked to this dynamic expression pattern across different blood cell lineages remains uncharacterized. Here, we performed an in-depth analysis of the relationships between gene expression, chromatin structure, 3D architecture, and trans-acting factors at PU.1 cis-regulatory elements (PCREs). By identifying phylogenetically conserved DNA elements within chromatin-accessible regions in primary human blood lineages, we discovered multiple novel candidate PCREs within the upstream region of the human PU.1 locus. A subset of these elements localizes within an 8-kb-wide cluster exhibiting enhancer features, including open chromatin, demethylated DNA, enriched enhancer histone marks, present enhancer RNAs, and PU.1 occupation, presumably mediating PU.1 autoregulation. Importantly, we revealed the presence of a common 35-kb-wide CTCF-flanked insulated neighborhood that contains the PCRE cluster (PCREC), forming a chromatin territory for lineage-specific and PCRE-mediated chromatin interactions. These include functional PCRE-promoter interactions in myeloid and B cells that are absent in erythroid and T cells. By correlating chromatin structure and 3D architecture with PU.1 expression in various lineages, we were able to attribute enhancer versus silencer functions to individual elements. Our findings provide mechanistic insights into the interplay between dynamic chromatin structure and 3D architecture in the chromatin regulation of PU.1 expression. This study lays crucial groundwork for additional experimental studies that validate and dissect the role of PCREs in epigenetic regulation of normal and malignant hematopoiesis.

Small DNA elements that act as both insulators and silencers in plants

Insulators are cis-regulatory elements that separate transcriptional units, whereas silencers are elements that repress transcription regardless of their position. In plants, these elements remain largely uncharacterized. Here, we use the massively parallel reporter assay Plant STARR-seq with short fragments of eight large insulators to identify more than 100 fragments that block enhancer activity. The short fragments can be combined to generate more powerful insulators that abolish the capacity of the strong viral 35S enhancer to activate the 35S minimal promoter. Unexpectedly, when tested upstream of weak enhancers, these fragments act as silencers and repress transcription. Thus, these elements are capable of both insulating or repressing transcription dependent upon regulatory context. We validate our findings in stable transgenic Arabidopsis, maize, and rice plants. The short elements identified here should be useful building blocks for plant biotechnology efforts.

Understanding variants of unknown significance: the computational frontier

The rapid advancement of sequencing technologies has led to the identification of numerous mutations in cancer genomes, many of which are variants of unknown significance (VUS). Computational models are increasingly being used to predict the functional impact of these mutations, in both coding and noncoding regions. Integration of these models with emerging genomic datasets will refine our understanding of mutation effects and guide clinical decision making. Future advancements in modeling protein interactions and transcriptional regulation will further enhance our ability to interpret VUS. Periodic incorporation of these developments into VUS reclassification practice has the potential to significantly improve personalized cancer care.

Genome-wide Cas9-mediated screening of essential non-coding regulatory elements via libraries of paired single-guide RNAs

The functions of non-coding regulatory elements (NCREs), which constitute a major fraction of the human genome, have not been systematically studied. Here we report a method involving libraries of paired single-guide RNAs targeting both ends of an NCRE as a screening system for the Cas9-mediated deletion of thousands of NCREs genome-wide to study their functions in distinct biological contexts. By using K562 and 293T cell lines and human embryonic stem cells, we show that NCREs can have redundant functions, and that many ultra-conserved elements have silencer activity and play essential roles in cell growth and in cellular responses to drugs (notably, the ultra-conserved element PAX6_Tarzan may be critical for heart development, as removing it from human embryonic stem cells led to defects in cardiomyocyte differentiation). The high-throughput screen, which is compatible with single-cell sequencing, may allow for the identification of druggable NCREs.

KSHV genome harbors both constitutive and lytically induced enhancers

Kaposi's sarcoma-associated herpesvirus (KSHV) belongs to the gamma-herpesvirus family and is a well-known human oncogenic virus. In infected cells, the viral genome of 165 kbp is circular DNA wrapped in chromatin. The tight control of gene expression is critical for latency, the transition into the lytic phase, and the development of viral-associated malignancies. Distal cis-regulatory elements, such as enhancers and silencers, can regulate gene expression in a position- and orientation-independent manner. Open chromatin is another characteristic feature of enhancers. To systematically search for enhancers, we cloned all the open chromatin regions in the KSHV genome downstream of the luciferase gene and tested their enhancer activity in infected and uninfected cells. A silencer was detected upstream of the latency-associated nuclear antigen promoter. Two constitutive enhancers were identified in the K12p-OriLyt-R and ORF29 Intron regions, where ORF29 Intron is a tissue-specific enhancer. The following promoters: OriLyt-L, PANp, ALTp, and the terminal repeats (TRs) acted as lytically induced enhancers. The expression of the replication and transcription activator (RTA), the master regulator of the lytic cycle, was sufficient to induce the activity of lytic enhancers in uninfected cells. We propose that the TRs that span about 24 kbp region serve as a "viral super-enhancer" that integrates the repressive effect of the latency-associated nuclear antigen (LANA) with the activating effect of RTA. Utilizing CRISPR activation and interference techniques, we determined the connections between these enhancers and their regulated genes. The silencer and enhancers described here provide an additional layer to the complex gene regulation of herpesviruses.IMPORTANCEIn this study, we performed a systematic functional assay to identify cis-regulatory elements within the genome of the oncogenic herpesvirus, Kaposi's sarcoma-associated herpesvirus (KSHV). Similar to other herpesviruses, KSHV presents both latent and lytic phases. Therefore, our assays were performed in uninfected cells, during latent infection, and under lytic conditions. We identified two constitutive enhancers, one of which seems to be a tissue-specific enhancer. In addition, four lytically induced enhancers, which are all responsive to the replication and transcription activator (RTA), were identified. Furthermore, a silencer was identified between the major latency promoter and the lytic gene locus. Utilizing CRISPR activation and interference techniques, we determined the connections between these enhancers and their regulated genes. The terminal repeats, spanning a region of about 24 kbp, seem like a "viral super-enhancer" that integrates the repressive effect of the latency-associated nuclear antigen (LANA) with the activating effect of RTA to regulate latency to lytic transition.

Evolution of a Human-Specific De Novo Open Reading Frame and Its Linked Transcriptional Silencer

In the human genome, two short open reading frames (ORFs) separated by a transcriptional silencer and a small intervening sequence stem from the gene SMIM45. The two ORFs show different translational characteristics, and they also show divergent patterns of evolutionary development. The studies presented here describe the evolution of the components of SMIM45. One ORF consists of an ultra-conserved 68 amino acid (aa) sequence, whose origins can be traced beyond the evolutionary age of divergence of the elephant shark, ~462 MYA. The silencer also has ancient origins, but it has a complex and divergent pattern of evolutionary formation, as it overlaps both at the 68 aa ORF and the intervening sequence. The other ORF consists of 107 aa. It develops during primate evolution but is found to originate de novo from an ancestral non-coding genomic region with root origins within the Afrothere clade of placental mammals, whose evolutionary age of divergence is ~99 MYA. The formation of the complete 107 aa ORF during primate evolution is outlined, whereby sequence development is found to occur through biased mutations, with disruptive random mutations that also occur but lead to a dead-end. The 107 aa ORF is of particular significance, as there is evidence to suggest it is a protein that may function in human brain development. Its evolutionary formation presents a view of a human-specific ORF and its linked silencer that were predetermined in non-primate ancestral species. The genomic position of the silencer offers interesting possibilities for the regulation of transcription of the 107 aa ORF. A hypothesis is presented with respect to possible spatiotemporal expression of the 107 aa ORF in embryonic tissues.

Iterative assay for transposase-accessible chromatin by sequencing to isolate functionally relevant neuronal subtypes

The Drosophila brain contains tens of thousands of distinct cell types. Thousands of different transgenic lines reproducibly target specific neuron subsets, yet most still express in several cell types. Furthermore, most lines were developed without a priori knowledge of where the transgenes would be expressed. To aid in the development of cell type-specific tools for neuronal identification and manipulation, we developed an iterative assay for transposase-accessible chromatin (ATAC) approach. Open chromatin regions (OCRs) enriched in neurons, compared to whole bodies, drove transgene expression preferentially in subsets of neurons. A second round of ATAC-seq from these specific neuron subsets revealed additional enriched OCR2s that further restricted transgene expression within the chosen neuron subset. This approach allows for continued refinement of transgene expression, and we used it to identify neurons relevant for sleep behavior. Furthermore, this approach is widely applicable to other cell types and to other organisms.

Pathogenic variants in CRX have distinct cis-regulatory effects on enhancers and silencers in photoreceptors

Dozens of variants in the gene for the homeodomain transcription factor (TF) cone-rod homeobox (CRX) are linked with human blinding diseases that vary in their severity and age of onset. How different variants in this single TF alter its function in ways that lead to a range of phenotypes is unclear. We characterized the effects of human disease-causing variants on CRX cis-regulatory function by deploying massively parallel reporter assays (MPRAs) in mouse retina explants carrying knock-ins of two variants, one in the DNA-binding domain (p.R90W) and the other in the transcriptional effector domain (p.E168d2). The degree of reporter gene dysregulation in these mutant Crx retinas corresponds with their phenotypic severity. The two variants affect similar sets of enhancers, and p.E168d2 has distinct effects on silencers. Cis-regulatory elements (CREs) near cone photoreceptor genes are enriched for silencers that are derepressed in the presence of p.E168d2. Chromatin environments of CRX-bound loci are partially predictive of episomal MPRA activity, and distal elements whose accessibility increases later in retinal development are enriched for CREs with silencer activity. We identified a set of potentially pleiotropic regulatory elements that convert from silencers to enhancers in retinas that lack a functional CRX effector domain. Our findings show that phenotypically distinct variants in different domains of CRX have partially overlapping effects on its cis-regulatory function, leading to misregulation of similar sets of enhancers while having a qualitatively different impact on silencers.

A screen for regeneration-associated silencer regulatory elements in zebrafish

Regeneration involves gene expression changes explained in part by context-dependent recruitment of transcriptional activators to distal enhancers. Silencers that engage repressive transcriptional complexes are less studied than enhancers and more technically challenging to validate, but they potentially have profound biological importance for regeneration. Here, we identified candidate silencers through a screening process that examined the ability of DNA sequences to limit injury-induced gene expression in larval zebrafish after fin amputation. A short sequence (s1) on chromosome 5 near several genes that reduce expression during adult fin regeneration could suppress promoter activity in stable transgenic lines and diminish nearby gene expression in knockin lines. High-resolution analysis of chromatin organization identified physical associations of s1 with gene promoters occurring preferentially during fin regeneration, and genomic deletion of s1 elevated the expression of these genes after fin amputation. Our study provides methods to identify "tissue regeneration silencer elements" (TRSEs) with the potential to reduce unnecessary or deleterious gene expression during regeneration.

Canalizing cell fate by transcriptional repression

Precision in the establishment and maintenance of cellular identities is crucial for the development of multicellular organisms and requires tight regulation of gene expression. While extensive research has focused on understanding cell type-specific gene activation, the complex mechanisms underlying the transcriptional repression of alternative fates are not fully understood. Here, we provide an overview of the repressive mechanisms involved in cell fate regulation. We discuss the molecular machinery responsible for suppressing alternative fates and highlight the crucial role of sequence-specific transcription factors (TFs) in this process. Depletion of these TFs can result in unwanted gene expression and increased cellular plasticity. We suggest that these TFs recruit cell type-specific repressive complexes to their cis-regulatory elements, enabling them to modulate chromatin accessibility in a context-dependent manner. This modulation effectively suppresses master regulators of alternative fate programs and their downstream targets. The modularity and dynamic behavior of these repressive complexes enables a limited number of repressors to canalize and maintain major and minor cell fate decisions at different stages of development.

Cis-regulatory modes of Ultrabithorax inactivation in butterfly forewings

Hox gene clusters encode transcription factors that drive regional specialization during animal development: for example the Hox factor Ubx is expressed in the insect metathoracic (T3) wing appendages and differentiates them from T2 mesothoracic identities. Hox transcriptional regulation requires silencing activities that prevent spurious activation and regulatory crosstalks in the wrong tissues, but this has seldom been studied in insects other than Drosophila, which shows a derived Hox dislocation into two genomic clusters that disjoined Antennapedia (Antp) and Ultrabithorax (Ubx). Here, we investigated how Ubx is restricted to the hindwing in butterflies, amidst a contiguous Hox cluster. By analysing Hi-C and ATAC-seq data in the butterfly Junonia coenia, we show that a Topologically Associated Domain (TAD) maintains a hindwing-enriched profile of chromatin opening around Ubx. This TAD is bordered by a Boundary Element (BE) that separates it from a region of joined wing activity around the Antp locus. CRISPR mutational perturbation of this BE releases ectopic Ubx expression in forewings, inducing homeotic clones with hindwing identities. Further mutational interrogation of two non-coding RNA encoding regions and one putative cis-regulatory module within the Ubx TAD cause rare homeotic transformations in both directions, indicating the presence of both activating and repressing chromatin features. We also describe a series of spontaneous forewing homeotic phenotypes obtained in Heliconius butterflies, and discuss their possible mutational basis. By leveraging the extensive wing specialization found in butterflies, our initial exploration of Ubx regulation demonstrates the existence of silencing and insulating sequences that prevent its spurious expression in forewings.

Chromatin structure and 3D architecture define differential functions of PU.1 cis regulatory elements in human blood cell lineages

The precise spatio-temporal expression of the hematopoietic ETS transcription factor PU.1 that determines the hematopoietic cell fates is tightly regulated at the chromatin level. However, it remains elusive as to how chromatin signatures are linked to this dynamic expression pattern of PU.1 across blood cell lineages. Here we performed an unbiased and in-depth analysis of the relationship between human PU.1 expression, the presence of trans-acting factors, and 3D architecture at various cis-regulatory elements (CRE) proximal to the PU.1 locus. We identified multiple novel CREs at the upstream region of the gene following an integrative inspection for conserved DNA elements at the chromatin-accessible regions in primary human blood lineages. We showed that a subset of CREs localize within a 10 kb-wide cluster that exhibits that exhibit molecular features of a myeloid-specific super-enhancer involved in mediating PU.1 autoregulation, including open chromatin, unmethylated DNA, histone enhancer marks, transcription of enhancer RNAs, and occupancy of the PU.1 protein itself. Importantly, we revealed the presence of common 35-kb-wide CTCF-bound insulated neighborhood that contains the CRE cluster, forming the chromatin territory for lineage-specific and CRE-mediated chromatin interactions. These include functional CRE-promoter interactions in myeloid and B cells but not in erythroid and T cells. Our findings also provide mechanistic insights into the interplay between dynamic chromatin structure and 3D architecture in defining certain CREs as enhancers or silencers in chromatin regulation of PU.1 expression. The study lays the groundwork for further examination of PU.1 CREs as well as epigenetic regulation in malignant hematopoiesis.

Pathogenic variants in Crx have distinct cis-regulatory effects on enhancers and silencers in photoreceptors

Dozens of variants in the photoreceptor-specific transcription factor (TF) CRX are linked with human blinding diseases that vary in their severity and age of onset. It is unclear how different variants in this single TF alter its function in ways that lead to a range of phenotypes. We examined the effects of human disease-causing variants on CRX cis-regulatory function by deploying massively parallel reporter assays (MPRAs) in live mouse retinas carrying knock-ins of two variants, one in the DNA binding domain (p.R90W) and the other in the transcriptional effector domain (p.E168d2). The degree of reporter gene dysregulation caused by the variants corresponds with their phenotypic severity. The two variants affect similar sets of enhancers, while p.E168d2 has stronger effects on silencers. Cis-regulatory elements (CREs) near cone photoreceptor genes are enriched for silencers that are de-repressed in the presence of p.E168d2. Chromatin environments of CRX-bound loci were partially predictive of episomal MPRA activity, and silencers were notably enriched among distal elements whose accessibility increases later in retinal development. We identified a set of potentially pleiotropic regulatory elements that convert from silencers to enhancers in retinas that lack a functional CRX effector domain. Our findings show that phenotypically distinct variants in different domains of CRX have partially overlapping effects on its cis-regulatory function, leading to misregulation of similar sets of enhancers, while having a qualitatively different impact on silencers.

Methylation-directed regulatory networks determine enhancing and silencing of mutation disease driver genes and explain inter-patient expression variation

BACKGROUND

Common diseases manifest differentially between patients, but the genetic origin of this variation remains unclear. To explore possible involvement of gene transcriptional-variation, we produce a DNA methylation-oriented, driver-gene-wide dataset of regulatory elements in human glioblastomas and study their effect on inter-patient gene expression variation.

RESULTS

In 175 of 177 analyzed gene regulatory domains, transcriptional enhancers and silencers are intermixed. Under experimental conditions, DNA methylation induces enhancers to alter their enhancing effects or convert into silencers, while silencers are affected inversely. High-resolution mapping of the association between DNA methylation and gene expression in intact genomes reveals methylation-related regulatory units (average size = 915.1 base-pairs). Upon increased methylation of these units, their target-genes either increased or decreased in expression. Gene-enhancing and silencing units constitute cis-regulatory networks of genes. Mathematical modeling of the networks highlights indicative methylation sites, which signified the effect of key regulatory units, and add up to make the overall transcriptional effect of the network. Methylation variation in these sites effectively describe inter-patient expression variation and, compared with DNA sequence-alterations, appears as a major contributor of gene-expression variation among glioblastoma patients.

CONCLUSIONS

We describe complex cis-regulatory networks, which determine gene expression by summing the effects of positive and negative transcriptional inputs. In these networks, DNA methylation induces both enhancing and silencing effects, depending on the context. The revealed mechanism sheds light on the regulatory role of DNA methylation, explains inter-individual gene-expression variation, and opens the way for monitoring the driving forces behind deferential courses of cancer and other diseases.

Three-dimensional genome architecture coordinates key regulators of lineage specification in mammary epithelial cells

Although lineage-specific genes have been identified in the mammary gland, little is known about the contribution of the 3D genome organization to gene regulation in the epithelium. Here, we describe the chromatin landscape of the three major epithelial subsets through integration of long- and short-range chromatin interactions, accessibility, histone modifications, and gene expression. While basal genes display exquisite lineage specificity via distal enhancers, luminal-specific genes show widespread promoter priming in basal cells. Cell specificity in luminal progenitors is largely mediated through extensive chromatin interactions with super-enhancers in gene-body regions in addition to interactions with polycomb silencer elements. Moreover, lineage-specific transcription factors appear to be controlled through cell-specific chromatin interactivity. Finally, chromatin accessibility rather than interactivity emerged as a defining feature of the activation of quiescent basal stem cells. This work provides a comprehensive resource for understanding the role of higher-order chromatin interactions in cell-fate specification and differentiation in the adult mouse mammary gland.

Exploring plant cis-regulatory elements at single-cell resolution: overcoming biological and computational challenges to advance plant research

Cis-regulatory elements (CREs) are important sequences for gene expression and for plant biological processes such as development, evolution, domestication, and stress response. However, studying CREs in plant genomes has been challenging. The totipotent nature of plant cells, coupled with the inability to maintain plant cell types in culture and the inherent technical challenges posed by the cell wall has limited our understanding of how plant cell types acquire and maintain their identities and respond to the environment via CRE usage. Advances in single-cell epigenomics have revolutionized the field of identifying cell-type-specific CREs. These new technologies have the potential to significantly advance our understanding of plant CRE biology, and shed light on how the regulatory genome gives rise to diverse plant phenomena. However, there are significant biological and computational challenges associated with analyzing single-cell epigenomic datasets. In this review, we discuss the historical and foundational underpinnings of plant single-cell research, challenges, and common pitfalls in the analysis of plant single-cell epigenomic data, and highlight biological challenges unique to plants. Additionally, we discuss how the application of single-cell epigenomic data in various contexts stands to transform our understanding of the importance of CREs in plant genomes.

Active learning of enhancer and silencer regulatory grammar in photoreceptors

Cis-regulatory elements (CREs) direct gene expression in health and disease, and models that can accurately predict their activities from DNA sequences are crucial for biomedicine. Deep learning represents one emerging strategy to model the regulatory grammar that relates CRE sequence to function. However, these models require training data on a scale that exceeds the number of CREs in the genome. We address this problem using active machine learning to iteratively train models on multiple rounds of synthetic DNA sequences assayed in live mammalian retinas. During each round of training the model actively selects sequence perturbations to assay, thereby efficiently generating informative training data. We iteratively trained a model that predicts the activities of sequences containing binding motifs for the photoreceptor transcription factor Cone-rod homeobox (CRX) using an order of magnitude less training data than current approaches. The model's internal confidence estimates of its predictions are reliable guides for designing sequences with high activity. The model correctly identified critical sequence differences between active and inactive sequences with nearly identical transcription factor binding sites, and revealed order and spacing preferences for combinations of motifs. Our results establish active learning as an effective method to train accurate deep learning models of cis-regulatory function after exhausting naturally occurring training examples in the genome.

Theory of chromatin organization maintained by active loop extrusion

The active loop extrusion hypothesis proposes that chromatin threads through the cohesin protein complex into progressively larger loops until reaching specific boundary elements. We build upon this hypothesis and develop an analytical theory for active loop extrusion which predicts that loop formation probability is a nonmonotonic function of loop length and describes chromatin contact probabilities. We validate our model with Monte Carlo and hybrid Molecular Dynamics-Monte Carlo simulations and demonstrate that our theory recapitulates experimental chromatin conformation capture data. Our results support active loop extrusion as a mechanism for chromatin organization and provide an analytical description of chromatin organization that may be used to specifically modify chromatin contact probabilities.

Cis-regulatory atlas of primary human CD4+ T cells

Cis-regulatory elements (CRE) are critical for coordinating gene expression programs that dictate cell-specific differentiation and homeostasis. Recently developed self-transcribing active regulatory region sequencing (STARR-Seq) has allowed for genome-wide annotation of functional CREs. Despite this, STARR-Seq assays are only employed in cell lines, in part, due to difficulties in delivering reporter constructs. Herein, we implemented and validated a STARR-Seq-based screen in human CD4+ T cells using a non-integrating lentiviral transduction system. Lenti-STARR-Seq is the first example of a genome-wide assay of CRE function in human primary cells, identifying thousands of functional enhancers and negative regulatory elements (NREs) in human CD4+ T cells. We find an unexpected difference in nucleosome organization between enhancers and NRE: enhancers are located between nucleosomes, whereas NRE are occupied by nucleosomes in their endogenous locations. We also describe chromatin modification, eRNA production, and transcription factor binding at both enhancers and NREs. Our findings support the idea of silencer repurposing as enhancers in alternate cell types. Collectively, these data suggest that Lenti-STARR-Seq is a successful approach for CRE screening in primary human cell types, and provides an atlas of functional CREs in human CD4+ T cells.

Changes in locus wide repression underlie the evolution of Drosophila abdominal pigmentation

Changes in gene regulation represent an important path to generate developmental differences affecting anatomical traits. Interspecific divergence in gene expression often results from changes in transcription-stimulating enhancer elements. While gene repression is crucial for precise spatiotemporal expression patterns, the relative contribution of repressive transcriptional silencers to regulatory evolution remains to be addressed. Here, we show that the Drosophila pigmentation gene ebony has mainly evolved through changes in the spatial domains of silencers patterning its abdominal expression. By precisely editing the endogenous ebony locus of D. melanogaster, we demonstrate the requirement of two redundant abdominal enhancers and three silencers that repress the redundant enhancers in a patterned manner. We observe a role for changes in these silencers in every case of ebony evolution observed to date. Our findings suggest that negative regulation by silencers likely has an under-appreciated role in gene regulatory evolution.

Gene silencing dynamics are modulated by transiently active regulatory elements

Transcriptional enhancers have been extensively characterized, but cis-regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread.

Mechanisms of Interaction between Enhancers and Promoters in Three Drosophila Model Systems

In higher eukaryotes, the regulation of developmental gene expression is determined by enhancers, which are often located at a large distance from the promoters they regulate. Therefore, the architecture of chromosomes and the mechanisms that determine the functional interaction between enhancers and promoters are of decisive importance in the development of organisms. Mammals and the model animal Drosophila have homologous key architectural proteins and similar mechanisms in the organization of chromosome architecture. This review describes the current progress in understanding the mechanisms of the formation and regulation of long-range interactions between enhancers and promoters at three well-studied key regulatory loci in Drosophila.

Deciphering the multi-scale, quantitative cis-regulatory code

Uncovering the cis-regulatory code that governs when and how much each gene is transcribed in a given genome and cellular state remains a central goal of biology. Here, we discuss major layers of regulation that influence how transcriptional outputs are encoded by DNA sequence and cellular context. We first discuss how transcription factors bind specific DNA sequences in a dosage-dependent and cooperative manner and then proceed to the cofactors that facilitate transcription factor function and mediate the activity of modular cis-regulatory elements such as enhancers, silencers, and promoters. We then consider the complex and poorly understood interplay of these diverse elements within regulatory landscapes and its relationships with chromatin states and nuclear organization. We propose that a mechanistically informed, quantitative model of transcriptional regulation that integrates these multiple regulatory layers will be the key to ultimately cracking the cis-regulatory code.

Enhancer-promoter entanglement explains their transcriptional interdependence

Enhancers not only activate target promoters to stimulate messenger RNA (mRNA) synthesis, but they themselves also undergo transcription to produce enhancer RNAs (eRNAs), the significance of which is not well understood. Transcription at the participating enhancer-promoter pair appears coordinated, but it is unclear why and how. Here, we employ cell-free transcription assays using constructs derived from the human GREB1 locus to demonstrate that transcription at an enhancer and its target promoter is interdependent. This interdependence is observable under conditions where direct enhancer-promoter contact (EPC) takes place. We demonstrate that transcription activation at a participating enhancer-promoter pair is dependent on i) the mutual availability of the enhancer and promoter, ii) the state of transcription at both the enhancer and promoter, iii) local abundance of both eRNA and mRNA, and iv) direct EPC. Our results suggest transcriptional interdependence between the enhancer and the promoter as the basis of their transcriptional concurrence and coordination throughout the genome. We propose a model where transcriptional concurrence, coordination and interdependence are possible if the participating enhancer and promoter are entangled in the form of EPC, reside in a proteinaceous bubble, and utilize shared transcriptional resources and regulatory inputs.

Whole-genome functional characterization of RE1 silencers using a modified massively parallel reporter assay

Both upregulation and downregulation by cis-regulatory elements help modulate precise gene expression. However, our understanding of repressive elements is far more limited than activating elements. To address this gap, we characterized RE1, a group of transcriptional silencers bound by REST, at genome-wide scale using a modified massively parallel reporter assay (MPRAduo). MPRAduo empirically defined a minimal binding strength of REST (REST motif-intrinsic value [m-value]), above which cofactors colocalize and silence transcription. We identified 1,500 human variants that alter RE1 silencing and found that their effect sizes are predictable when they overlap with REST-binding sites above the m-value. Additionally, we demonstrate that non-canonical REST-binding motifs exhibit silencer function only if they precisely align half sites with specific spacer lengths. Our results show mechanistic insights into RE1, which allow us to predict its activity and effect of variants on RE1, providing a paradigm for performing genome-wide functional characterization of transcription-factor-binding sites.

Characterizing cis-regulatory elements using single-cell epigenomics

Cell type-specific gene expression patterns and dynamics during development or in disease are controlled by cis-regulatory elements (CREs), such as promoters and enhancers. Distinct classes of CREs can be characterized by their epigenomic features, including DNA methylation, chromatin accessibility, combinations of histone modifications and conformation of local chromatin. Tremendous progress has been made in cataloguing CREs in the human genome using bulk transcriptomic and epigenomic methods. However, single-cell epigenomic and multi-omic technologies have the potential to provide deeper insight into cell type-specific gene regulatory programmes as well as into how they change during development, in response to environmental cues and through disease pathogenesis. Here, we highlight recent advances in single-cell epigenomic methods and analytical tools and discuss their readiness for human tissue profiling.

A global high-density chromatin interaction network reveals functional long-range and trans-chromosomal relationships

BACKGROUND

Chromatin contacts are essential for gene-expression regulation; however, obtaining a high-resolution genome-wide chromatin contact map is still prohibitively expensive owing to large genome sizes and the quadratic scale of pairwise data. Chromosome conformation capture (3C)-based methods such as Hi-C have been extensively used to obtain chromatin contacts. However, since the sparsity of these maps increases with an increase in genomic distance between contacts, long-range or trans-chromatin contacts are especially challenging to sample.

RESULTS

Here, we create a high-density reference genome-wide chromatin contact map using a meta-analytic approach. We integrate 3600 human, 6700 mouse, and 500 fly Hi-C experiments to create species-specific meta-Hi-C chromatin contact maps with 304 billion, 193 billion, and 19 billion contacts in respective species. We validate that meta-Hi-C contact maps are uniquely powered to capture functional chromatin contacts in both cis and trans. We find that while individual dataset Hi-C networks are largely unable to predict any long-range coexpression (median 0.54 AUC), meta-Hi-C networks perform comparably in both cis and trans (0.65 AUC vs 0.64 AUC). Similarly, for long-range expression quantitative trait loci (eQTL), meta-Hi-C contacts outperform all individual Hi-C experiments, providing an improvement over the conventionally used linear genomic distance-based association. Assessing between species, we find patterns of chromatin contact conservation in both cis and trans and strong associations with coexpression even in species for which Hi-C data is lacking.

CONCLUSIONS

We have generated an integrated chromatin interaction network which complements a large number of methodological and analytic approaches focused on improved specificity or interpretation. This high-depth "super-experiment" is surprisingly powerful in capturing long-range functional relationships of chromatin interactions, which are now able to predict coexpression, eQTLs, and cross-species relationships. The meta-Hi-C networks are available at https://labshare.cshl.edu/shares/gillislab/resource/HiC/ .

Building integrative functional maps of gene regulation

Every cell in the human body inherits a copy of the same genetic information. The three billion base pairs of DNA in the human genome, and the roughly 50 000 coding and non-coding genes they contain, must thus encode all the complexity of human development and cell and tissue type diversity. Differences in gene regulation, or the modulation of gene expression, enable individual cells to interpret the genome differently to carry out their specific functions. Here we discuss recent and ongoing efforts to build gene regulatory maps, which aim to characterize the regulatory roles of all sequences in a genome. Many researchers and consortia have identified such regulatory elements using functional assays and evolutionary analyses; we discuss the results, strengths and shortcomings of their approaches. We also discuss new techniques the field can leverage and emerging challenges it will face while striving to build gene regulatory maps of ever-increasing resolution and comprehensiveness.

Origin, evolution, and tissue-specific functions of the porcine repetitive element 1

Background

The porcine repetitive element 1 (PRE1) is the most abundant short interspersed nuclear element (SINE) in the Sus scrofa genome and it has been suggested that some PRE1 can have regulatory functions. The million copies of PRE1 in the porcine genome have accumulated abundant CpG dinucleotides and unique structural variations, such as direct repeats and patterns of sequence degeneration. The aims of this study were to analyse these structural variations to trace the origin and evolutionary pattern of PRE1 and to investigate potential methylation-related functions of PRE1 based on methylation patterns of PRE1 CpG dinucleotides in different tissues.

Results

We investigated the evolutionary trajectory of PRE1 and found that PRE1 originated from the ancestral CHRS-S1 family through three main successive partial duplications. We found that the partial duplications and deletions of PRE1 were likely due to RNA splicing events during retrotransposition. Functionally, correlation analysis showed that the methylation levels of 103 and 261 proximal PRE1 were, respectively, negatively and positively correlated with the expression levels of neighboring genes (Spearman correlation, P < 0.01). Further epigenomic analysis revealed that, in the testis, demethylation of proximal PRE1 in the HORMAD1 and HACD3 genes had tissue-specific enhancer and promoter functions, while in the muscle, methylation of proximal PRE1 repeats in the TCEA3 gene had an enhancer function.

Conclusions

The characteristic sequences of PRE1 reflect unique patterns of origin and evolution and provide a structural basis for diverse regulatory functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00745-3.

REDfly: An Integrated Knowledgebase for Insect Regulatory Genomics

We provide here an updated description of the REDfly (Regulatory Element Database for Fly) database of transcriptional regulatory elements, a unique resource that provides regulatory annotation for the genome of Drosophila and other insects. The genomic sequences regulating insect gene expression-transcriptional cis-regulatory modules (CRMs, e.g., "enhancers") and transcription factor binding sites (TFBSs)-are not currently curated by any other major database resources. However, knowledge of such sequences is important, as CRMs play critical roles with respect to disease as well as normal development, phenotypic variation, and evolution. Characterized CRMs also provide useful tools for both basic and applied research, including developing methods for insect control. REDfly, which is the most detailed existing platform for metazoan regulatory-element annotation, includes over 40,000 experimentally verified CRMs and TFBSs along with their DNA sequences, their associated genes, and the expression patterns they direct. Here, we briefly describe REDfly's contents and data model, with an emphasis on the new features implemented since 2020. We then provide an illustrated walk-through of several common REDfly search use cases.

Identification and classification of cis-regulatory elements in the amphipod crustacean Parhyale hawaiensis

Emerging research organisms enable the study of biology that cannot be addressed using classical 'model' organisms. New data resources can accelerate research in such animals. Here, we present new functional genomic resources for the amphipod crustacean Parhyale hawaiensis, facilitating the exploration of gene regulatory evolution using this emerging research organism. We use Omni-ATAC-seq to identify accessible chromatin genome-wide across a broad time course of Parhyale embryonic development. This time course encompasses many major morphological events, including segmentation, body regionalization, gut morphogenesis and limb development. In addition, we use short- and long-read RNA-seq to generate an improved Parhyale genome annotation, enabling deeper classification of identified regulatory elements. We discover differential accessibility, predict nucleosome positioning, infer transcription factor binding, cluster peaks based on accessibility dynamics, classify biological functions and correlate gene expression with accessibility. Using a Minos transposase reporter system, we demonstrate the potential to identify novel regulatory elements using this approach. This work provides a platform for the identification of novel developmental regulatory elements in Parhyale, and offers a framework for performing such experiments in other emerging research organisms.

Analysis of the landscape of human enhancer sequences in biological databases

The process of gene regulation extends as a network in which both genetic sequences and proteins are involved. The levels of regulation and the mechanisms involved are multiple. Transcription is the main control mechanism for most genes, being the downstream steps responsible for refining the transcription patterns. In turn, gene transcription is mainly controlled by regulatory events that occur at promoters and enhancers. Several studies are focused on analyzing the contribution of enhancers in the development of diseases and their possible use as therapeutic targets. The study of regulatory elements has advanced rapidly in recent years with the development and use of next generation sequencing techniques. All this information has generated a large volume of information that has been transferred to a growing number of public repositories that store this information. In this article, we analyze the content of those public repositories that contain information about human enhancers with the aim of detecting whether the knowledge generated by scientific research is contained in those databases in a way that could be computationally exploited. The analysis will be based on three main aspects identified in the literature: types of enhancers, type of evidence about the enhancers, and methods for detecting enhancer-promoter interactions. Our results show that no single database facilitates the optimal exploitation of enhancer data, most types of enhancers are not represented in the databases and there is need for a standardized model for enhancers. We have identified major gaps and challenges for the computational exploitation of enhancer data.

Harnessing changes in open chromatin determined by ATAC-seq to generate insulin-responsive reporter constructs

BACKGROUND

Gene regulation is critical for proper cellular function. Next-generation sequencing technology has revealed the presence of regulatory networks that regulate gene expression and essential cellular functions. Studies investigating the epigenome have begun to uncover the complex mechanisms regulating transcription. Assay for transposase-accessible chromatin by sequencing (ATAC-seq) is quickly becoming the assay of choice for many epigenomic investigations. However, whether intervention-mediated changes in accessible chromatin determined by ATAC-seq can be harnessed to generate intervention-inducible reporter constructs has not been systematically assayed.

RESULTS

We used the insulin signaling pathway as a model to investigate chromatin regions and gene expression changes using ATAC- and RNA-seq in insulin-treated Drosophila S2 cells. We found correlations between ATAC- and RNA-seq data, especially when stratifying differentially-accessible chromatin regions by annotated feature type. In particular, our data demonstrated a weak but significant correlation between chromatin regions annotated to enhancers (1-2 kb from the transcription start site) and downstream gene expression. We cloned candidate enhancer regions upstream of luciferase and demonstrate insulin-inducibility of several of these reporters.

CONCLUSIONS

Insulin-induced chromatin accessibility determined by ATAC-seq reveals enhancer regions that drive insulin-inducible reporter gene expression.

Long-Distance Repression by Human Silencers: Chromatin Interactions and Phase Separation in Silencers

Three-dimensional genome organization represents an additional layer in the epigenetic regulation of gene expression. Active transcription controlled by enhancers or super-enhancers has been extensively studied. Enhancers or super-enhancers can recruit activators or co-activators to activate target gene expression through long-range chromatin interactions. Chromatin interactions and phase separation play important roles in terms of enhancer or super-enhancer functioning. Silencers are another major type of cis-regulatory element that can mediate gene regulation by turning off or reducing gene expression. However, compared to active transcription, silencer studies are still in their infancy. This review covers the current knowledge of human silencers, especially the roles of chromatin interactions and phase separation in silencers. This review also proposes future directions for human silencer studies.

Structural Variation at a Disease Mutation Hotspot: Strategies to Investigate Gene Regulation and the 3D Genome

A rare form of X-linked Charcot-Marie-Tooth neuropathy, CMTX3, is caused by an interchromosomal insertion occurring at chromosome Xq27.1. Interestingly, eight other disease phenotypes have been associated with insertions (or insertion-deletions) occurring at the same genetic locus. To date, the pathogenic mechanism underlying most of these diseases remains unsolved, although local gene dysregulation has clearly been implicated in at least two phenotypes. The challenges of accessing disease-relevant tissue and modelling these complex genomic rearrangements has led to this research impasse. We argue that recent technological advancements can overcome many of these challenges, particularly induced pluripotent stem cells (iPSC) and their capacity to provide access to patient-derived disease-relevant tissue. However, to date these valuable tools have not been utilized to investigate the disease-associated insertions at chromosome Xq27.1. Therefore, using CMTX3 as a reference disease, we propose an experimental approach that can be used to explore these complex mutations, as well as similar structural variants located elsewhere in the genome. The mutational hotspot at Xq27.1 is a valuable disease paradigm with the potential to improve our understanding of the pathogenic consequences of complex structural variation, and more broadly, refine our knowledge of the multifaceted process of long-range gene regulation. Intergenic structural variation is a critically understudied class of mutation, although it is likely to contribute significantly to unsolved genetic disease.

The role of the epidermis enhancer element in positive and negative transcriptional regulation of ebony in Drosophila melanogaster

The spatiotemporal regulation of gene expression is essential to ensure robust phenotypic outcomes. Pigmentation patterns in Drosophila are determined by pigments biosynthesized in the developing epidermis and the cis-regulatory elements of the genes involved in this process are well-characterized. Here, we report that the known primary epidermal enhancer is dispensable for the transcriptional activation of ebony (involved in light-colored pigment synthesis) in the developing epidermis of Drosophila melanogaster. The evidence was obtained by introducing an approximately 1 kbp deletion at the primary epidermal enhancer by genome editing. The effect of the primary epidermal enhancer deletion on pigmentation and on the endogenous expression pattern of a mCherry-fused ebony allele was examined in the abdomen. The expression levels of the mCherry-fused ebony in the primary epidermal enhancer-deleted strains were slightly higher than that of the control strain, indicating that the sequences outside the primary epidermal enhancer have an ability to drive an expression of this gene in the epidermis. Interestingly, the primary epidermal enhancer deletion resulted in a derepression of this gene in the dorsal midline of the abdominal tergites, where dark pigmentation is present in the wild-type individuals. This indicated that the primary epidermal enhancer fragment contains a silencer. Furthermore, the endogenous expression pattern of ebony in the 2 additional strains with partially deleted primary epidermal enhancer revealed that the silencer resides within a 351-bp fragment in the 5' portion of the primary epidermal enhancer. These results demonstrated that deletion assays combined with reporter assays are highly effective in detecting the presence of positively and negatively regulating sequences within and outside the focal cis-regulatory elements.

Sex-specific evolution of a Drosophila sensory system via interacting cis- and trans-regulatory changes

The evolution of gene expression via cis-regulatory changes is well established as a major driver of phenotypic evolution. However, relatively little is known about the influence of enhancer architecture and intergenic interactions on regulatory evolution. We address this question by examining chemosensory system evolution in Drosophila. Drosophila prolongata males show a massively increased number of chemosensory bristles compared to females and males of sibling species. This increase is driven by sex-specific transformation of ancestrally mechanosensory organs. Consistent with this phenotype, the Pox neuro transcription factor (Poxn), which specifies chemosensory bristle identity, shows expanded expression in D. prolongata males. Poxn expression is controlled by nonadditive interactions among widely dispersed enhancers. Although some D. prolongata Poxn enhancers show increased activity, the additive component of this increase is slight, suggesting that most changes in Poxn expression are due to epistatic interactions between Poxn enhancers and trans-regulatory factors. Indeed, the expansion of D. prolongata Poxn enhancer activity is only observed in cells that express doublesex (dsx), the gene that controls sexual differentiation in Drosophila and also shows increased expression in D. prolongata males due to cis-regulatory changes. Although expanded dsx expression may contribute to increased activity of D. prolongata Poxn enhancers, this interaction is not sufficient to explain the full expansion of Poxn expression, suggesting that cis-trans interactions between Poxn, dsx, and additional unknown genes are necessary to produce the derived D. prolongata phenotype. Overall, our results demonstrate the importance of epistatic gene interactions for evolution, particularly when pivotal genes have complex regulatory architecture.

Regulating specificity in enhancer-promoter communication

Enhancers are cis-regulatory elements that can activate transcription remotely to regulate a specific pattern of a gene's expression. Genes typically have many enhancers that are often intermingled in the loci of other genes. To regulate expression, enhancers must therefore activate their correct promoter while ignoring others that may be in closer linear proximity. In this review, we discuss mechanisms by which enhancers engage with promoters, including recent findings on the role of cohesin and the Mediator complex, and how this specificity in enhancer-promoter communication is encoded. Genetic dissection of model loci, in addition to more recent findings using genome-wide approaches, highlight the core promoter sequence, its accessibility, cofactor-promoter preference, in addition to the surrounding genomic context, as key components.