Genome-wide mapping of in vivo targets of the Drosophila transcription factor Kruppel

Activation of Ftz-F1-Responsive Genes through Ftz/Ftz-F1 Dependent Enhancers

The orphan nuclear receptor Ftz-F1 is expressed in all somatic nuclei in Drosophila embryos, but mutations result in a pair-rule phenotype. This was explained by the interaction of Ftz-F1 with the homeodomain protein Ftz that is expressed in stripes in the primordia of segments missing in either ftz-f1 or ftz mutants. Ftz-F1 and Ftz were shown to physically interact and coordinately activate the expression of ftz itself and engrailed by synergistic binding to composite Ftz-F1/Ftz binding sites. However, attempts to identify additional target genes on the basis of Ftz-F1/ Ftz binding alone has met with only limited success. To discern rules for Ftz-F1 target site selection in vivo and to identify additional target genes, a microarray analysis was performed comparing wildtype and ftz-f1 mutant embryos. Ftz-F1-responsive genes most highly regulated included engrailed and nine additional genes expressed in patterns dependent on both ftz and ftz-f1. Candidate enhancers for these genes were identified by combining BDTNP Ftz ChIP-chip data with a computational search for Ftz-F1 binding sites. Of eight enhancer reporter genes tested in transgenic embryos, six generated expression patterns similar to the corresponding endogenous gene and expression was lost in ftz mutants. These studies identified a new set of Ftz-F1 targets, all of which are co-regulated by Ftz. Comparative analysis of enhancers containing Ftz/Ftz-F1 binding sites that were or were not bona fide targets in vivo suggested that GAF negatively regulates enhancers that contain Ftz/Ftz-F1 binding sites but are not actually utilized. These targets include other regulatory factors as well as genes involved directly in morphogenesis, providing insight into how pair-rule genes establish the body pattern.

Constraints, independence, and evolution of thermal plasticity: probing genetic architecture of long- and short-term thermal acclimation

Seasonal and daily thermal variation can limit species distributions because of physiological tolerances. Low temperatures are particularly challenging for ectotherms, which use both basal thermotolerance and acclimation, an adaptive plastic response, to mitigate thermal stress. Both basal thermotolerance and acclimation are thought to be important for local adaptation and persistence in the face of climate change. However, the evolutionary independence of basal and plastic tolerances remains unclear. Acclimation can occur over longer (seasonal) or shorter (hours to days) time scales, and the degree of mechanistic overlap is unresolved. Using a midlatitude population of Drosophila melanogaster, we show substantial heritable variation in both short- and long-term acclimation. Rapid cold hardening (short-term plasticity) and developmental acclimation (long-term plasticity) are positively correlated, suggesting shared mechanisms. However, there are independent components of these traits, because developmentally acclimated flies respond positively to short-term acclimation. A strong negative correlation between basal cold tolerance and developmental acclimation suggests that basal cold tolerance may constrain developmental acclimation, whereas a weaker negative correlation between basal cold tolerance and short-term acclimation suggests less constraint. Using genome-wide association mapping, we show the genetic architecture of rapid cold hardening and developmental acclimation responses are nonoverlapping at the SNP and corresponding gene level. However, genes associated with each trait share functional similarities, including genes involved in apoptosis and autophagy, cytoskeletal and membrane structural components, and ion binding and transport. These results indicate substantial opportunity for short-term and long-term acclimation responses to evolve separately from each other and for short-term acclimation to evolve separately from basal thermotolerance.

CisMiner: genome-wide in-silico cis-regulatory module prediction by fuzzy itemset mining

Eukaryotic gene control regions are known to be spread throughout non-coding DNA sequences which may appear distant from the gene promoter. Transcription factors are proteins that coordinately bind to these regions at transcription factor binding sites to regulate gene expression. Several tools allow to detect significant co-occurrences of closely located binding sites (cis-regulatory modules, CRMs). However, these tools present at least one of the following limitations: 1) scope limited to promoter or conserved regions of the genome; 2) do not allow to identify combinations involving more than two motifs; 3) require prior information about target motifs. In this work we present CisMiner, a novel methodology to detect putative CRMs by means of a fuzzy itemset mining approach able to operate at genome-wide scale. CisMiner allows to perform a blind search of CRMs without any prior information about target CRMs nor limitation in the number of motifs. CisMiner tackles the combinatorial complexity of genome-wide cis-regulatory module extraction using a natural representation of motif combinations as itemsets and applying the Top-Down Fuzzy Frequent- Pattern Tree algorithm to identify significant itemsets. Fuzzy technology allows CisMiner to better handle the imprecision and noise inherent to regulatory processes. Results obtained for a set of well-known binding sites in the S. cerevisiae genome show that our method yields highly reliable predictions. Furthermore, CisMiner was also applied to putative in-silico predicted transcription factor binding sites to identify significant combinations in S. cerevisiae and D. melanogaster, proving that our approach can be further applied genome-wide to more complex genomes. CisMiner is freely accesible at: http://genome2.ugr.es/cisminer. CisMiner can be queried for the results presented in this work and can also perform a customized cis-regulatory module prediction on a query set of transcription factor binding sites provided by the user.

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.

Chromatin immunoprecipitation reveals a novel role for the Drosophila SoxNeuro transcription factor in axonal patterning

In all metazoans, the expression of group B HMG domain Sox transcription factors is associated with the earliest stages of CNS development. In Drosophila, SoxNeuro (SoxN) is involved in dorso-ventral patterning of the neuroectoderm, and in the formation and segregation of neuroblasts. In this report, we show that SoxN expression persists in a subset of neurons and glial cells of the ventral nerve cord at embryonic stages 15/16. In an attempt to address SoxN function in late stages of CNS development, we have used a chromatin immunoprecipitation approach to isolate genomic regions bound in vivo by SoxN. We identified several genes involved in the regulation of axon scaffolding as potential direct target genes of SoxN, including beat1a, semaphorin2a, fasciclin2, longitudinal lacking and tailup/islet. We present genetic evidence for a direct involvement of SoxN in axonal patterning. Indeed, overexpressing a transcriptionally hyperactive mutated SoxN protein in neurons results in specific defects in axon scaffolding, which are also observed in transheterozygous combinations of SoxN null mutation and mutations in its target genes.

Mapping Dmef2-binding regulatory modules by using a ChIP-enriched in silico targets approach

Mapping the regulatory modules to which transcription factors bind in vivo is a key step toward understanding of global gene expression programs. We have developed a chromatin immunoprecipitation (ChIP)-chip strategy for identifying factor-specific regulatory regions acting in vivo. This method, called the ChIP-enriched in silico targets (ChEST) approach, combines immunoprecipitation of cross-linked protein-DNA complexes (X-ChIP) with in silico prediction of targets and generation of computed DNA microarrays. We report the use of ChEST in Drosophila to identify several previously unknown targets of myocyte enhancer factor 2 (MEF2), a key regulator of myogenic differentiation. Our approach was validated by demonstrating that the identified sequences act as enhancers in vivo and are able to drive reporter gene expression specifically in MEF2-positive muscle cells. Presented here, the ChEST strategy was originally designed to identify regulatory modules in Drosophila, but it can be adapted for any sequenced and annotated genome.

Prediction of cis-regulatory elements using binding site matrices--the successes, the failures and the reasons for both

Protein-DNA interactions control many aspects of animal development and cellular responses to the environment. Although profiling of individual transcription factor binding sites is not a reliable guide for predicting the position of cis-regulatory elements in large genomes, modelling the evolution and the organization of regulatory elements has provided enough information to make some successful predictions. For vertebrate genomes, the field is limited by the lack of sufficient experimental data upon which to build reliable models. Nonetheless, a combination of experimental, computational and comparative data is likely to reveal aspects of complex regulatory networks in vertebrates, just as it has already done for simple eukaryotic genomes.

Identification and analysis of vnd/NK-2 homeodomain binding sites in genomic DNA

Vnd/NK-2 homeodomain affinity column chromatography was used to purify Drosophila DNA fragments bound by the vnd/NK-2 homeodomain. Sequencing the selected genomic DNA fragments led to the identification of 77 Drosophila DNA fragments that were grouped into 42 vnd/NK-2 homeodomain-binding loci. Most loci were within upstream or intronic regions, especially first introns. Nineteen of the Drosophila DNA fragments cloned correspond to one locus, termed Clone A, which is 312 bp in length and contains five vnd/NK-2 homeodomain core consensus binding sites, 5'-AAGTG, and is part of the first intron of the Beadex gene. We further analyzed the interactions between Clone A and vnd/NK-2 homeodomain protein by mobility-shift assay, DNase I footprinting, methylation interference, and ethylation interference. The DNase I footprinting analysis of Clone A with vnd/NK-2 homeodomain protein revealed three strong binding sites and one weak binding site between 15 and 130 bp of Clone A. We also analyzed binding of the vnd/NK-2 homeodomain to the 5'-flanking sequence of vnd/NK-2 genomic DNA. The DNase I footprinting result showed that there are two strong binding sites and five weak binding sites in the fragment between -385 and -675 bp from the transcription start site of the vnd/NK-2 gene.