Regulation of mating and filamentation genes by two distinct Ste12 complexes in Saccharomyces cerevisiae

The transcription factor STE12 influences growth on several carbon sources and production of dehydroacetic acid (DHAA) in Trichoderma reesei

The filamentous ascomycete Trichoderma reesei, known for its prolific cellulolytic enzyme production, recently also gained attention for its secondary metabolite synthesis. Both processes are intricately influenced by environmental factors like carbon source availability and light exposure. Here, we explore the role of the transcription factor STE12 in regulating metabolic pathways in T. reesei in terms of gene regulation, carbon source utilization and biosynthesis of secondary metabolites. We show that STE12 is involved in regulating cellulase gene expression and growth on carbon sources associated with iron homeostasis. STE12 impacts gene regulation in a light dependent manner on cellulose with modulation of several CAZyme encoding genes as well as genes involved in secondary metabolism. STE12 selectively influences the biosynthesis of the sorbicillinoid trichodimerol, while not affecting the biosynthesis of bisorbibutenolide, which was recently shown to be regulated by the MAPkinase pathway upstream of STE12 in the signaling cascade. We further report on the biosynthesis of dehydroacetic acid (DHAA) in T. reesei, a compound known for its antimicrobial properties, which is subject to regulation by STE12. We conclude, that STE12 exerts functions beyond development and hence contributes to balance the energy distribution between substrate consumption, reproduction and defense.

A network of transcription factors in complex with a regulating cell cycle cyclin orchestrates fungal oxidative stress responses

BACKGROUND

Response to oxidative stress is universal in almost all organisms and the mitochondrial membrane protein, BbOhmm, negatively affects oxidative stress responses and virulence in the insect fungal pathogen, Beauveria bassiana. Nothing further, however, is known concerning how BbOhmm and this phenomenon is regulated.

RESULTS

Three oxidative stress response regulating Zn2Cys6 transcription factors (BbOsrR1, 2, and 3) were identified and verified via chromatin immunoprecipitation (ChIP)-qPCR analysis as binding to the BbOhmm promoter region, with BbOsrR2 showing the strongest binding. Targeted gene knockout of BbOsrR1 or BbOsrR3 led to decreased BbOhmm expression and consequently increased tolerances to free radical generating compounds (H2O2 and menadione), whereas the ΔBbOsrR2 strain showed increased BbOhmm expression with concomitant decreased tolerances to these compounds. RNA and ChIP sequencing analysis revealed that BbOsrR1 directly regulated a wide range of antioxidation and transcription-associated genes, negatively affecting the expression of the BbClp1 cyclin and BbOsrR2. BbClp1 was shown to localize to the cell nucleus and negatively mediate oxidative stress responses. BbOsrR2 and BbOsrR3 were shown to feed into the Fus3-MAPK pathway in addition to regulating antioxidation and detoxification genes. Binding motifs for the three transcription factors were found to partially overlap in the promoter region of BbOhmm and other target genes. Whereas BbOsrR1 appeared to function independently, co-immunoprecipitation revealed complex formation between BbClp1, BbOsrR2, and BbOsrR3, with BbClp1 partially regulating BbOsrR2 phosphorylation.

CONCLUSIONS

These findings reveal a regulatory network mediated by BbOsrR1 and the formation of a BbClp1-BbOsrR2-BbOsrR3 complex that orchestrates fungal oxidative stress responses.

The serine/threonine protein kinase MpSTE1 directly governs hyphal branching in Monascus spp

Monascus spp. are commercially important fungi due to their ability to produce beneficial secondary metabolites such as the cholesterol-lowering agent lovastatin and natural food colorants azaphilone pigments. Although hyphal branching intensively influenced the production of these secondary metabolites, the pivotal regulators of hyphal development in Monascus spp. remain unclear. To identify these important regulators, we developed an artificial intelligence (AI)-assisted image analysis tool for quantification of hyphae-branching and constructed a random T-DNA insertion library. High-throughput screening revealed that a STE kinase, MpSTE1, was considered as a key regulator of hyphal branching based on the hyphal phenotype. To further validate the role of MpSTE1, we generated an mpSTE1 gene knockout mutant, a complemented mutant, and an overexpression mutant (OE::mpSTE1). Microscopic observations revealed that overexpression of mpSTE1 led to a 63% increase in branch number while deletion of mpSTE1 reduced the hyphal branching by 68% compared to the wild-type strain. In flask cultures, the strain OE::mpSTE1 showed accelerated growth and glucose consumption. More importantly, the strain OE::mpSTE1 produced 9.2 mg/L lovastatin and 17.0 mg/L azaphilone pigments, respectively, 47.0% and 30.1% higher than those of the wild-type strain. Phosphoproteomic analysis revealed that MpSTE1 directly phosphorylated 7 downstream signal proteins involved in cell division, cytoskeletal organization, and signal transduction. To our best knowledge, MpSTE1 is reported as the first characterized regulator for tightly regulating the hyphal branching in Monascus spp. These findings significantly expanded current understanding of the signaling pathway governing the hyphal branching and development in Monascus spp. Furthermore, MpSTE1 and its analogs were demonstrated as promising targets for improving production of valuable secondary metabolites. KEY POINTS: • MpSTE1 is the first characterized regulator for tightly regulating hyphal branching • Overexpression of mpSTE1 significantly improves secondary metabolite production • A high-throughput image analysis tool was developed for counting hyphal branching.

Mutations in Podospora anserina MCM1 and VelC Trigger Spontaneous Development of Barren Fruiting Bodies

The ascomycete Podospora anserina is a heterothallic filamentous fungus found mainly on herbivore dung. It is commonly used in laboratories as a model system, and its complete life cycle lasting eight days is well mastered in vitro. The main objective of our team is to understand better the global process of fruiting body development, named perithecia, induced normally in this species by fertilization. Three allelic mutants, named pfd3, pfd9, and pfd23 (for "promoting fruiting body development") obtained by UV mutagenesis, were selected in view of their abilities to promote barren perithecium development without fertilization. By complete genome sequencing of pfd3 and pfd9, and mutant complementation, we identified point mutations in the mcm1 gene as responsible for spontaneous perithecium development. MCM1 proteins are MADS box transcription factors that control diverse developmental processes in plants, metazoans, and fungi. We also identified using the same methods a mutation in the VelC gene as responsible for spontaneous perithecium development in the vacua mutant. The VelC protein belongs to the velvet family of regulators involved in the control of development and secondary metabolite production. A key role of MCM1 and VelC in coordinating the development of P. anserina perithecia with gamete formation and fertilization is highlighted.

The mitogen-activated protein kinase module CcSte11-CcSte7-CcPmk1 regulates pathogenicity via the transcription factor CcSte12 in Cytospora chrysosperma

The pathogen Cytospora chrysosperma is the causal agent of poplar canker disease and causes considerable economic losses in China. Mitogen-activated protein kinase (MAPK) cascades play a crucial role in mediating cellular responses and Pmk1-MAPKs are indispensable for pathogenic related processes in plant pathogenic fungi. In previous studies, we demonstrated that the CcPmk1 acts as a core regulator of fungal pathogenicity by modulating a small number of master downstream targets, such as CcSte12. In this study, we identified and characterized two upstream components of CcPmk1: MAPKKK CcSte11 and MAPKK CcSte7. Deletion of CcSte11 and CcSte7, resulted in slowed growth, loss of sporulation and virulence, similar to the defects observed in the CcPmk1 deletion mutant. In addition, CcSte11, CcSte7 and CcPmk1 interact with each other, and the upstream adaptor protein CcSte50 interact with CcSte11 and CcSte7. Moreover, we explored the global regulation network of CcSte12 by transcriptional analysis between CcSte12 deletion mutants and wild-type during the simulated infection process. Two hydrolase activity GO terms (GO:0004553 and GO:0016798) and starch and sucrose metabolism (mgr00500) KEGG pathway were significantly enriched in the down-regulated genes of CcSte12 deletion mutants. In addition, a subset of glycosyl hydrolase genes and putative effector genes were significantly down-regulated in the CcSte12 deletion mutant, which might be important for fungal pathogenicity. Especially, CcSte12 bound to the CcSp84 promoter region containing the TGAAACA motif. Moreover, comparison of CcSte12-regulated genes with CcPmk1-regulated genes revealed 116 overlapping regulated genes in both CcSte12 and CcPmk1, including some virulence-associated genes. Taken together, the protein complexes CcSte11-CcSte7-CcPmk1 receive signals transmitted by upstream CcSte50 and transmit signals to downstream CcSte12, which regulates hydrolase, effectors and other genes to promote virulence. Overall, these results indicate that the CcPmk1-MAPK signaling pathway of C. chrysosperma plays a key role in the pathogenicity.

Transcription Factors Tec1 and Tec2 Play Key Roles in the Hyphal Growth and Virulence of Mucor lusitanicus Through Increased Mitochondrial Oxidative Metabolism

Mucormycosis is a lethal and difficult-to-treat fungal infection caused by fungi of the order Mucorales. Mucor lusitanicus, a member of Mucorales, is commonly used as a model to understand disease pathogenesis. However, transcriptional control of hyphal growth and virulence in Mucorales is poorly understood. This study aimed to investigate the role of Tec proteins, which belong to the TEA/ATTS transcription factor family, in the hyphal development and virulence of M. lusitanicus. Unlike in the genome of Ascomycetes and Basidiomycetes, which have a single Tec homologue, in the genome of Mucorales, two Tec homologues, Tec1 and Tec2, were found, except in that of Phycomyces blakesleeanus, with only one Tec homologue. tec1 and tec2 overexpression in M. lusitanicus increased mycelial growth, mitochondrial content and activity, expression of the rhizoferrin synthetase-encoding gene rfs, and virulence in nematodes and wax moth larvae but decreased cAMP levels and protein kinase A (PKA) activity. Furthermore, tec1- and tec2-overexpressing strains required adequate mitochondrial metabolism to promote the virulent phenotype. The heterotrimeric G beta subunit 1-encoding gene deletant strain (Δgpb1) increased cAMP-PKA activity, downregulation of both tec genes, decreased both virulence and hyphal development, but tec1 and tec2 overexpression restored these defects. Overexpression of allele-mutated variants of Tec1(S332A) and Tec2(S168A) in the putative phosphorylation sites for PKA increased both virulence and hyphal growth of Δgpb1. These findings suggest that Tec homologues promote mycelial development and virulence by enhancing mitochondrial metabolism and rhizoferrin accumulation, providing new information for the rational control of the virulent phenotype of M. lusitanicus.

Ecological inducers of the yeast filamentous growth pathway reveal environment-dependent roles for pathway components

Signaling modules, such as mitogen-activated protein kinase (MAPK) pathways, are evolutionarily conserved drivers of cell differentiation and stress responses. In many fungal species including pathogens, MAPK pathways control filamentous growth, where cells differentiate into an elongated cell type. The convenient model budding yeast Saccharomyces cerevisiae undergoes filamentous growth by the filamentous growth (fMAPK) pathway; however, the inducers of the pathway remain unclear, perhaps because pathway activity has been mainly studied in laboratory conditions. To address this knowledge gap, an ecological framework was used, which uncovered new fMAPK pathway inducers, including pectin, a material found in plants, and the metabolic byproduct ethanol. We also show that induction by a known inducer of the pathway, the non-preferred carbon source galactose, required galactose metabolism and induced the pathway differently than glucose limitation or other non-preferred carbon sources. By exploring fMAPK pathway function in fruit, we found that induction of the pathway led to visible digestion of fruit rind through a known target, PGU1, which encodes a pectolytic enzyme. Combinations of inducers (galactose and ethanol) stimulated the pathway to near-maximal levels, which showed dispensability of several fMAPK pathway components (e.g., mucin sensor, p21-activated kinase), but not others (e.g., adaptor, MAPKKK) and required the Ras2-protein kinase A pathway. This included a difference between the transcription factor binding partners for the pathway, as Tec1p, but not Ste12p, was partly dispensable for fMAPK pathway activity. Thus, by exploring ecologically relevant stimuli, new modes of MAPK pathway signaling were uncovered, perhaps revealing how a pathway can respond differently to specific environments. IMPORTANCE Filamentous growth is a cell differentiation response and important aspect of fungal biology. In plant and animal fungal pathogens, filamentous growth contributes to virulence. One signaling pathway that regulates filamentous growth is an evolutionarily conserved MAPK pathway. The yeast Saccharomyces cerevisiae is a convenient model to study MAPK-dependent regulation of filamentous growth, although the inducers of the pathway are not clear. Here, we exposed yeast cells to ecologically relevant compounds (e.g., plant compounds), which identified new inducers of the MAPK pathway. In combination, the inducers activated the pathway to near-maximal levels but did not cause detrimental phenotypes associated with previously identified hyperactive alleles. This context allowed us to identify conditional bypass for multiple pathway components. Thus, near-maximal induction of a MAPK pathway by ecologically relevant inducers provides a powerful tool to assess cellular signaling during a fungal differentiation response.

Fus3/Kss1-MAP kinase and Ste12-like control distinct biocontrol-traits besides regulation of insect cuticle penetration via phosphorylation cascade in a filamentous fungal pathogen

BACKGROUND

Homolog of the yeast Fus3/Kss1 mitogen-activated protein kinase (MAPK) pathway and its target transcription factor, Ste12-like, are involved in penetration of host cuticle/pathogenicity in many ascomycete pathogens. However, details of their interaction during fungal infection, as well as their controlled other virulence-associated traits, are unclear.

RESULTS

Ste12-like (BbSte12) and Fus3/Kss1 MAPK homolog (Bbmpk1) interacted in nucleus, and phosphorylation of BbSte12 by Bbmpk1 was essential for penetration of insect cuticle in an insect fungal pathogen, Beauveria bassiana. However, some distinct biocontrol-traits were found to be mediated by Ste12 and Bbmpk1. In contrast to ΔBbmpk1 colony that grew more rapid than wild-type strain, inactivation of BbSte12 resulted in the opposite phenotype, which was consistent with their different proliferation rates in insect hemocoel after direct injection of conidia bypass the cuticle. Reduced conidial yield with decreased hydrophobicity was examined in both mutants, however they displayed distinct conidiogenesis, accompanying with differently altered cell cycle, distinct hyphal branching and septum formation. Moreover, ΔBbmpk1 showed increased tolerance to oxidative agent, whereas the opposite phenotype was seen for ΔBbSte12 strain. RNA sequencing analysis revealed that Bbmpk1 controlled 356 genes depending on BbSte12 during cuticle penetration, but 1077 and 584 genes were independently controlled by Bbmpk1 and BbSte12.

CONCLUSION

BbSte12 and Bbmpk1 separately participate in additional pathways for control of conidiation, growth and hyphal differentiation, as well as oxidative stress response besides regulating cuticle penetration via phosphorylation cascade. © 2023 Society of Chemical Industry.

Management of mung bean leaf spot disease caused by Phoma herbarum through Penicillium janczewskii metabolites mediated by MAPK signaling cascade

Vigna radiata L., an imperative legume crop of Pakistan, faces hordes of damage due to fungi; infecting host tissues by the appressorium. The use of natural compounds is an innovative concern to manage mung-bean fungal diseases. The bioactive secondary metabolites of Penicillium species are well documented for their strong fungi-static ability against many pathogens. Presently, one-month-old aqueous culture filtrates of Penicillium janczewskii, P. digitatum, P. verrucosum, P. crustosum, and P. oxalicum were evaluated to check the antagonistic effect of different dilutions (0, 10, 20, … and 60%). There was a significant reduction of around 7-38%, 46-57%, 46-58%, 27-68%, and 21-51% in Phoma herbarum dry biomass production due to P. janczewskii, P. digitatum, P. verrucosum, P. crustosum, and P. oxalicum, respectively. Inhibition constants determined by a regression equation demonstrated the most significant inhibition by P. janczewskii. Finally, using real-time reverse transcription PCR (qPCR) the effect of P. Janczewskii metabolites was determined on the transcript level of StSTE12 gene involved in the development and penetration of appressorium. The expression pattern of the StSTE12 gene was determined by percent Knockdown (%KD) expression that was found to be decreased i.e. 51.47, 43.22, 40.67, 38.01, 35.97, and 33.41% for P. herbarum with an increase in metabolites concentrations viz., 10, 20, 30, 40, 50 and 60% metabolites, respectively. In silico studies were conducted to analyze the role of Ste12 a transcriptional factor in the MAPK signaling pathway. The present study concludes a strong fungicidal potential of Penicillium species against P. herbarum. Further studies to isolate the effective fungicidal constituents of Penicillium species through GCMS analysis and determination of their role in signaling pathways are requisite.

Polarization and cell-fate decision facilitated by the adaptor Ste50p in Saccharomyces cerevisiae

In response to pheromone, many proteins localize on the plasma membrane of yeast cell to reform it into a polarized shmoo structure. The adaptor protein Ste50p, known as a pheromone signal enhancer critical for shmoo polarization, has never been explored systematically for its localization and function in the polarization process. Time-lapse single-cell imaging and quantitation shown here characterizes Ste50p involvement in the establishment of cell polarity. We found that Ste50p patches on the cell cortex mark the point of shmoo initiation, these patches could move, and remain associated with the growing shmoo tip in a pheromone concentration time-dependent manner until shmoo maturation. A Ste50p mutant impaired in patch localization suffers a delay in polarization. By quantitative analysis we show that polarization correlates with the rising levels of Ste50p, enabling rapid cell responses to pheromone that correspond to a critical level of Ste50p at the initial G1 phase. We exploited the quantitative differences in the pattern of Ste50p expression to correlate with the cell-cell phenotypic heterogeneity, showing Ste50p involvement in the cellular differentiation choice. Taken together, these findings present Ste50p to be part of the early shmoo development phase, suggesting that Ste50p may be involved with the polarisome in the initiation of polarization, and plays a role in regulating the polarized growth of shmoo during pheromone response.

Prey sensing and response in a nematode-trapping fungus is governed by the MAPK pheromone response pathway

Detection of surrounding organisms in the environment plays a major role in the evolution of interspecies interactions, such as predator-prey relationships. Nematode-trapping fungi (NTF) are predators that develop specialized trap structures to capture, kill, and consume nematodes when food sources are limited. Despite the identification of various factors that induce trap morphogenesis, the mechanisms underlying the differentiation process have remained largely unclear. Here, we demonstrate that the highly conserved pheromone-response MAPK pathway is essential for sensing ascarosides, a conserved molecular signature of nemaotdes, and is required for the predatory lifestyle switch in the NTF Arthrobotrys oligospora. Gene deletion of STE7 (MAPKK) and FUS3 (MAPK) abolished nematode-induced trap morphogenesis and conidiation and impaired the growth of hyphae. The conserved transcription factor Ste12 acting downstream of the pheromone-response pathway also plays a vital role in the predation of A. oligospora. Transcriptional profiling of a ste12 mutant identified a small subset of genes with diverse functions that are Ste12 dependent and could trigger trap differentiation. Our work has revealed that A. oligospora perceives and interprets the ascarosides produced by nematodes via the conserved pheromone signaling pathway in fungi, providing molecular insights into the mechanisms of communication between a fungal predator and its nematode prey.

Turning Inside Out: Filamentous Fungal Secretion and Its Applications in Biotechnology, Agriculture, and the Clinic

Filamentous fungi are found in virtually every marine and terrestrial habitat. Vital to this success is their ability to secrete a diverse range of molecules, including hydrolytic enzymes, organic acids, and small molecular weight natural products. Industrial biotechnologists have successfully harnessed and re-engineered the secretory capacity of dozens of filamentous fungal species to make a diverse portfolio of useful molecules. The study of fungal secretion outside fermenters, e.g., during host infection or in mixed microbial communities, has also led to the development of novel and emerging technological breakthroughs, ranging from ultra-sensitive biosensors of fungal disease to the efficient bioremediation of polluted environments. In this review, we consider filamentous fungal secretion across multiple disciplinary boundaries (e.g., white, green, and red biotechnology) and product classes (protein, organic acid, and secondary metabolite). We summarize the mechanistic understanding for how various molecules are secreted and present numerous applications for extracellular products. Additionally, we discuss how the control of secretory pathways and the polar growth of filamentous hyphae can be utilized in diverse settings, including industrial biotechnology, agriculture, and the clinic.

New Aspects of Invasive Growth Regulation Identified by Functional Profiling of MAPK Pathway Targets in Saccharomyces cerevisiae

MAPK pathways are drivers of morphogenesis and stress responses in eukaryotes. A major function of MAPK pathways is the transcriptional induction of target genes, which produce proteins that collectively generate a cellular response. One approach to comprehensively understand how MAPK pathways regulate cellular responses is to characterize the individual functions of their transcriptional targets. Here, by examining uncharacterized targets of the MAPK pathway that positively regulates filamentous growth in Saccharomyces cerevisiae (fMAPK pathway), we identified a new role for the pathway in negatively regulating invasive growth. Specifically, four targets were identified that had an inhibitory role in invasive growth: RPI1, RGD2, TIP1, and NFG1/YLR042cNFG1 was a highly induced unknown open reading frame that negatively regulated the filamentous growth MAPK pathway. We also identified SFG1, which encodes a transcription factor, as a target of the fMAPK pathway. Sfg1p promoted cell adhesion independently from the fMAPK pathway target and major cell adhesion flocculin Flo11p, by repressing genes encoding presumptive cell-wall-degrading enzymes. Sfg1p also contributed to FLO11 expression. Sfg1p and Flo11p regulated different aspects of cell adhesion, and their roles varied based on the environment. Sfg1p also induced an elongated cell morphology, presumably through a cell-cycle delay. Thus, the fMAPK pathway coordinates positive and negative regulatory proteins to fine-tune filamentous growth resulting in a nuanced response. Functional analysis of other pathways' targets may lead to a more comprehensive understanding of how signaling cascades generate biological responses.

Functional Properties of the MAP Kinase UeKpp2 in Ustilago esculenta

Ustilago esculenta undergoes an endophytic life cycle in Zizania latifolia. It induces the stem of its host to swell, forming the edible galls called jiaobai in China, which are the second most commonly cultivated aquatic vegetable in China. Z. latifolia raised for jiaobai can only reproduce asexually because the U. esculenta infection completely inhibits flowering. The infection and proliferation in the host plants during the formation of edible gall differ from those of conventional pathogens. Previous studies have shown a close relationship between mitogen-activated protein kinase (MAPK) and fungal pathogenesis. In this study, we explored the functional properties of the MAPK UeKpp2. Cross-species complementation assays were carried out, which indicated a functional complementation between the UeKpp2 of U. esculenta and the Kpp2 of Ustilago maydis. Next, UeKpp2 mutants of the UeT14 and the UeT55 sporidia background were generated; these showed an aberrant morphology of budding cells, and attenuated mating and filamentous growth in vitro, in the context of normal pathogenicity. Interestingly, we identified another protein kinase, UeUkc1, which acted downstream of UeKpp2 and may participate in the regulation of cell shape. We also found a defect of filamentous growth in UeKpp2 mutants that was not related to a defect of the induction of mating-type genes but was directly related to a defect in UeRbf1 induction. Overall, our results indicate an important role for UeKpp2 in U. esculenta that is slightly different from those reported for other smut fungi.

Transcription factor regulatory modules provide the molecular mechanisms for functional redundancy observed among transcription factors in yeast

BACKGROUND

Current technologies for understanding the transcriptional reprogramming in cells include the transcription factor (TF) chromatin immunoprecipitation (ChIP) experiments and the TF knockout experiments. The ChIP experiments show the binding targets of TFs against which the antibody directs while the knockout techniques find the regulatory gene targets of the knocked-out TFs. However, it was shown that these two complementary results contain few common targets. Researchers have used the concept of TF functional redundancy to explain the low overlap between these two techniques. But the detailed molecular mechanisms behind TF functional redundancy remain unknown. Without knowing the possible molecular mechanisms, it is hard for biologists to fully unravel the cause of TF functional redundancy.

RESULTS

To mine out the molecular mechanisms, a novel algorithm to extract TF regulatory modules that help explain the observed TF functional redundancy effect was devised and proposed in this research. The method first searched for candidate TF sets from the TF binding data. Then based on these candidate sets the method utilized the modified Steiner Tree construction algorithm to construct the possible TF regulatory modules from protein-protein interaction data and finally filtered out the noise-induced results by using confidence tests. The mined-out regulatory modules were shown to correlate to the concept of functional redundancy and provided testable hypotheses of the molecular mechanisms behind functional redundancy. And the biological significance of the mined-out results was demonstrated in three different biological aspects: ontology enrichment, protein interaction prevalence and expression coherence. About 23.5% of the mined-out TF regulatory modules were literature-verified. Finally, the biological applicability of the proposed method was shown in one detailed example of a verified TF regulatory module for pheromone response and filamentous growth in yeast.

CONCLUSION

In this research, a novel method that mined out the potential TF regulatory modules which elucidate the functional redundancy observed among TFs is proposed. The extracted TF regulatory modules not only correlate the molecular mechanisms to the observed functional redundancy among TFs, but also show biological significance in inferring TF functional binding target genes. The results provide testable hypotheses for biologists to further design subsequent research and experiments.

MADS-box transcription factor Mcm1 controls cell cycle, fungal development, cell integrity and virulence in the filamentous insect pathogenic fungus Beauveria bassiana

MADS-box transcription factor Mcm1 plays crucial roles in regulating mating processes and pathogenesis in some fungi. However, its roles are varied in fungal species, and its function remains unclear in entomopathogenic fungi. Here, Mcm1 orthologue, Bbmcm1, was characterized in a filamentous entomopathogenic fungus, Beauveria bassiana. Disruption of Bbmcm1 resulted in a distinct reduction in growth with abnormal conidiogenesis, and a significant decrease in conidial viability with abnormal germination. ΔBbmcm1 displayed impaired cell integrity, with distorted cell wall structure and altered cell wall component. Abnormal cell cycle was detected in ΔBbmcm1 with longer G₂ /M phase but shorter G₁ /G₀ and S phases in unicellular blastospores, and sparser septa in multicellular hyphae, which might be responsible for defects in development and differentiation due to the regulation of cell cycle-involved genes, as well as the corresponding cellular events-associated genes. Dramatically decreased virulence was examined in ΔBbmcm1, with impaired ability to escape haemocyte encapsulation, which was consistent with markedly reduced cuticle-degrading enzyme production by repressing their coding genes, and downregulated fungal effector protein-coding genes, suggesting a novel role of Mcm1 in interaction with host insect. These data indicate that Mcm1 is a crucial regulator of development, cell integrity, cell cycle and virulence in insect fungal pathogens.

Filamentation Regulatory Pathways Control Adhesion-Dependent Surface Responses in Yeast

Signaling pathways can regulate biological responses by the transcriptional regulation of target genes. In yeast, multiple signaling pathways control filamentous growth, a morphogenetic response that occurs in many species including fungal pathogens. Here, we examine the role of signaling pathways that control filamentous growth in regulating adhesion-dependent surface responses, including mat formation and colony patterning. Expression profiling and mutant phenotype analysis showed that the major pathways that regulate filamentous growth [filamentous growth MAPK (fMAPK), RAS, retrograde (RTG), RIM101, RPD3, ELP, SNF1, and PHO85] also regulated mat formation and colony patterning. The chromatin remodeling complex, SAGA, also regulated these responses. We also show that the RAS and RTG pathways coregulated a common set of target genes, and that SAGA regulated target genes known to be controlled by the fMAPK, RAS, and RTG pathways. Analysis of surface growth-specific targets identified genes that respond to low oxygen, high temperature, and desiccation stresses. We also explore the question of why cells make adhesive contacts in colonies. Cell adhesion contacts mediated by the coregulated target and adhesion molecule, Flo11p, deterred entry into colonies by macroscopic predators and impacted colony temperature regulation. The identification of new regulators (e.g., SAGA), and targets of surface growth in yeast may provide insights into fungal pathogenesis in settings where surface growth and adhesion contributes to virulence.

Transcription Factor SsSte12 Was Involved in Mycelium Growth and Development in Sclerotinia sclerotiorum

Sclerotinia sclerotiorum is a challenging agricultural pathogen for management, causing large global economic losses annually. The sclerotia and infection cushions are critical for its long-term survival and successful penetration on a wide spectrum of hosts. The mitogen-activated protein kinase (MAPK) cascades serve as central signaling complexes that are involved in various aspects of sclerotia development and infection. In this study, the putative downstream transcription factor of MAPK pathway, SsSte12, was analyzed in S. sclerotiorum. Silencing SsSte12 in S. sclerotiorum resulted in phenotypes of delayed vegetative growth, reduced size of sclerotia, and fewer appressoria formation. Consequently, the SsSte12 RNAi mutants showed attenuated pathogenicity on the host plants due to the defect compound appressorium. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation assays demonstrated that the SsSte12 interacts with SsMcm1. However, the SsMcm1 expression is independent of the regulation of SsSte12 as revealed by qRT-PCR analysis in SsSte12 RNAi mutants. Together with high accumulation of SsSte12 transcripts in the early development of S. sclerotiorum, our results demonstrated that SsSte12 function was essential in the vegetative mycelial growth, sclerotia development, appressoria formation and penetration-dependent pathogenicity. Moreover, the SsSte12-SsMcm1 interaction might play a critical role in the regulation of the genes encoding these traits in S. sclerotiorum.

Preferences in a trait decision determined by transcription factor variants

Few mechanisms are known that explain how transcription factors can adjust phenotypic outputs to accommodate differing environments. In Saccharomyces cerevisiae, the decision to mate or invade relies on environmental cues that converge on a shared transcription factor, Ste12. Specificity toward invasion occurs via Ste12 binding cooperatively with the cofactor Tec1. Here, we determine the range of phenotypic outputs (mating vs. invasion) of thousands of DNA-binding domain variants in Ste12 to understand how preference for invasion may arise. We find that single amino acid changes in the DNA-binding domain can shift the preference of yeast toward either mating or invasion. These mutations define two distinct regions of this domain, suggesting alternative modes of DNA binding for each trait. We characterize the DNA-binding specificity of wild-type Ste12 to identify a strong preference for spacing and orientation of both homodimeric and heterodimeric sites. Ste12 mutants that promote hyperinvasion in a Tec1-independent manner fail to bind cooperative sites with Tec1 and bind to unusual dimeric Ste12 sites composed of one near-perfect and one highly degenerate site. We propose a model in which Ste12 alone may have evolved to activate invasion genes, which could explain how preference for invasion arose in the many fungal pathogens that lack Tec1.

Timing of gene expression in a cell-fate decision system

During development, morphogens provide extracellular cues allowing cells to select a specific fate by inducing complex transcriptional programs. The mating pathway in budding yeast offers simplified settings to understand this process. Pheromone secreted by the mating partner triggers the activity of a MAPK pathway, which results in the expression of hundreds of genes. Using a dynamic expression reporter, we quantified the kinetics of gene expression in single cells upon exogenous pheromone stimulation and in the physiological context of mating. In both conditions, we observed striking differences in the timing of induction of mating‐responsive promoters. Biochemical analyses and generation of synthetic promoter variants demonstrated how the interplay between transcription factor binding and nucleosomes contributes to determine the kinetics of transcription in a simplified cell‐fate decision system.

The transcription factor PstSTE12 is required for virulence of Puccinia striiformis f. sp. tritici

Puccinia striiformis f. sp. tritici (Pst) is an obligate biotrophic fungus that causes extensive damage in wheat. The pathogen is now known to be a heteroecious fungus with an intricate life cycle containing sexual and asexual stages. Orthologues of the STE12 transcription factor that regulate mating and filamentation in Saccharomyces cerevisiae, as well as virulence in other fungi, have been extensively described. Because reliable transformation and gene disruption methods are lacking for Pst, knowledge about the function of its STE12 orthologue is limited. In this study, we identified a putative orthologue of STE12 from Pst in haustoria-enriched transcripts and designated it as PstSTE12. The gene encodes a protein of 879 amino acids containing three helices in the homeodomain, conserved phenylalanine and tryptophan sites, and two C₂ /H₂ -Zn²⁺ finger domains. Real-time reverse transcription-polymerase chain reaction (RT-PCR) analyses revealed that the expression of PstSTE12 was highly induced during the early infection stages and peaked during haustorium formation and the pycniospore stage in the aecial host barberry. Subcellular localization assays indicated that PstSTE12 is localized in the nucleus and functions as a transcriptional activator. Yeast one-hybrid assays revealed that PstSTE12 exhibits transcriptional activity, and that its C-terminus is necessary for the activation of transcription. PstSTE12 complemented the mating defect in an α ste12 mutant of S. cerevisiae. In addition, it partially complemented the defects of the Magnaporthe oryzae mst12 mutant in plant infection. Knocking down PstSTE12 via host-induced gene silencing (HIGS) mediated by Barley stripe mosaic virus (BSMV) resulted in a substantial reduction in the growth and spread of hyphae in Pst and weakened the virulence of Pst on wheat. Our results suggest that PstSTE12 probably acts at an intersection participating in the invasion and mating processes of Pst, and provide new insights into the comprehension of the variation of virulence in cereal rust fungi.

Analysis of the thresholds for transcriptional activation by the yeast MAP kinases Fus3 and Kss1

Signaling in the pheromone response pathway of budding yeast activates two distinct MAP kinases (MAPKs), Fus3 and Kss1. Either MAPK alone can mediate pheromone-induced transcription, but it has been unclear to what degree each one contributes to transcriptional output in wild-type cells. Here, we report that transcription reflects the ratio of active to inactive MAPK, and not simply the level of active MAPK. For Kss1 the majority of MAPK molecules must be converted to the active form, whereas for Fus3 only a small minority must be activated. These different activation thresholds reflect two opposing effects of each MAPK, in which the inactive forms inhibit transcription, whereas the active forms promote transcription. Moreover, negative feedback from Fus3 limits activation of Kss1 so that it does not meet its required threshold in wild-type cells but does so only when hyperactivated in cells lacking Fus3. The results suggest that the normal transcriptional response involves asymmetric contributions from the two MAPKs, in which pheromone signaling reduces the negative effect of Kss1 while increasing the positive effect of Fus3. These findings reveal new functional distinctions between these MAPKs, and help illuminate how inhibitory functions shape positive pathway outputs in both pheromone and filamentation pathways.

From Genes to Networks: The Regulatory Circuitry Controlling Candida albicans Morphogenesis

Candida albicans is a commensal yeast of most healthy individuals, but also one of the most prevalent human fungal pathogens. During adaptation to the mammalian host, C. albicans encounters different niches where it is exposed to several types of stress, including oxidative, nitrosative (e.g., immune system), osmotic (e.g., kidney and oral cavity) stresses and pH variation (e.g., gastrointestinal (GI) tract and vagina). C. albicans has developed the capacity to respond to the environmental changes by modifying its morphology, which comprises the yeast-to-hypha transition, white-opaque switching, and chlamydospore formation. The yeast-to-hypha transition has been very well characterized and was shown to be modulated by several external stimuli that mimic the host environment. For instance, temperature above 37 ℃, serum, alkaline pH, and CO₂ concentration are all reported to enhance filamentation. The transition is characterized by the activation of an intricate regulatory network of signaling pathways, involving many transcription factors. The regulatory pathways that control either the stress response or morphogenesis are required for full virulence and promote survival of C. albicans in the host. Many of these transcriptional circuitries have been characterized, highlighting the complexity and the interconnections between the different pathways. Here, we present the major signaling pathways and the main transcription factors involved in the yeast-to-hypha transition. Furthermore, we describe the role of heat shock transcription factors in the morphogenetic transition, providing an edifying example of the complex cross talk between pathways involved in morphogenesis and stress response.

Negative regulation of filamentous growth in Candida albicans by Dig1p

Transcriptional regulation involves both positive and negative regulatory elements. The Dig1 negative regulators are part of a fungal-specific module that includes a transcription factor (a Ste12 family member) and a Dig1 family member. In Saccharomyces cerevisiae, the post-genome-duplication Dig1/Dig2 proteins regulate MAP kinase controlled signalling pathways involved in mating and filamentous growth. We have identified the single Dig1 orthologue in the fungal pathogen Candida albicans. Genetic studies and transcriptional profiling experiments show that this single protein is implicated in the regulation of MAP kinase-controlled processes involved in mating, filamentous growth and biofilm formation, and also influences cAMP-regulated processes. This suggests that the multiple cellular roles of the Dig1 protein are ancestral and predate the sub-functionalization apparent in S. cerevisiae after the genome duplication. Intriguingly, even though loss of Dig1 function in C. albicans enhances filamentous growth and biofilm formation, colonization of the murine gastrointestinal tract is reduced in the mutant. The complexity of the processes influenced by Dig1 in C. albicans, and the observation that Dig1 is one of the few regulatory proteins that were retained in the duplicated state after the whole genome duplication event in yeast, emphasizes the important role of these negative regulators in fungal transcriptional control.

Fungal Sex: The Ascomycota

This article provides an overview of sexual reproduction in the ascomycetes, a phylum of fungi that is named after the specialized sacs or "asci" that hold the sexual spores. They have therefore also been referred to as the Sac Fungi due to these characteristic structures that typically contain four to eight ascospores. Ascomycetes are morphologically diverse and include single-celled yeasts, filamentous fungi, and more complex cup fungi. The sexual cycles of many species, including those of the model yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe and the filamentous saprobes Neurospora crassa, Aspergillus nidulans, and Podospora anserina, have been examined in depth. In addition, sexual or parasexual cycles have been uncovered in important human pathogens such as Candida albicans and Aspergillus fumigatus, as well as in plant pathogens such as Fusarium graminearum and Cochliobolus heterostrophus. We summarize what is known about sexual fecundity in ascomycetes, examine how structural changes at the mating-type locus dictate sexual behavior, and discuss recent studies that reveal that pheromone signaling pathways can be repurposed to serve cellular roles unrelated to sex.

Transcriptional Profiling of Biofilm Regulators Identified by an Overexpression Screen in Saccharomyces cerevisiae

Biofilm formation by microorganisms is a major cause of recurring infections and removal of biofilms has proven to be extremely difficult given their inherent drug resistance . Understanding the biological processes that underlie biofilm formation is thus extremely important and could lead to the development of more effective drug therapies, resulting in better infection outcomes. Using the yeast Saccharomyces cerevisiae as a biofilm model, overexpression screens identified DIG1, SFL1, HEK2, TOS8, SAN1, and ROF1/YHR177W as regulators of biofilm formation. Subsequent RNA-seq analysis of biofilm and nonbiofilm-forming strains revealed that all of the overexpression strains, other than DIG1 and TOS8, were adopting a single differential expression profile, although induced to varying degrees. TOS8 adopted a separate profile, while the expression profile of DIG1 reflected the common pattern seen in most of the strains, plus substantial DIG1-specific expression changes. We interpret the existence of the common transcriptional pattern seen across multiple, unrelated overexpression strains as reflecting a transcriptional state, that the yeast cell can access through regulatory signaling mechanisms, allowing an adaptive morphological change between biofilm-forming and nonbiofilm states.

Dissecting Gene Expression Changes Accompanying a Ploidy-Based Phenotypic Switch

Aneuploidy, a state in which the chromosome number deviates from a multiple of the haploid count, significantly impacts human health. The phenotypic consequences of aneuploidy are believed to arise from gene expression changes associated with the altered copy number of genes on the aneuploid chromosomes. To dissect the mechanisms underlying altered gene expression in aneuploids, we used RNA-seq to measure transcript abundance in colonies of the haploid Saccharomyces cerevisiae strain F45 and two aneuploid derivatives harboring disomies of chromosomes XV and XVI. F45 colonies display complex “fluffy” morphologies, while the disomic colonies are smooth, resembling laboratory strains. Our two disomes displayed similar transcriptional profiles, a phenomenon not driven by their shared smooth colony morphology nor simply by their karyotype. Surprisingly, the environmental stress response (ESR) was induced in F45, relative to the two disomes. We also identified genes whose expression reflected a nonlinear interaction between the copy number of a transcriptional regulatory gene on chromosome XVI, DIG1, and the copy number of other chromosome XVI genes. DIG1 and the remaining chromosome XVI genes also demonstrated distinct contributions to the effect of the chromosome XVI disome on ESR gene expression. Expression changes in aneuploids appear to reflect a mixture of effects shared between different aneuploidies and effects unique to perturbing the copy number of particular chromosomes, including nonlinear copy number interactions between genes. The balance between these two phenomena is likely to be genotype- and environment-specific.

Spatial and temporal signal processing and decision making by MAPK pathways

Recent studies show that MAPK pathways perform exquisite spatial and temporal signal processing. This review discusses the mechanisms that process dynamic inputs into graded output responses, the role of positive and negative feedbacks, and feedforward regulation.

Mitogen-activated protein kinase (MAPK) pathways are conserved from yeast to man and regulate a variety of cellular processes, including proliferation and differentiation. Recent developments show how MAPK pathways perform exquisite spatial and temporal signal processing and underscores the importance of studying the dynamics of signaling pathways to understand their physiological response. The importance of dynamic mechanisms that process input signals into graded downstream responses has been demonstrated in the pheromone-induced and osmotic stress–induced MAPK pathways in yeast and in the mammalian extracellular signal-regulated kinase MAPK pathway. Particularly, recent studies in the yeast pheromone response have shown how positive feedback generates switches, negative feedback enables gradient detection, and coherent feedforward regulation underlies cellular memory. More generally, a new wave of quantitative single-cell studies has begun to elucidate how signaling dynamics determine cell physiology and represents a paradigm shift from descriptive to predictive biology.

Protein expression patterns of the yeast mating response

Microfluidics, in combination with time-lapse microscopy, is a transformative technology that significantly enhances our ability to monitor and probe biological processes in living cells. However, high-throughput microfluidic devices mostly require sophisticated preparatory and setup work and are thus hard to adopt by non-experts. In this work, we designed an easy-to-use microfluidic chip, which enables tracking of 48 GFP-tagged yeast strains, with each strain under two different stimulus conditions, in a single experiment. We used this technology to investigate the dynamic pattern of protein expression during the yeast mating differentiation response. High doses of pheromone induce cell cycle arrest and the shmoo morphology, whereas low doses of pheromone lead to elongation and chemotrophic growth. By systematically analyzing the protein dynamics of 156 pheromone-regulated genes, we identified groups of genes that are preferentially induced in response to low-dose pheromone (elongation during growth) or high-dose pheromone (shmoo formation and cell cycle arrest). The protein dynamics of these genes may provide insights into the mechanisms underlying the differentiation switch induced by different doses of pheromone.

The transcriptional regulators SteA and StuA contribute to keratin degradation and sexual reproduction of the dermatophyte Arthroderma benhamiae

Most superficial fungal infections are caused by dermatophytes, a specialized group of filamentous fungi which exclusively infect keratinized host structures such as hair, skin and nails. Since little is known about the molecular basis of pathogenicity and sexual reproduction in dermatophytes, here we functionally addressed two central transcriptional regulators, SteA and StuA. In the zoophilic species Arthroderma benhamiae a strategy for targeted genetic manipulation was recently established, and moreover, the species is teleomorphic and thus allows performing assays based on mating. By comparative genome analysis homologs of the developmental regulators SteA and StuA were identified in A. benhamiae. Knock-out mutants of the corresponding genes as well as complemented strains were generated and phenotypically characterized. In contrast to A. benhamiae wild type and complemented strains, both mutants failed to produce sexual reproductive structures in mating experiments. Analysis of growth on keratin substrates indicated that loss of steA resulted in the inability of ΔsteA mutants to produce hair perforation organs, but did not affect mycelia formation during growth on hair and nails. By contrast, ΔstuA mutants displayed a severe growth defect on these substrates, but were still able to produce hair perforations. Hence, formation of hair perforation organs and fungal growth on hair per se are differentially regulated processes. Our findings on the major role of SteA and StuA during sexual development and keratin degradation in A. benhamiae provide insights into their role in dermatophytes and further enhance our knowledge of basic biology and pathogenicity of these fungi.

cDREM: inferring dynamic combinatorial gene regulation

Genes are often combinatorially regulated by multiple transcription factors (TFs). Such combinatorial regulation plays an important role in development and facilitates the ability of cells to respond to different stresses. While a number of approaches have utilized sequence and ChIP-based datasets to study combinational regulation, these have often ignored the combinational logic and the dynamics associated with such regulation. Here we present cDREM, a new method for reconstructing dynamic models of combinatorial regulation. cDREM integrates time series gene expression data with (static) protein interaction data. The method is based on a hidden Markov model and utilizes the sparse group Lasso to identify small subsets of combinatorially active TFs, their time of activation, and the logical function they implement. We tested cDREM on yeast and human data sets. Using yeast we show that the predicted combinatorial sets agree with other high throughput genomic datasets and improve upon prior methods developed to infer combinatorial regulation. Applying cDREM to study human response to flu, we were able to identify several combinatorial TF sets, some of which were known to regulate immune response while others represent novel combinations of important TFs.

Integrative random forest for gene regulatory network inference

Motivation: Gene regulatory network (GRN) inference based on genomic data is one of the most actively pursued computational biological problems. Because different types of biological data usually provide complementary information regarding the underlying GRN, a model that integrates big data of diverse types is expected to increase both the power and accuracy of GRN inference. Towards this goal, we propose a novel algorithm named iRafNet: integrative random forest for gene regulatory network inference.

Results: iRafNet is a flexible, unified integrative framework that allows information from heterogeneous data, such as protein–protein interactions, transcription factor (TF)-DNA-binding, gene knock-down, to be jointly considered for GRN inference. Using test data from the DREAM4 and DREAM5 challenges, we demonstrate that iRafNet outperforms the original random forest based network inference algorithm (GENIE3), and is highly comparable to the community learning approach. We apply iRafNet to construct GRN in Saccharomyces cerevisiae and demonstrate that it improves the performance in predicting TF-target gene regulations and provides additional functional insights to the predicted gene regulations.

Availability and implementation: The R code of iRafNet implementation and a tutorial are available at: http://research.mssm.edu/tulab/software/irafnet.html

Contact:zhidong.tu@mssm.edu

Supplementary information: Supplementary data are available at Bioinformatics online.

Pheromone-induced morphogenesis and gradient tracking are dependent on the MAPK Fus3 binding to Gα

Unique roles are found for the MAPK Fus3 during the mating response of yeast. In particular, the interaction of Fus3 with the G-protein α-subunit is required for morphogenesis and gradient tracking and suppresses cell-to-cell variability between mating and chemotropic fates in a population of pheromone-responding cells.

Mitogen-activated protein kinase (MAPK) pathways control many cellular processes, including differentiation and proliferation. These pathways commonly activate MAPK isoforms that have redundant or overlapping function. However, recent studies have revealed circumstances in which MAPK isoforms have specialized, nonoverlapping roles in differentiation. The mechanisms that underlie this specialization are not well understood. To address this question, we sought to establish regulatory mechanisms that are unique to the MAPK Fus3 in pheromone-induced mating and chemotropic fate transitions of the budding yeast Saccharomyces cerevisiae. Our investigations reveal a previously unappreciated role for inactive Fus3 as a potent negative regulator of pheromone-induced chemotropism. We show that this inhibitory role is dependent on inactive Fus3 binding to the α-subunit of the heterotrimeric G-protein. Further analysis revealed that the binding of catalytically active Fus3 to the G-protein is required for gradient tracking and serves to suppress cell-to-cell variability between mating and chemotropic fates in a population of pheromone-responding cells.

Fine-tuning of histone H3 Lys4 methylation during pseudohyphal differentiation by the CDK submodule of RNA polymerase II

Transcriptional regulation is dependent upon the interactions between the RNA pol II holoenzyme complex and chromatin. RNA pol II is part of a highly conserved multiprotein complex that includes the core mediator and CDK8 subcomplex. In Saccharomyces cerevisiae, the CDK8 subcomplex, composed of Ssn2p, Ssn3p, Ssn8p, and Srb8p, is thought to play important roles in mediating transcriptional control of stress-responsive genes. Also central to transcriptional control are histone post-translational modifications. Lysine methylation, dynamically balanced by lysine methyltransferases and demethylases, has been intensively studied, uncovering significant functions in transcriptional control. A key question remains in understanding how these enzymes are targeted during stress response. To determine the relationship between lysine methylation, the CDK8 complex, and transcriptional control, we performed phenotype analyses of yeast lacking known lysine methyltransferases or demethylases in isolation or in tandem with SSN8 deletions. We show that the RNA pol II CDK8 submodule components SSN8/SSN3 and the histone demethylase JHD2 are required to inhibit pseudohyphal growth-a differentiation pathway induced during nutrient limitation-under rich conditions. Yeast lacking both SSN8 and JHD2 constitutively express FLO11, a major regulator of pseudohyphal growth. Interestingly, deleting known FLO11 activators including FLO8, MSS11, MFG1, TEC1, SNF1, KSS1, and GCN4 results in a range of phenotypic suppression. Using chromatin immunoprecipitation, we found that SSN8 inhibits H3 Lys4 trimethylation independently of JHD2 at the FLO11 locus, suggesting that H3 Lys4 hypermethylation is locking FLO11 into a transcriptionally active state. These studies implicate the CDK8 subcomplex in fine-tuning H3 Lys4 methylation levels during pseudohyphal differentiation.

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.

Global regulation of a differentiation MAPK pathway in yeast

Cell differentiation requires different pathways to act in concert to produce a specialized cell type. The budding yeast Saccharomyces cerevisiae undergoes filamentous growth in response to nutrient limitation. Differentiation to the filamentous cell type requires multiple signaling pathways, including a mitogen-activated protein kinase (MAPK) pathway. To identify new regulators of the filamentous growth MAPK pathway, a genetic screen was performed with a collection of 4072 nonessential deletion mutants constructed in the filamentous (Σ1278b) strain background. The screen, in combination with directed gene-deletion analysis, uncovered 97 new regulators of the filamentous growth MAPK pathway comprising 40% of the major regulators of filamentous growth. Functional classification extended known connections to the pathway and identified new connections. One function for the extensive regulatory network was to adjust the activity of the filamentous growth MAPK pathway to the activity of other pathways that regulate the response. In support of this idea, an unregulated filamentous growth MAPK pathway led to an uncoordinated response. Many of the pathways that regulate filamentous growth also regulated each other's targets, which brings to light an integrated signaling network that regulates the differentiation response. The regulatory network characterized here provides a template for understanding MAPK-dependent differentiation that may extend to other systems, including fungal pathogens and metazoans.

Identification of important nodes in directed biological networks: a network motif approach

Identification of important nodes in complex networks has attracted an increasing attention over the last decade. Various measures have been proposed to characterize the importance of nodes in complex networks, such as the degree, betweenness and PageRank. Different measures consider different aspects of complex networks. Although there are numerous results reported on undirected complex networks, few results have been reported on directed biological networks. Based on network motifs and principal component analysis (PCA), this paper aims at introducing a new measure to characterize node importance in directed biological networks. Investigations on five real-world biological networks indicate that the proposed method can robustly identify actually important nodes in different networks, such as finding command interneurons, global regulators and non-hub but evolutionary conserved actually important nodes in biological networks. Receiver Operating Characteristic (ROC) curves for the five networks indicate remarkable prediction accuracy of the proposed measure. The proposed index provides an alternative complex network metric. Potential implications of the related investigations include identifying network control and regulation targets, biological networks modeling and analysis, as well as networked medicine.

The transcription factors Tec1 and Ste12 interact with coregulators Msa1 and Msa2 to activate adhesion and multicellular development

In Saccharomyces cerevisiae and related yeast species, the TEA transcription factor Tec1, together with a second transcription factor, Ste12, controls development, including cell adhesion and filament formation. Tec1-Ste12 complexes control target genes through Tec1 binding sites (TEA consensus sequences [TCSs]) that can be further combined with Ste12 binding sites (pheromone response elements [PREs]) for cooperative DNA binding. The activity of Tec1-Ste12 complexes is known to be under negative control of the Dig1 and Dig2 (Dig1/2) transcriptional corepressors that confer regulation by upstream signaling pathways. Here, we found that Tec1 and Ste12 can associate with the transcriptional coregulators Msa1 and Msa2 (Msa1/2), which were previously found to associate with the cell cycle transcription factor complexes SBF (Swi4/Swi6 cell cycle box binding factor) and MBF (Mbp1/Swi6 cell cycle box binding factor) to control G1-specific transcription. We further show that Tec1-Ste12-Msa1/2 complexes (i) do not contain Swi4 or Mbp1, (ii) assemble at single TCSs or combined TCS-PREs in vitro, and (iii) coregulate genes involved in adhesive and filamentous growth by direct promoter binding in vivo. Finally, we found that, in contrast to Dig proteins, Msa1/2 seem to act as coactivators that enhance the transcriptional activity of Tec1-Ste12. Taken together, our findings add an additional layer of complexity to our understanding of the control mechanisms exerted by the evolutionarily conserved TEA domain and Ste12-like transcription factors.

Proteomic identification of potential target proteins regulated by the SCF(F) (bp1) -mediated proteolysis pathway in Fusarium oxysporum

F-box proteins function in the recruitment of proteins for SCF ubiquitination and proteasome degradation. Here, we studied the role of Fbp1, a nonessential F-box protein of the tomato pathogen Fusarium oxysporum f. sp. lycopersici. The Δfbp1 mutant showed a significant delay in the production of wilt symptoms on tomato plants and was impaired in invasive growth on cellophane membranes and on living plant tissue. To search for target proteins recruited by Fbp1, a combination of sodium dodecylsulphate-polyacrylamide gel electrophoresis (SDS-PAGE) and matrix-assisted laser desorption/ionization time-of-flight/time-of-flight (MALDI-TOF/TOF) was used to compare proteins in mycelia of the wild-type and Δfbp1 mutant. The proteomic approach identified 41 proteins differing significantly in abundance between the two strains, 17 of which were more abundant in the Δfbp1 mutant, suggesting a possible regulation by proteasome degradation. Interestingly, several of the identified proteins were related to vesicle trafficking. Microscopic analysis revealed an impairment of the Δfbp1 strain in directional growth and in the structure of the Spitzenkörper, suggesting a role of Fbp1 in hyphal orientation. Our results indicate that Fbp1 regulates protein turnover and pathogenicity in F. oxysporum.

Mating type gene homologues and putative sex pheromone-sensing pathway in arbuscular mycorrhizal fungi, a presumably asexual plant root symbiont

The fungal kingdom displays a fascinating diversity of sex-determination systems. Recent advances in genomics provide insights into the molecular mechanisms of sex, mating type determination, and evolution of sexual reproduction in many fungal species in both ancient and modern phylogenetic lineages. All major fungal groups have evolved sexual differentiation and recombination pathways. However, sexuality is unknown in arbuscular mycorrhizal fungi (AMF) of the phylum Glomeromycota, an ecologically vital group of obligate plant root symbionts. AMF are commonly considered an ancient asexual lineage dating back to the Ordovician, approximately 460 M years ago. In this study, we used genomic and transcriptomic surveys of several AMF species to demonstrate the presence of conserved putative sex pheromone-sensing mitogen-activated protein (MAP) kinases, comparable to those described in Ascomycota and Basidiomycota. We also find genes for high mobility group (HMG) transcription factors, homologous to SexM and SexP genes in the Mucorales. The SexM genes show a remarkable sequence diversity among multiple copies in the genome, while only a single SexP sequence was detected in some isolates of Rhizophagus irregularis. In the Mucorales and Microsporidia, the sexM gene is flanked by genes for a triosephosphate transporter (TPT) and a RNA helicase, but we find no evidence for synteny in the vicinity of the Sex locus in AMF. Nonetheless, our results, together with previous observations on meiotic machinery, suggest that AMF could undergo a complete sexual reproduction cycle.

Multiple MAPK cascades regulate the transcription of IME1, the master transcriptional activator of meiosis in Saccharomyces cerevisiae

The choice between alternative developmental pathways is primarily controlled at the level of transcription. Induction of meiosis in budding yeasts in response to nutrient levels provides a system to investigate the molecular basis of cellular decision-making. In Saccharomyces cerevisiae, entry into meiosis depends on multiple signals converging upon IME1, the master transcriptional activator of meiosis. Here we studied the regulation of the cis-acting regulatory element Upstream Activation Signal (UAS)ru, which resides within the IME1 promoter. Guided by our previous data acquired using a powerful high-throughput screening system, here we provide evidence that UASru is regulated by multiple stimuli that trigger distinct signal transduction pathways as follows: (i) The glucose signal inhibited UASru activity through the cyclic AMP (cAMP/protein kinase A (PKA) pathway, targeting the transcription factors (TFs), Com2 and Sko1; (ii) high osmolarity activated UASru through the Hog1/mitogen-activated protein kinase (MAPK) pathway and its corresponding TF Sko1; (iii) elevated temperature increased the activity of UASru through the cell wall integrity pathway and the TFs Swi4/Mpk1 and Swi4/Mlp1; (iv) the nitrogen source repressed UASru activity through Sum1; and (v) the absence of a nitrogen source was detected and transmitted to UASru by the Kss1 and Fus3 MAPK pathways through their respective downstream TFs, Ste12/Tec1 and Ste12/Ste12 as well as by their regulators Dig1/2. These signaling events were specific to UASru; they did not affect the mating and filamentation response elements that are regulated by MAPK pathways. The complex regulation of UASru through all the known vegetative MAPK pathways is unique to S. cerevisiae and is specific for IME1, likely because it is the master regulator of gametogenesis.

Dynamic ubiquitination of the mitogen-activated protein kinase kinase (MAPKK) Ste7 determines mitogen-activated protein kinase (MAPK) specificity

Ubiquitination is a post-translational modification that tags proteins for proteasomal degradation. In addition, there is a growing appreciation that ubiquitination can influence protein activity and localization. Ste7 is a prototype MAPKK in yeast that participates in both the pheromone signaling and nutrient deprivation/invasive growth pathways. We have shown previously that Ste7 is ubiquitinated upon pheromone stimulation. Here, we show that the Skp1/Cullin/F-box ubiquitin ligase SCF(Cdc4) and the ubiquitin protease Ubp3 regulate Ste7 ubiquitination and signal specificity. Using purified components, we demonstrate that SCF(Cdc4) ubiquitinates Ste7 directly. Using gene deletion mutants, we show that SCF(Cdc4) and Ubp3 have opposing effects on Ste7 ubiquitination. Although SCF(Cdc4) is necessary for proper activation of the pheromone MAPK Fus3, Ubp3 is needed to limit activation of the invasive growth MAPK Kss1. Finally, we show that Fus3 phosphorylates Ubp3 directly and that phosphorylation of Ubp3 is necessary to limit Kss1 activation. These results reveal a feedback loop wherein one MAPK limits the ubiquitination of an upstream MAPKK and thereby prevents spurious activation of a second competing MAPK.

Assigning quantitative function to post-translational modifications reveals multiple sites of phosphorylation that tune yeast pheromone signaling output

Cell signaling systems transmit information by post-translationally modifying signaling proteins, often via phosphorylation. While thousands of sites of phosphorylation have been identified in proteomic studies, the vast majority of sites have no known function. Assigning functional roles to the catalog of uncharacterized phosphorylation sites is a key research challenge. Here we present a general approach to address this challenge and apply it to a prototypical signaling pathway, the pheromone response pathway in Saccharomyces cerevisiae. The pheromone pathway includes a mitogen activated protein kinase (MAPK) cascade activated by a G-protein coupled receptor (GPCR). We used published mass spectrometry-based proteomics data to identify putative sites of phosphorylation on pheromone pathway components, and we used evolutionary conservation to assign priority to a list of candidate MAPK regulatory sites. We made targeted alterations in those sites, and measured the effects of the mutations on pheromone pathway output in single cells. Our work identified six new sites that quantitatively tuned system output. We developed simple computational models to find system architectures that recapitulated the quantitative phenotypes of the mutants. Our results identify a number of putative phosphorylation events that contribute to adjust the input-output relationship of this model eukaryotic signaling system. We believe this combined approach constitutes a general means not only to reveal modification sites required to turn a pathway on and off, but also those required for more subtle quantitative effects that tune pathway output. Our results suggest that relatively small quantitative influences from individual phosphorylation events endow signaling systems with plasticity that evolution may exploit to quantitatively tailor signaling outcomes.

RNA preparation of Saccharomyces cerevisiae using the digestion method may give misleading results

Zymolyase (lyticase) is used for cell wall digestion in yeast experiments and is needed for incubation processes under moderate experimental conditions. This has been thought to cause unfavorable effects, and many researchers are aware that the enzyme method is unsuitable for RNA preparation following several reports of stress responses to the enzyme process. However, RNA preparation with enzyme digestion continues to be used. This may be because there have been insufficient data directly comparing RNA preparation conditions with previous studies. We investigated the influence of enzyme processes in RNA preparation using a DNA microarray, and compared superoxide dismutase (SOD) activities with a non-treated control and the results of previous research. Gene expressions were commonly changed by enzyme processes, and SOD activities increased only during short-term incubation. Meanwhile, both SOD gene expressions and SOD activity during RNA preparation indicated different results than gained under conditions of long-term incubation. These results suggest that zymolyase treatment surely influences gene expressions and enzyme activity, although the effect assumed by previous studies is not necessarily in agreement with that of RNA preparation.

Role of transcription at centromeres in budding yeast

Centromeres are specialized chromosomal loci that are essential for proper chromosome segregation. Recent data show that a certain level of active transcription, regulated by transcription factors Cbf1 and Ste12, makes a direct contribution to centromere function in Saccharomyces cerevisiae. Here, we discuss the requirement and function of transcription at centromeres.

Positive roles for negative regulators in the mating response of yeast

All cells must detect and respond to changes in their environment, often through changes in gene expression. The yeast pheromone pathway has been extensively characterized, and is an ideal system for studying transcriptional regulation. Here we combine computational and experimental approaches to study transcriptional regulation mediated by Ste12, the key transcription factor in the pheromone response. Our mathematical model is able to explain multiple counterintuitive experimental results and led to several novel findings. First, we found that the transcriptional repressors Dig1 and Dig2 positively affect transcription by stabilizing Ste12. This stabilization through protein-protein interactions creates a large pool of Ste12 that is rapidly activated following pheromone stimulation. Second, we found that protein degradation follows saturating kinetics, explaining the long half-life of Ste12 in mutants expressing elevated amounts of Ste12. Finally, our model reveals a novel mechanism for robust perfect adaptation through protein-protein interactions that enhance complex stability. This mechanism allows the transcriptional response to act on a shorter time scale than upstream pathway activity.

The regulation of filamentous growth in yeast

Filamentous growth is a nutrient-regulated growth response that occurs in many fungal species. In pathogens, filamentous growth is critical for host-cell attachment, invasion into tissues, and virulence. The budding yeast Saccharomyces cerevisiae undergoes filamentous growth, which provides a genetically tractable system to study the molecular basis of the response. Filamentous growth is regulated by evolutionarily conserved signaling pathways. One of these pathways is a mitogen activated protein kinase (MAPK) pathway. A remarkable feature of the filamentous growth MAPK pathway is that it is composed of factors that also function in other pathways. An intriguing challenge therefore has been to understand how pathways that share components establish and maintain their identity. Other canonical signaling pathways-rat sarcoma/protein kinase A (RAS/PKA), sucrose nonfermentable (SNF), and target of rapamycin (TOR)-also regulate filamentous growth, which raises the question of how signals from multiple pathways become integrated into a coordinated response. Together, these pathways regulate cell differentiation to the filamentous type, which is characterized by changes in cell adhesion, cell polarity, and cell shape. How these changes are accomplished is also discussed. High-throughput genomics approaches have recently uncovered new connections to filamentous growth regulation. These connections suggest that filamentous growth is a more complex and globally regulated behavior than is currently appreciated, which may help to pave the way for future investigations into this eukaryotic cell differentiation behavior.

A framework for mapping, visualisation and automatic model creation of signal-transduction networks

An intuitive formalism for reconstructing cellular networks from empirical data is presented, and used to build a comprehensive yeast MAP kinase network. The accompanying rxncon software tool can convert networks to a range of standard graphical formats and mathematical models.

Network mapping at the granularity of empirical data that largely avoids combinatorial complexity

Automatic visualisation and model generation with the rxncon open source software tool

Visualisation in a range of formats, including all three SBGN formats, as well as contingency matrix or regulatory graph

Comprehensive and completely references map of the yeast MAP kinase network in the rxncon format

Intracellular signalling systems are highly complex. This complexity makes handling, analysis and visualisation of available knowledge a major challenge in current signalling research. Here, we present a novel framework for mapping signal-transduction networks that avoids the combinatorial explosion by breaking down the network in reaction and contingency information. It provides two new visualisation methods and automatic export to mathematical models. We use this framework to compile the presently most comprehensive map of the yeast MAP kinase network. Our method improves previous strategies by combining (I) more concise mapping adapted to empirical data, (II) individual referencing for each piece of information, (III) visualisation without simplifications or added uncertainty, (IV) automatic visualisation in multiple formats, (V) automatic export to mathematical models and (VI) compatibility with established formats. The framework is supported by an open source software tool that facilitates integration of the three levels of network analysis: definition, visualisation and mathematical modelling. The framework is species independent and we expect that it will have wider impact in signalling research on any system.

Curated collection of yeast transcription factor DNA binding specificity data reveals novel structural and gene regulatory insights

Background

Transcription factors (TFs) play a central role in regulating gene expression by interacting with cis-regulatory DNA elements associated with their target genes. Recent surveys have examined the DNA binding specificities of most Saccharomyces cerevisiae TFs, but a comprehensive evaluation of their data has been lacking.

Results

We analyzed in vitro and in vivo TF-DNA binding data reported in previous large-scale studies to generate a comprehensive, curated resource of DNA binding specificity data for all characterized S. cerevisiae TFs. Our collection comprises DNA binding site motifs and comprehensive in vitro DNA binding specificity data for all possible 8-bp sequences. Investigation of the DNA binding specificities within the basic leucine zipper (bZIP) and VHT1 regulator (VHR) TF families revealed unexpected plasticity in TF-DNA recognition: intriguingly, the VHR TFs, newly characterized by protein binding microarrays in this study, recognize bZIP-like DNA motifs, while the bZIP TF Hac1 recognizes a motif highly similar to the canonical E-box motif of basic helix-loop-helix (bHLH) TFs. We identified several TFs with distinct primary and secondary motifs, which might be associated with different regulatory functions. Finally, integrated analysis of in vivo TF binding data with protein binding microarray data lends further support for indirect DNA binding in vivo by sequence-specific TFs.

Conclusions

The comprehensive data in this curated collection allow for more accurate analyses of regulatory TF-DNA interactions, in-depth structural studies of TF-DNA specificity determinants, and future experimental investigations of the TFs' predicted target genes and regulatory roles.