Abstract
Budding yeast cells are assumed to trigger Start and enter the cell cycle only after they attain a critical size set by external conditions. However, arguing against deterministic models of cell size control, cell volume at Start displays great individual variability even under constant conditions. Here we show that cell size at Start is robustly set at a single-cell level by the volume growth rate in G1, which explains the observed variability. We find that this growth-rate-dependent sizer is intimately hardwired into the Start network and the Ydj1 chaperone is key for setting cell size as a function of the individual growth rate. Mathematical modelling and experimental data indicate that a growth-rate-dependent sizer is sufficient to ensure size homeostasis and, as a remarkable advantage over a rigid sizer mechanism, it reduces noise in G1 length and provides an immediate solution for size adaptation to external conditions at a population level.
Similar content being viewed by others
Introduction
Cell size can be considered as the outcome of complex parallel and interconnected processes that convey different sorts of intrinsic and extrinsic information. Of all these processes, cell growth and cell cycle are particularly relevant as they likely apply to all organisms, and their efficient and precise coordination is a fundamental biological problem1,2,3,4,5,6,7,8,9,10. As a unicellular eukaryote, Saccharomyces cerevisiae has been widely used to study the mechanisms that determine cell size, mainly because this budding yeast displays two of the most universal properties regarding cell size control, a constant mass/ploidy ratio11 and a critical size threshold for cell cycle progression12,13, which is mainly exerted at Start, a key event in late G1. Two heteromeric transcription factors14, SBF (Swi6–Swi4) and MBF (Mbp1–Swi6), drive a transcriptional wave of ~200 genes15 at Start to initiate budding and trigger S-phase entry. The most upstream activator of Start is a complex formed by Cln316,17,18,19,20, a G1 cyclin, and Cdc2821, the cyclin-dependent kinase that controls the cell cycle in budding yeast. The Cdc28–Cln3 complex phosphorylates Whi5 (and presumably Swi6) at multiple residues to activate SBF- and MBF-dependent transcription22,23 in a scheme homologous to the onset of the mammalian cell cycle, where Cdk4,6–cyclin D complexes phosphorylate retinoblastoma protein (pRB) to activate E2F–DP transcription factors24. Two other G1 cyclins, Cln1 and Cln2, are upregulated by SBF and have a key role with Cdc28 in a positive feedback loop essential for coherent and irreversible execution of Start25,26,27.
The Cln3 cyclin is present at constant levels throughout G119,28 and it is regulated by a retention/release mechanism that prevents its unscheduled accumulation in the nucleus29. Whi3, first identified as a negative regulator of Cln3 activity30, binds the CLN3 mRNA and recruits Cdc28 to help retain newly formed Cdc28–Cln3 complexes at the endoplasmic reticulum (ER) during early and mid G131,32,33. In late G1, the Ydj1 chaperone has an essential role in releasing Cln3 or the Cdc28–Cln3 complex from the ER33, which allows its accumulation in the nucleus to phosphorylate Whi5 and trigger Start.
Screens for small-cell-size mutants have not only identified key proteins of the Start network, such as Cln3, Whi3 and Whi5, but also central regulators of ribosome biogenesis34, namely Sfp1 and Sch9, indicating a close functional interaction between the machineries that control growth and proliferation4. On the other hand, and in addition to its role at Start, Ydj1 is heavily involved in protein translocation across the ER membrane35, a key process for cell growth. These lines of evidence suggest that the critical size could be set by mechanisms intimately linked to cell growth4,29.
Although it is generally accepted that budding yeast cells must grow to reach a critical size before cell cycle entry4,10,12,13, cell volume at Start displays a very large variability under constant conditions36, which strongly argues against the existence of a deterministic mode of cell size control. This key question has been analysed recently at a single-cell level, and it was proposed that cell size variability at Start is mainly due to molecular noise intrinsic to the mechanisms that activate the G1/S regulon37. Here we show that most of the large variation in cell size at Start is not because of an inherent stochastic behaviour. Instead, cell size at Start is an individual parameter robustly set by the rate at which each individual cell grows during G1. We also show that the growth-rate-dependent sizer is intimately hardwired into the Start network, the Ydj1 chaperone being the most important factor for setting the critical size as a function of the individual growth potential. We also show that a growth-rate-dependent sizer ensures size homeostasis and gives a simple solution to cell size adaptation to different external conditions as long as they have an impact on growth rate. Moreover, as a notable advantage over the constant critical size postulated up to now, a growth-rate-dependent sizer reduces noise in G1 length and equalizes the generation time of fast- and slow-growing cells, thus maintaining cell diversity at a population level.
Results
Discontinuous growth and cell-to-cell variability
As a first step to investigate interdependencies between growth and cell cycle entry at a single-cell level, we first determined volume growth kinetics of newborn cells by time-lapse microscopy (Fig. 1a, Supplementary Movie 1). Accurate volumetric data were obtained by means of a densitometric profile along radial axes to iteratively detect the optical boundary of the cell (Fig. 1b, Supplementary Movie 2) and obtain its volume assuming a prolate spheroid as shape. Volume distributions and average growth rates obtained by time-lapse microscopy were similar to those obtained with a Coulter Counter in batch cultures (Supplementary Fig. S1). Using this methodological approach, we observed a clearly discontinuous pattern of growth in volume during the cell cycle (Fig. 1c). First, cell volume increased in a fairly constant manner throughout G1, with only a slight decrease in the rate immediately before budding. Of particular relevance to the hypothesis of this work (see below), we observed a large cell-to-cell variability in volume growth rate during G1 similar to that observed in other biosynthetic processes38,39,40 (Supplementary Fig. S2). On the other hand, a strong increase in volume growth rate was observed approximately 20 min after budding, which decreased again to initial values by the end of the cell cycle (Fig. 1d, e). Thus, average volume growth rate in G2–M increased 2.3-fold compared with G1 phase, but individual rates in the two periods were not correlated significantly (Supplementary Fig. S2). This pattern of growth in volume was very similar in a different strain background (Supplementary Fig. S3). Overall, our data support previous studies that have proposed a rate-change point (RCP) for growth in volume shortly after budding41,42.
Cell size at Start is dictated by volume growth rate in G1
We next extended our analysis to determine cell volume at Start, when nuclear export of the transcriptional repressor Whi5 takes place23,43. Elutriated and asynchronously growing daughter cells were analysed using time-lapse microscopy (Supplementary Movies 3 and 4, respectively), and a Whi5—green fluorescent protein (GFP) nuclear index was quantified to determine cell volume at Start (Vs), time at Start (T1) and volume growth rate in G1 (α; Fig. 2a,b; Supplementary Fig. S4). Although we assumed linear kinetics to obtain volume growth rates in G1 as the simplest approximation to our data (Fig. 1c–e), very similar results were obtained if exponential kinetics were assumed (data not shown). Although there was no correlation between volume and time at Start, volume at Start was closely correlated with volume growth rate in G1 (Fig. 2c,e; Supplementary Tables S1 and S2), giving support to the hypothesis that individual growth rate could be a key determinant of cell size at Start4,29. If G1 length were a constant, this finding could have an obvious explanation, but this is not the case. First, time at Start fluctuated over a tenfold range, whereas volume at Start only varied across a threefold range. Second, time at Start did not correlate with growth rate (Fig. 2g; Supplementary Table S3). Finally, relative s.e. for volume and time at Start estimated from growth rate were 7.5% and 44.7%, respectively, which makes time a very unlikely determinant of volume at Start as a function of growth rate. As nutrients are potent modulators of growth rate in yeast, we analysed cells growing asynchronously in different carbon sources, and we found that volume at Start correlated with volume growth rate in G1 not only for each nutritional condition considered separately, but also when all results were pooled (Fig. 2d; Supplementary Table S1). As discussed above, it is very unlikely that time at Start mediates this correlation (Fig. 2h). Time from Start to budding (T2) showed less variability43 and it did not correlate with growth rate (Supplementary Fig. S5). Overall, our data suggest that individual growth rate determines cell volume at Start independently of the nutritional environment.
Growth rate during G1 governs cell size variability at Start
One approach to studying whether a cell cycle transition is controlled by a sizer is to test whether smaller cells gain more size before the transition44, and this can be effected using an αT-plot, where increase in volume (that is, αT, the product of volume growth rate and time) is plotted against initial volume37,44,45. If cell volume at Start (Vs) were a constant determined by external conditions, αT1 should correlate negatively with the initial volume (Vi) with a slope of −1 (Fig. 3a). Although we observed a negative correlation between αT1 and Vi in asynchronously growing cells, the slope was significantly different from −1 and the coefficient of determination was only 0.114, with data points being notably more scattered at small initial volumes (Fig. 3b; Supplementary Table S4). Very similar patterns have been interpreted to be the result of molecular noise inherent to the Start network37. However, we show here (Fig. 2c, d) that volume at Start (Vs) can be expressed as a linear function of the individual volume growth rate (α) by Vs=k.α+V0, where V0 is the minimum volume that a cell would attain at Start when its growth rate approaches zero, and k is a constant inherent to the sizer mechanism or, more formally, the time that a cell born with volume V0 would always spend in G1 independently of its growth rate. Once this function was incorporated into the equation for αT1, our data produced a completely different outcome (Fig. 3c). First, the slope was −0.91, a value not significantly different from −1 (Supplementary Table S5). Second, the coefficient of determination increased to 0.783, with an evident decrease in the scatter for small newborn cells. Despite the fact that Vi and T1 for elutriated cells could only be estimated, the results obtained were comparable to those obtained from asynchronously growing cells (Fig. 3b, c). More importantly, the same correlation was obtained with cells growing in different carbon sources (Fig. 3d, e), and a single regression line could be fitted independently of the carbon source used as nutrient. Thus, we conclude that variation in the individual growth rate and not inherent molecular noise is what causes the main variability in cell volume at Start.
Start network sets critical size as a function of growth rate
Cells execute Start through a network of proteins (Fig. 4a) that ultimately activate the G1/S regulon. Accordingly, we tested whether these proteins have a direct impact on cell size regulation by growth rate (Fig. 4, Supplementary Tables S1-S5). First, we analysed the role of two well-characterized inhibitors of Start: Whi3 and Whi5. Whi3 binds the CLN3 mRNA and contributes to retain the Cdc28–Cln3 complex in the cytoplasm in G1 cells31,32. On the other hand, Whi522,23 binds G1/S transcription factors14 to inhibit their activity in the nucleus during G1. Thus, we analysed cells deficient for Whi3 (Fig. 4b) or expressing whi5N (Fig. 4c), which produces an N-terminal fragment of Whi5 unable to inhibit the G1/S transition but still subject to cell-cycle-regulated export46. A priori, we expected that smaller cells at Start could be obtained by decreasing either k, the slope inherent to the sizer mechanism that defines the linear growth rate dependency, or V0, the minimal volume at Start for a cell with a growth rate approaching zero. We found that size at Start for both whi3 and whi5N cells was still dependent on the individual volume growth rate, but with a much smaller slope k compared with wild-type cells, indicating that both inhibitors participate in the sizer mechanism by linking cell size at Start to volume growth rate in G1. The small cell size phenotype is particularly important at high growth rates, which indicates that both negative regulators are needed in the mechanisms by which growth inhibits Start4,47.
We next analysed the role of Cln3, the most upstream G1 cyclin that activates Start16,17,18,19,20. We observed a clear dependence of cell size at Start on growth rate in Cln3-deficient cells (Fig. 4d), and similar results were obtained with cells lacking both Cln1 and Cln2 (Fig. 4f), the two G1 cyclins that contribute to Start in a positive feedback loop25,26,27. Thus, rather than being mediators, individual G1 cyclins would be target effectors of the sizing mechanism. However, Cln3- and Cln1,2-deficient cells showed a lower correlation between volume at Start and growth rate, indicating that their cooperative presence is important to achieve a robust growth-rate-dependent size control.
Cells deficient for Ydj1, the J-chaperone that releases the Cdc28–Cln3 complex from the ER to trigger Start33, were much larger than wild-type as expected but, more important, there was no correlation between cell volume at Start and volume growth rate in G1 (Fig. 4e). This result indicates that this chaperone is essential to determine cell size at Start as a function of growth rate.
Ribosome biogenesis activators Sfp1 and Sch9 also act as important negative regulators of Start34,47. As described, we found that Sfp1-deficient cells were smaller than wild-type34, but cell volume at Start still showed a very strong correlation with growth rate in G1 (Fig. 4g), suggesting that ribosome biogenesis regulates size as an upstream determinant of growth rate. In addition, as the correlation function was steeper in sfp1 than wild-type cells, Sfp1 could also have a positive role at Start, possibly by increasing translation rates and Cln3 levels4.
Cln3 and Ydj1 are both rate limiting for cell cycle entry19,33 and, as expected, their overexpression decreased the average cell volume at Start (Supplementary Table S6). Both the slope of the correlation function and the coefficient of determination were also smaller in oCLN3 (Fig. 4h) and oYDJ1 (Fig. 4i) cells, indicating that unusually high levels of Cln3 or Ydj1 produce a small cell size at Start because they weaken its dependence on growth rate in a similar manner to the lack of negative regulators Whi3 or Whi5.
Cells simultaneously lacking activating (Cln3) and inhibiting (Whi5, Stb1) components of the Start network display a budding size very similar to wild-type cells48. Accordingly, when we analysed cln3 whi5N stb1 cells by time-lapse microscopy, we found that average cell volume at Start was very similar to wild-type cells (Supplementary Table S6). However, correlation of individual cell volume at Start to growth rate in G1 was almost lost (Fig. 4j).
In summary, although deletion or overexpression of components of the Start network causes small or large cell phenotypes consistent with their positive or negative role, all these genetic alterations loose the strong dependency of the critical size on growth rate (Fig. 4k). Thus, a key role of the Start network in regulating cell size is by setting it as a function of growth rate for each individual cell.
Role of a growth-rate-dependent sizer
To test whether a growth-rate-dependent sizer is able to ensure size homeostasis at a population level, we used our experimental data to model growth in volume during the yeast cell cycle (Fig. 5a; a detailed description can be found in Supplementary Methods). We assumed that volume at Start is a stochastic variable that depends on growth rate, and we used two different approaches for modelling, a germline (Fig. 5b and Supplementary Data 1) and a pedigree algorithm (Fig. 6a and Supplementary Software 1). Validating their design, both modelling approaches predicted the observed independence of time to Start (T1) on growth rate (α), and produced αT-plots very similar to the experimental data (Fig. 5c,d and Supplementary Fig. S6). Predicted values for all dependent variables, that is, G1 length and volumes at the beginning of G1, Start, budding and division, were not significantly different from experimental data (Fig. 5e). Starting from a single cell, the pedigree approach generated a stable variability in daughter cell volume at Start in less than 10 generations (Fig. 6b) and predicted that mother cells should increase their volume in each generation (Fig. 6c) as experimentally observed13,49. Moreover, the pedigree analysis produced steady-state cell size distributions very similar to the experimental data obtained in different nutritional conditions (Fig. 6d). Germline modelling predicted a stable pattern of fluctuations in daughter cell volume at Start over 1,000 generations (Fig. 7c). As expected, much larger fluctuations were predicted by our model if volume growth rate in a cell cycle phase was made to depend on the initial volume at that phase (Supplementary Fig. S7). However, we found that growth rates in G1 (α) and S–G2–M (β) do not correlate with the respective initial volumes (Supplementary Fig. S7). Overall, both modelling approaches indicate that a single growth-rate-dependent size control within the cell cycle is sufficient to maintain cell size homeostasis.
Finally, we took advantage of the germline modelling approach to compare the growth-rate-dependent sizer with the steady sizer generally assumed up to now (Fig. 7a,b). Notably, G1 length fluctuations were much smaller for growth-rate-dependent (Fig. 7c) versus constant (Fig. 7d) sizer mechanisms. The predicted difference is so important that the average generation time for the population would be approximately 15% longer if critical size did not depend on growth rate. This effect is mainly due to slow-growing cells, because they have to grow for a longer time when the critical size is not reduced as a function of growth rate (Fig. 7e, f). Thus, G1 length variability is insensitive to growth rate variation only for the growth-rate-dependent sizer (Fig. 7g). As they showed no significant correlation between volume at Start and growth rate, ydj1 and whi5N stb1 cln3 cells would behave as if they were subject to a constant sizer mechanism. Indeed, slow-growing cells of both mutant strains displayed a clear increase in G1 length compared with wild-type cells with same growth rate (Fig. 7h).
Discussion
Our methodological approach to accurately determine the optical cell boundary has allowed us to observe a clearly discontinuous pattern of growth in volume during the cell cycle of budding yeast. Growth rate in volume was rather constant during G1 and it only showed a slight decrease immediately before budding. Notably, average volume growth rate in G2–M increased 2.3-fold compared with G1 phase. Mitchison41 had shown the existence of an RCP for growth in volume shortly after budding. On the other hand, although cell volumes were estimated from protein concentration, a recent report also has given support to this notion42. As protein synthesis seems to increase exponentially during the cell cycle37,42,50,51, a discontinuous pattern of growth in volume should cause changes in cell density over the cell cycle, consistent with the fact that yeast cells reach a maximum density shortly after budding52,53. All these observations point to a discontinuous mode of volume growth during the cell cycle of budding yeast, a property that has been well established in the fission yeast Schizosaccharomyces pombe44,45,54 and may also apply to mammalian cells55. Finally, from an opposite but complementary point of view, it is clear that cell growth is in turn governed by cell cycle progression56; so, in addition to gene-dosage effects as a result of DNA replication, shifting growth from apical to isotropic modes in G2 could give raise to the observed RCP.
It is generally accepted that budding yeast cells must reach a critical size before entering the cell cycle12,13. Nevertheless, individual cells offer a large variability in their size at Start under constant external and nutritional conditions36,37, and this fact has seriously challenged the existence of a deterministic mechanism setting the critical size in a dedicated manner. We have found that yeast cells set their size at Start as a robust function of the individual growth rate, the observed population variability in size at Start being a consequence of the intrinsic variation of growth rate in volume among cells. Our findings clearly support a deterministic mode of cell size control and rule out the notion of a universal critical size that, once defined by external conditions, had to be satisfied by all cells in the population to enter the cell cycle. Concerning the variability in volume growth rate, similar variations in different biosynthetic processes have already been reported in yeast38,39,40. However, we have found no (or very weak) correlations between volume growth rate and either ribosomal protein levels, GFP synthesis rates or ER contents (F.F., N.C. and M.A, unpublished results). Growth in volume is a very complex process, with many different factors likely suffering from some stochastic behaviour, which altogether could cause a large variability in the overall rate.
The dependency of cell size at Start on growth rate in G1 was extremely dependent on the proteins that control entry into the cell cycle and thus define the Start network. Small cells lacking negative regulators (Whi3, Whi5) or overexpressing positive effectors (Cln3, Ydj1) still showed a growth-rate-dependent volume at Start, but with lower coefficients of correlation compared with wild-type cells. Loss of G1 cyclins Cln3 or Cln1,2, or G1/S transcription factors, which produce a large-cell phenotype, also caused a decrease in the correlation between volume at Start and growth rate in G1. Moreover, cells simultaneously lacking activating (Cln3) and inhibiting (Whi5, Stb1) components of the Start network displayed a much weaker correlation of cell volume at Start to growth rate in G1, indicating that the Start network is deeply involved in the growth-rate-dependent sizer mechanism. Notably, loss of Ydj1 caused the most severe effects in the dependency of the critical size on growth rate. As Ydj1 is also involved in cell wall growth processes, such as protein translocation at the ER35, our data suggest that the key molecular device operating this correlation is at the interface between the cell cycle and cell growth machineries. In this regard, ribosome biogenesis activators Sfp1 and Sch9 also act as important negative regulators of Start47. However, we found that Sfp1-deficient cells still showed a very strong correlation between cell size at Start and growth rate in G1, which suggests that ribosome biogenesis would determine the critical size through its immediate effects on growth rate.
On the basis of a mathematical model of volume growth during the yeast cell cycle, we have shown that a growth-rate-dependent sizer is able to provide size homeostasis only if growth rate is made independent of the initial cell volume. Notably, we have found that growth rates in G1 and G2–M do not depend on the respective initial volumes, which confirms similar results obtained from α-factor-arrested cells in G1 (Supplementary Fig. S2b–d in ref. 39). An intrinsic control that sets the critical size as a function of growth rate also provides a broad solution for adaptation, thus making unnecessary the existence of a plethora of sensing mechanisms transmitting different external cues to the cell cycle machinery. Our experimental data suggest that, in addition to provide with homeostasis and adaptation, a growth-rate-dependent sizer would reduce noise in G1 length, which would in turn prevent loss of slow-growing cells and shorten the generation time at the population level (Fig. 7i). As slow-growing cells show a longer lifespan, these effects would provide the population with a significant evolutionary advantage.
Although robustness of the mechanism requires all Start network components, we have found that Ydj1 has the most essential role in dictating the critical size as a function of growth rate. We have proposed that high volume growth rates should require higher levels of chaperones at ribosome biogenesis, protein synthesis and protein translocation, thus decreasing available chaperone levels for other processes, such as Cln3 release from the ER29. As a consequence, a fast-growing cell would delay cell cycle entry until it has grown large enough and has accumulated sufficient chaperone levels to trigger Start. Indeed, our work demonstrates that budding yeast cells set their critical size at a single-cell level as a function of growth rate in volume and supports a model whereby competition for molecular chaperones could subordinate proliferation to growth and help cells to attain an individually tailored size.
Methods
Strain constructions
BY4741 WHI5::GFP57 (supplied by Invitrogen) and CML128 (a S288C derivative58) were used as wild-type strains of Saccharomyces cerevisiae. Standard PCR-based gene-replacement methods were used to obtain the indicated deletions, C-terminal translational fusions to GFP and mCherry, or transcriptional fusions to the GAL1 promoter at the respective chromosomal loci. To analyse Whi5-deficient cells, we used a GFP fusion of the N-terminal fragment of Whi5 (whi5N) unable to inhibit the G1/S transition but still subject to cell-cycle-regulated export46.
Growth conditions
Cells were grown under exponential conditions for at least 14 h in synthetic complete (SC) medium59 with either 2% glucose, 2% galactose, 2% raffinose or 3% ethanol at 30 °C before microscopy. CLN3 and YDJ1 overexpression was attained adding 2% galactose to cells that contained transcriptional fusions to the GAL1 promoter at the respective chromosomal loci. CLN3 overexpression (oCLN3) was evaluated in asynchronous daughter cells growing in SC plus galactose. YDJ1 overexpression (oYDJ1) was assessed in daughter cells elutriated 3 h after adding galactose to a culture grown in raffinose. As they contained a MET3p-CLN2 construct to avoid accumulation of aberrantly shaped cells, cln1,2 cells were grown in SC medium without methionine, which was added immediately before the time-lapse experiment27. Newly born cells were obtained by elutriation19. Coverslips (25 mm) were precoated at a small central area with 0.1 ml 0.1mg ml−1 concanavalin A type V (Sigma-Aldrich) in water for 60 min in a humid chamber, washed twice with water and kept at 4 °C in water for up to 24 h before use39. Cells (0.1 ml, optical density at 600 nm (OD600)=0.5) were added to the coverslip and left at room temperature for 2 min to obtain 20–40 bound cells per field. After two washes with 0.1 ml SC medium to remove unbound cells, the coverslip was assembled in a prewarmed RC-43C chamber (Warner Instruments, www.warneronline.com) with ca. 0.3 ml SC medium, such that concentration of bound cells within the chamber would correspond to an OD600 lower than 0.1. When elutriated cells were analysed, they were first suspended in SC medium at OD600=2 and bound to the centre of a concanavalin A-coated round 25-mm coverslip. Temperature at the centre of the chamber was maintained at 30±1 °C with the help of a DH-40i micro-incubator system and a Single Channel Heater Controller TC-324B (Warner Instruments).
Time-lapse microscopy
Cells were imaged by time-lapse microscopy under a fully motorized Axio Observer Z1 inverted microscope (Zeiss) equipped with a Plan-Apochromat ×100/1.40 objective, differential interference contrast (DIC) optics, an XY-motorized stage, a 100 W mercury lamp controlled by a FluoArc unit (Zeiss), adequate filters for GFP and mCherry detection, and an AxioCam MRm camera (Zeiss). Typically, 14 fields were automatically and repeatedly analysed every 3 min for 180 min using AxioVision 4.1 software (Zeiss) to obtain DIC, bright-field and fluorescence images. Autofocus was applied to the DIC image with a 50% threshold, and an offset of 1 μm was required to correct for an invariable deviation and obtain a precise focus of the cells before imaging. An additional offset of 1 μm was applied before taking the bright-field image to obtain the highest contrast at the cell boundary. To minimize irradiation of the cells, which were exposed for 1 s to obtain fluorescence images, illumination from the mercury lamp was reduced to 5% with the help of an appropriate neutral filter and the FluoArc unit. A diaphragm slider was used to limit illumination to the imaged field. A small percentage of cells (less than 2%) failed to grow and were not analysed further. Under these conditions, the average growth rate in cell volume was comparable to that obtained with a Coulter Counter from a batch culture (Supplementary Fig. S1). In addition, average growth rates from fields imaged or not for fluorescence were not significantly different (Supplementary Fig. S1).
Image analysis and cell measurements
We obtained cell dimensions from bright-field images (see ref. 60 for a discussion on the use of fluorescence or bright-field images for identifying cell boundaries). Specific procedures in Visual Basic using object libraries supplied by AxioVision (Zeiss) software were written to detect cell boundaries as pixels markedly darker compared with both the surrounding background and the cell interior. The algorithm first obtains a densitometric profile along a radial axis from a seed point within the cell (initially provided by the user) to find two consecutive maximum and minimum pixel values produced around the cell boundary, which is arbitrarily established at an intermediate pixel having their arithmetic grayscale mean as value (Fig. 1b; Supplementary Movie 2). Then, cell boundary pixels are iteratively defined by full rotation of the radial axis and a preliminary ellipse object is fitted. Outliers are first eliminated by their distance to the centre of the ellipse when it is larger than a fixed value. Second, additional outliers are eliminated when their distance to neighbour boundary pixels is larger than a fixed value. Once outliers have been removed, missing boundary pixels are estimated by fitting ellipsoidal segments. This process was essential to analyse daughter cells before cell separation. A final ellipse is then fitted to the obtained boundary pixel array, and major and minor axes are used to calculate the cell volume assuming a prolate as shape. The cell was followed through consecutive time-lapse images by using the centre of the ellipse as seed point to obtain radial profiles in the following time-point image. Volume distributions generated from time-lapse microscopy data were similar to those obtained with the Coulter Counter in batch cultures (Supplementary Fig. S1). These procedures have been rewritten in Java as BudJ, an ImageJ plugin that can be obtained from www.ibmb.csic.es\home\maldea or, see Supplementary Software 2.
Quantification of Whi5 nuclear localization
To avoid colocalization with a compatible fluorescent nuclear marker, which would necessarily increase light damage during growth, Whi5 localization in the nucleus was assessed as a relative index inversely proportional to the mean distance of the brightest pixels to their gravity centre (Supplementary Fig. S4). Once cell boundaries are determined as described above, brightest pixels within the cell are first selected having a grey value larger than the mean plus two s.d. Next, the mean distance to their gravity centre is calculated and made relative to 8r/15, the theoretical mean distance to the centre from random points in a sphere projected to a circle with radius r. The nuclear index thus obtained when Whi5 has been already exported to the cytoplasm was very close to one (Fig. 2b). In addition, this approach produced results that correlated well with those obtained by colocalization with Htb2-mCherry as a nuclear marker (Supplementary Fig. S4). Briefly, the Htb2-mCherry image was thresholded to obtain a binary mask and delimit nuclear and cytoplasmic regions within each cell identified in the bright-field image as described above. Next, colocalization was calculated as the mean brightness of pixels within the nuclear region made relative to that of pixels within the cytoplasmic compartment.
Determination of cell parameters
Growth rates in volume during G1 phase (α) were obtained by linear regression of volume data during this time period. Occasionally, a change in rate was clearly observed within G1, and in this case, individual growth rate in volume was obtained by linear regression of volume data during the 60 min before Whi5 nuclear export. Instantaneous growth rates as a function of time during G1 were estimated from data ±6 min around every time point. In asynchronous cultures, T1 was the time period from Whi5 import to Whi5 export (Start), and T2 was the time period from Start to budding. Once nuclear index values were smoothed by a sliding mean over three consecutive points, we defined Whi5 nuclear entry as the first time point with a nuclear index significantly different from 1, and Whi5 nuclear export as the first time point with a nuclear index not significantly different from 1. In elutriation experiments, we collected less than 0.2% of the cells loaded in the rotor, being the smallest cells among all newly born cells in the asynchronous culture. We could only approximate T1 in elutriated cells as the period from time 0 (when we placed the newborn cells in the chamber with prewarmed medium) to Whi5 nuclear export. Thus, true Vi and T1 values for elutriated cells are expected to be slightly smaller and longer, respectively, than those measured. Initial cell volumes (Vi) and volumes at Start (Vs) and budding (Vb) were obtained directly from raw data of the corresponding time points. T3 was the time at division and was approximated as the time from budding to Whi5 import to the nucleus. Similarly, volume at division (Vd) was also obtained at Whi5 import. Although growth rates for the G2–M period were estimated by linear regression of volume data from 30 to 60 min after budding, growth rates for the S-G2–M period (β) were obtained as (Vd–Vb)/T3. Finally, a few cells (less than 5% of all cells analysed) did not trigger Start during the course of the experiment or displayed growth rate differences to the average larger than two s.d. and were discarded as outliers.
Additional information
How to cite this article: Ferrezuelo, F. et al. The critical size is set at a single-cell level by growth rate to attain homeostasis and adaptation. Nat. Commun. 3:1012 doi: 10.1038/ncomms2015 (2012).
References
Donachie, W. D. & Blakely, G. W. Coupling the initiation of chromosome replication to cell size in Escherichia coli. Curr. Opin. Microbiol. 6, 146–150 (2003).
Mitchison, J.M. Growth during the cell cycle. Int. Rev. Cytol. 226, 165–258 (2003).
Conlon, I. & Raff, M. Control and maintenance of mammalian cell size: Response. BMC Cell Biol. 5, 36 (2004).
Jorgensen, P. & Tyers, M. How cells coordinate growth and division. Curr. Biol. 14, R1014–1027 (2004).
Cavalier-Smith, T. Economy, speed and size matter: Evolutionary forces driving nuclear genome miniaturization and expansion. Ann. Bot. 95, 147–175 (2005).
Grebien, F., Dolznig, H., Beug, H. & Mullner, E.W. Cell size control: new evidence for a general mechanism. Cell Cycle 4, 418–421 (2005).
Baserga, R. Is cell size important? Cell Cycle 6, 814–816 (2007).
Zaritsky, A., Vischer, N. & Rabinovitch, A. Changes of initiation mass and cell dimensions by the 'eclipse'. Mol. Microbiol. 63, 15–22 (2007).
Cook, M. & Tyers, M. Size control goes global. Curr. Opin. Biotechnol. 18, 341–350 (2007).
Leslie, M. How does a cell know its size? Science 334, 1047–1048 (2011).
Mortimer, R. K. Radiobiological and genetic studies on a polyploid series haploid to hexaploid of Saccharomyces cerevisiae. Rad. Res. 9, 312–326 (1958).
Hartwell, L. H. & Unger, M. W. Unequal division in Saccharomyces cerevisiae and its implications for the control of cell division. J. Cell Biol. 1, 422–435 (1977).
Johnston, G. C., Pringle, J. R. & Hartwell, L. H. Coordination of growth with cell division in the yeast Saccharomyces cerevisiae. Exp. Cell Res. 105, 79–98 (1977).
Nasmyth, K. At the heart of the budding yeast cell cycle. Trends Genet. 12, 405–412 (1996).
Futcher, B. Transcriptional regulatory networks and the yeast cell cycle. Curr. Opin. Cell Biol. 14, 676–683 (2002).
Sudbery, P. E., Goodey, A. R. & Carter, B. L. Genes which control cell proliferation in the yeast Saccharomyces cerevisiae. Nature 288, 401–404 (1980).
Cross, F. R. DAF1, a mutant gene affecting size control, pheromone arrest, and cell cycle kinetics of Saccharomyces cerevisiae. Mol. Cell Biol. 8, 4675–4684 (1988).
Nash, R., Tokiwa, G., Anand, S., Erickson, K. & Futcher, B. The WHI1+ gene of Saccharomyces cerevisiae tethers cell division to cell size and is a cyclin homolog. EMBO J. 7, 4335–4346 (1988).
Tyers, M., Tokiwa, G. & Futcher, B. Comparison of the Saccharomyces cerevisiae G1 cyclins: Cln3 may be an upstream activator of Cln1, Cln2 and other cyclins. EMBO J. 12, 1955–1968 (1993).
Dirick, L., Bohm, T. & Nasmyth, K. Roles and regulation of Cln/Cdc28 kinases at the start of the cell cycle of Saccharomyces cerevisiae. EMBO J. 14, 4803–4813 (1995).
Mendenhall, M. D. & Hodge, A. E. Regulation of Cdc28 cyclin-dependent protein kinase activity during the cell cycle of the yeast Saccharomyces cerevisiae. Microbiol. Mol. Biol. Rev. 62, 1191–1243 (1998).
de Bruin, R. A. M., McDonald, W. H., Kalashnikova, T. I., Yates, J. III & Wittenberg, C. Cln3 activates G1-specific transcription via phosphorylation of the SBF-bound repressor Whi5. Cell 117, 887–898 (2004).
Costanzo, M. et al. CDK activity antagonizes Whi5, an inhibitor of G1/S transcription in yeast. Cell 117, 899–913 (2004).
Kaelin Jr, W. G. Functions of the retinoblastoma protein. Bioessays 21, 950–958 (1999).
Cross, F. R. & Tinkelenberg, A. H. A potential positive feedback loop controlling CLN1 and CLN2 gene expression at the start of the yeast cell cycle. Cell 65, 875–883 (1991).
Dirick, L. & Nasmyth, K. Positive feedback in the activation of G1 cyclins in yeast. Nature 351, 754–757 (1991).
Skotheim, J. M., Di Talia, S., Siggia, E. D. & Cross, F. R. Positive feedback of G1 cyclins ensures coherent cell cycle entry. Nature 454, 291–296 (2008).
McInerny, C. J., Partridge, J. F., Mikesell, G. E., Creemer, D. P. & Breeden, L. L. A novel Mcm1–dependent element in the SWI4, CLN3, CDC6, and CDC47 promoters activates M/G1-specific transcription. Genes Dev. 11, 1277–1288 (1997).
Aldea, M., Garí, E. & Colomina, N. Control of cell cycle and cell growth by molecular chaperones. Cell Cycle 6, 2599–2603 (2007).
Nash, R. S., Volpe, T. & Futcher, B. Isolation and characterization of WHI3, a size–control gene of Saccharomyces cerevisiae. Genetics 157, 1469–1480 (2001).
Garí, E. et al. Whi3 binds the mRNA of the G1 cyclin CLN3 to modulate cell fate in budding yeast. Genes Dev. 15, 2803–2808 (2001).
Wang, H., Garí, E., Vergés, E., Gallego, C. & Aldea, M. Recruitment of Cdc28 by Whi3 restricts nuclear accumulation of the G1 cyclin-CDK complex to late G1. EMBO J. 23, 180–190 (2004).
Vergés, E., Colomina, N., Garí, E., Gallego, C. & Aldea, M. Cyclin Cln3 is retained at the ER and released by the J chaperone Ydj1 in late G1 to trigger cell cycle entry. Mol. Cell 26, 649–662 (2007).
Jorgensen, P., Nishikawa, J. L., Breitkreutz, B. J. & Tyers, M. Systematic identification of pathways that couple cell growth and division in yeast. Science 297, 395–400 (2002).
Caplan, A. J., Cyr, D. M. & Douglas, M. G. Ydj1p facilitates polypeptide translocation across different intracellular membranes by a conserved mechanism. Cell 71, 1143–1155 (1992).
Lord, P. G. & Wheals, A. E. Variability in individual cell cycles of Saccharomyces cerevisiae. J. Cell Sci. 50, 361–376 (1981).
Di Talia, S., Skotheim, J. M., Bean, J. M., Siggia, E. D. & Cross, F. R. The effects of molecular noise and size control on variability in the budding yeast cell cycle. Nature 448, 947–951 (2007).
Blake, W. J., Kaern, M., Cantor, C. R. & Collins, J. J. Noise in eukaryotic gene expression. Nature 422, 633–637 (2003).
Colman-Lerner, A. et al. Regulated cell-to-cell variation in a cell-fate decision system. Nature 437, 699–706 (2005).
Volfson, D. et al. Origins of extrinsic variability in eukaryotic gene expression. Nature 439, 861–864 (2006).
Mitchison, J. M. The growth of single cells. II. Saccharomyces cerevisiae. Exp. Cell Res. 15, 214–221 (1958).
Cookson, N. A., Cookson, S. W., Tsimring, L. S. & Hasty, J. Cell cycle-dependent variations in protein concentration. Nucleic Acids Res. 38, 2676–2681 (2010).
Bean, J. M., Siggia, E. D. & Cross, F. R. Coherence and timing of cell cycle start examined at single-cell resolution. Mol. Cell 21, 3–14 (2006).
Fantes, P. A. Control of cell size and cycle time in Schizosaccharomyces pombe. J. Cell Sci. 24, 51–67 (1977).
Sveiczer, A., Novak, B. & Mitchison, J. M. The size control of fission yeast revisited. J. Cell Sci. 109, 2947–2957 (1996).
Taberner, F. J., Quilis, I. & Igual, J. C. Spatial regulation of the start repressor Whi5. Cell Cycle 8, 3010–3018 (2009).
Jorgensen, P. et al. A dynamic transcriptional network communicates growth potential to ribosome synthesis and critical cell size. Genes Dev. 18, 2491–2505 (2004).
Wang, H., Carey, L. B., Cai, Y., Wijnen, H. & Futcher, B. Recruitment of Cln3 cyclin to promoters controls cell cycle entry via histone deacetylase and other targets. PLoS Biol. 7, e1000189 (2009).
Yang, J. et al. Cell size and growth rate are major determinants of replicative lifespan. Cell Cycle 10, 144–155 (2011).
Elliott, S. G. & McLaughlin, C. S. Rate of macromolecular synthesis through the cell cycle of the yeast Saccharomyces cerevisiae. Proc. Natl Acad. Sci. USA 75, 4384–4388 (1978).
Godin, M. et al. Using buoyant mass to measure the growth of single cells. Nat. Methods 7, 387–390 (2010).
Hartwell, L. H. Periodic density fluctuation during the yeast cell cycle and the selection of synchronous cultures. J. Bacteriol. 104, 1280–1285 (1970).
Bryan, A. K., Goranov, A., Amon, A. & Manalis, S. R. Measurement of mass, density, and volume during the cell cycle of yeast. Proc. Natl Acad. Sci. USA 107, 999–1004 (2010).
Mitchison, J. M. & Nurse, P. Growth in cell length in the fission yeast Schizosaccharomyces pombe. J. Cell Sci. 75, 357–376 (1985).
Mir, M. et al. Optical measurement of cycle-dependent cell growth. Proc. Natl Acad. Sci. USA 108, 13124–13129 (2011).
Goranov, A. I. et al. The rate of cell growth is governed by cell cycle stage. Genes Dev. 23, 1408–1422 (2009).
Huh, W. et al. Global analysis of protein localization in budding yeast. Nature 425, 686–691 (2003).
Gallego, C., Garí, E., Colomina, N., Herrero, E. & Aldea, M. The Cln3 cyclin is downregulated by translational repression and degradation during the G1 arrest caused by nitrogen deprivation in budding yeast. EMBO J. 16, 7196–7206 (1997).
Sambrook, J. & Russell, D. W. Molecular Cloning: A Laboratory Manual 3rd edn (CSHL: Cold Spring Harbor, NY), (2001).
Gordon, A. et al. Single-cell quantification of molecules and rates using open-source microscope-based cytometry. Nat. Methods 4, 175–181 (2007).
Acknowledgements
We thank A. Cornadó and S. Rius for technical assistance, J.C. Igual for communicating unpublished results and R. Y. Tsien for the mCherry fusion vector. We also thank A. Aubanell, T. Mazza and J. Sánchez-Matamala for key preliminary contributions, and B. Novak, P. Nurse and J. Torres-Rosell for stimulating discussions. This work was funded by the Ministry of Science and Innovation of Spain (BFU2010-20205), Consolider-Ingenio 2010 (CSD2007-15), the Italian Research Fund FIRB (RBPR0523C3) and the European Union (FEDER). F.F. and N.C. are researchers of the Ramón y Cajal programme.
Author information
Authors and Affiliations
Contributions
F.F., N.C. and M.A. performed the experiments. E.G., C.G. and M.A. conceived and designed the project and developed the required algorithms for image analysis, A.P., A. C.-N. and M.A. developed and analysed the model. M.A. wrote the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Figures, Tables and Methods
Supplementary Figures S1-S7, Supplementary Tables S1-S7 and Supplementary Methods (PDF 471 kb)
Supplementary Data 1
Germline simulation spreadsheet (XLSX 1742 kb)
Supplementary Movie 1
Time-lapse microscopy of a new-born daughter cell during G1 (MOV 130 kb)
Supplementary Movie 2
Determination of cell boundaries for cell size computation (MOV 1043 kb)
Supplementary Movie 3
Cell size analysis and Whi5-GFP localization of an elutriated cell during cell cycle entry (MOV 225 kb)
Supplementary Movie 4
Cell size analysis and Whi5-GFP localization during cell cycle entry of a newlyborn daughter cell in an asynchronous culture (MOV 243 kb)
Supplementary Software 1
Pedigree simulation software (DOC 60 kb)
Supplementary Software 2
BudJ analysis software (DOC 201 kb)
Rights and permissions
About this article
Cite this article
Ferrezuelo, F., Colomina, N., Palmisano, A. et al. The critical size is set at a single-cell level by growth rate to attain homeostasis and adaptation. Nat Commun 3, 1012 (2012). https://doi.org/10.1038/ncomms2015
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms2015
This article is cited by
-
Temporal segregation of biosynthetic processes is responsible for metabolic oscillations during the budding yeast cell cycle
Nature Metabolism (2023)
-
A high-throughput microfluidic diploid yeast long-term culturing (DYLC) chip capable of bud reorientation and concerted daughter dissection for replicative lifespan determination
Journal of Nanobiotechnology (2022)
-
Determination of protoplast growth properties using quantitative single-cell tracking analysis
Plant Methods (2022)
-
High-resolution mass measurements of single budding yeast reveal linear growth segments
Nature Communications (2022)
-
Sublinear scaling of the cellular proteome with ploidy
Nature Communications (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.