m:Explorer: multinomial regression models reveal positive and negative regulators of longevity in yeast quiescence

Jüri Reimand
Anu Aun
Jaak Vilo
Juan M Vaquerizas
Juhan Sedman
Nicholas M Luscombe
m:Explorer: multinomial regression models reveal positive and negative regulators of longevity in yeast quiescence

We developed m:Explorer for identifying process-specific transcription factors (TFs) from multiple genome-wide sources, including transcriptome, DNA-binding and chromatin data. m:Explorer robustly outperforms similar techniques in finding cell cycle TFs in Saccharomyces cerevisiae. We predicted and experimentally tested regulators of quiescence (G 0 ), a model of ageing, over a six-week time-course. We validated nine of top-12 predictions as novel G 0 TFs, including Δmga2, Δcst6, Δbas1 with higher viability and G 0 -essential TFs Tup1, Swi3. Pathway analysis associates longevity to reduced growth, reprogrammed metabolism and cell wall remodeling. m:Explorer (http:// biit.cs.ut.ee/mexplorer/) is instrumental in interrogating eukaryotic regulatory systems using heterogeneous data.

Background

Eukaryotic transcriptional regulation is a core cellular process that governs the expression of genes. Understanding gene expression is crucial in explaining complex biological processes including development, disease and cancer. Transcription factors (TF) are key proteins that activate or repress transcription by binding sequence-specifically to DNA in promoter regions of target genes. Mapping such regulatory networks and TF functions is therefore an important goal of current biomedical research. In complex vertebrate organisms like human, this task is hindered by enormous genomic space, numerous cell types, and distinct experimental procedures with data that is often unsuitable for direct comparison. The relatively simple unicellular model organism budding yeast (S. cerevisiae) serves as a platform for regulatory genomics. Multiple types of globalscale data of yeast gene regulation are available to date, including microarrays with TF deletion (ΔTF) strains [1,2], predictions of TF binding sites (TFBS) [3][4][5], and measurements of chromatin state such as nucleosome positioning [6]. These data appear to be complete, however the agreement between transcript expression and TF binding events remains modest [2,7]. While part of this controversy can be attributed to experimental and statistical noise, we may still lack significant details regarding the biological relationships among such heterogeneous information. Consequently high-throughput data constitute less reliable evidence and much functional knowledge is extracted from careful and expensive focused studies. Most TFs and their exact roles in cellular processes remain poorly understood. Therefore biologically meaningful computational analysis is an important challenge in deciphering cellular regulatory networks.

Computational prediction of TF function from gene expression and DNA binding data is an active area of research. Numerous algorithms have been published elsewhere, albeit few have been validated experimentally. Earliest approaches focused on a specific class of data and used alternative types of evidence for computational validation. For instance, microarray clustering followed by DNA motif discovery in gene promoters helped establish the genome-scale link between mRNA expression profiles and TF binding [8,9]. Similarly, analysis of cell cycle expression patterns of TF-bound genes led to recovery of cell cycle TFs [10]. More recent methods use statistical modeling to integrate multiple types of evidence. For example, ARACNE extracts transcriptional networks from numeric microarray data using mutual information [11], and MARINA is a down-stream method that identifies master regulators of these networks through association tests with TF binding target genes [12]. The SAMBA biclustering algorithm studies matrices of regulators and target genes, and highlights regulatory relationships between genes and TFs that co-occur in clusters [13]. The linear regression method REDUCE integrates numeric microarray data, DNA sequence and TF affinity matrices by modeling the linear relationship between gene expression levels and TF-DNA interactions [14]. The GeneClass algorithm additionally integrates information about gene function, as it constructs decision trees of discrete microarray profiles and TF binding sites to select predictors of process-specific genes [15]. While this method provides direct modeling of gene function, TFs and gene expression data are studied as independent predictors. Notably, none of the above methods take advantage of recent ΔTF microarrays that reveal regulator target genes [1,2]. Nested effects models are designed to extract regulatory networks from perturbation data [16], although integration of TFBS and gene annotations is not supported. Nucleosome positioning measurements also remain unexplored in all above approaches. In summary, additional computational efforts are required for meaningful integration of versatile biological data.

Here we propose a method m:Explorer that uses multinomial logistic regression models to predict process-specific transcription factors. We aim to provide the following improvements in comparison to earlier methods. First, our method allows simultaneous analysis of four classes of data: (i) gene expression data, including perturbation screens, (ii) TF binding sites, (iii) chromatin state in gene promoters, and (iv) functional gene classification. The model is based on the assumption that TF target genes from perturbation screens and TF binding assays are equally informative about TF process specificity. Second, we reduce noise by including only high-confidence regulatory relationships, and do not assume linear relationships between regulators and target genes. Third, we integrate detailed information to better reflect underlying biology: multiple subprocesses may be studied in a single model, and chromatin state data are incorporated into TF binding site analysis. TF target genes with simultaneous evidence from gene expression and TFBS data are highlighted separately. Fourth, our analysis is robust to highly redundant biological networks, as statistical independence is not required. We use univariate models to study all TFs independently and avoid overfitting that is characteristic to many model-based approaches. This is statistically valid under the assumption that a complex model may be understood by examining its components.

To test our method, we compiled a comprehensive dataset covering most TFs of the budding yeast. We benchmarked m:Explorer in a well-studied biological system and establish its improved performance in comparison to several similar methods. Then we used the tool to discover regulators of quiescence (G 0 , stationary phase), a cellular resting state that serves as a model of chronological ageing. Experimental validations of our predictions revealed nine TFs with significant impact on G 0 viability. Besides demonstrating the applicability of our computational method, these findings are of great potential interest to yeast biologists and researchers of G 0 -related processes like ageing, development and cancer.

Results m:Explorer -multinomial logistic regression for inferring process-specific gene regulation

Here we tackle the problem of identifying transcription factors that regulate process-specific genes ( Figure 1). Our model m:Explorer uses three types of independent regulatory information to characterize target genes of TFs: gene expression measurements from TF perturbation screens, TF binding sites in gene promoters and DNA nucleosome occupancy in binding sites. The fourth input is a list of process-specific genes for which potential transcriptional regulators are sought.

The first stage of our analysis involves data preprocessing and discretization in which high-confidence TF target genes are identified from multiple sources ( Figure  1A). We assumed that genes responding to TF perturbation are likely targets of the regulator. We previously analyzed a large collection of ΔTF microarrays, extracted genes with significant up or down-regulation (moderated t-test, FDR p ≤ 0.05), and assigned these to perturbed regulators (Step 1, methods described in [2]). We also followed the assumption that TF binding in promoters is likely to indicate regulation of downstream genes, and binding sites in low nucleosome occupancy regions are more likely targets of TFs. We collected TF-DNA interactions from multiple datasets and classified genes as TF-bound if at least one dataset showed significant binding in 600 bp promoters (Step 2). We further categorized our TFBS collection into nucleosomedepleted TFBS (onesided t-test, FDR p ≤ 0.05) and sites with no nucleosome depletion. Next we integrated TF target genes into a genome-wide matrix, by assigning non-related genes to a baseline class and creating extra classes for genes with multiple evidence (Steps 3, 5).

Besides regulatory targets of transcription factors, our method requires a list of process-specific genes for which potential regulators are predicted. These may originate from literature, additional microarray datasets, pathway databases or biomedical ontologies. Several non-overlapping lists of genes may be provided to integrate further information about sub-process specificity, sample treatment or differential expression. These genes are organized similarly to TF targets (Steps 4, 5).

The second stage of our analysis involves multinomial regression analysis of process-specific genes and TF targets ( Figure 1B). It is a generalization of linear regression that associates a multi-class categorical response (process-specific genes) with one or more predictors (TF target genes). Through the logistic transformation, each gene is assigned a log-odds probability of being process-specific given its relation to a particular TF, as
g(y i ) = log p i,c p i,C = β 0,c + K k=1 β k,c x i,k ,
where y i is the process annotation of the i-th gene, and p i,c is the probability that gene i is part of sub-process c, given a linear combination of K types of evidence x X regarding TF target genes. All probabilities are computed relative to the baseline genes denoted by class C. The TF relation to process genes is quantified through regression coefficients b such that positive coefficients reflect a higher probability of TF target genes involving in the given process. Coefficients b are sought iteratively in maximum likelihood estimation. Likelihood reflects the estimated probabilities of all N genes belonging to their actual class, and thus provides a measure for model evaluation: For every TF, assign non-associated genes to the baseline class (no regulatory data ). Apply same treatment to genes unrelated to the process .   Figure 1A: Data preprocessing. High-confidence TF target genes of four primary classes are collected from several datasets (Steps 1-2) and merged into composite lists with four extra classes for multiple lines of evidence (Step 3). Process-specific gene lists (Step 4) and TF target genes are assembled into a regulatory matrix (Step 5) such that unrelated genes are assigned to an additional "baseline" class. Figure 1B: TF significance tests with multinomial logistic regression models. For a given TF, the alternative model H 1 is a univariate multinomial regression model that associates response (process genes) and one predictor (TF target genes), such that TF targets are linearly associated to probabilities of process gene classes (Step 6). The null model H 0 associates response (process genes) to their relative frequency in the dataset (Step 7). Log-likelihood ratio test measures if H 1 provides a better fit to data than the simpler H 0 model (Step 8). All TFs are subject to independent testing (Step 9) and subsequent multiple testing correction (Step 10). TF, transcription factor knockout strain; ChIP, chromatin immunoprecipitation; TF, transcription factor; TFBS, transcription factor binding site.
L(β|Y) = N i=1 C−1 c=1 ( p i,c p i,C ) y i,
where y i,c = 1 if y i is of class c and 0 otherwise, and the probability of gene-class relationship is computed as
p i,c = exp(β 0,c + K k=1 β k,c x i,k ) 1 + C−1 c=1 exp(β 0,c + K k=1 β k,c x i,k )

Maximising the log likelihood l leads to optimal regression coefficientsβ and the corresponding likelihood valueˆl :
(β,l) = arg maxβ β N i=1 C−1 c=1 (yi,c(β0 + K k=1 βk,cxi,k)) − yi,C log(1 + C−1 c=1 exp(β0,c + K k=1

Here we implemented a statistical test to assess the process specificity of a given TF by comparing two multinomial regression models. The null model H 0 : g(Y) = b 0 is an intercept-only model where process-specific genes are predicted solely based on their frequency in the full dataset (Step 7). The alternative model H 1 : g(Y) = b 0 + b k X k is a univariate model in which TF targets are also considered as predictors of process genes (Step 6). We use the likelihood ratio (LR) test with the chi-square distribution to compare the likelihoods of the two models, and decide if adding TF information substantially improves fit to data given its additional complexity (Step 8), as
P(H 0 ) = P χ 2 (−2(l(H 0 ) −l(H 1 )), ν 1 − ν 0 ),
where ν corresponds to degrees of freedom and reflects number of model parameters. To predict all regulators to a process of interest, we test all TFs independently, correct for multiple testing and find TFs with significant chi-square p-values (Benjamini-Yekutieli FDR, p ≤ 0.05).

In summary, m:Explorer uses the multinomial regression framework to associate process genes with TF regulatory targets from TFBS maps, gene expression patterns and nucleosome positioning data. Our method finds candidate TFs whose targets are especially informative of process genes, and thus may regulate their expression.

Yeast TF dataset with perturbation targets, DNA binding sites and nucleosome positioning

We used m:Explorer to study transcriptional regulation and TF function in yeast, as it has the widest collection of relevant genome-wide evidence. First we compiled a dataset of 285 regulators that contains carefully selected target genes for nearly all yeast TFs from microarrays, DNAbinding assays and nucleosome positioning measurements. Statistically significant target genes from regulator deletion experiments originate from our recent reanalysis [2] of an earlier study [1]. High-confidence TFBS targets were assembled from earlier chromatin immunoprecipitation (ChIP) assays by Harbison et al. [3], in silico TFBS predictions [4,17], and recent refinements with protein-binding microarrays by Zhu et al. [5]. The data were further processed with in vivo nucleosome positioning measurements [6] to distinguish binding sites where lower nucleosome occupancy reflects open chromatin structure.

Our dataset of 285 regulators contains 128,656 significant associations between regulators and target genes. Statistically reasoned cutoffs render our dataset sparse: it comprises high-confidence signals to 7.2% of approximately 1.8 million potential TF-gene pairs. The dataset includes 107 TF target sets with knockout data, 16 TFs with TFBS predictions and 162 TFs with both types of evidence. The majority of all gene-regulator associations (84%) originate from TF perturbation arrays ( Figure 2A). As observed previously, the agreement between binding sites and ΔTF targets is low: only 1.5% of all high-confidence targets constitute both types of evidence. Along with 170 confirmed or putative DNA-binding TFs, our dataset covers cofactors, chromatin modifiers and other regulatory proteins ( Figure 2B).

In conclusion, the yeast TF dataset is a useful resource for studying gene regulation.

High-confidence recovery of cell cycle regulators

First we tested m:Explorer in a well-defined biological context. Cell cycle is a thoroughly described regulatory system with four consecutive phases: gap-1 (G1), synthesis (S), gap-2 (G2) and mitosis (M). Some of the earliest microarray experiments identified cell cycle-regulated yeast genes [18,19], and a computational analysis organized these into phase-specific groups [20]. Several focused studies have investigated the roles of individual cell cycle TFs [21][22][23][24][25], and a genome-wide experiment outlined the underlying regulatory network in its interconnected, circular nature [26]. Altogether, the core cell cycle network comprises nine transcriptional regulators (Swi4, Swi6, Mbp1, Ndd1, Fkh1, Fkh2, Swi5, Ace2, Mcm1, Additional file 1, Table s1).

Here we applied m:Explorer and the TF dataset to select regulators to cell cycle genes. We focused on a recent tiling array study that measured genome-wide transcription during cell cycle at five minute resolution [27]. We used the list of 600 periodically expressed genes that contains specific groups for the four cell cycle phases and two checkpoints (G1, S, G2, G2/M, M, M/G1; 41-257 genes). This structured list of genes was then analyzed in a single m:Explorer run. We identified 46 statistically significant TFs (Benjamini-Yekutieli FDR p ≤ 0.05, LR test from m:Explorer) including all nine core TFs ( Figure 3A). Our results are ordered meaningfully, as eight of nine core TFs are ranked first (all p ≤ 10 −9 ). Besides core TFs, our results include at least four regulators that interact directly with the core TFs or act as secondary regulators. Notably, Stb1 forms a complex with G1/S TFs to affect gene expression in G1 [28], whereas Yox1 cooperates with Mcm1 to repress the expression of M/G1 specific genes [29]. The negative cell cycle regulator Ste12 is known to interact with Mcm1 in a specific pheromoneinduced response [30]. In addition to cell cycle regulators, we found components of the transcriptional machinery, including the general transcription factor Taf14 and multiple subunits of the Mediator complex (Ssn2, Cse2, Srb2, Srb8, Gal11). Several chromatin modifiers are also present, e.g. the silent information regulators (Sir2, Sir3) carry out genome silencing and are related to replicative cell ageing [31]. We expected to see such regulators among our predictions, since their disruption is likely to affect any process that involves transcription.

Our method reveals additional details about cell cycle regulation. First, as we model all cell cycle phases in one run, relative TF phase activities can be quantified through regression coefficients ( Figure 3B). For instance Swi4, Swi6 and Mbp1 make up the G1-S specific TF complexes MBF and SBF [21], and m:Explorer correctly highlights the phases with the strongest signal of regulatory activity. Second, we can assess the relative contribution of different kinds of regulatory evidence, and show that combined TFBS and ΔTF evidence are most informative of cell cycle regulation ( Figure 3C). Third, simultaneous analysis of multiple sub-processes in a single multinomial model is advantageous to separate logistic models for each related subprocess, since the latter approach is more prone to false positive predictions (Additional file 1, Figure s1). We performed m:Explorer analysis for four cell cycle phases and two checkpoints separately and recovered all cell cycle TFs found by the multinomial model, however also retrieved a large number (28) of additional false positive TFs not associated to cell cycle. Despite the above, analysis of sub-processes showed that m:Explorer is applicable to relatively small gene lists, for instance Mcm1 and Yox1 are correctly recovered as regulators of M-phase through only 55 informative genes.

Next we compared m:Explorer with eight similar methods for predicting TF function in regulatory networks (Additional file 1). As no other method allows exact replication of m:Explorer models, we used combinations of discretized and numeric gene expression, TF binding and cell cycle data as required (Table 1). Method performance evaluation was carried out with the Area Under Curve (AUC) statistic that accounted for 18 cell cycle TFs (Additional File 1, Tables s1-s2). To measure performance robustness, we also conducted  Figure 2A: Distribution of TF target genes from gene expression data (red, green), TF binding data (yellow, orange), and combined evidence (blue). Figure 2B: Distribution of regulator classes in the TF dataset. TF, transcription factor knockout strain; TF, transcription factor; TFBS, transcription factor binding site. a benchmark in which random subsets of input data were presented to each method (30, 50, 70 and 90% of yeast genes, 100 subsets each). The simulation shows that m:Explorer substantially outperforms all tested methods in recovering cell cycle regulators ( Figure 3D, AUC = 0.835 for 18 TFs, AUC = 0.996 for nine core TFs). Our method is reasonably accurate even when 50% of genes are discarded from the analysis (mean AUC = 0.747). The only method with comparable performance is the Fisher's exact test, a standard statistic for detecting significant biases in frequency tables. Comparison of m:Explorer and Fisher's test shows that our method is less prone to false positive discovery from randomly shuffled data (Additional File 1, Figure s2), and less dependent on microarray discretization parameters (Additional File 1, Figure s3). Fisher's test also prohibits the combined use of multiple features like gene expression, TF binding, nucleosome occupancy, and cell cycle phases. Simultaneous modeling of all data types in m:Explorer is likely to contribute to the demonstrated advantage over other approaches.

In conclusion, the cell cycle analysis showed that our approach successfully recovers a well-characterized regulatory system from multiple lines of high-throughput data. m:Explorer greatly outperformed several similar methods and showed robustness to incomplete data.  Figure 3A: Significance scores for 46 predicted TFs for cell cycle genes (FDR p ≤ 0.05). Red bars represent nine core cell cycle TFs and blue bars show secondary cell cycle TFs. Figure 3B: Predicted phase specificity for cell cycle TFs. Cell cycle phases are shown clockwise from top right, and TF activities across phases are shown in concentric circles (not drawn to scale). Intensity of red tone shows predicted proportion of TF activity from scaled regression coefficients. TF names are coloured according to their known phase of action, confirming the agreement between known and predicted functions of TFs. Figure 3C: Contribution of regulatory evidence in recovering 13 known cell cycle TFs. Bar color denotes class of evidence and bar height shows proportion of weight assigned to this class. Figure 3D: Performance and robustness comparison of m:Explorer and other methods. X-axis corresponds to proportion of genes presented to method, and Y-axis shows performance with the Area Under Curve statistic for 18 cell cycle TFs (AUC, 95% confidence intervals for 100 runs). m:Explorer values are plotted in wider, black lines. Numbers at the right end of curves reflect method performance with the full yeast dataset of 6253 genes. AUC, area under curve; G1, gap-1; G2, gap-2; M, mitosis; S, synthesis; TFs, transcription factors; TFBS, transcription factor binding site.

Computational prediction of TFs for quiescence entry and maintenance

Next, we applied m:Explorer in a less familiar biological context to create experimentally verifiable hypotheses about TF function. We focused on the transcriptional mechanisms that govern cell quiescence (G 0 , reviewed by Gray [32] and Kaeberlein [33]). G 0 is a cellular resting state with no proliferation, silenced genomes, reduced metabolism and translation, and greater stress resistance. Studying G 0 has proven difficult and related regulatory programs remain elusive.

Quiescence of yeast cells can be experimentally induced as a response to prolonged starvation ( Figure 4A). When glucose is depleted in exponentially growing cultures, growth rate is reduced as cells pass diauxic shift in which metabolic reprogramming initiates respiration of non-optimal carbon sources. Nutrients are depleted in post-diauxic phase, resulting in halted growth and differentiation to quiescent and non-quiescent cell populations [34]. The quiescent fraction of homogeneous cells may survive for extended periods of time, while the ageing heterogeneous non-quiescent fraction dies on further starvation. Consequently, culture viability starts decreasing rapidly in later stages of G 0 . Induction and inhibition of quiescence has been associated to several highly conserved signalling pathways, including protein kinases A and C (PKA, PKC), TOR and Snf1 [32].

Here we studied two public microarray datasets and executed m:Explorer in two independent rounds. First, we retrieved 207 diauxic shift genes in three distinct subgroups of early, transient and late expression from the dataset by Radonjic et al. [35]. Second, we used 594 genes and 676 genes characteristic of quiescent and non-quiescent cells from the study by Aragon et al. [36] ( Figure 4B). We identified 29 and 82 statistically significant candidate TFs in the two runs, log-transformed the scores and produced a final list of 97 G 0 regulators ( Figure 4C). A large number of regulators is expected, as G 0 entry is thought to comprise large-scale cellular reprogramming [35]. Several top-ranking TFs have high scores in both m:Explorer predictions. This ranking is not an artifact of the overlap between diauxic shift and quiescence genes. Although the two lists comprise a considerable number of common genes (n = 62, p = 0.0005, Fisher's exact test), these were not sufficient for predicting a similar collection of G 0 TFs, as m:Explorer analysis with the 62 genes only provided in a single significant TF (Mga2, p = 0.0005, LR test from m: Explorer).

In summary, the result of this analysis is an inclusive, prioritized list of candidate G 0 TFs that serves as a resource for hypothesis generation and experimental testing.

Experimental validation reveals super-wildtype and essential G 0 TFs Next we selected top 12 high-scoring TFs from our predictions for experimental testing. In total, 17 different strains were grown to G 0 and assessed for viability in six consecutive weekly measurements ( Figure 5A). We included deletion strains of candidate TFs (Δmga2, Δcst6, Δswi3, Δsds3, Δspt10, Δsin3, Δbas1, Δsnf2, Δspt20, Δhaa1, Δtup1, Δsnf11), positive controls (Δard1, Δmip1), negative controls (Δpdr3, Δgal3) and wildtype strains (Additional file 2). The viability of some strains was additionally monitored in five measurements over the first 72 hours of growth ( Figure 5B, Additional file 3). To confirm the timeframe of exponential growth and diauxic shift, we measured culture density and glucose levels of wildtype strains during 48 hours of growth (Additional file 1, Figure s4). To distinguish TFs with significant viability deviations, we used a linear error model that accounted for viability in wildtype and negative control strains as well as experimental batch effects.

All tested strains showed significant deviances from background viability at different stages of the quiescence time-course ( Figure 5A). The deletion strains of Bas1, Sds3, cst6, Mga2, and Spt10 show consistently greater viability in G 0 , indicating that their normal presence in wildtype cells suppresses viability and hastens cell ageing ( Figure 5C). We refer to these knockout phenotypes as super-wildtypes (WT+). In particular, Δbas1 strains are on average 1.7-4.5 times more viable than wildtype in weeks 3-6 of quiescence (all FDR p ≤ 10 −4 , LR test from error model). The transcription factor Bas1 is involved in the regulation of amino acid and nucleic acid metabolic pathways [37], and cst6 is related to chromosome stability and non-optimal carbon source regulation [38,39]. Spt10 and Sds3 are chromatin modifiers involved in genome silencing [40,41], and Mga2 regulates fatty acid metabolism, transcriptional silencing and response to low oxygen [42][43][44]. Deletion of Sds3 of the Sin3-Rpd3 histone deacetylase complex has been associated to increased chronological cell ageing [45]. The deletion strains Δtup1, Δswi3, Δhaa1 are significantly less viable than wildtype in quiescence ( Figure 5A). In particular, Δtup1 and Δswi3 strains become inviable in later stages of G 0 (viability ≤ 0.005) and can be considered essential for survival in this cell state ( Figure 5D). Two further strains Δspt20 and Δsnf2 are less viable in early quiescence, while Δsin3 shows later deviations. With the exceptions of Sin3 and Haa1, corresponding null mutants Figure 4 Quiescence (G 0 ) regulation and prediction of candidate TFs. Figure 4A: G 0 in S. cerevisiae occurs when a saturated culture is depleted of nutrients. As exponential growth stops with the exhaustion of glucose, cells pass diauxic shift and switch to slower respiratory growth. Growth is halted in post-diauxic phase, cells differentiate into quiescent and non-quiescent populations and enter G 0 . Culture viability starts decreasing rapidly as G 0 progresses. Two sets of process-specific genes were used in independent m:Explorer runs: genes expressed during diauxic shift (yellow box), and genes expressed in G 0 cultures (blue box, G 0 maintenance). Figure 4B: Venn diagram for cell cycle genes, diauxic shift genes, and G 0 maintenance genes. Statistically significant enrichment is observed between diauxic shift and quiescence genes (Fisher's exact test, p = 0.0005, 62 genes). Figure 4C: Log-scaled TF scores from diauxic shift and G 0 maintenance predictions (top -TFs 1-50, bottom -TFs 51-97). Bar color represents fraction of final score attributed to G 0 profile (diauxic shift -yellow; G 0 maintenance -blue). First 12 TFs highlighted with asterisks were selected for experimental testing. TF, transcription factor. are previously known for decreased or absent respiratory growth. Tup1 is a general inhibitor of transcription that establishes repressive chromatin structure [46]. Other factors are also involved in regulation of chromatin, transcription and genome stability, such as Swi3 and Snf2 of the SWI-SNF complex [47], Sin3 of Sin3-Rpd3 complex [48] and Spt20 of the SAGA complex [49]. While the factors have not been specifically described in the context of quiescence, disruption of their global functions is likely to affect this cellular state. Besides the above, the reduced G 0 viability of Δ haa1 potentially relates to its role in regulating cell wall proteins [50].

Before entering quiescence, most tested TF strains have similar viability to wildtype strains, suggesting that their function in regulating viability is specific to G 0 ( Figure 5B). During exponential growth at seven hours after inoculation, only three strains including the positive control Δ ard1 are significantly less viable. Ard1 encodes an N-terminal acetyltransferase subunit that guides genome silencing, and Δard1 fails to enter G 0 as observed previously [51]. In contrast, the other positive control Δmip1 is as viable as wildtype in exponential phase, and more viable in post-diauxic phase. Mip1 encodes a mitochondrial DNA polymerase subunit required for cell respiration [52], and Δmip1 loses viability in a similar manner to Δ tup1 ( Figure 5E). Curiously, Δspt10 is less viable in exponential growth phase and early quiescence, while its viability exceeds wildtype after week three of our time-course. The negative control strains Δgal3 and Δpdr3 expectedly show no major  Figure 5A: Viability of ΔTF strains in contrast to wildtype and negative controls, shown as a hierarchically clustered heat-map. Strains are shown vertically and time-points in days horizontally, coloured cells denote log fold change of viability (red -increased viability, blue -decreased viability). Asterisks denote statistical significance of fold change (linear regression, LR test, FDR p ≤ 0.05). Inviable strains with viability less than CFU = 0.005 are labeled with the symbol #. Time-point on day one corresponds to exponentially growing cells, measured seven hours after inoculation. Figure 5B: High-resolution viability curves of the first 72 hours of the time-course covering exponential growth, diauxic shift and post-diauxic phase. A subset of ΔTF strains (coloured solid lines), as well as wildtype strains (black solid lines) and controls (dashed and dotted lines) are shown. Figure 5C: G 0 viability curve of super-wildtype strain Δbas1. Figure 5D: Viability curve of G 0 -essential strain Δ tup1. Figure 5E: Viability curve of positive control strain Δmip1. Figure 5F: Viability curve of negative control strain Δtup1. TF, transcription factor knockout strain; FDR, false discovery rate.

deviations from wildtype viability ( Figure 5F). The TFs are related to alternative carbon metabolism and drug resistance, respectively [53,54], and show non-significant scores in m:Explorer predictions of G 0 TFs. Finally, our glycerol growth assays confirm the respiratory properties of tested strains (Additional file 1, Table s3) and mostly agree with previous studies [55,56]. However, in contrast to those reports, our data indicate that Δcst6 is viable on glycerol and indeed displays increased G 0 viability.

According to our knowledge, most of our predicted TFs are not recognized as quiescence regulators. However previous functional evidence refers to processes important in quiescence, and hence lends confidence to our experimental observations. Besides uncovering novel regulators of viability in G 0 , our experiments show that m:Explorer provides biologically meaningful prediction of regulator function.

Functional enrichment analysis explains roles of G 0 TFs

To gain insight into G 0 gene regulation of validated TFs, we performed a functional enrichment analysis of their G 0 target genes. We focused on quiescence genes defined by Aragon et al. [36] and identified the subset of genes that were bound by at least one WT+ TF or showed differential gene expression in at least one WT+ ΔTF microarray [2]. Target genes were then scored by product of differential expression p-values across all WT+ ΔTF microarrays and ranked such that genes with most dramatic transcriptional changes were prioritized. The target gene list for viability-deficient TF strains was complied in a similar fashion. We expect that ΔTF differential expression is informative of regulatory relationships in quiescence. The strains underlying microarray profiling are genetically identical to the strains in our G 0 experiments, although the former assays were performed with exponentially growing cells. Intersection of known quiescence genes with target genes of validated G 0 TFs, and subsequent prioritization according to differential expression, is therefore likely to highlight high-confidence TF targets and functional relationships. To investigate this in detail, we then used the ordered gene list analysis of g:Profiler [57] to study the functional importance of significanceranked target genes of WT+ and viability-deficient TFs.

Our analysis revealed 62 non-redundant Gene Ontology categories and KEGG and Reactome pathways with statistically significant enrichment in quiescence-related targets of G 0 TFs (FDR p ≤ 0.05, hypergeometric test, Figure 6). A number of functions were found to be enriched in TF targets corresponding to both viability phenotypes, suggesting that improved and reduced viability in quiescence may involve common regulatory pathways. The most significant results include the KEGG pathway of ribosome (p = 10 −15 ), proteolysis (p = 10 −11 ), reproduction (p = 10 −9 ) and oxidation-reduction process (p = 10 −10 ). Other functions are informative of TFs responsible for reduced G 0 viability. For instance, metabolic and catabolic genes (p = 0.0070 and p = 0.0035) are mostly up-regulated, while genes related to cell wall organization are inhibited (p = 0.030). In contrast, WT+ TFs with increased G 0 viability associate to down-regulation of protein metabolic genes (p = 10 −7 ) and modulation of alternative energy pathways such as fatty acid catabolism (p = 0.034) and glutamine metabolism (p = 0.047).

Taken together, the above results associate to known mechanisms of quiescence and provide clues of the regulatory programs of predicted G 0 TFs. Inhibition of growth through down-regulation of ribosome genes has been linked to increased replicative lifespan [58]. Efficient cell wall remodeling and response to increased oxidative stress are essential prerequisites of quiescence entry and survival [32]. Expectedly, increased viability appears to correlate with reduced metabolism, as related genes show opposite expression patterns in corresponding strains. Further discussion on G 0 TFs and related pathways can be found below.

Discussion

Function of G 0 regulators

It is tempting to speculate about the role of identified quiescence TFs in modulating quiescence signalling, as links between the factors and global G 0 -related pathways are apparent in our dataset. Our findings of WT+ regulators are especially intriguing, since their normal presence in wildtype cells reduces viability in quiescence and causes increased chronological ageing. From the perspective of evolutionary maintenance, WT+ regulators should engage in significant cellular functions that compensate for such negative properties.

As an example of G 0 regulation, protein kinase A (PKA) mediates nutritional signals to the cell and is known as an inhibitor of quiescence [32]. Its primary regulatory subunit Bcy1 acts as an inhibitor of the pathway, and mutations in Bcy1 cause viability loss and death in G 0 [59,60]. This double negative regulatory mechanism provides a potential explanation to observed viability phenotypes. In our TF dataset, Δmga2 has significantly higher levels of Bcy1, potentially allowing more starving cells to pass into quiescence. The G 0 -essential Tup1 and Swi3 knockout strains have depleted levels of Bcy1 and as a possible consequence, we observe reduction and loss of viability. As another example, protein kinase C (PKC) guides cell wall remodeling in response to starvation and its activity is required for G 0 entry [61]. The cell wall biosynthesis enzyme Gsc2 is a downstream target of PKC [62] and part of the gene expression signature of quiescent cells [36]. In ΔTF microarrays, Δmga2 and Δcst6 strains have elevated levels of Gsc2, while Δswi3 and Δtup1 show inhibition of PKC upstream of Gsc2 (Additional file 1, Figure s5).

Other genes with known function in G 0 appear to be regulated by WT+ and viability-deficient TFs. Notably, the conserved superoxide dismutase (SOD) genes are responsible for neutralizing oxidative damage of mitochondrial respiration. In yeast, SOD genes are required for G 0 survival and extend chronological lifespan when over-expressed [63,64]. Induced levels of Sod2 expression in Δcst6 may explain our observations of increased G 0 viability.

Several confirmed G 0 TFs are also associated to mammalian gene regulation.Cst6 carries the DNA-binding domain of CREB, an extensively studied TF that regulates a variety of processes, including cell survival and proliferation, cellular metabolism, and synaptic plasticity of long-term memory [65]. Bas1 is homologous to the MYB TF that regulates stem and progenitor cells and appears as an oncogene in multiple tumour types [66]. Chromatin modifier complexes Swi/Snf, Sin3/Rpd3 and SAGA are also broadly conserved, for instance Swi3 homolog SMARCC1 is involved in versatile functions, including neural stem cell renewal and differentiation [67]. As the yeast quiescence model associates to hallmark cancer properties of cell cycle control, proliferation and differentiation, further analysis of our Findings may reveal intriguing links to cancer biology.

Applicability and validity of m:Explorer

Here we present the robust computational method m: Explorer for predicting functions of gene regulators from high-throughput data. We applied a model that probabilistically accounts for multiple types of regulatory signals and Figure 6 Gene Ontology and pathway analysis of G 0 TFs.Quiescence genes with differential expression in G 0 TF knockout microarrays were studied with ranked enrichment analysis in g:Profiler. Statistically significant non-redundant functional Gene Ontology categories are shown (hypergeometric test, FDR p ≤ 0.05). Black bars correspond to total number of G 0 genes in a given category, and coloured bars show the number of times these genes were up-regulated (red) or down-regulated (green) in related G 0 ΔTFs, according to knockout data. Enrichment pvalue is shown in vertical colour strips on the left. Top: functional enrichments of quiescence genes associated to all G 0 TFs; middle: functional enrichments specific to TFs with reduced viability phenotype; bottom: Functional enrichments specific to WT+ TFs with increased viability phenotype. GO, Gene Ontology; TF, transcription factor. functional gene annotations. To take advantage of abundant genome-wide data and powerful experimental approaches, we present a case study for predicting transcription factors (TF) in the unicellular budding yeast. However, our method is not restricted to yeast and even not to these classes of data and regulators, being easily scalable to more complex regulatory systems of vertebrate organisms. Our method is also applicable to data such as protein-protein and genetic interactions that are categorical in nature. As shown here, m:Explorer is particularly useful in investigating sparse, high-confidence sets of data that may be controversial and not entirely comparable. For instance, we envisage large-scale characterization of human pathways in the context of heterogeneous tumours, utilizing sequence mutations, gene expression and chromatin modification data that are collected in cancer genomics projects.

In our model benchmarks, we demonstrate the advantage of univariate multinomial models in m:Explorer over similar multivariate models (AUC = 0.83 vs AUC = 0.67 for recovering cell cycle TFs). Briefly, the former models treat each TF independently in process gene classification, while the latter models include a non redundant collection of TFs as predictors. However, TF redundancy is an inherent property of robust biological networks that have evolved through gene and genome duplication [68]. In our case, the core cell cycle system involves three pairs of homologous TFs (Swi4 and Swi6; Fkh1 and Fkh2; Ace2 and Swi5) that have strikingly similar TFBS and expression patterns. Due to redundancy, such TFs are not treated as significant predictors in the multivariate setting. This is evident in our simulations: none of the tested multivariate models included both TFs of homologous pairs as significant predictors.

This analysis provides multiple lines of evidence to establish m:Explorer among other methods with similar goals. First, we carried out a highly detailed reconstruction of the known cell cycle regulatory system and proved the validity of our approach through existing knowledge. Second, we repeated the same analysis using eight alternative computational methods and random samples of input data, and provided quantitative proof to the robustness and better performance of our method. Third, we predicted regulators to the enigmatic cellular state of quiescence and validated our top-ranking candidate TFs in follow-up experiments. Nine of twelve tested TFs were confirmed to have consistent and significant G 0 viability deviations in gene knockout screens, while the remaining three factors showed differences in subsections of our time-course. Thus we proved a high success rate given our relatively simple experimental assays. Besides demonstrating the biological validity of our method, our findings reveal novel, previously unrecognized regulators of quiescence. m:Explorer web server and data availability m:Explorer is available as an R package on our web site [69] and elsewhere. The yeast TF dataset may prove to be a useful resource for the community and is also provided. We have established a web server at [69], allowing online prediction of regulator function using the yeast TF dataset.

Conclusions m:Explorer is a generally applicable method for inferring transcription factor function from heterogeneous highthroughput datasets. Our approach outperforms similar state-of-the art tools in recovering regulatory relationships in a well-studied eukaryotic system. Furthermore, the algorithm helps explore uncharacterized regulatory networks and propose valuable hypotheses for detailed assays. Our case study of quiescence G 0 and subsequent experimental validations revealed nine novel regulators that enhance or reduce cellular longevity, providing insights to investigators of this cryptic cellular state. In conclusion, our computational and experimental analyses provide strong support to the validity and usefulness of m:Explorer.

Materials and methods

Data processing

The yeast transcription factor dataset of 6253 genes and 285 transcription factors was compiled from gene expression, TF binding and nucleosome positioning data. Perturbation microarrays for 269 regulators were originally produced by Hu et al. [1], while our recently reanalyzed dataset [2] was used here for discretized, high-confidence values of up-and down-regulation (moderated t-test, FDR p ≤ 0.05). Further details on microarray preprocessing are available in the related publication [2]. TF binding site data for 178 TFs were compiled from multiple datasets of ChIP-chip [3], protein-binding microarrays [5] and computational predictions [4,17], using custom filtering and significance cutoffs proposed by the authors. Each promoter of 600 bp was considered to be bound by a TF if at least one binding site occurred in the dataset, and the TFBS was considered nucleosome-depleted (NDTFBS) if nucleosome occupancy [6] at the site was considerably below normalized genome-wide average (t-test, FDR p ≤ 0.05). Finally, gene expression and TF binding targets for each regulator were integrated and split into eight classes (up, down, TFBS, NDTFBS, up +TFBS, down+TFBS, up+NDTFBS, down+NDTFBS). All other genes except TF targets were assigned to the baseline class (not regulated).

Process-specific gene lists originate from previous highthroughput gene expression experiments. 600 cell cycle specific genes were retrieved from the tiling array experiment by Granovskaia et al. [27] and split into six sublists (G1, S, G2, G2/M, M, M/G1) according to authors' instructions. Three classes of diauxic shift genes (early, transient and late expression) originate from the G 0 time series [35], and genes specific to quiescent and nonquiescent cell cultures were first mapped in the analysis by Aragon et al. [36].

Computational methods

m:Explorer is based on univariate multinomial regression and implements the functionality of the R NNET package [70] for model fitting. We use a list of process-specific genes as categorical model response, and TF target genes as predictors. Briefly, m:Explorer compares two models: the null intercept-only model classifies process gene through their frequency, and the alternative univariate model additionally incorporates TF regulatory targets as predictors. We apply the log likelihood ratio test with null and alternative models to decide if TF target genes are significantly informative of process-related genes. Detailed description of the model is available in Additional file 1.

Yeast cell cycle TFs were predicted from a single structured gene list and directly ranked according to log pvalues from m:Explorer. G 0 TFs were predicted in two independent m:Explorer runs using genes from two datasets. TF p-values from LR tests were log-transformed, scaled to unit range and summed across the two runs to create unbiased composite scores for final ranking. Unitscaled positive regression coefficients were used to assess the relative phase specificity of cell cycle TFs, since these indicate over-represented regulatory targets in contrast to baseline genes. Relative contribution of regulatory evidence was computed in a similar way.

Linear regression was used to assess the significance of mutant strain viability deviations from control and wildtype strains. With viability as model response v, three types of variance were included as model predictors for assessing each mutant/time-point combination across all related replicas, as the alternative model H 1 : v~i + c + b + m. The above reflect global variance i, variance of negative controls c, variance between two batches of independent time-courses b, and additional variance of the tested strain m. Significance of viability deviation was assessed with a LR test, similarly to the m:Explorer algorithm. Specifically, the null model comprised only global variance, negative control variance and batch variance as H 0 : v~i + c + b, and null and alternative models were compared using the chi-square distribution. Resulting p-values were corrected for multiple testing with FDR.

Fisher's exact tests were used in multiple cases to evaluate the correlation of two binary variables. In the case of TF target genes and cell cycle genes, we applied the Fisher's test to assess whether the proportion of TF-regulated genes was statistically unexpected in the set of cell cycle genes. The Fisher's probability of observing a particular configuration in a two-way contingency table is computed as P(g CT , g Ct , g cT , g ct ) = g CT + g Ct g CT g cT + g ct g cT n g CT + g cT , where g denotes the number of genes in a particular set, C indicates cell cycle genes, T indicates TF targets, c shows genes unrelated to cell cycle, t shows genes not regulated by the particular TF, and n = g CT +g Ct +g cT +g ct reflects the number of all yeast genes. As Fisher's test does not support large contingency tables of multi-level variables, different types of TF regulatory targets were treated as the first category and non-regulated genes were assigned to second category, and cell cycle phasespecific genes were similarly merged into a bivariate discrete variable. A similar analysis was carried out to compare the overlap between diauxic shift genes and quiescence genes, using the set of all yeast genes as statistical background.

Gene Ontology (GO) and pathway enrichment analysis for G 0 TFs was carried out with with g:Profiler software [57]. We defined two ranked gene lists: G 0 genes [36] that were differentially expressed in WT+ TF knockout strains (Mga2, Cst6, Sds3, Spt10, Bas1), and G 0 genes that were differentially expressed in viability-deficient TF strains (Swi3, Sin3, Snf2, Spt20, Tup1, Haa1), according to TF knockout microarrays [2]. The gene lists were ordered according to statistical significance in TF knockout data [2], computed as products of p-values across WT+ and RD strains for every gene. We used the ordered enrichment analysis of g:Profiler to find GO functions and pathways in ranked gene lists and applied statistical filtering to find significant enrichments (FDR p ≤ 0.05).

The one-tailed hypergeometric tests calculated by g: Profiler assess the significance of observing k or more genes of a certain functional category in a list of n genes, as
P(x ≥ k) = n x=k K x N−K n−x N n ,
given that there are N genes in total and K of which are part of the functional category. As ordered enrichment analysis assumes that genes with stronger signals are ranked first, it consequently tests different subsets of the top list and returns the portion of top genes with the strongest p-value for a particular functional category [71]. Resulting G 0 functional categories were grouped into three classes: enriched G 0 categories associating to WT+ TF targets, categories of viability-deficient TF targets, and categories with statistical enrichment in both groups of targets. Enrichment p-values were corrected for multiple testing with the FDR procedure. To rank the third class of common functional categories, we multiplied corresponding p-values of WT+ target genes and viability-deficient TF target genes. After functional enrichment analysis, redundant categories whose genes formed a subset of some other category were removed. To quantify each GO category and function, we also counted up-regulated and down-regulated G 0 genes across all related TF strains.

Experimental procedures

Regulator knockout strains were selected as 12 top-ranking candidates from m:Explorer results. S. cerevisiae deletion strains originate from the EUROSCARF deletion collection in the BY4741 strain (MATa his3Δ1 leu2Δ0 met15Δ0 ura3Δ0). Liquid cultures were grown in triplicate at 30°C with aeration in YPD (1% yeast extract, 2% peptone, 2% glucose) for 28 days and subsequently shifted to room temperature without aeration. Viability measurements of the six-week time-course were taken in eight time-points: 7h after colony initiation, 48h after colony initiation, followed by six weekly measurements on days 7, 14, 21, 28, 35 and 42. Two independent batches involved distinct sets of tested strains, while wildtypes and controls were covered in both batches. A shorter, independent time-course covered the first three days of growth and involved viability measurements at 7h, 11h, 24h, 48h, and 72h. Cell density was measured at 600 nm. Colony forming units (CFU/ml) were determined by plating cells on YPD agar and counting colonies after three days of growth at 30°C. Culture viability was determined by dividing CFU/ml with total cell number per milliliter in corresponding culture (OD600 units ×10 7 ). Growth on glycerol was determined by streaking strains onto YPG plates (1% yeast extract, 2% peptone, 3% glycerol, 2% agar). Glucose concentration was determined by measuring NADPH production in hexokinase and glucose-6phosphate dehydrogenase coupled reactions provided by Roche.

Figure 1
1Method summary.

Figure 2
2Overview of yeast TF dataset.

Figure 3
3Cell cycle TF prediction.

Figure 5
5Experimental validation of G 0 TFs.

TF and TFBS target gene lists merged, 4 base classes and 4 combined classes of regulatory evidence used.c

p

y i,C

i,C ,

STEP 4
Pr
genes

TFBS1
TFBS2
TFBS3

TF1
TF3
TF2

TF1
TF3
TF2

STEP 1
TF knockout
microarrays

STEP 2
TF-DNA
interactions

STEP 3
Integrate
regulatory
data

+

Likelihood ratio test
2 (deviance(H1,H0))

H1 : p ~ TF 3
TF 3
redictor
of genes in p

H 0 : p ~ (intercept)
TF 3 predictive power
similar to random guess

P<0.05
P>0.05

(a) Data preprocessing and integration
(b) Regression analysis and

STEP 9: Repeat steps 6 and 8 for each TF

STEP 5
Compose genome-wide
regulatory matrix

One or more non-overlapping pr
gene lists (p1, p2, p3
by the user

no differentially
expressed genes

no binding sites

STEP 10: Multiple testing correction for tests

STEP 6: Fit model H1 with
process genes and TF targets

STEP 7: Fit model H0
with process genes

STEP 8: compar
1 and H 0

TF3 F

S1
TFBS2

TF1
TF2
TF3

All genes

r

g1
g2
g3
g4
g5
g6
g7
g8
g9
g10
g11
g12
g13
g14
g15
g16
g17
g18
g19
g20
g21
g22
g23

g1
g5

g7

g15

g10

g15

g10
g1

g7

g10

g5
g14
g3

g11
g21

g18
g19

g2

g15
g10

l
i

genes

g1 1
g2 2 2 2
g3 3 3 3
g4 4 4 4
g5 5 5 5
g6 6 6 6
g7 7 7 7
g8 8 8 8
g9 9 9 9
g10 0 0 0
g11 1
g12 2 2

p3
p2
p1

TF1

TF2
TF3

g5 g18 g19 g1 g7 g23 g10 g15

Induced and inhibited TF target genes from TF microarrays [P<0.05]

ChIP-chip and motif discovery data from multiple sources [P<0.001]
Nucleosome occupancy analysis shows accessible TFBS [P<0.05]

g18
g19

g2

g15
g10

g14
g3

g11
g21

Table 1
1Summary of comparison with similar methodsMethod
Software
Reference
ΔTF expression
TFBS
Gene function
Nucleosomes

Univariate multinomial regression
m:Explorer
Discretized
Discretized
Discretized
Discretized

Multivariate multinomial regression
Discretized
Discretized
Discretized
Discretized

Decision tree
GeneClass
[15]
Discretized
Discretized
Discretized

Kolmogorov Smirnov test
[10]
Discretized
Numeric 1

Mutual information
ARACNE
[11]
Numeric
Discretized 4

Fisher's exact test
Discretized 3
Discretized 3
Discretized
Discretized 3

Biclustering
SAMBA
[13]
Discretized
Discretized
Discretized 4
Discretized

Univariate linear regression
REDUCE
[14]
Numeric 2
Numeric 1
Multivariate linear regression
REDUCE
[14]
Numeric 2
Numeric 1

1. Gene expression values from cell cycle time-course. 2. Promoter sequence data and TF-DNA binding affinity matrices. 3. TF target genes treated as union of
gene expression and TFBS targets. 4. Fisher's test used for finding cell cycle TFs in constructed network or clustering. TF, transcription factor knockout strain;
TFBS, transcription factor binding site.

Reimand et al. Genome Biology 2012, 13:R55 http://genomebiology.com/2012/13/6/R55
Reimand et al. Genome Biology 2012, 13:R55 http://genomebiology.com/2012/13/6/R55 Page 7 of 16
Reimand et al. Genome Biology 2012, 13:R55 http://genomebiology.com/2012/13/6/R55 Page 10 of 16
Reimand et al. Genome Biology 2012, 13:R55 http://genomebiology.com/2012/13/6/R55 Page 11 of 16
Reimand et al. Genome Biology 2012, 13:R55 http://genomebiology.com/2012/13/6/R55 Page 16 of 16
AcknowledgementsWe are grateful to Gary Bader, Richard Bourgon, John Marioni, Gabriella Rustici, and Annabel Todd for their valuable advice, and to Tambet Arak for his work on the m:Explorer website. JR would like to acknowledge financial support from IKTDK, EXCS, ENFIN, EITSA, Marie Curie Biostar, the Ustus Agur Foundation, and the Artur Lind Foundation. Computational resources were provided by the University of Tartu HPC Centre.Additional materialAdditional file 2: Additional file 2 contains normalized colony forming unit measurements for tested ΔTF strains, wildtypes and controls of the six-week quiescence time-course.Additional file 3: Additional file 3 contains normalized colony forming unit measurements for tested ΔTF strains, wildtypes and controls of the 72-hour quiescence time-course. Abbreviations. ΔTFtranscription factor knockout strain; AUC -area under curve; CFU -colony forming units; ChIP -chromatin immunoprecipitation; DNA -deoxyribonucleic acid; FDR -false discovery rate; G 0 -stationary phase, quiescence; G1 -gap-1; G2 -gap-2; GO -gene ontology; KEGG -Kyoto encyclopedia of genes and genomes; LR -likelihood ratio; M -mitosis; MBF -mlu1-box binding factor; NADPH -nicotinamide adenine dinucleotide phosphate; NDTFBS -nucleosome-depleted transcription factor binding site; OD -optical density; PKA -protein kinase A; PKCprotein kinase C; S -synthesis; SAGA -Spt-Ada-Gcn5-acetyltransferase; SBF -SCB binding factor; TF -transcription factor; TFBS -transcription factor binding site; TOR -target of rapamycin; WT+ -super-wildtype; YPD -yeast extract peptone dextrose; YPG -yeast extract peptone glycerol.Authors' contributionsCompeting interestsThe authors declare that they have no competing interests.
Genetic reconstruction of a functional transcriptional regulatory network. Z Hu, P J Killion, V R Iyer, Nat Genet. 39Hu Z, Killion PJ, Iyer VR: Genetic reconstruction of a functional transcriptional regulatory network. Nat Genet 2007, 39:683-687.

Comprehensive reanalysis of transcription factor knockout expression data in Saccharomyces cerevisiae reveals many new targets. J Reimand, J M Vaquerizas, A E Todd, J Vilo, N M Luscombe, Nucleic Acids Res. 38Reimand J, Vaquerizas JM, Todd AE, Vilo J, Luscombe NM: Comprehensive reanalysis of transcription factor knockout expression data in Saccharomyces cerevisiae reveals many new targets. Nucleic Acids Res 2010, 38:4768-4777.

Transcriptional regulatory code of a eukaryotic genome. C T Harbison, D B Gordon, T I Lee, N J Rinaldi, K D Macisaac, T W Danford, N M Hannett, J B Tagne, D B Reynolds, J Yoo, E G Jennings, J Zeitlinger, D K Pokholok, M Kellis, P A Rolfe, K T Takusagawa, E S Lander, D K Gifford, E Fraenkel, R A Young, Nature. 431Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J, Jennings EG, Zeitlinger J, Pokholok DK, Kellis M, Rolfe PA, Takusagawa KT, Lander ES, Gifford DK, Fraenkel E, Young RA: Transcriptional regulatory code of a eukaryotic genome. Nature 2004, 431:99-104.

An improved map of conserved regulatory sites for Saccharomyces cerevisiae. K D Macisaac, T Wang, D B Gordon, D K Gifford, G D Stormo, E Fraenkel, BMC Bioinformatics. 7113MacIsaac KD, Wang T, Gordon DB, Gifford DK, Stormo GD, Fraenkel E: An improved map of conserved regulatory sites for Saccharomyces cerevisiae. BMC Bioinformatics 2006, 7:113.

High-resolution DNAbinding specificity analysis of yeast transcription factors. C Zhu, K J Byers, R P Mccord, Z Shi, M F Berger, D E Newburger, K Saulrieta, Z Smith, M V Shah, M Radhakrishnan, A A Philippakis, Y Hu, De Masi, F Pacek, M Rolfs, A Murthy, T Labaer, J Bulyk, M L , Genome Res. 19Zhu C, Byers KJ, McCord RP, Shi Z, Berger MF, Newburger DE, Saulrieta K, Smith Z, Shah MV, Radhakrishnan M, Philippakis AA, Hu Y, De Masi F, Pacek M, Rolfs A, Murthy T, Labaer J, Bulyk ML: High-resolution DNA- binding specificity analysis of yeast transcription factors. Genome Res 2009, 19:556-566.

The DNA-encoded nucleosome organization of a eukaryotic genome. N Kaplan, I K Moore, Y Fondufe-Mittendorf, A J Gossett, D Tillo, Y Field, E M Leproust, T R Hughes, J D Lieb, J Widom, E Segal, Nature. 458Kaplan N, Moore IK, Fondufe-Mittendorf Y, Gossett AJ, Tillo D, Field Y, LeProust EM, Hughes TR, Lieb JD, Widom J, Segal E: The DNA-encoded nucleosome organization of a eukaryotic genome. Nature 2009, 458:362-366.

Comparison of transcription regulatory interactions inferred from high-throughput methods: what do they reveal?. S Balaji, L M Iyer, M M Babu, L Aravind, Trends Genet. 24Balaji S, Iyer LM, Babu MM, Aravind L: Comparison of transcription regulatory interactions inferred from high-throughput methods: what do they reveal? Trends Genet 2008, 24:319-323.

. Reimand, Genome Biology. 1355Reimand et al. Genome Biology 2012, 13:R55

Cluster analysis and display of genome-wide expression patterns. M B Eisen, P T Spellman, P O Brown, D Botstein, Proc Natl Acad Sci. 95Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA 1998, 95:14863-14868.

Predicting gene regulatory elements in silico on a genomic scale. A Brazma, I Jonassen, J Vilo, E Ukkonen, Genome Res. 8Brazma A, Jonassen I, Vilo J, Ukkonen E: Predicting gene regulatory elements in silico on a genomic scale. Genome Res 1998, 8:1202-1215.

Statistical methods for identifying yeast cell cycle transcription factors. H K Tsai, H H Lu, W H Li, Proc Natl Acad Sci. 102Tsai HK, Lu HH, Li WH: Statistical methods for identifying yeast cell cycle transcription factors. Proc Natl Acad Sci USA 2005, 102:13532-13537.

ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. A A Margolin, I Nemenman, K Basso, C Wiggins, G Stolovitzky, Dalla Favera, R Califano, A , BMC Bioinformatics. 77Suppl 1Margolin AA, Nemenman I, Basso K, Wiggins C, Stolovitzky G, Dalla Favera R, Califano A: ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinformatics 2006, 7(Suppl 1):S7.

A human B-cell interactome identifies MYB and FOXM1 as master regulators of proliferation in germinal centers. C Lefebvre, P Rajbhandari, M J Alvarez, P Bandaru, W K Lim, M Sato, K Wang, P Sumazin, M Kustagi, B C Bisikirska, K Basso, P Beltrao, N Krogan, J Gautier, R Dalla-Favera, A Califano, Mol Syst Biol. 6377Lefebvre C, Rajbhandari P, Alvarez MJ, Bandaru P, Lim WK, Sato M, Wang K, Sumazin P, Kustagi M, Bisikirska BC, Basso K, Beltrao P, Krogan N, Gautier J, Dalla-Favera R, Califano A: A human B-cell interactome identifies MYB and FOXM1 as master regulators of proliferation in germinal centers. Mol Syst Biol 2010, 6:377.

Revealing modularity and organization in the yeast molecular network by integrated analysis of highly heterogeneous genomewide data. A Tanay, R Sharan, M Kupiec, R Shamir, Proc Natl Acad Sci. 1019Tanay A, Sharan R, Kupiec M, Shamir R: Revealing modularity and organization in the yeast molecular network by integrated analysis of highly heterogeneous genomewide data. Proc Natl Acad Sci USA 2004, 101(9):2981-2986.

Regulatory element detection using correlation with expression. H J Bussemaker, H Li, E D Siggia, Nat Genet. 27Bussemaker HJ, Li H, Siggia ED: Regulatory element detection using correlation with expression. Nat Genet 2001, 27:167-171.

Predicting genetic regulatory response using classification. M Middendorf, A Kundaje, C Wiggins, Y Freund, C Leslie, Bioinformatics. 201SupplMiddendorf M, Kundaje A, Wiggins C, Freund Y, Leslie C: Predicting genetic regulatory response using classification. Bioinformatics 2004, 20(Suppl 1):i232-240.

Non-transcriptional pathway features reconstructed from secondary effects of RNA interference. F Markowetz, J Bloch, R Spang, Bioinformatics. 21Markowetz F, Bloch J, Spang R: Non-transcriptional pathway features reconstructed from secondary effects of RNA interference. Bioinformatics 2005, 21:4026-4032.

Statistical features of yeast's transcriptional regulatory code. I Erb, E Van Nimwegen, IEEE Proceedings ICCSB. Erb I, van Nimwegen E: Statistical features of yeast's transcriptional regulatory code. IEEE Proceedings ICCSB 2006.

A genome-wide transcriptional analysis of the mitotic cell cycle. R J Cho, M J Campbell, E A Winzeler, L Steinmetz, A Conway, L Wodicka, T G Wolfsberg, A E Gabrielian, D Landsman, D J Lockhart, R W Davis, Mol Cell. 2Cho RJ, Campbell MJ, Winzeler EA, Steinmetz L, Conway A, Wodicka L, Wolfsberg TG, Gabrielian AE, Landsman D, Lockhart DJ, Davis RW: A genome-wide transcriptional analysis of the mitotic cell cycle. Mol Cell 1998, 2:65-73.

Comprehensive identification of cell cycleregulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. P T Spellman, G Sherlock, M Q Zhang, V R Iyer, K Anders, M B Eisen, P O Brown, D Botstein, B Futcher, Mol Biol Cell. 9Spellman PT, Sherlock G, Zhang MQ, Iyer VR, Anders K, Eisen MB, Brown PO, Botstein D, Futcher B: Comprehensive identification of cell cycle- regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. Mol Biol Cell 1998, 9:3273-3297.

Comparison of computational methods for the identification of cell cycle-regulated genes. U De Lichtenberg, L J Jensen, A Fausbøll, T S Jensen, P Bork, S Brunak, Bioinformatics. 21de Lichtenberg U, Jensen LJ, Fausbøll A, Jensen TS, Bork P, Brunak S: Comparison of computational methods for the identification of cell cycle-regulated genes. Bioinformatics 2005, 21:1164-1171.

A role for the transcription factors Mbp1 and Swi4 in progression from G1 to S phase. C Koch, T Moll, M Neuberg, H Ahorn, K Nasmyth, Science. 261Koch C, Moll T, Neuberg M, Ahorn H, Nasmyth K: A role for the transcription factors Mbp1 and Swi4 in progression from G1 to S phase. Science 1993, 261:1551-1557.

Forkhead transcription factors, Fkh1p and Fkh2p, collaborate with Mcm1p to control transcription required for M-phase. R Kumar, D M Reynolds, A Shevchenko, A Shevchenko, S D Goldstone, S Dalton, Curr Biol. 10Kumar R, Reynolds DM, Shevchenko A, Shevchenko A, Goldstone SD, Dalton S: Forkhead transcription factors, Fkh1p and Fkh2p, collaborate with Mcm1p to control transcription required for M-phase. Curr Biol 2000, 10:896-906.

Distinct regions of the Swi5 and Ace2 transcription factors are required for specific gene activation. H J Mcbride, Y Yu, D J Stillman, J Biol Chem. 274McBride HJ, Yu Y, Stillman DJ: Distinct regions of the Swi5 and Ace2 transcription factors are required for specific gene activation. J Biol Chem 1999, 274:21029-21036.

NDD1, a high-dosage suppressor of cdc28-1N, is essential for expression of a subset of late-S-phase-specific genes in Saccharomyces cerevisiae. C J Loy, D Lydall, U Surana, Mol Cell Biol. 19Loy CJ, Lydall D, Surana U: NDD1, a high-dosage suppressor of cdc28-1N, is essential for expression of a subset of late-S-phase-specific genes in Saccharomyces cerevisiae. Mol Cell Biol 1999, 19:3312-3327.

Mcm1 is required to coordinate G2-specific transcription in Saccharomyces cerevisiae. H Althoefer, A Schleiffer, K Wassmann, A Nordheim, G Ammerer, Mol Cell Biol. 15Althoefer H, Schleiffer A, Wassmann K, Nordheim A, Ammerer G: Mcm1 is required to coordinate G2-specific transcription in Saccharomyces cerevisiae. Mol Cell Biol 1995, 15:5917-5928.

Serial regulation of transcriptional regulators in the yeast cell cycle. I Simon, J Barnett, N Hannett, C T Harbison, N J Rinaldi, T L Volkert, J J Wyrick, J Zeitlinger, D K Gifford, T S Jaakkola, R A Young, Cell. 106Simon I, Barnett J, Hannett N, Harbison CT, Rinaldi NJ, Volkert TL, Wyrick JJ, Zeitlinger J, Gifford DK, Jaakkola TS, Young RA: Serial regulation of transcriptional regulators in the yeast cell cycle. Cell 2001, 106:697-708.

High-resolution transcription atlas of the mitotic cell cycle in budding yeast. M V Granovskaia, L J Jensen, M E Ritchie, J Toedling, Y Ning, P Bork, W Huber, L M Steinmetz, Genome Biol. 1124Granovskaia MV, Jensen LJ, Ritchie ME, Toedling J, Ning Y, Bork P, Huber W, Steinmetz LM: High-resolution transcription atlas of the mitotic cell cycle in budding yeast. Genome Biol 2010, 11:R24.

G1 transcription factors are differentially regulated in Saccharomyces cerevisiae by the Swi6-binding protein Stb1. M Costanzo, O Schub, B Andrews, Mol Cell Biol. 23Costanzo M, Schub O, Andrews B: G1 transcription factors are differentially regulated in Saccharomyces cerevisiae by the Swi6-binding protein Stb1. Mol Cell Biol 2003, 23:5064-5077.

Conserved homeodomain proteins interact with MADS box protein Mcm1 to restrict ECB-dependent transcription to the M/G1 phase of the cell cycle. T Pramila, S Miles, D Guhathakurta, D Jemiolo, L L Breeden, Genes Dev. 16Pramila T, Miles S, GuhaThakurta D, Jemiolo D, Breeden LL: Conserved homeodomain proteins interact with MADS box protein Mcm1 to restrict ECB-dependent transcription to the M/G1 phase of the cell cycle. Genes Dev 2002, 16:3034-3045.

Ste12 and Mcm1 regulate cell cycledependent transcription of FAR1. L J Oehlen, J D Mckinney, F R Cross, Mol Cell Biol. 16Oehlen LJ, McKinney JD, Cross FR: Ste12 and Mcm1 regulate cell cycle- dependent transcription of FAR1. Mol Cell Biol 1996, 16:2830-2837.

The SIR2/3/4 complex and SIR2 alone promote longevity in Saccharomyces cerevisiae by two different mechanisms. M Kaeberlein, M Mcvey, L Guarente, Genes Dev. 13Kaeberlein M, McVey M, Guarente L: The SIR2/3/4 complex and SIR2 alone promote longevity in Saccharomyces cerevisiae by two different mechanisms. Genes Dev 1999, 13:2570-2580.

J V Gray, G A Petsko, G C Johnston, D Ringe, R A Singer, M Werner-Washburne, Sleeping beauty":quiescence in Saccharomyces cerevisiae. 68Gray JV, Petsko GA, Johnston GC, Ringe D, Singer RA, Werner-Washburne M: Sleeping beauty":quiescence in Saccharomyces cerevisiae. Microbiol Mol Biol Rev 2004, 68:187-206.

Recent developments in yeast aging. M Kaeberlein, C R Burtner, B K Kennedy, PLoS Genet. 384Kaeberlein M, Burtner CR, Kennedy BK: Recent developments in yeast aging. PLoS Genet 2007, 3:e84.

Isolation of quiescent and nonquiescent cells from yeast stationary-phase cultures. C Allen, S Buttner, A D Aragon, J A Thomas, O Meirelles, J E Jaetao, D Benn, S W Ruby, M Veenhuis, F Madeo, M Werner-Washburne, J Cell Biol. 174Allen C, Buttner S, Aragon AD, Thomas JA, Meirelles O, Jaetao JE, Benn D, Ruby SW, Veenhuis M, Madeo F, Werner-Washburne M: Isolation of quiescent and nonquiescent cells from yeast stationary-phase cultures. J Cell Biol 2006, 174:89-100.

Genome-wide analyses reveal RNA polymerase II located upstream of genes poised for rapid response upon S. cerevisiae stationary phase exit. M Radonjic, J C Andrau, P Lijnzaad, P Kemmeren, T T Kockelkorn, D Van Leenen, N L Van Berkum, F C Holstege, Mol Cell. 18Radonjic M, Andrau JC, Lijnzaad P, Kemmeren P, Kockelkorn TT, van Leenen D, van Berkum NL, Holstege FC: Genome-wide analyses reveal RNA polymerase II located upstream of genes poised for rapid response upon S. cerevisiae stationary phase exit. Mol Cell 2005, 18:171-183.

Characterization of differentiated quiescent and nonquiescent cells in yeast stationaryphase cultures. A D Aragon, A L Rodriguez, O Meirelles, S Roy, G S Davidson, P H Tapia, C Allen, R Joe, D Benn, M Werner-Washburne, Mol Biol Cell. 19Aragon AD, Rodriguez AL, Meirelles O, Roy S, Davidson GS, Tapia PH, Allen C, Joe R, Benn D, Werner-Washburne M: Characterization of differentiated quiescent and nonquiescent cells in yeast stationary- phase cultures. Mol Biol Cell 2008, 19:1271-1280.

Role of the myb-like protein bas1p in Saccharomyces cerevisiae: a proteome analysis. V Denis, H Boucherie, C Monribot, B Daignan-Fornier, Mol Microbiol. 30Denis V, Boucherie H, Monribot C, Daignan-Fornier B: Role of the myb-like protein bas1p in Saccharomyces cerevisiae: a proteome analysis. Mol Microbiol 1998, 30:557-566.

New yeast genes important for chromosome integrity and segregation identified by dosage effects on genome stability. I I Ouspenski, S J Elledge, B R Brinkley, Nucleic Acids Res. 27Ouspenski II, Elledge SJ, Brinkley BR: New yeast genes important for chromosome integrity and segregation identified by dosage effects on genome stability. Nucleic Acids Res 1999, 27:3001-3008.

Aca1 and Aca2, ATF/CREB activators in Saccharomyces cerevisiae, are important for carbon source utilization but not the response to stress. M A Garcia-Gimeno, K Struhl, Mol Cell Biol. 20Garcia-Gimeno MA, Struhl K: Aca1 and Aca2, ATF/CREB activators in Saccharomyces cerevisiae, are important for carbon source utilization but not the response to stress. Mol Cell Biol 2000, 20:4340-4349.

Global regulation by the yeast Spt10 protein is mediated through chromatin structure and the histone upstream activating sequence elements. P R Riksson, G Mendiratta, N B Mclaughlin, T G Wolfsberg, L Marino-Ramirez, T A Pompa, M Jainerin, D Landsman, C H Shen, D J Clark, Mol Cell Biol. 25riksson PR, Mendiratta G, McLaughlin NB, Wolfsberg TG, Marino-Ramirez L, Pompa TA, Jainerin M, Landsman D, Shen CH, Clark DJ: Global regulation by the yeast Spt10 protein is mediated through chromatin structure and the histone upstream activating sequence elements. Mol Cell Biol 2005, 25:9127-9137.

Evidence that the transcriptional regulators SIN3 and RPD3, and a novel gene (SDS3) with similar functions, are involved in transcriptional silencing in S. cerevisiae. D Vannier, D Balderes, D Shore, Genetics. 144Vannier D, Balderes D, Shore D: Evidence that the transcriptional regulators SIN3 and RPD3, and a novel gene (SDS3) with similar functions, are involved in transcriptional silencing in S. cerevisiae. Genetics 1996, 144:1343-1353.

MGA2 and SPT23 are modifiers of transcriptional silencing in yeast. M L Dula, S G Holmes, Genetics. 156Dula ML, Holmes SG: MGA2 and SPT23 are modifiers of transcriptional silencing in yeast. Genetics 2000, 156:933-941.

The membrane proteins, Spt23p and Mga2p, play distinct roles in the activation of Saccharomyces cerevisiae OLE1 gene expression. Fatty acid-mediated regulation of Mga2p activity is independent of its proteolytic processing into a soluble transcription activator. R Chellappa, P Kandasamy, C S Oh, Y Jiang, M Vemula, C E Martin, J Biol Chem. 276Chellappa R, Kandasamy P, Oh CS, Jiang Y, Vemula M, Martin CE: The membrane proteins, Spt23p and Mga2p, play distinct roles in the activation of Saccharomyces cerevisiae OLE1 gene expression. Fatty acid-mediated regulation of Mga2p activity is independent of its proteolytic processing into a soluble transcription activator. J Biol Chem 2001, 276:43548-43556.

MGA2 is involved in the low-oxygen response element-dependent hypoxic induction of genes in Saccharomyces cerevisiae. Y Jiang, M J Vasconcelles, S Wretzel, A Light, C E Martin, M A Goldberg, Mol Cell Biol. 21Jiang Y, Vasconcelles MJ, Wretzel S, Light A, Martin CE, Goldberg MA: MGA2 is involved in the low-oxygen response element-dependent hypoxic induction of genes in Saccharomyces cerevisiae. Mol Cell Biol 2001, 21:6161-6169.

Induction of autophagy by spermidine promotes longevity. T Eisenberg, H Knauer, A Schauer, S Buttner, C Ruckenstuhl, D Carmona-Gutierrez, J Ring, S Schroeder, C Magnes, L Antonacci, H Fussi, L Deszcz, R Hartl, E Schraml, A Criollo, E Megalou, D Weiskopf, P Laun, G Heeren, M Breitenbach, B Grubeck-Loebenstein, E Herker, B Fahrenkrog, K U Frohlich, F Sinner, N Tavernarakis, N Minois, G Kroemer, F Madeo, Nat Cell Biol. 11Eisenberg T, Knauer H, Schauer A, Buttner S, Ruckenstuhl C, Carmona- Gutierrez D, Ring J, Schroeder S, Magnes C, Antonacci L, Fussi H, Deszcz L, Hartl R, Schraml E, Criollo A, Megalou E, Weiskopf D, Laun P, Heeren G, Breitenbach M, Grubeck-Loebenstein B, Herker E, Fahrenkrog B, Frohlich KU, Sinner F, Tavernarakis N, Minois N, Kroemer G, Madeo F: Induction of autophagy by spermidine promotes longevity. Nat Cell Biol 2009, 11:1305-1314.

Transcriptional repression: the long and the short of it. A J Courey, S Jia, Genes Dev. 15Courey AJ, Jia S: Transcriptional repression: the long and the short of it. Genes Dev 2001, 15:2786-2796.

The Swi/Snf family nucleosome-remodeling complexes and transcriptional control. P Sudarsanam, F Winston, Trends Genet. 16Sudarsanam P, Winston F: The Swi/Snf family nucleosome-remodeling complexes and transcriptional control. Trends Genet 2000, 16:345-351.

Sin3: a flexible regulator of global gene expression and genome stability. R A Silverstein, K Ekwall, Curr Genet. 47Silverstein RA, Ekwall K: Sin3: a flexible regulator of global gene expression and genome stability. Curr Genet 2005, 47:1-17.

Functional organization of the yeast SAGA complex: distinct components involved in structural integrity, nucleosome acetylation, and TATA-binding protein interaction. D E Sterner, P A Grant, S M Roberts, L J Duggan, R Belotserkovskaya, L A Pacella, F Winston, J L Workman, S L Berger, Mol Cell Biol. 19Sterner DE, Grant PA, Roberts SM, Duggan LJ, Belotserkovskaya R, Pacella LA, Winston F, Workman JL, Berger SL: Functional organization of the yeast SAGA complex: distinct components involved in structural integrity, nucleosome acetylation, and TATA-binding protein interaction. Mol Cell Biol 1999, 19:86-98.

. Reimand, Genome Biology. 1355Reimand et al. Genome Biology 2012, 13:R55

Sa-Correia I: The SPI1 gene, encoding a glycosylphosphatidylinositol-anchored cell wall protein, plays a prominent role in the development of yeast resistance to lipophilic weak-acid food preservatives. T Simoes, N P Mira, A R Fernandes, Appl Environ Microbiol. 72Simoes T, Mira NP, Fernandes AR, Sa-Correia I: The SPI1 gene, encoding a glycosylphosphatidylinositol-anchored cell wall protein, plays a prominent role in the development of yeast resistance to lipophilic weak-acid food preservatives. Appl Environ Microbiol 2006, 72:7168-7175.

Identification and characterization of genes and mutants for an N-terminal acetyltransferase from yeast. J R Mullen, P S Kayne, R P Moerschell, S Tsunasawa, M Gribskov, M Colavito-Shepanski, M Grunstein, F Sherman, R Sternglanz, EMBO J. 8Mullen JR, Kayne PS, Moerschell RP, Tsunasawa S, Gribskov M, Colavito- Shepanski M, Grunstein M, Sherman F, Sternglanz R: Identification and characterization of genes and mutants for an N-terminal acetyltransferase from yeast. EMBO J 1989, 8:2067-2075.

Cloning and sequencing of the nuclear gene MIP1 encoding the catalytic subunit of the yeast mitochondrial DNA polymerase. F Foury, J Biol Chem. 264Foury F: Cloning and sequencing of the nuclear gene MIP1 encoding the catalytic subunit of the yeast mitochondrial DNA polymerase. J Biol Chem 1989, 264:20552-20560.

Genetic and molecular analysis of the GAL3 gene in the expression of the galactose/melibiose regulon of Saccharomyces cerevisiae. T E Torchia, J E Hopper, Genetics. 113Torchia TE, Hopper JE: Genetic and molecular analysis of the GAL3 gene in the expression of the galactose/melibiose regulon of Saccharomyces cerevisiae. Genetics 1986, 113:229-246.

PDR3, a new yeast regulatory gene, is homologous to PDR1 and controls the multidrug resistance phenomenon. T Delaveau, A Delahodde, E Carvajal, J Subik, C Jacq, Mol Gen Genet. 244Delaveau T, Delahodde A, Carvajal E, Subik J, Jacq C: PDR3, a new yeast regulatory gene, is homologous to PDR1 and controls the multidrug resistance phenomenon. Mol Gen Genet 1994, 244:501-511.

Systematic screen for human disease genes in yeast. L M Steinmetz, C Scharfe, A M Deutschbauer, D Mokranjac, Z S Herman, T Jones, A M Chu, G Giaever, H Prokisch, P J Oefner, R W Davis, Nat Genet. 31Steinmetz LM, Scharfe C, Deutschbauer AM, Mokranjac D, Herman ZS, Jones T, Chu AM, Giaever G, Prokisch H, Oefner PJ, Davis RW: Systematic screen for human disease genes in yeast. Nat Genet 2002, 31:400-404.

Genetic basis of mitochondrial function and morphology in Saccharomyces cerevisiae. K S Dimmer, S Fritz, F Fuchs, M Messerschmitt, N Weinbach, W Neupert, B Westermann, Mol Biol Cell. 13Dimmer KS, Fritz S, Fuchs F, Messerschmitt M, Weinbach N, Neupert W, Westermann B: Genetic basis of mitochondrial function and morphology in Saccharomyces cerevisiae. Mol Biol Cell 2002, 13:847-853.

Profiler-a web server for functional interpretation of gene lists. J Reimand, T Arak, J Vilo, Nucleic Acids Res. 39Reimand J, Arak T, Vilo J: g:Profiler-a web server for functional interpretation of gene lists (2011 update). Nucleic Acids Res 2011, 39: W307-315.

Regulation of yeast replicative life span by TOR and Sch9 in response to nutrients. M Kaeberlein, R W Powers, K K Steffen, E A Westman, D Hu, N Dang, E O Kerr, K T Kirkland, S Fields, B K Kennedy, Science. 310Kaeberlein M, Powers RW, Steffen KK, Westman EA, Hu D, Dang N, Kerr EO, Kirkland KT, Fields S, Kennedy BK: Regulation of yeast replicative life span by TOR and Sch9 in response to nutrients. Science 2005, 310:1193-1196.

Glucose signaling in Saccharomyces cerevisiae. G M Santangelo, Microbiol Mol Biol Rev. 70Santangelo GM: Glucose signaling in Saccharomyces cerevisiae. Microbiol Mol Biol Rev 2006, 70:253-282.

Yeast bcy1 mutants with stationary phase-specific defects. V M Peck, E K Fuge, P A Padilla, M A Gomez, M Werner-Washburne, Curr Genet. 32Peck VM, Fuge EK, Padilla PA, Gomez MA, Werner-Washburne M: Yeast bcy1 mutants with stationary phase-specific defects. Curr Genet 1997, 32:83-92.

The protein kinase C pathway is required for viability in quiescence in Saccharomyces cerevisiae. S A Krause, J V Gray, Curr Biol. 12Krause SA, Gray JV: The protein kinase C pathway is required for viability in quiescence in Saccharomyces cerevisiae. Curr Biol 2002, 12:588-593.

Temperature-induced expression of yeast FKS2 is under the dual control of protein kinase C and calcineurin. C Zhao, U S Jung, P Garrett-Engele, T Roe, M S Cyert, D E Levin, Mol Cell Biol. 18Zhao C, Jung US, Garrett-Engele P, Roe T, Cyert MS, Levin DE: Temperature-induced expression of yeast FKS2 is under the dual control of protein kinase C and calcineurin. Mol Cell Biol 1998, 18:1013-1022.

Superoxide dismutase activity is essential for stationary phase survival in Saccharomyces cerevisiae. Mitochondrial production of toxic oxygen species in vivo. V D Longo, E B Gralla, J S Valentine, J Biol Chem. 271Longo VD, Gralla EB, Valentine JS: Superoxide dismutase activity is essential for stationary phase survival in Saccharomyces cerevisiae. Mitochondrial production of toxic oxygen species in vivo. J Biol Chem 1996, 271:12275-12280.

P Fabrizio, L L Liou, V N Moy, A Diaspro, J S Valentine, E B Gralla, V D Longo, SOD2 functions downstream of Sch9 to extend longevity in yeast. 163Fabrizio P, Liou LL, Moy VN, Diaspro A, Valentine JS, Gralla EB, Longo VD: SOD2 functions downstream of Sch9 to extend longevity in yeast. Genetics 2003, 163:35-46.

Transcriptional regulation by the phosphorylationdependent factor CREB. B Mayr, M Montminy, Nat Rev Mol Cell Biol. 2Mayr B, Montminy M: Transcriptional regulation by the phosphorylation- dependent factor CREB. Nat Rev Mol Cell Biol 2001, 2:599-609.

MYB function in normal and cancer cells. R G Ramsay, T J Gonda, Nat Rev Cancer. 8Ramsay RG, Gonda TJ: MYB function in normal and cancer cells. Nat Rev Cancer 2008, 8:523-534.

An essential switch in subunit composition of a chromatin remodeling complex during neural development. J Lessard, J I Wu, J A Ranish, M Wan, M M Winslow, B T Staahl, H Wu, R Aebersold, I A Graef, G R Crabtree, Neuron. 55Lessard J, Wu JI, Ranish JA, Wan M, Winslow MM, Staahl BT, Wu H, Aebersold R, Graef IA, Crabtree GR: An essential switch in subunit composition of a chromatin remodeling complex during neural development. Neuron 2007, 55:201-215.

H: Biological Kitano, Robustness, 69. m:Explorer website. 5Kitano H: Biological robustness. Nat Rev Genet 2004, 5:826-837. 69. m:Explorer website. [http://biit.cs.ut.ee/mexplorer/].

W N Venables, B D Ripley, Modern Applied Statistics with S. 4 edition. New York, USASpringerVenables WN, Ripley BD: Modern Applied Statistics with S. 4 edition. New York, USA: Springer; 2003.

Profiler-a web-based toolset for functional profiling of gene lists from large-scale experiments. J Reimand, M Kull, H Peterson, J Hansen, J Vilo, 10.1186/gb-2012-13-6-r55Nucleic Acids Res. 35Reimand J, Kull M, Peterson H, Hansen J, Vilo J: g:Profiler-a web-based toolset for functional profiling of gene lists from large-scale experiments. Nucleic Acids Res 2007, 35:193-200. doi:10.1186/gb-2012-13-6-r55

m:Explorer: multinomial regression models reveal positive and negative regulators of longevity in yeast quiescence. Reimand, Genome Biology. 1355Cite this article asCite this article as: Reimand et al.: m:Explorer: multinomial regression models reveal positive and negative regulators of longevity in yeast quiescence. Genome Biology 2012 13:R55.