Characterization of the DNA-binding domains from the yeast cell-cycle transcription factors Mbp1 and Swi4

FvMbp1-Swi6 complex regulates vegetative growth, stress tolerance, and virulence in Fusarium verticillioides

The mycotoxigenic fungus Fusarium verticillioides is a common pathogen of grain and medicine that contaminates the host with fumonisin B1 (FB1) mycotoxin, poses serious threats to human and animal health. Therefore, it is crucial to unravel the regulatory mechanisms of growth, and pathogenicity of F. verticillioides. Mbp1 is a component of the MluI cell cycle box binding factor complex and acts as an APSES-type transcription factor that regulates cell cycle progression. However, no information is available regarding its role in F. verticillioides. In this study, we demonstrate that FvMbp1 interacts with FvSwi6 that acts as the cell cycle transcription factor, to form the heteromeric transcription factor complexes in F. verticillioides. Our results show that ΔFvMbp1 and ΔFvSwi6 both cause a severe reduction of vegetative growth, conidiation, and increase tolerance to diverse environmental stresses. Moreover, ΔFvMbp1 and ΔFvSwi6 dramatically decrease the virulence of the pathogen on the stalk and ear of maize. Transcriptome profiling show that FvMbp1-Swi6 complex co-regulates the expression of genes associated with multiple stress responses. These results indicate the functional importance of the FvMbp1-Swi6 complex in the filamentous fungi F. verticillioides and reveal a potential target for the effective prevention and control of Fusarium diseases.

Control of Gene Expression via the Yeast CWI Pathway

Living cells exposed to stressful environmental situations can elicit cellular responses that guarantee maximal cell survival. Most of these responses are mediated by mitogen-activated protein kinase (MAPK) cascades, which are highly conserved from yeast to humans. Cell wall damage conditions in the yeast Saccharomyces cerevisiae elicit rescue mechanisms mainly associated with reprogramming specific transcriptional responses via the cell wall integrity (CWI) pathway. Regulation of gene expression by this pathway is coordinated by the MAPK Slt2/Mpk1, mainly via Rlm1 and, to a lesser extent, through SBF (Swi4/Swi6) transcription factors. In this review, we summarize the molecular mechanisms controlling gene expression upon cell wall stress and the role of chromatin structure in these processes. Some of these mechanisms are also discussed in the context of other stresses governed by different yeast MAPK pathways. Slt2 regulates both transcriptional initiation and elongation by interacting with chromatin at the promoter and coding regions of CWI-responsive genes but using different mechanisms for Rlm1- and SBF-dependent genes. Since MAPK pathways are very well conserved in eukaryotic cells and are essential for controlling cellular physiology, improving our knowledge regarding how they regulate gene expression could impact the future identification of novel targets for therapeutic intervention.

Regulation of Yeast-to-Hyphae Transition in Yarrowia lipolytica

Many yeasts undergo a morphological transition from yeast-to-hyphal growth in response to environmental conditions. We used forward and reverse genetic techniques to identify genes regulating this transition in Yarrowia lipolytica. We confirmed that the transcription factor Ylmsn2 is required for the transition to hyphal growth and found that signaling by the histidine kinases Ylchk1 and Ylnik1 as well as the MAP kinases of the HOG pathway (Ylssk2, Ylpbs2, and Ylhog1) regulates the transition to hyphal growth. These results suggest that Y. lipolytica transitions to hyphal growth in response to stress through multiple kinase pathways. Intriguingly, we found that a repetitive portion of the genome containing telomere-like and rDNA repeats may be involved in the transition to hyphal growth, suggesting a link between this region and the general stress response.

The yeast Yarrowia lipolytica undergoes a morphological transition from yeast-to-hyphal growth in response to environmental conditions. A forward genetic screen was used to identify mutants that reliably remain in the yeast phase, which were then assessed by whole-genome sequencing. All the smooth mutants identified, so named because of their colony morphology, exhibit independent loss of DNA at a repetitive locus made up of interspersed ribosomal DNA and short 10- to 40-mer telomere-like repeats. The loss of repetitive DNA is associated with downregulation of genes with stress response elements (5′-CCCCT-3′) and upregulation of genes with cell cycle box (5′-ACGCG-3′) motifs in their promoter region. The stress response element is bound by the transcription factor Msn2p in Saccharomyces cerevisiae. We confirmed that the Y. lipolyticamsn2 (Ylmsn2) ortholog is required for hyphal growth and found that overexpression of Ylmsn2 enables hyphal growth in smooth strains. The cell cycle box is bound by the Mbp1p/Swi6p complex in S. cerevisiae to regulate G₁-to-S phase progression. We found that overexpression of either the Ylmbp1 or Ylswi6 homologs decreased hyphal growth and that deletion of either Ylmbp1 or Ylswi6 promotes hyphal growth in smooth strains. A second forward genetic screen for reversion to hyphal growth was performed with the smooth-33 mutant to identify additional genetic factors regulating hyphal growth in Y. lipolytica. Thirteen of the mutants sequenced from this screen had coding mutations in five kinases, including the histidine kinases Ylchk1 and Ylnik1 and kinases of the high-osmolarity glycerol response (HOG) mitogen-activated protein (MAP) kinase cascade Ylssk2, Ylpbs2, and Ylhog1. Together, these results demonstrate that Y. lipolytica transitions to hyphal growth in response to stress through multiple signaling pathways.

IMPORTANCE Many yeasts undergo a morphological transition from yeast-to-hyphal growth in response to environmental conditions. We used forward and reverse genetic techniques to identify genes regulating this transition in Yarrowia lipolytica. We confirmed that the transcription factor Ylmsn2 is required for the transition to hyphal growth and found that signaling by the histidine kinases Ylchk1 and Ylnik1 as well as the MAP kinases of the HOG pathway (Ylssk2, Ylpbs2, and Ylhog1) regulates the transition to hyphal growth. These results suggest that Y. lipolytica transitions to hyphal growth in response to stress through multiple kinase pathways. Intriguingly, we found that a repetitive portion of the genome containing telomere-like and rDNA repeats may be involved in the transition to hyphal growth, suggesting a link between this region and the general stress response.

The Inner Nuclear Membrane Protein Bqt4 in Fission Yeast Contains a DNA-Binding Domain Essential for Telomere Association with the Nuclear Envelope

Telomeres, the protective caps at the end of the chromosomes, are often associated with the nuclear envelope (NE). Telomere positioning to the NE is dynamically regulated during mitosis and meiosis. One inner nuclear membrane protein, Bqt4, in Schizosaccharomyces pombe plays essential roles in connecting telomeres to the NE. However, the structural basis of Bqt4 in mediating telomere-NE association is not clear. Here, we report the crystal structure of the N-terminal domain of Bqt4. The N-terminal domain of Bqt4 structurally resembles the APSES-family DNA-binding domain and has a moderate double-stranded DNA-binding activity. Disruption of Bqt4-DNA interaction results in telomere detachment from the NE. These data suggest that the DNA-binding activity of Bqt4 may function to prime the chromosome onto the NE and promote telomere-NE association.

Using single-index ODEs to study dynamic gene regulatory network

With the development of biotechnology, high-throughput studies on protein-protein, protein-gene, and gene-gene interactions become possible and attract remarkable attention. To explore the interactions in dynamic gene regulatory networks, we propose a single-index ordinary differential equation (ODE) model and develop a variable selection procedure. We employ the smoothly clipped absolute deviation penalty (SCAD) penalized function for variable selection. We analyze a yeast cell cycle gene expression data set to illustrate the usefulness of the single-index ODE model. In real data analysis, we group genes into functional modules using the smoothing spline clustering approach. We estimate state functions and their first derivatives for functional modules using penalized spline-based nonparametric mixed-effects models and the spline method. We substitute the estimates into the single-index ODE models, and then use the penalized profile least-squares procedure to identify network structures among the models. The results indicate that our model fits the data better than linear ODE models and our variable selection procedure identifies the interactions that may be missed by linear ODE models but confirmed in biological studies. In addition, Monte Carlo simulation studies are used to evaluate and compare the methods.

The evolution of a G1/S transcriptional network in yeasts

The G1-to-S cell cycle transition is promoted by the periodic expression of a large set of genes. In Saccharomyces cerevisiae G1/S gene expression is regulated by two transcription factor (TF) complexes, the MBF and SBF, which bind to specific DNA sequences, the MCB and SCB, respectively. Despite extensive research little is known regarding the evolution of the G1/S transcription regulation including the co-evolution of the DNA binding domains with their respective DNA binding sequences. We have recently examined the co-evolution of the G1/S TF specificity through the systematic generation and examination of chimeric Mbp1/Swi4 TFs containing different orthologue DNA binding domains in S. cerevisiae (Hendler et al. in PLoS Genet 13:e1006778. doi: 10.1371/journal.pgen.1006778 , 2017). Here, we review the co-evolution of G1/S transcriptional network and discuss the evolutionary dynamics and specificity of the MBF-MCB and SBF-SCB interactions in different fungal species.

Gene duplication and co-evolution of G1/S transcription factor specificity in fungi are essential for optimizing cell fitness

Transcriptional regulatory networks play a central role in optimizing cell survival. How DNA binding domains and cis-regulatory DNA binding sequences have co-evolved to allow the expansion of transcriptional networks and how this contributes to cellular fitness remains unclear. Here we experimentally explore how the complex G1/S transcriptional network evolved in the budding yeast Saccharomyces cerevisiae by examining different chimeric transcription factor (TF) complexes. Over 200 G1/S genes are regulated by either one of the two TF complexes, SBF and MBF, which bind to specific DNA binding sequences, SCB and MCB, respectively. The difference in size and complexity of the G1/S transcriptional network across yeast species makes it well suited to investigate how TF paralogs (SBF and MBF) and DNA binding sequences (SCB and MCB) co-evolved after gene duplication to rewire and expand the network of G1/S target genes. Our data suggests that whilst SBF is the likely ancestral regulatory complex, the ancestral DNA binding element is more MCB-like. G1/S network expansion took place by both cis- and trans- co-evolutionary changes in closely related but distinct regulatory sequences. Replacement of the endogenous SBF DNA-binding domain (DBD) with that from more distantly related fungi leads to a contraction of the SBF-regulated G1/S network in budding yeast, which also correlates with increased defects in cell growth, cell size, and proliferation.

Structural basis of DNA recognition by PCG2 reveals a novel DNA binding mode for winged helix-turn-helix domains

The MBP1 family proteins are the DNA binding subunits of MBF cell-cycle transcription factor complexes and contain an N terminal winged helix-turn-helix (wHTH) DNA binding domain (DBD). Although the DNA binding mechanism of MBP1 from Saccharomyces cerevisiae has been extensively studied, the structural framework and the DNA binding mode of other MBP1 family proteins remains to be disclosed. Here, we determined the crystal structure of the DBD of PCG2, the Magnaporthe oryzae orthologue of MBP1, bound to MCB-DNA. The structure revealed that the wing, the 20-loop, helix A and helix B in PCG2-DBD are important elements for DNA binding. Unlike previously characterized wHTH proteins, PCG2-DBD utilizes the wing and helix-B to bind the minor groove and the major groove of the MCB-DNA whilst the 20-loop and helix A interact non-specifically with DNA. Notably, two glutamines Q89 and Q82 within the wing were found to recognize the MCB core CGCG sequence through making hydrogen bond interactions. Further in vitro assays confirmed essential roles of Q89 and Q82 in the DNA binding. These data together indicate that the MBP1 homologue PCG2 employs an unusual mode of binding to target DNA and demonstrate the versatility of wHTH domains.

Gene expression analyses to explore the biomarkers and therapeutic targets for gliomas

To improve treatment strategies of glioma, microarray data were applied to screen target molecules that were regulated by microRNAs (miRNAs). GSE31262 was downloaded from Gene Expression Omnibus, including five neural stem cells samples from normal human and nine stem cells samples from glioma patients. Differentially expressed genes (DEGs) were identified with Multtest package and Limma package of R language, and false discovery rate < 0.05 and |log2FC (fold change)| >1 were chosen as cut-off criterion. Hierarchical clustering and pathway enrichment analysis of DEGs were performed using pheatmap package of R language and KOBAS software, respectively. miRNAs related to up- and down-regulated DEGs were, respectively, predicted by WebGestalt software and its miRNAs-target DEGs interaction network were, respectively, constructed by STRING database. Bingo plug-in in Cytoscape software was applied to analyze Gene Ontology functional enrichment analysis for up- and down-regulated DEGs in network, respectively. A total of 428 DEGs were selected, including 331 down-regulated and 97 up-regulated DEGs. Hierarchical clustering analysis showed that glioma samples and normal samples were completely separated. Pathway analysis indicated that CDK2 and WEE1 participated in the cell cycle. miR-124 could simultaneously regulate up-regulated (ELAVL1 and EZH2) and down-regulated (BACE1) DEGs. The down-regulated genes (KIF23, WEE1 and CDK2) were associated with cell division, while the up-regulated genes (PLP1 and MBP) were related to myelination of neurons. miR-124 might participate in development of glioma by regulating BACE1, ELAVL1 and EZH2. The biomarkers (KIF23, WEE1, CDK2, PLP1 and MBP) were considered as therapeutic targets of glioma.

Topology and control of the cell-cycle-regulated transcriptional circuitry

Nearly 20% of the budding yeast genome is transcribed periodically during the cell division cycle. The precise temporal execution of this large transcriptional program is controlled by a large interacting network of transcriptional regulators, kinases, and ubiquitin ligases. Historically, this network has been viewed as a collection of four coregulated gene clusters that are associated with each phase of the cell cycle. Although the broad outlines of these gene clusters were described nearly 20 years ago, new technologies have enabled major advances in our understanding of the genes comprising those clusters, their regulation, and the complex regulatory interplay between clusters. More recently, advances are being made in understanding the roles of chromatin in the control of the transcriptional program. We are also beginning to discover important regulatory interactions between the cell-cycle transcriptional program and other cell-cycle regulatory mechanisms such as checkpoints and metabolic networks. Here we review recent advances and contemporary models of the transcriptional network and consider these models in the context of eukaryotic cell-cycle controls.

A unique spumavirus Gag N-terminal domain with functional properties of orthoretroviral matrix and capsid

The Spumaretrovirinae, or foamyviruses (FVs) are complex retroviruses that infect many species of monkey and ape. Although FV infection is apparently benign, trans-species zoonosis is commonplace and has resulted in the isolation of the Prototypic Foamy Virus (PFV) from human sources and the potential for germ-line transmission. Despite little sequence homology, FV and orthoretroviral Gag proteins perform equivalent functions, including genome packaging, virion assembly, trafficking and membrane targeting. In addition, PFV Gag interacts with the FV Envelope (Env) protein to facilitate budding of infectious particles. Presently, there is a paucity of structural information with regards FVs and it is unclear how disparate FV and orthoretroviral Gag molecules share the same function. Therefore, in order to probe the functional overlap of FV and orthoretroviral Gag and learn more about FV egress and replication we have undertaken a structural, biophysical and virological study of PFV-Gag. We present the crystal structure of a dimeric amino terminal domain from PFV, Gag-NtD, both free and in complex with the leader peptide of PFV Env. The structure comprises a head domain together with a coiled coil that forms the dimer interface and despite the shared function it is entirely unrelated to either the capsid or matrix of Gag from other retroviruses. Furthermore, we present structural, biochemical and virological data that reveal the molecular details of the essential Gag-Env interaction and in addition we also examine the specificity of Trim5α restriction of PFV. These data provide the first information with regards to FV structural proteins and suggest a model for convergent evolution of gag genes where structurally unrelated molecules have become functionally equivalent.

Structural models of protein-DNA complexes based on interface prediction and docking

Protein-DNA interactions are the physical basis of gene expression and DNA modification. Structural models that reveal these interactions are essential for their understanding. As only a limited number of structures for protein-DNA complexes have been determined by experimental methods, computation methods provide a potential way to fill the need. We have developed the DISPLAR method to predict DNA binding sites on proteins. Predicted binding sites have been used to assist the building of structural models by docking, either by guiding the docking or by selecting near-native candidates from the docked poses. Here we applied the DISPLAR method to predict the DNA binding sites for 20 DNA-binding proteins, which have had their DNA binding sites characterized by NMR chemical shift perturbation. For two of these proteins, the structures of their complexes with DNA have also been determined. With the help of the DISPLAR predictions, we built structural models for these two complexes. Evaluations of both the DNA binding sites for 20 proteins and the structural models of the two protein-DNA complexes against experimental results demonstrate the significant promise of our model-building approach.

Regulation of cell wall biogenesis in Saccharomyces cerevisiae: the cell wall integrity signaling pathway

The yeast cell wall is a strong, but elastic, structure that is essential not only for the maintenance of cell shape and integrity, but also for progression through the cell cycle. During growth and morphogenesis, and in response to environmental challenges, the cell wall is remodeled in a highly regulated and polarized manner, a process that is principally under the control of the cell wall integrity (CWI) signaling pathway. This pathway transmits wall stress signals from the cell surface to the Rho1 GTPase, which mobilizes a physiologic response through a variety of effectors. Activation of CWI signaling regulates the production of various carbohydrate polymers of the cell wall, as well as their polarized delivery to the site of cell wall remodeling. This review article centers on CWI signaling in Saccharomyces cerevisiae through the cell cycle and in response to cell wall stress. The interface of this signaling pathway with other pathways that contribute to the maintenance of cell wall integrity is also discussed.

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.

Extensive characterization of NF-κB binding uncovers non-canonical motifs and advances the interpretation of genetic functional traits

BACKGROUND

Genetic studies have provided ample evidence of the influence of non-coding DNA polymorphisms on trait variance, particularly those occurring within transcription factor binding sites. Protein binding microarrays and other platforms that can map these sites with great precision have enhanced our understanding of how a single nucleotide polymorphism can alter binding potential within an in vitro setting, allowing for greater predictive capability of its effect on a transcription factor binding site.

RESULTS

We have used protein binding microarrays and electrophoretic mobility shift assay-sequencing (EMSA-Seq), a deep sequencing based method we developed to analyze nine distinct human NF-κB dimers. This family of transcription factors is one of the most extensively studied, but our understanding of its DNA binding preferences has been limited to the originally described consensus motif, GGRRNNYYCC. We highlight differences between NF-κB family members and also put under the spotlight non-canonical motifs that have so far received little attention. We utilize our data to interpret the binding of transcription factors between individuals across 1,405 genomic regions laden with single nucleotide polymorphisms. We also associated binding correlations made using our data with risk alleles of disease and demonstrate its utility as a tool for functional studies of single nucleotide polymorphisms in regulatory regions.

CONCLUSIONS

NF-κB dimers bind specifically to non-canonical motifs and these can be found within genomic regions in which a canonical motif is not evident. Binding affinity data generated with these different motifs can be used in conjunction with data from chromatin immunoprecipitation-sequencing (ChIP-Seq) to enable allele-specific analyses of expression and transcription factor-DNA interactions on a genome-wide scale.

Spt10 and Swi4 control the timing of histone H2A/H2B gene activation in budding yeast

The expression of the histone genes is regulated during the cell cycle to provide histones for nucleosome assembly during DNA replication. In budding yeast, histones H2A and H2B are expressed from divergent promoters at the HTA1-HTB1 and HTA2-HTB2 loci. Here, we show that the major activator of HTA1-HTB1 is Spt10, a sequence-specific DNA binding protein with a putative histone acetyltransferase (HAT) domain. Spt10 binds to two pairs of upstream activation sequence (UAS) elements in the HTA1-HTB1 promoter: UAS1 and UAS2 drive HTA1 expression, and UAS3 and UAS4 drive HTB1 expression. UAS3 and UAS4 also contain binding sites for the cell cycle regulator SBF (an Swi4-Swi6 heterodimer), which overlap the Spt10 binding sites. The binding of Spt10 and binding of SBF to UAS3 and UAS4 are mutually exclusive in vitro. Both SBF and Spt10 are bound in cells arrested with α-factor, apparently awaiting a signal to activate transcription. Soon after the removal of α-factor, SBF initiates a small, early peak of HTA1 and HTB1 transcription, which is followed by a much larger peak due to Spt10. Both activators dissociate from the HTA1-HTB1 promoter after expression has been activated. Thus, SBF and Spt10 cooperate to control the timing of HTA1-HTB1 expression.

Transcriptional reporters for genes activated by cell wall stress through a non-catalytic mechanism involving Mpk1 and SBF

The Mpk1 MAP kinase of the cell wall integrity (CWI) signalling pathway induces transcription of the FKS2 gene in response to cell wall stress through a non-catalytic mechanism that involves stable association of Mpk1 with the Swi4 transcription factor. This dimeric complex binds to a Swi4 recognition site in the FKS2 promoter. The Swi6 transcription factor is also required to bind this ternary complex for transcription initiation to ensue. In this context, the Mlp1 pseudokinase serves a redundant function with Mpk1. We have identified three additional genes, CHA1, YLR042c and YKR013w, that are induced by cell wall stress through the same mechanism. We report on the behaviour of several promoter-lacZ reporter plasmids designed to detect cell wall stress transcription through this pathway.

Identifying functional relationships within sets of co-expressed genes by combining upstream regulatory motif analysis and gene expression information

BACKGROUND

Existing clustering approaches for microarray data do not adequately differentiate between subsets of co-expressed genes. We devised a novel approach that integrates expression and sequence data in order to generate functionally coherent and biologically meaningful subclusters of genes. Specifically, the approach clusters co-expressed genes on the basis of similar content and distributions of predicted statistically significant sequence motifs in their upstream regions.

RESULTS

We applied our method to several sets of co-expressed genes and were able to define subsets with enrichment in particular biological processes and specific upstream regulatory motifs.

CONCLUSIONS

These results show the potential of our technique for functional prediction and regulatory motif identification from microarray data.

MAPK cell-cycle regulation in Saccharomyces cerevisiae and Candida albicans

The cell cycle is the sequential set of events that living cells undergo in order to duplicate. This process must be tightly regulated as alterations may lead to diseases such as cancer. The molecular events that control the cell cycle are directional and involve regulatory molecules such as cyclins and cyclin-dependent kinases (CDKs). The budding yeast Saccharomyces cerevisiae has become a model to study this complex system since it shares several mechanisms with higher eukaryotes. Signal transduction pathways are biochemical mechanisms that sense environmental changes and there is recent evidence that they control the progression through the cell cycle in response to several stimuli. In response to pheromone, the budding yeast arrests the cell cycle in the G1 phase at the START stage. Activation of the pheromone response pathway leads to the phosphorylation of Far1, which inhibits the function of complexes formed by G1 cyclins (Cln1 and Cln2) and the CDK (Cdc28), blocking the transition to the S phase. This response prepares the cells to fuse cytoplasms and nuclei to generate a diploid cell. Activation of the Hog1 MAP kinase in response to osmotic stress or arsenite leads to the transient arrest of the cell cycle in G1 phase, which is mediated by direct phosphorylation of the CDK inhibitor, Sic1, and by downregulation of cyclin expression. Osmotic stress also induces a delay in G2 phase by direct phosphorylation of Hsl7 via Hog1, which results in the accumulation of Swe1. As a consequence, cell cycle arrest allows cells to survive upon stress. Finally, cell wall damage can induce cell cycle arrest at G2 via the cell integrity MAPK Slt2. By linking MAPK signal transduction pathways to the cell cycle machinery, a tight and precise control of the cell division takes place in response to environmental changes. Research into similar MAPK-mediated cell cycle regulation in the opportunistic pathogen Candida albicans may result in the development of new antifungal therapies.

The transcriptional network activated by Cln3 cyclin at the G1-to-S transition of the yeast cell cycle

Background

The G1-to-S transition of the cell cycle in the yeast Saccharomyces cerevisiae involves an extensive transcriptional program driven by transcription factors SBF (Swi4-Swi6) and MBF (Mbp1-Swi6). Activation of these factors ultimately depends on the G1 cyclin Cln3.

Results

To determine the transcriptional targets of Cln3 and their dependence on SBF or MBF, we first have used DNA microarrays to interrogate gene expression upon Cln3 overexpression in synchronized cultures of strains lacking components of SBF and/or MBF. Secondly, we have integrated this expression dataset together with other heterogeneous data sources into a single probabilistic model based on Bayesian statistics. Our analysis has produced more than 200 transcription factor-target assignments, validated by ChIP assays and by functional enrichment. Our predictions show higher internal coherence and predictive power than previous classifications. Our results support a model whereby SBF and MBF may be differentially activated by Cln3.

Conclusions

Integration of heterogeneous genome-wide datasets is key to building accurate transcriptional networks. By such integration, we provide here a reliable transcriptional network at the G1-to-S transition in the budding yeast cell cycle. Our results suggest that to improve the reliability of predictions we need to feed our models with more informative experimental data.

Yeast Mpk1 cell wall integrity mitogen-activated protein kinase regulates nucleocytoplasmic shuttling of the Swi6 transcriptional regulator

The yeast SBF transcription factor is a heterodimer comprised of Swi4 and Swi6 that has a well defined role in cell cycle-specific transcription. SBF serves a second function in the transcriptional response to cell wall stress in which activated Mpk1 mitogen-activated protein kinase of the cell wall integrity signaling pathway forms a complex with Swi4, the DNA binding subunit of SBF, conferring upon Swi4 the ability to bind DNA and activate transcription of FKS2. Although Mpk1-Swi4 complex formation and transcriptional activation of FKS2 does not require Mpk1 catalytic activity, Swi6 is phosphorylated by Mpk1 and must be present in the Mpk1-Swi4 complex for transcriptional activation of FKS2. Here, we find that Mpk1 regulates Swi6 nucleocytoplasmic shuttling in a biphasic manner. First, formation of the Mpk1-Swi4 complex recruits Swi6 to the nucleus for transcriptional activation. Second, Mpk1 negatively regulates Swi6 by phosphorylation on Ser238, which inhibits nuclear entry. Ser238 neighbors a nuclear localization signal (NLS) whose function is blocked by phosphorylation at Ser238 in a manner similar to the regulation by Cdc28 of another Swi6 NLS, revealing a mechanism for the integration of multiple signals to a single endpoint. Finally, the Kap120 beta-importin binds the Mpk1-regulated Swi6 NLS but not the Cdc28-regulated NLS.

Structural analysis of the DNA target site and its interaction with Mbp1

The solution structure of a 14 base-pair non-self complementary DNA duplex containing the consensus-binding site of the yeast transcription factor Mbp1 has been determined by NMR using a combination of scalar coupling analysis, time-dependent NOEs, residual dipolar couplings and 13C-edited NMR spectroscopy of a duplex prepared with one strand uniformly labeled with 13C-nucleotides. As expected, the free DNA duplex is within the B-family of structures, and within experimental limits is straight. However, there are clear local structural variations associated with the consensus CGCG element in the binding sequence that are important for sequence recognition. In the complex, the DNA bends around the protein, which also undergoes some conformational rearrangement in the C-terminal region. Structural constraints derived from paramagnetic perturbation experiments with spin-labeled DNA, chemical shift perturbation experiments of the DNA, previous cross-saturation, chemical shift perturbation experiments on the protein, information from mutational analysis, and electrostatics calculations have been used to produce a detailed docked structure using the known solution conformation of the free protein and other spectroscopic information about the Mbp1:DNA complex. A Monte Carlo-based docking procedure with restrained MD in a fully solvated system subjected to available experimental constraints produced models that account for the available structural data, and can rationalize the extensive thermodynamic data about the Mbp1:DNA complex. The protein:DNA interface is closely packed and is associated with a small number of specific contacts. The structure shows an extensive positively charged surface that accounts for the high polyelectrolyte contribution to binding.

Mechanism of Mpk1 mitogen-activated protein kinase binding to the Swi4 transcription factor and its regulation by a novel caffeine-induced phosphorylation

The Mpk1 mitogen-activated protein kinase (MAPK) of the cell wall integrity signaling pathway uses a noncatalytic mechanism to activate the SBF (Swi4/Swi6) transcription factor. Active Mpk1 forms a complex with Swi4, the DNA-binding subunit of SBF, conferring the ability to bind DNA. Because SBF activation is independent of Mpk1 catalytic activity but requires Mpk1 to be in an active conformation, we sought to understand how Mpk1 interacts with Swi4. Mutational analysis revealed that binding and activation of Swi4 by Mpk1 requires an intact D-motif-binding site, a docking surface common to MAPKs that resides distal to the phosphorylation loop but does not require the substrate-binding site, revealing a novel mechanism for MAPK target regulation. Additionally, we found that Mpk1 binds near the autoinhibitory C terminus of Swi4, suggesting an activation mechanism in which Mpk1 substitutes for Swi6 in promoting Swi4 DNA binding. Finally, we show that caffeine is an atypical activator of cell wall integrity signaling, because it induces phosphorylation of the Mpk1 C-terminal extension at Ser423 and Ser428. These phosphorylations were dependent on the DNA damage checkpoint kinases, Mec1/Tel1 and Rad53. Phosphorylation of Ser423 specifically blocked SBF activation by preventing Mpk1 association with Swi4, revealing a novel mechanism for regulating MAPK target specificity.

Distinguishing direct versus indirect transcription factor-DNA interactions

Transcriptional regulation is largely enacted by transcription factors (TFs) binding DNA. Large numbers of TF binding motifs have been revealed by ChIP-chip experiments followed by computational DNA motif discovery. However, the success of motif discovery algorithms has been limited when applied to sequences bound in vivo (such as those identified by ChIP-chip) because the observed TF-DNA interactions are not necessarily direct: Some TFs predominantly associate with DNA indirectly through protein partners, while others exhibit both direct and indirect binding. Here, we present the first method for distinguishing between direct and indirect TF-DNA interactions, integrating in vivo TF binding data, in vivo nucleosome occupancy data, and motifs from in vitro protein binding microarray experiments. When applied to yeast ChIP-chip data, our method reveals that only 48% of the data sets can be readily explained by direct binding of the profiled TF, while 16% can be explained by indirect DNA binding. In the remaining 36%, none of the motifs used in our analysis was able to explain the ChIP-chip data, either because the data were too noisy or because the set of motifs was incomplete. As more in vitro TF DNA binding motifs become available, our method could be used to build a complete catalog of direct and indirect TF-DNA interactions. Our method is not restricted to yeast or to ChIP-chip data, but can be applied in any system for which both in vivo binding data and in vitro DNA binding motifs are available.

Yeast Mpk1 mitogen-activated protein kinase activates transcription through Swi4/Swi6 by a noncatalytic mechanism that requires upstream signal

The cell wall integrity mitogen-activated protein kinase (MAPK) cascade of Saccharomyces cerevisiae drives changes in gene expression in response to cell wall stress. We show that the MAPK of this pathway (Mpk1) and its pseudokinase paralog (Mlp1) use a noncatalytic mechanism to activate transcription of the FKS2 gene. Transcriptional activation of FKS2 was dependent on the Swi4/Swi6 (SBF) transcription factor and on an activating signal to Mpk1 but not on protein kinase activity. Activated (phosphorylated) Mpk1 and Mlp1 were detected in a complex with Swi4 and Swi6 at the FKS2 promoter. Mpk1 association with Swi4 in vivo required phosphorylation of Mpk1. Promoter association of Mpk1 and the Swi4 DNA-binding subunit of SBF were codependent but did not require Swi6, indicating that the MAPK confers DNA-binding ability to Swi4. Based on these data, we propose a model in which phosphorylated Mpk1 or Mlp1 forms a dimeric complex with Swi4 that is competent to associate with the FKS2 promoter. This complex then recruits Swi6 to activate transcription. Finally, we show that human ERK5, a functional ortholog of Mpk1, is similarly capable of driving FKS2 expression in the absence of protein kinase activity, suggesting that this mammalian MAPK may also have a noncatalytic function in vivo.

Unique 31P spectral response to the formation of a specific restriction enzyme-DNA complex

Protein-induced distortion is a dramatic but not universally observed feature of sequence-specific DNA interactions. This is illustrated by the crystal structures of restriction enzyme-DNA complexes: While some of these structures exhibit DNA distortion, others do not. Among the latter is Pvull endonuclease, a small enzyme that is also amenable to NMR spectroscopic studies. Here 31P NMR spectroscopy is applied to demonstrate the unique spectral response of DNA to sequence-specific protein interactions. The 31P NMR spectrum of a noncognate DNA exhibits only spectral broadening upon the addition of enzyme. However, when enzyme is added to target DNA, a number of 31P resonances shift dramatically. The magnitudes of the chemical shifts (2-3 ppm) are among the largest observed. Site-specific substitution with phosphoramidates and phosphorothioates are used analyze these effects. While such spectral features have been interpreted as indicative of DNA backbone distortions, FRET analysis indicates that this does not occur in PvuII-cognate DNA complexes in solution. The distinct 31P spectral signature observed for cognate DNA mirrors that observed for the enzyme, underscoring the unique features of cognate complex formation.

Mapping and conformational characterization of the DNA-binding region of the breast cancer susceptibility protein BRCA1

The breast cancer susceptibility gene, BRCA1, encodes a large nuclear phosphoprotein, the major isoform of which is 1863 amino acids in size. Structure-function studies have been largely restricted to the only two domains identified by homology searches: the RING (really interesting new gene) and BRCT (BRCA1 C-terminus) domains. However, we have recently reported the identification of a large central soluble region of BRCA1 (residues 230-534) that binds specifically to four-way junction DNA, a property that potentially facilitates its role in the repair of DNA lesions by homologous recombination. We have now used a combination of limited proteolysis and extension cloning to identify more accurately the DNA-binding region of BRCA1. Limited trypsinolysis of BRCA1-(230-534) resulted in the production of a soluble domain identified as residues 230-339. However, after cloning, expression and purification of this region, studies revealed that it was unable to bind to four-way junctions, suggesting that the DNA-binding activity, in part, resides within residues 340-534. A series of fragments extending from residue 340 were produced, and each was tested for its ability to bind to four-way junction DNA in gel retardation assays. In these experiments, residues 340-554 of BRCA1 were identified as the minimal DNA-binding region. We then went on to characterize the conformation of this region using CD spectroscopy and analytical centrifugation.

The cell cycle-regulated genes of Schizosaccharomyces pombe

Many genes are regulated as an innate part of the eukaryotic cell cycle, and a complex transcriptional network helps enable the cyclic behavior of dividing cells. This transcriptional network has been studied in Saccharomyces cerevisiae (budding yeast) and elsewhere. To provide more perspective on these regulatory mechanisms, we have used microarrays to measure gene expression through the cell cycle of Schizosaccharomyces pombe (fission yeast). The 750 genes with the most significant oscillations were identified and analyzed. There were two broad waves of cell cycle transcription, one in early/mid G2 phase, and the other near the G2/M transition. The early/mid G2 wave included many genes involved in ribosome biogenesis, possibly explaining the cell cycle oscillation in protein synthesis in S. pombe. The G2/M wave included at least three distinctly regulated clusters of genes: one large cluster including mitosis, mitotic exit, and cell separation functions, one small cluster dedicated to DNA replication, and another small cluster dedicated to cytokinesis and division. S. pombe cell cycle genes have relatively long, complex promoters containing groups of multiple DNA sequence motifs, often of two, three, or more different kinds. Many of the genes, transcription factors, and regulatory mechanisms are conserved between S. pombe and S. cerevisiae. Finally, we found preliminary evidence for a nearly genome-wide oscillation in gene expression: 2,000 or more genes undergo slight oscillations in expression as a function of the cell cycle, although whether this is adaptive, or incidental to other events in the cell, such as chromatin condensation, we do not know.

A comprehensive examination of gene expression throughout the cell cycle of fission yeast is compared with recent related studies to highlight robust transcriptional patterns.

High functional overlap between MluI cell-cycle box binding factor and Swi4/6 cell-cycle box binding factor in the G1/S transcriptional program in Saccharomyces cerevisiae

In budding yeast, many genes are induced early in the cell cycle. Induction of these genes has been predominantly attributed to two transcription factors, Swi4-Swi6 (SBF) and Mbp1-Swi6 (MBF). Swi4 and Mbp1 are related DNA-binding proteins with dissimilar target sequences. For most G1/S-regulated genes that we tested in a cdc20 block-release protocol for cell-cycle synchronization, removal of both Swi4 and Mbp1 was necessary and sufficient to essentially eliminate cell-cycle-regulated expression. Detectable SBF or MBF binding sites (SCBs or MCBs) in the promoters or available genome-wide promoter occupancy data do not consistently explain this functional overlap. The overlapping ability of these transcription factors to regulate many promoters with very similar cell-cycle kinetics may provide robustness to the G1/S transcriptional response, but poses a puzzle with respect to promoter-transcription factor specificity. In addition, for some genes, deletion of Mbp1 or Swi4 enhances transcription, suggesting that these factors can also function as transcriptional repressors. Finally, we observe residual G1/S transcriptional regulation in the absence of Swi4 and Mbp1.

The yeast pafl-rNA polymerase II complex is required for full expression of a subset of cell cycle-regulated genes

We have previously described an alternative form of RNA polymerase II in yeast lacking the Srb and Med proteins but including Pafl, Cdc73, Hprl, and Ccr4. The Pafl-RNA polymerase II complex (Paf1 complex) acts in the same pathway as the Pkc1-mitogen-activated protein kinase cascade and is required for full expression of many cell wall biosynthetic genes. The expression of several of these cell integrity genes, as well as many other Paf1-requiring genes identified by differential display and microarray analyses, is regulated during the cell cycle. To determine whether the Paf1 complex is required for basal or cyclic expression of these genes, we assayed transcript abundance throughout the cell cycle. We found that transcript abundance for a subset of cell cycle-regulated genes, including CLN1, HO, RNR1, and FAR1, is reduced from 2- to 13-fold in a paf1delta strain, but that this reduction is not promoter dependent. Despite the decreased expression levels, cyclic expression is still observed. We also examined the possibility that the Paf1 complex acts in the same pathway as either SBF (Swi4/Swi6) or MBF (Mbp1/Swi6), the partially redundant cell cycle transcription factors. Consistent with the possibility that they have overlapping essential functions, we found that loss of Paf1 is lethal in combination with loss of Swi4 or Swi6. In addition, overexpression of either Swi4 or Mbp1 suppresses some paf1delta phenotypes. These data establish that the Paf1 complex plays an important role in the essential regulatory pathway controlled by SBF and MBF.

Structure of the retinal determination protein Dachshund reveals a DNA binding motif

The Dachshund proteins are essential components of a regulatory network controlling cell fate determination. They have been implicated in eye, limb, brain, and muscle development. These proteins cannot be assigned to any recognizable structural or functional class based on amino acid sequence analysis. The 1.65 A crystal structure of the most conserved domain of human DACHSHUND is reported here. The protein forms an alpha/beta structure containing a DNA binding motif similar to that found in the winged helix/forkhead subgroup of the helix-turn-helix family. This unexpected finding alters the previously proposed molecular models for the role of Dachshund in the eye determination pathway. Furthermore, it provides a rational framework for future mechanistic analyses of the Dachshund proteins in several developmental contexts.

Transcriptional coregulation by the cell integrity mitogen-activated protein kinase Slt2 and the cell cycle regulator Swi4

In Saccharomyces cerevisiae, the heterodimeric transcription factor SBF (for SCB binding factor) is composed of Swi4 and Swi6 and activates gene expression at the G(1)/S-phase transition of the mitotic cell cycle. Cell cycle commitment is associated not only with major alterations in gene expression but also with highly polarized cell growth; the mitogen-activated protein kinase (MAPK) Slt2 is required to maintain cell wall integrity during periods of polarized growth and cell wall stress. We describe experiments aimed at defining the regulatory pathway involving the cell cycle transcription factor SBF and Slt2-MAPK. Gene expression assays and chromatin immunoprecipitation experiments revealed Slt2-dependent recruitment of SBF to the promoters of the G(1) cyclins PCL1 and PCL2 after activation of the Slt2-MAPK pathway. We performed DNA microarray analysis and identified other genes whose expression was reduced in both SLT2 and SWI4 deletion strains. Genes that are sensitive to both Slt2 and Swi4 appear to be uniquely regulated and reveal a role for Swi4, the DNA-binding component of SBF, which is independent of the regulatory subunit Swi6. Some of the Swi4- and Slt2-dependent genes do not require Swi6 for either their expression or for Swi4 localization to their promoters. Consistent with these results, we found a direct interaction between Swi4 and Slt2. Our results establish a new Slt2-dependent mode of Swi4 regulation and suggest roles for Swi4 beyond its prominent role in controlling cell cycle transcription.

Assessing clusters and motifs from gene expression data

Large-scale gene expression studies and genomic sequencing projects are providing vast amounts of information that can be used to identify or predict cellular regulatory processes. Genes can be clustered on the basis of the similarity of their expression profiles or function and these clusters are likely to contain genes that are regulated by the same transcription factors. Searches for cis-regulatory elements can then be undertaken in the noncoding regions of the clustered genes. However, it is necessary to assess the efficiency of both the gene clustering and the postulated regulatory motifs, as there are many difficulties associated with clustering and determining the functional relevance of matches to sequence motifs. We have developed a method to assess the potential functional significance of clusters and motifs based on the probability of finding a certain number of matches to a motif in all of the gene clusters. To avoid problems with threshold scores for a match, the top matches to a motif are taken in several sample sizes. Genes from a sample are then counted by the cluster in which they appear. The probability of observing these counts by chance is calculated using the hypergeometric distribution. Because of the multiple sample sizes, strong and weak matching motifs can be detected and refined and significant matches to motifs across cluster boundaries are observed as all clusters are considered. By applying this method to many motifs and to a cluster set of yeast genes, we detected a similarity between Swi Five Factor and forkhead proteins and suggest that the currently unidentified Swi Five Factor is one of the yeast forkhead proteins.