Long-term microfluidic tracking of coccoid cyanobacterial cells reveals robust control of division timing
© The Author(s). 2017
Received: 31 July 2016
Accepted: 10 December 2016
Published: 14 February 2017
Cyanobacteria are important agents in global carbon and nitrogen cycling and hold great promise for biotechnological applications. Model organisms such as Synechocystis sp. and Synechococcus sp. have advanced our understanding of photosynthetic capacity and circadian behavior, mostly using population-level measurements in which the behavior of individuals cannot be monitored. Synechocystis sp. cells are small and divide slowly, requiring long-term experiments to track single cells. Thus, the cumulative effects of drift over long periods can cause difficulties in monitoring and quantifying cell growth and division dynamics.
To overcome this challenge, we enhanced a microfluidic cell-culture device and developed an image analysis pipeline for robust lineage reconstruction. This allowed simultaneous tracking of many cells over multiple generations, and revealed that cells expand exponentially throughout their cell cycle. Generation times were highly correlated for sister cells, but not between mother and daughter cells. Relationships between birth size, division size, and generation time indicated that cell-size control was inconsistent with the “sizer” rule, where division timing is based on cell size, or the “timer” rule, where division occurs after a fixed time interval. Instead, single cell growth statistics were most consistent with the “adder” rule, in which division occurs after a constant increment in cell volume. Cells exposed to light-dark cycles exhibited growth and division only during the light period; dark phases pause but do not disrupt cell-cycle control.
Our analyses revealed that the “adder” model can explain both the growth-related statistics of single Synechocystis cells and the correlation between sister cell generation times. We also observed rapid phenotypic response to light-dark transitions at the single cell level, highlighting the critical role of light in cyanobacterial cell-cycle control. Our findings suggest that by monitoring the growth kinetics of individual cells we can build testable models of circadian control of the cell cycle in cyanobacteria.
Cyanobacteria are ancient oxygenic photoautotrophs with important roles in global carbon and nitrogen cycles, and hold promise as chassis organisms for products such as biofuels . Cyanobacteria possess a circadian clock and cell-cycle regulation that allow them to robustly respond to diel cycles. Synchronized populations of the unicellular cyanobacterium Synechococcus elongatus PCC7942 have been used to identify the main components responsible for circadian oscillations . Another model species, Synechocystis sp. PCC6803 (hereafter Synechocystis), has played an important role in elucidating photosynthetic pathways  and phototaxis [4–6], in addition to providing insight into circadian cycle regulation [7, 8]. Synechocystis can be engineered to produce many biomolecules . However, it remains unknown how the cell cycle is coupled with growth (here referring to volume expansion) in single cells and across generations and how this coupling is influenced by diel cycles. A detailed understanding of the phenotypic heterogeneity across populations and how environmental factors such as rapid changes in light affect growth may provide insight into how cells integrate external stimuli with internal mechanisms of cell-cycle and cell-size regulation. This understanding will also be required for optimizing the efficiency of large-scale Synechocystis bioreactors.
Bacteria typically maintain a size and shape that is characteristic of the species, suggesting that cell-size control is fundamental across the kingdom. Most studies of bacterial growth have focused on fast-growing heterotrophs such as Escherichia coli , Caulobacter crescentus , Bacillus subtilis , and Pseudomonas aeruginosa , which differ in many respects from slow-growing cells such as Synechocystis. Recently, microscopy has been used to track single fast-growing cells on agar pads or in microfluidic devices and to characterize correlations between cell size and generation time (defined as the time between cell birth and cell division). For several organisms, studies have demonstrated that size homeostasis is maintained via an adder rule whereby cells increase by a constant volume each generation regardless of birth size . These studies have focused almost entirely on rod-shaped bacteria with short generation times of less than 1 h; it remains to be seen whether similar homeostatic behaviors are exhibited by cells with other morphologies and/or much longer doubling times.
Several technical challenges complicate the single-cell microscopy-based analysis of slow-growing cocci such as Synechocystis. Although their small size (1–2 μm) is typical of many model bacteria, Synechocystis and other cyanobacteria require light and carbon dioxide for photosynthesis. Evaporation makes hydrogel surfaces unfit for long-term tracking of slow-growing cells. Microfluidics alleviates problems associated with evaporation, but devices can be difficult to use, particularly in high throughput, due to lack of automation and system-level integration of a comprehensively controlled microfluidic system including microscope, stage, image acquisition, and actuation of microfluidic valves. In addition, some microfluidic devices have been designed to exploit the elongation of rod-shaped cells along only one direction [14, 15]; such one-dimensional expansion is unlikely to be the case for many non-rod-shaped organisms and hence mechanical constraint within a micron-sized channel would not reflect normal growth. To address these issues, we modified a microfluidic cell-culture system for monitoring Synechocystis growth and division over several generations in continuous illumination or with light-dark cycling . We determined that cells undergo exponential growth during times of illumination, with expansion and division almost completely inhibited in the dark. Sister-cell pairs exhibited highly correlated generation times, even maintaining synchrony throughout dark periods. By comparing our experimental data to simulations of various cell-size control models, we found that Synechocystis cells are unlikely to follow the ‘sizer’ or ‘timer’ models; instead, the ‘adder’ rule of constant volume increment better explains the observed trends. In summary, our analyses reveal how light plays a critical role and is tightly integrated with the Synechocystis cell cycle.
Microfluidics and probabilistic image analysis facilitate long-term quantification of growth behavior
The coccoid shape and small size of Synechocystis cells make robust identification of cell division events challenging. To address this, we developed an automated image analysis pipeline to track cell positions and to identify newly divided sister cells in a set of time-lapse frames (Fig. 1b, Additional file 2: Figure S2). The key advantage of our analysis method is a probabilistic framework specifically trained on Synechocystis morphologies (Additional file 3: Figure S3, Additional file 4). This framework avoids hard thresholds that define cell boundaries and division events, and allows for correction of classification errors using information from the changes in cell shape over time. Moreover, in cases where a pair of cells is not accurately segmented, the algorithm still classifies the cluster as distinct from a single cell, avoiding lineages with artefactually high division times due to missing the division event. Our image analysis method can operate solely on bright-field or phase-contrast microscopy images, eliminating the dependence on fluorescence images for cell segmentation. This aspect is particularly important for cyanobacteria, which exhibit high levels of auto-fluorescence. In general, removing the requirement of fluorescence imaging also avoids potential inhibition of cell growth due to fluorescence excitation , or frees up the fluorescence channel for other applications.
To determine the growth dynamics of Synechocystis cells over multiple generations, we estimated the volume of individual cells by assuming rotational symmetry of the cell contour (Additional file 5: Figure S4) and tracked cell lineages from the single-cell stage for 60 h in 20 different chambers (Fig. 1c, Methods). We observed that all cells grew, though at different rates (Additional file 6: Figure S5A). Total volume of all lineages, normalized to the volume of the initial cell in the first frame, increased approximately exponentially for the first 40 h (Fig. 1c). Mean residuals after fitting two separate sections of the lineage growth curve further confirmed exponential growth (Fig. 1d). At later times, lineage growth rate slowed down, presumably reflecting the consumption of nutrients in the medium as cell density increased over the course of the experiment. To determine whether rates of division were coordinated with lineage growth rates, we automatically counted the number of cells in each lineage over time and found that mean cell number increased at the same rate as the mean lineage volume (Fig. 1c), suggesting cell-size homeostasis. The deviation between mean cell number and mean volume in the final 20 h is due, at least in part, to the presence of clusters with many cells in which accurate number quantification is challenging. Regardless, the combination of our experimental and analysis platforms enables rapid and robust quantification of bacterial growth and division across multiple days, empowering long-term single-cell analyses of slow-growing species and ellipsoidal cells such as Synechocystis.
Synechocystis cell volume expands exponentially under continuous light
Additional file 7: Movie S1. Time-lapse video showing three sets of Synechocystis sister cells, grown under continuous illumination, moving apart by a few microns after division. This movement usually takes place in 10-20 min, which is short compared to the mean generation time (16.9 h), and is more prevalent for isolated doublets than for larger clusters of cells.
The exponential growth of a microcolony (Fig. 1c, d) does not automatically imply exponential growth of individual cells over the cell cycle. To examine whether single cells also expanded their biomass exponentially or underwent distinct growth phases during their cell cycle, we quantified the volume of single cells for which boundaries could be confidently identified throughout their cell cycle (n = 140, Additional file 8: Movie S2). Most cells continuously increased in volume exponentially throughout the cell cycle under continuous illumination (Fig. 2e), even though lineage growth eventually slowed during the experiment, suggesting that they were growing in a relatively constant environment throughout their cell cycle. Therefore, our microfluidic device supports exponential expansion of cells and cell populations over multiple days and multiple cell-division events.
Additional file 8: Movie S2. Time-lapse video tracking two cell cycles of growth and division under continuous illumination starting from a single Synechocystis cell. Red boundaries outline single cells extracted by our image analysis pipeline.
Growth and division of Synechocystis cells are rapidly inhibited in the dark
Unlike most heterotrophic fast-growing bacterial species whose growth has been characterized at the single-cell level, cyanobacteria divide relatively slowly, rely on photosynthesis for energy, possess a robust circadian cycle, and respond to environmental light stimuli . Thus, it is important to determine the growth dynamics of Synechocystis cells under light-dark cycles that are similar to conditions encountered in the environment. Most previous studies have entrained cyanobacteria using light-dark cycles and then observed free-running behavior under continuous illumination ; however, this strategy does not reveal how quickly cells respond to changes in light conditions or if there is heterogeneity in cellular responses. Our microfluidic culture system has the advantage of allowing direct observation of Synechocystis cells during the dark phase, using short (millisecond) pulses of low-intensity light to record bright-field images (Additional file 4).
Additional file 10: Movie S3. Time-lapse video comparing Synechocystis movement during an illumination period (left) with a dark period (right). Left and right panels represent the same area of the same microfluidic chamber. Growth and division are apparent under illumination, while increased movement occurs in the dark period.
In the dark, cell division also halted, even in those cells with substantial constriction prior to the LED being switched off (Fig. 3c, Additional file 11: Movie S4). In the subsequent illumination period, cells completed cytokinesis. Only 6/547 (1%) of division events were observed in the dark, all of which occurred within 30 min after the light was turned off (Fig. 3d). The timing of division events displayed no preference for the beginning or end of the illuminated intervals. There was a peak of division events in the middle of the first interval, while the distribution was approximately uniform in the second and third intervals (Fig. 3d). We observed an increase in the number of divisions in the first light interval that stabilized by the second light interval. The initial increase was largely due to a burst of divisions that occurred once cells began incubation in the device. We do not know the origin of this synchronization, but we note that the first division event for each cell does not contribute to our generation time statistics because we can only measure birth time after the first division has taken place. Taken together, our results indicate that light is necessary for both growth and division of Synechocystis cells.
Additional file 11: Movie S4. Time-lapse video of a single Synechocystis cell grown under light-dark cycles. Cell growth and division do not occur during dark periods, but resume rapidly after light is restored.
Sister cells have similar generation times whether grown under continuous light or light-dark cycles
To extract single-cell growth related parameters from experiments under light-dark cycles, we ignored intervals of single-cell growth curves in the dark, in which neither growth nor division was observed (Methods). While the distributions of growth rates and generation times were similar under light-dark cycles compared with continuous illumination, the distribution of cell sizes was slightly smaller under light-dark cycles (Additional file 14: Figure S9). Cell birth and division volume distributions had coefficients of variation of 0.12 and 0.13 in continuous light and 0.15 and 0.17 in light-dark cycles, respectively, in close agreement with the coefficients of variation reported for other bacterial species . Sister cell birth volumes were also highly correlated, indicating that cells generally divided symmetrically, in both continuous light (R = 0.95, P < 1 × 10–69, Fig. 3e) and light-dark cycles (R = 0.83, P < 1 × 10–11, Fig. 4f). Nonetheless, there were a few cells that divided asymmetrically (8/139 cell pairs with birth volume asymmetry > 7%) (Fig. 4e). Interestingly, the resulting daughter cell pairs exhibited large differences in division timing (Fig. 4c), indicating that division asymmetry may influence the ability of daughter cells to maintain their otherwise synchronized generation timing. In summary, the striking similarities between sister cell generation times under continuous light and light-dark growth conditions suggest that the underlying regulatory mechanism is suspended in the dark but otherwise unaffected by light input.
Synechocystis cell-cycle statistics are not consistent with regulation of division timing based on fixed division size or cell-cycle interval
Under ideal conditions (constant mean cell size and normally distributed growth rates that are independent of cell size), the sizer model predicts that division volume should be independent of birth volume, while the adder and timer models predict slopes of +1 and +2, respectively. However, the expected values of these slopes are altered somewhat due to experimental noise and deviations from ideal conditions. To incorporate how distributions of our measured quantities modify the predicted slopes, we extended a governing set of equations to take into consideration imperfect distributions of various single-cell growth parameters (Additional file 4) [25, 29]. Then, we simulated exponentially growing cells using the three models with noise distributions extracted from our experimental data (Fig. 4a and b, Additional file 14: Figure S9). To make the simulations more comparable to our experiments, we also used our experimentally measured distributions of growth rates and birth sizes. Our measurements of cells under continuous illumination revealed a significant correlation between division and birth volume with a slope of 0.75 (Fig. 5b, P < 1 × 10–7). Compared to the sizer and timer models, this slope most closely mimicked simulations of the adder model (Fig. 5b, slope m = 0.97 ± 0.12) and was inconsistent with simulations of either the sizer or timer models (Fig. 5b). Also consistent with only the adder model, increment volume was uncorrelated with birth volume (Fig. 5c, m = –0.25, P = 0.07). The timer model predicts in ideal conditions that generation time is independent of birth volume, whereas the adder and sizer models predict similar inverse relationships. The experimentally determined negative slope of –0.79 for generation time with respect to birth volume indicated that smaller cells take longer to divide than larger cells, and was in reasonable agreement with our simulations of the adder model (Fig. 5d, m = –0.48 ± 0.05, P = 7 × 10–4). Finally, normalized differences in generation times between sister cells were negatively correlated to the asymmetry in birth volumes (Additional file 15: Figure S10A, slope = –1.04, P = 4 × 10–8), indicating that the smaller of the sister cells tended to spend a longer time growing before dividing. The slope was closest to that of simulations based on the adder model (Additional file 15: Figure S10A, slope = –0.72 ± 0.29, P = 0.07). Thus, while it remains possible that Synechocystis cell-size regulation follows a rule that differs subtly from the adder model, Synechocystis growth under continuous illumination is clearly inconsistent with the sizer or timer models.
To determine whether light-dark cycles altered the regulation of cell-cycle timing, we computed generation times ignoring the dark periods (as in Fig. 4c, d, Methods) and performed simulations of each control model, sampling birth volumes and growth rates from our light-dark cycle experiment. As with continuous illumination, slopes of division (Fig. 5e), increment volume (Fig. 5f), and generation time (Fig. 5g) as a function of birth volume were more consistent with the adder model compared to the sizer and timer models. The data for generation time asymmetry and birth volume asymmetry were too noisy to determine the significance of the relationships (Additional file 15: Figure S10B). Thus, Synechocystis cell growth and division behaviors under light-dark cycles provide further support against the sizer and timer models, independent of intervening dark intervals.
Cyanobacteria are significantly impacted by light and nutrient status. Hence, studying and modeling their growth kinetics provide a useful paradigm for how complex environmental inputs are integrated into cell-cycle control in photosynthetic microorganisms. To determine growth behaviors, size-control mechanisms, and the role of light in cell-cycle progression, we tracked single-cell growth kinetics of Synechocystis in a modified microfluidic cell culture system under continuous illumination and light-dark cycles (Fig. 1). With features such as integrated LED lighting and automated refocusing and image acquisition, our microfluidic cell-culture device allows facile multiplexing and long-term tracking of single cells for days, enabling the study of slow-growing organisms such as Synechocystis. Moreover, our device does not constrain the movement or growth directions of cells. This aspect is critical for Synechocystis cells, whose division planes rotate by 90° every generation (Fig. 2c), and is in contrast to “mother machine” devices  that exploit the one-dimensional elongation of rod-shaped organisms to track cells.
Most previous studies of circadian control in cyanobacteria have used the rod-shaped Synechococcus elongatus sp. PCC7942, for which batch cultures were entrained over several light-dark cycles, followed by fluorescence imaging of circadian-clock proteins under continuous illumination [30, 31]. In such experiments, expression levels of circadian genes have been observed to oscillate during intervals classified subjectively as “light” and “dark” [2, 32], suggesting that a direct light input can entrain the system and that expression of circadian genes may gate cell division . However, recent studies have shown that clock genes also respond to the ADP/ATP ratio within the cell, which is a read out of metabolic status determined by rates of photosynthesis during the light period . Thus, cyanobacterial growth and division can also be affected by light through metabolism, and cell behaviors after entrainment but under continuous illumination are likely distinct from phenotypes that emerge after transfer to a dark environment in which energetics also change dramatically. Our microfluidic platform provides the ability to directly observe the growth behavior of single Synechocystis cells during the dark phase, with short, low-intensity light exposures. The level of light used is sufficient for accurate cell tracking and demonstrably does not induce any cell growth in the dark (Fig. 3a, b). Furthermore, our custom image analysis pipeline does not require fluorescence labeling of the cell periphery for cell-size quantification, thus reducing stress imposed on cells during imaging. In future experiments, our device would also permit the localization of fluorescently tagged proteins in concert with bright-field imaging.
Under continuous illumination, Synechocystis cells followed exponential expansion kinetics at low cell density, which was previously observed in fast-growing coccoid Staphylococcus aureus cells . On average, cell size increased slightly over time, which may be due to the transition from batch culture to a surface-associated mode of growth. Twelve-hour dark periods simply suspended growth and division (Fig. 3a, b), but did not alter exponential growth (Additional file 9: Figure S6B), generation times (Fig. 4b), sister-cell generation time correlation (Fig. 4c and d), division symmetry (Fig. 4e, f), or cell-size control (Fig. 5) as compared to cells grown in continuous light. Cells showed no obvious signs at the gross level of growth of anticipating transitions into or out of the dark periods, even after three dark phases. The rapid cessation and resumption of growth when transitioning from light to dark and vice versa, respectively, suggest that light affects biomass accumulation through rapid metabolic control rather than via changes mediated by transcriptional/translational mechanisms, which are typically on the timescale of hours.
Despite substantial variation in growth rates (~30%), sister cell generation times were strikingly similar; for some sisters, the variation in generation times was only a few percent. The positive correlation between sister generation times argues against the uneven partition of molecules (mRNA, proteins, metabolites) as the source of generation time variation because such mechanisms would yield a negative correlation between sister generation times. Sister cells with different generation times tended to result from an asymmetric division (Fig. 4c, e), suggesting that the maintenance of generation times between sisters requires similarity in cellular composition between the two sister cells produced by a symmetric division and that generation times are then determined relatively deterministically (and similarly) in the two sisters. Another study has observed a positive correlation between sister generation times in mammalian cells . One potential explanation for the high degree of correlation between sisters, as compared with that between mother and daughters (Additional file 13: Figure S8C), involves deterministic components shared by sisters that are not inherited. One study argues that an underlying nonlinear process affecting generation time would produce such a correlation, whereby cell divisions occurring during a particular phase of the nonlinear process would produce daughter cells with generation times corresponding to that inherited phase . On the other hand, mother and daughter cells are unlikely to inherit the same phase, and hence would have uncorrelated generation times. In Synechocystis cells, a likely candidate for such a nonlinear effect on generation time is the circadian cycle. Although it is possible that phases of the circadian cycle influence cell-cycle duration, resulting in the generation time patterns that we observe, we did not observe any evidence of circadian regulation in our single-cell growth data. Instead, size-based cell-cycle regulation alone tends to produce correlated sister generation times.
By comparing our data with simulations of three size-control models, we determined that compared with the sizer or timer models, Synechocystis follows more closely to an adder principle whereby a constant volume is added each cell cycle. This model explains both the strong correlation between generation times of sisters resulting from a symmetric division (Fig. 4c, d), given that their similar size implies that a similar time period is required to accumulate the appropriate volume increment, and the difference in generation times between sisters resulting from an asymmetric division (Fig. 4c, e), with the smaller of the two cells requiring longer to accumulate the volume increment during exponential growth. The adder model also recapitulates many other growth statistics better than the sizer and timer models, including sister asymmetries in both continuous illumination (Additional file 4: Table S1) and during light-dark cycles (Additional file 4: Table S2), although in some cases measurement noise precludes determination of the nature of the correlation and in other cases there were small deviations between the predicted and experimental slopes (Additional file 4: Tables S1 and S2). Molecular mechanisms underlying size regulation via any of the three models have not been determined in any bacterial species. It is possible that certain (perhaps all) species actually implement a combination of cell-size regulation methods, which are in turn controlled by translational and/or metabolic processes. It has been proposed that size control is affected by the dilution of transcription factors or the initiation of DNA replication rather than upon cell division, and that the regulated quantity is cell size per genome or replication origin rather than cell size per se [36, 37]. If this is indeed the case, the fact that Synechocystis cells are considered to be polyploid  may underlie inconsistencies between our experimental data and simulations based on the adder model.
Size and growth control are fundamental physiological features of all cells, and tools such as microfluidics and automated image analysis make possible the careful quantification of these parameters with great precision and can be combined with statistical analyses. The ability to image lineages for multiple generations, over several days, potentiates studies in other slow growing cyanobacteria such as Synechococcus to address the generality of the behaviors we have uncovered, particularly the immediate responses to changes in light. How cells respond to changes in the environment such as nutrient starvation is not generally understood, and cyanobacteria experience daily light cycles that likely require adaptation of their size and growth; single-cell imaging of such transitions can be a powerful tool to shed light on the underlying control mechanisms . Given the likely commonality of adder-based cell-size control in Synechocystis with fast-growing heterotrophs such as E. coli and eukaryotes such as S. cerevisiae, it is tempting to speculate about the generality of the adder rule in other walled organisms that exhibit size homeostasis such as the shoot apical meristem of plants (or even in wall-less eukaryotes). The diversity of mechanisms for cell-size determination and maintenance from bacteria to eukaryotes suggests that the basis of the adder principle is either a common molecular component across kingdoms, such as DNA, or is based on common physical constraints. Uncovering the molecular mechanisms of size homeostasis, particularly during environmental perturbations, represents an exciting future challenge.
The PDMS microfluidic cell culture chip was fabricated using standard multilayer soft lithography techniques  (Additional file 16: Figure S11). Briefly, molds were created using SU8 (MicroChem, Massachusetts, USA) or AZ (EMD Performance Materials, Darmstadt, Germany) photoresist. SU8 was used to achieve rectangular-channel cross sections while AZ photoresist was used to create rounded-channel cross-sections that facilitate valve formation. After mold fabrication, PDMS mixture was poured onto the top-layer mold to form a thick block and spun onto the bottom-layer mold to form a thin film. Finally, the two PDMS layers were bonded to each other and to a glass slide. The dimensions of each microfluidic chamber are 1850 × 900 μm, and the height is 25 μm (Additional file 1: Figure S1A). The heights of control and flow channels are 25 and 20 μm, respectively.
Bulk Synechocystis culturing
Synechocystis cell culture stocks were grown in glass culture tubes inside a shaking incubator at 30 °C. Tubes were secured at an angle for improved agitation. Cells were grown in BG11 medium (Sigma) under white light (80 μE m–2 s–1). Cells were not entrained beforehand.
Field emission scanning electron microscopy (FE-SEM)
Synechocystis cells were spotted on polycarbonate membrane filters, after which they were fixed and gently processed using a Vacuum Microfiltration apparatus. Disks were sputter-coated with gold palladium and visualization was performed on a Zeiss Sigma FE-SEM. Images were processed using a noise-reduction algorithm.
Microfluidic-based Synechocystis growth
Before each experiment, the microfluidic chip was flushed with water to passivate uncured PDMS bonds. Flow channels were treated with 2% Pluronic F-127 solution (Sigma-Aldrich, Missouri, USA) for 1 h to prevent cells from sticking to the surfaces during the loading process, which would cause blockage. Pluronic was used to passivate only channel walls and was rinsed away prior to cell loading; the microfluidic chambers used for culturing cells were not treated with Pluronic solution. Chambers were treated with 0.01% polylysine (PLL; Sigma-Aldrich, Missouri, USA) for 2 h to promote attachment of cells to PDMS surfaces and to facilitate imaging . PLL was only used to treat PDMS culturing chambers and was washed away after 2 h of incubation. Due to the push-up valve structure of the microfluidic device, culturing chambers have PDMS surfaces on all sides, which likely also results in less PLL adsorption compared with a glass surface.
Exponentially growing cells cultured in a flask under continuous illumination were seeded into the chambers at low cell density with each field of view initially containing approximately 50–100 cells, so that cells would experience an approximately constant environment for a long period of time; indeed, we observed that growth rate remained constant for the first ~30 h (Fig. 1c). Even though the capacity of each chamber is approximately 40 nL, only the center portion (400 × 400 μm) was imaged at the higher magnification (20×) used to image bacterial cells. Cultures were not entrained before introduction into the microfluidic device. During all experiments, 5% CO2 was continuously added into the microscope’s environmental chamber through a water bubbler, which helped to maintain humidity and CO2 concentration. For light-dark cycle experiments, LEDs were switched on and off every 12 h. Toward the end of the experiments, cell growth, as measured by combining all lineages, slowed down (Fig. 1c), but growth was not significantly affected within the cell cycle of individual cells (Fig. 2e). The slowdown in growth is likely due to nutrient depletion; we did not continuously flush the channels with new medium to avoid washing out cells from the chambers, which would have made it challenging to track cells for several generations. Interestingly, we found that sister cell generation times were highly correlated in both the first and last 30 h of the experiment (Additional file 13: Figure S8A, B), indicating that the mechanism driving this correlation was unaffected by the slowdown in growth.
Image acquisition and analysis
Bright-field images were acquired every 10 min from 20 chambers of the microfluidic chip with 20× magnification (Leica DMI 6000, Illinois, USA). The image exposure time (16 ms) did not affect cell growth, as demonstrated by the lack of growth during the dark phase of light-dark experiments (Fig. 3a, b). A custom MATLAB-based script coordinated all components of the cell culture system including the microscope, stage, camera, and LEDs. Every hour, contrast-based autofocusing was performed to correct for vertical stage drift. Automated image analysis was performed with custom MATLAB scripts. Details of acquisition and analysis methods can be found in the Additional file 4 .
Cell volume calculation
To reconstruct the volume of the cell from cell contours extracted from time-lapse images, we assumed rotational symmetry with respect to the major axis of the cell or doublet (Additional file 5: Figure S4). For each cell, we first extracted the orientation of the major axis. Total cell volume was then calculated by summing the volumes of a series of circular disks with radii equal to the perpendicular distances between the major axis and the cell boundary.
Extraction of growth parameters from single-cell growth curves
For single-cell growth curves obtained from experiments under continuous illumination, we extract birth size, division size, and generation time following Additional file 12: Figure S7. In order to analyze single-cell growth curves obtained from experiments under light-dark cycles, we removed segments corresponding to dark periods of the experiment. The fact that neither growth nor division is observed in dark allowed us to piece together single-cell growth curve segments from only the illuminated periods and still obtain continuous growth curves. In this case, generation time corresponds to the total amount of time a cell spends under light.
To determine how much of the observed correlation between sister cell generation times (Fig. 4c, d) could be attributed to their time of birth, we bootstrapped the generation-time differences, randomizing generation times based on cells having either similar birth times or similar division times (as determined by binning with a 5-h bin width). The standard deviation of unrandomized sister cell generation time differences was 3% of their mean generation time, whereas the standard deviation was 8.3% for randomly selected cell pairs in a similar growth window. The distributions after randomizing differed significantly from our data, indicating that sister cell generation times are highly correlated.
Simulations of cell division models
Distributions of birth volumes, cell division asymmetry, and growth rates were used as inputs for simulations of the sizer, adder, and timer models of cell size regulation. We observed an increase in cell size (Fig. 2e) that occurred largely in the first full cell cycle; cell birth and cell division sizes were approximately consistent for later cell cycles (Additional file 14: Figure S9). We then directly incorporated measurements of the size increase into our simulations using the model described in Additional file 4. Details of our simulation methodology are included in Additional file 4.
The authors would like to thank Michael Robles and Sean Liu from the Stanford Microfluidics Foundry for help in device fabrication.
This work was supported by an NSF Graduate Research Fellowship and Stanford Graduate Fellowship (to FBY), a Bio-X Interdisciplinary Initiatives Program seed grant (to DB and KCH), a National Academies Keck Future Initiatives seed grant (to KCH), NSF grant EF-1038697 (to KCH), NSF CAREER Award MCB-1149328 (to KCH), and NSF grant MCB-1331151 (to DB). The funders had no role in the design of the study, nor in the collection, analysis, or interpretation of data or writing of the manuscript.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files, or is available upon request.
Conceptualization, FBY, RMWC, AZ, DB, KCH, and SRQ. Methodology, FBY, LW, RMWC, AZ; Software, FBY and LW. Investigation, FBY, LW, RMWC, and AZ. Writing, FBY, LW, MH, DB, KCH, and SRQ. Supervision, MH, DB, KCH, and SRQ. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Logan BE. Exoelectrogenic bacteria that power microbial fuel cells. Nat Rev Microbiol. 2009;7(5):375–81.View ArticlePubMedGoogle Scholar
- Chabot JR, Pedraza JM, Luitel P, van Oudenaarden A. Stochastic gene expression out-of-steady-state in the cyanobacterial circadian clock. Nature. 2007;450(7173):1249–52.View ArticlePubMedGoogle Scholar
- Nogales J, Gudmundsson S, Knight EM, Palsson BO, Thiele I. Detailing the optimality of photosynthesis in cyanobacteria through systems biology analysis. Proc Natl Acad Sci U S A. 2012;109(7):2678–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Chau RM, Ursell T, Wang S, Huang KC, Bhaya D. Maintenance of motility bias during cyanobacterial phototaxis. Biophys J. 2015;108(7):1623–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Ng WO, Grossman AR, Bhaya D. Multiple light inputs control phototaxis in Synechocystis sp. strain PCC6803. J Bacteriol. 2003;185(5):1599–607.View ArticlePubMedPubMed CentralGoogle Scholar
- Ursell T, Chau RM, Wisen S, Bhaya D, Huang KC. Motility enhancement through surface modification is sufficient for cyanobacterial community organization during phototaxis. PLoS Comput Biol. 2013;9(9):e1003205.View ArticlePubMedPubMed CentralGoogle Scholar
- Dorrich AK, Mitschke J, Siadat O, Wilde A. Deletion of the Synechocystis sp. PCC 6803 kaiAB1C1 gene cluster causes impaired cell growth under light-dark conditions. Microbiology. 2014;160(Pt 11):2538–50.View ArticlePubMedGoogle Scholar
- Wiegard A, Dorrich AK, Deinzer HT, Beck C, Wilde A, Holtzendorff J, et al. Biochemical analysis of three putative KaiC clock proteins from Synechocystis sp. PCC 6803 suggests their functional divergence. Microbiology. 2013;159(Pt 5):948–58.View ArticlePubMedGoogle Scholar
- Panoff J-M, Priem B, Morvan H, Joset F. Sulphated exopolysaccharides produced by two unicellular strains of cyanobacteria, Synechocystis PCC 6803 and 6714. Arch Microbiol. 1988;150(6):558–63.View ArticleGoogle Scholar
- Wang P, Robert L, Pelletier J, Dang WL, Taddei F, Wright A, et al. Robust growth of Escherichia coli. Curr Biol. 2010;20(12):1099–103.View ArticlePubMedPubMed CentralGoogle Scholar
- Campos M, Surovtsev IV, Kato S, Paintdakhi A, Beltran B, Ebmeier SE, et al. A constant size extension drives bacterial cell size homeostasis. Cell. 2014;159(6):1433–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Jones LJ, Carballido-Lopez R, Errington J. Control of cell shape in bacteria: helical, actin-like filaments in Bacillus subtilis. Cell. 2001;104(6):913–22.View ArticlePubMedGoogle Scholar
- Connell JL, Wessel AK, Parsek MR, Ellington AD, Whiteley M, Shear JB. Probing prokaryotic social behaviors with bacterial “lobster traps”. MBio. 2010;1(4):e00202–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Long Z, Nugent E, Javer A, Cicuta P, Sclavi B, Cosentino Lagomarsino M, et al. Microfluidic chemostat for measuring single cell dynamics in bacteria. Lab Chip. 2013;13(5):947–54.View ArticlePubMedGoogle Scholar
- Moffitt JR, Lee JB, Cluzel P. The single-cell chemostat: an agarose-based, microfluidic device for high-throughput, single-cell studies of bacteria and bacterial communities. Lab Chip. 2012;12(8):1487–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Gomez-Sjoberg R, Leyrat AA, Pirone DM, Chen CS, Quake SR. Versatile, fully automated, microfluidic cell culture system. Anal Chem. 2007;79(22):8557–63.View ArticlePubMedGoogle Scholar
- Taheri-Araghi S, Brown SD, Sauls JT, McIntosh DB, Jun S. Single-cell physiology. Annu Rev Biophys. 2015;44:123–42.View ArticlePubMedGoogle Scholar
- Mazouni K, Domain F, Cassier-Chauvat C, Chauvat F. Molecular analysis of the key cytokinetic components of cyanobacteria: FtsZ. ZipN and MinCDE Mol Microbiol. 2004;52(4):1145–58.View ArticlePubMedGoogle Scholar
- Tzagoloff H, Novick R. Geometry of cell division in Staphylococcus aureus. J Bacteriol. 1977;129(1):343–50.PubMedPubMed CentralGoogle Scholar
- Westling-Haggstrom B, Elmros T, Normark S, Winblad B. Growth pattern and cell division in Neisseria gonorrhoeae. J Bacteriol. 1977;129(1):333–42.PubMedPubMed CentralGoogle Scholar
- Pattanayak G, Rust MJ. The cyanobacterial clock and metabolism. Curr Opin Microbiol. 2014;18:90–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Johnson CH, Stewart PL, Egli M. The cyanobacterial circadian system: from biophysics to bioevolution. Annu Rev Biophys. 2011;40:143–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Taheri-Araghi S, Bradde S, Sauls JT, Hill NS, Levin PA, Paulsson J, et al. Cell-size control and homeostasis in bacteria. Curr Biol. 2015;25(3):385–91.View ArticlePubMedGoogle Scholar
- Osella M, Nugent E, Cosentino LM. Concerted control of Escherichia coli cell division. Proc Natl Acad Sci U S A. 2014;111(9):3431–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Amir A. Cell size regulation in bacteria. Phys Rev Lett. 2014;112(20):208102.View ArticleGoogle Scholar
- Ghusinga KR, Vargas-Garcia CA, Singh A. A mechanistic stochastic framework for regulating bacterial cell division. Sci Rep. 2016;6:30229.View ArticlePubMedPubMed CentralGoogle Scholar
- Sauls JT, Li D, Jun S. Adder and a coarse-grained approach to cell size homeostasis in bacteria. Curr Opin Cell Biol. 2016;38:38–44.View ArticlePubMedGoogle Scholar
- Soifer I, Robert L, Amir A. Single-cell analysis of growth in budding yeast and bacteria reveals a common size regulation strategy. Curr Biol. 2016;26(3):356–61.View ArticlePubMedGoogle Scholar
- Tanouchi Y, Pai A, Park H, Huang SQ, Stamatov R, Buchler NE, et al. A noisy linear map underlies oscillations in cell size and gene expression in bacteria. Nature. 2015;523(7560):357–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Q, Pando BF, Dong G, Golden SS, van Oudenaarden A. Circadian gating of the cell cycle revealed in single cyanobacterial cells. Science. 2010;327(5972):1522–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Pando BF, van Oudenaarden A. Coupling cellular oscillators--circadian and cell division cycles in cyanobacteria. Curr Opin Genet Dev. 2010;20(6):613–8.View ArticlePubMedGoogle Scholar
- Teng SW, Mukherji S, Moffitt JR, de Buyl S, O’Shea EK. Robust circadian oscillations in growing cyanobacteria require transcriptional feedback. Science. 2013;340(6133):737–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Pattanayak GK, Lambert G, Bernat K, Rust MJ. Controlling the cyanobacterial clock by synthetically rewiring metabolism. Cell Rep. 2015;13(11):2362–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou X, Halladin DK, Rojas ER, Koslover EF, Lee TK, Huang KC, et al. Mechanical crack propagation drives millisecond daughter cell separation in Staphylococcus aureus. Science. 2015;348(6234):574–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Sandler O, Mizrahi SP, Weiss N, Agam O, Simon I, Balaban NQ. Lineage correlations of single cell division time as a probe of cell-cycle dynamics. Nature. 2015;519(7544):468–71.View ArticlePubMedGoogle Scholar
- Ho PY, Amir A. Simultaneous regulation of cell size and chromosome replication in bacteria. Front Microbiol. 2015;6:662.View ArticlePubMedPubMed CentralGoogle Scholar
- Taheri-Araghi S. Self-consistent examination of Donachie’s constant initiation size at the single-cell level. Front Microbiol. 2015;6:1349.View ArticlePubMedPubMed CentralGoogle Scholar
- Griese M, Lange C, Soppa J. Ploidy in cyanobacteria. FEMS Microbiol Lett. 2011;323(2):124–31.View ArticlePubMedGoogle Scholar
- Sutterlin HA, Shi H, May KL, Miguel A, Khare S, Huang KC, et al. Disruption of lipid homeostasis in the Gram-negative cell envelope activates a novel cell death pathway. Proc Natl Acad Sci U S A. 2016;113(11):E1565–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou J, Ellis AV, Voelcker NH. Recent developments in PDMS surface modification for microfluidic devices. Electrophoresis. 2010;31(1):2–16.View ArticlePubMedGoogle Scholar
- Donoser M, Bischof H, editors. Efficient Maximally Stable Extremal Region (MSER) Tracking. Computer Vision and Pattern Recognition, 2006, IEEE Computer Society Conference on 17–22 June 2006. doi:10.1109/CVPR.2006.107.