Two yeast forkhead genes regulate the cell cycle and pseudohyphal growth

Knock-out of the major regulator Flo8 in Komagataella phaffii results in unique host strain performance for methanol-free recombinant protein production

Flo8 is a main transcriptional regulator of flocculation and pseudohyphal growth in yeast. Disruption of FLO8 in the popular recombinant protein production host Komagataella phaffii (Pichia pastoris) prevents pseudohyphal growth and reduces cell-to-surface adherence, making the mutant an interesting platform for research and industry. However, knowledge of the physiological impact of the mutation remained scarce. In-depth analysis of transcriptome data from FLO8-deficient K. phaffii revealed that Flo8 affects genes involved in cell cycle, mating, respiration, and catabolite repression additionally to flocculation targets. One gene with considerably increased expression in flo8 was GTH1, encoding a high-affinity glucose transporter in K. phaffii. Its promoter (PG1) was previously established as a strong, glucose-regulatable alternative to methanol-induced promoters. PG1 and its improved derivatives PG1-3, D-PGS4 and D-PGS5, proved to be promising candidates for controlling recombinant protein production in the FLO8-deficient background. In small-scale screenings, PG13-controlled intracellular EGFP levels were 2.8-fold higher, and yields of different secreted recombinant proteins were up to 4.8-fold increased. The enhanced productivity of the flo8 mutant in combination with the PG1 variants was transferrable to glucose-limited fed-batch processes and could largely be attributed to higher transcriptional activity of the promoter, leading to a much higher productivity per chromosomally integrated gene copy. K. phaffii flo8 has many advantageous characteristics, such as reduced surface growth and increased transcriptional strength of glucose-regulatable promoters. These features turn the flo8 strain into a valuable new base strain for various experimental designs and establish flo8 as an excellent strain background for methanol-free recombinant protein production processes.

A winged-helix DNA-binding protein is essential for self-fertility during sexual development of the homothallic fungus Fusarium graminearum

Sexual reproduction is crucial for increasing the genetic diversity of populations and providing overwintering structures, such as perithecia and associated tissue, in the destructive plant pathogenic fungus Fusarium graminearum. While mating-type genes serve as master regulators in fungal sexual reproduction, the molecular mechanisms underlying this process remain elusive. Winged-helix DNA-binding proteins are key regulators of embryogenesis and cell differentiation in higher eukaryotes. These proteins are implicated in the morphogenesis and development of several fungal species. However, their involvement in sexual reproduction remains largely unexplored in F. graminearum. Here, we investigated the function of winged-helix DNA-binding proteins in vegetative growth, conidiation, and sexual reproduction, with a specific focus on the FgWING27, which is highly conserved among Fusarium species. Deletion of FgWING27 resulted in an abnormal pattern characterized by a gradual increase in the expression of mating-type genes during sexual development, indicating its crucial role in the stage-specific genetic regulation of MAT genes in the late stages of sexual development. Furthermore, using chromatin immunoprecipitation followed by sequencing analysis, we identified Fg17056 as a downstream gene of Fgwing27, which is essential for sexual reproduction. These findings underscore the significance of winged-helix DNA-binding proteins in fungal development and reproduction in F. graminearum, and highlight the pivotal role of Fgwing27 as a core genetic factor in the intricate genetic regulatory network governing sexual reproduction.IMPORTANCEFusarium graminearum is a devastating plant pathogenic fungus causing significant economic losses due to reduced crop yields. In Fusarium Head Blight epidemics, spores produced through sexual and asexual reproduction serve as inoculum, making it essential to understand the fungal reproduction process. Here, we focus on winged-helix DNA-binding proteins, which have been reported to play crucial roles in cell cycle regulation and differentiation, and address their requirement in the sexual reproduction of F. graminearum. Furthermore, we identified a highly conserved protein in Fusarium as a key factor in self-fertility, along with the discovery of its direct downstream genes. This provides crucial information for constructing the complex genetic regulatory network of sexual reproduction and significantly contribute to further research on sexual reproduction in Fusarium species.

The budding yeast Fkh1 Forkhead associated (FHA) domain promotes a G1-chromatin state and the activity of chromosomal DNA replication origins

In Saccharomyces cerevisiae, the forkhead (Fkh) transcription factor Fkh1 (forkhead homolog) enhances the activity of many DNA replication origins that act in early S-phase (early origins). Current models posit that Fkh1 acts directly to promote these origins' activity by binding to origin-adjacent Fkh1 binding sites (FKH sites). However, the post-DNA binding functions that Fkh1 uses to promote early origin activity are poorly understood. Fkh1 contains a conserved FHA (forkhead associated) domain, a protein-binding module with specificity for phosphothreonine (pT)-containing partner proteins. At a small subset of yeast origins, the Fkh1-FHA domain enhances the ORC (origin recognition complex)-origin binding step, the G1-phase event that initiates the origin cycle. However, the importance of the Fkh1-FHA domain to either chromosomal replication or ORC-origin interactions at genome scale is unclear. Here, S-phase SortSeq experiments were used to compare genome replication in proliferating FKH1 and fkh1-R80A mutant cells. The Fkh1-FHA domain promoted the activity of ≈ 100 origins that act in early to mid- S-phase, including the majority of centromere-associated origins, while simultaneously inhibiting ≈ 100 late origins. Thus, in the absence of a functional Fkh1-FHA domain, the temporal landscape of the yeast genome was flattened. Origins are associated with a positioned nucleosome array that frames a nucleosome depleted region (NDR) over the origin, and ORC-origin binding is necessary but not sufficient for this chromatin organization. To ask whether the Fkh1-FHA domain had an impact on this chromatin architecture at origins, ORC ChIPSeq data generated from proliferating cells and MNaseSeq data generated from G1-arrested and proliferating cell populations were assessed. Origin groups that were differentially regulated by the Fkh1-FHA domain were characterized by distinct effects of this domain on ORC-origin binding and G1-phase chromatin. Thus, the Fkh1-FHA domain controlled the distinct chromatin architecture at early origins in G1-phase and regulated origin activity in S-phase.

Interactions between Fkh1 monomers stabilize its binding to DNA replication origins

Eukaryotic DNA replication is initiated from multiple genomic origins, which can be broadly categorized as firing early or late in the S phase. Several factors can influence the temporal usage of origins to determine the timing of their firing. In budding yeast, the Forkhead family proteins Fkh1 and Fkh2 bind to a subset of replication origins and activate them at the beginning of the S phase. In these origins, the Fkh1/2 binding sites are arranged in a strict configuration, suggesting that Forkhead factors must bind the origins in a specific manner. To explore these binding mechanisms in more detail, we mapped the domains of Fkh1 that were required for its role in DNA replication regulation. We found that a short region of Fkh1 near its DNA binding domain was essential for the protein to bind and activate replication origins. Analysis of purified Fkh1 proteins revealed that this region mediates dimerization of Fkh1, suggesting that intramolecular contacts of Fkh1 are required for efficient binding and regulation of DNA replication origins. We also show that the Sld3-Sld7-Cdc45 complex is recruited to Forkhead-regulated origins already in the G1 phase and that Fkh1 is constantly required to keep these factors bound on origins before the onset of the S phase. Together, our results suggest that dimerization-mediated stabilization of DNA binding by Fkh1 is crucial for its ability to activate DNA replication origins.

SIR telomere silencing depends on nuclear envelope lipids and modulates sensitivity to a lysolipid

The nuclear envelope (NE) is important in maintaining genome organization. The role of lipids in communication between the NE and telomere regulation was investigated, including how changes in lipid composition impact gene expression and overall nuclear architecture. Yeast was treated with the non-metabolizable lysophosphatidylcholine analog edelfosine, known to accumulate at the perinuclear ER. Edelfosine induced NE deformation and disrupted telomere clustering but not anchoring. Additionally, the association of Sir4 at telomeres decreased. RNA-seq analysis showed altered expression of Sir-dependent genes located at sub-telomeric (0-10 kb) regions, consistent with Sir4 dispersion. Transcriptomic analysis revealed that two lipid metabolic circuits were activated in response to edelfosine, one mediated by the membrane sensing transcription factors, Spt23/Mga2, and the other by a transcriptional repressor, Opi1. Activation of these transcriptional programs resulted in higher levels of unsaturated fatty acids and the formation of nuclear lipid droplets. Interestingly, cells lacking Sir proteins displayed resistance to unsaturated-fatty acids and edelfosine, and this phenotype was connected to Rap1.

Broadly Applicable Control Approaches Improve Accuracy of ChIP-Seq Data

Chromatin ImmunoPrecipitation (ChIP) is a widely used method for the analysis of protein-DNA interactions in vivo; however, ChIP has pitfalls, particularly false-positive signal enrichment that permeates the data. We have developed a new approach to control for non-specific enrichment in ChIP that involves the expression of a non-genome-binding protein targeted in the IP alongside the experimental target protein due to the sharing of epitope tags. ChIP of the protein provides a "sensor" for non-specific enrichment that can be used for the normalization of the experimental data, thereby correcting for non-specific signals and improving data quality as validated against known binding sites for several proteins that we tested, including Fkh1, Orc1, Mcm4, and Sir2. We also tested a DNA-binding mutant approach and showed that, when feasible, ChIP of a site-specific DNA-binding mutant of the target protein is likely an ideal control. These methods vastly improve our ChIP-seq results in S. cerevisiae and should be applicable in other systems.

New roles of DNA-binding and forkhead-associated domains of Fkh1 and Fkh2 in cellular functions

Two yeast forkhead transcription factors Fkh1 and Fkh2 regulate the transcription of CLB2 cluster genes important for mitosis. Both proteins contain a DNA-binding domain (DBD) and a forkhead-associated domain (FHAD), which are essential for ternary complex formation with transcription factor Mcm1, the transcription of CLB2 cluster genes and the physical interaction with Ndd1 and Clb2. Fkh2 also contains an additional C' domain that contains six consensus Cdk phosphorylation sites, but the function of this domain is dispensable. Here, we found new roles of the DBD, the FHAD, and the C' domain of Fkh1 and Fkh2 in cellular functions. The Fkh2 DBD determines the genetic interaction with NDD1, while both the FHAD and DBD of Fkh1 or Fkh2 determine cell morphology and stability of their own transcripts. Both HFADs, but not DBDs, also mediate physical interaction between Fkh1 and Fkh2. DBD and HFAD of Fkh1 and DBD, but not HFAD, of Fkh2 are also fundamental for nuclear localization. However, the Fkh2-specific C' domain has no role in these aspects except in the stability of some fkh mutant transcripts, which is either increased or decreased in the presence of this domain. These findings reveal that Fkh1 and Fkh2 have multiple cellular functions and function mainly via their DBD and FHAD through a domain-controlled feedback regulation mechanism.

The Polo kinase Cdc5 is regulated at multiple levels in the adaptation response to telomere dysfunction

Telomere dysfunction activates the DNA damage checkpoint to induce a cell cycle arrest. After an extended period of time, however, cells can bypass the arrest and undergo cell division despite the persistence of the initial damage, a process called adaptation to DNA damage. The Polo kinase Cdc5 in Saccharomyces cerevisiae is essential for adaptation and for many other cell cycle processes. How the regulation of Cdc5 in response to telomere dysfunction relates to adaptation is not clear. Here, we report that Cdc5 protein level decreases after telomere dysfunction in a Mec1-, Rad53- and Ndd1-dependent manner. This regulation of Cdc5 is important to maintain long-term cell cycle arrest but not for the initial checkpoint arrest. We find that both Cdc5 and the adaptation-deficient mutant protein Cdc5-ad are heavily phosphorylated and several phosphorylation sites modulate adaptation efficiency. The PP2A phosphatases are involved in Cdc5-ad phosphorylation status and contribute to adaptation mechanisms. We finally propose that Cdc5 orchestrates multiple cell cycle pathways to promote adaptation.

Activation of the Anaphase Promoting Complex Reverses Multiple Drug Resistant Cancer in a Canine Model of Multiple Drug Resistant Lymphoma

Simple Summary

Multiple drug resistant cancers develop all too soon in patients who received successful cancer treatment. A lack of treatment options often leaves palliative care as the last resort. We tested whether the insulin sensitizer, metformin, known to have anti-cancer activity, could impact canines with drug resistant lymphoma when added to chemotherapy. All canines in the study expressed protein markers of drug resistance and within weeks of receiving metformin, the markers were decreased. A microarray was performed, and from four canines assessed, a common set of 290 elevated genes were discovered in tumor cells compared to control cells. This cluster was enriched with genes that stall the cell cycle, with a large component representing substrates of the Anaphase Promoting Complex (APC), which degrades proteins. One canine entered partial remission. RNAs from this canine showed that APC substrates were decreased during remission and elevated again during relapse, suggesting that the APC was impaired in drug resistant canines and restored when remission occurred. We validated our results in cell lines using APC inhibitors and activators. We conclude that the APC may be a vital guardian of the genome and could delay the onset of multiple drug resistance when activated.

Abstract

Like humans, canine lymphomas are treated by chemotherapy cocktails and frequently develop multiple drug resistance (MDR). Their shortened clinical timelines and tumor accessibility make canines excellent models to study MDR mechanisms. Insulin-sensitizers have been shown to reduce the incidence of cancer in humans prescribed them, and we previously demonstrated that they also reverse and delay MDR development in vitro. Here, we treated canines with MDR lymphoma with metformin to assess clinical and tumoral responses, including changes in MDR biomarkers, and used mRNA microarrays to determine differential gene expression. Metformin reduced MDR protein markers in all canines in the study. Microarrays performed on mRNAs gathered through longitudinal tumor sampling identified a 290 gene set that was enriched in Anaphase Promoting Complex (APC) substrates and additional mRNAs associated with slowed mitotic progression in MDR samples compared to skin controls. mRNAs from a canine that went into remission showed that APC substrate mRNAs were decreased, indicating that the APC was activated during remission. In vitro validation using canine lymphoma cells selected for resistance to chemotherapeutic drugs confirmed that APC activation restored MDR chemosensitivity, and that APC activity was reduced in MDR cells. This supports the idea that rapidly pushing MDR cells that harbor high loads of chromosome instability through mitosis, by activating the APC, contributes to improved survival and disease-free duration.

Unveiling Forkhead-mediated regulation of yeast cell cycle and metabolic networks

•

Findings from genome-wide ChIP studies on budding yeast Forkheads are interpreted.

•

Power, challenges and limitation of ChIP studies are presented by target gene analysis.

•

Forkheads regulate metabolic targets through which cell division may be coordinated.

Transcription factors are regulators of the cell’s genomic landscape. By switching single genes or entire molecular pathways on or off, transcription factors modulate the precise timing of their activation. The Forkhead (Fkh) transcription factors are evolutionarily conserved to regulate organismal physiology and cell division. In addition to molecular biology and biochemical efforts, genome-wide studies have been conducted to characterize the genomic landscape potentially regulated by Forkheads in eukaryotes. Here, we discuss and interpret findings reported in six genome-wide Chromatin ImmunoPrecipitation (ChIP) studies, with a particular focus on ChIP-chip and ChIP-exo. We highlight their power and challenges to address Forkhead-mediated regulation of the cellular landscape in budding yeast. Expression changes of the targets identified in the binding assays are investigated by taking expression data for Forkhead deletion and overexpression into account. Forkheads are revealed as regulators of the metabolic network through which cell cycle dynamics may be temporally coordinated further, in addition to their well-known role as regulators of the gene cluster responsible for cell division.

Cyclin/Forkhead-mediated coordination of cyclin waves: an autonomous oscillator rationalizing the quantitative model of Cdk control for budding yeast

Networks of interacting molecules organize topology, amount, and timing of biological functions. Systems biology concepts required to pin down 'network motifs' or 'design principles' for time-dependent processes have been developed for the cell division cycle, through integration of predictive computer modeling with quantitative experimentation. A dynamic coordination of sequential waves of cyclin-dependent kinases (cyclin/Cdk) with the transcription factors network offers insights to investigate how incompatible processes are kept separate in time during the eukaryotic cell cycle. Here this coordination is discussed for the Forkhead transcription factors in light of missing gaps in the current knowledge of cell cycle control in budding yeast. An emergent design principle is proposed where cyclin waves are synchronized by a cyclin/Cdk-mediated feed-forward regulation through the Forkhead as a transcriptional timer. This design is rationalized by the bidirectional interaction between mitotic cyclins and the Forkhead transcriptional timer, resulting in an autonomous oscillator that may be instrumental for a well-timed progression throughout the cell cycle. The regulation centered around the cyclin/Cdk-Forkhead axis can be pivotal to timely coordinate cell cycle dynamics, thereby to actuate the quantitative model of Cdk control.

Granger-causal testing for irregularly sampled time series with application to nitrogen signalling in Arabidopsis

MOTIVATION

Identification of system-wide causal relationships can contribute to our understanding of long-distance, intercellular signalling in biological organisms. Dynamic transcriptome analysis holds great potential to uncover coordinated biological processes between organs. However, many existing dynamic transcriptome studies are characterized by sparse and often unevenly spaced time points that make the identification of causal relationships across organs analytically challenging. Application of existing statistical models, designed for regular time series with abundant time points, to sparse data may fail to reveal biologically significant, causal relationships. With increasing research interest in biological time series data, there is a need for new statistical methods that are able to determine causality within and between time series data sets. Here, a statistical framework was developed to identify (Granger) causal gene-gene relationships of unevenly spaced, multivariate time series data from two different tissues of Arabidopsis thaliana in response to a nitrogen signal.

RESULTS

This work delivers a statistical approach for modelling irregularly sampled bivariate signals which embeds functions from the domain of engineering that allow to adapt the model's dependence structure to the specific sampling time. Using maximum-likelihood to estimate the parameters of this model for each bivariate time series, it is then possible to use bootstrap procedures for small samples (or asymptotics for large samples) in order to test for Granger-Causality. When applied to the A.thaliana data, the proposed approach produced 3078 significant interactions, in which 2012 interactions have root causal genes and 1066 interactions have shoot causal genes. Many of the predicted causal and target genes are known players in local and long-distance nitrogen signalling, including genes encoding transcription factors, hormones and signalling peptides. Of the 1007 total causal genes (either organ), 384 are either known or predicted mobile transcripts, suggesting that the identified causal genes may be directly involved in long-distance nitrogen signalling through intercellular interactions. The model predictions and subsequent network analysis identified nitrogen-responsive genes that can be further tested for their specific roles in long-distance nitrogen signalling.

AVAILABILITY AND IMPLEMENTATION

The method was developed with the R statistical software and is made available through the R package 'irg' hosted on the GitHub repository https://github.com/SMAC-Group/irg where also a running example vignette can be found (https://smac-group.github.io/irg/articles/vignette.html). A few signals from the original data set are made available in the package as an example to apply the method and the complete A.thaliana data can be found at: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE97500.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Quantitative model of eukaryotic Cdk control through the Forkhead CONTROLLER

In budding yeast, synchronization of waves of mitotic cyclins that activate the Cdk1 kinase occur through Forkhead transcription factors. These molecules act as controllers of their sequential order and may account for the separation in time of incompatible processes. Here, a Forkhead-mediated design principle underlying the quantitative model of Cdk control is proposed for budding yeast. This design rationalizes timing of cell division, through progressive and coordinated cyclin/Cdk-mediated phosphorylation of Forkhead, and autonomous cyclin/Cdk oscillations. A "clock unit" incorporating this design that regulates timing of cell division is proposed for both yeast and mammals, and has a DRIVER operating the incompatible processes that is instructed by multiple CLOCKS. TIMERS determine whether the clocks are active, whereas CONTROLLERS determine how quickly the clocks shall function depending on external MODULATORS. This "clock unit" may coordinate temporal waves of cyclin/Cdk concentration/activity in the eukaryotic cell cycle making the driver operate the incompatible processes, at separate times.

Granger-Causal Testing for Irregularly Sampled Time Series with Application to Nitrogen Signaling in Arabidopsis

Motivation

Identification of system-wide causal relationships can contribute to our understanding of long-distance, intercellular signalling in biological organisms. Dynamic transcriptome analysis holds great potential to uncover coordinated biological processes between organs. However, many existing dynamic transcriptome studies are characterized by sparse and often unevenly spaced time points that make the identification of causal relationships across organs analytically challenging. Application of existing statistical models, designed for regular time series with abundant time points, to sparse data may fail to reveal biologically significant, causal relationships. With increasing research interest in biological time series data, there is a need for new statistical methods that are able to determine causality within and between time series data sets. Here, a statistical framework was developed to identify (Granger) causal gene-gene relationships of unevenly spaced, multivariate time series data from two different tissues of Arabidopsis thaliana in response to a nitrogen signal.

Results

This work delivers a statistical approach for modelling irregularly sampled bivariate signals which embeds functions from the domain of engineering that allow to adapt the model’s dependence structure to the specific sampling time. Using maximum-likelihood to estimate the parameters of this model for each bivariate time series, it is then possible to use bootstrap procedures for small samples (or asymptotics for large samples) in order to test for Granger-Causality. When applied to the A.thaliana data, the proposed approach produced 3078 significant interactions, in which 2012 interactions have root causal genes and 1066 interactions have shoot causal genes. Many of the predicted causal and target genes are known players in local and long-distance nitrogen signalling, including genes encoding transcription factors, hormones and signalling peptides. Of the 1007 total causal genes (either organ), 384 are either known or predicted mobile transcripts, suggesting that the identified causal genes may be directly involved in long-distance nitrogen signalling through intercellular interactions. The model predictions and subsequent network analysis identified nitrogen-responsive genes that can be further tested for their specific roles in long-distance nitrogen signalling.

Availability and implementation

The method was developed with the R statistical software and is made available through the R package ‘irg’ hosted on the GitHub repository https://github.com/SMAC-Group/irg where also a running example vignette can be found (https://smac-group.github.io/irg/articles/vignette.html). A few signals from the original data set are made available in the package as an example to apply the method and the complete A.thaliana data can be found at: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE97500.

Supplementary information

Supplementary data are available at Bioinformatics online.

Forkhead transcription factor Fkh1: insights into functional regulatory domains crucial for recruitment of Sin3 histone deacetylase complex

Transcription factors are inextricably linked with histone deacetylases leading to compact chromatin. The Forkhead transcription factor Fkh1 is mainly a negative transcriptional regulator which affects cell cycle control, silencing of mating-type cassettes and induction of pseudohyphal growth in the yeast Saccharomyces cerevisiae. Markedly, Fkh1 impinges chromatin architecture by recruiting large regulatory complexes. Implication of Fkh1 with transcriptional corepressor complexes remains largely unexplored. In this work we show that Fkh1 directly recruits corepressors Sin3 and Tup1 (but not Cyc8), providing evidence for its influence on epigenetic regulation. We also identified the specific domain of Fkh1 mediating Sin3 recruitment and substantiated that amino acids 51–125 of Fkh1 bind PAH2 of Sin3. Importantly, this part of Fkh1 overlaps with its Forkhead-associated domain (FHA). To analyse this domain in more detail, selected amino acids were replaced by alanine, revealing that hydrophobic amino acids L74 and I78 are important for Fkh1-Sin3 binding. In addition, we could prove Fkh1 recruitment to promoters of cell cycle genes CLB2 and SWI5. Notably, Sin3 is also recruited to these promoters but only in the presence of functional Fkh1. Our results disclose that recruitment of Sin3 to Fkh1 requires precisely positioned Fkh1/Sin3 binding sites which provide an extended view on the genetic control of cell cycle genes CLB2 and SWI5 and the mechanism of transcriptional repression by modulation of chromatin architecture at the G2/M transition.

The Amazing Acrobat: Yeast's Histone H3K56 Juggles Several Important Roles While Maintaining Perfect Balance

Acetylation on lysine 56 of histone H3 of the yeast Saccharomyces cerevisiae has been implicated in many cellular processes that affect genome stability. Despite being the object of much research, the complete scope of the roles played by K56 acetylation is not fully understood even today. The acetylation is put in place at the S-phase of the cell cycle, in order to flag newly synthesized histones that are incorporated during DNA replication. The signal is removed by two redundant deacetylases, Hst3 and Hst4, at the entry to G2/M phase. Its crucial location, at the entry and exit points of the DNA into and out of the nucleosome, makes this a central modification, and dictates that if acetylation and deacetylation are not well concerted and executed in a timely fashion, severe genomic instability arises. In this review, we explore the wealth of information available on the many roles played by H3K56 acetylation and the deacetylases Hst3 and Hst4 in DNA replication and repair.

Accumulation of cis- and trans-regulatory variations is associated with phenotypic divergence of a complex trait between yeast species

Gene regulatory variations accumulate during evolution and alter gene expression. While the importance of expression variation in phenotypic evolution is well established, the molecular basis remains largely unknown. Here, we examine two closely related yeast species, Saccharomyces cerevisiae and Saccharomyces paradoxus, which show phenotypical differences in morphology and cell cycle progression when grown in the same environment. By profiling the cell cycle transcriptome and binding of key transcription factors (TFs) in the two species and their hybrid, we show that changes in expression levels and dynamics of oscillating genes are dominated by upstream trans-variations. We find that multiple cell cycle regulators show both cis- and trans-regulatory variations, which alters their expression in favor of the different cell cycle phenotypes. Moreover, we show that variations in the cell cycle TFs, Fkh1, and Fkh2 affect both the expression of target genes, and the binding specificity of an interacting TF, Ace2. Our study reveals how multiple variations accumulate and propagate through the gene regulatory network, alter TFs binding, contributing to phenotypic changes in cell cycle progression.

Cell cycle transcriptomics of Capsaspora provides insights into the evolution of cyclin-CDK machinery

Progression through the cell cycle in eukaryotes is regulated on multiple levels. The main driver of the cell cycle progression is the periodic activity of cyclin-dependent kinase (CDK) complexes. In parallel, transcription during the cell cycle is regulated by a transcriptional program that ensures the just-in-time gene expression. Many core cell cycle regulators are widely conserved in eukaryotes, among them cyclins and CDKs; however, periodic transcriptional programs are divergent between distantly related species. In addition, many otherwise conserved cell cycle regulators have been lost and independently evolved in yeast, a widely used model organism for cell cycle research. For a better understanding of the evolution of the cell cycle regulation in opisthokonts, we investigated the transcriptional program during the cell cycle of the filasterean Capsaspora owczarzaki, a unicellular species closely related to animals. We developed a protocol for cell cycle synchronization in Capsaspora cultures and assessed gene expression over time across the entire cell cycle. We identified a set of 801 periodic genes that grouped into five clusters of expression over time. Comparison with datasets from other eukaryotes revealed that the periodic transcriptional program of Capsaspora is most similar to that of animal cells. We found that orthologues of cyclin A, B and E are expressed at the same cell cycle stages as in human cells and in the same temporal order. However, in contrast to human cells where these cyclins interact with multiple CDKs, Capsaspora cyclins likely interact with a single ancestral CDK1-3. Thus, the Capsaspora cyclin-CDK system could represent an intermediate state in the evolution of animal-like cyclin-CDK regulation. Overall, our results demonstrate that Capsaspora could be a useful unicellular model system for animal cell cycle regulation.

Improved recovery of cell-cycle gene expression in Saccharomyces cerevisiae from regulatory interactions in multiple omics data

BACKGROUND

Gene expression is regulated by DNA-binding transcription factors (TFs). Together with their target genes, these factors and their interactions collectively form a gene regulatory network (GRN), which is responsible for producing patterns of transcription, including cyclical processes such as genome replication and cell division. However, identifying how this network regulates the timing of these patterns, including important interactions and regulatory motifs, remains a challenging task.

RESULTS

We employed four in vivo and in vitro regulatory data sets to investigate the regulatory basis of expression timing and phase-specific patterns cell-cycle expression in Saccharomyces cerevisiae. Specifically, we considered interactions based on direct binding between TF and target gene, indirect effects of TF deletion on gene expression, and computational inference. We found that the source of regulatory information significantly impacts the accuracy and completeness of recovering known cell-cycle expressed genes. The best approach involved combining TF-target and TF-TF interactions features from multiple datasets in a single model. In addition, TFs important to multiple phases of cell-cycle expression also have the greatest impact on individual phases. Important TFs regulating a cell-cycle phase also tend to form modules in the GRN, including two sub-modules composed entirely of unannotated cell-cycle regulators (STE12-TEC1 and RAP1-HAP1-MSN4).

CONCLUSION

Our findings illustrate the importance of integrating both multiple omics data and regulatory motifs in order to understand the significance regulatory interactions involved in timing gene expression. This integrated approached allowed us to recover both known cell-cycles interactions and the overall pattern of phase-specific expression across the cell-cycle better than any single data set. Likewise, by looking at regulatory motifs in the form of TF-TF interactions, we identified sets of TFs whose co-regulation of target genes was important for cell-cycle expression, even when regulation by individual TFs was not. Overall, this demonstrates the power of integrating multiple data sets and models of interaction in order to understand the regulatory basis of established biological processes and their associated gene regulatory networks.

The Aspergillus fumigatus Phosphoproteome Reveals Roles of High-Osmolarity Glycerol Mitogen-Activated Protein Kinases in Promoting Cell Wall Damage and Caspofungin Tolerance

Aspergillus fumigatus is an opportunistic human pathogen causing allergic reactions or systemic infections, such as invasive pulmonary aspergillosis in immunocompromised patients. The mitogen-activated protein kinase (MAPK) signaling pathways are essential for fungal adaptation to the human host. Fungal cell survival, fungicide tolerance, and virulence are highly dependent on the organization, composition, and function of the cell wall. Upon cell wall stress, MAPKs phosphorylate multiple target proteins involved in the remodeling of the cell wall. Here, we investigate the global phosphoproteome of the ΔsakA and ΔmpkCA. fumigatus and high-osmolarity glycerol (HOG) pathway MAPK mutants upon cell wall damage. This showed the involvement of the HOG pathway and identified novel protein kinases and transcription factors, which were confirmed by fungal genetics to be involved in promoting tolerance of cell wall damage. Our results provide understanding of how fungal signal transduction networks modulate the cell wall. This may also lead to the discovery of new fungicide drug targets to impact fungal cell wall function, fungicide tolerance, and virulence.

The filamentous fungus Aspergillus fumigatus can cause a distinct set of clinical disorders in humans. Invasive aspergillosis (IA) is the most common life-threatening fungal disease of immunocompromised humans. The mitogen-activated protein kinase (MAPK) signaling pathways are essential to the adaptation to the human host. Fungal cell survival is highly dependent on the organization, composition, and function of the cell wall. Here, an evaluation of the global A. fumigatus phosphoproteome under cell wall stress caused by the cell wall-damaging agent Congo red (CR) revealed 485 proteins potentially involved in the cell wall damage response. Comparative phosphoproteome analyses with the ΔsakA, ΔmpkC, and ΔsakA ΔmpkC mutant strains from the osmotic stress MAPK cascades identify their additional roles during the cell wall stress response. Our phosphoproteomics allowed the identification of novel kinases and transcription factors (TFs) involved in osmotic stress and in the cell wall integrity (CWI) pathway. Our global phosphoproteome network analysis showed an enrichment for protein kinases, RNA recognition motif domains, and the MAPK signaling pathway. In contrast to the wild-type strain, there is an overall decrease of differentially phosphorylated kinases and phosphatases in ΔsakA, ΔmpkC, and ΔsakA ΔmpkC mutants. We constructed phosphomutants for the phosphorylation sites of several proteins differentially phosphorylated in the wild-type and mutant strains. For all the phosphomutants, there is an increase in the sensitivity to cell wall-damaging agents and a reduction in the MpkA phosphorylation upon CR stress, suggesting these phosphosites could be important for the MpkA modulation and CWI pathway regulation.

The Emerging Roles of Fox Family Transcription Factors in Chromosome Replication, Organization, and Genome Stability

The forkhead box (Fox) transcription factors (TFs) are widespread from yeast to humans. Their mutations and dysregulation have been linked to a broad spectrum of malignant neoplasias. They are known as critical players in DNA repair, metabolism, cell cycle control, differentiation, and aging. Recent studies, especially those from the simple model eukaryotes, revealed unexpected contributions of Fox TFs in chromosome replication and organization. More importantly, besides functioning as a canonical TF in cell signaling cascades and gene expression, Fox TFs can directly participate in DNA replication and determine the global replication timing program in a transcription-independent mechanism. Yeast Fox TFs preferentially recruit the limiting replication factors to a subset of early origins on chromosome arms. Attributed to their dimerization capability and distinct DNA binding modes, Fkh1 and Fkh2 also promote the origin clustering and assemblage of replication elements (replication factories). They can mediate long-range intrachromosomal and interchromosomal interactions and thus regulate the four-dimensional chromosome organization. The novel aspects of Fox TFs reviewed here expand their roles in maintaining genome integrity and coordinating the multiple essential chromosome events. These will inevitably be translated to our knowledge and new treatment strategies of Fox TF-associated human diseases including cancer.

ChIP-exo analysis highlights Fkh1 and Fkh2 transcription factors as hubs that integrate multi-scale networks in budding yeast

The understanding of the multi-scale nature of molecular networks represents a major challenge. For example, regulation of a timely cell cycle must be coordinated with growth, during which changes in metabolism occur, and integrate information from the extracellular environment, e.g. signal transduction. Forkhead transcription factors are evolutionarily conserved among eukaryotes, and coordinate a timely cell cycle progression in budding yeast. Specifically, Fkh1 and Fkh2 are expressed during a lengthy window of the cell cycle, thus are potentially able to function as hubs in the multi-scale cellular environment that interlocks various biochemical networks. Here we report on a novel ChIP-exo dataset for Fkh1 and Fkh2 in both logarithmic and stationary phases, which is analyzed by novel and existing software tools. Our analysis confirms known Forkhead targets from available ChIP-chip studies and highlights novel ones involved in the cell cycle, metabolism and signal transduction. Target genes are analyzed with respect to their function, temporal expression during the cell cycle, correlation with Fkh1 and Fkh2 as well as signaling and metabolic pathways they occur in. Furthermore, differences in targets between Fkh1 and Fkh2 are presented. Our work highlights Forkhead transcription factors as hubs that integrate multi-scale networks to achieve proper timing of cell division in budding yeast.

A Co-Expression Network in Hexaploid Wheat Reveals Mostly Balanced Expression and Lack of Significant Gene Loss of Homeologous Meiotic Genes Upon Polyploidization

Polyploidization has played an important role in plant evolution. However, upon polyploidization, the process of meiosis must adapt to ensure the proper segregation of increased numbers of chromosomes to produce balanced gametes. It has been suggested that meiotic gene (MG) duplicates return to a single copy following whole genome duplication to stabilize the polyploid genome. Therefore, upon the polyploidization of wheat, a hexaploid species with three related (homeologous) genomes, the stabilization process may have involved rapid changes in content and expression of MGs on homeologous chromosomes (homeologs). To examine this hypothesis, sets of candidate MGs were identified in wheat using co-expression network analysis and orthology informed approaches. In total, 130 RNA-Seq samples from a range of tissues including wheat meiotic anthers were used to define co-expressed modules of genes. Three modules were significantly correlated with meiotic tissue samples but not with other tissue types. These modules were enriched for GO terms related to cell cycle, DNA replication, and chromatin modification and contained orthologs of known MGs. Overall, 74.4% of genes within these meiosis-related modules had three homeologous copies which was similar to other tissue-related modules. Amongst wheat MGs identified by orthology, rather than co-expression, the majority (93.7%) were either retained in hexaploid wheat at the same number of copies (78.4%) or increased in copy number (15.3%) compared to ancestral wheat species. Furthermore, genes within meiosis-related modules showed more balanced expression levels between homeologs than genes in non-meiosis-related modules. Taken together, our results do not support extensive gene loss nor changes in homeolog expression of MGs upon wheat polyploidization. The construction of the MG co-expression network allowed identification of hub genes and provided key targets for future studies.

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.

Variation in Filamentous Growth and Response to Quorum-Sensing Compounds in Environmental Isolates of Saccharomyces cerevisiae

In fungi, filamentous growth is a major developmental transition that occurs in response to environmental cues. In diploid Saccharomyces cerevisiae, it is known as pseudohyphal growth and presumed to be a foraging mechanism. Rather than unicellular growth, multicellular filaments composed of elongated, attached cells spread over and into surfaces. This morphogenetic switch can be induced through quorum sensing with the aromatic alcohols phenylethanol and tryptophol. Most research investigating pseudohyphal growth has been conducted in a single lab background, Σ1278b. To investigate the natural variation in this phenotype and its induction, we assayed the diverse 100-genomes collection of environmental isolates. Using computational image analysis, we quantified the production of pseudohyphae and observed a large amount of variation. Population origin was significantly associated with pseudohyphal growth, with the West African population having the most. Surprisingly, most strains showed little or no response to exogenous phenylethanol or tryptophol. We also investigated the amount of natural genetic variation in pseudohyphal growth using a mapping population derived from a highly-heterozygous clinical isolate that contained as much phenotypic variation as the environmental panel. A bulk-segregant analysis uncovered five major peaks with candidate loci that have been implicated in the Σ1278b background. Our results indicate that the filamentous growth response is a generalized, highly variable phenotype in natural populations, while response to quorum sensing molecules is surprisingly rare. These findings highlight the importance of coupling studies in tractable lab strains with natural isolates in order to understand the relevance and distribution of well-studied traits.

The cell-penetrating FOXM1 N-terminus (M1-138) demonstrates potent inhibitory effects on cancer cells by targeting FOXM1 and FOXM1-interacting factor SMAD3

Transcription factor FOXM1 is involved in stimulating cell proliferation, enhancing DNA damage repair, promoting metastasis of cancer cells, and the inhibition of FOXM1 has been shown to prevent the initiation and progression of multiple cancers and FOXM1 is considered to be an effective target for tumor therapeutic drug development. The N-terminus of FOXM1 has been found to prevent transcriptional activities of FOXM1 and to mediate the interaction between FOXM1 and SMAD3.

Methods: A recombinant FOXM1 N-terminal domain (1-138aa) fused with a nine arginine cell-penetrating peptide is produced with an E. coli expression system and named as M1-138. The effects of M1-138 on the proliferation, migration, and tumorigenic ability of cancer cells are analyzed in vitro with cell counting, transwell assays, and colony formation assays. Electrophoretic mobility shift assays (EMSAs) and Luciferase activity assays are used to test the DNA binding ability and transcriptional activity of transcription factors. The levels of mRNAs and proteins are measured by quantitative-PCR, Western blotting or Immunohistochemistry. The interactions among proteins are analyzed with Pull-down and Co-immunoprecipitation (Co-IP) assays. The nude mouse engrafted tumor models are used to test the inhibitory effects of M1-138 in vivo.

Results: M1-138 diminishes the proliferation and migration abilities of cancer cells through binding to FOXM1 and FOXM1-interacting factor SMAD3, and consequently attenuating FOXM1 transcriptional activities from both direct and indirect FOXM1-promoter binding mechanisms and interfering with the interaction between FOXM1 and SMAD3. Treatment of M1-138 prevents tumorigenicity of cancer cells and inhibits tumor growth in nude mouse xenograft models with no obvious signs of toxicity.

Conclusion: M1-138 is a promising drug candidate for the development of anti-cancer therapeutics targeting FOXM1 and SMAD3.

FOXK2 Transcription Factor and Its Emerging Roles in Cancer

Forkhead box (FOX) transcription factors compose a large family of regulators of key biological processes within a cell. FOXK2 is a member of FOX family, whose biological functions remain relatively unexplored, despite its description in the early nineties. More recently, growing evidence has been pointing towards a role of FOXK2 in cancer, which is likely to be context-dependent and tumour-specific. Here, we provide an overview of important aspects concerning the mechanisms of regulation of FOXK2 expression and function, as well as its complex interactions at the chromatin level, which orchestrate how it differentially regulates the expression of gene targets in pathophysiology. Particularly, we explore the emerging functions of FOXK2 as a regulator of a broad range of cancer features, such as cell proliferation and survival, DNA damage, metabolism, migration, invasion and metastasis. Finally, we discuss the prognostic value of assessing FOXK2 expression in cancer patients and how it can be potentially targeted for future anticancer interventions.

The cell-cycle transcriptional network generates and transmits a pulse of transcription once each cell cycle

Multiple studies have suggested the critical roles of cyclin-dependent kinases (CDKs) as well as a transcription factor (TF) network in generating the robust cell-cycle transcriptional program. However, the precise mechanisms by which these components function together in the gene regulatory network remain unclear. Here we show that the TF network can generate and transmit a "pulse" of transcription independently of CDK oscillations. The premature firing of the transcriptional pulse is prevented by early G1 inhibitors, including transcriptional corepressors and the E3 ubiquitin ligase complex APCCdh1. We demonstrate that G1 cyclin-CDKs facilitate the activation and accumulation of TF proteins in S/G2/M phases through inhibiting G1 transcriptional corepressors (Whi5 and Stb1) and APCCdh1, thereby promoting the initiation and propagation of the pulse by the TF network. These findings suggest a unique oscillatory mechanism in which global phase-specific transcription emerges from a pulse-generating network that fires once-and-only-once at the start of the cycle.

Activating the Anaphase Promoting Complex to Enhance Genomic Stability and Prolong Lifespan

In aging cells, genomic instability is now recognized as a hallmark event. Throughout life, cells encounter multiple endogenous and exogenous DNA damaging events that are mostly repaired, but inevitably DNA mutations, chromosome rearrangements, and epigenetic deregulation begins to mount. Now that people are living longer, more and more late life time is spent suffering from age-related disease, in which genomic instability plays a critical role. However, several major questions remain heavily debated, such as the following: When does aging start? How long can we live? In order to minimize the impact of genomic instability on longevity, it is important to understand when aging starts, and to ensure repair mechanisms remain optimal from the very start to the very end. In this review, the interplay between the stress and nutrient response networks, and the regulation of homeostasis and genomic stability, is discussed. Mechanisms that link these two networks are predicted to be key lifespan determinants. The Anaphase Promoting Complex (APC), a large evolutionarily conserved ubiquitin ligase, can potentially serve this need. Recent work demonstrates that the APC maintains genomic stability, mounts a stress response, and increases longevity in yeast. Furthermore, inhibition of APC activity by glucose and nutrient response factors indicates a tight link between the APC and the stress/nutrient response networks.

Rapid Nuclear Exclusion of Hcm1 in Aging Saccharomyces cerevisiae Leads to Vacuolar Alkalization and Replicative Senescence

The yeast, Saccharomyces cerevisiae, like other higher eukaryotes, undergo a finite number of cell divisions before exiting the cell cycle due to the effects of aging. Here, we show that yeast aging begins with the nuclear exclusion of Hcm1 in young cells, resulting in loss of acidic vacuoles. Autophagy is required for healthy aging in yeast, with proteins targeted for turnover by autophagy directed to the vacuole. Consistent with this, vacuolar acidity is necessary for vacuolar function and yeast longevity. Using yeast genetics and immunofluorescence microscopy, we confirm that vacuolar acidity plays a critical role in cell health and lifespan, and is potentially maintained by a series of Forkhead Box (Fox) transcription factors. An interconnected transcriptional network involving the Fox proteins (Fkh1, Fkh2 and Hcm1) are required for transcription of v-ATPase subunits and vacuolar acidity. As cells age, Hcm1 is rapidly excluded from the nucleus in young cells, blocking the expression of Hcm1 targets (Fkh1 and Fkh2), leading to loss of v-ATPase gene expression, reduced vacuolar acidification, increased α-syn-GFP vacuolar accumulation, and finally, diminished replicative lifespan (RLS). Loss of vacuolar acidity occurs about the same time as Hcm1 nuclear exclusion and is conserved; we have recently demonstrated that lysosomal alkalization similarly contributes to aging in C. elegans following a transition from progeny producing to post-reproductive life. Our data points to a molecular mechanism regulating vacuolar acidity that signals the end of RLS when acidification is lost.

Identification of Fkh1 and Fkh2 binding site variants associated with dynamically bound DNA elements including replication origins

Forkhead Box (Fox) DNA binding proteins control multiple genome activities, including transcription, replication, and repair. These activities are organized spatially and temporally in the nucleus, and Fox proteins Fkh1 and Fkh2 have emerged as regulators of long-range chromosomal interactions involved with these activities, such as the clustering of replication origins programmed for early initiation. Fkh1 and Fkh2 bind a subset of replication origins and are thought to dimerize to mediate long-range chromosomal contacts between these origins. The binding of Fkh1 and/or Fkh2 (Fkh1/2) to replication origins and the recombination enhancer (RE), which is involved in DNA repair required for mating-type switching, is cell cycle-regulated and thus appears to be more dynamic than Fkh1/2 binding at regulated target genes. Here we report the identification of Fkh1/2 binding sequence variants at replication origins and the RE compared with Fkh1/2 binding sequences found at target genes of the CLB2 group. These different binding sequences have previously been characterized as weak and strong, respectively, suggesting that the presence of weak sites contributes to more dynamic interactions at replication origins and RE, possibly facilitated by Fkh1/2 dimerization and cooperative interactions with accessory proteins. We discuss the wealth of regulatory potential imbued in these features of the DNA and its binding proteins.

Reconciling conflicting models for global control of cell-cycle transcription

Models for the control of global cell-cycle transcription have advanced from a CDK-APC/C oscillator, a transcription factor (TF) network, to coupled CDK-APC/C and TF networks. Nonetheless, current models were challenged by a recent study that concluded that the cell-cycle transcriptional program is primarily controlled by a CDK-APC/C oscillator in budding yeast. Here we report an analysis of the transcriptome dynamics in cyclin mutant cells that were not queried in the previous study. We find that B-cyclin oscillation is not essential for control of phase-specific transcription. Using a mathematical model, we demonstrate that the function of network TFs can be retained in the face of significant reductions in transcript levels. Finally, we show that cells arrested at mitotic exit with non-oscillating levels of B-cyclins continue to cycle transcriptionally. Taken together, these findings support a critical role of a TF network and a requirement for CDK activities that need not be periodic.

Flexible model-based clustering of mixed binary and continuous data: application to genetic regulation and cancer

Clustering is used widely in ‘omics’ studies and is often tackled with standard methods, e.g. hierarchical clustering. However, the increasing need for integration of multiple data sets leads to a requirement for clustering methods applicable to mixed data types, where the straightforward application of standard methods is not necessarily the best approach. A particularly common problem involves clustering entities characterized by a mixture of binary data (e.g. presence/absence of mutations, binding, motifs and epigenetic marks) and continuous data (e.g. gene expression, protein abundance, metabolite levels). Here, we present a generic method based on a probabilistic model for clustering this type of data, and illustrate its application to genetic regulation and the clustering of cancer samples. We show that the resulting clusters lead to useful hypotheses: in the case of genetic regulation these concern regulation of groups of genes by specific sets of transcription factors and in the case of cancer samples combinations of gene mutations are related to patterns of gene expression. The clusters have potential mechanistic significance and in the latter case are significantly linked to survival. The method is available as a stand-alone software package (GNU General Public Licence) from http://github.com/BioToolsLeeds/FlexiCoClusteringPackage.git.

Design of experiment for nonlinear dynamic gene regulatory network identification

The gene regulatory network (GRN) is critical for understanding the regulatory interaction between genes. Time-course microarray experiments provide ample information for constructing GRN. The designs for microarray experiments serve different purposes. However, the experiment design specifically for GRN identification is still sparse. In this article, we use a simulation-based approach to deal with design problems in the framework of nonparametric differential equations. We investigate a number of feasible designs. In particular, we evaluate whether earlier samplings can result in more useful information for GRN identification. We also evaluate the effectiveness of two strategies: more frequent samplings per replicate with fewer replicates versus fewer samplings per replicate with more replicates while keeping the total number of samplings constant. The results of our investigation provide quantitative guidance for designing and selecting microarray experiments for the purpose of GRN identification.

Form and function of topologically associating genomic domains in budding yeast

The genome of metazoan cells is organized into topologically associating domains (TADs) that have similar histone modifications, transcription level, and DNA replication timing. Although similar structures appear to be conserved in fission yeast, computational modeling and analysis of high-throughput chromosome conformation capture (Hi-C) data have been used to argue that the small, highly constrained budding yeast chromosomes could not have these structures. In contrast, herein we analyze Hi-C data for budding yeast and identify 200-kb scale TADs, whose boundaries are enriched for transcriptional activity. Furthermore, these boundaries separate regions of similarly timed replication origins connecting the long-known effect of genomic context on replication timing to genome architecture. To investigate the molecular basis of TAD formation, we performed Hi-C experiments on cells depleted for the Forkhead transcription factors, Fkh1 and Fkh2, previously associated with replication timing. Forkhead factors do not regulate TAD formation, but do promote longer-range genomic interactions and control interactions between origins near the centromere. Thus, our work defines spatial organization within the budding yeast nucleus, demonstrates the conserved role of genome architecture in regulating DNA replication, and identifies a molecular mechanism specifically regulating interactions between pericentric origins.

Conserved forkhead dimerization motif controls DNA replication timing and spatial organization of chromosomes in S. cerevisiae

Forkhead Box (Fox) proteins share the Forkhead domain, a winged-helix DNA binding module, which is conserved among eukaryotes from yeast to humans. These sequence-specific DNA binding proteins have been primarily characterized as transcription factors regulating diverse cellular processes from cell cycle control to developmental fate, deregulation of which contributes to developmental defects, cancer, and aging. We recently identified Saccharomyces cerevisiae Forkhead 1 (Fkh1) and Forkhead 2 (Fkh2) as required for the clustering of a subset of replication origins in G₁ phase and for the early initiation of these origins in the ensuing S phase, suggesting a mechanistic role linking the spatial organization of the origins and their activity. Here, we show that Fkh1 and Fkh2 share a unique structural feature of human FoxP proteins that enables FoxP2 and FoxP3 to form domain-swapped dimers capable of bridging two DNA molecules in vitro. Accordingly, Fkh1 self-associates in vitro and in vivo in a manner dependent on the conserved domain-swapping region, strongly suggestive of homodimer formation. Fkh1- and Fkh2-domain-swap-minus (dsm) mutations are functional as transcription factors yet are defective in replication origin timing control. Fkh1-dsm binds replication origins in vivo but fails to cluster them, supporting the conclusion that Fkh1 and Fkh2 dimers perform a structural role in the spatial organization of chromosomal elements with functional importance.

Transcriptional and post-transcriptional regulation of Cdc20 during the spindle assembly checkpoint in S. cerevisiae

The anaphase-promoting complex (APC) is a ubiquitin ligase responsible for promoting the degradation of many cell cycle regulators. One of the activators and substrate-binding proteins for the APC is Cdc20. It has been shown previously that Cdc20 can promote its own degradation by the APC in normal cycling cells mainly through a cis-degradation mode (i.e. via an intramolecular mechanism). However, how Cdc20 is degraded during the spindle assembly checkpoint (SAC) is still not fully clear. In this study, we used a dual-Cdc20 system to investigate this issue and found that the cis-degradation mode is also the major pathway responsible for Cdc20 degradation during the SAC. In addition, we found that there is an inverse relationship between APC^Cdc20 activity and the transcriptional activity of the CDC20 promoter, which likely occurs through feedback regulation by APC^Cdc20 substrates, such as the cyclins Clb2 and Clb5. These findings contribute to our understanding of how the inhibition of APC^Cdc20 activity and enhanced Cdc20 degradation are required for proper spindle checkpoint arrest.

Inference of gene regulation functions from dynamic transcriptome data

To quantify gene regulation, a function is required that relates transcription factor binding to DNA (input) to the rate of mRNA synthesis from a target gene (output). Such a ‘gene regulation function’ (GRF) generally cannot be measured because the experimental titration of inputs and simultaneous readout of outputs is difficult. Here we show that GRFs may instead be inferred from natural changes in cellular gene expression, as exemplified for the cell cycle in the yeast S. cerevisiae. We develop this inference approach based on a time series of mRNA synthesis rates from a synchronized population of cells observed over three cell cycles. We first estimate the functional form of how input transcription factors determine mRNA output and then derive GRFs for target genes in the CLB2 gene cluster that are expressed during G2/M phase. Systematic analysis of additional GRFs suggests a network architecture that rationalizes transcriptional cell cycle oscillations. We find that a transcription factor network alone can produce oscillations in mRNA expression, but that additional input from cyclin oscillations is required to arrive at the native behaviour of the cell cycle oscillator.

DOI:http://dx.doi.org/10.7554/eLife.12188.001

Living cells rely on networks of genes to control their behavior, including how they grow, develop and respond to stress. Genes encode instructions needed to make proteins and other molecules, and much of the control is exerted at the first stage of protein production, known as transcription. During this process, a gene is copied to make molecules known as transcripts that may later be used as templates to make proteins.

Many genes encode proteins that act to regulate transcription. Therefore, an individual gene may receive inputs from other genes, and these inputs affect how much transcript the gene produces, which can be considered as the gene’s output. While these inputs and outputs can often be wired together to form a network, it is less clear exactly how all the different inputs at a gene interact to determine its output. These interactions are known as “gene regulation functions”, and knowing them would be an important step towards understanding gene networks, which would help us to predict how cells will behave in different situations.

Gene regulation functions are difficult to measure directly, so researchers would like to find other ways to assess them indirectly. A recently developed experimental technique called “dynamic transcriptome analysis” seemed promising as it measures both the inputs and outputs of all genes in a cell over time. Hillenbrand et al. used this technique to infer gene regulation functions with one or two inputs in yeast cells.

Comparing these estimates with experimental data from previous studies showed that these inferred gene regulation functions could successfully predict the output of a gene based on its inputs. Hillenbrand et al. then used these estimates to search and model a well-known genetic network that is thought to be part of the molecular clockwork that controls the timing of events that cause a cell to divide.

Currently, the approach used by Hillenbrand et al. treats gene regulation functions like “black boxes”. This means that, while an output can be predicted if the inputs are known, it cannot reveal all of the detailed mechanisms behind it. Gaining insights into the inner workings of these black boxes will require taking more data into account, such as how abundant the proteins that regulate transcription are, where they are located within cells or whether they are active or not. Therefore, the next challenge is to incorporate these kinds of data to gain a fuller picture of how gene networks operate within cells.

DOI:http://dx.doi.org/10.7554/eLife.12188.002

Decoding common and divergent cellular functions of the domains of forkhead transcription factors Fkh1 and Fkh2

Forkhead transcription factors play a key role in embryonic patterning during development. In Saccharomyces cerevisiae, two forkhead transcription factors, Fkh1 and Fkh2, regulate the transcription of CLB2 cluster genes important for mitosis. Fkh1 reduces, whereas Fkh2 elevates, the transcription of CLB2 cluster genes. However, the mechanism for this observation remains unclear. Fkh1 and Fkh2 each contain a forkhead domain (DNA-binding domain, DBD) and a forkhead-associated domain (FHAD), whereas Fkh2 possesses an extra C' domain containing six consensus cyclin-dependent kinase phosphorylation sites. In the present study, roles of these domains in protein complexes, the regulation of cell growth and CLB2 cluster genes and protein interactions were investigated using various domain mutants. The result showed that the DBD was vital for ternary complex formation with Mcm1, whereas the FHAD was central for the regulation of cell growth and CLB2 cluster transcription and for interactions with Ndd1 and Clb2. However, the Fkh2 C' domain was dispensable for the above functions. Both DBDs and FHADs had functional divergences in the cell, and Ndd1 functioned via its phosphorylated form. These data provide important insights into the functional mechanism of Fkh1 and Fkh2 in cell cycle control.

Binding of the Fkh1 Forkhead Associated Domain to a Phosphopeptide within the Mph1 DNA Helicase Regulates Mating-Type Switching in Budding Yeast

The Saccharomyces cerevisiae Fkh1 protein has roles in cell-cycle regulated transcription as well as a transcription-independent role in recombination donor preference during mating-type switching. The conserved FHA domain of Fkh1 regulates donor preference by juxtaposing two distant regions on chromosome III to promote their recombination. A model posits that this Fkh1-mediated long-range chromosomal juxtaposition requires an interaction between the FHA domain and a partner protein(s), but to date no relevant partner has been described. In this study, we used structural modeling, 2-hybrid assays, and mutational analyses to show that the predicted phosphothreonine-binding FHA domain of Fkh1 interacted with multiple partner proteins. The Fkh1 FHA domain was important for its role in cell-cycle regulation, but no single interaction partner could account for this role. In contrast, Fkh1’s interaction with the Mph1 DNA repair helicase regulated donor preference during mating-type switching. Using 2-hybrid assays, co-immunoprecipitation, and fluorescence anisotropy, we mapped a discrete peptide within the regulatory Mph1 C-terminus required for this interaction and identified two threonines that were particularly important. In vitro binding experiments indicated that at least one of these threonines had to be phosphorylated for efficient Fkh1 binding. Substitution of these two threonines with alanines (mph1-2TA) specifically abolished the Fkh1-Mph1 interaction in vivo and altered donor preference during mating-type switching to the same degree as mph1Δ. Notably, the mph1-2TA allele maintained other functions of Mph1 in genome stability. Deletion of a second Fkh1-interacting protein encoded by YMR144W also resulted in a change in Fkh1-FHA-dependent donor preference. We have named this gene FDO1 for Forkhead one interacting protein involved in donor preference. We conclude that a phosphothreonine-mediated protein-protein interface between Fkh1-FHA and Mph1 contributes to a specific long-range chromosomal interaction required for mating-type switching, but that Fkh1-FHA must also interact with several other proteins to achieve full functionality in this process.

Specific chromosomal interactions between distal regions of the genome allow for DNA transactions necessary for normal cell function, but the protein-protein interfaces that regulate such interactions remain largely unknown. The budding yeast Fkh1 protein uses its evolutionarily conserved phosphothreonine-binding FHA domain to regulate a long-range DNA transaction called mating-type switching that allows yeast cells to switch their sexual phenotype. In this study, another conserved nuclear protein, the Mph1 DNA repair helicase, was shown to interact directly with the FHA domain of Fkh1 to regulate mating-type switching. The Fkh1-Mph1 interaction required two phosphorylated threonines on Mph1 that were dispensable for many other Mph1-protein interactions and other Mph1 chromosomal functions. Thus a discrete protein-protein interface between two multifunctional chromosomal proteins helps define a long-range chromosomal interaction important for controlling cell behavior.

Regulation of wee1(+) expression during meiosis in fission yeast

In eukaryotes, the cyclin-dependent kinase Cdk1p (Cdc2p) plays a central role in entry into and progression through nuclear division during mitosis and meiosis. Cdk1p is activated during meiotic nuclear divisions by dephosphorylation of its tyrosine-15 residue. The phosphorylation status of this residue is largely determined by the Wee1p kinase and the Cdc25p phosphatase. In fission yeast, the forkhead-type transcription factor Mei4p is essential for entry into the first meiotic nuclear division. We recently identified cdc25⁺ as an essential target of Mei4p in the control of entry into meiosis I. Here, we show that wee1⁺ is another important target of Mei4p in the control of entry into meiosis I. Mei4p bound to the upstream region of wee1⁺ in vivo and in vitro and inhibited expression of wee1⁺, whereas Mei4p positively regulated expression of the adjacent pseudogene. Overexpression of Mei4p inhibited expression of wee1⁺ and induced that of the pseudogene. Conversely, deletion of Mei4p did not decrease expression of wee1⁺ but inhibited that of the pseudogene. In addition, deletion of Mei4p-binding regions delayed repression of wee1⁺ expression as well as induction of expression of the pseudogene. These results suggest that repression of wee1⁺ expression is primarily owing to Mei4p-mediated transcriptional interference.

Decoupling of divergent gene regulation by sequence-specific DNA binding factors

Divergent gene pairs (DGPs) are abundant in eukaryotic genomes. Since two genes in a DGP potentially share the same regulatory sequence, one might expect that they should be co-regulated. However, an inspection of yeast DGPs containing cell-cycle or stress response genes revealed that most DGPs are differentially-regulated. The mechanism underlying DGP differential regulation is not understood. Here, we showed that co- versus differential regulation cannot be explained by genetic features including promoter length, binding site orientation, TATA elements, nucleosome distribution, or presence of non-coding RNAs. Using time-lapse fluorescence microscopy, we carried out an in-depth study of a differentially regulated DGP, PFK26-MOB1. We found that their differential regulation is mainly achieved through two DNA-binding factors, Tbf1 and Mcm1. Similar to 'enhancer-blocking insulators' in higher eukaryotes, these factors shield the proximal promoter from the action of more distant transcription regulators. We confirmed the blockage function of Tbf1 using synthetic promoters. We further presented evidence that the blockage mechanism is widely used among genome-wide DGPs. Besides elucidating the DGP regulatory mechanism, our work revealed a novel class of insulators in yeast.

Behavior of replication origins in Eukaryota - spatio-temporal dynamics of licensing and firing

Although every organism shares some common features of replication, this process varies greatly among eukaryotic species. Current data show that mathematical models of the organization of origins based on possibility theory may be applied (and remain accurate) in every model organism i.e. from yeast to humans. The major differences lie within the dynamics of origin firing and the regulation mechanisms that have evolved to meet new challenges throughout the evolution of the organism. This article elaborates on the relations between chromatin structure, organization of origins, their firing times and the impact that these features can have on genome stability, showing both differences and parallels inside the eukaryotic domain.

Degradation of Ndd1 by APC/C(Cdh1) generates a feed forward loop that times mitotic protein accumulation

Ndd1 activates the Mcm1-Fkh2 transcription factor to transcribe mitotic regulators. The anaphase-promoting complex/cyclosome activated by Cdh1 (APC/C(Cdh1)) mediates the degradation of proteins throughout G1. Here we show that the APC/C(Cdh1) ubiquitinates Ndd1 and mediates its degradation, and that APC/C(Cdh1) activity suppresses accumulation of Ndd1 targets. We confirm putative Ndd1 targets and identify novel ones, many of them APC/C(Cdh1) substrates. The APC/C(Cdh1) thus regulates these proteins in a dual manner—both pretranscriptionally and post-translationally, forming a multi-layered feedforward loop (FFL). We predict by mathematical modelling and verify experimentally that this FFL introduces a lag between APC/C(Cdh1) inactivation at the end of G1 and accumulation of genes transcribed by Ndd1 in G2. This regulation generates two classes of APC/C(Cdh1) substrates, early ones that accumulate in S and late ones that accumulate in G2. Our results show how the dual state APC/C(Cdh1) activity is converted into multiple outputs by interactions between its substrates.

Transcriptional control of fungal cell cycle and cellular events by Fkh2, a forkhead transcription factor in an insect pathogen

Transcriptional control of the cell cycle by forkhead (Fkh) transcription factors is likely associated with fungal adaptation to host and environment. Here we show that Fkh2, an ortholog of yeast Fkh1/2, orchestrates cell cycle and many cellular events of Beauveria bassiana, a filamentous fungal insect pathogen. Deletion of Fkh2 in B. bassiana resulted in dramatic down-regulation of the cyclin-B gene cluster and hence altered cell cycle (longer G₂/M and S, but shorter G₀/G₁, phases) in unicellular blastospores. Consequently, ΔFkh2 produced twice as many, but smaller, blastospores than wild-type under submerged conditions, and formed denser septa and shorter/broader cells in aberrantly branched hyphae. In these hyphae, clustered genes required for septation and conidiation were remarkedly up-regulated, followed by higher yield and slower germination of aerial conidia. Moreover, ΔFkh2 displayed attenuated virulence and decreased tolerance to chemical and environmental stresses, accompanied with altered transcripts and activities of phenotype-influencing proteins or enzymes. All the changes in ΔFkh2 were restored by Fkh2 complementation. All together, Fkh2-dependent transcriptional control is vital for the adaptation of B. bassiana to diverse habitats of host insects and hence contributes to its biological control potential against arthropod pests.

A new transcription factor for mitosis: in Schizosaccharomyces pombe, the RFX transcription factor Sak1 works with forkhead factors to regulate mitotic expression

Mitotic genes are one of the most strongly oscillating groups of genes in the eukaryotic cell cycle. Understanding the regulation of mitotic gene expression is a key issue in cell cycle control but is poorly understood in most organisms. Here, we find a new mitotic transcription factor, Sak1, in the fission yeast Schizosaccharomyces pombe. Sak1 belongs to the RFX family of transcription factors, which have not previously been connected to cell cycle control. Sak1 binds upstream of mitotic genes in close proximity to Fkh2, a forkhead transcription factor previously implicated in regulation of mitotic genes. We show that Sak1 is the major activator of mitotic gene expression and also confirm the role of Fkh2 as the opposing repressor. Sep1, another forkhead transcription factor, is an activator for a small subset of mitotic genes involved in septation. From yeasts to humans, forkhead transcription factors are involved in mitotic gene expression and it will be interesting to see whether RFX transcription factors may also be involved in other organisms.

Ndd1 turnover by SCF(Grr1) is inhibited by the DNA damage checkpoint in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, Ndd1 is the dedicated transcriptional activator of the mitotic gene cluster, which includes thirty-three genes that encode key mitotic regulators, making Ndd1 a hub for the control of mitosis. Previous work has shown that multiple kinases, including cyclin-dependent kinase (Cdk1), phosphorylate Ndd1 to regulate its activity during the cell cycle. Previously, we showed that Ndd1 was inhibited by phosphorylation in response to DNA damage. Here, we show that Ndd1 is also subject to regulation by protein turnover during the mitotic cell cycle: Ndd1 is unstable during an unperturbed cell cycle, but is strongly stabilized in response to DNA damage. We find that Ndd1 turnover in metaphase requires Cdk1 activity and the ubiquitin ligase SCF^Grr1. In response to DNA damage, Ndd1 stabilization requires the checkpoint kinases Mec1/Tel1 and Swe1, the S. cerevisiae homolog of the Wee1 kinase. In both humans and yeast, the checkpoint promotes Wee1-dependent inhibitory phosphorylation of Cdk1 following exposure to DNA damage. While this is critical for checkpoint-induced arrest in most organisms, this is not true in budding yeast, where the function of damage-induced inhibitory phosphorylation is less well understood. We propose that the DNA damage checkpoint stabilizes Ndd1 by inhibiting Cdk1, which we show is required for targeting Ndd1 for destruction.

All cells must regulate cell division in response to extracellular and intracellular cues, and one of the most critical steps to regulate is the process of cell division, or mitosis. In response to DNA damage in budding yeast, cells activate a checkpoint that promotes DNA repair and arrests the cell cycle before division to give the cell time to repair the lesion. One of the key regulators of mitosis is an essential transcription factor called Ndd1. Ndd1 is known to be regulated by transcription and phosphorylation, both in unperturbed cells and following exposure to DNA damage. Here, we show that Ndd1 protein turnover is also regulated in both situations. Ndd1 is degraded quickly during an unperturbed cell cycle, but is strongly stabilized following exposure to DNA damage. We identify the machinery that targets Ndd1 for turnover and the signaling pathways required to stabilize Ndd1 in response to DNA damage. Maintaining high levels of Ndd1 after exposure to DNA damage may allow the cell to reactivate Ndd1 after the damage has been repaired to promote mitosis.

Cell cycle-independent phospho-regulation of Fkh2 during hyphal growth regulates Candida albicans pathogenesis

The opportunistic human fungal pathogen, Candida albicans, undergoes morphological and transcriptional adaptation in the switch from commensalism to pathogenicity. Although previous gene-knockout studies have identified many factors involved in this transformation, it remains unclear how these factors are regulated to coordinate the switch. Investigating morphogenetic control by post-translational phosphorylation has generated important regulatory insights into this process, especially focusing on coordinated control by the cyclin-dependent kinase Cdc28. Here we have identified the Fkh2 transcription factor as a regulatory target of both Cdc28 and the cell wall biosynthesis kinase Cbk1, in a role distinct from its conserved function in cell cycle progression. In stationary phase yeast cells 2D gel electrophoresis shows that there is a diverse pool of Fkh2 phospho-isoforms. For a short window on hyphal induction, far before START in the cell cycle, the phosphorylation profile is transformed before reverting to the yeast profile. This transformation does not occur when stationary phase cells are reinoculated into fresh medium supporting yeast growth. Mass spectrometry and mutational analyses identified residues phosphorylated by Cdc28 and Cbk1. Substitution of these residues with non-phosphorylatable alanine altered the yeast phosphorylation profile and abrogated the characteristic transformation to the hyphal profile. Transcript profiling of the phosphorylation site mutant revealed that the hyphal phosphorylation profile is required for the expression of genes involved in pathogenesis, host interaction and biofilm formation. We confirmed that these changes in gene expression resulted in corresponding defects in pathogenic processes. Furthermore, we identified that Fkh2 interacts with the chromatin modifier Pob3 in a phosphorylation-dependent manner, thereby providing a possible mechanism by which the phosphorylation of Fkh2 regulates its specificity. Thus, we have discovered a novel cell cycle-independent phospho-regulatory event that subverts a key component of the cell cycle machinery to a role in the switch from commensalism to pathogenicity.

The fungus Candida albicans is a commensal in the human microbiota, responsible for superficial infections such as oral and vaginal thrush. However, it can become highly virulent, causing life-threatening systemic candidemia in severely immunocompromised patients, including those taking immunosuppressive drugs for transplantation, sufferers of AIDS and neutropenia, and individuals undergoing chemotherapy or at extremes of age. With a rapidly increasing ageing population worldwide, C. albicans and other fungal pathogens will become more prevalent, demanding a greater understanding of their pathogenesis for the development of effective therapeutics. Fungal pathogenicity requires a coordinated change in the pattern of gene expression orchestrated by a set of transcription factors. Here we have discovered that a transcription factor, Fkh2, is modified by phosphorylation under the control of the kinases Cdc28 and Cbk1 in response to conditions that activate virulence factor expression. Fkh2 is involved in a wide variety of cellular processes including cell proliferation, but this phosphorylation endows it with a specialized function in promoting the expression of genes required for tissue invasion, biofilm formation, and pathogenesis in the host. This study highlights the role of protein phosphorylation in regulating pathogenesis and furthers our understanding of the pathogenic switch in this important opportunistic fungal pathogen.

Dominant Epistasis Between Two Quantitative Trait Loci Governing Sporulation Efficiency in Yeast Saccharomyces cerevisiae

Sporulation efficiency in the yeast Saccharomyces cerevisiae is a well-established model for studying quantitative traits. A variety of genes and nucleotides causing different sporulation efficiencies in laboratory, as well as in wild strains, has already been extensively characterised (mainly by reciprocal hemizygosity analysis and nucleotide exchange methods). We applied a different strategy in order to analyze the variation in sporulation efficiency of laboratory yeast strains. Coupling classical quantitative genetic analysis with simulations of phenotypic distributions (a method we call phenotype modelling) enabled us to obtain a detailed picture of the quantitative trait loci (QTLs) relationships underlying the phenotypic variation of this trait. Using this approach, we were able to uncover a dominant epistatic inheritance of loci governing the phenotype. Moreover, a molecular analysis of known causative quantitative trait genes and nucleotides allowed for the detection of novel alleles, potentially responsible for the observed phenotypic variation. Based on the molecular data, we hypothesise that the observed dominant epistatic relationship could be caused by the interaction of multiple quantitative trait nucleotides distributed across a 60--kb QTL region located on chromosome XIV and the RME1 locus on chromosome VII. Furthermore, we propose a model of molecular pathways which possibly underlie the phenotypic variation of this trait.