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

Gene networks in Drosophila melanogaster: integrating experimental data to predict gene function

BACKGROUND

Discovering the functions of all genes is a central goal of contemporary biomedical research. Despite considerable effort, we are still far from achieving this goal in any metazoan organism. Collectively, the growing body of high-throughput functional genomics data provides evidence of gene function, but remains difficult to interpret.

RESULTS

We constructed the first network of functional relationships for Drosophila melanogaster by integrating most of the available, comprehensive sets of genetic interaction, protein-protein interaction, and microarray expression data. The complete integrated network covers 85% of the currently known genes, which we refined to a high confidence network that includes 20,000 functional relationships among 5,021 genes. An analysis of the network revealed a remarkable concordance with prior knowledge. Using the network, we were able to infer a set of high-confidence Gene Ontology biological process annotations on 483 of the roughly 5,000 previously unannotated genes. We also show that this approach is a means of inferring annotations on a class of genes that cannot be annotated based solely on sequence similarity. Lastly, we demonstrate the utility of the network through reanalyzing gene expression data to both discover clusters of coregulated genes and compile a list of candidate genes related to specific biological processes.

CONCLUSIONS

Here we present the the first genome-wide functional gene network in D. melanogaster. The network enables the exploration, mining, and reanalysis of experimental data, as well as the interpretation of new data. The inferred annotations provide testable hypotheses of previously uncharacterized genes.

The complex spatio-temporal regulation of the Drosophila myoblast attractant gene duf/kirre

A key early player in the regulation of myoblast fusion is the gene dumbfounded (duf, also known as kirre). Duf must be expressed, and function, in founder cells (FCs). A fixed number of FCs are chosen from a pool of equivalent myoblasts and serve to attract fusion-competent myoblasts (FCMs) to fuse with them to form a multinucleate muscle-fibre. The spatial and temporal regulation of duf expression and function are important and play a deciding role in choice of fibre number, location and perhaps size. We have used a combination of bioinformatics and functional enhancer deletion approaches to understand the regulation of duf. By transgenic enhancer-reporter deletion analysis of the duf regulatory region, we found that several distinct enhancer modules regulate duf expression in specific muscle founders of the embryo and the adult. In addition to existing bioinformatics tools, we used a new program for analysis of regulatory sequence, PhyloGibbs-MP, whose development was largely motivated by the requirements of this work. The results complement our deletion analysis by identifying transcription factors whose predicted binding regions match with our deletion constructs. Experimental evidence for the relevance of some of these TF binding sites comes from available ChIP-on-chip from the literature, and from our analysis of localization of myogenic transcription factors with duf enhancer reporter gene expression. Our results demonstrate the complex regulation in each founder cell of a gene that is expressed in all founder cells. They provide evidence for transcriptional control—both activation and repression—as an important player in the regulation of myoblast fusion. The set of enhancer constructs generated will be valuable in identifying novel trans-acting factor-binding sites and chromatin regulation during myoblast fusion in Drosophila. Our results and the bioinformatics tools developed provide a basis for the study of the transcriptional regulation of other complex genes.

The role of the visceral mesoderm in the development of the gastrointestinal tract

The gastrointestinal (GI) tract forms from the endoderm (which gives rise to the epithelium) and the mesoderm (which develops into the smooth muscle layer, the mesenchyme, and numerous other cell types). Much of what is known of GI development has been learned from studies of the endoderm and its derivatives, because of the importance of epithelial biology in understanding and treating human diseases. Although the necessity of epithelial-mesenchymal cross talk for GI development is uncontested, the role of the mesoderm remains comparatively less well understood. The transformation of the visceral mesoderm during development is remarkable; it differentiates from a very thin layer of cells into a complex tissue comprising smooth muscle cells, myofibroblasts, neurons, immune cells, endothelial cells, lymphatics, and extracellular matrix molecules, all contributing to the form and function of the digestive system. Understanding the molecular processes that govern the development of these cell types and elucidating their respective contribution to GI patterning could offer insight into the mechanisms that regulate cell fate decisions in the intestine, which has the unique property of rapid cell renewal for the maintenance of epithelial integrity. In reviewing evidence from both mammalian and nonmammalian models, we reveal the important role of the visceral mesoderm in the ontogeny of the GI tract.

Cardiac gene regulatory networks in Drosophila

The Drosophila system has proven a powerful tool to help unlock the regulatory processes that occur during specification and differentiation of the embryonic heart. In this review, we focus upon a temporal analysis of the molecular events that result in heart formation in Drosophila, with a particular emphasis upon how genomic and other cutting-edge approaches are being brought to bear upon the subject. We anticipate that systems-level approaches will contribute greatly to our comprehension of heart development and disease in the animal kingdom.

A systematic analysis of Tinman function reveals Eya and JAK-STAT signaling as essential regulators of muscle development

Nk-2 proteins are essential developmental regulators from flies to humans. In Drosophila, the family member tinman is the major regulator of cell fate within the dorsal mesoderm, including heart, visceral, and dorsal somatic muscle. To decipher Tinman's direct regulatory role, we performed a time course of ChIP-on-chip experiments, revealing a more prominent role in somatic muscle specification than previously anticipated. Through the combination of transgenic enhancer-reporter assays, colocalization studies, and phenotypic analyses, we uncovered two additional factors within this myogenic network: by activating eyes absent, Tinman's regulatory network extends beyond developmental stages and tissues where it is expressed; by regulating stat92E expression, Tinman modulates the transcriptional readout of JAK/STAT signaling. We show that this pathway is essential for somatic muscle development in Drosophila and for myotome morphogenesis in zebrafish. Taken together, these data uncover a conserved requirement for JAK/STAT signaling and an important component of the transcriptional network driving myogenesis.

The Xenopus MEF2 gene family: evidence of a role for XMEF2C in larval tendon development

MEF2 transcription factors are well-established regulators of muscle development. In this report, we describe the cloning of multiple splicing isoforms of the XMEF2A and XMEF2C encoding genes, differentially expressed during Xenopus development. Using whole-mount in situ hybridization, we found that the accumulation of XMEF2C mRNA in the tadpole stages was restricted to intersomitic regions and to the peripheral edges of hypaxial and cranial muscle masses in contrast to XMEF2A and XMEF2D, characterized by a continuous muscle cell expression. The XMEF2C positive cells express the bHLH transcription factor, Xscleraxis, known as a specific marker for tendons. Gain of function experiments revealed that the use of a hormone-inducible XMEF2C construct is able to induce Xscleraxis expression. Furthermore, XMEF2C specifically cooperates with Xscleraxis to induce tenascin C and betaig-h3, two genes preferentially expressed in Xenopus larval tendons. These findings 1) highlight a previously unappreciated and specific role for XMEF2C in tendon development and 2) identify a novel gene transactivation pathway where MEF2C cooperates with the bHLH protein, Xscleraxis, to activate specific gene expression.

Toward a systems-level understanding of developmental regulatory networks

Developmental regulatory networks constitute all the interconnections among molecular components that guide embryonic development. Developmental transcriptional regulatory networks (TRNs) are circuits of transcription factors and cis-acting DNA elements that control expression of downstream regulatory and effector genes. Developmental networks comprise functional subnetworks that are deployed sequentially in requisite spatiotemporal patterns. Here, we discuss integrative genomics approaches for elucidating TRNs, with an emphasis on those involved in Drosophila mesoderm development and mammalian embryonic stem cell maintenance and differentiation. As examples of regulatory subnetworks, we consider the transcriptional and signaling regulation of genes that interact to control cell morphology and migration. Finally, we describe integrative experimental and computational strategies for defining the entirety of molecular interactions underlying developmental regulatory networks.

Alternative requirements for Vestigial, Scalloped, and Dmef2 during muscle differentiation in Drosophila melanogaster

Vertebrate development requires the activity of the myocyte enhancer factor 2 (mef2) gene family for muscle cell specification and subsequent differentiation. Additionally, several muscle-specific functions of MEF2 family proteins require binding additional cofactors including members of the Transcription Enhancing Factor-1 (TEF-1) and Vestigial-like protein families. In Drosophila there is a single mef2 (Dmef2) gene as well single homologues of TEF-1 and vestigial-like, scalloped (sd), and vestigial (vg), respectively. To clarify the role(s) of these factors, we examined the requirements for Vg and Sd during Drosophila muscle specification. We found that both are required for muscle differentiation as loss of sd or vg leads to a reproducible loss of a subset of either cardiac or somatic muscle cells in developing embryos. This muscle requirement for Sd or Vg is cell specific, as ubiquitous overexpression of either or both of these proteins in muscle cells has a deleterious effect on muscle differentiation. Finally, using both in vitro and in vivo binding assays, we determined that Sd, Vg, and Dmef2 can interact directly. Thus, the muscle-specific phenotypes we have associated with Vg or Sd may be a consequence of alternative binding of Vg and/or Sd to Dmef2 forming alternative protein complexes that modify Dmef2 activity.

Gaussian process modelling of latent chemical species: applications to inferring transcription factor activities

MOTIVATION

Inference of latent chemical species in biochemical interaction networks is a key problem in estimation of the structure and parameters of the genetic, metabolic and protein interaction networks that underpin all biological processes. We present a framework for Bayesian marginalization of these latent chemical species through Gaussian process priors.

RESULTS

We demonstrate our general approach on three different biological examples of single input motifs, including both activation and repression of transcription. We focus in particular on the problem of inferring transcription factor activity when the concentration of active protein cannot easily be measured. We show how the uncertainty in the inferred transcription factor activity can be integrated out in order to derive a likelihood function that can be used for the estimation of regulatory model parameters. An advantage of our approach is that we avoid the use of a coarsegrained discretization of continuous time functions, which would lead to a large number of additional parameters to be estimated. We develop exact (for linear regulation) and approximate (for non-linear regulation) inference schemes, which are much more efficient than competing sampling-based schemes and therefore provide us with a practical toolkit for model-based inference.

AVAILABILITY

The software and data for recreating all the experiments in this paper is available in MATLAB from http://www.cs.man. ac.uk/~neill/gpsim.

Global gene expression analysis identifies PDEF transcriptional networks regulating cell migration during cancer progression

Prostate derived ETS factor (PDEF) is an ETS (epithelial-specific E26 transforming sequence) family member that has been identified as a potential tumor suppressor. In multiple invasive breast cancer cells, PDEF expression inhibits cell migration by preventing the acquisition of directional morphological polarity conferred by changes in cytoskeleton organization. In this study, microarray analysis was used to identify >200 human genes that displayed a common differential expression pattern in three invasive breast cancer cell lines after expression of exogenous PDEF protein. Gene ontology associations and data mining analysis identified focal adhesion, adherens junctions, cell adhesion, and actin cytoskeleton regulation as cell migration-associated interaction pathways significantly impacted by PDEF expression. Validation experiments confirmed the differential expression of four cytoskeleton-associated genes with known functional associations with these pathways: uPA, uPAR, LASP1, and VASP. Significantly, chromatin immunoprecipitation studies identified PDEF as a direct negative regulator of the metastasis-associated gene uPA and phenotypic rescue experiments demonstrate that exogenous urokinase plasminogen activator (uPA) expression can restore the migratory ability of invasive breast cancer cells expressing PDEF. Furthermore, immunofluorescence studies identify the subcellular relocalization of urokinase plasminogen activator receptor (uPAR), LIM and SH3 protein (LASP1), and vasodilator-stimulated protein (VASP) as a possible mechanism accounting for the loss of morphological polarity observed upon PDEF expression.

RNA interference screening in Drosophila primary cells for genes involved in muscle assembly and maintenance

To facilitate the genetic analysis of muscle assembly and maintenance, we have developed a method for efficient RNA interference (RNAi) in Drosophila primary cells using double-stranded RNAs (dsRNAs). First, using molecular markers, we confirm and extend the observation that myogenesis in primary cultures derived from Drosophila embryonic cells follows the same developmental course as that seen in vivo. Second, we apply this approach to analyze 28 Drosophila homologs of human muscle disease genes and find that 19 of them, when disrupted, lead to abnormal muscle phenotypes in primary culture. Third, from an RNAi screen of 1140 genes chosen at random, we identify 49 involved in late muscle differentiation. We validate our approach with the in vivo analyses of three genes. We find that Fermitin 1 and Fermitin 2, which are involved in integrin-containing adhesion structures, act in a partially redundant manner to maintain muscle integrity. In addition, we characterize CG2165, which encodes a plasma membrane Ca2+-ATPase, and show that it plays an important role in maintaining muscle integrity. Finally, we discuss how Drosophila primary cells can be manipulated to develop cell-based assays to model human diseases for RNAi and small-molecule screens.

Daughterless dictates Twist activity in a context-dependent manner during somatic myogenesis

Somatic myogenesis in Drosophila relies on the reiterative activity of the basic helix-loop-helix transcriptional regulator, Twist (Twi). How Twi directs multiple cell fate decisions over the course of mesoderm and muscle development is unclear. Previous work has shown that Twi is regulated by its dimerization partner: Twi homodimers activate genes necessary for somatic myogenesis, whereas Twi/Daughterless (Da) heterodimers lead to the repression of these genes. Here, we examine the nature of Twi/Da heterodimer repressive activity. Analysis of the Da protein structure revealed a Da repression (REP) domain, which is required for Twi/Da-mediated repression of myogenic genes, such as Dmef2, both in tissue culture and in vivo. This domain is crucial for the allocation of mesodermal cells to distinct fates, such as heart, gut and body wall muscle. By contrast, the REP domain is not required in vivo during later stages of myogenesis, even though Twi activity is required for muscles to achieve their final pattern and morphology. Taken together, we present evidence that the repressive activity of the Twi/Da dimer is dependent on the Da REP domain and that the activity of the REP domain is sensitive to tissue context and developmental timing.

Transcription factor functionality and transcription regulatory networks

Now that numerous high-quality complete genome sequences are available, many efforts are focusing on the "second genomic code", namely the code that determines how the precise temporal and spatial expression of each gene in the genome is achieved. In this regard, the elucidation of transcription regulatory networks that describe combined transcriptional circuits for an organism of interest has become valuable to our understanding of gene expression at a systems level. Such networks describe physical and regulatory interactions between transcription factors (TFs) and the target genes they regulate under different developmental, physiological, or pathological conditions. The mapping of high-quality transcription regulatory networks depends not only on the accuracy of the experimental or computational method chosen, but also relies on the quality of TF predictions. Moreover, the total repertoire of TFs is not only determined by the protein-coding capacity of the genome, but also by different protein properties, including dimerization, co-factor interactions and post-translational modifications. Here, we discuss the factors that influence TF functionality and, hence, the functionality of the networks in which they operate.

CF2 activity and enhancer integration are required for proper muscle gene expression in Drosophila

The creation of the contractile apparatus in muscle involves the co-activation of a group of genes encoding muscle-specific proteins and the production of high levels of protein in a short period of time. We have studied the transcriptional control of six Drosophila muscle genes that have similar expression profiles and we have compared these mechanisms with those employed to control the distinct expression profiles of other Drosophila genes. The regulatory elements controlling the transcription of co-expressed muscle genes share an Upstream Regulatory Element and an Intronic Regulatory Element. Moreover, similar clusters of MEF2 and CF2 binding sites are present in these elements. Here, we demonstrate that CF2 depletion alters the relative expression of thin and thick filament components. We propose that the appropriate rapid gene expression responses during muscle formation and the maintenance of each muscle type is guaranteed in Drosophila by equivalent duplicate enhancer-like elements. This mechanism may be exceptional and restricted to muscle genes, reflecting the specific requirement to mediate rapid muscle responses. However, it may also be a more general mechanism to control the correct levels of gene expression during development in each cell type.

Erk kinases link pre-B cell receptor signaling to transcriptional events required for early B cell expansion

The pre-B cell receptor (pre-BCR) plays a crucial role in the development of immature B cells. Although certain aspects of proximal pre-BCR signaling have been studied, the intermediate signal transducers and the distal transcription modulators are poorly characterized. Here, we demonstrate that deletion of both Erk1 and Erk2 kinases was associated with defective pre-BCR-mediated cell expansion as well as a block in the transition of pro-B to pre-B cells. Phosphorylation of transcription factors Elk1 and CREB was mediated by Erk, and a dominant-negative mutation in the Erk-mediated phosphorylation sites of Elk1 or CREB suppressed pre-BCR-mediated cell expansion as well as expression of genes including Myc, which is involved in the cell-cycle progression. Together, our results identify a crucial role for Erk kinases in regulating B cell development by initiating transcriptional regulatory network and thereby pre-BCR-mediated cell expansion.

Myocyte enhancer factor 2 and chorion factor 2 collaborate in activation of the myogenic program in Drosophila

The process of myogenesis requires the coordinated activation of many structural genes whose products are required for myofibril assembly, function, and regulation. Although numerous reports have documented the importance of the myogenic regulator myocyte enhancer factor 2 (MEF2) in muscle differentiation, the interaction of MEF2 with cofactors is critical to the realization of muscle fate. We identify here a genomic region required for full MEF2-mediated activation of actin gene expression in Drosophila, and we identify the zinc finger transcriptional regulator chorion factor 2 (CF2) as a factor functioning alongside MEF2 via this region. Furthermore, although both MEF2 and CF2 can individually activate actin gene expression, we demonstrate that these two factors collaborate in regulating the Actin57B target gene in vitro and in vivo. More globally, MEF2 and CF2 synergistically activate the enhancers of a number of muscle-specific genes, and loss of CF2 function in vivo results in reductions in the levels of several muscle structural gene transcripts. These findings validate a general importance of CF2 alongside MEF2 as a critical regulator of the myogenic program, identify a new regulator functioning with MEF2 to control cell fate, and provide insight into the network of regulatory events that shape the developing musculature.

The conserved microRNA miR-8 tunes atrophin levels to prevent neurodegeneration in Drosophila

microRNAs (miRNAs) bind to specific messenger RNA targets to posttranscriptionally modulate their expression. Understanding the regulatory relationships between miRNAs and targets remains a major challenge. Many miRNAs reduce expression of their targets to inconsequential levels. It has also been proposed that miRNAs might adjust target expression to an optimal level. Here we analyze the consequences of mutating the conserved miRNA miR-8 in Drosophila. We identify atrophin as a direct target of miR-8. miR-8 mutant phenotypes are attributable to elevated atrophin activity, resulting in elevated apoptosis in the brain and in behavioral defects. Reduction of atrophin levels in miR-8-expressing cells to below the level generated by miR-8 regulation is detrimental, providing evidence for a "tuning target" relationship between them. Drosophila atrophin is related to the atrophin family of mammalian transcriptional regulators, implicated in the neurodegenerative disorder DRPLA. The regulatory relationship between miR-8 and atrophin orthologs is conserved in mammals.

mef2 activity levels differentially affect gene expression during Drosophila muscle development

Cell differentiation is controlled by key transcription factors, and a major question is how they orchestrate cell-type-specific genetic programs. Muscle differentiation is a well studied paradigm in which the conserved Mef2 transcription factor plays a pivotal role. Recent genomic studies have identified a large number of mef2-regulated target genes with distinct temporal expression profiles during Drosophila myogenesis. However, the question remains as to how a single transcription factor can control such diverse patterns of gene expression. In this study we used a strategy combining genomics and developmental genetics to address this issue in vivo during Drosophila muscle development. We found that groups of mef2-regulated genes respond differently to changes in mef2 activity levels: some require higher levels for their expression than others. Furthermore, this differential requirement correlates with when the gene is first expressed during the muscle differentiation program. Genes that require higher levels are activated later. These results implicate mef2 in the temporal regulation of muscle gene expression, and, consistent with this, we show that changes in mef2 activity levels can alter the start of gene expression in a predictable manner. Together these results indicate that Mef2 is not an all-or-none regulator; rather, its action is more subtle, and levels of its activity are important in the differential expression of muscle genes. This suggests a route by which mef2 can orchestrate the muscle differentiation program and contribute to the stringent regulation of gene expression during myogenesis.

Genome-wide view of cell fate specification: ladybird acts at multiple levels during diversification of muscle and heart precursors

Correct diversification of cell types during development ensures the formation of functional organs. The evolutionarily conserved homeobox genes from ladybird/Lbx family were found to act as cell identity genes in a number of embryonic tissues. A prior genetic analysis showed that during Drosophila muscle and heart development ladybird is required for the specification of a subset of muscular and cardiac precursors. To learn how ladybird genes exert their cell identity functions we performed muscle and heart-targeted genome-wide transcriptional profiling and a chromatin immunoprecipitation (ChIP)-on-chip search for direct Ladybird targets. Our data reveal that ladybird not only contributes to the combinatorial code of transcription factors specifying the identity of muscle and cardiac precursors, but also regulates a large number of genes involved in setting cell shape, adhesion, and motility. Among direct ladybird targets, we identified bric-a-brac 2 gene as a new component of identity code and inflated encoding alphaPS2-integrin playing a pivotal role in cell-cell interactions. Unexpectedly, ladybird also contributes to the regulation of terminal differentiation genes encoding structural muscle proteins or contributing to muscle contractility. Thus, the identity gene-governed diversification of cell types is a multistep process involving the transcriptional control of genes determining both morphological and functional properties of cells.

The Him gene reveals a balance of inputs controlling muscle differentiation in Drosophila

Tissue development requires the controlled regulation of cell-differentiation programs. In muscle, the Mef2 transcription factor binds to and activates the expression of many genes and has a major positive role in the orchestration of differentiation [1–4]. However, little is known about how Mef2 activity is regulated in vivo during development. Here, we characterize a gene, Holes in muscle (Him), which our results indicate is part of this control in Drosophila. Him expression rapidly declines as embryonic muscle differentiates, and consistent with this, Him overexpression inhibits muscle differentiation. This inhibitory effect is suppressed by mef2, implicating Him in the mef2 pathway. We then found that Him downregulates the transcriptional activity of Mef2 in both cell culture and in vivo. Furthermore, Him protein binds Groucho, a conserved, transcriptional corepressor, through a WRPW motif and requires this motif and groucho function to inhibit both muscle differentiation and Mef2 activity during development. Together, our results identify a mechanism that can inhibit muscle differentiation in vivo. We conclude that a balance of positive and negative inputs, including Mef2, Him, and Groucho, controls muscle differentiation during Drosophila development and suggest that one outcome is to hold developing muscle cells in a state with differentiation genes poised to be expressed.

RNA polymerase stalling at developmental control genes in the Drosophila melanogaster embryo

It is widely assumed that the key rate-limiting step in gene activation is the recruitment of RNA polymerase II (Pol II) to the core promoter. Although there are well-documented examples in which Pol II is recruited to a gene but stalls, a general role for Pol II stalling in development has not been established. We have carried out comprehensive Pol II chromatin immunoprecipitation microarray (ChIP-chip) assays in Drosophila embryos and identified three distinct Pol II binding behaviors: active (uniform binding across the entire transcription unit), no binding, and stalled (binding at the transcription start site). The notable feature of the approximately 10% genes that are stalled is that they are highly enriched for developmental control genes, which are either repressed or poised for activation during later stages of embryogenesis. We propose that Pol II stalling facilitates rapid temporal and spatial changes in gene activity during development.

Reliable prediction of regulator targets using 12 Drosophila genomes

Gene expression is regulated pre- and post-transcriptionally via cis-regulatory DNA and RNA motifs. Identification of individual functional instances of such motifs in genome sequences is a major goal for inferring regulatory networks yet has been hampered due to the motifs' short lengths that lead to many chance matches and poor signal-to-noise ratios. In this paper, we develop a general methodology for the comparative identification of functional motif instances across many related species, using a phylogenetic framework that accounts for the evolutionary relationships between species, allows for motif movements, and is robust against missing data due to artifacts in sequencing, assembly, or alignment. We also provide a robust statistical framework for evaluating motif confidence, which enables us to translate evolutionary conservation into a confidence measure for each motif instance, correcting for varying motif length, composition, and background conservation of the target regions. We predict targets of fly transcription factors and miRNAs in alignments of 12 recently sequenced Drosophila species. When compared to extensive genome-wide experimental data, predicted targets are of high quality, matching and surpassing ChIP-chip microarrays and recovering miRNA targets with high sensitivity. The resulting regulatory network suggests significant redundancy between pre- and post-transcriptional regulation of gene expression.

Temporal ChIP-on-chip reveals Biniou as a universal regulator of the visceral muscle transcriptional network

Smooth muscle plays a prominent role in many fundamental processes and diseases, yet our understanding of the transcriptional network regulating its development is very limited. The FoxF transcription factors are essential for visceral smooth muscle development in diverse species, although their direct regulatory role remains elusive. We present a transcriptional map of Biniou (a FoxF transcription factor) and Bagpipe (an Nkx factor) activity, as a first step to deciphering the developmental program regulating Drosophila visceral muscle development. A time course of chromatin immunoprecipitatation followed by microarray analysis (ChIP-on-chip) experiments and expression profiling of mutant embryos reveal a dynamic map of in vivo bound enhancers and direct target genes. While Biniou is broadly expressed, it regulates enhancers driving temporally and spatially restricted expression. In vivo reporter assays indicate that the timing of Biniou binding is a key trigger for the time span of enhancer activity. Although bagpipe and biniou mutants phenocopy each other, their regulatory potential is quite different. This network architecture was not apparent from genetic studies, and highlights Biniou as a universal regulator in all visceral muscle, regardless of its developmental origin or subsequent function. The regulatory connection of a number of Biniou target genes is conserved in mice, suggesting an ancient wiring of this developmental program.

REDfly 2.0: an integrated database of cis-regulatory modules and transcription factor binding sites in Drosophila

The identification and study of the cis-regulatory elements that control gene expression are important areas of biological research, but few resources exist to facilitate large-scale bioinformatics studies of cis-regulation in metazoan species. Drosophila melanogaster, with its well-annotated genome, exceptional resources for comparative genomics and long history of experimental studies of transcriptional regulation, represents the ideal system for regulatory bioinformatics. We have merged two existing Drosophila resources, the REDfly database of cis-regulatory modules and the FlyReg database of transcription factor binding sites (TFBSs), into a single integrated database containing extensive annotation of empirically validated cis-regulatory modules and their constituent binding sites. With the enhanced functionality made possible through this integration of TFBS data into REDfly, together with additional improvements to the REDfly infrastructure, we have constructed a one-stop portal for Drosophila cis-regulatory data that will serve as a powerful resource for both computational and experimental studies of transcriptional regulation. REDfly is freely accessible at http://redfly.ccr.buffalo.edu.

Collier transcription in a single Drosophila muscle lineage: the combinatorial control of muscle identity

Specification of muscle identity in Drosophila is a multistep process: early positional information defines competence groups termed promuscular clusters, from which muscle progenitors are selected, followed by asymmetric division of progenitors into muscle founder cells (FCs). Each FC seeds the formation of an individual muscle with morphological and functional properties that have been proposed to reflect the combination of transcription factors expressed by its founder. However, it is still unclear how early patterning and muscle-specific differentiation are linked. We addressed this question, using Collier (Col; also known as Knot) expression as both a determinant and read-out of DA3 muscle identity. Characterization of the col upstream region driving DA3 muscle specific expression revealed the existence of three separate phases of cis-regulation, correlating with conserved binding sites for different mesodermal transcription factors. Examination of col transcription in col and nautilus (nau) loss-of-function and gain-of-function conditions showed that both factors are required for col activation in the ;naïve' myoblasts that fuse with the DA3 FC, thereby ensuring that all DA3 myofibre nuclei express the same identity programme. Together, these results indicate that separate sets of cis-regulatory elements control the expression of identity factors in muscle progenitors and myofibre nuclei and directly support the concept of combinatorial control of muscle identity.

Discovery of functional elements in 12 Drosophila genomes using evolutionary signatures

Sequencing of multiple related species followed by comparative genomics analysis constitutes a powerful approach for the systematic understanding of any genome. Here, we use the genomes of 12 Drosophila species for the de novo discovery of functional elements in the fly. Each type of functional element shows characteristic patterns of change, or 'evolutionary signatures', dictated by its precise selective constraints. Such signatures enable recognition of new protein-coding genes and exons, spurious and incorrect gene annotations, and numerous unusual gene structures, including abundant stop-codon readthrough. Similarly, we predict non-protein-coding RNA genes and structures, and new microRNA (miRNA) genes. We provide evidence of miRNA processing and functionality from both hairpin arms and both DNA strands. We identify several classes of pre- and post-transcriptional regulatory motifs, and predict individual motif instances with high confidence. We also study how discovery power scales with the divergence and number of species compared, and we provide general guidelines for comparative studies.

ChIP-on-chip protocol for genome-wide analysis of transcription factor binding in Drosophila melanogaster embryos

This protocol describes a method to detect in vivo associations between proteins and DNA in developing Drosophila embryos. It combines formaldehyde crosslinking and immunoprecipitation of protein-bound sequences with genome-wide analysis using microarrays. After crosslinking, nuclei are enriched using differential centrifugation and the chromatin is sheared by sonication. Antibodies specifically recognizing wild-type protein or, alternatively, a genetically encoded epitope tag are used to enrich for specifically bound DNA sequences. After purification and polymerase chain reaction-based amplification, the samples are fluorescently labeled and hybridized to genomic tiling microarrays. This protocol has been successfully used to study different tissue-specific transcription factors, and is generally applicable to in vivo analysis of any DNA-binding proteins in Drosophila embryos. The full protocol, including the collection of embryos and the collection of raw microarray data, can be completed within 10 days.

MEF2: a central regulator of diverse developmental programs

The myocyte enhancer factor 2 (MEF2) transcription factor acts as a lynchpin in the transcriptional circuits that control cell differentiation and organogenesis. The spectrum of genes activated by MEF2 in different cell types depends on extracellular signaling and on co-factor interactions that modulate MEF2 activity. Recent studies have revealed MEF2 to form an intimate partnership with class IIa histone deacetylases, which together function as a point of convergence of multiple epigenetic regulatory mechanisms. We review the myriad roles of MEF2 in development and the mechanisms through which it couples developmental, physiological and pathological signals with programs of cell-specific transcription.

Mef2s are required for thick filament formation in nascent muscle fibres

During skeletal muscle differentiation, the actomyosin motor is assembled into myofibrils, multiprotein machines that generate and transmit force to cell ends. How expression of muscle proteins is coordinated to build the myofibril is unknown. Here we show that zebrafish Mef2d and Mef2c proteins are required redundantly for assembly of myosin-containing thick filaments in nascent muscle fibres, but not for the earlier steps of skeletal muscle fibre differentiation, elongation, fusion or thin filament gene expression. mef2d mRNA and protein is present in myoblasts, whereas mef2c expression commences in muscle fibres. Knockdown of both Mef2s with antisense morpholino oligonucleotides or in mutant fish blocks muscle function and prevents sarcomere assembly. Cell transplantation and heat-shock-driven rescue reveal a cell-autonomous requirement for Mef2 within fibres. In nascent fibres, Mef2 drives expression of genes encoding thick, but not thin, filament proteins. Among genes analysed, myosin heavy and light chains and myosin-binding protein C require Mef2 for normal expression, whereas actin, tropomyosin and troponin do not. Our findings show that Mef2 controls skeletal muscle formation after terminal differentiation and define a new maturation step in vertebrate skeletal muscle development at which thick filament gene expression is controlled.

A core transcriptional network for early mesoderm development in Drosophila melanogaster

Embryogenesis is controlled by large gene-regulatory networks, which generate spatially and temporally refined patterns of gene expression. Here, we report the characteristics of the regulatory network orchestrating early mesodermal development in the fruitfly Drosophila, where the transcription factor Twist is both necessary and sufficient to drive development. Through the integration of chromatin immunoprecipitation followed by microarray analysis (ChIP-on-chip) experiments during discrete time periods with computational approaches, we identified >2000 Twist-bound cis-regulatory modules (CRMs) and almost 500 direct target genes. Unexpectedly, Twist regulates an almost complete cassette of genes required for cell proliferation in addition to genes essential for morophogenesis and cell migration. Twist targets almost 25% of all annotated Drosophila transcription factors, which may represent the entire set of regulators necessary for the early development of this system. By combining in vivo binding data from Twist, Mef2, Tinman, and Dorsal we have constructed an initial transcriptional network of early mesoderm development. The network topology reveals extensive combinatorial binding, feed-forward regulation, and complex logical outputs as prevalent features. In addition to binary activation and repression, we suggest that Twist binds to almost all mesodermal CRMs to provide the competence to integrate inputs from more specialized transcription factors.

From functional genomics to systems biology

This review discusses the talks presented at the third EMBL Biennial Symposium, From functional genomics to systems biology, held in Heidelberg, Germany, 14-17 October 2006. Current issues and trends in various subfields of functional genomics and systems biology are considered, including analysis of regulatory elements, signalling networks, transcription networks, protein-protein interaction networks, genetic interaction networks, medical applications of DNA microarrays, and metagenomics. Several technological advances in the fields of DNA microarrays, identification of regulatory elements in the genomes of higher eukaryotes, and MS for detection of protein interactions are introduced. Major directions of future systems biology research are also discussed.

Genomic characterization of Gli-activator targets in sonic hedgehog-mediated neural patterning

Sonic hedgehog (Shh) acts as a morphogen to mediate the specification of distinct cell identities in the ventral neural tube through a Gli-mediated (Gli1-3) transcriptional network. Identifying Gli targets in a systematic fashion is central to the understanding of the action of Shh. We examined this issue in differentiating neural progenitors in mouse. An epitope-tagged Gli-activator protein was used to directly isolate cis-regulatory sequences by chromatin immunoprecipitation (ChIP). ChIP products were then used to screen custom genomic tiling arrays of putative Hedgehog (Hh) targets predicted from transcriptional profiling studies, surveying 50-150 kb of non-transcribed sequence for each candidate. In addition to identifying expected Gli-target sites, the data predicted a number of unreported direct targets of Shh action. Transgenic analysis of binding regions in Nkx2.2, Nkx2.1 (Titf1) and Rab34 established these as direct Hh targets. These data also facilitated the generation of an algorithm that improved in silico predictions of Hh target genes. Together, these approaches provide significant new insights into both tissue-specific and general transcriptional targets in a crucial Shh-mediated patterning process.

Identification and mapping of protein kinase A binding sites in the costameric protein myospryn

Recently we identified a novel target gene of MEF2A named myospryn that encodes a large, muscle-specific, costamere-restricted alpha-actinin binding protein. Myospryn belongs to the tripartite motif (TRIM) superfamily of proteins and was independently identified as a dysbindin-interacting protein. Dysbindin is associated with alpha-dystrobrevin, a component of the dystrophin-glycoprotein complex (DGC) in muscle. Apart from these initial findings little else is known regarding the potential function of myospryn in striated muscle. Here we reveal that myospryn is an anchoring protein for protein kinase A (PKA) (or AKAP) whose closest homolog is AKAP12, also known as gravin/AKAP250/SSeCKS. We demonstrate that myospryn co-localizes with RII alpha, a type II regulatory subunit of PKA, at the peripheral Z-disc/costameric region in striated muscle. Myospryn interacts with RII alpha and this scaffolding function has been evolutionarily conserved as the zebrafish ortholog also interacts with PKA. Moreover, myospryn serves as a substrate for PKA. These findings point to localized PKA signaling at the muscle costamere.

CoCo: a web application to display, store and curate ChIP-on-chip data integrated with diverse types of gene expression data

MOTIVATION

CoCo, ChIP-on-Chip online, is an open-source web application that supports the annotation and curation of regulatory regions and associated target genes discovered in ChIP-on-chip experiments. CoCo integrates ChIP-on-chip results with diverse types of gene expression data (expression profiling, in situ hybridization) and displays them within a genomic context. Regulatory relationships between the transcription factor-bound regions and putative target genes can be stored and expanded throughout different sessions.

AVAILABILITY

http://furlonglab.embl.de/methods/tools/coco.

Multiple upstream modules regulate zebrafish myf5 expression

Background

Myf5 is one member of the basic helix-loop-helix family of transcription factors, and it functions as a myogenic factor that is important for the specification and differentiation of muscle cells. The expression of myf5 is somite- and stage-dependent during embryogenesis through a delicate regulation. However, this complex regulatory mechanism of myf5 is not clearly understood.

Results

We isolated a 156-kb bacterial artificial chromosome clone that includes an upstream 80-kb region and a downstream 70-kb region of zebrafish myf5 and generated a transgenic line carrying this 156-kb segment fused to a green fluorescent protein (GFP) reporter gene. We find strong GFP expression in the most rostral somite and in the presomitic mesoderm during segmentation stages, similar to endogenous myf5 expression. Later, the GFP signals persist in caudal somites near the tail bud but are down-regulated in the older, rostral somites. During the pharyngula period, we detect GFP signals in pectoral fin buds, dorsal rostral myotomes, hypaxial myotomes, and inferior oblique and superior oblique muscles, a pattern that also corresponds well with endogenous myf5 transcripts. To characterize the specific upstream cis-elements that regulate this complex and dynamic expression pattern, we also generated several transgenic lines that harbor various lengths within the upstream 80-kb segment. We find that (1) the -80 kb/-9977 segment contains a fin and cranial muscle element and a notochord repressor; (2) the -9977/-6213 segment contains a strong repressive element that does not include the notochord-specific repressor; (3) the -6212/-2938 segment contains tissue-specific elements for bone and spinal cord; (4) the -2937/-291 segment contains an eye enhancer, and the -2937/-2457 segment is required for notochord and myocyte expression; and (5) the -290/-1 segment is responsible for basal transcription in somites and the presomitic mesoderm.

Conclusion

We suggest that the cell lineage-specific expression of myf5 is delicately orchestrated by multiple modules within the distal upstream region. This study provides an insight to understand the molecular control of myf5 and myogenesis in the zebrafish.

Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers

Restriction site associated DNA (RAD) tags are a genome-wide representation of every site of a particular restriction enzyme by short DNA tags. Most organisms segregate large numbers of DNA sequence polymorphisms that disrupt restriction sites, which allows RAD tags to serve as genetic markers spread at a high density throughout the genome. Here, we demonstrate the applicability of RAD markers for both individual and bulk-segregant genotyping. First, we show that these markers can be identified and typed on pre-existing microarray formats. Second, we present a method that uses RAD marker DNA to rapidly produce a low-cost microarray genotyping resource that can be used to efficiently identify and type thousands of RAD markers. We demonstrate the utility of the former approach by using a tiling path array for the fruit fly to map a recombination breakpoint, and the latter approach by creating and using an enriched RAD marker array for the threespine stickleback. The high number of RAD markers enabled localization of a previously identified region, as well as a second region also associated with the lateral plate phenotype. Taken together, our results demonstrate that RAD markers, and the method to develop a RAD marker microarray resource, allow high-throughput, high-resolution genotyping in both model and nonmodel systems.

Mes2, a MADF-containing transcription factor essential for Drosophila development

The development of the Drosophila mesoderm is initiated by the basic helix-loop-helix transcription factor twist. We identified a gene encoding a putative transcription factor, mes2, in a screen for essential mesoderm-expressed genes that function downstream of twist. Mes2 protein belongs to a family of 48 Drosophila proteins containing MADF domains. MADF domains exist in worms, flies, and fish. Mes2 is a nuclear protein first produced in trunk and head mesoderm during late gastrulation. At later embryonic stages, mes2 is expressed in glia of the central and peripheral nervous systems, and in tissues derived from the head mesoderm. We have identified a null mutation of mes2 that leads to developmental arrest in first instar larvae. Increased production of Mes2 in multiple embryonic and larval tissues almost always causes lethality. The ubiquitous or epidermal misexpression of mes2 in the embryo causes a dramatic loss of epidermal integrity resulting in the failure of dorsal closure. Our data show that the precise regulation of mes2 expression is critical for normal development in Drosophila and implicate Mes2 in the regulation of essential target genes.

Genomics and development: Taking developmental biology to new heights

The 2006 Arolla meeting brought together scientists from around the globe to discuss how genomic scale analyses can enhance progress in understanding developmental biology.

The bHLH transcription factor Hand is regulated by Alk in the Drosophila embryonic gut

During embryonic development the midgut visceral muscle is formed by fusion of cells within the visceral mesoderm, a process initiated by the specification of a specialised cell type, the founder cell, within this tissue. Activation of the receptor tyrosine kinase Anaplastic lymphoma kinase (Alk) in the developing visceral muscle of Drosophila melanogaster initiates a signal transduction pathway required for muscle fusion. In this paper, we have investigated downstream components which are regulated by this novel signalling pathway. Here we show that Alk-mediated signal transduction drives the expression of the bHLH transcription factor Hand in vivo. Loss of Alk function results in a complete lack of Hand expression in this tissue, whereas Alk gain of function results in an expansion of Hand expression. Finally, we have investigated the process of muscle fusion in the gut of Hand mutant animals and can find no obvious defects in this process, suggesting that Hand is not critical for visceral muscle fusion per se.