Uncovering cis-regulatory sequence requirements for context-specific transcription factor binding

Systematic dissection of sequence features affecting binding specificity of a pioneer factor reveals binding synergy between FOXA1 and AP-1

Despite the unique ability of pioneer factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called chromatin immunoprecipitation with integrated synthetic oligonucleotides (ChIP-ISO) to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1, in human A549 cells. Combining ChIP-ISO with in vitro and neural network analyses, we find that (1) FOXA1 binding is strongly affected by co-binding transcription factors (TFs) AP-1 and CEBPB; (2) FOXA1 and AP-1 show binding cooperativity in vitro; (3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin; and (4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

Context dependent activity of p63-bound gene regulatory elements

The p53 family of transcription factors regulate numerous organismal processes including the development of skin and limbs, ciliogenesis, and preservation of genetic integrity and tumor suppression. p53 family members control these processes and gene expression networks through engagement with DNA sequences within gene regulatory elements. Whereas p53 binding to its cognate recognition sequence is strongly associated with transcriptional activation, p63 can mediate both activation and repression. How the DNA sequence of p63-bound gene regulatory elements is linked to these varied activities is not yet understood. Here, we use massively parallel reporter assays (MPRA) in a range of cellular and genetic contexts to investigate the influence of DNA sequence on p63-mediated transcription. Most regulatory elements with a p63 response element motif (p63RE) activate transcription, with those sites bound by p63 more frequently or adhering closer to canonical p53 family response element sequences driving higher transcriptional output. The most active regulatory elements are those also capable of binding p53. Elements uniquely bound by p63 have varied activity, with p63RE-mediated repression associated with lower overall GC content in flanking sequences. Comparison of activity across cell lines suggests differential activity of elements may be regulated by a combination of p63 abundance or context-specific cofactors. Finally, changes in p63 isoform expression dramatically alters regulatory element activity, primarily shifting inactive elements towards a strong p63-dependent activity. Our analysis of p63-bound gene regulatory elements provides new insight into how sequence, cellular context, and other transcription factors influence p63-dependent transcription. These studies provide a framework for understanding how p63 genomic binding locally regulates transcription. Additionally, these results can be extended to investigate the influence of sequence content, genomic context, chromatin structure on the interplay between p63 isoforms and p53 family paralogs.

A fine kinetic balance of interactions directs transcription factor hubs to genes

Eukaryotic gene regulation relies on the binding of sequence-specific transcription factors (TFs). TFs bind chromatin transiently yet occupy their target sites by forming high-local concentration microenvironments (hubs and condensates) that increase the frequency of binding events. Despite their ubiquity, such microenvironments have been difficult to study in endogenous contexts due to technical limitations. Here, we overcome these limitations and investigate how hubs drive TF occupancy at their targets. Using a DNA binding perturbation to a hub-forming TF, Zelda, in Drosophila embryos, we find that hub properties, including the stability and frequencies of associations to targets, are key determinants of TF occupancy. Our data suggest that the targeting of these hubs is driven not just by specific DNA motif recognition, but also by a fine-tuned kinetic balance of interactions between TFs and their co-binding partners.

Protein-intrinsic properties and context-dependent effects regulate pioneer factor binding and function

Chromatin is a barrier to the binding of many transcription factors. By contrast, pioneer factors access nucleosomal targets and promote chromatin opening. Despite binding to target motifs in closed chromatin, many pioneer factors display cell-type-specific binding and activity. The mechanisms governing pioneer factor occupancy and the relationship between chromatin occupancy and opening remain unclear. We studied three Drosophila transcription factors with distinct DNA-binding domains and biological functions: Zelda, Grainy head and Twist. We demonstrated that the level of chromatin occupancy is a key determinant of pioneering activity. Multiple factors regulate occupancy, including motif content, local chromatin and protein concentration. Regions outside the DNA-binding domain are required for binding and chromatin opening. Our results show that pioneering activity is not a binary feature intrinsic to a protein but occurs on a spectrum and is regulated by a variety of protein-intrinsic and cell-type-specific features.

Multiomics-integrated deep language model enables in silico genome-wide detection of transcription factor binding site in unexplored biosamples

MOTIVATION

Transcription factor binding sites (TFBS) are regulatory elements that have significant impact on transcription regulation and cell fate determination. Canonical motifs, biological experiments, and computational methods have made it possible to discover TFBS. However, most existing in silico TFBS prediction models are solely DNA-based, and are trained and utilized within the same biosample, which fail to infer TFBS in experimentally unexplored biosamples.

RESULTS

Here, we propose TFBS prediction by modified TransFormer (TFTF), a multimodal deep language architecture which integrates multiomics information in epigenetic studies. In comparison to existing computational techniques, TFTF has state-of-the-art accuracy, and is also the first approach to accurately perform genome-wide detection for cell-type and species-specific TFBS in experimentally unexplored biosamples. Compared to peak calling methods, TFTF consistently discovers true TFBS in threshold tuning-free way, with higher recalled rates. The underlying mechanism of TFTF reveals greater attention to the targeted TF's motif region in TFBS, and general attention to the entire peak region in non-TFBS. TFTF can benefit from the integration of broader and more diverse data for improvement and can be applied to multiple epigenetic scenarios.

AVAILABILITY AND IMPLEMENTATION

We provide a web server (https://tftf.ibreed.cn/) for users to utilize TFTF model. Users can train TFTF model and discover TFBS with their own data.

Chromatin accessibility in the Drosophila embryo is determined by transcription factor pioneering and enhancer activation

Chromatin accessibility is integral to the process by which transcription factors (TFs) read out cis-regulatory DNA sequences, but it is difficult to differentiate between TFs that drive accessibility and those that do not. Deep learning models that learn complex sequence rules provide an unprecedented opportunity to dissect this problem. Using zygotic genome activation in Drosophila as a model, we analyzed high-resolution TF binding and chromatin accessibility data with interpretable deep learning and performed genetic validation experiments. We identify a hierarchical relationship between the pioneer TF Zelda and the TFs involved in axis patterning. Zelda consistently pioneers chromatin accessibility proportional to motif affinity, whereas patterning TFs augment chromatin accessibility in sequence contexts where they mediate enhancer activation. We conclude that chromatin accessibility occurs in two tiers: one through pioneering, which makes enhancers accessible but not necessarily active, and the second when the correct combination of TFs leads to enhancer activation.

Setting the stage for development: the maternal-to-zygotic transition in Drosophila

The zygote has a daunting task ahead of itself; it must develop from a single cell (fertilized egg) into a fully functioning adult with a multitude of different cell types. In the beginning, the zygote has help from its mother, in the form of gene products deposited into the egg, but eventually, it must rely on its own resources to proceed through development. The transfer of developmental control from the mother to the embryo is called the maternal-to-zygotic transition (MZT). All animals undergo this transition, which is defined by two main processes-the degradation of maternal RNAs and the synthesis of new RNAs from the zygote's own genome. Here, we review the regulation of the MZT in Drosophila, but given the broad conservation of this essential process, much of the regulation is shared among metazoans.

Prediction accuracy of regulatory elements from sequence varies by functional sequencing technique

Introduction

Various sequencing based approaches are used to identify and characterize the activities of cis-regulatory elements in a genome-wide fashion. Some of these techniques rely on indirect markers such as histone modifications (ChIP-seq with histone antibodies) or chromatin accessibility (ATAC-seq, DNase-seq, FAIRE-seq), while other techniques use direct measures such as episomal assays measuring the enhancer properties of DNA sequences (STARR-seq) and direct measurement of the binding of transcription factors (ChIP-seq with transcription factor-specific antibodies). The activities of cis-regulatory elements such as enhancers, promoters, and repressors are determined by their sequence and secondary processes such as chromatin accessibility, DNA methylation, and bound histone markers.

Methods

Here, machine learning models are employed to evaluate the accuracy with which cis-regulatory elements identified by various commonly used sequencing techniques can be predicted by their underlying sequence alone to distinguish between cis-regulatory activity that is reflective of sequence content versus secondary processes.

Results and discussion

Models trained and evaluated on D. melanogaster sequences identified through DNase-seq and STARR-seq are significantly more accurate than models trained on sequences identified by H3K4me1, H3K4me3, and H3K27ac ChIP-seq, FAIRE-seq, and ATAC-seq. These results suggest that the activity detected by DNase-seq and STARR-seq can be largely explained by underlying DNA sequence, independent of secondary processes. Experimentally, a subset of DNase-seq and H3K4me1 ChIP-seq sequences were tested for enhancer activity using luciferase assays and compared with previous tests performed on STARR-seq sequences. The experimental data indicated that STARR-seq sequences are substantially enriched for enhancer-specific activity, while the DNase-seq and H3K4me1 ChIP-seq sequences are not. Taken together, these results indicate that the DNase-seq approach identifies a broad class of regulatory elements of which enhancers are a subset and the associated data are appropriate for training models for detecting regulatory activity from sequence alone, STARR-seq data are best for training enhancer-specific sequence models, and H3K4me1 ChIP-seq data are not well suited for training and evaluating sequence-based models for cis-regulatory element prediction.

Protein-intrinsic properties and context-dependent effects regulate pioneer-factor binding and function

Chromatin is a barrier to the binding of many transcription factors. By contrast, pioneer factors access nucleosomal targets and promote chromatin opening. Despite binding to target motifs in closed chromatin, many pioneer factors display cell-type specific binding and activity. The mechanisms governing pioneer-factor occupancy and the relationship between chromatin occupancy and opening remain unclear. We studied three Drosophila transcription factors with distinct DNA-binding domains and biological functions: Zelda, Grainy head, and Twist. We demonstrated that the level of chromatin occupancy is a key determinant of pioneering activity. Multiple factors regulate occupancy, including motif content, local chromatin, and protein concentration. Regions outside the DNA-binding domain are required for binding and chromatin opening. Our results show that pioneering activity is not a binary feature intrinsic to a protein but occurs on a spectrum and is regulated by a variety of protein-intrinsic and cell-type-specific features.

The Drosophila embryo as a tabula rasa for the epigenome

The control of gene expression in eukaryotes relies on how transcription factors and RNA polymerases manipulate the structure of chromatin. These interactions are especially important in development as gene expression programs change. Chromatin generally limits the accessibility of DNA, and thus exposing sequences at regulatory elements is critical for gene expression. However, it is challenging to understand how transcription factors manipulate chromatin structure and the sequence of regulatory events. The Drosophila embryo has provided a powerful setting to directly observe the establishment and elaboration of chromatin features and experimentally test the causality of transcriptional events that are shared among many metazoans. The large embryo is tractable by live imaging, and a variety of well-developed tools allow the manipulation of factors during early development. The early embryo develops as a syncytium with rapid nuclear divisions and no zygotic transcription, with largely featureless chromatin. Thus, studies in this system have revealed the progression of genome activation triggered by pioneer factors that initiate DNA exposure at regulatory elements and the establishment of chromatin domains, including heterochromatin, the nucleolus, and nuclear bodies. The de novo emergence of nuclear structures in the early embryo reveals features of chromatin dynamics that are likely to be central to transcriptional regulation in all cells.

DeepSTARR predicts enhancer activity from DNA sequence and enables the de novo design of synthetic enhancers

Enhancer sequences control gene expression and comprise binding sites (motifs) for different transcription factors (TFs). Despite extensive genetic and computational studies, the relationship between DNA sequence and regulatory activity is poorly understood, and de novo enhancer design has been challenging. Here, we built a deep-learning model, DeepSTARR, to quantitatively predict the activities of thousands of developmental and housekeeping enhancers directly from DNA sequence in Drosophila melanogaster S2 cells. The model learned relevant TF motifs and higher-order syntax rules, including functionally nonequivalent instances of the same TF motif that are determined by motif-flanking sequence and intermotif distances. We validated these rules experimentally and demonstrated that they can be generalized to humans by testing more than 40,000 wildtype and mutant Drosophila and human enhancers. Finally, we designed and functionally validated synthetic enhancers with desired activities de novo.

Analysis of long and short enhancers in melanoma cell states

Understanding how enhancers drive cell-type specificity and efficiently identifying them is essential for the development of innovative therapeutic strategies. In melanoma, the melanocytic (MEL) and the mesenchymal-like (MES) states present themselves with different responses to therapy, making the identification of specific enhancers highly relevant. Using massively parallel reporter assays (MPRAs) in a panel of patient-derived melanoma lines (MM lines), we set to identify and decipher melanoma enhancers by first focusing on regions with state-specific H3K27 acetylation close to differentially expressed genes. An in-depth evaluation of those regions was then pursued by investigating the activity of overlapping ATAC-seq peaks along with a full tiling of the acetylated regions with 190 bp sequences. Activity was observed in more than 60% of the selected regions, and we were able to precisely locate the active enhancers within ATAC-seq peaks. Comparison of sequence content with activity, using the deep learning model DeepMEL2, revealed that AP-1 alone is responsible for the MES enhancer activity. In contrast, SOX10 and MITF both influence MEL enhancer function with SOX10 being required to achieve high levels of activity. Overall, our MPRAs shed light on the relationship between long and short sequences in terms of their sequence content, enhancer activity, and specificity across melanoma cell states.

Spatial organization of transcribing loci during early genome activation in Drosophila

The early Drosophila embryo provides unique experimental advantages for addressing fundamental questions of gene regulation at multiple levels of organization, from individual gene loci to the entire genome. Using 1.5-h-old Drosophila embryos undergoing the first wave of genome activation,¹ we detected ∼110 discrete "speckles" of RNA polymerase II (RNA Pol II) per nucleus, two of which were larger and localized to the histone locus bodies (HLBs).^2,3 In the absence of the primary driver of Drosophila genome activation, the pioneer factor Zelda (Zld),^1,4,5 70% fewer speckles were present; however, the HLBs tended to be larger than wild-type (WT) HLBs, indicating that RNA Pol II accumulates at the HLBs in the absence of robust early-gene transcription. We observed a uniform distribution of distances between active genes in the nuclei of both WT and zld mutant embryos, indicating that early co-regulated genes do not cluster into nuclear sub-domains. However, in instances whereby transcribing genes did come into close 3D proximity (within 400 nm), they were found to have distinct RNA Pol II speckles. In contrast to the emerging model whereby active genes are clustered to facilitate co-regulation and sharing of transcriptional resources, our data support an "individualist" model of gene control at early genome activation in Drosophila. This model is in contrast to a "collectivist" model, where active genes are spatially clustered and share transcriptional resources, motivating rigorous tests of both models in other experimental systems.

The non-coding genome in genetic brain disorders: new targets for therapy?

The non-coding genome, consisting of more than 98% of all genetic information in humans and once judged as ‘Junk DNA’, is increasingly moving into the spotlight in the field of human genetics. Non-coding regulatory elements (NCREs) are crucial to ensure correct spatio-temporal gene expression. Technological advancements have allowed to identify NCREs on a large scale, and mechanistic studies have helped to understand the biological mechanisms underlying their function. It is increasingly becoming clear that genetic alterations of NCREs can cause genetic disorders, including brain diseases. In this review, we concisely discuss mechanisms of gene regulation and how to investigate them, and give examples of non-coding alterations of NCREs that give rise to human brain disorders. The cross-talk between basic and clinical studies enhances the understanding of normal and pathological function of NCREs, allowing better interpretation of already existing and novel data. Improved functional annotation of NCREs will not only benefit diagnostics for patients, but might also lead to novel areas of investigations for targeted therapies, applicable to a wide panel of genetic disorders. The intrinsic complexity and precision of the gene regulation process can be turned to the advantage of highly specific treatments. We further discuss this exciting new field of ‘enhancer therapy’ based on recent examples.

CLAMP regulates zygotic genome activation in Drosophila embryos

Embryonic patterning is critically dependent on zygotic genome activation (ZGA). In Drosophila melanogaster embryos, the pioneer factor Zelda directs ZGA, possibly in conjunction with other factors. Here, we have explored the novel involvement of Chromatin-Linked Adapter for MSL Proteins (CLAMP) during ZGA. CLAMP binds thousands of sites genome-wide throughout early embryogenesis. Interestingly, CLAMP relocates to target promoter sequences across the genome when ZGA is initiated. Although there is a considerable overlap between CLAMP and Zelda binding sites, the proteins display distinct temporal dynamics. To assess whether CLAMP occupancy affects gene expression, we analyzed transcriptomes of embryos zygotically compromised for either clamp or zelda and found that transcript levels of many zygotically activated genes are similarly affected. Importantly, compromising either clamp or zelda disrupted the expression of critical segmentation and sex determination genes bound by CLAMP (and Zelda). Furthermore, clamp knockdown embryos recapitulate other phenotypes observed in Zelda-depleted embryos, including nuclear division defects, centrosome aberrations, and a disorganized actomyosin network. Based on these data, we propose that CLAMP acts in concert with Zelda to regulate early zygotic transcription.

Extremely rapid and reversible optogenetic perturbation of nuclear proteins in living embryos

Many developmental regulators have complex and context-specific roles in different tissues and stages, making the dissection of their function extremely challenging. As regulatory processes often occur within minutes, perturbation methods that match these dynamics are needed. Here, we present the improved light-inducible nuclear export system (iLEXY), an optogenetic loss-of-function approach that triggers translocation of proteins from the nucleus to the cytoplasm. By introducing a series of mutations, we substantially increased LEXY's efficiency and generated variants with different recovery times. iLEXY enables rapid (t1/2 < 30 s), efficient, and reversible nuclear protein depletion in embryos, and is generalizable to proteins of diverse sizes and functions. Applying iLEXY to the Drosophila master regulator Twist, we phenocopy loss-of-function mutants, precisely map the Twist-sensitive embryonic stages, and investigate the effects of timed Twist depletions. Our results demonstrate the power of iLEXY to dissect the function of pleiotropic factors during embryogenesis with unprecedented temporal precision.

Pioneering the developmental frontier

Coordinated changes in gene expression allow a single fertilized oocyte to develop into a complex multi-cellular organism. These changes in expression are controlled by transcription factors that gain access to discrete cis-regulatory elements in the genome, allowing them to activate gene expression. Although nucleosomes present barriers to transcription factor occupancy, pioneer transcription factors have unique properties that allow them to bind DNA in the context of nucleosomes, define cis-regulatory elements, and facilitate the subsequent binding of additional factors that determine gene expression. In this capacity, pioneer factors act at the top of gene-regulatory networks to control developmental transitions. Developmental context also influences pioneer factor binding and activity. Here we discuss the interplay between pioneer factors and development, their role in driving developmental transitions, and the influence of the cellular environment on pioneer factor binding and activity.

GAF is essential for zygotic genome activation and chromatin accessibility in the early Drosophila embryo

Following fertilization, the genomes of the germ cells are reprogrammed to form the totipotent embryo. Pioneer transcription factors are essential for remodeling the chromatin and driving the initial wave of zygotic gene expression. In Drosophila melanogaster, the pioneer factor Zelda is essential for development through this dramatic period of reprogramming, known as the maternal-to-zygotic transition (MZT). However, it was unknown whether additional pioneer factors were required for this transition. We identified an additional maternally encoded factor required for development through the MZT, GAGA Factor (GAF). GAF is necessary to activate widespread zygotic transcription and to remodel the chromatin accessibility landscape. We demonstrated that Zelda preferentially controls expression of the earliest transcribed genes, while genes expressed during widespread activation are predominantly dependent on GAF. Thus, progression through the MZT requires coordination of multiple pioneer-like factors, and we propose that as development proceeds control is gradually transferred from Zelda to GAF.

AP-1 subunits converge promiscuously at enhancers to potentiate transcription

The AP-1 transcription factor (TF) dimer contributes to many biological processes and environmental responses. AP-1 can be composed of many interchangeable subunits. Unambiguously determining the binding locations of these subunits in the human genome is challenging because of variable antibody specificity and affinity. Here, we definitively establish the genome-wide binding patterns of five AP-1 subunits by using CRISPR to introduce a common antibody tag on each subunit. We find limited evidence for strong dimerization preferences between subunits at steady state and find that, under a stimulus, dimerization patterns reflect changes in the transcriptome. Further, our analysis suggests that canonical AP-1 motifs indiscriminately recruit all AP-1 subunits to genomic sites, which we term AP-1 hotspots. We find that AP-1 hotspots are predictive of cell type–specific gene expression and of genomic responses to glucocorticoid signaling (more so than super-enhancers) and are significantly enriched in disease-associated genetic variants. Together, these results support a model where promiscuous binding of many AP-1 subunits to the same genomic location play a key role in regulating cell type–specific gene expression and environmental responses.

ConnecTF: A platform to integrate transcription factor-gene interactions and validate regulatory networks

Deciphering gene regulatory networks (GRNs) is both a promise and challenge of systems biology. The promise lies in identifying key transcription factors (TFs) that enable an organism to react to changes in its environment. The challenge lies in validating GRNs that involve hundreds of TFs with hundreds of thousands of interactions with their genome-wide targets experimentally determined by high-throughput sequencing. To address this challenge, we developed ConnecTF, a species-independent, web-based platform that integrates genome-wide studies of TF-target binding, TF-target regulation, and other TF-centric omic datasets and uses these to build and refine validated or inferred GRNs. We demonstrate the functionality of ConnecTF by showing how integration within and across TF-target datasets uncovers biological insights. Case study 1 uses integration of TF-target gene regulation and binding datasets to uncover TF mode-of-action and identify potential TF partners for 14 TFs in abscisic acid signaling. Case study 2 demonstrates how genome-wide TF-target data and automated functions in ConnecTF are used in precision/recall analysis and pruning of an inferred GRN for nitrogen signaling. Case study 3 uses ConnecTF to chart a network path from NLP7, a master TF in nitrogen signaling, to direct secondary TF2s and to its indirect targets in a Network Walking approach. The public version of ConnecTF (https://ConnecTF.org) contains 3,738,278 TF-target interactions for 423 TFs in Arabidopsis, 839,210 TF-target interactions for 139 TFs in maize (Zea mays), and 293,094 TF-target interactions for 26 TFs in rice (Oryza sativa). The database and tools in ConnecTF will advance the exploration of GRNs in plant systems biology applications for model and crop species.

Epigenetic Regulation of Kupffer Cell Function in Health and Disease

Kupffer cells, the resident macrophages of the liver, comprise the largest pool of tissue macrophages in the body. Within the liver sinusoids Kupffer cells perform functions common across many tissue macrophages including response to tissue damage and antigen presentation. They also engage in specialized activities including iron scavenging and the uptake of opsonized particles from the portal blood. Here, we review recent studies of the epigenetic pathways that establish Kupffer cell identity and function. We describe a model by which liver-environment specific signals induce lineage determining transcription factors necessary for differentiation of Kupffer cells from bone-marrow derived monocytes. We conclude by discussing how these lineage determining transcription factors (LDTFs) drive Kupffer cell behavior during both homeostasis and disease, with particular focus on the relevance of Kupffer cell LDTF pathways in the setting of non-alcoholic fatty liver disease and non-alcoholic steatohepatitis.

Evolution of vertebrate gill covers via shifts in an ancient Pou3f3 enhancer

Whereas the gill chambers of jawless vertebrates open directly into the environment, jawed vertebrates evolved skeletal appendages that drive oxygenated water unidirectionally over the gills. A major anatomical difference between the two jawed vertebrate lineages is the presence of a single large gill cover in bony fishes versus separate covers for each gill chamber in cartilaginous fishes. Here, we find that these divergent patterns correlate with the pharyngeal arch expression of Pou3f3 orthologs. We identify a deeply conserved Pou3f3 arch enhancer present in humans through sharks but undetectable in jawless fish. Minor differences between the bony and cartilaginous fish enhancers account for their restricted versus pan-arch expression patterns. In zebrafish, mutation of Pou3f3 or the conserved enhancer disrupts gill cover formation, whereas ectopic pan-arch Pou3f3b expression generates ectopic skeletal elements resembling the multimeric covers of cartilaginous fishes. Emergence of this Pou3f3 arch enhancer >430 Mya and subsequent modifications may thus have contributed to the acquisition and diversification of gill covers and respiratory strategies during gnathostome evolution.

Transcriptional Enhancers in Drosophila

Key discoveries in Drosophila have shaped our understanding of cellular "enhancers." With a special focus on the fly, this chapter surveys properties of these adaptable cis-regulatory elements, whose actions are critical for the complex spatial/temporal transcriptional regulation of gene expression in metazoa. The powerful combination of genetics, molecular biology, and genomics available in Drosophila has provided an arena in which the developmental role of enhancers can be explored. Enhancers are characterized by diverse low- or high-throughput assays, which are challenging to interpret, as not all of these methods of identifying enhancers produce concordant results. As a model metazoan, the fly offers important advantages to comprehensive analysis of the central functions that enhancers play in gene expression, and their critical role in mediating the production of phenotypes from genotype and environmental inputs. A major challenge moving forward will be obtaining a quantitative understanding of how these cis-regulatory elements operate in development and disease.

Odd-paired is a pioneer-like factor that coordinates with Zelda to control gene expression in embryos

Pioneer factors such as Zelda (Zld) help initiate zygotic transcription in Drosophila early embryos, but whether other factors support this dynamic process is unclear. Odd-paired (Opa), a zinc-finger transcription factor expressed at cellularization, controls the transition of genes from pair-rule to segmental patterns along the anterior-posterior axis. Finding that Opa also regulates expression through enhancer sog_Distal along the dorso-ventral axis, we hypothesized Opa's role is more general. Chromatin-immunoprecipitation (ChIP-seq) confirmed its in vivo binding to sog_Distal but also identified widespread binding throughout the genome, comparable to Zld. Furthermore, chromatin assays (ATAC-seq) demonstrate that Opa, like Zld, influences chromatin accessibility genome-wide at cellularization, suggesting both are pioneer factors with common as well as distinct targets. Lastly, embryos lacking opa exhibit widespread, late patterning defects spanning both axes. Collectively, these data suggest Opa is a general timing factor and likely late-acting pioneer factor that drives a secondary wave of zygotic gene expression.

Regulation of spatiotemporal limits of developmental gene expression via enhancer grammar

The regulatory specificity of a gene is determined by the structure of its enhancers, which contain multiple transcription factor binding sites. A unique combination of transcription factor binding sites in an enhancer determines the boundary of target gene expression, and their disruption often leads to developmental defects. Despite extensive characterization of binding motifs in an enhancer, it is still unclear how each binding site contributes to overall transcriptional activity. Using live imaging, quantitative analysis, and mathematical modeling, we measured the contribution of individual binding sites in transcriptional regulation. We show that binding site arrangement within the Rho-GTPase component t48 enhancer mediates the expression boundary by mainly regulating the timing of transcriptional activation along the dorsoventral axis of Drosophila embryos. By tuning the binding affinity of the Dorsal (Dl) and Zelda (Zld) sites, we show that single site modulations are sufficient to induce significant changes in transcription. Yet, no one site seems to have a dominant role; rather, multiple sites synergistically drive increases in transcriptional activity. Interestingly, Dl and Zld demonstrate distinct roles in transcriptional regulation. Dl site modulations change spatial boundaries of t48, mostly by affecting the timing of activation and bursting frequency rather than transcriptional amplitude or bursting duration. However, modulating the binding site for the pioneer factor Zld affects both the timing of activation and amplitude, suggesting that Zld may potentiate higher Dl recruitment to target DNAs. We propose that such fine-tuning of dynamic gene control via enhancer structure may play an important role in ensuring normal development.

From mother to embryo: A molecular perspective on zygotic genome activation

In animals, the early embryo is mostly transcriptionally silent and development is fueled by maternally supplied mRNAs and proteins. These maternal products are important not only for survival, but also to gear up the zygote's genome for activation. Over the last three decades, research with different model organisms and experimental approaches has identified molecular factors and proposed mechanisms for how the embryo transitions from being transcriptionally silent to transcriptionally competent. In this chapter, we discuss the molecular players that shape the molecular landscape of ZGA and provide insights into their mode of action in activating the transcription program in the developing embryo.

Eukaryotic transcription factors can track and control their target genes using DNA antennas

Eukaryotic transcription factors (TF) function by binding to short 6-10 bp DNA recognition sites located near their target genes, which are scattered through vast genomes. Such process surmounts enormous specificity, efficiency and celerity challenges using a molecular mechanism that remains poorly understood. Combining biophysical experiments, theory and bioinformatics, we dissect the interplay between the DNA-binding domain of Engrailed, a Drosophila TF, and the regulatory regions of its target genes. We find that Engrailed binding affinity is strongly amplified by the DNA regions flanking the recognition site, which contain long tracts of degenerate recognition-site repeats. Such DNA organization operates as an antenna that attracts TF molecules in a promiscuous exchange among myriads of intermediate affinity binding sites. The antenna ensures a local TF supply, enables gene tracking and fine control of the target site's basal occupancy. This mechanism illuminates puzzling gene expression data and suggests novel engineering strategies to control gene expression.

The Role of Chromatin Accessibility in cis-Regulatory Evolution

Transcription factor (TF) binding is determined by sequence as well as chromatin accessibility. Although the role of accessibility in shaping TF-binding landscapes is well recorded, its role in evolutionary divergence of TF binding, which in turn can alter cis-regulatory activities, is not well understood. In this work, we studied the evolution of genome-wide binding landscapes of five major TFs in the core network of mesoderm specification, between Drosophila melanogaster and Drosophila virilis, and examined its relationship to accessibility and sequence-level changes. We generated chromatin accessibility data from three important stages of embryogenesis in both Drosophila melanogaster and Drosophila virilis and recorded conservation and divergence patterns. We then used multivariable models to correlate accessibility and sequence changes to TF-binding divergence. We found that accessibility changes can in some cases, for example, for the master regulator Twist and for earlier developmental stages, more accurately predict binding change than is possible using TF-binding motif changes between orthologous enhancers. Accessibility changes also explain a significant portion of the codivergence of TF pairs. We noted that accessibility and motif changes offer complementary views of the evolution of TF binding and developed a combined model that captures the evolutionary data much more accurately than either view alone. Finally, we trained machine learning models to predict enhancer activity from TF binding and used these functional models to argue that motif and accessibility-based predictors of TF-binding change can substitute for experimentally measured binding change, for the purpose of predicting evolutionary changes in enhancer activity.

Continued Activity of the Pioneer Factor Zelda Is Required to Drive Zygotic Genome Activation

Reprogramming cell fate during the first stages of embryogenesis requires that transcriptional activators gain access to the genome and remodel the zygotic transcriptome. Nonetheless, it is not clear whether the continued activity of these pioneering factors is required throughout zygotic genome activation or whether they are only required early to establish cis-regulatory regions. To address this question, we developed an optogenetic strategy to rapidly and reversibly inactivate the master regulator of genome activation in Drosophila, Zelda. Using this strategy, we demonstrate that continued Zelda activity is required throughout genome activation. We show that Zelda binds DNA in the context of nucleosomes and suggest that this allows Zelda to occupy the genome despite the rapid division cycles in the early embryo. These data identify a powerful strategy to inactivate transcription factor function during development and suggest that reprogramming in the embryo may require specific, continuous pioneering functions to activate the genome.

Evolutionary Potential of Cis-Regulatory Mutations to Cause Rapid Changes in Transcription Factor Binding

Transcriptional regulation is crucial for all biological processes and well investigated at the molecular level for a wide range of organisms. However, it is quite unclear how innovations, such as the activity of a novel regulatory element, evolve. In the case of transcription factor (TF) binding, both a novel TF and a novel-binding site would need to evolve concertedly. Since promiscuous functions have recently been identified as important intermediate steps in creating novel specific functions in many areas such as enzyme evolution and protein-protein interactions, we ask here how promiscuous binding of TFs to TF-binding sites (TFBSs) affects the robustness and evolvability of this tightly regulated system. Specifically, we investigate the binding behavior of several hundred TFs from different species at unprecedented breadth. Our results illustrate multiple aspects of TF-binding interactions, ranging from correlations between the strength of the interaction bond and specificity, to preferences regarding TFBS nucleotide composition in relation to both domains and binding specificity. We identified a subset of high A/T binding motifs. Motifs in this subset had many functionally neutral one-error mutants, and were bound by multiple different binding domains. Our results indicate that, especially for some TF-TFBS associations, low binding specificity confers high degrees of evolvability, that is that few mutations facilitate rapid changes in transcriptional regulation, in particular for large and old TF families. In this study we identify binding motifs exhibiting behavior indicating high evolutionary potential for innovations in transcriptional regulation.

Regulatory principles governing the maternal-to-zygotic transition: insights from Drosophila melanogaster

The onset of metazoan development requires that two terminally differentiated germ cells, a sperm and an oocyte, become reprogrammed to the totipotent embryo, which can subsequently give rise to all the cell types of the adult organism. In nearly all animals, maternal gene products regulate the initial events of embryogenesis while the zygotic genome remains transcriptionally silent. Developmental control is then passed from mother to zygote through a process known as the maternal-to-zygotic transition (MZT). The MZT comprises an intimately connected set of molecular events that mediate degradation of maternally deposited mRNAs and transcriptional activation of the zygotic genome. This essential developmental transition is conserved among metazoans but is perhaps best understood in the fruit fly, Drosophila melanogaster. In this article, we will review our understanding of the events that drive the MZT in Drosophila embryos and highlight parallel mechanisms driving this transition in other animals.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

An efficient method to transcription factor binding sites imputation via simultaneous completion of multiple matrices with positional consistency

Transcription factors (TFs) are DNA-binding proteins that have a central role in regulating gene expression. Identification of DNA-binding sites of TFs is a key task in understanding transcriptional regulation, cellular processes and disease. Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) enables genome-wide identification of in vivo TF binding sites. However, it is still difficult to map every TF in every cell line owing to cost and biological material availability, which poses an enormous obstacle for integrated analysis of gene regulation. To address this problem, we propose a novel computational approach, TFBSImpute, for predicting additional TF binding profiles by leveraging information from available ChIP-seq TF binding data. TFBSImpute fuses the dataset to a 3-mode tensor and imputes missing TF binding signals via simultaneous completion of multiple TF binding matrices with positional consistency. We show that signals predicted by our method achieve overall similarity with experimental data and that TFBSImpute significantly outperforms baseline approaches, by assessing the performance of imputation methods against observed ChIP-seq TF binding profiles. Besides, motif analysis shows that TFBSImpute preforms better in capturing binding motifs enriched in observed data compared with baselines, indicating that the higher performance of TFBSImpute is not simply due to averaging related samples. We anticipate that our approach will constitute a useful complement to experimental mapping of TF binding, which is beneficial for further study of regulation mechanisms and disease.

Integrated analysis of motif activity and gene expression changes of transcription factors

The ability to predict transcription factors based on sequence information in regulatory elements is a key step in systems-level investigation of transcriptional regulation. Here, we have developed a novel tool, IMAGE, for precise prediction of causal transcription factors based on transcriptome profiling and genome-wide maps of enhancer activity. High precision is obtained by combining a near-complete database of position weight matrices (PWMs), generated by compiling public databases and systematic prediction of PWMs for uncharacterized transcription factors, with a state-of-the-art method for PWM scoring and a novel machine learning strategy, based on both enhancers and promoters, to predict the contribution of motifs to transcriptional activity. We applied IMAGE to published data obtained during 3T3-L1 adipocyte differentiation and showed that IMAGE predicts causal transcriptional regulators of this process with higher confidence than existing methods. Furthermore, we generated genome-wide maps of enhancer activity and transcripts during human mesenchymal stem cell commitment and adipocyte differentiation and used IMAGE to identify positive and negative transcriptional regulators of this process. Collectively, our results demonstrate that IMAGE is a powerful and precise method for prediction of regulators of gene expression.

Local sequence features that influence AP-1 cis-regulatory activity

In the genome, most occurrences of transcription factor binding sites (TFBS) have no cis-regulatory activity, which suggests that flanking sequences contain information that distinguishes functional from nonfunctional TFBS. We interrogated the role of flanking sequences near Activator Protein 1 (AP-1) binding sites that reside in DNase I Hypersensitive Sites (DHS) and regions annotated as Enhancers. In these regions, we found that sequence features directly adjacent to the core motif distinguish high from low activity AP-1 sites. Some nearby features are motifs for other TFs that genetically interact with the AP-1 site. Other features are extensions of the AP-1 core motif, which cause the extended sites to match motifs of multiple AP-1 binding proteins. Computational models trained on these data distinguish between sequences with high and low activity AP-1 sites and also predict changes in cis-regulatory activity due to mutations in AP-1 core sites and their flanking sequences. Our results suggest that extended AP-1 binding sites, together with adjacent binding sites for additional TFs, encode part of the information that governs TFBS activity in the genome.

Uncoupling evolutionary changes in DNA sequence, transcription factor occupancy and enhancer activity

Sequence variation within enhancers plays a major role in both evolution and disease, yet its functional impact on transcription factor (TF) occupancy and enhancer activity remains poorly understood. Here, we assayed the binding of five essential TFs over multiple stages of embryogenesis in two distant Drosophila species (with 1.4 substitutions per neutral site), identifying thousands of orthologous enhancers with conserved or diverged combinatorial occupancy. We used these binding signatures to dissect two properties of developmental enhancers: (1) potential TF cooperativity, using signatures of co-associations and co-divergence in TF occupancy. This revealed conserved combinatorial binding despite sequence divergence, suggesting protein-protein interactions sustain conserved collective occupancy. (2) Enhancer in-vivo activity, revealing orthologous enhancers with conserved activity despite divergence in TF occupancy. Taken together, we identify enhancers with diverged motifs yet conserved occupancy and others with diverged occupancy yet conserved activity, emphasising the need to functionally measure the effect of divergence on enhancer activity.

Assessing the model transferability for prediction of transcription factor binding sites based on chromatin accessibility

Background

Computational prediction of transcription factor (TF) binding sites in different cell types is challenging. Recent technology development allows us to determine the genome-wide chromatin accessibility in various cellular and developmental contexts. The chromatin accessibility profiles provide useful information in prediction of TF binding events in various physiological conditions. Furthermore, ChIP-Seq analysis was used to determine genome-wide binding sites for a range of different TFs in multiple cell types. Integration of these two types of genomic information can improve the prediction of TF binding events.

Results

We assessed to what extent a model built upon on other TFs and/or other cell types could be used to predict the binding sites of TFs of interest. A random forest model was built using a set of cell type-independent features such as specific sequences recognized by the TFs and evolutionary conservation, as well as cell type-specific features derived from chromatin accessibility data. Our analysis suggested that the models learned from other TFs and/or cell lines performed almost as well as the model learned from the target TF in the cell type of interest. Interestingly, models based on multiple TFs performed better than single-TF models. Finally, we proposed a universal model, BPAC, which was generated using ChIP-Seq data from multiple TFs in various cell types.

Conclusion

Integrating chromatin accessibility information with sequence information improves prediction of TF binding.The prediction of TF binding is transferable across TFs and/or cell lines suggesting there are a set of universal “rules”. A computational tool was developed to predict TF binding sites based on the universal “rules”.

Electronic supplementary material

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

Short DNA sequence patterns accurately identify broadly active human enhancers

Background

Enhancers 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.

Results

We 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.

Conclusions

Our 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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3934-9) contains supplementary material, which is available to authorized users.

ConTra v3: a tool to identify transcription factor binding sites across species, update 2017

Transcription factors are important gene regulators with distinctive roles in development, cell signaling and cell cycling, and they have been associated with many diseases. The ConTra v3 web server allows easy visualization and exploration of predicted transcription factor binding sites (TFBSs) in any genomic region surrounding coding or non-coding genes. In this updated version, with a completely re-implemented user interface using latest web technologies, users can choose from nine reference organisms ranging from human to yeast. ConTra v3 can analyze promoter regions, 5΄-UTRs, 3΄-UTRs and introns or any other genomic region of interest. Thousands of position weight matrices are available to choose from for detecting specific binding sites. Besides this visualization option, additional new exploration functionality is added to the tool that will automatically detect TFBSs having at the same time the highest regulatory potential, the highest conservation scores of the genomic regions covered by the predicted TFBSs and strongest co-localizations with genomic regions exhibiting regulatory activity. The ConTra v3 web server is freely available at http://bioit2.irc.ugent.be/contra/v3.

The origins of developmental gene regulation

The leap from simple unicellularity to complex multicellularity remains one of life's major enigmas. The origins of metazoan developmental gene regulatory mechanisms are sought by analyzing gene regulation in extant eumetazoans, sponges, and unicellular organisms. The main hypothesis of this manuscript is that, developmental enhancers evolved from unicellular inducible promoters that diversified the expression of regulatory genes during metazoan evolution. Promoters and enhancers are functionally similar; both can regulate the transcription of distal promoters and both direct local transcription. Additionally, enhancers have experimentally characterized structural features that reveal their origin from inducible promoters. The distal co-operative regulation among promoters identified in unicellular opisthokonts possibly represents the precursor of distal regulation of promoters by enhancers. During metazoan evolution, constitutive-type promoters of regulatory genes would have acquired novel receptivity to distal regulatory inputs from promoters of inducible genes that eventually specialized as enhancers. The novel regulatory interactions would have caused constitutively expressed genes controlling differential gene expression in unicellular organisms to become themselves differentially expressed. The consequence of the novel regulatory interactions was that regulatory pathways of unicellular organisms became interlaced and ultimately evolved into the intricate developmental gene regulatory networks (GRNs) of extant metazoans.

Combinatorial function of transcription factors and cofactors

Differential gene expression gives rise to the many cell types of complex organisms. Enhancers regulate transcription by binding transcription factors (TFs), which in turn recruit cofactors to activate RNA Polymerase II at core promoters. Transcriptional regulation is typically mediated by distinct combinations of TFs, enabling a relatively small number of TFs to generate a large diversity of cell types. However, how TFs achieve combinatorial enhancer control and how enhancers, enhancer-bound TFs, and the cofactors they recruit regulate RNA Polymerase II activity is not entirely clear. Here, we review how TF synergy is mediated at the level of DNA binding and after binding, the role of cofactors and the post-translational modifications they catalyze, and discuss different models of enhancer-core-promoter communication.

Stable Binding of the Conserved Transcription Factor Grainy Head to its Target Genes Throughout Drosophila melanogaster Development

It has been suggested that transcription factor binding is temporally dynamic, and that changes in binding determine transcriptional output. Nonetheless, this model is based on relatively few examples in which transcription factor binding has been assayed at multiple developmental stages. The essential transcription factor Grainy head (Grh) is conserved from fungi to humans, and controls epithelial development and barrier formation in numerous tissues. Drosophila melanogaster, which possess a single grainy head (grh) gene, provide an excellent system to study this conserved factor. To determine whether temporally distinct binding events allow Grh to control cell fate specification in different tissue types, we used a combination of ChIP-seq and RNA-seq to elucidate the gene regulatory network controlled by Grh during four stages of embryonic development (spanning stages 5-17) and in larval tissue. Contrary to expectations, we discovered that Grh remains bound to at least 1146 genomic loci over days of development. In contrast to this stable DNA occupancy, the subset of genes whose expression is regulated by Grh varies. Grh transitions from functioning primarily as a transcriptional repressor early in development to functioning predominantly as an activator later. Our data reveal that Grh binds to target genes well before the Grh-dependent transcriptional program commences, suggesting it sets the stage for subsequent recruitment of additional factors that execute stage-specific Grh functions.

EP-DNN: A Deep Neural Network-Based Global Enhancer Prediction Algorithm

We present EP-DNN, a protocol for predicting enhancers based on chromatin features, in different cell types. Specifically, we use a deep neural network (DNN)-based architecture to extract enhancer signatures in a representative human embryonic stem cell type (H1) and a differentiated lung cell type (IMR90). We train EP-DNN using p300 binding sites, as enhancers, and TSS and random non-DHS sites, as non-enhancers. We perform same-cell and cross-cell predictions to quantify the validation rate and compare against two state-of-the-art methods, DEEP-ENCODE and RFECS. We find that EP-DNN has superior accuracy with a validation rate of 91.6%, relative to 85.3% for DEEP-ENCODE and 85.5% for RFECS, for a given number of enhancer predictions and also scales better for a larger number of enhancer predictions. Moreover, our H1 → IMR90 predictions turn out to be more accurate than IMR90 → IMR90, potentially because H1 exhibits a richer signature set and our EP-DNN model is expressive enough to extract these subtleties. Our work shows how to leverage the full expressivity of deep learning models, using multiple hidden layers, while avoiding overfitting on the training data. We also lay the foundation for exploration of cross-cell enhancer predictions, potentially reducing the need for expensive experimentation.

Drosophila poised enhancers are generated during tissue patterning with the help of repression

Histone modifications are frequently used as markers for enhancer states, but how to interpret enhancer states in the context of embryonic development is not clear. The poised enhancer signature, involving H3K4me1 and low levels of H3K27ac, has been reported to mark inactive enhancers that are poised for future activation. However, future activation is not always observed, and alternative reasons for the widespread occurrence of this enhancer signature have not been investigated. By analyzing enhancers during dorsal-ventral (DV) axis formation in the Drosophila embryo, we find that the poised enhancer signature is specifically generated during patterning in the tissue where the enhancers are not induced, including at enhancers that are known to be repressed by a transcriptional repressor. These results suggest that, rather than serving exclusively as an intermediate step before future activation, the poised enhancer state may be a mark for spatial regulation during tissue patterning. We discuss the possibility that the poised enhancer state is more generally the result of repression by transcriptional repressors.

Transcription factors involved in drought tolerance and their possible role in developing drought tolerant cultivars with emphasis on wheat (Triticum aestivum L.)

TFs involved in drought tolerance in plants may be utilized in future for developing drought tolerant cultivars of wheat and some other crops. Plants have developed a fairly complex stress response system to deal with drought and other abiotic stresses. These response systems often make use of transcription factors (TFs); a gene encoding a specific TF together with -its target genes constitute a regulon, and take part in signal transduction to activate/silence genes involved in response to drought. Since, five specific families of TFs (out of >80 known families of TFs) have gained widespread attention on account of their significant role in drought tolerance in plants, TFs and regulons belonging to these five multi-gene families (AP2/EREBP, bZIP, MYB/MYC, NAC and WRKY) have been described and their role in improving drought tolerance discussed in this brief review. These TFs often undergo reversible phosphorylation to perform their function, and are also involved in complex networks. Therefore, some details about reversible phosphorylation of TFs by different protein kinases/phosphatases and the co-regulatory networks, which involve either only TFs or TFs with miRNAs, have also been discussed. Literature on transgenics involving genes encoding TFs and that on QTLs and markers associated with TF genes involved in drought tolerance has also been reviewed. Throughout the review, there is a major emphasis on wheat as an important crop, although examples from the model cereal rice (sometimes maize also), and the model plant Arabidopsis have also been used. This knowledge base may eventually allow the use of TF genes for development of drought tolerant cultivars, particularly in wheat.

An Intronic Enhancer Element Regulates Angiotensin II Type 2 Receptor Expression during Satellite Cell Differentiation, and Its Activity Is Suppressed in Congestive Heart Failure

Patients with advanced congestive heart failure (CHF) or chronic kidney disease often have increased angiotensin II (Ang II) levels and cachexia. We previously demonstrated that Ang II, via its type 1 receptor, causes muscle protein breakdown and apoptosis and inhibits satellite cell (SC) proliferation and muscle regeneration, likely contributing to cachexia in CHF and chronic kidney disease. In contrast, Ang II, via its type 2 receptor (AT2R) expression, is robustly induced during SC differentiation, and it potentiates muscle regeneration. To understand the mechanisms regulating AT2R expression and its potential role in muscle regeneration in chronic diseases, we used a mouse model of CHF and found that muscle regeneration was markedly reduced and that this was accompanied by blunted increase of AT2R expression. We performed AT2R promoter reporter analysis during satellite cell differentiation and found that the 70 bp upstream of the AT2R transcription start site contain a core promoter region, and regions upstream of 70 bp to 3 kbp are dispensable for AT2R induction. Instead, AT2R intron 2 acts as a transcriptional enhancer during SC differentiation. Further deletion/mutation analysis revealed that multiple transcription factor binding sites in the +286/+690 region within intron 2 coordinately regulate AT2R transcription. Importantly, +286/+690 enhancer activity was suppressed in CHF mouse skeletal muscle, suggesting that AT2R expression is suppressed in CHF via inhibition of AT2R intronic enhancer activity, leading to lowered muscle regeneration. Thus targeting intron 2 enhancer element could lead to the development of a novel intervention to increase AT2R expression in SCs and potentiate skeletal muscle regenerative capacity in chronic diseases.

Genome-Wide Ultrabithorax Binding Analysis Reveals Highly Targeted Genomic Loci at Developmental Regulators and a Potential Connection to Polycomb-Mediated Regulation

Hox homeodomain transcription factors are key regulators of animal development. They specify the identity of segments along the anterior-posterior body axis in metazoans by controlling the expression of diverse downstream targets, including transcription factors and signaling pathway components. The Drosophila melanogaster Hox factor Ultrabithorax (Ubx) directs the development of thoracic and abdominal segments and appendages, and loss of Ubx function can lead for example to the transformation of third thoracic segment appendages (e.g. halters) into second thoracic segment appendages (e.g. wings), resulting in a characteristic four-wing phenotype. Here we present a Drosophila melanogaster strain with a V5-epitope tagged Ubx allele, which we employed to obtain a high quality genome-wide map of Ubx binding sites using ChIP-seq. We confirm the sensitivity of the V5 ChIP-seq by recovering 7/8 of well-studied Ubx-dependent cis-regulatory regions. Moreover, we show that Ubx binding is predictive of enhancer activity as suggested by comparison with a genome-scale resource of in vivo tested enhancer candidates. We observed densely clustered Ubx binding sites at 12 extended genomic loci that included ANTP-C, BX-C, Polycomb complex genes, and other regulators and the clustered binding sites were frequently active enhancers. Furthermore, Ubx binding was detected at known Polycomb response elements (PREs) and was associated with significant enrichments of Pc and Pho ChIP signals in contrast to binding sites of other developmental TFs. Together, our results show that Ubx targets developmental regulators via strongly clustered binding sites and allow us to hypothesize that regulation by Ubx might involve Polycomb group proteins to maintain specific regulatory states in cooperative or mutually exclusive fashion, an attractive model that combines two groups of proteins with prominent gene regulatory roles during animal development.

Opening up the blackbox: an interpretable deep neural network-based classifier for cell-type specific enhancer predictions

Background

Gene expression is mediated by specialized cis-regulatory modules (CRMs), the most prominent of which are called enhancers. Early experiments indicated that enhancers located far from the gene promoters are often responsible for mediating gene transcription. Knowing their properties, regulatory activity, and genomic targets is crucial to the functional understanding of cellular events, ranging from cellular homeostasis to differentiation. Recent genome-wide investigation of epigenomic marks has indicated that enhancer elements could be enriched for certain epigenomic marks, such as, combinatorial patterns of histone modifications.

Methods

Our efforts in this paper are motivated by these recent advances in epigenomic profiling methods, which have uncovered enhancer-associated chromatin features in different cell types and organisms. Specifically, in this paper, we use recent state-of-the-art Deep Learning methods and develop a deep neural network (DNN)-based architecture, called EP-DNN, to predict the presence and types of enhancers in the human genome. It uses as features, the expression levels of the histone modifications at the peaks of the functional sites as well as in its adjacent regions. We apply EP-DNN to four different cell types: H1, IMR90, HepG2, and HeLa S3. We train EP-DNN using p300 binding sites as enhancers, and TSS and random non-DHS sites as non-enhancers. We perform EP-DNN predictions to quantify the validation rate for different levels of confidence in the predictions and also perform comparisons against two state-of-the-art computational models for enhancer predictions, DEEP-ENCODE and RFECS.

Results

We find that EP-DNN has superior accuracy and takes less time to make predictions. Next, we develop methods to make EP-DNN interpretable by computing the importance of each input feature in the classification task. This analysis indicates that the important histone modifications were distinct for different cell types, with some overlaps, e.g., H3K27ac was important in cell type H1 but less so in HeLa S3, while H3K4me1 was relatively important in all four cell types. We finally use the feature importance analysis to reduce the number of input features needed to train the DNN, thus reducing training time, which is often the computational bottleneck in the use of a DNN.

Conclusions

In this paper, we developed EP-DNN, which has high accuracy of prediction, with validation rates above 90 % for the operational region of enhancer prediction for all four cell lines that we studied, outperforming DEEP-ENCODE and RFECS. Then, we developed a method to analyze a trained DNN and determine which histone modifications are important, and within that, which features proximal or distal to the enhancer site, are important.