Bacteria adopt a wide variety of sizes and shapes, with many species exhibiting stereotypical morphologies. How morphology changes, and over what timescales, is less clear. Previous work examining cell morphology in an experiment with Escherichia coli showed that populations evolved larger cells and, in some cases, cells that were less rod-like. That experiment has now run for over two more decades. Meanwhile, genome sequence data are available for these populations, and new computational methods enable high-throughput microscopic analyses. Here, we measured stationary-phase cell volumes for the ancestor and 12 populations at 2,000, 10,000, and 50,000 generations, including measurements during exponential growth at the last timepoint. We measured the distribution of cell volumes for each sample using a Coulter counter and microscopy, the latter of which also provided data on cell shape. Our data confirm the trend toward larger cells, while also revealing substantial variation in size and shape across replicate populations. Most populations first evolved wider cells, but later reverted to the ancestral length-to-width ratio. All but one population evolved mutations in rod-shape maintenance genes. We also observed many ghost-like cells in the only population that evolved the novel ability to grow on citrate, supporting the hypothesis that this lineage struggles with maintaining balanced growth. Lastly, we show that cell size and fitness remain correlated across 50,000 generations. Our results suggest larger cells are beneficial in the experimental environment, while the reversion toward ancestral length-to-width ratios suggests partial compensation for the less favorable surface area-to-volume ratios of the evolved cells. IMPORTANCE Bacteria exhibit great morphological diversity, yet we have only a limited understanding of how their cell sizes and shapes evolve, and of how these features affect organismal fitness. This knowledge gap reflects, in part, the paucity of the fossil record for bacteria. Here, we revive and analyze samples extending over 50,000 generations from 12 populations of experimentally evolving Escherichia coli to investigate the relation between cell size, shape, and fitness. Using this “frozen fossil record” we show that all 12 populations evolved larger cells concomitant with increased fitness, with substantial heterogeneity in cell size and shape across the replicate lines. Our work demonstrates that cell morphology can readily evolve and diversify, even among populations living in identical environments.

Starving Myxococcus xanthus bacteria use short-range C-signaling to coordinate their movements and construct multicellular mounds, which mature into fruiting bodies as rods differentiate into spherical spores. Differentiation requires efficient C-signaling to drive the expression of developmental genes, but how the arrangement of cells within nascent fruiting bodies (NFBs) affects C-signaling is not fully understood. Here, we used confocal microscopy and cell segmentation to visualize and quantify the arrangement, morphology, and gene expression of cells near the bottom of NFBs at much higher resolution than previously achieved. We discovered that “transitioning cells” (TCs), intermediate in morphology between rods and spores, comprised 10 to 15% of the total population. Spores appeared midway between the center and the edge of NFBs early in their development and near the center as maturation progressed. The developmental pattern, as well as C-signal–dependent gene expression in TCs and spores, were correlated with cell density, the alignment of neighboring rods, and the tangential orientation of rods early in the development of NFBs. These dynamic radial patterns support a model in which the arrangement of cells within the NFBs affects C-signaling efficiency to regulate precisely the expression of developmental genes and cellular differentiation in space and time. Developmental patterns in other bacterial biofilms may likewise rely on short-range signaling to communicate multiple aspects of cellular arrangement, analogous to juxtacrine and paracrine signaling during animal development.

The cell morphology of rod-shaped bacteria is determined by the rigid net of peptidoglycan forming the cell wall. Alterations to the rod shape, such as the curved rod, occur through manipulating the process of cell wall synthesis. The human pathogen Vibrio cholerae typically exists as a curved rod, but straight rods have been observed under certain conditions. While this appears to be a regulated process, the regulatory pathways controlling cell shape transitions in V. cholerae and the benefits of switching between rod and curved shape have not been determined. We demonstrate that cell shape in V. cholerae is regulated by the bacterial second messenger cyclic dimeric guanosine monophosphate (c-di-GMP) by posttranscriptionally repressing expression of crvA, a gene encoding an intermediate filament-like protein necessary for curvature formation in V. cholerae. This regulation is mediated by the transcriptional cascade that also induces production of biofilm matrix components, indicating that cell shape is coregulated with V. cholerae’s induction of sessility. During microcolony formation, wild-type V. cholerae cells tended to exist as straight rods, while genetically engineering cells to maintain high curvature reduced microcolony formation and biofilm density. Conversely, straight V. cholerae mutants have reduced swimming speed when using flagellar motility in liquid. Our results demonstrate regulation of cell shape in bacteria is a mechanism to increase fitness in planktonic and biofilm lifestyles.

The cell morphology of rod-shaped bacteria is determined by the rigid net of peptidoglycan forming the cell wall. While V. cholerae grows into a curved shape under most conditions, straight rods have been observed. However, the signals and regulatory pathways controlling cell shape transitions in V. cholerae and the benefits of switching between rod and curved shape have not been determined. We demonstrate that cell shape in V. cholerae is regulated by the bacterial second messenger cyclic dimeric guanosine monophosphate (c-di-GMP) by repressing expression of crvA, a gene encoding an intermediate filament-like protein necessary for curvature formation in V. cholerae. This regulation is mediated by the transcriptional cascade that also induces production of biofilm matrix components, indicating that cell shape is coregulated with V. cholerae’s induction of sessility. Wild-type V. cholerae cells adhering to a surface lose their characteristic curved shape to become as straight as cells lacking crvA while genetically engineering cells to maintain high curvature reduced microcolony formation and biofilm density. Conversely, straight V. cholerae mutants have reduced speed when swimming using flagellar motility in liquid. Our results demonstrate regulation of cell shape in bacteria is a mechanism to increase fitness in planktonic or biofilm lifestyles.

Populations exposed to a high mutation rate harbor abundant deleterious genetic variation, leading to depressed mean fitness. This reduction in mean fitness presents an opportunity for selection to restore fitness through the evolution of mutational robustness. In extreme cases, selection for mutational robustness can lead to flat genotypes (with low fitness but high robustness) outcompeting fit genotypes (with high fitness but low robustness)—a phenomenon known as survival of the flattest. While this effect was previously explored using the digital evolution system Avida, a complete analysis of the local fitness landscapes of fit and flat genotypes has been lacking, leading to uncertainty about the genetic basis of the survival-of-the-flattest effect. Here, we repeated the survival-of-the-flattest study and analyzed the mutational neighborhoods of fit and flat genotypes. We found that the flat genotypes, compared to the fit genotypes, had a reduced likelihood of deleterious mutations as well as an increased likelihood of neutral and, surprisingly, of lethal mutations. This trend holds for mutants one to four substitutions away from the wild-type sequence. We also found that flat genotypes have, on average, no epistasis between mutations, while fit genotypes have, on average, positive epistasis. Our results demonstrate that the genetic causes of mutational robustness on complex fitness landscapes are multifaceted. While the traditional idea of the survival of the flattest emphasized the evolution of increased neutrality, others have argued for increased mutational sensitivity in response to strong mutational loads. Our results show that both increased neutrality and increased lethality can lead to the evolution of mutational robustness. Furthermore, strong negative epistasis is not required for mutational sensitivity to lead to mutational robustness. Overall, these results suggest that mutational robustness is achieved by minimizing heritable deleterious variation.