A Temporal Map of Transcription Factor Activity: Mef2 Directly Regulates Target Genes at All Stages of Muscle Development

A matter of time - How transient transcription factor interactions create dynamic gene regulatory networks

Dynamic reprogramming of transcriptional networks enables cells to adapt to a changing environment. Thus, it is crucial not only to understand what gene targets are regulated by a transcription factor (TF) but also when. This review explores the way TFs function with respect to time, paying particular attention to discoveries made in plants - where coordinated, genome-wide responses to environmental change is crucial to the survival of these sessile organisms. We investigate the molecular mechanisms that mediate transient TF-DNA binding, and assess how these rapid and dynamic interactions translate to long-term temporal regulation of genomes. We also discuss how current molecular techniques can catch, and sometimes miss, transient TF-target interactions that underlie dynamic cellular responses. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

The Function of the MEF2 Family of Transcription Factors in Cardiac Development, Cardiogenomics, and Direct Reprogramming

Proper formation of the mammalian heart requires precise spatiotemporal transcriptional regulation of gene programs in cardiomyocytes. Sophisticated regulatory networks have evolved to not only integrate the activities of distinct transcription factors to control tissue-specific gene programs but also, in many instances, to incorporate multiple members within these transcription factor families to ensure accuracy and specificity in the system. Unsurprisingly, perturbations in this elaborate transcriptional circuitry can lead to severe cardiac abnormalities. Myocyte enhancer factor–2 (MEF2) transcription factor belongs to the evolutionarily conserved cardiac gene regulatory network. Given its central role in muscle gene regulation and its evolutionary conservation, MEF2 is considered one of only a few core cardiac transcription factors. In addition to its firmly established role as a differentiation factor, MEF2 regulates wide variety of, sometimes antagonistic, cellular processes such as cell survival and death. Vertebrate genomes encode multiple MEF2 family members thereby expanding the transcriptional potential of this core transcription factor in the heart. This review highlights the requirement of the MEF2 family and their orthologs in cardiac development in diverse animal model systems. Furthermore, we describe the recently characterized role of MEF2 in direct reprogramming and genome-wide cardiomyocyte gene regulation. A thorough understanding of the regulatory functions of the MEF2 family in cardiac development and cardiogenomics is required in order to develop effective therapeutic strategies to repair the diseased heart.

Genome-Wide Chromatin Immunoprecipitation Sequencing Analysis of the Penicillium chrysogenum Velvet Protein PcVelA Identifies Methyltransferase PcLlmA as a Novel Downstream Regulator of Fungal Development

Penicillium chrysogenum is the sole industrial producer of the β-lactam antibiotic penicillin, which is the most commonly used drug for treating bacterial infections. In P. chrysogenum and other filamentous fungi, secondary metabolism and morphogenesis are controlled by the highly conserved multisubunit velvet complex. Here we present the first chromatin immunoprecipitation next-generation sequencing (ChIP-seq) analysis of a fungal velvet protein, providing experimental evidence that a velvet homologue in P. chrysogenum (PcVelA) acts as a direct transcriptional regulator at the DNA level in addition to functioning as a regulator at the protein level in P. chrysogenum, which was previously described. We identified many target genes that are related to processes known to be dependent on PcVelA, e.g., secondary metabolism as well as asexual and sexual development. We also identified seven PcVelA target genes that encode putative methyltransferases. Yeast two-hybrid and bimolecular fluorescence complementation analyses showed that one of the putative methyltransferases, PcLlmA, directly interacts with PcVelA. Furthermore, functional characterization of PcLlmA demonstrated that this protein is involved in the regulation of conidiosporogenesis, pellet formation, and hyphal morphology, all traits with major biotechnological relevance. IMPORTANCE Filamentous fungi are of major interest for biotechnological and pharmaceutical applications. This is due mainly to their ability to produce a wide variety of secondary metabolites, many of which are relevant as antibiotics. One of the most prominent examples is penicillin, a β-lactam antibiotic that is produced on the industrial scale by fermentation of P. chrysogenum. In recent years, the multisubunit protein complex velvet has been identified as one of the key regulators of fungal secondary metabolism and development. However, until recently, only a little has been known about how velvet mediates regulation at the molecular level. To address this issue, we performed ChIP-seq (chromatin immunoprecipitation in combination with next-generation sequencing) on and follow-up analysis of PcVelA, the core component of the velvet complex in P. chrysogenum. We demonstrate direct involvement of velvet in transcriptional control and present the putative methyltransferase PcLlmA as a new downstream factor and interaction partner of PcVelA.

Loss of Frataxin induces iron toxicity, sphingolipid synthesis, and Pdk1/Mef2 activation, leading to neurodegeneration

Mutations in Frataxin (FXN) cause Friedreich's ataxia (FRDA), a recessive neurodegenerative disorder. Previous studies have proposed that loss of FXN causes mitochondrial dysfunction, which triggers elevated reactive oxygen species (ROS) and leads to the demise of neurons. Here we describe a ROS independent mechanism that contributes to neurodegeneration in fly FXN mutants. We show that loss of frataxin homolog (fh) in Drosophila leads to iron toxicity, which in turn induces sphingolipid synthesis and ectopically activates 3-phosphoinositide dependent protein kinase-1 (Pdk1) and myocyte enhancer factor-2 (Mef2). Dampening iron toxicity, inhibiting sphingolipid synthesis by Myriocin, or reducing Pdk1 or Mef2 levels, all effectively suppress neurodegeneration in fh mutants. Moreover, increasing dihydrosphingosine activates Mef2 activity through PDK1 in mammalian neuronal cell line suggesting that the mechanisms are evolutionarily conserved. Our results indicate that an iron/sphingolipid/Pdk1/Mef2 pathway may play a role in FRDA.

Hidden Markov induced Dynamic Bayesian Network for recovering time evolving gene regulatory networks

Dynamic Bayesian Networks (DBN) have been widely used to recover gene regulatory relationships from time-series data in computational systems biology. Its standard assumption is ‘stationarity’, and therefore, several research efforts have been recently proposed to relax this restriction. However, those methods suffer from three challenges: long running time, low accuracy and reliance on parameter settings. To address these problems, we propose a novel non-stationary DBN model by extending each hidden node of Hidden Markov Model into a DBN (called HMDBN), which properly handles the underlying time-evolving networks. Correspondingly, an improved structural EM algorithm is proposed to learn the HMDBN. It dramatically reduces searching space, thereby substantially improving computational efficiency. Additionally, we derived a novel generalized Bayesian Information Criterion under the non-stationary assumption (called BWBIC), which can help significantly improve the reconstruction accuracy and largely reduce over-fitting. Moreover, the re-estimation formulas for all parameters of our model are derived, enabling us to avoid reliance on parameter settings. Compared to the state-of-the-art methods, the experimental evaluation of our proposed method on both synthetic and real biological data demonstrates more stably high prediction accuracy and significantly improved computation efficiency, even with no prior knowledge and parameter settings.

Expression of the Troponin C at 41C Gene in Adult Drosophila Tubular Muscles Depends upon Both Positive and Negative Regulatory Inputs

Most animals express multiple isoforms of structural muscle proteins to produce tissues with different physiological properties. In Drosophila, the adult muscles include tubular-type muscles and the fibrillar indirect flight muscles. Regulatory processes specifying tubular muscle fate remain incompletely understood, therefore we chose to analyze the transcriptional regulation of TpnC41C, a Troponin C gene expressed in the tubular jump muscles, but not in the fibrillar flight muscles. We identified a 300-bp promoter fragment of TpnC41C sufficient for the fiber-specific reporter expression. Through an analysis of this regulatory element, we identified two sites necessary for the activation of the enhancer. Mutations in each of these sites resulted in 70% reduction of enhancer activity. One site was characterized as a binding site for Myocyte Enhancer Factor-2. In addition, we identified a repressive element that prevents activation of the enhancer in other muscle fiber types. Mutation of this site increased jump muscle-specific expression of the reporter, but more importantly reporter expression expanded into the indirect flight muscles. Our findings demonstrate that expression of the TpnC41C gene in jump muscles requires integration of multiple positive and negative transcriptional inputs. Identification of the transcriptional regulators binding the cis-elements that we identified will reveal the regulatory pathways controlling muscle fiber differentiation.

Immune-metabolic interaction in Drosophila

Immune responses and metabolic regulation are tightly coupled in all animals, but the underlying mechanistic connections are nowhere completely clear. In flies and in humans, prolonged or excessive immune activation can drive metabolic disruption and cause loss of metabolic stores. Conversely, disruptions of metabolic homeostasis, such as periods of malnutrition, can have significant impacts on immune function. We have recently identified the transcription factor MEF2 as a critical switch between anabolic and immune function in the adult Drosophila fat body. A conserved phosphorylation determines the affinity of MEF2 for the TATA-binding protein, effecting a choice between energy storage and immune function. The goal of this review is to place this molecular event in the broader context of metabolic-immune interaction in Drosophila, exploring what is and is not known about the ties between these 2 critical physiological functions.

Transcription Factor Activity Mapping of a Tissue-Specific in vivo Gene Regulatory Network

A wealth of physical interaction data between transcription factors (TFs) and DNA has been generated, but these interactions often do not have apparent regulatory consequences. Thus, equating physical interaction data with gene regulatory networks (GRNs) is problematic. Here, we comprehensively assay TF activity, rather than binding, to construct a network of gene regulatory interactions in the C. elegans intestine. By manually observing the in vivo tissue-specific knockdown of 921 TFs on a panel of 19 fluorescent transcriptional reporters, we identified a GRN of 411 interactions between 19 promoters and 177 TFs. This GRN shows only modest overlap with physical interactions, indicating that many regulatory interactions are indirect. We applied nested effects modeling to uncover information flow between TFs in the intestine that converges on a small set of physical TF-promoter interactions. We found numerous cell nonautonomous regulatory interactions, illustrating tissue-to-tissue communication. Altogether, our study illuminates the complexity of gene regulation in the context of a living animal.

Genome-Wide Mapping of Collier In Vivo Binding Sites Highlights Its Hierarchical Position in Different Transcription Regulatory Networks

Collier, the single Drosophila COE (Collier/EBF/Olf-1) transcription factor, is required in several developmental processes, including head patterning and specification of muscle and neuron identity during embryogenesis. To identify direct Collier (Col) targets in different cell types, we used ChIP-seq to map Col binding sites throughout the genome, at mid-embryogenesis. In vivo Col binding peaks were associated to 415 potential direct target genes. Gene Ontology analysis revealed a strong enrichment in proteins with DNA binding and/or transcription-regulatory properties. Characterization of a selection of candidates, using transgenic CRM-reporter assays, identified direct Col targets in dorso-lateral somatic muscles and specific neuron types in the central nervous system. These data brought new evidence that Col direct control of the expression of the transcription regulators apterous and eyes-absent (eya) is critical to specifying neuronal identities. They also showed that cross-regulation between col and eya in muscle progenitor cells is required for specification of muscle identity, revealing a new parallel between the myogenic regulatory networks operating in Drosophila and vertebrates. Col regulation of eya, both in specific muscle and neuronal lineages, may illustrate one mechanism behind the evolutionary diversification of Col biological roles.

Myocyte enhancer factor 2D regulates ectoderm specification and adhesion properties of animal cap cells in the early Xenopus embryo

In Xenopus, animal cap (AC) cells give rise to ectoderm and its derivatives: epidermis and the central nervous system. Ectoderm has long been considered a default pathway of embryonic development, with cells that are not under the influence of vegetal Nodal signaling adopting an ectodermal program of gene expression. In the present study, we describe the involvement of the animally-localized maternal transcription factor myocyte enhancer factor (Mef) 2D in regulating the identity of AC cells. We find that Mef2D is required for the formation of both ectodermal lineages: neural and epidermis. Gain and loss of function experiments indicate that Mef2D regulates early gastrula expression of key ectodermal/epidermal genes in the animal region. Mef2D controls the activity of zygotic bone morphogenetic protein (BMP) signaling known to dictate the epidermal differentiation program. Exogenous expression of Mef2D in vegetal blastomeres was sufficient to induce ectopic expression of ectoderm/epidermal genes in the vegetal half of the embryo, when Nodal signaling was inhibited. Depletion of Mef2D caused a loss of AC cell adhesion that was rescued by the expression of E-cadherin or bone morphogenetic protein 4. In addition, expression of Mef2D in the prospective endoderm caused unusual aggregation of vegetal cells with animal cells in vitro and inappropriate segregation to other germ layers in vivo. Mef2D cooperates with another animally-expressed transcription factor, FoxI1e. Together, they regulate the expression of genes encoding signaling proteins and the transcription factors that control the regional identity of animal cells. Therefore, we describe a new role for the animally-localized Mef2D protein in early ectoderm specification, which is similar to that of the vegetally-localized VegT in endoderm and mesoderm formation.

Identification of singles bar as a direct transcriptional target of Drosophila Myocyte enhancer factor-2 and a regulator of adult myoblast fusion

In Drosophila, myoblast fusion is a conserved process in which founder cells (FCs) and fusion competent myoblasts (FCMs) fuse to form a syncytial muscle fiber. Mutants for the myogenic regulator Myocyte enhancer factor-2 (MEF2) show a failure of myoblast fusion, indicating that MEF2 regulates the fusion process. Indeed, chromatin immunoprecipitation studies show that several genes involved in myoblast fusion are bound by MEF2 during embryogenesis. Of these, the MARVEL domain gene singles bar (sing), is down-regulated in MEF2 knockdown pupae, and has five consensus MEF2 binding sites within a 9000-bp region. To determine if MEF2 is an essential and direct regulator of sing during pupal muscle development, we identified a 315-bp myoblast enhancer of sing. This enhancer was active during myoblast fusion, and mutation of two MEF2 sites significantly decreased enhancer activity. We show that lack of sing expression resulted in adult lethality and muscle loss, due to a failure of fusion during the pupal stage. Additionally, we sought to determine if sing was required in either FCs or FCMs to support fusion. Interestingly, knockdown of sing in either population did not significantly affect fusion, however, knockdown in both FCs and FCMs resulted in muscles with significantly reduced nuclei numbers, provisionally indicating that sing function is required in either cell type, but not both. Finally, we found that MEF2 regulated sing expression at the embryonic stage through the same 315-bp enhancer, indicating that sing is a MEF2 target at both critical stages of myoblast fusion. Our studies define for the first time how MEF2 directly controls fusion at multiple stages of the life cycle, and provide further evidence that the mechanisms of fusion characterized in Drosophila embryos is also used in the formation of the more complex adult muscles.

MEF2D drives photoreceptor development through a genome-wide competition for tissue-specific enhancers

Organismal development requires the precise coordination of genetic programs to regulate cell fate and function. MEF2 transcription factors (TFs) play essential roles in this process but how these broadly expressed factors contribute to the generation of specific cell types during development is poorly understood. Here we show that despite being expressed in virtually all mammalian tissues, in the retina MEF2D binds to retina-specific enhancers and controls photoreceptor cell development. MEF2D achieves specificity by cooperating with a retina-specific factor CRX, which recruits MEF2D away from canonical MEF2 binding sites and redirects it to retina-specific enhancers that lack the consensus MEF2-binding sequence. Once bound to retina-specific enhancers, MEF2D and CRX co-activate the expression of photoreceptor-specific genes that are critical for retinal function. These findings demonstrate that broadly expressed TFs acquire specific functions through competitive recruitment to enhancers by tissue-specific TFs and through selective activation of these enhancers to regulate tissue-specific genes.

A cis-regulatory mutation in troponin-I of Drosophila reveals the importance of proper stoichiometry of structural proteins during muscle assembly

Rapid and high wing-beat frequencies achieved during insect flight are powered by the indirect flight muscles, the largest group of muscles present in the thorax. Any anomaly during the assembly and/or structural impairment of the indirect flight muscles gives rise to a flightless phenotype. Multiple mutagenesis screens in Drosophila melanogaster for defective flight behavior have led to the isolation and characterization of mutations that have been instrumental in the identification of many proteins and residues that are important for muscle assembly, function, and disease. In this article, we present a molecular-genetic characterization of a flightless mutation, flightless-H (fliH), originally designated as heldup-a (hdp-a). We show that fliH is a cis-regulatory mutation of the wings up A (wupA) gene, which codes for the troponin-I protein, one of the troponin complex proteins, involved in regulation of muscle contraction. The mutation leads to reduced levels of troponin-I transcript and protein. In addition to this, there is also coordinated reduction in transcript and protein levels of other structural protein isoforms that are part of the troponin complex. The altered transcript and protein stoichiometry ultimately culminates in unregulated acto-myosin interactions and a hypercontraction muscle phenotype. Our results shed new insights into the importance of maintaining the stoichiometry of structural proteins during muscle assembly for proper function with implications for the identification of mutations and disease phenotypes in other species, including humans.

The control operated by the cell cycle machinery on MEF2 stability contributes to the downregulation of CDKN1A and entry into S phase

MEF2s are pleiotropic transcription factors (TFs) which supervise multiple cellular activities. During the cell cycle, MEF2s are activated at the G0/G1 transition to orchestrate the expression of the immediate early genes in response to growth factor stimulation. Here we show that, in human and murine fibroblasts, MEF2 activities are downregulated during late G1. MEF2C and MEF2D interact with the E3 ligase F-box protein SKP2, which mediates their subsequent degradation through the ubiquitin-proteasome system. The cyclin-dependent kinase 4 (CDK4)/cyclin D1 complex phosphorylates MEF2D on serine residues 98 and 110, and phosphorylation of these residues is an important determinant for SKP2 binding. Unscheduled MEF2 transcription during the cell cycle reduces cell proliferation, whereas its containment sustains DNA replication. The CDK inhibitor p21/CDKN1A gene is a MEF2 target gene required to exert this antiproliferative influence. MEF2C and MEF2D bind a region within the first intron of CDKN1A, presenting epigenetic markers of open chromatin. Importantly, H3K27 acetylation within this regulative region depends on the presence of MEF2D. We propose that following the initial engagement in the G0/G1 transition, MEF2C and MEF2D must be polyubiquitylated and degraded during G1 progression to diminish the transcription of the CDKN1A gene, thus favoring entry into S phase.

The myogenic repressor gene Holes in muscles is a direct transcriptional target of Twist and Tinman in the Drosophila embryonic mesoderm

Understanding the regulatory circuitry controlling myogenesis is critical to understanding developmental mechanisms and developmentally-derived diseases. We analyzed the transcriptional regulation of a Drosophila myogenic repressor gene, Holes in muscles (Him). Previously, Him was shown to inhibit Myocyte enhancer factor-2 (MEF2) activity, and is expressed in myoblasts but not differentiating myotubes. We demonstrate that different phases of Him embryonic expression arises through the actions of different enhancers, and we characterize the enhancer required for its early mesoderm expression. This Him early mesoderm enhancer contains two conserved binding sites for the basic helix-loop-helix regulator Twist, and one binding site for the NK homeodomain protein Tinman. The sites for both proteins are required for enhancer activity in early embryos. Twist and Tinman activate the enhancer in tissue culture assays, and ectopic expression of either factor is sufficient to direct ectopic expression of a Him-lacZ reporter, or of the endogenous Him gene. Moreover, sustained expression of twist in the mesoderm up-regulates mesodermal Him expression in late embryos. Our findings provide a model to define mechanistically how Twist can both promotes myogenesis through direct activation of Mef2, and can place a brake on myogenesis, through direct activation of Him.

Gene coexpression networks reveal key drivers of phenotypic divergence in porcine muscle

Background

Domestication of the wild pig has led to obese and lean phenotype breeds, and evolutionary genome research has sought to identify the regulatory mechanisms underlying this phenotypic diversity. However, revealing the molecular mechanisms underlying muscle phenotype variation based on differentially expressed genes has proved to be difficult. To characterize the mechanisms regulating muscle phenotype variation under artificial selection, we aimed to provide an integrated view of genome organization by weighted gene coexpression network analysis.

Results

Our analysis was based on 20 publicly available next-generation sequencing datasets of lean and obese pig muscle generated from 10 developmental stages. The evolution of the constructed coexpression modules was examined using the genome resequencing data of 37 domestic pigs and 11 wild boars. Our results showed the regulation of muscle development might be more complex than had been previously acknowledged, and is regulated by the coordinated action of muscle, nerve and immunity related genes. Breed-specific modules that regulated muscle phenotype divergence were identified, and hundreds of hub genes with major roles in muscle development were determined to be responsible for key functional distinctions between breeds. Our evolutionary analysis showed that the role of changes in the coding sequence under positive selection in muscle phenotype divergence was minor.

Conclusions

Muscle phenotype divergence was found to be regulated by the divergence of coexpression network modules under artificial selection, and not by changes in the coding sequence of genes. Our results present multiple lines of evidence suggesting links between modules and muscle phenotypes, and provide insights into the molecular bases of genome organization in muscle development and phenotype variation.

Electronic supplementary material

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

MEF2 transcription factors regulate distinct gene programs in mammalian skeletal muscle differentiation

Skeletal muscle differentiation requires precisely coordinated transcriptional regulation of diverse gene programs that ultimately give rise to the specialized properties of this cell type. In Drosophila, this process is controlled, in part, by MEF2, the sole member of an evolutionarily conserved transcription factor family. By contrast, vertebrate MEF2 is encoded by four distinct genes, Mef2a, -b, -c, and -d, making it far more challenging to link this transcription factor to the regulation of specific muscle gene programs. Here, we have taken the first step in molecularly dissecting vertebrate MEF2 transcriptional function in skeletal muscle differentiation by depleting individual MEF2 proteins in myoblasts. Whereas MEF2A is absolutely required for proper myoblast differentiation, MEF2B, -C, and -D were found to be dispensable for this process. Furthermore, despite the extensive redundancy, we show that mammalian MEF2 proteins regulate a significant subset of nonoverlapping gene programs. These results suggest that individual MEF2 family members are able to recognize specific targets among the entire cohort of MEF2-regulated genes in the muscle genome. These findings provide opportunities to modulate the activity of MEF2 isoforms and their respective gene programs in skeletal muscle homeostasis and disease.

Genetic networks governing heart development

Animal genomes contain a code for construction of the body plan from a fertilized egg. Understanding how genome information is deciphered to create the complex multilayered regulatory systems that drive organismal development, and which become altered in disease, is one of the greatest challenges in the biological sciences. The development of methods that effectively represent and communicate the complexity inherent in gene regulatory networks remains a major barrier. This review introduces the philosophy of systems biology and discusses recent progress in understanding the development of the heart at a systems biology level.

Mef2 interacts with the Notch pathway during adult muscle development in Drosophila melanogaster

Myogenesis of indirect flight muscles (IFMs) in Drosophila melanogaster follows a well-defined cellular developmental scheme. During embryogenesis, a set of cells, the Adult Muscle Precursors (AMPs), are specified. These cells will become proliferating myoblasts during the larval stages which will then give rise to the adult IFMs. Although the cellular aspect of this developmental process is well studied, the molecular biology behind the different stages is still under investigation. In particular, the interactions required during the transition from proliferating myoblasts to differentiated myoblasts ready to fuse to the muscle fiber. It has been previously shown that the Notch pathway is active in proliferating myoblasts, and that this pathway is inhibited in developing muscle fibers. Furthermore, the Myocyte Enhancing Factor 2 (Mef2), Vestigial (Vg) and Scalloped (Sd) transcription factors are necessary for IFM development and that Vg is required for Notch pathway repression in differentiating fibers. Here we examine the interactions between Notch and Mef2 and mechanisms by which the Notch pathway is inhibited during differentiation. We show that Mef2 is capable of inhibiting the Notch pathway in non myogenic cells. A previous screen for Mef2 potential targets identified Delta a component of the Notch pathway. Dl is expressed in Mef2 and Sd-positive developing fibers. Our results show that Mef2 and possibly Sd regulate a Dl enhancer specifically expressed in the developing IFMs and that Mef2 is required for Dl expression in developing IFMs.

Global MEF2 target gene analysis in cardiac and skeletal muscle reveals novel regulation of DUSP6 by p38MAPK-MEF2 signaling

MEF2 plays a profound role in the regulation of transcription in cardiac and skeletal muscle lineages. To define the overlapping and unique MEF2A genomic targets, we utilized ChIP-exo analysis of cardiomyocytes and skeletal myoblasts. Of the 2783 and 1648 MEF2A binding peaks in skeletal myoblasts and cardiomyocytes, respectively, 294 common binding sites were identified. Genomic targets were compared to differentially expressed genes in RNA-seq analysis of MEF2A depleted myogenic cells, revealing two prominent genetic networks. Genes largely associated with muscle development were down-regulated by loss of MEF2A while up-regulated genes reveal a previously unrecognized function of MEF2A in suppressing growth/proliferative genes. Several up-regulated (Tprg, Mctp2, Kitl, Prrx1, Dusp6) and down-regulated (Atp1a2, Hspb7, Tmem182, Sorbs2, Lmod3) MEF2A target genes were chosen for further investigation. Interestingly, siRNA targeting of the MEF2A/D heterodimer revealed a somewhat divergent role in the regulation of Dusp6, a MAPK phosphatase, in cardiac and skeletal myogenic lineages. Furthermore, MEF2D functions as a p38MAPK-dependent repressor of Dusp6 in myoblasts. These data illustrate that MEF2 orchestrates both common and non-overlapping programs of signal-dependent gene expression in skeletal and cardiac muscle lineages.

Coordinated repression and activation of two transcriptional programs stabilizes cell fate during myogenesis

Molecular models of cell fate specification typically focus on the activation of specific lineage programs. However, the concurrent repression of unwanted transcriptional networks is also essential to stabilize certain cellular identities, as shown in a number of diverse systems and phyla. Here, we demonstrate that this dual requirement also holds true in the context of Drosophila myogenesis. By integrating genetics and genomics, we identified a new role for the pleiotropic transcriptional repressor Tramtrack69 in myoblast specification. Drosophila muscles are formed through the fusion of two discrete cell types: founder cells (FCs) and fusion-competent myoblasts (FCMs). When tramtrack69 is removed, FCMs appear to adopt an alternative muscle FC-like fate. Conversely, ectopic expression of this repressor phenocopies muscle defects seen in loss-of-function lame duck mutants, a transcription factor specific to FCMs. This occurs through Tramtrack69-mediated repression in FCMs, whereas Lame duck activates a largely distinct transcriptional program in the same cells. Lineage-specific factors are therefore not sufficient to maintain FCM identity. Instead, their identity appears more plastic, requiring the combination of instructive repressive and activating programs to stabilize cell fate.

The transient expression of miR-203 and its inhibiting effects on skeletal muscle cell proliferation and differentiation

Previous studies have shown that miR-203 is a skin-specific microRNA (miRNA) with a profound role in skin cell differentiation. However, emerging microarray and deep sequencing data revealed that miR-203 is also expressed in embryonic skeletal muscle and myoblasts. In this study, we found that miR-203 was transiently upregulated in chicken embryos on days 10 to 16 (E10-E16) and was sharply downregulated and even not expressed after E16 in chicken embryonic skeletal muscle. Histological profiles and weight variations of embryo skeletal muscle revealed that miR-203 expression is correlated with muscle development. In vitro experiments showed that miR-203 exhibited downregulated expression during myoblast differentiation into myotubes. miR-203 overexpression inhibited myoblast proliferation and differentiation, whereas its loss-of-function increased myoblast proliferation and differentiation. During myogenesis, miR-203 can target and inhibit the expression of c-JUN and MEF2C, which were important for cell proliferation and muscle development, respectively. The overexpression of c-JUN significantly promoted myoblast proliferation. Conversely, knockdown of c-JUN by siRNA suppressed myoblast proliferation. In addition, the knockdown of MEF2C by siRNA significantly inhibited myoblast differentiation. Altogether, these data not only suggested that the expression of miR-203 is transitory during chicken skeletal muscle development but also showed a novel role of miR-203 in inhibiting skeletal muscle cell proliferation and differentiation by repressing c-JUN and MEF2C, respectively.

Emergence of tissue sensitivity to Hox protein levels underlies the evolution of an adaptive morphological trait

Growth control scales morphological attributes and, therefore, provides a critical contribution to the evolution of adaptive traits. Yet, the genetic mechanisms underlying growth in the context of specific ecological adaptations are poorly understood. In water striders, adaptation to locomotion on the water surface is associated with allometric and functional changes in thoracic appendages, such that T2-legs, used as propelling oars, are longer than T3-legs, used as steering rudders. The Hox gene Ubx establishes this derived morphology by elongating T2-legs but shortening T3-legs. Using gene expression assays, RNAi knockdown, and comparative transcriptomics, we demonstrate that the evolution of water surface rowing as a novel means of locomotion is associated with the evolution of a dose-dependent promoting-repressing effect of Ubx on leg growth. In the water strider Limnoporus dissortis, T3-legs express six to seven times higher levels of Ubx compared to T2-legs. Ubx RNAi shortens T2-legs and the severity of this phenotype increases with increased depletion of Ubx protein. Conversely, Ubx RNAi lengthens T3-legs but this phenotype is partially rescued when Ubx protein is further depleted. This dose-dependent effect of Ubx on leg growth is absent in non-rowing relatives that retain the ancestral relative leg length. We also show that the spatial patterns of expression of dpp, wg, hh, egfr, dll, exd, hth, and dac are unchanged in Ubx RNAi treatments. This indicates that the dose-dependent opposite effect of Ubx on T2- and T3-legs operates without any apparent effect on the spatial expression of major leg patterning genes. Our data suggest that scaling of adaptive allometries can evolve through changes in the levels of expression of Hox proteins early during ontogeny, and in the sensitivity of the tissues that express them, without any major effects on pattern formation.

•

Ubx is generally expressed at higher levels in T3- relative to T2-legs in semi-aquatic insects.

•

It is only in the derived Gerridae where the high levels of Ubx result in reduced T3-leg length.

•

In the Gerridae, the response of leg tissues to Ubx levels is bimodal.

•

Changes in Ubx regulation and function have evolved in Limnoporus without disrupting patterning hierarchies.

•

Changes in Hox protein levels and emergence of tissue sensitivity to these levels can shape adaptive morphological traits.

The genome-wide molecular signature of transcription factors in leukemia

Transcription factors control expression of genes essential for the normal functioning of the hematopoietic system and regulate development of distinct blood cell types. During leukemogenesis, aberrant regulation of transcription factors such as RUNX1, CBFβ, MLL, C/EBPα, SPI1, GATA, and TAL1 is central to the disease. Here, we will discuss the mechanisms of transcription factor deregulation in leukemia and how in recent years next-generation sequencing approaches have helped to elucidate the molecular role of many of these aberrantly expressed transcription factors. We will focus on the complexes in which these factors reside, the role of posttranslational modification of these factors, their involvement in setting up higher order chromatin structures, and their influence on the local epigenetic environment. We suggest that only comprehensive knowledge on all these aspects will increase our understanding of aberrant gene expression in leukemia as well as open new entry points for therapeutic intervention.

Two enhancers control transcription of Drosophila muscleblind in the embryonic somatic musculature and in the central nervous system

The phylogenetically conserved family of Muscleblind proteins are RNA-binding factors involved in a variety of gene expression processes including alternative splicing regulation, RNA stability and subcellular localization, and miRNA biogenesis, which typically contribute to cell-type specific differentiation. In humans, sequestration of Muscleblind-like proteins MBNL1 and MBNL2 has been implicated in degenerative disorders, particularly expansion diseases such as myotonic dystrophy type 1 and 2. Drosophila muscleblind was previously shown to be expressed in embryonic somatic and visceral muscle subtypes, and in the central nervous system, and to depend on Mef2 for transcriptional activation. Genomic approaches have pointed out candidate gene promoters and tissue-specific enhancers, but experimental confirmation of their regulatory roles was lacking. In our study, luciferase reporter assays in S2 cells confirmed that regions P1 (515 bp) and P2 (573 bp), involving the beginning of exon 1 and exon 2, respectively, were able to initiate RNA transcription. Similarly, transgenic Drosophila embryos carrying enhancer reporter constructs supported the existence of two regulatory regions which control embryonic expression of muscleblind in the central nerve cord (NE, neural enhancer; 830 bp) and somatic (skeletal) musculature (ME, muscle enhancer; 3.3 kb). Both NE and ME were able to boost expression from the Hsp70 heterologous promoter. In S2 cell assays most of the ME enhancer activation could be further narrowed down to a 1200 bp subregion (ME.3), which contains predicted binding sites for the Mef2 transcription factor. The present study constitutes the first characterization of muscleblind enhancers and will contribute to a deeper understanding of the transcriptional regulation of the gene.

Practical guidelines for the comprehensive analysis of ChIP-seq data

Mapping the chromosomal locations of transcription factors, nucleosomes, histone modifications, chromatin remodeling enzymes, chaperones, and polymerases is one of the key tasks of modern biology, as evidenced by the Encyclopedia of DNA Elements (ENCODE) Project. To this end, chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) is the standard methodology. Mapping such protein-DNA interactions in vivo using ChIP-seq presents multiple challenges not only in sample preparation and sequencing but also for computational analysis. Here, we present step-by-step guidelines for the computational analysis of ChIP-seq data. We address all the major steps in the analysis of ChIP-seq data: sequencing depth selection, quality checking, mapping, data normalization, assessment of reproducibility, peak calling, differential binding analysis, controlling the false discovery rate, peak annotation, visualization, and motif analysis. At each step in our guidelines we discuss some of the software tools most frequently used. We also highlight the challenges and problems associated with each step in ChIP-seq data analysis. We present a concise workflow for the analysis of ChIP-seq data in Figure 1 that complements and expands on the recommendations of the ENCODE and modENCODE projects. Each step in the workflow is described in detail in the following sections.

Transcriptional networks regulating the costamere, sarcomere, and other cytoskeletal structures in striated muscle

Structural abnormalities in striated muscle have been observed in numerous transcription factor gain- and loss-of-function phenotypes in animal and cell culture model systems, indicating that transcription is important in regulating the cytoarchitecture. While most characterized cytoarchitectural defects are largely indistinguishable by histological and ultrastructural criteria, analysis of dysregulated gene expression in each mutant phenotype has yielded valuable information regarding specific structural gene programs that may be uniquely controlled by each of these transcription factors. Linking the formation and maintenance of each subcellular structure or subset of proteins within a cytoskeletal compartment to an overlapping but distinct transcription factor cohort may enable striated muscle to control cytoarchitectural function in an efficient and specific manner. Here we summarize the available evidence that connects transcription factors, those with established roles in striated muscle such as MEF2 and SRF, as well as other non-muscle transcription factors, to the regulation of a defined cytoskeletal structure. The notion that genes encoding proteins localized to the same subcellular compartment are coordinately transcriptionally regulated may prompt rationally designed approaches that target specific transcription factor pathways to correct structural defects in muscle disease.

A Cdx4-Sall4 regulatory module controls the transition from mesoderm formation to embryonic hematopoiesis

Deletion of caudal/cdx genes alters hox gene expression and causes defects in posterior tissues and hematopoiesis. Yet, the defects in hox gene expression only partially explain these phenotypes. To gain deeper insight into Cdx4 function, we performed chromatin immunoprecipitation sequencing (ChIP-seq) combined with gene-expression profiling in zebrafish, and identified the transcription factor spalt-like 4 (sall4) as a Cdx4 target. ChIP-seq revealed that Sall4 bound to its own gene locus and the cdx4 locus. Expression profiling showed that Cdx4 and Sall4 coregulate genes that initiate hematopoiesis, such as hox, scl, and lmo2. Combined cdx4/sall4 gene knockdown impaired erythropoiesis, and overexpression of the Cdx4 and Sall4 target genes scl and lmo2 together rescued the erythroid program. These findings suggest that auto- and cross-regulation of Cdx4 and Sall4 establish a stable molecular circuit in the mesoderm that facilitates the activation of the blood-specific program as development proceeds.

•

Cdx4 and Sall4 bind to each other’s genomic loci

•

Cdx4 and Sall4 coregulate genes responsible for the mesoderm-to-blood transition

•

Scl and Lmo2 overexpression rescues blood defects in cdx4/sall4 double morphants

Modulation of lysine methylation in myocyte enhancer factor 2 during skeletal muscle cell differentiation

Myocyte enhancer factor 2 (MEF2) is a family of transcription factors that regulates many processes, including muscle differentiation. Due to its many target genes, MEF2D requires tight regulation of transcription activity over time and by location. Epigenetic modifiers have been suggested to regulate MEF2-dependent transcription via modifications to histones and MEF2. However, the modulation of MEF2 activity by lysine methylation, an important posttranslational modification that alters the activities of transcription factors, has not been studied. We report the reversible lysine methylation of MEF2D by G9a and LSD1 as a regulatory mechanism of MEF2D activity and skeletal muscle differentiation. G9a methylates lysine-267 of MEF2D and represses its transcriptional activity, but LSD1 counteracts it. This residue is highly conserved between MEF2 members in mammals. During myogenic differentiation of C2C12 mouse skeletal muscle cells, the methylation of MEF2D by G9a decreased, on which MEF2D-dependent myogenic genes were upregulated. We have also identified lysine-267 as a methylation/demethylation site and demonstrate that the lysine methylation state of MEF2D regulates its transcriptional activity and skeletal muscle cell differentiation.

MEF2 is an in vivo immune-metabolic switch

Infections disturb metabolic homeostasis in many contexts, but the underlying connections are not completely understood. To address this, we use paired genetic and computational screens in Drosophila to identify transcriptional regulators of immunity and pathology and their associated target genes and physiologies. We show that Mef2 is required in the fat body for anabolic function and the immune response. Using genetic and biochemical approaches, we find that MEF2 is phosphorylated at a conserved site in healthy flies and promotes expression of lipogenic and glycogenic enzymes. Upon infection, this phosphorylation is lost, and the activity of MEF2 changes—MEF2 now associates with the TATA binding protein to bind a distinct TATA box sequence and promote antimicrobial peptide expression. The loss of phosphorylated MEF2 contributes to loss of anabolic enzyme expression in Gram-negative bacterial infection. MEF2 is thus a critical transcriptional switch in the adult fat body between metabolism and immunity.

•

Mef2 is required in Drosophila for immune function and storage of fat and glycogen

•

MEF2 is phosphorylated in vivo at a conserved site (T20) to promote anabolism

•

Infection reduces phospho-T20, allowing MEF2 to bind TBP and an immune TATA box

•

MEF2 dephosphorylation leads to metabolic dysfunction in Gram-negative infection

Polyethyleneimine mediated DNA transfection in schistosome parasites and regulation of the WNT signaling pathway by a dominant-negative SmMef2

Schistosomiasis is a serious global problem and the second most devastating parasitic disease following malaria. Parasitic worms of the genus Schistosoma are the causative agents of schistosomiasis and infect more than 240 million people worldwide. The paucity of molecular tools to manipulate schistosome gene expression has made an understanding of genetic pathways in these parasites difficult, increasing the challenge of identifying new potential drugs for treatment. Here, we describe the use of a formulation of polyethyleneimine (PEI) as an alternative to electroporation for the efficacious transfection of genetic material into schistosome parasites. We show efficient expression of genes from a heterologous CMV promoter and from the schistosome Sm23 promoter. Using the schistosome myocyte enhancer factor 2 (SmMef2), a transcriptional activator critical for myogenesis and other developmental pathways, we describe the development of a dominant-negative form of the schistosome Mef2. Using this mutant, we provide evidence that SmMef2 may regulate genes in the WNT pathway. We also show that SmMef2 regulates its own expression levels. These data demonstrate the use of PEI to facilitate effective transfection of nucleic acids into schistosomes, aiding in the study of schistosome gene expression and regulation, and development of genetic tools for the characterization of molecular pathways in these parasites.

The transcription factor Mef2 links the Drosophila core clock to Fas2, neuronal morphology, and circadian behavior

The transcription factor Mef2 regulates activity-dependent neuronal plasticity and morphology in mammals, and clock neurons are reported to experience activity-dependent circadian remodeling in Drosophila. We show here that Mef2 is required for this daily fasciculation-defasciculation cycle. Moreover, the master circadian transcription complex CLK/CYC directly regulates Mef2 transcription. ChIP-Chip analysis identified numerous Mef2 target genes implicated in neuronal plasticity, including the cell-adhesion gene Fas2. Genetic epistasis experiments support this transcriptional regulatory hierarchy, CLK/CYC- > Mef2- > Fas2, indicate that it influences the circadian fasciculation cycle within pacemaker neurons, and suggest that this cycle also contributes to circadian behavior. Mef2 therefore transmits clock information to machinery involved in neuronal remodeling, which contributes to locomotor activity rhythms.

Transcriptional regulation of mesoderm genes by MEF2D during early Xenopus development

In Xenopus, specification of the three germ layers is one of the earliest developmental decisions occurring prior to gastrulation. The maternally-expressed vegetally-localized transcription factor VegT has a central role in cell autonomous specification of endoderm and in the generation of mesoderm-inducing signals. Yet, marginally-expressed transcription factors that cooperate with mesoderm-inducing signals are less investigated. Here we report that the transcription factors MEF2A and MEF2D are expressed in the animal hemisphere before mid-blastula transition. At the initiation of zygotic transcription, expression of MEF2D expands into the marginal region that gives rise to mesoderm. Knockdown of MEF2D delayed gastrulation movements, prevented embryo elongation at the subsequent tailbud stage and caused severe defects in axial tissues. At the molecular level, MEF2D knockdown reduced the expression of genes involved in mesoderm formation and patterning. We also report that MEF2D functions with FGF signaling in a positive feedback loop; each augments the expression of the other in the marginal region and both are necessary for mesodermal gene expression. One target of MEF2D is the Nodal-related 1 gene (Xnr1) that mediates some of MEF2D mesodermal activities. Chromatin immunoprecipitation analysis revealed that MEF2D associates with transcriptional regulatory sequences of the Xnr1 gene. Several MEF2 binding sites within the proximal promoter region of Xnr1 were identified by their in vitro association with MEF2D protein. The same promoter region was necessary but not sufficient to mediate MEF2D activity in a reporter gene assay. In sum, our results indicate that the MEF2D protein is a key transcription factor in the marginal zone acting in a positive feedback loop with FGF signaling that promotes mesoderm specification at late blastula stages.

Conserved non-coding elements and cis regulation: actions speak louder than words

It is a truth (almost) universally acknowledged that conserved non-coding genomic sequences function in the cis regulation of neighbouring genes. But is this a misconception? The literature is strewn with examples of conserved non-coding sequences being able to drive reporter expression, but the extent to which such sequences are actually used endogenously in vivo is only now being rigorously explored using unbiased genome-scale approaches. Here, we review the emerging picture, examining the extent to which conserved non-coding sequences equivalently regulate gene expression in different species, or at different developmental stages, and how genomics approaches are revealing the relationship between sequence conservation and functional use of cis-regulatory elements.

Mef2A, a homologue of animal Mef2 transcription factors, regulates cell differentiation in Dictyostelium discoideum

Background

Transcription factors from the MADS-box family play a relevant role in cell differentiation and development and include the animal SRF (serum response factor) and MEF2 (myocyte enhancer factor 2) proteins. The social amoeba Dictyostelium discoideum contains four genes coding for MADS-box transcription factors, two of these genes code for proteins that are more similar to SRF, and the other two code for proteins that are more similar to MEF2 animal factors.

Results

The biological function of one of the two genes that codes for MEF2-related proteins, a gene known as mef2A, is described in this article. This gene is expressed under the transcriptional control of two alternative promoters in growing cells, and its expression is induced during development in prespore cells. Mutant strains where the mef2A gene has been partially deleted were generated to study its biological function. The mutant strains showed reduced growth when feeding on bacteria and were able to develop and form fruiting bodies, but spore production was significantly reduced. A study of developmental markers showed that prespore cells differentiation was impaired in the mutant strains. When mutant and wild-type cells were set to develop in chimeras, mutant spores were underrepresented in the fruiting bodies. The mutant cells were also unable to form spores in vitro. In addition, mutant cells also showed a poor contribution to the formation of the tip-organizer and the upper region of slugs and culminant structures. In agreement with these observations, a comparison of the genes transcribed by mutant and wild-type strains during development indicated that prestalk gene expression was enhanced, while prespore gene expression decreased in the mef2A^- strain.

Conclusions

Our data shows that mef2A plays a role in cell differentiation in D. discoideum and modulates the expression of prespore and prestalk genes.

Temporal mapping of CEBPA and CEBPB binding during liver regeneration reveals dynamic occupancy and specific regulatory codes for homeostatic and cell cycle gene batteries

Dynamic shifts in transcription factor binding are central to the regulation of biological processes by allowing rapid changes in gene transcription. However, very few genome-wide studies have examined how transcription factor occupancy is coordinated temporally in vivo in higher animals. Here, we quantified the genome-wide binding patterns of two key hepatocyte transcription factors, CEBPA and CEBPB (also known as C/EBPalpha and C/EBPbeta), at multiple time points during the highly dynamic process of liver regeneration elicited by partial hepatectomy in mouse. Combining these profiles with RNA polymerase II binding data, we find three temporal classes of transcription factor binding to be associated with distinct sets of regulated genes involved in the acute phase response, metabolic/homeostatic functions, or cell cycle progression. Moreover, we demonstrate a previously unrecognized early phase of homeostatic gene expression prior to S-phase entry. By analyzing the three classes of CEBP bound regions, we uncovered mutually exclusive sets of sequence motifs, suggesting temporal codes of CEBP recruitment by differential cobinding with other factors. These findings were validated by sequential ChIP experiments involving a panel of central transcription factors and/or by comparison to external ChIP-seq data. Our quantitative investigation not only provides in vivo evidence for the involvement of many new factors in liver regeneration but also points to similarities in the circuitries regulating self-renewal of differentiated cells. Taken together, our work emphasizes the power of global temporal analyses of transcription factor occupancy to elucidate mechanisms regulating dynamic biological processes in complex higher organisms.

Hedgehog signaling regulates MyoD expression and activity

The inhibition of MyoD expression is important for obtaining muscle progenitors that can replenish the satellite cell niche during muscle repair. Progenitors could be derived from either embryonic stem cells or satellite cells. Hedgehog (Hh) signaling is important for MyoD expression during embryogenesis and adult muscle regeneration. To date, the mechanistic understanding of MyoD regulation by Hh signaling is unclear. Here, we demonstrate that the Hh effector, Gli2, regulates MyoD expression and associates with MyoD gene elements. Gain- and loss-of-function experiments in pluripotent P19 cells show that Gli2 activity is sufficient and required for efficient MyoD expression during skeletal myogenesis. Inhibition of Hh signaling reduces MyoD expression during satellite cell activation in vitro. In addition to regulating MyoD expression, Hh signaling regulates MyoD transcriptional activity, and MyoD activates Hh signaling in myogenic conversion assays. Finally, Gli2, MyoD, and MEF2C form a protein complex, which enhances MyoD activity on skeletal muscle-related promoters. We therefore link Hh signaling to the function and expression of MyoD protein during myogenesis in stem cells.

Differential RNAi screening provides insights into the rewiring of signalling networks during oxidative stress

Reactive Oxygen Species (ROS) are a natural by-product of cellular growth and proliferation, and are required for fundamental processes such as protein-folding and signal transduction. However, ROS accumulation, and the onset of oxidative stress, can negatively impact cellular and genomic integrity. Signalling networks have evolved to respond to oxidative stress by engaging diverse enzymatic and non-enzymatic antioxidant mechanisms to restore redox homeostasis. The architecture of oxidative stress response networks during periods of normal growth, and how increased ROS levels dynamically reconfigure these networks are largely unknown. In order to gain insight into the structure of signalling networks that promote redox homeostasis we first performed genome-scale RNAi screens to identify novel suppressors of superoxide accumulation. We then infer relationships between redox regulators by hierarchical clustering of phenotypic signatures describing how gene inhibition affects superoxide levels, cellular viability, and morphology across different genetic backgrounds. Genes that cluster together are likely to act in the same signalling pathway/complex and thus make "functional interactions". Moreover we also calculate differential phenotypic signatures describing the difference in cellular phenotypes following RNAi between untreated cells and cells submitted to oxidative stress. Using both phenotypic signatures and differential signatures we construct a network model of functional interactions that occur between components of the redox homeostasis network, and how such interactions become rewired in the presence of oxidative stress. This network model predicts a functional interaction between the transcription factor Jun and the IRE1 kinase, which we validate in an orthogonal assay. We thus demonstrate the ability of systems-biology approaches to identify novel signalling events.

Mef2d acts upstream of muscle identity genes and couples lateral myogenesis to dermomyotome formation in Xenopus laevis

Xenopus myotome is formed by a first medial and lateral myogenesis directly arising from the presomitic mesoderm followed by a second myogenic wave emanating from the dermomyotome. Here, by a series of gain and loss of function experiments, we showed that Mef2d, a member of the Mef2 family of MADS-box transcription factors, appeared as an upstream regulator of lateral myogenesis, and as an inducer of dermomyotome formation at the beginning of neurulation. In the lateral presomitic cells, we showed that Mef2d transactivates Myod expression which is necessary for lateral myogenesis. In the most lateral cells of the presomitic mesoderm, we showed that Mef2d and Paraxis (Tcf15), a member of the Twist family of transcription factors, were co-localized and activate directly the expression of Meox2, which acts upstream of Pax3 expression during dermomyotome formation. Cell tracing experiments confirm that the most lateral Meox2 expressing cells of the presomitic mesoderm correspond to the dermomyotome progenitors since they give rise to the most dorsal cells of the somitic mesoderm. Thus, Xenopus Mef2d couples lateral myogenesis to dermomyotome formation before somite segmentation. These results together with our previous works reveal striking similarities between dermomyotome and tendon formation in Xenopus: both develop in association with myogenic cells and both involve a gene transactivation pathway where one member of the Mef2 family, Mef2d or Mef2c, cooperates with a bHLH protein of the Twist family, Paraxis or Scx (Scleraxis) respectively. We propose that these shared characteristics in Xenopus laevis reflect the existence of a vertebrate ancestral mechanism which has coupled the development of the myogenic cells to the formation of associated tissues during somite compartmentalization.

Predicting spatial and temporal gene expression using an integrative model of transcription factor occupancy and chromatin state

Precise patterns of spatial and temporal gene expression are central to metazoan complexity and act as a driving force for embryonic development. While there has been substantial progress in dissecting and predicting cis-regulatory activity, our understanding of how information from multiple enhancer elements converge to regulate a gene's expression remains elusive. This is in large part due to the number of different biological processes involved in mediating regulation as well as limited availability of experimental measurements for many of them. Here, we used a Bayesian approach to model diverse experimental regulatory data, leading to accurate predictions of both spatial and temporal aspects of gene expression. We integrated whole-embryo information on transcription factor recruitment to multiple cis-regulatory modules, insulator binding and histone modification status in the vicinity of individual gene loci, at a genome-wide scale during Drosophila development. The model uses Bayesian networks to represent the relation between transcription factor occupancy and enhancer activity in specific tissues and stages. All parameters are optimized in an Expectation Maximization procedure providing a model capable of predicting tissue- and stage-specific activity of new, previously unassayed genes. Performing the optimization with subsets of input data demonstrated that neither enhancer occupancy nor chromatin state alone can explain all gene expression patterns, but taken together allow for accurate predictions of spatio-temporal activity. Model predictions were validated using the expression patterns of more than 600 genes recently made available by the BDGP consortium, demonstrating an average 15-fold enrichment of genes expressed in the predicted tissue over a naïve model. We further validated the model by experimentally testing the expression of 20 predicted target genes of unknown expression, resulting in an accuracy of 95% for temporal predictions and 50% for spatial. While this is, to our knowledge, the first genome-wide approach to predict tissue-specific gene expression in metazoan development, our results suggest that integrative models of this type will become more prevalent in the future.

Development is a complex process in which a single cell gives rise to a multi-cellular organism comprised of diverse cell types and well-organized tissues. This transformation requires tightly coordinated expression, both spatially and temporally, of hundreds to thousands of genes specific to any given tissue. To orchestrate these patterns, gene expression is regulated at multiple steps, from TF binding to cis-regulatory modules, general transcription factor and RNA polymerase II recruitment to promoters, chromatin remodeling, and three-dimensional looping interactions. Despite this level of complexity, the regulation of gene expression is typically modeled in the context of transcription factor binding and a single enhancer's activity as this is where the majority of experimental data is available. Recent advances in the measurement of chromatin modifications and insulator binding during embryogenesis provide new datasets that can be used for modeling gene expression. Here we use a Bayesian approach to integrate all three levels of information to combine the activity of multiple regulatory elements into a single model of a gene's expression, implementing an expectation maximization strategy to overcome the problem of missing data. Importantly, while the data for histone modifications and insulator binding represents merged signals from all cells in the embryo, the model can extract cell type specific and stage-specific predictions on gene expression for hundreds of genes of unknown expression.

Integrative analysis of the zinc finger transcription factor Lame duck in the Drosophila myogenic gene regulatory network

Contemporary high-throughput technologies permit the rapid identification of transcription factor (TF) target genes on a genome-wide scale, yet the functional significance of TFs requires knowledge of target gene expression patterns, cooperating TFs, and cis-regulatory element (CRE) structures. Here we investigated the myogenic regulatory network downstream of the Drosophila zinc finger TF Lame duck (Lmd) by combining both previously published and newly performed genomic data sets, including ChIP sequencing (ChIP-seq), genome-wide mRNA profiling, cell-specific expression patterns of putative transcriptional targets, analysis of histone mark signatures, studies of TF cooccupancy by additional mesodermal regulators, TF binding site determination using protein binding microarrays (PBMs), and machine learning of candidate CRE motif compositions. Our findings suggest that Lmd orchestrates an extensive myogenic regulatory network, a conclusion supported by the identification of Lmd-dependent genes, histone signatures of Lmd-bound genomic regions, and the relationship of these features to cell-specific gene expression patterns. The heterogeneous cooccupancy of Lmd-bound regions with additional mesodermal regulators revealed that different transcriptional inputs are used to mediate similar myogenic gene expression patterns. Machine learning further demonstrated diverse combinatorial motif patterns within tissue-specific Lmd-bound regions. PBM analysis established the complete spectrum of Lmd DNA binding specificities, and site-directed mutagenesis of Lmd and additional newly discovered motifs in known enhancers demonstrated the critical role of these TF binding sites in supporting full enhancer activity. Collectively, these findings provide insights into the transcriptional codes regulating muscle gene expression and offer a generalizable approach for similar studies in other systems.

miR-92b regulates Mef2 levels through a negative-feedback circuit during Drosophila muscle development

Mef2 is the key transcription factor for muscle development and differentiation in Drosophila. It activates hundreds of downstream target genes, including itself. Precise control of Mef2 levels is essential for muscle development as different Mef2 protein levels activate distinct sets of muscle genes, but how this is achieved remains unclear. Here, we have identified a novel heart- and muscle-specific microRNA, miR-92b, which is activated by Mef2 and subsequently downregulates Mef2 through binding to its 3'UTR, forming a negative regulatory circuit that fine-tunes the level of Mef2. Deletion of miR-92b caused abnormally high Mef2 expression, leading to muscle defects and lethality. Blocking miR-92b function using microRNA sponge techniques also increased Mef2 levels and caused muscle defects similar to those seen with the miR-92b deletion. Additionally, overexpression of miR-92b reduced Mef2 levels and caused muscle defects similar to those seen in Mef2 RNAi, and Mef2 overexpression led to reversal of these defects. Our results suggest that the negative feedback circuit between miR-92b and Mef2 efficiently maintains the stable expression of both components that is required for homeostasis during Drosophila muscle development.

Transcription factors: from enhancer binding to developmental control

Developmental progression is driven by specific spatiotemporal domains of gene expression, which give rise to stereotypically patterned embryos even in the presence of environmental and genetic variation. Views of how transcription factors regulate gene expression are changing owing to recent genome-wide studies of transcription factor binding and RNA expression. Such studies reveal patterns that, at first glance, seem to contrast with the robustness of the developmental processes they encode. Here, we review our current knowledge of transcription factor function from genomic and genetic studies and discuss how different strategies, including extensive cooperative regulation (both direct and indirect), progressive priming of regulatory elements, and the integration of activities from multiple enhancers, confer specificity and robustness to transcriptional regulation during development.

Studying and modelling dynamic biological processes using time-series gene expression data

Biological processes are often dynamic, thus researchers must monitor their activity at multiple time points. The most abundant source of information regarding such dynamic activity is time-series gene expression data. These data are used to identify the complete set of activated genes in a biological process, to infer their rates of change, their order and their causal effects and to model dynamic systems in the cell. In this Review we discuss the basic patterns that have been observed in time-series experiments, how these patterns are combined to form expression programs, and the computational analysis, visualization and integration of these data to infer models of dynamic biological systems.

i-cisTarget: an integrative genomics method for the prediction of regulatory features and cis-regulatory modules

The field of regulatory genomics today is characterized by the generation of high-throughput data sets that capture genome-wide transcription factor (TF) binding, histone modifications, or DNAseI hypersensitive regions across many cell types and conditions. In this context, a critical question is how to make optimal use of these publicly available datasets when studying transcriptional regulation. Here, we address this question in Drosophila melanogaster for which a large number of high-throughput regulatory datasets are available. We developed i-cisTarget (where the 'i' stands for integrative), for the first time enabling the discovery of different types of enriched 'regulatory features' in a set of co-regulated sequences in one analysis, being either TF motifs or 'in vivo' chromatin features, or combinations thereof. We have validated our approach on 15 co-expressed gene sets, 21 ChIP data sets, 628 curated gene sets and multiple individual case studies, and show that meaningful regulatory features can be confidently discovered; that bona fide enhancers can be identified, both by in vivo events and by TF motifs; and that combinations of in vivo events and TF motifs further increase the performance of enhancer prediction.

Notch and Mef2 synergize to promote proliferation and metastasis through JNK signal activation in Drosophila

Genetic analyses in Drosophila revealed a synergy between Notch and the pleiotropic transcription factor Mef2 (myocyte enhancer factor 2), which profoundly influences proliferation and metastasis. We show that these hyperproliferative and invasive Drosophila phenotypes are attributed to upregulation of eiger, a member of the tumour necrosis factor superfamily of ligands, and the consequent activation of Jun N-terminal kinase signalling, which in turn triggers the expression of the invasive marker MMP1. Expression studies in human breast tumour samples demonstrate correlation between Notch and Mef2 paralogues and support the notion that Notch-MEF2 synergy may be significant for modulating human mammary oncogenesis.

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

The regulation of gene expression is mediated at the transcriptional level by enhancer regions that are bound by sequence-specific transcription factors (TFs). Recent studies have shown that the in vivo binding sites of single TFs differ between developmental or cellular contexts. How this context-specific binding is encoded in the cis-regulatory DNA sequence has, however, remained unclear. We computationally dissect context-specific TF binding sites in Drosophila, Caenorhabditis elegans, mouse, and human and find distinct combinations of sequence motifs for partner factors, which are predictive and reveal specific motif requirements of individual binding sites. We predict that TF binding in the early Drosophila embryo depends on motifs for the early zygotic TFs Vielfaltig (also known as Zelda) and Tramtrack. We validate experimentally that the activity of Twist-bound enhancers and Twist binding itself depend on Vielfaltig motifs, suggesting that Vielfaltig is more generally important for early transcription. Our finding that the motif content can predict context-specific binding and that the predictions work across different Drosophila species suggests that characteristic motif combinations are shared between sites, revealing context-specific motif codes (cis-regulatory signatures), which appear to be conserved during evolution. Taken together, this study establishes a novel approach to derive predictive cis-regulatory motif requirements for individual TF binding sites and enhancers. Importantly, the method is generally applicable across different cell types and organisms to elucidate cis-regulatory sequence determinants and the corresponding trans-acting factors from the increasing number of tissue- and cell-type-specific TF binding studies.