GibbsST: a Gibbs sampling method for motif discovery with enhanced resistance to local optima

BACKGROUND: Computational discovery of transcription factor binding sites (TFBS) is a challenging but important problem of bioinformatics. In this study, improvement of a Gibbs sampling based technique for TFBS discovery is attempted through an approach that is widely known, but which has never been investigated before: reduction of the effect of local optima. RESULTS: To alleviate the vulnerability of Gibbs sampling to local optima trapping, we propose to combine a thermodynamic method, called simulated tempering, with Gibbs sampling. The resultant algorithm, GibbsST, is then validated using synthetic data and actual promoter sequences extracted from Saccharomyces cerevisiae. It is noteworthy that the marked improvement of the efficiency presented in this paper is attributable solely to the improvement of the search method. CONCLUSION: Simulated tempering is a powerful solution for local optima problems found in pattern discovery. Extended application of simulated tempering for various bioinformatic problems is promising as a robust solution against local optima problems

The Stress Response Factors Yap6, Cin5, Phd1, and Skn7 Direct Targeting of the Conserved Co-Repressor Tup1-Ssn6 in S. cerevisiae

Maintaining the proper expression of the transcriptome during development or in response to a changing environment requires a delicate balance between transcriptional regulators with activating and repressing functions. The budding yeast transcriptional co-repressor Tup1-Ssn6 is a model for studying similar repressor complexes in multicellular eukaryotes. Tup1-Ssn6 does not bind DNA directly, but is directed to individual promoters by one or more DNA-binding proteins, referred to as Tup1 recruiters. This functional architecture allows the Tup1-Ssn6 to modulate the expression of genes required for the response to a variety of cellular stresses. To understand the targeting or the Tup1-Ssn6 complex, we determined the genomic distribution of Tup1 and Ssn6 by ChIP-chip. We found that most loci bound by Tup1-Ssn6 could not be explained by co-occupancy with a known recruiting cofactor and that deletion of individual known Tup1 recruiters did not significantly alter the Tup1 binding profile. These observations suggest that new Tup1 recruiting proteins remain to be discovered and that Tup1 recruitment typically depends on multiple recruiting cofactors. To identify new recruiting proteins, we computationally screened for factors with binding patterns similar to the observed Tup1-Ssn6 genomic distribution. Four top candidates, Cin5, Skn7, Phd1, and Yap6, all known to be associated with stress response gene regulation, were experimentally confirmed to physically interact with Tup1 and/or Ssn6. Incorporating these new recruitment cofactors with previously characterized cofactors now explains the majority of Tup1 targeting across the genome, and expands our understanding of the mechanism by which Tup1-Ssn6 is directed to its targets

Identification of the CRE-1 Cellulolytic Regulon in Neurospora crassa

Background: In filamentous ascomycete fungi, the utilization of alternate carbon sources is influenced by the zinc finger transcription factor CreA/CRE-1, which encodes a carbon catabolite repressor protein homologous to Mig1 from Saccharomyces cerevisiae. In Neurospora crassa, deletion of cre-1 results in increased secretion of amylase and b-galactosidase. Methodology/Principal Findings: Here we show that a strain carrying a deletion of cre-1 has increased cellulolytic activity and increased expression of cellulolytic genes during growth on crystalline cellulose (Avicel). Constitutive expression of cre-1 complements the phenotype of a N. crassa Dcre-1 strain grown on Avicel, and also results in stronger repression of cellulolytic protein secretion and enzyme activity. We determined the CRE-1 regulon by investigating the secretome and transcriptome of a Dcre-1 strain as compared to wild type when grown on Avicel versus minimal medium. Chromatin immunoprecipitation-PCR of putative target genes showed that CRE-1 binds to only some adjacent 59-SYGGRG-39 motifs, consistent with previous findings in other fungi, and suggests that unidentified additional regulatory factors affect CRE-1 binding to promoter regions. Characterization of 30 mutants containing deletions in genes whose expression level increased in a Dcre-1 strain under cellulolytic conditions identified novel genes that affect cellulase activity and protein secretion

Interplay of Nkx3.2, Sox9 and Pax3 Regulates Chondrogenic Differentiation of Muscle Progenitor Cells

Muscle satellite cells make up a stem cell population that is capable of differentiating into myocytes and contributing to muscle regeneration upon injury. In this work we investigate the mechanism by which these muscle progenitor cells adopt an alternative cell fate, the cartilage fate. We show that chick muscle satellite cells that normally would undergo myogenesis can be converted to express cartilage matrix proteins in vitro when cultured in chondrogenic medium containing TGFß3 or BMP2. In the meantime, the myogenic program is repressed, suggesting that muscle satellite cells have undergone chondrogenic differentiation. Furthermore, ectopic expression of the myogenic factor Pax3 prevents chondrogenesis in these cells, while chondrogenic factors Nkx3.2 and Sox9 act downstream of TGFß or BMP2 to promote this cell fate transition. We found that Nkx3.2 and Sox9 repress the activity of the Pax3 promoter and that Nkx3.2 acts as a transcriptional repressor in this process. Importantly, a reverse function mutant of Nkx3.2 blocks the ability of Sox9 to both inhibit myogenesis and induce chondrogenesis, suggesting that Nkx3.2 is required for Sox9 to promote chondrogenic differentiation in satellite cells. Finally, we found that in an in vivo mouse model of fracture healing where muscle progenitor cells were lineage-traced, Nkx3.2 and Sox9 are significantly upregulated while Pax3 is significantly downregulated in the muscle progenitor cells that give rise to chondrocytes during fracture repair. Thus our in vitro and in vivo analyses suggest that the balance of Pax3, Nkx3.2 and Sox9 may act as a molecular switch during the chondrogenic differentiation of muscle progenitor cells, which may be important for fracture healing