Genome-wide identification of the regulatory targets of a transcription factor using biochemical characterization and computational genomic analysis

Characterization of Schistosome Sox Genes and Identification of a Flatworm Class of Sox Regulators

Schistosome helminths infect over 200 million people across 78 countries and are responsible for nearly 300,000 deaths annually. However, our understanding of basic genetic pathways crucial for schistosome development is limited. The sex determining region Y-box 2 (Sox2) protein is a Sox B type transcriptional activator that is expressed prior to blastulation in mammals and is necessary for embryogenesis. Sox expression is associated with pluripotency and stem cells, neuronal differentiation, gut development, and cancer. Schistosomes express a Sox-like gene expressed in the schistosomula after infecting a mammalian host when schistosomes have about 900 cells. Here, we characterized and named this Sox-like gene SmSOXS1. SmSoxS1 protein is a developmentally regulated activator that localizes to the anterior and posterior ends of the schistosomula and binds to Sox-specific DNA elements. In addition to SmSoxS1, we have also identified an additional six Sox genes in schistosomes, two Sox B, one SoxC, and three Sox genes that may establish a flatworm-specific class of Sox genes with planarians. These data identify novel Sox genes in schistosomes to expand the potential functional roles for Sox2 and may provide interesting insights into early multicellular development of flatworms.

Developmentally regulated internal transcription initiation during meiosis in budding yeast

Sporulation of budding yeast is a developmental process in which cells undergo meiosis to generate stress-resistant progeny. The dynamic nature of the budding yeast meiotic transcriptome has been well established by a number of genome-wide studies. Here we develop an analysis pipeline to systematically identify novel transcription start sites that reside internal to a gene. Application of this pipeline to data from a synchronized meiotic time course reveals over 40 genes that display specific internal initiations in mid-sporulation. Consistent with the time of induction, motif analysis on upstream sequences of these internal transcription start sites reveals a significant enrichment for the binding site of Ndt80, the transcriptional activator of middle sporulation genes. Further examination of one gene, MRK1, demonstrates the Ndt80 binding site is necessary for internal initiation and results in the expression of an N-terminally truncated protein isoform. When the MRK1 paralog RIM11 is downregulated, the MRK1 internal transcript promotes efficient sporulation, indicating functional significance of the internal initiation. Our findings suggest internal transcriptional initiation to be a dynamic, regulated process with potential functional impacts on development.

Gene regulatory network plasticity predates a switch in function of a conserved transcription regulator

The rewiring of gene regulatory networks can generate phenotypic novelty. It remains an open question, however, how the large number of connections needed to form a novel network arise over evolutionary time. Here, we address this question using the network controlled by the fungal transcription regulator Ndt80. This conserved protein has undergone a dramatic switch in function-from an ancestral role regulating sporulation to a derived role regulating biofilm formation. This switch in function corresponded to a large-scale rewiring of the genes regulated by Ndt80. However, we demonstrate that the Ndt80-target gene connections were undergoing extensive rewiring prior to the switch in Ndt80's regulatory function. We propose that extensive drift in the Ndt80 regulon allowed for the exploration of alternative network structures without a loss of ancestral function, thereby facilitating the formation of a network with a new function.

Functional Profiling of Human MeCP2 by Automated Data Comparison Analysis and Computerized Expression Pathway Modeling

OBJECTIVES

Methyl-CpG binding protein 2 (MeCP2) is a ubiquitous epigenetic factor that represses gene expression by modifying chromatin. Mutations in the MeCP2 gene cause Rett syndrome, a progressive neurodevelopmental disorder. Recent studies also have shown that MeCP2 plays a role in carcinogenesis. Specifically, functional ablation of MeCP2 suppresses cell growth and leads to the proliferation of cancer cells. However, MeCP2's function in adult tissues remains poorly understood. We utilized a weight matrix-based comparison software to identify transcription factor binding site (TFBS) of MeCP2-regulated genes, which were recognized by cDNA microarray analysis.

METHODS

MeCP2 expression was silenced using annealed siRNA in HEK293 cells, and then a cDNA microarray analysis was performed. Functional analysis was carried out, and transcriptional levels in target genes regulated by MeCP2 were investigated. TFBS analysis was done within genes selected by the cDNA microarray analysis, using a weight matrix-based program and the TRANSFAC 6.0 database.

RESULTS

Among the differentially expressed genes with a change in expression greater than two-fold, 189 genes were up-regulated and 91 genes were down-regulated. Genes related to apoptosis and cell proliferation (JUN, FOSL2, CYR61, SKIL, ATF3, BMABI, BMPR2, RERE, and FALZ) were highly up-regulated. Genes with anti-apoptotic and anti-proliferative functions (HNRPA0, HIS1, and FOXC1) were down-regulated. Using TFBS analysis within putative promoters of novel candidate target genes of MeCP2, disease-related transcription factors were identified.

CONCLUSIONS

The present results provide insights into the new target genes regulated by MeCP2 under epigenetic control. This information will be valuable for further studies aimed at clarifying the pathogenesis of Rett syndrome and neoplastic diseases.

Identification and characterization of a Mef2 transcriptional activator in schistosome parasites

Myocyte enhancer factor 2 protein (Mef2) is an evolutionarily conserved activator of transcription that is critical to induce and control complex processes in myogenesis and neurogenesis in vertebrates and insects, and osteogenesis in vertebrates. In Drosophila, Mef2 null mutants are unable to produce differentiated muscle cells, and in vertebrates, Mef2 mutants are embryonic lethal. Schistosome worms are responsible for over 200 million cases of schistosomiasis globally, but little is known about early development of schistosome parasites after infecting a vertebrate host. Understanding basic schistosome development could be crucial to delineating potential drug targets. Here, we identify and characterize Mef2 from the schistosome worm Schistosoma mansoni (SmMef2). We initially identified SmMef2 as a homolog to the yeast Mef2 homolog, Resistance to Lethality of MKK1P386 overexpression (Rlm1), and we show that SmMef2 is homologous to conserved Mef2 family proteins. Using a genetics approach, we demonstrate that SmMef2 is a transactivator that can induce transcription of four separate heterologous reporter genes by yeast one-hybrid analysis. We also show that Mef2 is expressed during several stages of schistosome development by quantitative PCR and that it can bind to conserved Mef2 DNA consensus binding sequences.

Schistosome parasites infect more than 200 million people worldwide and cause human schistosomiasis. Free-swimming schistosome larvae are highly mobile and invade and penetrate the host's skin to perpetuate their lifecycle in their human host, growing from 90–215 micrometers in length as a schistosomulum to a 7–20 millimeter long adult worm. Few molecular pathways have been identified in schistosome worms that are important for parasite early development. The myocyte enhancer factor protein 2 is a major regulator of muscle and nerve development in mammals and insects and is highly conserved from bread yeast to vertebrates. Here we identify and characterize the Mef2 activator from parasitic schistosome worms, the first described in any parasitic worm, and delineation of its function may be important to further understanding the basic biology of schistosome early development. Additionally, since schistosomes developed early evolutionarily, an investigation of schistosome Mef2 regulatory mechanisms could lead to a greater understanding of the development of early muscle and neurogenic development in animals.

An orthologue of the Myelin-gene Regulatory Transcription Factor regulates Dictyostelium prestalk differentiation

The prestalk region of the Dictyostelium slug is comprised of an anterior population of pstA cells and a posterior population of pstO cells. They are distinguished by their ability to utilize different parts of the promoter of the ecmA gene. We identify, by mutational analysis and DNA transformation, CA-rich sequence elements within the ecmA promoter that are essential for pstA-specific expression and sufficient to direct pstA-specific expression when multimerised. The CA-rich region was used in affinity chromatography with nuclear extracts and bound proteins were identified by mass spectrometry. The CA-rich elements purify MrfA, a protein with extensive sequence similarity to animal Myelin-gene Regulatory Factor (MRF)-like proteins. The MRF-like proteins and MrfA also display more spatially limited but significant sequence similarity with the DNA binding domain of the yeast Ndt80 sporulation-specific transcription factor. Furthermore, the ecmA CA-rich elements show sequence similarity to the core consensus Ndt80 binding site (the MSE) and point mutation of highly conserved arginine residues in MrfA, that in Ndt80 make critical contacts with the MSE, ablate binding of MrfA to its sites within the ecmA promoter. MrfA null strains are delayed in multicellular development and highly defective in pstA-specific gene expression. These results provide a first insight into the intracellular signaling pathway that directs pstA differentiation and identify a non-metazoan orthologue of a family of molecularly uncharacterised transcription factors.

Cdc7-Dbf4 is a gene-specific regulator of meiotic transcription in yeast

Meiosis divides the chromosome number of the cell in half by having two rounds of chromosome segregation follow a single round of chromosome duplication. The first meiotic division is unique in that homologous pairs of sister chromatids segregate to opposite poles. Recent work in budding and fission yeast has shown that the cell cycle kinase, Cdc7-Dbf4, is required for many meiosis-specific chromosomal functions necessary for proper disjunction at meiosis I. This work reveals another role for Cdc7 in meiosis as a gene-specific regulator of the global transcription factor, Ndt80, which is required for exit from pachytene and entry into the meiotic divisions in budding yeast. Cdc7-Dbf4 promotes NDT80 transcription by relieving repression mediated by a complex of Sum1, Rfm1, and a histone deacetylase, Hst1. Sum1 exhibits meiosis-specific Cdc7-dependent phosphorylation, and mass spectrometry analysis reveals a dynamic and complex pattern of phosphorylation events, including four constitutive cyclin-dependent kinase (Cdk1) sites and 11 meiosis-specific Cdc7-Dbf4-dependent sites. Analysis of various phosphorylation site mutants suggests that Cdc7 functions with both Cdk1 and the meiosis-specific kinase Ime2 to control this critical transition point during meiosis.

The functional consequences of relative substrate specificity in complex biochemical systems

A biochemical activity, that is, enzymatic reaction or molecular interaction, frequently involves a molecule, for example, an enzyme, capable of interacting with numerous substrates or partners. Specificity is a fundamental property of biochemical activities, and relative specificity refers to the situation whereby a molecule interacts with multiple substrates or partners but with different affinities. Here, a hypothesis is proposed that any molecule, such as an enzyme, would have a range of preferences or relative specificity for its many native substrates, which differentially impacts the phenotypes of these substrates and hence shapes the relevant biological processes in vivo. While the mechanisms underlying the specific recognition between enzymes and individual substrates have been studied extensively, whether any enzyme exhibits intrinsic selectivity toward its ensemble of substrates is often overlooked, and whether this selectivity has any functional consequences is much less appreciated. There are, however, several lines of evidence in the literature that are consistent with the hypothesis and reviewed here. Furthermore, this hypothesis is supported by our analyses of a number of diverse biochemical systems at a large scale. Thus, the human microRNA processing machinery possesses relative specificity toward its hundreds of substrates, which might contribute to differential microRNA biogenesis; the promoter binding affinity of the transcription factor Ndt80 might regulate Ndt80 target mRNA expression in the budding yeast; Cdk1 kinase specificity might lead to variable substrate phosphorylation in vivo; and the density of HuR deposition to its thousands of RNA targets might partly explain differential RNA expression in human cells. It is proposed, therefore, that relative specificity is a universal property of complex biochemical systems and that the hypothesis could denote a general principle in biology.

Identification of functional TFAP2A and SP1 binding sites in new TFAP2A-modulated genes

Background

Different approaches have been developed to dissect the interplay between transcription factors (TFs) and their cis-acting sequences on DNA in order to identify TF target genes. Here we used a combination of computational and experimental approaches to identify novel direct targets of TFAP2A, a key TF for a variety of physiological and pathological cellular processes. Gene expression profiles of HeLa cells either silenced for TFAP2A by RNA interference or not were previously compared and a set of differentially expressed genes was revealed.

Results

The regulatory regions of 494 TFAP2A-modulated genes were analyzed for the presence of TFAP2A binding sites, employing the canonical TFAP2A Positional Weight Matrix (PWM) reported in Jaspar http://jaspar.genereg.net/. 264 genes containing at least 2 high score TFAP2A binding sites were identified, showing a central role in "Cellular Movement" and "Cellular Development". In an attempt to identify TFs that could cooperate with TFAP2A, a statistically significant enrichment for SP1 binding sites was found for TFAP2A-activated but not repressed genes. The direct binding of TFAP2A or SP1 to a random subset of TFAP2A-modulated genes was demonstrated by Chromatin ImmunoPrecipitation (ChIP) assay and the TFAP2A-driven regulation of DCBLD2/ESDN/CLCP1 gene studied in details.

Conclusions

We proved that our computational approaches applied to microarray selected genes are valid tools to identify functional TF binding sites in gene regulatory regions as confirmed by experimental validations. In addition, we demonstrated a fine-tuned regulation of DCBLD2/ESDN transcription by TFAP2A.

Mitotic expression of Spo13 alters M-phase progression and nucleolar localization of Cdc14 in budding yeast

Spo13 is a key meiosis-specific regulator required for centromere cohesion and coorientation, and for progression through two nuclear divisions. We previously reported that it causes a G2/M arrest and may delay the transition from late anaphase to G1, when overexpressed in mitosis. Yet its mechanism of action has remained elusive. Here we show that Spo13, which is phosphorylated and stabilized at G2/M in a Cdk/Clb-dependent manner, acts at two stages during mitotic cell division. Spo13 provokes a G2/M arrest that is reversible and largely independent of the Mad2 spindle checkpoint. Since mRNAs whose induction requires Cdc14 activation are reduced, we propose that its anaphase delay results from inhibition of Cdc14 function. Indeed, the Spo13-induced anaphase delay correlates with Cdc14 phosphatase retention in the nucleolus and with cyclin B accumulation, which both impede anaphase exit. At the onset of arrest, Spo13 is primarily associated with the nucleolus, where Cdc14 accumulates. Significantly, overexpression of separase (Esp1), which promotes G2/M and anaphase progression, suppresses Spo13 effects in mitosis, arguing that Spo13 acts upstream or parallel to Esp1. Given that Spo13 overexpression reduces Pds1 and cyclin B degradation, our findings are consistent with a role for Spo13 in regulating APC, which controls both G2/M and anaphase. Similar effects of Spo13 during meiotic MI may prevent cell cycle exit and initiation of DNA replication prior to MII, thereby ensuring two successive chromosome segregation events without an intervening S phase.

Analysis of plant regulatory DNA sequences by the yeast-one-hybrid assay

Regulatory DNA sequences harbor the essential information to control specific gene expression changes and integrate information derived from upstream signaling cascades. This regulatory potential is mediated by direct binding of proteins, e.g., transcription factors, to defined stretches of DNA motifs in regulatory regions. The analysis of these DNA regions, at which several signaling pathways could merge to orchestrate gene expression, is still a challenging task. To date, the combination of functional approaches in the laboratory and computer aided sequence evaluation is frequently used for regulatory sequence analysis. The yeast-one-hybrid method is a possible approach to test for direct binding of plant proteins to DNA in a heterologous system. Moreover, it is the most frequently used method for the identification of DNA-binding proteins targeting a given DNA sequence by screening a cDNA library.

An in silico strategy identified the target gene candidates regulated by dehydration responsive element binding proteins (DREBs) in Arabidopsis genome

Identification of downstream target genes of stress-relating transcription factors (TFs) is desirable in understanding cellular responses to various environmental stimuli. However, this has long been a difficult work for both experimental and computational practices. In this research, we presented a novel computational strategy which combined the analysis of the transcription factor binding site (TFBS) contexts and machine learning approach. Using this strategy, we conducted a genome-wide investigation into novel direct target genes of dehydration responsive element binding proteins (DREBs), the members of AP2-EREBPs transcription factor super family which is reported to be responsive to various abiotic stresses in Arabidopsis. The genome-wide searching yielded in total 474 target gene candidates. With reference to the microarray data for abiotic stresses-inducible gene expression profile, 268 target gene candidates out of the total 474 genes predicted, were induced during the 24-h exposure to abiotic stresses. This takes about 57% of total predicted targets. Furthermore, GO annotations revealed that these target genes are likely involved in protein amino acid phosphorylation, protein binding and Endomembrane sorting system. The results suggested that the predicted target gene candidates were adequate to meet the essential biological principle of stress-resistance in plants.

The Saccharomyces cerevisiae CLB5 promoter contains two middle sporulation elements (MSEs) that are differentially regulated during sporulation

The B-type cyclins Clb5 and Clb6 are essential activators of DNA replication during sporulation in Saccharomyces cerevisiae. The expression of CLB5 is maximally induced during the middle phase of sporulation by the transcription factor Ndt80. We have performed an analysis of the CLB5 promoter and have identified two middle sporulation elements (MSEs) that act as binding sites for Ndt80. Although both MSE sequences bind Ndt80 in vitro, they display differential effectiveness in their ability to function as cis-acting regulatory sequences in vivo. Mutation of both MSE sequences in the CLB5 promoter profoundly reduces the induction of CLB5 transcription during the middle phase of sporulation but results in no obvious defect in progression through meiosis and sporulation, implying that the Ndt80-dependent induction of CLB5 is not required for effective DNA replication or chromosome division.

Computational methods to dissect cis-regulatory transcriptional networks

The formation of diverse cell types from an invariant set of genes is governed by biochemical and molecular processes that regulate gene activity. A complete understanding of the regulatory mechanisms of gene expression is the major function of genomics. Computational genomics is a rapidly emerging area for deciphering the regulation of metazoan genes as well as interpreting the results of high-throughput screening. The integration of computer science with biology has expedited molecular modelling and processing of large-scale data inputs such as microarrays, analysis of genomes, transcriptomes and proteomes. Many bioinformaticians have developed various algorithms for predicting transcriptional regulatory mechanisms from the sequence, gene expression and interaction data. This review contains compiled information of various computational methods adopted to dissect gene expression pathways.

Locating mammalian transcription factor binding sites: a survey of computational and experimental techniques

Fields such as genomics and systems biology are built on the synergism between computational and experimental techniques. This type of synergism is especially important in accomplishing goals like identifying all functional transcription factor binding sites in vertebrate genomes. Precise detection of these elements is a prerequisite to deciphering the complex regulatory networks that direct tissue specific and lineage specific patterns of gene expression. This review summarizes approaches for in silico, in vitro, and in vivo identification of transcription factor binding sites. A variety of techniques useful for localized- and high-throughput analyses are discussed here, with emphasis on aspects of data generation and verification.

Gene expression from random libraries of yeast promoters

Genomewide techniques to assay gene expression and transcription factor binding are in widespread use, but are far from providing predictive rules for the function of regulatory DNA. To investigate more intensively the grammar rules for active regulatory sequence, we made libraries from random ligations of a very restricted set of sequences. Working with the yeast Saccharomyces cerevisiae, we developed a novel screen based on the sensitivity of ascospores lacking dityrosine to treatment with lytic enzymes. We tested two separate libraries built by random ligation of a single type of activator site either for a well-characterized sporulation factor, Ndt80, or for a new sporulation-specific regulatory site that we identified and several neutral spacer elements. This selective system achieved up to 1:10(4) enrichment of the artificial sequences that were active during sporulation, allowing a high-throughput analysis of large libraries of synthetic promoters. This is not practical with methods involving direct screening for expression, such as those based on fluorescent reporters. There were very few false positives, since active promoters always passed the screen when retested. The survival rate of our libraries containing roughly equal numbers of spacers and activators was a few percent that of libraries made from activators alone. The sequences of approximately 100 examples of active and inactive promoters could not be distinguished by simple binary rules; instead, the best model for the data was a linear regression fit of a quantitative measure of gene activity to multiple features of the regulatory sequence.