Dual-functioning transcription factors in the developmental gene network of Drosophila melanogaster

scANANSE gene regulatory network and motif analysis of single-cell clusters

The recent development of single-cell techniques is essential to unravel complex biological systems. By measuring the transcriptome and the accessible genome on a single-cell level, cellular heterogeneity in a biological environment can be deciphered. Transcription factors act as key regulators activating and repressing downstream target genes, and together they constitute gene regulatory networks that govern cell morphology and identity. Dissecting these gene regulatory networks is crucial for understanding molecular mechanisms and disease, especially within highly complex biological systems. The gene regulatory network analysis software ANANSE and the motif enrichment software GimmeMotifs were both developed to analyse bulk datasets. We developed scANANSE, a software pipeline for gene regulatory network analysis and motif enrichment using single-cell RNA and ATAC datasets. The scANANSE pipeline can be run from either R or Python. First, it exports data from standard single-cell objects. Next, it automatically runs multiple comparisons of cell cluster data. Finally, it imports the results back to the single-cell object, where the result can be further visualised, integrated, and interpreted. Here, we demonstrate our scANANSE pipeline on a publicly available PBMC multi-omics dataset. It identifies well-known cell type-specific hematopoietic factors. Importantly, we also demonstrated that scANANSE combined with GimmeMotifs is able to predict transcription factors with both activating and repressing roles in gene regulation.

Thermodynamics-based modeling reveals regulatory effects of indirect transcription factor-DNA binding

Transcription factors (TFs) influence gene expression by binding to DNA, yet experimental data suggests that they also frequently bind regulatory DNA indirectly by interacting with other DNA-bound proteins. Here, we used a data modeling approach to test if such indirect binding by TFs plays a significant role in gene regulation. We first incorporated regulatory function of indirectly bound TFs into a thermodynamics-based model for predicting enhancer-driven expression from its sequence. We then fit the new model to a rich data set comprising hundreds of enhancers and their regulatory activities during mesoderm specification in Drosophila embryogenesis and showed that the newly incorporated mechanism results in significantly better agreement with data. In the process, we derived the first sequence-level model of this extensively characterized regulatory program. We further showed that allowing indirect binding of a TF explains its localization at enhancers more accurately than with direct binding only. Our model also provided a simple explanation of how a TF may switch between activating and repressive roles depending on context.

•

Inclusion of indirect DNA binding of transcription factor improves enhancer function prediction

•

Context specific activating or repressive roles of TFs

•

Indirect binding improves fits to experimental TF-DNA binding data

•

Role of Tinman depends on its DNA-binding mode (direct or indirect)

Epigenetics, Enhancer Function and 3D Chromatin Organization in Reprogramming to Pluripotency

Genome architecture, epigenetics and enhancer function control the fate and identity of cells. Reprogramming to induced pluripotent stem cells (iPSCs) changes the transcriptional profile and chromatin landscape of the starting somatic cell to that of the pluripotent cell in a stepwise manner. Changes in the regulatory networks are tightly regulated during normal embryonic development to determine cell fate, and similarly need to function in cell fate control during reprogramming. Switching off the somatic program and turning on the pluripotent program involves a dynamic reorganization of the epigenetic landscape, enhancer function, chromatin accessibility and 3D chromatin topology. Within this context, we will review here the current knowledge on the processes that control the establishment and maintenance of pluripotency during somatic cell reprogramming.

Targeting the Transcriptome Through Globally Acting Components

Transcription is a step in gene expression that defines the identity of cells and its dysregulation is associated with diseases. With advancing technologies revealing molecular underpinnings of the cell with ever-higher precision, our ability to view the transcriptomes may have surpassed our knowledge of the principles behind their organization. The human RNA polymerase II (Pol II) machinery comprises thousands of components that, in conjunction with epigenetic and other mechanisms, drive specialized programs of development, differentiation, and responses to the environment. Parts of these programs are repurposed in oncogenic transformation. Targeting of cancers is commonly done by inhibiting general or broadly acting components of the cellular machinery. The critical unanswered question is how globally acting or general factors exert cell type specific effects on transcription. One solution, which is discussed here, may be among the events that take place at genes during early Pol II transcription elongation. This essay turns the spotlight on the well-known phenomenon of promoter-proximal Pol II pausing as a step that separates signals that establish pausing genome-wide from those that release the paused Pol II into the gene. Concepts generated in this rapidly developing field will enhance our understanding of basic principles behind transcriptome organization and hopefully translate into better therapies at the bedside.

The DREAM complex represses growth in response to DNA damage in Arabidopsis

The DNA of all organisms is constantly damaged by physiological processes and environmental conditions. Upon persistent damage, plant growth and cell proliferation are reduced. Based on previous findings that RBR1, the only Arabidopsis homolog of the mammalian tumor suppressor gene retinoblastoma, plays a key role in the DNA damage response in plants, we unravel here the network of RBR1 interactors under DNA stress conditions. This led to the identification of homologs of every DREAM component in Arabidopsis, including previously not recognized homologs of LIN52. Interestingly, we also discovered NAC044, a mediator of DNA damage response in plants and close homolog of the major DNA damage regulator SOG1, to directly interact with RBR1 and the DREAM component LIN37B. Consistently, not only mutants in NAC044 but also the double mutant of the two LIN37 homologs and mutants for the DREAM component E2FB showed reduced sensitivities to DNA-damaging conditions. Our work indicates the existence of multiple DREAM complexes that work in conjunction with NAC044 to mediate growth arrest after DNA damage.

Collaborative repressive action of the antagonistic ETS transcription factors Pointed and Yan fine-tunes gene expression to confer robustness in Drosophila

The acquisition of cellular identity during development depends on precise spatiotemporal regulation of gene expression, with combinatorial interactions between transcription factors, accessory proteins and the basal transcription machinery together translating complex signaling inputs into appropriate gene expression outputs. The opposing repressive and activating inputs of the Drosophila ETS family transcription factors Yan and Pointed orchestrate numerous cell fate transitions downstream of receptor tyrosine kinase signaling, providing one of the premier systems for studying this process. Current models describe the differentiative transition as a switch from Yan-mediated repression to Pointed-mediated activation of common target genes. We describe here a new layer of regulation whereby Yan and Pointed co-occupy regulatory elements to repress gene expression in a coordinated manner, with Pointed being unexpectedly required for the genome-wide occupancy of both Yan and the co-repressor Groucho. Using even skipped as a test-case, synergistic genetic interactions between Pointed, Groucho, Yan and components of the RNA polymerase II pausing machinery suggest that Pointed integrates multiple scales of repressive regulation to confer robustness. We speculate that this mechanism may be used broadly to fine-tune the expression of many genes crucial for development.

Profiling of Accessible Chromatin Regions across Multiple Plant Species and Cell Types Reveals Common Gene Regulatory Principles and New Control Modules

The transcriptional regulatory structure of plant genomes remains poorly defined relative to animals. It is unclear how many cis-regulatory elements exist, where these elements lie relative to promoters, and how these features are conserved across plant species. We employed the assay for transposase-accessible chromatin (ATAC-seq) in four plant species (Arabidopsis thaliana, Medicago truncatula, Solanum lycopersicum, and Oryza sativa) to delineate open chromatin regions and transcription factor (TF) binding sites across each genome. Despite 10-fold variation in intergenic space among species, the majority of open chromatin regions lie within 3 kb upstream of a transcription start site in all species. We find a common set of four TFs that appear to regulate conserved gene sets in the root tips of all four species, suggesting that TF-gene networks are generally conserved. Comparative ATAC-seq profiling of Arabidopsis root hair and non-hair cell types revealed extensive similarity as well as many cell-type-specific differences. Analyzing TF binding sites in differentially accessible regions identified a MYB-driven regulatory module unique to the hair cell, which appears to control both cell fate regulators and abiotic stress responses. Our analyses revealed common regulatory principles among species and shed light on the mechanisms producing cell-type-specific transcriptomes during development.

Profiling of Accessible Chromatin Regions across Multiple Plant Species and Cell Types Reveals Common Gene Regulatory Principles and New Control Modules

The transcriptional regulatory structure of plant genomes remains poorly defined relative to animals. It is unclear how many cis-regulatory elements exist, where these elements lie relative to promoters, and how these features are conserved across plant species. We employed the assay for transposase-accessible chromatin (ATAC-seq) in four plant species (Arabidopsis thaliana, Medicago truncatula, Solanum lycopersicum, and Oryza sativa) to delineate open chromatin regions and transcription factor (TF) binding sites across each genome. Despite 10-fold variation in intergenic space among species, the majority of open chromatin regions lie within 3 kb upstream of a transcription start site in all species. We find a common set of four TFs that appear to regulate conserved gene sets in the root tips of all four species, suggesting that TF-gene networks are generally conserved. Comparative ATAC-seq profiling of Arabidopsis root hair and non-hair cell types revealed extensive similarity as well as many cell-type-specific differences. Analyzing TF binding sites in differentially accessible regions identified a MYB-driven regulatory module unique to the hair cell, which appears to control both cell fate regulators and abiotic stress responses. Our analyses revealed common regulatory principles among species and shed light on the mechanisms producing cell-type-specific transcriptomes during development.

Ehrlichia chaffeensis TRP32 is a Nucleomodulin that Directly Regulates Expression of Host Genes Governing Differentiation and Proliferation

Ehrlichia chaffeensis is an obligately intracellular bacterium that reprograms the mononuclear phagocyte through diverse effector-host interactions to modulate numerous host cell processes, including transcription. In a previous study, we reported that E. chaffeensis TRP32, a type 1 secreted effector, interacts with multiple host nucleus-associated proteins and also auto-activates reporter gene expression in yeast. In this study, we demonstrate that TRP32 is a nucleomodulin that binds host DNA and alters host gene transcription. TRP32 enters the host cell nucleus via a noncanonical translocation mechanism that involves phosphorylation of Y179 located in a C-terminal tri-tyrosine motif. Both genistein and mutation of Y179 inhibited TRP32 nuclear entry. An electromobility shift assay (EMSA) demonstrated TRP32 host DNA binding via its tandem repeat domain. TRP32 DNA binding and motif preference were further confirmed by supershift assays, as well as competition and mutant probe analyses. Using ChIP-Seq, we determined that TRP32 binds a G-rich motif primarily within ±500 bp of the gene transcription start site. An ontology analysis identified genes involved in processes such as immune cell differentiation, chromatin remodeling, and RNA transcription and processing, as primary TRP32 targets. TRP32 bound genes (n=1223) were distributed on all chromosomes and included several global regulators of proliferation and inflammation such as FOS and JUN, AKT3 and NRAS, and non-coding RNA genes, miRNA 21 and miRNA 142. TRP32 target genes were differentially regulated during infection, the majority of which were repressed, and direct repression/activation of these genes by TRP32 was confirmed in vitro with a cellular luciferase reporter assay.

Mid-embryo patterning and precision in Drosophila segmentation: Krüppel dual regulation of hunchback

In early development, genes are expressed in spatial patterns which later define cellular identities and tissue locations. The mechanisms of such pattern formation have been studied extensively in early Drosophila (fruit fly) embryos. The gap gene hunchback (hb) is one of the earliest genes to be expressed in anterior-posterior (AP) body segmentation. As a transcriptional regulator for a number of downstream genes, the spatial precision of hb expression can have significant effects in the development of the body plan. To investigate the factors contributing to hb precision, we used fine spatial and temporal resolution data to develop a quantitative model for the regulation of hb expression in the mid-embryo. In particular, modelling hb pattern refinement in mid nuclear cleavage cycle 14 (NC14) reveals some of the regulatory contributions of simultaneously-expressed gap genes. Matching the model to recent data from wild-type (WT) embryos and mutants of the gap gene Krüppel (Kr) indicates that a mid-embryo Hb concentration peak important in thoracic development (at parasegment 4, PS4) is regulated in a dual manner by Kr, with low Kr concentration activating hb and high Kr concentration repressing hb. The processes of gene expression (transcription, translation, transport) are intrinsically random. We used stochastic simulations to characterize the noise generated in hb expression. We find that Kr regulation can limit the positional variability of the Hb mid-embryo border. This has been recently corroborated in experimental comparisons of WT and Kr- mutant embryos. Further, Kr regulation can decrease uncertainty in mid-embryo hb expression (i.e. contribute to a smooth Hb boundary) and decrease between-copy transcriptional variability within nuclei. Since many tissue boundaries are first established by interactions between neighbouring gene expression domains, these properties of Hb-Kr dynamics to diminish the effects of intrinsic expression noise may represent a general mechanism contributing to robustness in early development.

Dissection of thousands of cell type-specific enhancers identifies dinucleotide repeat motifs as general enhancer features

Gene expression is determined by genomic elements called enhancers, which contain short motifs bound by different transcription factors (TFs). However, how enhancer sequences and TF motifs relate to enhancer activity is unknown, and general sequence requirements for enhancers or comprehensive sets of important enhancer sequence elements have remained elusive. Here, we computationally dissect thousands of functional enhancer sequences from three different Drosophila cell lines. We find that the enhancers display distinct cis-regulatory sequence signatures, which are predictive of the enhancers’ cell type-specific or broad activities. These signatures contain transcription factor motifs and a novel class of enhancer sequence elements, dinucleotide repeat motifs (DRMs). DRMs are highly enriched in enhancers, particularly in enhancers that are broadly active across different cell types. We experimentally validate the importance of the identified TF motifs and DRMs for enhancer function and show that they can be sufficient to create an active enhancer de novo from a nonfunctional sequence. The function of DRMs as a novel class of general enhancer features that are also enriched in human regulatory regions might explain their implication in several diseases and provides important insights into gene regulation.

ISMARA: automated modeling of genomic signals as a democracy of regulatory motifs

Accurate reconstruction of the regulatory networks that control gene expression is one of the key current challenges in molecular biology. Although gene expression and chromatin state dynamics are ultimately encoded by constellations of binding sites recognized by regulators such as transcriptions factors (TFs) and microRNAs (miRNAs), our understanding of this regulatory code and its context-dependent read-out remains very limited. Given that there are thousands of potential regulators in mammals, it is not practical to use direct experimentation to identify which of these play a key role for a particular system of interest. We developed a methodology that models gene expression or chromatin modifications in terms of genome-wide predictions of regulatory sites and completely automated it into a web-based tool called ISMARA (Integrated System for Motif Activity Response Analysis). Given only gene expression or chromatin state data across a set of samples as input, ISMARA identifies the key TFs and miRNAs driving expression/chromatin changes and makes detailed predictions regarding their regulatory roles. These include predicted activities of the regulators across the samples, their genome-wide targets, enriched gene categories among the targets, and direct interactions between the regulators. Applying ISMARA to data sets from well-studied systems, we show that it consistently identifies known key regulators ab initio. We also present a number of novel predictions including regulatory interactions in innate immunity, a master regulator of mucociliary differentiation, TFs consistently disregulated in cancer, and TFs that mediate specific chromatin modifications.

Global SUMO Proteome Responses Guide Gene Regulation, mRNA Biogenesis, and Plant Stress Responses

Small Ubiquitin-like MOdifier (SUMO) is a key regulator of abiotic stress, disease resistance, and development in plants. The identification of >350 plant SUMO targets has revealed many processes modulated by SUMO and potential consequences of SUMO on its targets. Importantly, highly related proteins are SUMO-modified in plants, yeast, and metazoans. Overlapping SUMO targets include heat-shock proteins (HSPs), transcription regulators, histones, histone-modifying enzymes, proteins involved in DNA damage repair, but also proteins involved in mRNA biogenesis and nucleo-cytoplasmic transport. Proteomics studies indicate key roles for SUMO in gene repression by controlling histone (de)acetylation activity at genomic loci. The responsible heavily sumoylated transcriptional repressor complexes are recruited by plant transcription factors (TFs) containing an (ERF)-associated Amphiphilic Repression (EAR) motif. These TFs are not necessarily themselves a SUMO target. Conversely, SUMO acetylation (Ac) prevents binding of downstream partners by blocking binding of their SUMO-interaction peptide motifs to Ac-SUMO. In addition, SUMO acetylation has emerged as a mechanism to recruit specifically bromodomains. Bromodomains are generally linked with gene activation. These findings strengthen the idea of a bi-directional sumo-acetylation switch in gene regulation. Quantitative proteomics has highlighted that global sumoylation provides a dynamic response to protein damage involving SUMO chain-mediated protein degradation, but also SUMO E3 ligase-dependent transcription of HSP genes. With these insights in SUMO function and novel technical advancements, we can now study SUMO dynamics in responses to (a)biotic stress in plants.

Global SUMO Proteome Responses Guide Gene Regulation, mRNA Biogenesis, and Plant Stress Responses

Small Ubiquitin-like MOdifier (SUMO) is a key regulator of abiotic stress, disease resistance, and development in plants. The identification of >350 plant SUMO targets has revealed many processes modulated by SUMO and potential consequences of SUMO on its targets. Importantly, highly related proteins are SUMO-modified in plants, yeast, and metazoans. Overlapping SUMO targets include heat-shock proteins (HSPs), transcription regulators, histones, histone-modifying enzymes, proteins involved in DNA damage repair, but also proteins involved in mRNA biogenesis and nucleo-cytoplasmic transport. Proteomics studies indicate key roles for SUMO in gene repression by controlling histone (de)acetylation activity at genomic loci. The responsible heavily sumoylated transcriptional repressor complexes are recruited by plant transcription factors (TFs) containing an (ERF)-associated Amphiphilic Repression (EAR) motif. These TFs are not necessarily themselves a SUMO target. Conversely, SUMO acetylation (Ac) prevents binding of downstream partners by blocking binding of their SUMO-interaction peptide motifs to Ac-SUMO. In addition, SUMO acetylation has emerged as a mechanism to recruit specifically bromodomains. Bromodomains are generally linked with gene activation. These findings strengthen the idea of a bi-directional sumo-acetylation switch in gene regulation. Quantitative proteomics has highlighted that global sumoylation provides a dynamic response to protein damage involving SUMO chain-mediated protein degradation, but also SUMO E3 ligase-dependent transcription of HSP genes. With these insights in SUMO function and novel technical advancements, we can now study SUMO dynamics in responses to (a)biotic stress in plants.

Tramtrack is genetically upstream of genes controlling tracheal tube size in Drosophila

The Drosophila transcription factor Tramtrack (Ttk) is involved in a wide range of developmental decisions, ranging from early embryonic patterning to differentiation processes in organogenesis. Given the wide spectrum of functions and pleiotropic effects that hinder a comprehensive characterisation, many of the tissue specific functions of this transcription factor are only poorly understood. We recently discovered multiple roles of Ttk in the development of the tracheal system on the morphogenetic level. Here, we sought to identify some of the underlying genetic components that are responsible for the tracheal phenotypes of Ttk mutants. We therefore profiled gene expression changes after Ttk loss- and gain-of-function in whole embryos and cell populations enriched for tracheal cells. The analysis of the transcriptomes revealed widespread changes in gene expression. Interestingly, one of the most prominent gene classes that showed significant opposing responses to loss- and gain-of-function was annotated with functions in chitin metabolism, along with additional genes that are linked to cellular responses, which are impaired in ttk mutants. The expression changes of these genes were validated by quantitative real-time PCR and further functional analysis of these candidate genes and other genes also expected to control tracheal tube size revealed at least a partial explanation of Ttk's role in tube size regulation. The computational analysis of our tissue-specific gene expression data highlighted the sensitivity of the approach and revealed an interesting set of novel putatively tracheal genes.

Transcriptional activators and activation mechanisms

Transcriptional activators are required to turn on the expression of genes in a eukaryotic cell. Activators bound to the enhancer can facilitate either the recruitment of RNA polymerase II to the promoter or its elongation. This article examines a few selected issues in understanding activator functions and activation mechanisms.

Complex SUMO-1 regulation of cardiac transcription factor Nkx2-5

Reversible post-translational protein modifications such as SUMOylation add complexity to cardiac transcriptional regulation. The homeodomain transcription factor Nkx2-5/Csx is essential for heart specification and morphogenesis. It has been previously suggested that SUMOylation of lysine 51 (K51) of Nkx2-5 is essential for its DNA binding and transcriptional activation. Here, we confirm that SUMOylation strongly enhances Nkx2-5 transcriptional activity and that residue K51 of Nkx2-5 is a SUMOylation target. However, in a range of cultured cell lines we find that a point mutation of K51 to arginine (K51R) does not affect Nkx2-5 activity or DNA binding, suggesting the existence of additional Nkx2-5 SUMOylated residues. Using biochemical assays, we demonstrate that Nkx2-5 is SUMOylated on at least one additional site, and this is the predominant site in cardiac cells. The second site is either non-canonical or a "shifting" site, as mutation of predicted consensus sites and indeed every individual lysine in the context of the K51R mutation failed to impair Nkx2-5 transcriptional synergism with SUMO, or its nuclear localization and DNA binding. We also observe SUMOylation of Nkx2-5 cofactors, which may be critical to Nkx2-5 regulation. Our data reveal highly complex regulatory mechanisms driven by SUMOylation to modulate Nkx2-5 activity.