Rapid, precise quantification of bacterial cellular dimensions across a genomic-scale knockout library
- Tristan Ursell†1, 2,
- Timothy K. Lee†1,
- Daisuke Shiomi4, 5,
- Handuo Shi1,
- Carolina Tropini3, 9,
- Russell D. Monds1, 10,
- Alexandre Colavin3,
- Gabriel Billings6,
- Ilina Bhaya-Grossman1,
- Michael Broxton7,
- Bevan Emma Huang8,
- Hironori Niki4 and
- Kerwyn Casey Huang1, 9Email author
© Ursell et al. 2017
Received: 18 August 2016
Accepted: 6 January 2017
Published: 21 February 2017
The determination and regulation of cell morphology are critical components of cell-cycle control, fitness, and development in both single-cell and multicellular organisms. Understanding how environmental factors, chemical perturbations, and genetic differences affect cell morphology requires precise, unbiased, and validated measurements of cell-shape features.
Here we introduce two software packages, Morphometrics and BlurLab, that together enable automated, computationally efficient, unbiased identification of cells and morphological features. We applied these tools to bacterial cells because the small size of these cells and the subtlety of certain morphological changes have thus far obscured correlations between bacterial morphology and genotype. We used an online resource of images of the Keio knockout library of nonessential genes in the Gram-negative bacterium Escherichia coli to demonstrate that cell width, width variability, and length significantly correlate with each other and with drug treatments, nutrient changes, and environmental conditions. Further, we combined morphological classification of genetic variants with genetic meta-analysis to reveal novel connections among gene function, fitness, and cell morphology, thus suggesting potential functions for unknown genes and differences in modes of action of antibiotics.
Morphometrics and BlurLab set the stage for future quantitative studies of bacterial cell shape and intracellular localization. The previously unappreciated connections between morphological parameters measured with these software packages and the cellular environment point toward novel mechanistic connections among physiological perturbations, cell fitness, and growth.
Cell shape varies widely across bacterial species and has been linked to a diverse range of processes including adhesion, motility, pathogenicity, and differentiation . The cell wall, a polymer network of sugar strands crosslinked by short peptides, is both necessary and sufficient for defining a particular cell shape . The precise morphology and size of a cell is determined by the equilibrium between cell-wall mechanics and the turgor pressure caused by the high concentration of osmolytes inside the cell. Cell shape is maintained via feedback between the spatial pattern of cell-wall synthesis and the cell’s current geometry. In many rod-shaped organisms, especially those that grow by inserting cell wall material along the cylindrical region of the cell, the spatial pattern of growth is dictated by the cytoskeletal protein MreB , an actin homolog that forms filaments bound to the inner membrane . Cell volume is positively linked with fitness  and increases exponentially with growth rate on different nutrient sources that vary growth rate over a wide range . Moreover, in long-term evolution experiments with Escherichia coli, cell volume more than doubled in all evolved lines after 10,000 generations . Thus, an outstanding challenge in biology is to understand the mapping between genotype and morphological phenotypes.
To achieve such an understanding requires accurate quantification of cell morphology, particularly for very subtle changes such as the sub-micron curvature preference of MreB . Several computational tools were previously developed to quantify cell shape  in order to investigate intracellular organization and size homeostasis. The first software to interpolate cell contours at subpixel resolution was PSICIC , which has generally been applied to precisely quantify the subcellular localization of proteins. Simulations of point spread functions and their effects were combined with diffraction-limited imaging to achieve generational tracking and superior cell-division classification using ObjectJ [10, 11]. Another software package, Microbetracker, enabled segmentation of cells within a dense population ; Microbetracker and its successor Oufti  were recently used to investigate the relationships among growth rate, elongation, and division in E. coli  and Caulobacter crescentus [15, 16]. For rod-shaped bacteria, most quantitative studies involving cell size have essentially studied the dynamics of cell length, since cell width is generally maintained during elongation. However, E. coli B/r cells that experienced a nutrient upshift from minimal to rich medium increased in cell width progressively over a few doublings [17, 18], consistent with bulk measurements linking growth rate and cell volume . Further, mutations in MreB  and key cell-wall synthesis enzymes such as PBP2  have been identified that alter cell width, and sublethal doses of antibiotics such as A22, which depolymerizes MreB, or mecillinam, which inhibits PBP2, lead to cell-width increases in a concentration-dependent manner . Finally, osmotic shock subtly alters cell width , signifying a change in turgor pressure. These data are evidence that the cell’s ability to determine its width may be important for its regulation of cell growth and fitness. While powerful for many applications, packages such as PSICIC, Microbetracker, and MicrobeJ , the latter of which has an elegant interface for tracking lineages and measuring sub-cellular localization [22–24], require a relatively large number of parameters; measurements of cell width are sensitive to the values of these parameters. Critically, our ability to link these subtle shape changes to underlying genotypes and chemical environments relies on accurate, unbiased morphological characterization.
The Keio collection of single, nonessential gene deletions in E. coli BW25113 is a powerful resource for discovering the phenotypes of genes of unknown function . A visual screen of the qualitative shapes of the knockouts in this collection revealed only one mutant that was obviously non-rod-shaped . ∆rodZ cells are round, and it was subsequently found that RodZ interacts with MreB [26–28]. By profiling mutants from the Keio collection across hundreds of chemical treatments and environmental conditions, the functions of several genes have been discovered , such as the lipoprotein co-factors LpoA/B that activate the bifunctional penicillin binding proteins PBP1A/B, respectively . This chemical-genomics approach can be used to cluster genes whose functions are related by virtue of a common pathway. Given previous discoveries of close connections between cell size and growth rate  and size and fitness , measuring cell shape and size in distinct environments will likely reveal the mechanisms of growth regulation. Moreover, imaging data may constitute a phenotype vector for individual cells or populations of cells containing multiple morphological features such as cell width and length, curvature, and polar morphology . A preliminary analysis of cell shape classified mutants in the Keio collection as short, normal, long, or very long (https://shigen.nig.ac.jp/ecoli/strain/resource/keioCollection/list). However, detailed features such as cell width, size variability, or polar morphology have been difficult to accurately measure due to computational and software limitations.
To quantify various aspects of cell morphology, a software platform must accurately and robustly identify changes in cell width and curvature, ideally with high computational efficiency on imaging datasets from large libraries of strains. The focus of many existing software packages has been on defining a cell contour that can be used for comparing intracellular localization patterns or for computing the dynamics of a global parameter such as cell length. Datasets estimating local cell geometry with high accuracy can enable machine-learning tools to identify low-dimensional representations of cell shape and may reveal novel biological principles connecting cell shape to other behaviors. Principal Component Analysis (PCA) was previously harnessed to analyze the cell contours of populations of cells, leading to the identification of cell-shape modes in the bacterium C. crescentus  and in keratocytes . For C. crescentus, PCA enabled the clustering of MreB mutants , while for keratocytes, distinct PCA modes were strongly correlated with motility characteristics such as speed or turning . In previous studies, we developed software tools to analyze cell shape in a variety of contexts. Using phase-contrast images, we measured the changes in cell width and length resulting from point mutations in MreB ; we verified that the changes in cell width correlated with the distance between peaks in fluorescence of a membrane dye on opposite sides of the cell . In another study, we previously measured the curvature of E. coli cell contours and showed that MreB localized to concave regions of the cell , targeting new cell-wall growth to these locations and straightening the cell. When we measured the correlations among cell size, cytoskeletal dynamics, and cell twisting for cells with a range of sizes generated via genetic or chemical perturbations, we found that cells systematically altered cell-wall structure as cell width increased . In a morphological screen of the effects of depleting essential genes in B. subtilis, partial depletion led to cell-width outliers that highlighted both shape actuators (involved in cell-wall synthesis) and modulators (e.g., involved in DNA replication) . However, the molecular mechanisms that regulate cell size are currently relatively unknown, motivating an unbiased examination of a genomic-scale library.
Here we present cell-shape measurement tools in a user-friendly, computationally efficient MATLAB-based package called Morphometrics. This software segments cell contours from phase contrast images, fluorescence labeling of the cell surface, or cytoplasmic fluorescence, without assumptions about cell shape or size. Once cells have been identified, cell contours are determined via a straightforward, essentially parameter-free algorithm that yields robust measurements of cellular dimensions and contour curvature, allowing for automatic characterization of mutants with subtle variations in morphology. Since judging the accuracy of contour detection requires the ability to measure cell shape from images of cells with known size, we also present a separate software platform, BlurLab, that generates simulated fluorescence microscopy images . While other software tools have been developed to generate fluorescence images for predefined shapes and structures , BlurLab addresses arbitrary distributions of fluorescent molecules in space and time and has the capacity to mimic a wide variety of techniques, features, and sources of noise in light microscopy. In the current investigation, we used Morphometrics to measure cell size from phase contrast and fluorescence images, and BlurLab to validate our measurements and to determine the relative shifts in cellular dimensions between imaging modalities. We then applied Morphometrics to quantify cell shape and size across ~14,000 images of the Keio collection, revealing an inverse correlation between cell width and the robustness of cell-shape maintenance. Finally, we demonstrated that cell-shape parameters such as width and length correlate with particular chemical sensitivities. Ultimately, we envision that Morphometrics and Blurlab will provide fast, reproducible quantification of cell shape as well as the ability to test quantitative models, thus complementing canonical tools for biochemistry and cell biology.
Implementation and Results
Assessment of the consistency of cell contours determined by different imaging modalities
To facilitate quantitative analysis of contours from cells with a wide variety of shapes and sizes, we sought to implement an algorithm that made no assumptions about specific cell shapes and that extracted a parameter-free contour not subject to user biases. The Morphometrics algorithm can be conceptualized in two stages. First, with a small number of user-defined parameters, discrete ‘objects’ are detected as contiguous groups of pixels through watershed and distance-transform segmentation. Second, a smooth parameter-free contour defining the boundary of each object is calculated by treating the image intensity as a metric surface on which contour ‘energy’ can be minimized, with the segmented object boundary as an initial contour guess. Further information about the algorithm can be found in the Morphometrics user manual included with the software download. Three types of images can be used for contour detection: i) phase contrast images, in which the cell interior appears dark; ii) interior fluorescence images, e.g. from uniformly distributed cytoplasmic fluorescent proteins; and iii) peripheral fluorescence images, e.g. from membrane dye. Calculating the magnitude of the image gradient transforms the first two image types into an intensity map similar to that given by a fluorescence marker on the surface, which constitutes the common basis for calculating contours.
After optional image scaling, contrast adjustment, and background removal, the software offers multiple algorithms for segmenting contiguous groups of pixels (objects) that meet specified constraints on size and intensity. Each object is checked against criteria for false-positive detection based on the ratio of interior to boundary pixel intensities and, depending on user input, these objects may be linked with objects in other frames for tracking across a set of time-lapse images (see examples in the Morphometrics user manual). Parameters may be tested on individual images before being applied to the processing of large data sets. Ultimately, each segmented region serves as a seed to begin contour fitting. Once a contiguous object composed of discrete pixels is identified, a contour is calculated by treating the intensity features of the object as an interpolated energy landscape and then fitting a closed-loop contour in continuous coordinates to the minima of that energy landscape. From the cell contours, Morphometrics can be used to calculate an interior mesh; for rod-shaped cells, this mesh defines a cellular coordinate system with a midline and associated perpendicular meshlines that connect the two sides of the cell, thereby also measuring local cell width along the midline. The contour is also used to calculate one-dimensional profiles such as curvature (a measure of the radius of the circle that best fits the contour surrounding a particular point and whether the contour is concave or convex), fluorescence signals along the cell boundary (e.g., membrane dye, surface markers , or membrane-bound proteins ), or fluorescence signals along the interior centerline of the cell, among other features. Optical shifts between fluorescence imaging channels can be corrected by translating the contour coordinates by an amount appropriate to the particular imaging system; these translation values can be determined visually using the included post-processing contour viewing software.
From these cells and their associated meshline grids, we measured width profiles along the midline, obtaining values that increased from 0 at the two poles to ~1 μm near the middle of the cell (Fig. 1d). As previously reported , there was little intracellular variability in cell width along the midline away from the poles, and all three profiles showed similar variability (Fig. 1d). The contour measurements represented an approximately cylindrical body with hemispherical endcaps (Fig. 1d), although all three contour measurements displayed slight variations in curvature (Fig. 1e) that were previously shown to correlate with MreB localization [3, 40].
Simulated fluorescence microscopy for validation of quantitative analyses
Given the differences in cellular dimensions extracted from PC, IF, and PF data (Fig. 1), we wondered which modality accurately represented cell size. Other experimental methodologies with higher resolution such as electron microscopy still do not provide a “true” measure of cell size, since sample preparation likely perturbs the cell, for example by disrupting turgor pressure. To identify a strategy for comparing measurements with known geometric parameters, we developed BlurLab, a software package that generates simulated fluorescence images. BlurLab takes as input a set of locations of fluorescent molecules and convolves these locations with a point spread function (PSF) to generate a simulated image. The PSF can be directly measured for a particular microscope and objective using sub-diffraction-limited particles such as quantum dots or fluorescent beads; alternatively, BlurLab can generate a PSF for a given set of objective parameters (numerical aperture, wavelength, magnification, index of refraction, and pixel size). BlurLab can also mimic camera noise, thermal noise, and shot noise, yielding simulated images that are more realistic for head-to-head comparison with experimental images. Additional BlurLab functionalities include simulation of imaging at other focal planes for creating z-stacks, modeling total internal reflection fluorescence imaging, boxcar averaging of positions during simulated time-lapse imaging to account for particle motion during the exposure interval, simulating mean-field and stochastic photobleaching, and simulating fluorescence recovery after photobleaching. Detailed descriptions of these functionalities and examples can be found in the BlurLab manual.
Next, we examined the contours extracted from simulated images of in silico cells with different widths (Fig. 3c). The error in the extracted contour was practically zero for widths greater than 2 μm, but increased as cell width decreased (Fig. 3d). The extracted contours were smaller than the true midplane contours because the cell curvature introduced light from out-of-focus planes in which the cellular cross-section had a smaller width; for cell widths much larger than the wavelength of light used for imaging, this curvature became negligible (Fig. 3d). To determine the accuracy of other geometric measurements, we computed the curvature along the cell contour. The curvature along a 2-μm in silico cell was relatively accurate everywhere except in the transition region between the cylinder and the hemispherical end caps (Fig. 3e), where the PSF blurred the step-function curvature into a smooth transition. For a 1-μm in silico cell, the curvature remained accurate along the cylindrical portions of the cell (where the measured curvature was close to zero), although there was a slight overestimate of the curvature at the ends of the cell (Fig. 3e); this error was exacerbated as the cell width was further narrowed (Fig. 3e). We note that these systematic morphological errors with decreasing cell size result from the limitations of light imaging at wavelengths comparable to the cell size, not from imprecision in fluorescence simulation or contour detection.
To evaluate the overall bias, we used BlurLab to compute the difference between the actual and measured widths and areas for in silico cells 400 nm to 8 μm in width and a range of cell lengths in the peripheral fluorescence modality. The error in the width asymptotically approached zero, and was essentially undetectable in cells with width above ~1.5 μm (Fig. 3f). Error in area behaved in a similar fashion, regardless of cell length (Fig. 3g), indicating that PF measurements provide an accurate measure of cell length without systematic bias. This application illustrates the intended purpose and power of BlurLab: to assess the accuracy of imaging data in the context of a particular model by consistently comparing analyses of both experimental and simulated images (Fig. 3h).
Morphological analysis of a genomic library of nonessential gene deletions
To identify other morphological correlations across the library or in particular strains, we aligned 150 cell contours from each strain of the Keio collection for which at least 150 contours were available in order to calculate an average cell shape, and performed PCA on the covariance matrix to identify significant shape variations. The first two, three, and four PCA modes accounted for 97.3%, 99.0% and 99.3% of variation in cell shape, respectively (Fig. 4e–h), and by far the greatest amount of variation was accounted for by a mode that clearly captured elongation (Fig. 4e). This result, which is not surprising for rod-shaped growth, indicates that length changes are the most significant source of shape variation within the Keio library. Nonetheless, the next three modes representing cell bending, widening, and tapering, respectively (Fig. 4f–h), have potential for revealing cells or strains that are shape outliers. Mean width strongly correlated with PCA mode 3 (width, R = 0.998, Student’s t-test: p < 0.001) and with the projection from mode 1 (length, R = 0.43, Student’s t-test: p < 0.001) (Fig. 4e, bottom), as expected based on the correlation between length and width noted above (Fig. 4a). There was also a significant correlation between mean width and mode 4 (tapering, R = 0.49, Student’s t-test: p < 0.001) (Fig. 4h, bottom), potentially indicating a connection between cell-width determination and cell division. Our analysis demonstrates that most of the variation in cell shape is captured by length and width, although other morphological features such as tapering may be informative for characterizing certain outlier strains. PCA of the correlation matrix, which involves rescaling that avoids heavy skewing by cell length variation, led to an increased emphasis on cell width, bending, and tapering in the decomposition, with the mode corresponding to the largest eigenvalue representing tip morphology (Additional file 5: Figure S5). Taken together, correlations among features suggest underlying feedback between elongation and division in rod-shaped cells .
Identification of chemical sensitivities correlated with cellular dimensions
Given the range of cellular dimensions across the Keio library and the correlations between morphological observables, we wanted to systematically probe the physiological significance of cell size. We previously found that the MreBA53T mutation led to wider cells during growth in a variety of carbon sources; cells harboring this mutation had a large gain in fitness when competed against the parental strain in glucose-rich medium . However, the change in fitness was carbon-source dependent, with neutral fitness in lactose and a reduction in fitness in galactose . We also previously observed that sublethal treatment with the MreB inhibitor A22 led to a dose-dependent increase in cell width in wild-type E. coli MG1655 cells . Based on these data and the increase in width variability with increasing mean cell width across the Keio collection detected here (Fig. 4c), we hypothesized that wider cells may be more sensitive to A22 than thinner cells, and more generally that morphological observables may be predictive of the severity of phenotypes in certain environments or chemical treatments.
There were 15 unique conditions for which the knockouts had S-scores with significant positive correlation with cell width (Additional file 6: Table S2) , indicating that wider cells were less sensitive. These conditions target a broader range of cellular processes, including DNA/RNA, stress, and fatty acid metabolism in addition to the ribosome and cell-wall and membrane synthesis (Additional file 6: Table S2) . Carbonyl cyanide m-chlorophenyl hydrazone, a proton ionophore that inhibits oxidative phosphorylation, was positively correlated at all tested concentrations (Additional file 6: Table S2) ; interestingly, we observed in a separate study that this ionophore increases bending rigidity, which is an expected outcome of increased cell width based on the mechanics of a thin cylindrical shell. There was also a compound (theophylline) with an unknown target that exhibited a significant positive correlation with cell width (Additional file 6: Table S2) , indicating that cell-size profiling may be an effective tool for associating chemical exposure with the molecular pathways that control cell morphology.
We next considered correlations of S-scores with cell length. Six conditions were associated with significant negative correlations (Fig. 5cii): high iron, treatment with the cell-wall inhibitor ceftazidime, and four carbon-source limitations, which may reflect the known coupling of cell size with nutrient-dependent growth rate . Of the eight conditions positively correlated with length (Fig. 5cii), several involved ribosomal inhibitors (n = 3), inhibitors of fatty acid and membrane synthesis (n = 2), or DNA/RNA synthesis (n = 1). In some cases, correlations indicated differences in the underlying mode of drug action; for example, the cephalosporin cefaclor was positively correlated with cell length, while the structurally similar compound ceftazidime was negatively correlated with length.
Given that the observed correlations between cell size and S-scores involved chemicals with many target processes, we conjectured that the correlations for subsets of strains deleted for genes with similar functions may yield further insight into the relationships between cell size and cellular processes. We separated genes into 23 Clusters of Orthologous Groups, a common classification scheme . As compared to the correlations from the full library (Fig. 5c), we detected higher correlations between cell size and S-scores within some of these gene clusters (Fig. 5d,e; Additional file 6: Table S3) . A22 again was significantly (Student’s t-test) and negatively correlated with cell width for strains harboring deletions of genes associated with signal transduction mechanisms (Fig. 5ei), indicating that signaling pathways, such as those activated in response to stress, may respond to the changes in cell width caused by A22 treatment. Although the COG for cell-cycle control, cell division, and chromosome partitioning comprises only 34 genes, sensitivity to the division inhibitor verapamil nevertheless exhibited a strong and significant correlation of −0.92 with cell width (Fig. 5eii), while S-scores for minocycline sensitivity of knockouts of genes involved in nucleotide metabolism and transport were positivity correlated with cell width (R = 0.47, Fig. 5eiii). Surprisingly, for the large number of strains carrying deletions of genes of unknown function (COG class S, n = 918 genes), mean cell width was negatively correlated with sensitivity to low iron stress and the folic acid synthesis inhibitor sulfamethizole (Fig. 5d). Strains harboring deletions in genes involved in amino-acid transport and metabolism had mean cell lengths that were negatively correlated with sensitivity to four carbon-source starvations (Fig. 5d), while cells harboring deletions in genes involved in transcription and post-translational modification, protein turnover, and chaperones had lengths positively correlated with ribosomal inhibitor sensitivity (Fig. 5eiv). As with cell width, strains deleted for genes of unknown function had lengths both positively and negatively correlated with some treatments (Additional file 6: Table S4) , further supporting the importance of these genes to cell shape. Taken together, our correlation analyses based on the entire Keio collection and subdivided into Clusters of Orthologous Groups suggest that a variety of intracellular factors, beyond those controlling cell-wall synthesis or turgor, contribute to the cell’s determination of its size, and that responses to some extracellular perturbations have general connections with cell size (Fig. 5d). These correlations were revealed through precise and automated measurements of cellular dimensions that were enabled by our open-source software packages Morphometrics and BlurLab.
Discussion and Conclusions
Quantifying subtle connections among environmental conditions, cellular morphology, and genetics is a crucial step in uncovering new biological roles for genes and their corresponding phenotypes. Microscopy is a powerful workhorse for establishing these connections, but the spatial constraints of light microscopy and our computational assumptions limit the set of measurable morphological perturbations. We developed Morphometrics and BlurLab to expand the utility of multiple imaging modalities and to minimize the presence of bias in image-based analyses. Here, we applied our software to bacteria, which are a technically challenging class of organisms to explore with quantitative image processing due to their small size and strong connections among fitness, genotype, and morphology. We envision that combining the morphological measurement techniques of Morphometrics with the unbiased hypothesis testing of BlurLab will enable the quantitative and automated characterization of libraries targeting cell morphology (such as libraries constructed through error-prone PCR mutagenesis of genes that affect cell morphology  and the set of Keio mutants with a range of cell lengths and widths (Fig. 4)) as well as the growing collection of genomic libraries, including knockout libraries of Salmonella typhimurium  and transposon-based libraries of the pathogen Pseudomonas aeruginosa  and the gut commensal Bacteroides thetaiotaomicron . The relationships we discovered among cell length, cell width, and width variability (Fig. 4) across the genomic-scale Keio collection suggest general links between cellular physiology and cell size . Screens that used the Keio collection to reveal new phenotypes previously focused on growth  or envelope permeability , both of which naturally led to studies of cell shape. Since Morphometrics is sensitive to small morphological variations and can rapidly analyze tens of thousands of images, it crucially enables forward genetic screens for genes or mutations that affect a morphological phenotype.
The relatively high throughput and unbiased morphological characterization of Morphometrics make it well-suited to the study of morphological and functional connections in genomic libraries, to the screening of environmental conditions like carbon sources, osmolytes, and antibiotics, and to analyses of dynamic cell-cycle data via time-lapse imaging . Similarly, new rules of cellular homeostasis encoded by time-dependent morphology in many organisms [14–16, 49] will be accessible with our software, particularly when spatial resolution is limited. It remains to be seen to what extent detailed measurements of cell morphology will be sufficient to uncover relationships among genotypes, chemical treatments, and cell morphology (Fig. 5); in some cases, lack of shape variation or degeneracies in shape phenotypes may make identification of a morphological signature challenging. For instance, drug treatments and attendant mutations could lead to global changes in gene expression such as a stress response that result in similar non-specific changes in cell morphology (e.g. filamentation).
While we have focused on bacteria, Morphometrics is amenable to morphological analysis of any organism for which high-contrast images can be obtained, especially walled organisms like plants (Fig. 2a)  and fungi (Fig. 2e) . Morphometrics has relatively few parameters and makes no assumptions about object morphology, strengths that are important considerations when tuning algorithms to examine new organisms with poorly understood morphology and growth cycles and when developing custom scripts for post-processing, such as the detection of spatial and temporal morphological correlations , detailed localization studies , and testing of biophysical models . The small size of bacteria leads to a smaller visual dynamic range, which makes human observers more prone to apophenia (seeing meaningful patterns in random data). BlurLab has a wide array of potential applications in defining null hypotheses about fluorescence imaging data, such as the expected distribution of a uniform surface label at the poles versus the lateral walls of a rod-shaped cell, the integrated intensity of a homogeneous cytoplasmic label such as 4′,6-diamidino-2-phenylindole (DAPI) as a function of cell size (given that wider cells have a greater fraction of the cell volume at the edges of the focal plane), and the level of noise in cytoplasmic distributions as a function of molecule number for testing whether noise is spatially heterogeneous. BlurLab can be coupled to any analysis platform, and should have growing utility given recent genome-scale screens of fluorescence localization [48, 52]. We note that evaluations of null hypotheses should also be important for eukaryotic cells, for which it is equally critical to connect molecular models to diffraction-limited fluorescence images .
For large data sets, image processing must be efficient in terms of computational resources and speed. The data structures, algorithms, and graphical user interfaces of Morphometrics and Blurlab are designed to exploit optimized image processing algorithms in MATLAB and to allow users to easily construct seamless custom scripts for data analysis ; these packages also take advantage of the large body of user-generated algorithms online (e.g. the MATLAB File Exchange). Additionally, high-level programming languages reduce the barrier to acquire coding skills and ensure that the varied and general needs of the quantitative imaging community for custom analysis can be met. Speed limitations in MATLAB can be resolved by recoding computationally intensive subroutines in C, as we did in Morphometrics. While Morphometrics and BlurLab enable new levels of precision and quantitative morphological characterization, like any analysis software they have certain limitations. Ultimately, the precision of contours from Morphometrics is limited by i) image quality, specifically spatial resolution, signal-to-noise, evenness of illumination, dynamic range of intensities, and degree of saturation; and ii) sample characteristics such as object contrast, the proximity of objects and attendant overlap in the light fields, the degree of blurring due to motion, and, where fine segmentation is desired, the number and degree of construction points in the image outline. As a general rule of thumb, images that are visually difficult to segment will be difficult for Morphometrics to segment. BlurLab is a simulator for linear optical microscopy, but specific nuances of an optical system or objective, atypical or time-dependent noise sources, camera chip-specific noise, and sample-dependent effects (e.g. absorption, quenching, or a sample’s index of refraction) all reduce the accuracy of simulation. Similar limitations will very likely apply to any image analysis or simulation software.
The importance of quantitation in cell biology will continue to increase, and the small size of bacteria and lack of organelles make both Morphometrics and BlurLab especially important for identifying subtle localization and morphological phenotypes in these organisms. However, despite the breadth and versatility of Morphometrics and BlurLab, significant challenges for computational image processing remain, including reconstruction of three-dimensional morphology in both static and dynamic environments, segmentation and tracking of dense and/or highly dynamic groups of objects, and the development of algorithms to process different imaging modalities such as fluorescence recovery after photobleaching, total internal reflection fluorescence, and super-resolution imaging. Morphometrics and BlurLab should serve as a foundation for developing new software to address these challenges.
Precise quantification of cell morphology in bacteria and eukaryotes will undoubtedly be a valuable tool for mapping genotype-phenotype relationships. Extending the analysis of static images carried out here to the dynamic response of cells to perturbations, for instance in microfluidic chambers, can further reveal the physiological basis of a particular phenotype, such as the mechanism of cell death during entry into stationary phase in an E. coli mutant with disrupted lipid homeostasis . Excitingly, cell morphology can even serve as a diagnostic tool for assessing cellular states in diseases such as cancer. To fully exploit the information obtained through these studies, we must develop and implement computational tools with high levels of accuracy and precision and couple them to methods for validation and compelling visual display. As demonstrated here, Morphometrics and BlurLab constitute an important step toward meeting these goals.
Wild-type MG1655 E. coli cells expressing cytoplasmic GFP and a kanamycin resistance cassette from plasmid pZS21-GFP (gift from Tom Silhavy, Princeton University) were labeled with the N-acetylglucosamine- and sialic acid-specific lectin Wheat Germ Agglutinin conjugated to Alexa-594 (Life Technologies). Five milliliters of cells were grown in lysogeny broth (LB) with shaking at 37 °C to exponential phase (optical density at 620 nm ∼ 0.5). One milliliter of cells was washed with fresh LB via centrifugation (10,000 g) and resuspension in 1 mL of LB and subsequently diluted 1:10 into 1 mL of fresh LB. Twenty-five microliters of a once-frozen 1 mg/mL fluorescent Wheat Germ Agglutinin stock solution were added, and the sample was briefly vortexed. Cells were incubated with the lectin for 20 min (approximately one cell cycle) with shaking at 37 °C in the dark. After incubation, cells were washed twice with fresh LB to remove excess lectin, and 5 μL of labeled cells were deposited onto a LB + 1% agarose pad, allowed to air dry on the pad, and promptly sealed with a #1.5 coverslip in a 125-μL FastWell (Grace BioLabs).
Labeled cells were imaged on a Nikon Eclipse Ti-E inverted fluorescence microscope with a 100X (NA 1.40) oil-immersion objective (Nikon Instruments). Images were collected using an Andor DU885 EMCCD camera (Andor Technology). Cells were maintained at 37 °C during imaging with an active-control environmental chamber (HaisonTech). Images were collected using μManager v. 1.3 .
Imaging of the Keio collection
Images were obtained from the NBRP. In brief, to obtain these images, strains from the Keio collection were inoculated in LB with 30 μg/mL kanamycin and grown overnight in 96-well plates at 30 °C. Cells were then diluted in LB plus 30 μg/mL kanamycin and grown for 2 h at 37 °C. After reaching exponential phase, cells were harvested via centrifugation and resuspended in LB. These cells were mounted on poly-lysine-coated cover-slips, fixed with methanol, washed with water, and stained with (4′,6-diamidino-2-phenylindole).
Analysis of Keio images
Keio collection images from the NBRP (1–3 per strain) were analyzed using Morphometrics. Cells with segmentation errors were filtered by only including contours with two identifiable points of high curvature (corresponding to the poles). This filtering eliminated segmentation errors, which we evaluated with manual curation (data not shown). Cell length and width were calculated according to a mesh representation of the cell contour computed by Morphometrics. Subsequent analyses only included strains for which 150 cells passed the above filtering step (2465/4353 strains).
Shape-variation modes were calculated from cell contours using PCA. Briefly, the center of mass and principal axes were calculated from the cell contour and the coordinates were shifted and rotated to a common alignment. Then, 300 equally spaced points were sampled from the cell contour using linear interpolation. A mean cell contour was subtracted from each contour and principal components were calculated using the eigenvalue decomposition of the covariance matrix between contour coordinates.
For comparing cell morphology to chemical sensitivity (Fig. 5), Pearson’s correlation coefficient between the mean cell width or length against the S-score across all strains was calculated for each chemical condition. The statistical significance of each correlation was calculated with a Student’s t-distribution, Bonferroni-corrected with the number of conditions.
The authors thank Nassos Typas and the Huang lab for helpful discussions. We also thank members of the Shaevitz, Gitai, and Garner laboratories for software alpha testing and feedback. Images used in Fig. 2 were kindly provided by David Quint (Fig. 2a), Leigh Harris (Fig. 2b,c), Wikipedia under the Creative Commons License (Fig. 2d), Liam Holt (Fig. 2e), Cynthia Whitchurch (Fig. 2f), Joseph Dillard (Fig. 2g), Anthony Shiver (Fig. 2h), Tristan Ursell (Fig. 2i), and Handuo Shi (Fig. 2j).
This work was supported by a Bio-X Postdoctoral Fellowship and University of Oregon seed funds (to T.S.U.), a Stanford Interdisciplinary Graduate Fellowship (to C.T.), a Siebel Scholars Graduate Fellowship (to T.K.L.), support from a National Institutes of Health (NIH) Biotechnology Training Grant (to T.K.L.), an Agilent Fellowship and a Stanford Interdisciplinary Graduate Fellowship (to H.S.), a Bio-X Senior Postdoctoral Fellowship (to R.D.M.), a Stanford Graduate Fellowship and a Gerald J. Lieberman Fellowship (to A.C.), a National Science Foundation Graduate Research Fellowship (to G.B.), National Institutes of Health Director’s New Innovator Award DP2OD006466 (to K.C.H.), National Science Foundation CAREER award MCB-1149328 (to K.C.H.), and a Bio-X Interdisciplinary Initiatives Program seed grant (to K.C.H.). 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
Morphometrics, an associated manual, and the images used in Fig. 2 can be downloaded at https://simtk.org/projects/morphometrics. BlurLab and an associated manual can be downloaded at https://simtk.org/projects/blurlab. All data generated or analysed during this study are included in this published article and its supplementary information files, or are available upon request.
TU, TKL, and KCH conceived the project. TU was the major developer of both Morphometrics and BlurLab, and carried out experimental measurements and simulations to measure cell size through different modalities. Keio collection images were generated by DS and HN, and TKL analyzed the images and performed PCA. HS, CT, RDM, AC, and GB contributed to Morphometrics development. IBG contributed to Morphometrics validation. CT contributed to BlurLab validation. MB contributed the PSF generation algorithm in BlurLab. BEH provided input on statistics. The paper was written by TU, TKL, and KCH. 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.
- Young KD. The selective value of bacterial shape. Microbiol Mol Biol Rev. 2006;70(3):660–703.View ArticlePubMedPubMed CentralGoogle Scholar
- Holtje JV. Growth of the stress-bearing and shape-maintaining murein sacculus of Escherichia coli. Microbiol Mol Biol Rev. 1998;62(1):181–203.PubMedPubMed CentralGoogle Scholar
- Ursell TS, Nguyen J, Monds RD, Colavin A, Billings G, Ouzounov N, et al. Rod-like bacterial shape is maintained by feedback between cell curvature and cytoskeletal localization. Proc Natl Acad Sci U S A. 2014;111(11):E1025–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Salje J, van den Ent F, de Boer P, Lowe J. Direct Membrane Binding by Bacterial Actin MreB. Mol Cell. 2011;43(3):478–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Monds RD, Lee TK, Colavin A, Ursell T, Quan S, Cooper TF, et al. Systematic perturbation of cytoskeletal function reveals a linear scaling relationship between cell geometry and fitness. Cell Rep. 2014;9(4):1528–37.View ArticlePubMedGoogle Scholar
- Schaechter M, Maaloe O, Kjeldgaard NO. Dependency on medium and temperature of cell size and chemical composition during balanced grown of Salmonella typhimurium. J Gen Microbiol. 1958;19(3):592–606.View ArticlePubMedGoogle Scholar
- Lenski RE, Travisano M. Dynamics of adaptation and diversification: a 10,000-generation experiment with bacterial populations. Proc Natl Acad Sci. 1994;91(15):6808–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Young JW, Locke JC, Altinok A, Rosenfeld N, Bacarian T, Swain PS, et al. Measuring single-cell gene expression dynamics in bacteria using fluorescence time-lapse microscopy. Nat Protoc. 2012;7(1):80–8.View ArticleGoogle Scholar
- Guberman JM, Fay A, Dworkin J, Wingreen NS, Gitai Z. PSICIC: noise and asymmetry in bacterial division revealed by computational image analysis at sub-pixel resolution. PLoS Comput Biol. 2008;4(11):e1000233.View ArticlePubMedPubMed CentralGoogle Scholar
- Syvertsson S, Vischer NO, Gao Y, Hamoen LW. When Phase Contrast Fails: ChainTracer and NucTracer, Two ImageJ Methods for Semi-Automated Single Cell Analysis Using Membrane or DNA Staining. PLoS One. 2016;11(3):e0151267.View ArticlePubMedPubMed CentralGoogle Scholar
- Vischer NOE, Verheul J, Postma M, van den Berg van Saparoea B, Galli E, Natale P, Gerdes K, Luirink J, Vollmer W, Vicente M, den Blaauwen T. Cell age dependent concentration of Escherichia coli divisome proteins analyzed with ImageJ and ObjectJ. Front Microbiol. 2015;6:586.Google Scholar
- Sliusarenko O, Cabeen MT, Wolgemuth CW, Jacobs-Wagner C, Emonet T. Processivity of peptidoglycan synthesis provides a built-in mechanism for the robustness of straight-rod cell morphology. Proc Natl Acad Sci U S A. 2010;107(22):10086–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Paintdakhi A, Parry B, Campos M, Irnov I, Elf J, Surovtsev I, et al. Oufti: an integrated software package for high-accuracy, high-throughput quantitative microscopy analysis. Mol Microbiol. 2016;99(4):767–77.View ArticlePubMedGoogle 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
- 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
- Iyer-Biswas S, Wright CS, Henry JT, Lo K, Burov S, Lin Y, et al. Scaling laws governing stochastic growth and division of single bacterial cells. Proc Natl Acad Sci U S A. 2014;111(45):15912–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Grover NB, Zaritsky A, Woldringh CL, Rosenberger RF. Dimensional rearrangement of rod-shaped bacteria following nutritional shift-up. I. Theory. J Theor Biol. 1980;86(3):421–39.View ArticlePubMedGoogle Scholar
- Woldringh CL, Grover NB, Rosenberger RF, Zaritsky A. Dimensional rearrangement of rod-shaped bacteria following nutritional shift-up. II. Experiments with Escherichia coli B/r. J Theor Biol. 1980;86(3):441–54.View ArticlePubMedGoogle Scholar
- Philippe N, Alcaraz JP, Coursange E, Geiselmann J, Schneider D. Improvement of pCVD442, a suicide plasmid for gene allele exchange in bacteria. Plasmid. 2004;51(3):246–55.View ArticlePubMedGoogle Scholar
- Tropini C, Lee TK, Hsin J, Desmarais SM, Ursell T, Monds RD, et al. Principles of bacterial cell-size determination revealed by cell-wall synthesis perturbations. Cell Rep. 2014;9(4):1520–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Rojas E, Theriot JA, Huang KC. Response of Escherichia coli growth rate to osmotic shock. Proc Natl Acad Sci U S A. 2014;111(21):7807–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Ducret A, Quardokus EM, Brun YV. MicrobeJ, a tool for high throughput bacterial cell detection and quantitative analysis. Nat Microbiol. 2016;1:16077.View ArticlePubMedPubMed CentralGoogle Scholar
- Alberge F, Espinosa L, Seduk F, Sylvi L, Toci R, Walburger A, et al. Dynamic subcellular localization of a respiratory complex controls bacterial respiration. eLife. 2015;4:e05357.View ArticlePubMed CentralGoogle Scholar
- Jiang C, Brown PJ, Ducret A, Brun YV. Sequential evolution of bacterial morphology by co-option of a developmental regulator. Nature. 2014;506(7489):489–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, et al. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006;2:2006 0008.View ArticlePubMedPubMed CentralGoogle Scholar
- Shiomi D, Sakai M, Niki H. Determination of bacterial rod shape by a novel cytoskeletal membrane protein. EMBO J. 2008;27(23):3081–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Alyahya SA, Alexander R, Costa T, Henriques AO, Emonet T, Jacobs-Wagner C. RodZ, a component of the bacterial core morphogenic apparatus. Proc Natl Acad Sci U S A. 2009;106(4):1239–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Bendezu FO, Hale CA, Bernhardt TG, de Boer PA. RodZ (YfgA) is required for proper assembly of the MreB actin cytoskeleton and cell shape in E. coli. EMBO J. 2009;28(3):193–204.View ArticlePubMedGoogle Scholar
- Nichols RJ, Sen S, Choo YJ, Beltrao P, Zietek M, Chaba R, et al. Phenotypic landscape of a bacterial cell. Cell. 2011;144(1):143–56.View ArticlePubMedGoogle Scholar
- Typas A, Banzhaf M, van den Berg van Saparoea B, Verheul J, Biboy J, Nichols RJ, et al. Regulation of peptidoglycan synthesis by outer-membrane proteins. Cell. 2010;143(7):1097–109.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang KC. Applications of imaging for bacterial systems biology. Curr Opin Microbiol. 2015;27:114–20.View ArticlePubMedGoogle Scholar
- Pincus Z, Theriot JA. Comparison of quantitative methods for cell-shape analysis. J Microsc. 2007;227(Pt 2):140–56.View ArticlePubMedGoogle Scholar
- Keren K, Pincus Z, Allen GM, Barnhart EL, Marriott G, Mogilner A, et al. Mechanism of shape determination in motile cells. Nature. 2008;453(7194):475–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Dye NA, Pincus Z, Fisher IC, Shapiro L, Theriot JA. Mutations in the nucleotide binding pocket of MreB can alter cell curvature and polar morphology in Caulobacter. Mol Microbiol. 2011;81(2):368–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Peters JM, Colavin A, Shi H, Czarny TL, Larson MH, Wong S, et al. A Comprehensive, CRISPR-based Functional Analysis of Essential Genes in Bacteria. Cell. 2016;165(6):1493–506.View ArticlePubMedGoogle Scholar
- Ursell T, Huang K. Resolution limits of optical microscopy and the mind. Biomed Comput Rev. 2011;7:27–9.Google Scholar
- Gardner MK, Odde DJ, Bloom K. Hypothesis testing via integrated computer modeling and digital fluorescence microscopy. Methods. 2007;41(2):232–7.View ArticlePubMedGoogle Scholar
- Ursell TS, Trepagnier EH, Huang KC, Theriot JA. Analysis of surface protein expression reveals the growth pattern of the Gram-negative outer membrane. PLoS Comput Biol. 2012;8(9):e1002680.View ArticlePubMedPubMed CentralGoogle Scholar
- Furchtgott L, Wingreen NS, Huang KC. Mechanisms for maintaining cell shape in rod-shaped Gram-negative bacteria. Mol Microbiol. 2011;81(2):340–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Billings G, Ouzounov N, Ursell T, Desmarais SM, Shaevitz J, Gitai Z, et al. De novo morphogenesis in L-forms via geometric control of cell growth. Mol Microbiol. 2014;93(5):883–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Auer GK, Lee TK, Rajendram M, Cesar S, Miguel A, Huang KC, et al. Mechanical genomics: high-throughput identification of bacterial cell-stiffness modulators. Cell Syst. 2016;2:402–11.Google Scholar
- Harris LK, Theriot JA. Relative Rates of Surface and Volume Synthesis Set Bacterial Cell Size. Cell. 2016;165(6):1479–92.View ArticlePubMedGoogle Scholar
- Galperin MY, Makarova KS, Wolf YI, Koonin EV. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 2015;43(Database issue):D261–9.View ArticlePubMedGoogle Scholar
- Porwollik S, Santiviago CA, Cheng P, Long F, Desai P, Fredlund J, et al. Defined single-gene and multi-gene deletion mutant collections in Salmonella enterica sv Typhimurium. PLoS One. 2014;9(7):e99820.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee SA, Gallagher LA, Thongdee M, Staudinger BJ, Lippman S, Singh PK, et al. General and condition-specific essential functions of Pseudomonas aeruginosa. Proc Natl Acad Sci U S A. 2015;112(16):5189–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Cullen TW, Schofield WB, Barry NA, Putnam EE, Rundell EA, Trent MS, et al. Gut microbiota. Antimicrobial peptide resistance mediates resilience of prominent gut commensals during inflammation. Science. 2015;347(6218):170–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Paradis-Bleau C, Kritikos G, Orlova K, Typas A, Bernhardt TG. A genome-wide screen for bacterial envelope biogenesis mutants identifies a novel factor involved in cell wall precursor metabolism. PLoS Genet. 2014;10(1):e1004056.View ArticlePubMedPubMed CentralGoogle Scholar
- Kuwada NJ, Traxler B, Wiggins PA. Genome-scale quantitative characterization of bacterial protein localization dynamics throughout the cell cycle. Mol Microbiol. 2015;95(1):64–79.View ArticlePubMedGoogle Scholar
- Brown PJ, de Pedro MA, Kysela DT, Van der Henst C, Kim J, De Bolle X, et al. Polar growth in the Alphaproteobacterial order Rhizobiales. Proc Natl Acad Sci U S A. 2012;109(5):1697–701.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang KC, Ehrhardt DW, Shaevitz JW. The molecular origins of chiral growth in walled cells. Curr Opin Microbiol. 2012;15(6):707–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou Z, Munteanu EL, He J, Ursell T, Bathe M, Huang KC, et al. The contractile ring coordinates curvature-dependent septum assembly during fission yeast cytokinesis. Mol Biol Cell. 2015;26(1):78–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Werner JN, Chen EY, Guberman JM, Zippilli AR, Irgon JJ, Gitai Z. Quantitative genome-scale analysis of protein localization in an asymmetric bacterium. Proc Natl Acad Sci U S A. 2009;106(19):7858–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Rafelski SM, Viana MP, Zhang Y, Chan YH, Thorn KS, Yam P, et al. Mitochondrial network size scaling in budding yeast. Science. 2012;338(6108):822–4.View ArticlePubMedPubMed CentralGoogle 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: E1565–74.Google Scholar
- Edelstein A, Amodaj N, Hoover K, Vale R, Stuurman N. Computer control of microscopes using µManager. Current protocols in molecular biology/edited by Frederick M Ausubel [et al.]. 2010; Chapter 14: Unit14 20.Google Scholar