|Home | About | Journals | Submit | Contact Us | Français|
Cell cycle distribution of adherent cells is typically assessed using flow cytometry, which precludes the measurements of many cell properties and their cycle phase in the same environment. Here we develop and validate a microscopy system to quantitatively analyze the cell-cycle phase of thousands of adherent cells and their associated cell properties simultaneously. This assay demonstrates that population-averaged cell phenotypes can be written as a linear combination of cell-cycle fractions and phase-dependent phenotypes. By perturbing cell cycle through inhibition of cell-cycle regulators or changing nuclear morphology by depletion of structural proteins, our results reveal that cell cycle regulators and structural proteins can significantly interfere with each other’s prima facie functions. This study introduces a high-throughput method to simultaneously measure cell cycle and phenotypes at single-cell resolution, which reveals a complex functional interplay between cell cycle and cell phenotypes.
The cell cycle is the series of highly regulated steps that lead to controlled cell division. Typically, cells first prepare for DNA synthesis (G1 phase), replicate their DNA (S phase), prepare for mitosis (G2 phase), and undergo mitosis (M phase).1, 2 During this cell cycle, specific proteins serve as door guards at every phase to prevent cells from early entrance into the next stage of cell cycle.3 Misregulation of cell cycle in human and rodent cells has been implicated in a number of disease states.4, 5, 6 For example, mutated p53 causes cells to lose the function of the G1/S checkpoint, replicating defective DNA, and finally leading to cancer.4, 6 Flow cytometry (FC) is the instrument of predilection to measure cell-cycle distribution, particularly of adherent cells, and the effects of drug treatment or genetic alteration (knockdown, knockout, over-expression, etc.) on cell cycle.7, 8 A major advantage of FC is its ability to analyze a large number of cells in a short time. However, conventional FC analysis requires cells to be detached from their substrate and therefore cannot measure cell properties (e.g. nuclear shape, cell migration, cytoskeleton organization, etc.) at the same time in the same environment. Moreover, since the expression of a wide range of proteins greatly vary during cell cycle,9, 10, 11, 12 these cell properties may adopt significantly different values in different phases. Consequently, without simultaneous measurement of cell cycle phase and cell properties in the same cells, an observed change in cell properties following a forced change in protein expression does not necessarily mean that this protein is a regulator of the cell property of interest. Rather this protein could be a cell cycle regulator (Fig. 1A).
Here we use a microscope-based assay to measure both the cell cycle phase of a thousand of individual adherent cells and their associated cellular and nuclear properties rapidly and simultaneously. This assay demonstrates that population-averaged cell morphological properties strongly depend on cell-cycle phase and could be written as linear combinations of cell-cycle fractions and phase-dependent morphological properties. This assay reveals that key structural nuclear-envelope proteins (Nesprins, Lamin A/C) are regulators of nuclear size and nuclear shape partially because they affect cell cycle distribution; they are not bona fide (intrinsic)regulators of nuclear morphology.13, 14, 15 Vice versa, this assay indicates that inhibition of cell cycle regulator cyclin-dependent kinase 4/6, cdk4/6, are also nuclear morphological regulators, and that some commonly used cell cycle synchronization methods have significant and lasting effects on cell and nuclear morphology.16, 17
First we established and validated a microscopy-based method to measure the cell-cycle phase in individual adherent cells. Mouse myoblasts (c2c12) were plated on a glass-bottom dish. After 48h incubation, cells were fixed, permeabilized, and incubated with nuclear DNA stain Hoechst33342, a dye routinely used in FC analysis of cell cycle distribution,7, 18 as well as additional stains to detect actin filament structures at the cell cortex and delineate the cytoplasm for cell edge detection (Fig. 1B). The dish was placed on a customized scanning fluorescence light microscope equipped with a motorized stage (Fig. 1B). Importantly, to quantitatively determine cell cycle phases at single-cell resolution, we introduced a calibration step to correlate measured Hoechst light intensity in the nucleus to DNA content (see Materials and Methods). Light intensities from solutions of Hoechst molecules were measured using the same customized scanning microscope used subsequently to measure cell properties in the same cells (Fig. 1B). For each dish, reference illumination and dark images were collected and used to normalize the intensities in the illumination field (see additional details under Methods). This calibration step addressed potential non-uniform illumination of the samples by the light microscope and the non-uniform recording by the optical train of the microscope and the CCD camera.19
The cell cycle-phase distribution of c2c12 mouse myoblasts measured by this microscopy-based assay was carefully compared to the cycle-phase distribution obtained from conventional FC analysis (Fig. 1, C and D).7 For both methods, cells were fixed and their DNA was stained with the same nuclear DNA dye Hoechst 33342. Even though different numbers of cells were assessed (~10,000 cells for FC vs. 1,000–2,000 for our assay), cell cycle distributions were statistically indistinguishable (Fig 1, C and D). By zooming in on the G0/G1 phase peak, we found that the coefficients of variation (CV) of these cell sub-populations using Gaussian fits were close, and that the signal-to-noise ratios (1/CV) were highly similar (Fig 1D and Table 1).20
A direct comparison between FC analysis and the analysis provided by our assay also showed statistically indistinguishable cell cycle distributions for human breast carcinoma cells (MD-MB-231) (Fig. S1). Our assay was not limited to c2c12 and MD-MB-231cell lines. We tested seven additional types of adherent cells and found that our assay could readily measure their cell-cycle distributions (Fig. S2). These cells were chosen to determine whether our assay worked on both normal and diseases cells, immortalized and primary cells, as well as human and rodent cells. The tested cells include human pancreatic normal epithelial cells, patient-derived pancreatic cancer ductal adenocarcinoma cells harvested from the primary tumor and from liver metastatic sites, normal human breast epithelial cells (MCF10A), and primary mouse embryonic fibroblasts (Fig. S2).
Next, we used the same microscopy method to measure both cell-cycle phase and cell and nuclear properties simultaneously in the same adherent cells. Flow cytometry can simultaneously measure cell size and DNA content in suspension cells. However, other important morphological properties of adherent cells including cell shape and nuclear shape become effectively meaningless if adherent cells are (artificially) detached from their substrate. Moreover, spatial information is lost. We measured DNA content of cells in culture and color-coded nuclei according to the associated cell cycle phase (Fig 2A and B). In the same cells, we measured nuclear size and shape and cellular size and shape to determine not only population-averaged values of these properties, but also their phase-dependent values (Fig. 2, C–F). As would be expected, cell and nuclear size increased as cell cycled from the G0/G1, S, ad G2/M phases (Fig. 2D, F, H, and J). More unexpectedly, we found that nuclei became more elongated during cycle progression, an information lost using conventional approaches. Together these results show that our microscopy assay can measure the following six independent quantities simultaneously: the three percentages of cells in each phase of the cell cycle (fG0/G1, fS, and fG2/M; Fig. 1 and Fig. 2B) and the three associated values of the cell property x (e.g., cell shape, nuclear shape, etc.) in each phase, xG0/G1, xS, and xG2/M, which are typically different from each other.
To help clarify the potential influence of cell cycle-phase distribution on cell properties (such as nuclear size or cell shape),21 we wrote the mean value of a given cell property <x> as a linear combination of cell fractions in each phase and mean values of the cell property in each of the cell cycle phase:
Here, <x> is the cell population-averaged value of the cell property being considered (e.g. the cell population-averaged values of nuclear and cell sizes and shapes), xi are the mean values of this property in the cell-cycle phases i (i = G0/G1, S, and G2/M phases), and fi are the fractions of cells in each phase i, which is also the relative time cells spend in each phase of the cell cycle. The method presented here enables us to measure xi and fi separately and simultaneously in the same cells.
When assessing the role of the expression or activity of a protein in a given cell function, cells are typically subjected to a drug that specifically inhibits/activates the protein or the gene of interest is knocked down (KD), knocked out (KO) or over-expressed. It is then pervasively assumed than any measured change in mean cell property (i.e. a change in the population averaged value <x>) reflects the involvement of this protein in the cell function being assessed, without considering the possible redistribution of the cells along the cell cycle. For example, the nuclear lamina-associated structural proteins Nesprins are widely believed to be structural regulators of nuclear shape because of their role of physically connecting the nuclear lamina to the cytoskeleton and because the depletion of Nesprins changes nuclear roundness.13, 22, 23, 24
This approach to assess the role of a protein in cell functions is only legitimate if the fractions fi of cells in the different cell-cycle phases remain unchanged following application of the inhibitor/activator or genetic manipulation, i.e. that the protein of interest is not also a cell cycle regulator. Alternatively, when in doubt that it is actually correct and that cell properties could be cell-cycle dependent, then cells can be synchronized. To synchronize cells, cells are often subjected to serum-starvation or chemicals that arrest cells in a specific phase of the cell cycle (Fig. 3A).8, 25, 26, 27 While used routinely in cell biology, such chemical treatments can induce their own side effects, including the fact that they may not synchronize cells. For instance, thymidine arrests cells in the early S phase. However, thymidine treatment to assess cell phenotypes in the S phase is only correct if the cell properties, xi, are the same for S-phase-synchronized cells as the cell properties of non-synchronized cells that are happen to be in the S phase. The same assumption is typically made when cell synchronization is induced by serum starvation or cell treatment with the microtubule depolymerizing drug nocodazole, which enrich the G0/G1 phase and M phase, respectively.
Exploiting the advantage that our assay can measure cycle phase and cell properties in the same individual cells simultaneously, we tested whether these commonly used cell synchronization methods verified this assumption. Surprisingly, synchronization methods, such as serum starvation and thymidine/nocodazole treatments, affected not only cell cycle distributions, as expected, but also greatly changed cell and nuclear properties (Fig 3, B–F). Results from our analysis confirmed that, as expected, these methods of synchronization did enrich cells in the target cycle phase: the S phase for thymidine treatment (orange), the G1 phase following serum starvation (green), and the M phase for nocodazole treatment (purple, Fig. 3, B and C). However, these synchronization methods also had a significant and lasting effect on cell properties in every cell cycle phase (Fig 3, D–F). For instance, thymidine treatment greatly increased the size of cells in the enriched S phase by 120% compared to untreated cells in the same S phase, i.e. the mean size of treated cells in the S phase was 2.2 times the mean size of untreated cells in the S phase. The mean size of thymidine-treated cells in the enriched G0/G1and G2/M phases were similarly increased by 130% and 220% compared to untreated (asynchronized) cells in the same G0/G1 and G2/M phases, respectively. Serum starvation and nocodazole treatment fared somewhat better, but still changed measured cell properties by >30%, including a 75% increase in cell size for nocodazole-treated cells in the G2/M phase. If synchronization were non-invasive, these phenotypic changes should have been vanishingly small.
More importantly, we verified that adding back serum to serum-starved cells or washing cells to eliminate the synchronizing drug from the medium allowed cells to recover the cycle-phase distribution of untreated cells prior to forced synchronization within 24 h (Fig S3, A and B). However, changes in cell and nuclear properties did not recover their untreated values, i.e. all tested synchronization methods continued to affect cell/nuclear size and cell/nuclear shape, even 24h after release (Fig S3, C and D). Together these results indicate that commonly used synchronization methods are not appropriate to reduce or eliminate possible effects of cell cycle re-distribution on cell properties, as they do not keep cell properties constant in the enriched phase compared to the value of the property in untreated cells in that phase. One should use synchronization methods with caution, and should not ignore the effect of changes in phenotypic property associated with forced synchronization.
Next, we asked whether inhibition of well-characterized cell cycle regulators, such as Cdk4/6, could also inadvertently affect cellular phenotypes. Eq. 1 suggests that a change in the cell-population-averaged value of the cell phenotype <x> (say nuclear size) following protein inhibition could result from cell-cycle redistribution (i.e. changes in fi). As expected, treatment of cells with specific Cdk4/6 inhibitor IV, trans-4-((6-(ethylamino)-2-((1-(phenylmethyl)-1H-indol-5-yl)amino)-4-pyrimidinyl)amino)-cyclohexanol (CINK4), blocked cells in the G0/G1 phase and reduced the population of cells in the S and G2/M phase (Fig. 4A).28 Moreover, Cdk4/6 inhibition also significantly affected population-averaged nuclear and cell morphology (e.g. ~50% decrease in nuclear size accompanied by rounding of the nucleus) (Fig. 4, B–D), and continued to affect cell properties up to 6h after release (Fig. S3, E and F). However, our assay which can evaluate cell-cycle-dependent values of phenotypes simultaneously (i.e. changes in xi) revealed that inhibition of Cdk4/6 caused significant effects on the values of phenotypes in each phase. For instance, following Cdk4/6 inhibition, the nuclear size of cells in the G0/G1 phase decreased by >50% (P < 0.0001). Hence our assay suggests that besides its well-known central role in cell cycle, Cdk4/6 may also play a structural role by regulating nuclear morphology.
Having manipulated cell cycle distributions through various forced synchronization schemes or Cdk4/6 inhibition to determine their (unexpected) effect on cellular and nuclear properties in each phase, we next manipulated cell properties through shRNA-mediated depletion of seemingly well-established structural regulators of these properties, while simultaneously measuring potential changes in cell cycle distribution. To help clarify this test, we return to Eq. 1. For a protein to be an “intrinsic” regulator of cell/nuclear morphology, that protein should not have the side effect of significantly affecting cell cycle distribution (which in turn would affect cell properties, Fig. 1A): the depletion of that protein should mostly induce changes in cell properties xi, not changes in cell cycle fractions fi.
Nuclear lamina protein Lamin A/C, which forms a thin elastic filamentous meshwork underneath the nuclear envelope and penetrates the intranuclear space, and LINC complex molecules Nesprin2giant and Nesprin3, which physically connect the nuclear lamina to the cytoskeleton, have recently been established as major structural proteins that regulate nuclear size and nuclear shape.22, 29, 30, 31 Indeed, we found that shRNA-mediated depletion of Lamin A/C, Nesprin2giant and Nesprin3 all significantly deformed and decreased the size of the nucleus and cell in all phases of the cell cycle (“average bars” in Fig. 5, C–E).
However, the depletion of these nuclear envelope proteins also affected cycle-phase distribution (Fig. 5, A and B).14, 32, 33 For instance, shRNA-mediated depletion of Lamin A/C significantly enriched the G0/G1 phase, while reducing the fractions of cells in the S and G2/M phases (Fig. 4B). This result was confirmed in primary LMNA+/+ and LMNA−/− mouse embryonic fibroblasts (Fig S4): Lamin A/C deficiency enriched the G0/G1 phase, while reducing the number of cells in other phases. Depletion of Nesprin3 also significantly enriched the G0/G1 phase and slightly, but significantly, reduced the fractions of cells in the S and G2/M phases. Nesprin 3 depletion had a much larger effect on cell cycle redistribution than Nesprin2giant (Fig. 5B).
Our assay allowed us to assess changes in nuclear size, phase by phase (for example, by comparing the mean values of nuclear size in the G2/M phase of both control cells and shRNA-depleted cells directly) (Fig. 5, D–I). Furthermore, the assay revealed that the role of Lamin A/C and Nesprins in nuclear morphology and cell size strongly depended on cell-cycle phase (Fig. 5, C–E). For instance, depletion of Lamin A/C and Nesprins had a 3-fold greater effect on nuclear size for cells in the G0/G1 phase than cells in the G2/M phase (Fig. 5C). Together these results showed that while Lamin A/C, Nesprin2giant and Nesprin3 should be considered both regulators of nuclear morphology and cell cycle phase distribution, which indirectly affect nuclear morphology.
When cells are subjected to nocodazole or thymidine treatments, serum starvation, or shRNA-depletion of structural proteins such as Lamin A/C and Nesprins, the population-averaged mean values of nuclear size, cell size, and nuclear shape change (Figs. 3 and and5).5). We can quantify the contributions to changes in these mean values due to: (i) direct changes in intrinsic values of these properties (“phenotypic effect”, Fig. 6A), i.e. changes in the cell-cycle-independent values of these properties, (ii) indirect changes due to cell cycle re-distribution (“cell cycle effect”, Fig. 6A), and (iii) second-order effects due to coupled phenotypic changes and cell cycle distribution (“2nd order term”, Fig. 6A). Indeed, one can easily show that the global change in a cell property x can be written as:
Here Δx is the total change in the mean value of nuclear size, cell size, and nuclear shape caused by shRNA-mediated depletion (subscript KD) of Nesprin2giant, Nesprin3, or LaminA/C or by pharmacological synchronization or protein inhibition (subscript Drug) compared to control cells; Δxi are the same differences but evaluated for cells in each cell-cycle phase i; and Δfi are the changes in cell-cycle fractions for each phase i. Therefore, overall changes in cell properties following protein depletion or drug treatment may stem from three distinct contributions: (i) intrinsic changes in cell properties independent of changes in cell cycle-phase distribution (first term in Eq 2), indirect changes in cell properties due a cell cycle-phase redistribution (second term), and coupled changes in cell cycle and cell properties (third term), which are expected to be second-order in magnitude (Fig. 6A). For cells depleted of structural proteins such as nuclear-envelope-associated Nesprins, changes in cell properties should mostly be due to the first term, which means that Δfi ≈ 0 in Eq. 2 and, in turn, Δx ≈ ΣΔxificontrol. For cells subjected to synchronization by serum-starvation or drug treatment, changes in cell properties should mostly be due to the second term (i.e. due to cell-cycle redistribution), which means that Δxi ≈ 0 in Eq. 2 and Δx ≈ ΣΔfixicontrol. Here ficontrol and xicontrol are the fractions of cells and values of the cell property of interest in phases i for control cells.
Here, we found that overall changes in nuclear size upon depletion of Lamin A/C due to direct intrinsic changes in nucleus size (when properly compared phase-to-phase) was only 75%, and that the contribution due to cell-cycle redistribution was significant (>20%) (Fig. 6B). Depletion of Nesprin2giant and Nesprin3 still induced 9% contribution to global change in nuclear size from cell cycle re-distribution (Fig. 6B). Similar conclusions held for changes in cell size (Fig. 6C) and nuclear shape (Fig. 6D): changes in cell size and nuclear shape by depletion of Lamin A/C, Nesprin3, or Nesprin2giant were partly due to non-negligible contributions from cell cycle re-distribution.
Vice versa, we analyzed the contributions of (unwelcome) changes in cell/nuclear properties following forced synchronization by serum starvation, thymidine treatment, or nocodazole treatment (Fig. 6, B–D). Ideally, if synchronization were non-invasive, then the contribution from the first term in Eq. 2 to changes in nuclear shape and size would be 100% and the contribution from the other terms would be zero. Remarkably, we found that the contributions to changes in nuclear size due to forced cell-cycle redistribution was significantly smaller than the contributions from intrinsic changes in nuclear size: only 31% for serum starvation, 9% for thymidine treatment, and 35% for nocodazole treatment were due to changes in nuclear size induced by forced synchronization compared to contributions of 68%, 91%, and 41% that were due to (uncontrolled) changes in nuclear size in each phase (Fig. 6B). Similar conclusions held when assessing changes in nuclear size and nuclear shape (Fig. 6, C and D).
A literature survey indicates that the vast majority of biological studies that make use of shRNA-induced depletion or genetic knockout or over-expression of specific proteins as methods to assess the function of a protein in mammalian cells only use changes in population-averaged values of the cell property under study (the left hand-side of Eq. 2) as a way to quantitatively assess a protein function in cell physiology, without taking into account possible changes in phase distribution caused by these cell manipulations. The remainder of these studies used forced synchronization as a way to reduce cell-cycle-dependent effects, and none measured both changes in cell properties and cell cycle distribution at the same time, as presented in this study. Our results suggest that this general approach to assess the role of proteins in establishing cell and nuclear properties could lead to erroneous conclusions.
Here, we used a microscopy-based assay that allowed us to measure cell cycle phase and cell/nuclear properties simultaneously in the same cells. This assay: (i) distinguished bona fide cell-cycle-independent regulators of cell and nuclear properties from cell-cycle-dependent regulators that regulate cell properties only or partially though cell cycle redistribution, (ii) revealed new regulators of cell cycle distribution (e.g. Nesprin2giant and Nesprin3 were not known to be cell cycle regulators), and (iii) quantified the distinct contributions of specific proteins to cell properties due to direct/intrinsic regulation by these proteins and due to indirect changes caused by uncontrolled cell cycle redistribution.
This study illuminates the critical importance of measuring cell and nuclear properties in each cell phase, highlights and quantifies the danger of using commonly used synchronization methods to eliminate potential effects of cell cycle redistribution on cell phenotypes, and suggests that regulators of other cell functions (e.g. nuclear morphology, cell motility), previously identified through shRNA-mediated depletion or knockout studies, may indeed be mixed cell-property/cell-cycle regulators.
Mouse myoblasts (c2c12) were cultured in Dulbecco’s Modification of Eagle’s Medium (DMEM) (Invitrogen, Carlsbad, CA) supplemented with 10% of fetal bovine serum (FBS, Hyclone, Logan, UT) and 100U penicillin and 100μg streptomycin (Sigma, St. Louis, MO) and 0.1% of puromcyin. Human breast cancer cells, MDA-MB-231 (ATCC, Manasas, VA), and freshly harvested LMNA+/+ and LMNA−/− mouse embryonic fibroblasts were cultured in DMEM with 10% of FBS and 100U penicillin and 100μg streptomycin. Human transformed epithelial cells, MCF-10A, were cultured in DMEM F-12 medium (Invitrogen) with 5% of horse serum, 0.5 μg/ml of hydrocortisone, 20 ng/ml of hEGF, 10 μg/ml of bovine insulin (Sigma), and 100μg/ml cholera toxin (Sigma). Human pancreatic normal epithelial cells (HPNE) were cultured in DMEM (low glucose), M3 base medium, FBS, Gentamicin, EGF, P/S (Sigma). Patient-derived pancreatic ductal adenocarcinomas, liver metastasis and primary tumor cells were cultured in DMEM with 10% of FBS and 100U penicillin and 100μg streptomycin. The culture environment was maintained at 37° C and 5% CO2. Cells were passaged every three days.
For each type of cell, ~10,000 cells were plated on a glass bottom dish (World Precision Instruments, Sarasota, FL). After incubation for 48h, the treatment for synchronization was applied. To serum-starve cells, cells were washed three times with Hanks’ balanced salt solution (HBSS) (GBICO) and cultured in serum-free medium for 72h. For synchronization at early S phase, cells were treated with 2 mM of thymidine for 18h twice. Between these two treatments, cells were rinsed with HBSS three times and cultured in normal growth medium to release cell cycle for 9h. To synchronize cells in the M phase, cells were first treated with 2 mM of thymidine for 24h. Followed by rinsing with HBSS three times, cells were released with normal growth medium for 3h. After the short release, cells were treated with 100 ng/ml of microtubule-depolymerizing drug nocodazole for 12h. For inhibition of Cdk4/6, cells were treated with 10 nM trans-4-((6-(ethylamino)-2-((1-(phenylmethyl)-1H-indol-5-yl)amino)-4-pyrimidinyl)amino)-cyclohexanol (CINK4) Cdk4/6 inhibitor (EMD Chemical) for 24h, prior to cell cycle and phenotypic analysis.
Our protocol follows that from Current Protocol from Cytometry (Wiley).18 Briefly, cells were grown in a 10-cm dish for 48h, then trypsinized, spin down at 1,000 rpm for 5 min, and resuspended in fresh culture medium without any serum. 106 cells/ml were resuspended with a final concentration 5 μg/ml of Hoechst 33342 and incubated for 30 min at room temperature. After incubation, cells were spin down to remove the dye solution, resuspended in fresh culture medium, and then promptly subjected to flow cytometry. The raw data was extracted with FlowJo software (Tree Star Inc., Ashland, OR) and processed with Matlab.
Cells were fixed with 3.7% of formaldehyde (Sigma) for 15 min at room temperature (RT). After fixation, cells were permeabilized with 0.1% Triton X-100 (Sigma) for 10 min and blocked for nonspecific binding with phosphate-buffered saline (PBS) supplemented with 10% of goat serum for 30 min. Cells were incubated with specific concentration of dye diluted from stock solution for 1h. Nuclear DNA was stained with Hoechst 33342 (Sigma) at 1:40 dilution. Cytoplasm was stained with HCS CellMask Cy5 (Invitrogen) at 1:20000 dilution. Actin was stained with Phalloidin 488 (Invitrogen) at 1:40 dilution. PBS rinse was conducted three times between each step.
Fluorescent images were collected with a Nikon DS-QiMc camera installed on a customized Nikon TE300 microscope with a 10x Plan Fluor lens (N.A. 0.3, Nikon Melville, NY), a motorized stage, and motorized excitation and emission filters (Prior Scientific, Rockland, MA) controlled by Nikon NIS Elements. Eighty-one (9-by-9 grid) fields of views were consistently generated using the software Nikon NIS-Elements. The size of the image acquired from the camera was 1280 × 1024 pixels, and the pixel size was 0.57 μm for a 10x objective. Our ΔX and ΔY for image capture was 662 μm and 520 μm respectively to allow 10% overlay between adjacent fields. The total size of the scanning region (81 images) was ~28.5 mm2. Four different channels (UV, GFP, Cy5, and phase contrast) were collected for every field of view. The calibration glass-bottom dish contained dyes for three different fluorescent channels, UV, GFP, and Cy5. For each fluorescence channel, UV, GFP, and Cy5, two calibration images were acquired: with and without illumination. The calibration images were used to reduce the non-uniform illumination of fluorescent images. All of the phenotypic information was calculated with a custom high throughput program (written in Matlab) developed in our laboratory. Typically > 1,000 cells were analyzed per different biological repeat (i.e. different cells) for a total >3,000 cells for each tested condition.
To precisely segment individual cells and nuclei, we used slightly different approaches based on the same principle. For nuclear segmentation, because of the relatively circular shape and relatively even intensity of the Hoechst stain, we filtered calibrated images (as described in the previous section) with a 23 × 23 pixel normalized Gaussian filter (similar scale as the size of nuclei) and an averaging filter (same size) to obtain IG (Gaussian intensity) and IM (averaged intensity). Subtracting IM from IG gives IN, the nuclear intensity values without regional background. Empirical testing showed that a threshold setting of 10 was optimal.
Because cells are larger than their nuclei, larger size of filters is needed. However, one major limitation with spatial domain image filters is non-continuous edges. Increasing the size of the filter increases the size of this “non-trustable” region. The use of a spatial filter 2r+1 in size will lead to the loss of r+1 pixels from the edge because of incomplete information, which greatly reduces the usable image size. For the nucleus, there are only ~12 pixels lost, which is acceptable. For an object the size of the cell, the much larger number of lost pixels is unacceptable. Thus, we did not use any spatial filter to segment individual cells. Rather, images of cells were processed with a 3×3 averaging smoothing filter.
To properly threshold single cell boundaries, measurement and elimination of the background intensity of images is more critical. First, we measured the average, <IBG>, and associated standard deviation, IRBG, of the background intensity of the smoothed imaged to obtain a set of pixel intensities less than <IBG> + 3.5*IRBG = Inn. Then, we update the value of <IBG> with Inn, IRBG, and the standard deviation of Inn. Three to five iterations will generally result in stable values of IBG and IRBG, which represent the average background intensity value and associated noise in background intensity magnitude, respectively. Next, we use IBG and IRBG to select the signal region of the fluorescence-labeled cells. We define the threshold factor, thc, and select all the pixels in the image with an intensity value larger than IBG + thc + IRBG. We make the assumption that the background noise intensity can be described by a Gaussian distribution and then set thc > 2, which represents > 95% of the background signal will not be selected. Depending on the signal intensity level, the value of thc will range between 2 and 5.
Using the above approach, we determined the cell boundary using phalloidin-stained F-actin images. F-actin usually gives a stronger signal at the cell boundary than at the cell center, differentiating the boundary from the cytoplasm with less bias than a more homogenous dye (such as HCS cell mask) would allow for. In fact, HCS cell mask intensities concentrated around the nucleus - the ticker region of the cell - and decayed towards the edge of the cell; because of the low NA objective, the edge intensity values was blurred, making edge detection very sensitive to bias and sample-to-sample variation.
Fig. S1. Cell cycle distributions of MD-MB-231 cells were verified by flow cytometry (FC) analysis. A. Cell cycle distribution of human breast cancer cells (MDA-MB-231) were measured by both FC and the microscopy assay. The number of analyzed cells is indicated in each panel. B. Fractions of cells in each phase of the cell cycle quantified with our microscopy-based assay and compared to FC.
Fig. S2. Cell cycle analysis applied to a wide range of cell types. To demonstrate the versatility of our method, seven additional types of cells were analyzed. These cells include human pancreatic normal epithelial cells (HPNE), human breast normal epithelial cells (MCF10A), patient-derived pancreatic cancer ductal adenocarcinoma cells that metastasized to the liver, LMNA+/+ mouse embryonic fibroblasts (MEFs), LMNA−/− MEFs, metastatic breast cancer cells (MD-MB-231), and patient-derived pancreatic cancer ductal adenocarcinoma cells from the primary tumor. At least 1,000 cells were analyzed for each cell line.
Fig. S3. Alteration of cell cycle causes long term effects on cell properties. A and B. Asynchronized cell cycle distribution is recovered after 24h release of from synchronizing drug. Orange represents thymidine treatment (A); purple represents nocodazole treatment (B). C and D. Mean changes in nuclear size and cell size induced 0h, 6h, and 24h after release of thymidine (C, orange) and nocodazole (D, purple) compared to control cells (blue). ***: P< 0.0001 (one-way ANOVA). E. Cell-cycle distributions of control cells and cells treated with Cdk4/6 inhibitor IV. Blue represents control. Black represents Cdk4/6 inhibitor IV treatment. F. Mean changes in nuclear size and cell size induced 0h, 6h, and 24h after removing Cdk4/6 inhibitor IV. Three biological repeats on different cells were analyzed for a total of >3,000 cells for each tested condition.
Fig. S4. Lamin deficiency enriches cells in the G0/G1 phase. A. Cell cycle distributions of LMNA−/− MEFs and Lamin A/C KD c2c12 cells compared to LMNA+/+ MEFs and control c2c12 cells transfected with a scrambled construct, respectively, obtained by the microscopy assay. B. Quantitative analysis demonstrates significant (P<0.0001) increase in the G0/G1 fractions and decrease in the G2/M fractions of both LMNA−/− MEFs and Lamin A/C KD c2c12 cells compared to control cells. C. Relative contributions to population-averaged changes in nuclear size of LMNA−/− MEFs and Lamin A/C KD c2c12 cells show only ~75% contribution from direct intrinsic nuclear size changes. The rest of ~25% is contributed from re-distribution of cell cycle and the second order term. At least 1,000 cells were analyzed in three biological repeats (different cells are analyzed each time) for a total of >3,000 cells for each cell line.
The authors thank Dr. Jeffrey Stirman and Prof. Gregory D. Longmore for their careful reading of the manuscript.