Search

Many species of pathogenic bacteria harbor critical plasmid-encoded virulence factors, and yet the regulation of plasmid replication is often poorly understood despite playing a critical role in plasmid-encoded gene expression. Human pathogenic Yersinia , including the plague agent Y. pestis and its close relative Y. pseudotuberculosis , require the type III secretion system (T3SS) virulence factor to subvert host defense mechanisms and colonize host tissues. The Yersinia T3SS is encoded on the IncFII plasmid for Yersinia virulence (pYV). Several layers of gene regulation enables a large increase in expression of Yersinia T3SS genes at mammalian body temperature. Surprisingly, T3SS expression is also controlled at the level of gene dosage. The number of pYV molecules relative to the number of chromosomes per cell, referred to as plasmid copy number, increases with temperature. The ability to increase and maintain elevated pYV plasmid copy number, and therefore T3SS gene dosage, at 37°C is important for Yersinia virulence. In addition, pYV is highly stable in Yersinia at all temperatures, despite being dispensable for growth outside the host. Yet how Yersinia reinforces elevated plasmid replication and plasmid stability remains unclear. In this study, we show that the chromosomal gene pcnB encoding the polyadenylase PAP I is required for regulation of pYV plasmid copy number (PCN), maintenance of pYV in the bacterial population outside the host, robust T3SS activity, and Yersinia virulence in a mouse infection model. Likewise, pcnB /PAP I is also required for robust expression of the Shigella flexneri virulence plasmid-encoded T3SS. Furthermore, Yersinia and Shigella pcnB /PAP I is required for maintaining normal PCN of model antimicrobial resistance (AMR) plasmids whose replication is regulated by sRNA, thereby increasing antibiotic resistance by ten-fold. These data suggest that pcnB /PAP I contributes to the spread and stabilization of virulence and AMR plasmids in bacterial pathogens, and is essential in maintaining the gene dosage required to mediate plasmid-encoded traits. Importantly pcnB /PAP I has been bioinformatically identified in many species of bacteria despite being studied in only a few species to date. Our work highlights the potential importance of pcnB /PAP I in antibiotic resistance, and shows for the first time that pcnB /PAP I reinforces PCN and virulence plasmid stability in natural pathogenic hosts with a direct impact on bacterial virulence.

In proliferating bacteria, growth rate is often assumed to be similar between daughter cells. However, most of our knowledge of cell growth derives from studies on symmetrically dividing bacteria. In many α-proteobacteria, asymmetric division is a normal part of the life cycle, with each division producing daughter cells with different sizes and fates. Here, we demonstrate that the functionally distinct swarmer and stalked daughter cells produced by the model α-proteobacterium Caulobacter crescentus can have different average growth rates under nutrient-replete conditions despite sharing an identical genome and environment. The discrepancy in growth rate is due to a growth slowdown associated with the cell cycle stage preceding DNA replication (the G1 phase), which initiates in the late predivisional mother cell before daughter cell separation. Both progenies experience a G1-associated growth slowdown, but the effect is more severe in swarmer cells because they have a longer G1 phase. Activity of SpoT, which produces the (p)ppGpp alarmone and extends the G1 phase, accentuates the cell cycle-dependent growth slowdown. Collectively, our data identify a coupling between cell growth, the G1 phase, and asymmetric division that C. crescentus may exploit for environmental adaptation through SpoT activity. This coupling differentially modulates the growth rate of functionally distinct daughter cells, thereby altering the relative abundance of ecologically important G1-specific traits within the population., Competing Interests: Competing interests statement:The authors declare no competing interest.

The interplay between the human innate immune system and bacterial cell wall components is pivotal in understanding diseases such as Crohn's disease and Lyme arthritis. Lyme disease, caused by Borrelia burgdorferi , is the most prevalent tick-borne illness in the United States, with a substantial number of cases reported annually. While antibiotic treatments are generally effective, approximately 10% of Lyme disease cases develop persistent arthritis, suggesting a dysregulated host immune response. We have previously identified a link between the immunogenic B. burgdorferi peptidoglycan (PG) and Lyme arthritis and showed that this pathogen sheds significant amounts of PG fragments during growth. Here, we synthesize these PG fragments, including ornithine-containing monosaccharides and disaccharides, to mimic the unique composition of Borrelia cell walls, using reproducible and rigorous synthetic methods. This synthetic approach allows for the modular preparation of PG derivatives, providing a diverse library of well-defined fragments. These fragments will serve as valuable tools for investigating the role of PG-mediated innate immune response in Lyme disease and aid in the development of improved diagnostic methods and treatment strategies.

Conventional antimicrobial discovery relies on targeting essential enzymes in pathogenic organisms, contributing to a paucity of new antibiotics to address resistant strains. Here, by targeting a non-essential enzyme, Borrelia burgdorferi HtpG, to deliver lethal payloads, we expand what can be considered druggable within any pathogen. We synthesized HS-291, an HtpG inhibitor tethered to the photoactive toxin verteporfin. Reactive oxygen species, generated by light, enables HS-291 to sterilize Borrelia cultures by causing oxidation of HtpG, and a discrete subset of proteins in proximity to the chaperone. This caused irreversible nucleoid collapse and membrane blebbing. Tethering verteporfin to the HtpG inhibitor was essential, since free verteporfin was not retained by Borrelia in contrast to HS-291. For this reason, we liken HS-291 to a berserker, wreaking havoc upon the pathogen's biology once selectively absorbed and activated. This strategy expands the druggable pathogenic genome and offsets antibiotic resistance by targeting non-essential proteins., Competing Interests: Declaration of interests T.A.J.H. and P.F.H. have multiple patents issued or disclosed with Duke University around the tethering technology described in this study., (Copyright © 2023 The Authors. Published by Elsevier Ltd.. All rights reserved.)

This personal story recounts the accidental observation, the struggles, the breakthroughs, and the collaborative spirit of a few individuals that led to the discovery that bacterial cells expend energy to effectively fluidize their otherwise "glass-like" cytoplasm and promote the dispersal of large cytoplasmic components. This adventure, which led us into an uncharted world at the intersection of cell biology and condensed matter physics about ten years ago, forever transformed the way I view cells and conduct research., Competing Interests: Declaration of interests The author declares no competing interests., (Copyright © 2023 Elsevier Inc. All rights reserved.)

To examine how bacteria achieve robust cell proliferation across diverse conditions, we developed a method that quantifies 77 cell morphological, cell cycle, and growth phenotypes of a fluorescently labeled Escherichia coli strain and >800 gene deletion derivatives under multiple nutrient conditions. This approach revealed extensive phenotypic plasticity and deviating mutant phenotypes were often nutrient dependent. From this broad phenotypic landscape emerged simple and robust unifying rules (laws) that connect DNA replication initiation, nucleoid segregation, FtsZ ring formation, and cell constriction to specific aspects of cell size (volume, length, or added length) at the population level. Furthermore, completion of cell division followed the initiation of cell constriction after a constant time delay across strains and nutrient conditions, identifying cell constriction as a key control point for cell size determination. Our work provides a population-level description of the governing principles by which E. coli integrates cell cycle processes and growth rate with cell size to achieve its robust proliferative capability. A record of this paper's transparent peer review process is included in the supplemental information., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

Borrelia burgdorferi, a causative agent of Lyme disease, contains the most segmented bacterial genome known to date, with one linear chromosome and over twenty plasmids. How this unusually complex genome is organized, and whether and how the different replicons interact are unclear. We recently demonstrated that B. burgdorferi is polyploid and that the copies of the chromosome and plasmids are regularly spaced in each cell, which is critical for faithful segregation of the genome to daughter cells. Regular spacing of the chromosome is controlled by two separate partitioning systems that involve the protein pairs ParA/ParZ and ParB/Smc. Here, using chromosome conformation capture (Hi-C), we characterized the organization of the B. burgdorferi genome and the interactions between the replicons. We uncovered that although the linear chromosome lacks contacts between the two replication arms, the two telomeres are in frequent contact. Moreover, several plasmids specifically interact with the chromosome oriC region, and a subset of plasmids interact with each other more than with others. We found that Smc and the Smc-like MksB protein mediate long-range interactions on the chromosome, but they minimally affect plasmid-chromosome or plasmid-plasmid interactions. Finally, we found that disruption of the two partition systems leads to chromosome restructuring, correlating with the mis-positioning of chromosome oriC. Altogether, this study revealed the conformation of a complex genome and analyzed the contribution of the partition systems and SMC family proteins to this organization. This work expands the understanding of the organization and maintenance of multipartite bacterial genomes., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: C.J.-W. is an investigator of the Howard Hughes Medical Institute. H.B.B. is an employee of Illumina, Inc., (Copyright: © 2023 Ren et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The spirochete Borrelia burgdorferi, which causes Lyme disease, has the most segmented genome among known bacteria. In addition to a linear chromosome, the B. burgdorferi genome contains over 20 linear and circular endogenous plasmids. While many of these plasmids are dispensable under in vitro culture conditions, they are maintained during the natural life cycle of the pathogen. Plasmid-encoded functions are required for colonization of the tick vector, transmission to the vertebrate host, and evasion of host immune defenses. Different Borrelia strains can vary substantially in the type of plasmids they carry. The gene composition within the same type of plasmid can also differ from strain to strain, impeding the inference of plasmid function from one strain to another. To facilitate the investigation of the role of specific B. burgdorferi plasmids, we developed a Cas9-based approach that targets a plasmid for removal. As a proof-of-principle, we showed that targeting wild-type Cas9 to several loci on the endogenous plasmids lp25 or lp28-1 of the B. burgdorferi type strain B31 results in sgRNA-specific plasmid loss even when homologous sequences (i.e., potential sequence donors for DNA recombination) are present nearby. Cas9 nickase versions, Cas9D10A or Cas9H840A, also cause plasmid loss, though not as robustly. Thus, sgRNA-directed Cas9 DNA cleavage provides a highly efficient way to eliminate B. burgdorferi endogenous plasmids that are non-essential in axenic culture., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2022 Takacs et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Borrelia burgdorferi, the tick-transmitted spirochete agent of Lyme disease, has a highly segmented genome with a linear chromosome and various linear or circular plasmids. Here, by imaging several chromosomal loci and 16 distinct plasmids, we show that B. burgdorferi is polyploid during growth in culture and that the number of genome copies decreases during stationary phase. B. burgdorferi is also polyploid inside fed ticks and chromosome copies are regularly spaced along the spirochete's length in both growing cultures and ticks. This patterning involves the conserved DNA partitioning protein ParA whose localization is controlled by a potentially phage-derived protein, ParZ, instead of its usual partner ParB. ParZ binds its own coding region and acts as a centromere-binding protein. While ParA works with ParZ, ParB controls the localization of the condensin, SMC. Together, the ParA/ParZ and ParB/SMC pairs ensure faithful chromosome inheritance. Our findings underscore the plasticity of cellular functions, even those as fundamental as chromosome segregation., (© 2022. The Author(s).)

Adenosine triphosphate (ATP) is an abundant and essential metabolite that cells consume and regenerate in large amounts to support growth. Although numerous studies have inferred the intracellular concentration of ATP in bacterial cultures, what happens in individual bacterial cells under stable growth conditions is less clear. Here, we use the QUEEN-2m biosensor to quantify ATP dynamics in single Escherichia coli cells in relation to their growth rate, metabolism, cell cycle, and cell lineage. We find that ATP dynamics are more complex than expected from population studies and are associated with growth-rate variability. Under stable nutrient-rich condition, cells can display large fluctuations in ATP level that are partially coordinated with the cell cycle. Abrogation of aerobic acetate fermentation (overflow metabolism) through genetic deletion considerably reduces both the amplitude of ATP level fluctuations and the cell-cycle trend. Similarly, growth in media in which acetate fermentation is lower or absent results in the reduction of ATP level fluctuation and cell-cycle trend. This suggests that overflow metabolism exhibits temporal dynamics, which contributes to fluctuating ATP levels during growth. Remarkably, at the single-cell level, growth rate negatively correlates with the amplitude of ATP fluctuation for each tested condition, linking ATP dynamics to growth-rate heterogeneity in clonal populations. Our work highlights the importance of single-cell analysis in studying metabolism and its implication to phenotypic diversity and cell growth., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Microtubules are polarized intracellular polymers that play key roles in the cell, including in transport, polarity, and cell division. Across eukaryotic cell types, microtubules adopt diverse intracellular organization to accommodate these distinct functions coordinated by specific cellular sites called microtubule-organizing centers (MTOCs). Over 50 years of research on MTOC biology has focused mainly on the centrosome; however, most differentiated cells employ non-centrosomal MTOCs (ncMTOCs) to organize their microtubules into diverse arrays, which are critical to cell function. To identify essential ncMTOC components, we developed the biotin ligase-based, proximity-labeling approach TurboID for use in C. elegans. We identified proteins proximal to the microtubule minus end protein PTRN-1/Patronin at the apical ncMTOC of intestinal epithelial cells, focusing on two conserved proteins: spectraplakin protein VAB-10B/MACF1 and WDR-62, a protein we identify as homologous to vertebrate primary microcephaly disease protein WDR62. VAB-10B and WDR-62 do not associate with the centrosome and instead specifically regulate non-centrosomal microtubules and the apical targeting of microtubule minus-end proteins. Depletion of VAB-10B resulted in microtubule mislocalization and delayed localization of a microtubule nucleation complex ɣ-tubulin ring complex (γ-TuRC), while loss of WDR-62 decreased the number of dynamic microtubules and abolished γ-TuRC localization. This regulation occurs downstream of cell polarity and in conjunction with actin. As this is the first report for non-centrosomal roles of WDR62 family proteins, we expand the basic cell biological roles of this important disease protein. Our studies identify essential ncMTOC components and suggest a division of labor where microtubule growth and localization are distinctly regulated., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Interview with Christine Jacobs-Wagner, who studies the replication of bacterial cells at Stanford University., (Copyright © 2021.)

All cells fold their genomes, including bacterial cells, where the chromosome is compacted into a domain-organized meshwork called the nucleoid. How compaction and domain organization arise is not fully understood. Here, we describe a method to estimate the average mesh size of the nucleoid in Escherichia coli. Using nucleoid mesh size and DNA concentration estimates, we find that the cytoplasm behaves as a poor solvent for the chromosome when the cell is considered as a simple semidilute polymer solution. Monte Carlo simulations suggest that a poor solvent leads to chromosome compaction and DNA density heterogeneity (i.e., domain formation) at physiological DNA concentration. Fluorescence microscopy reveals that the heterogeneous DNA density negatively correlates with ribosome density within the nucleoid, consistent with cryoelectron tomography data. Drug experiments, together with past observations, suggest the hypothesis that RNAs contribute to the poor solvent effects, connecting chromosome compaction and domain formation to transcription and intracellular organization., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 Elsevier Inc. All rights reserved.)

The spirochete Borrelia burgdorferi causes Lyme disease, an increasingly prevalent infection. While previous studies have provided important insight into B. burgdorferi biology, many aspects, including basic cellular processes, remain underexplored. To help speed up the discovery process, we adapted a CRISPR interference (CRISPRi) platform for use in B. burgdorferi For efficiency and flexibility of use, we generated various CRISPRi template constructs that produce different basal and induced levels of dcas9 and carry different antibiotic resistance markers. We characterized the effectiveness of our CRISPRi platform by targeting the motility and cell morphogenesis genes flaB, mreB, rodA, and ftsI, whose native expression levels span two orders of magnitude. For all four genes, we obtained gene repression efficiencies of at least 95%. We showed by darkfield microscopy and cryo-electron tomography that flagellin (FlaB) depletion reduced the length and number of periplasmic flagella, which impaired cellular motility and resulted in cell straightening. Depletion of FtsI caused cell filamentation, implicating this protein in cell division in B. burgdorferi Finally, localized cell bulging in MreB- and RodA-depleted cells matched the locations of new peptidoglycan insertion specific to spirochetes of the Borrelia genus. These results therefore implicate MreB and RodA in the particular mode of cell wall elongation of these bacteria. Collectively, our results demonstrate the efficiency and ease of use of our B. burgdorferi CRISPRi platform, which should facilitate future genetic studies of this important pathogen. IMPORTANCE Gene function studies are facilitated by the availability of rapid and easy-to-use genetic tools. Homologous recombination-based methods traditionally used to genetically investigate gene function remain cumbersome to perform in B. burgdorferi , as they often are relatively inefficient. In comparison, our CRISPRi platform offers an easy and fast method to implement as it only requires a single plasmid transformation step and IPTG addition to obtain potent (>95%) downregulation of gene expression. To facilitate studies of various genes in wild-type and genetically modified strains, we provide over 30 CRISPRi plasmids that produce distinct levels of dcas9 expression and carry different antibiotic resistance markers. Our CRISPRi platform represents a useful and efficient complement to traditional genetic and chemical methods to study gene function in B. burgdorferi ., (Copyright © 2020 Takacs et al.)

Lyme disease, the most common vector-borne illness in North America, is caused by the spirochete Borrelia burgdorferi. Infection begins in the skin following a tick bite and can spread to the hearts, joints, nervous system, and other organs. Diverse host responses influence the level of B. burgdorferi infection in mice and humans. Using a systems biology approach, we examined potential molecular interactions between human extracellular and secreted proteins and B. burgdorferi. A yeast display library expressing 1031 human extracellular proteins was probed against 36 isolates of B. burgdorferi sensu lato. We found that human Peptidoglycan Recognition Protein 1 (PGLYRP1) interacted with the vast majority of B. burgdorferi isolates. In subsequent experiments, we demonstrated that recombinant PGLYRP1 interacts with purified B. burgdorferi peptidoglycan and exhibits borreliacidal activity, suggesting that vertebrate hosts may use PGLYRP1 to identify B. burgdorferi. We examined B. burgdorferi infection in mice lacking PGLYRP1 and observed an increased spirochete burden in the heart and joints, along with splenomegaly. Mice lacking PGLYRP1 also showed signs of immune dysregulation, including lower serum IgG levels and higher levels of IFNγ, CXCL9, and CXCL10.Taken together, our findings suggest that PGLYRP1 plays a role in the host's response to B. burgdorferi and further demonstrate the utility of expansive yeast display screening in capturing biologically relevant interactions between spirochetes and their hosts., Competing Interests: No authors have competing interests.

Exponentially growing systems are prevalent in nature, spanning all scales from biochemical reaction networks in single cells to food webs of ecosystems. How exponential growth emerges in nonlinear systems is mathematically unclear. Here, we describe a general theoretical framework that reveals underlying principles of long-term growth: scalability of flux functions and ergodicity of the rescaled systems. Our theory shows that nonlinear fluxes can generate not only balanced growth but also oscillatory or chaotic growth modalities, explaining nonequilibrium dynamics observed in cell cycles and ecosystems. Our mathematical framework is broadly useful in predicting long-term growth rates from natural and synthetic networks, analyzing the effects of system noise and perturbations, validating empirical and phenomenological laws on growth rate, and studying autocatalysis and network evolution., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

In scientific research, we often rely on well-established model systems to tackle important questions. In this context, extensive characterization of specific bacterial species such as Escherichia coli and Bacillus subtilis has provided a vast amount of knowledge that extends well beyond the biology of these two organisms. However, the bacterial world is large and extremely diverse, necessitating the development of additional models that complement the classical rod-shaped and symmetrically dividing systems. Caulobacter crescentus is a species that has met this need effectively, as its dimorphic lifestyle showcases distinctive features, including cellular asymmetry and differentiation during the cell cycle. Studying C. crescentus has reformed our understanding of bacterial intracellular organization, cellular development, and cell-cycle regulation. These findings have, in turn, stimulated studies in other bacteria, shedding light on how protein function and cell morphology can evolve and diversify. Studies in C. crescentus have also deepened our knowledge of other topics (e.g. cell mechanosensing, motility, and bacterial aging), while opening the door to biotechnological innovations. In this Primer, we provide some general background to this peculiar bacterium and highlight specific features that have contributed to its rise as a versatile bacterial model. This Primer is not meant to be exhaustive on any topic and is instead intended to provide a taste of the power of C. crescentus as a model system to explore a diverse range of topics., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Genes are often transcribed by multiple RNA polymerases (RNAPs) at densities that can vary widely across genes and environmental conditions. Here, we provide in vitro and in vivo evidence for a built-in mechanism by which co-transcribing RNAPs display either collaborative or antagonistic dynamics over long distances (>2 kb) through transcription-induced DNA supercoiling. In Escherichia coli, when the promoter is active, co-transcribing RNAPs translocate faster than a single RNAP, but their average speed is not altered by large variations in promoter strength and thus RNAP density. Environmentally induced promoter repression reduces the elongation efficiency of already-loaded RNAPs, causing premature termination and quick synthesis arrest of no-longer-needed proteins. This negative effect appears independent of RNAP convoy formation and is abrogated by topoisomerase I activity. Antagonistic dynamics can also occur between RNAPs from divergently transcribed gene pairs. Our findings may be broadly applicable given that transcription on topologically constrained DNA is the norm across organisms., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Lyme disease is a multisystem disorder caused by the spirochete Borrelia burgdorferi A common late-stage complication of this disease is oligoarticular arthritis, often involving the knee. In ∼10% of cases, arthritis persists after appropriate antibiotic treatment, leading to a proliferative synovitis typical of chronic inflammatory arthritides. Here, we provide evidence that peptidoglycan (PG), a major component of the B. burgdorferi cell envelope, may contribute to the development and persistence of Lyme arthritis (LA). We show that B. burgdorferi has a chemically atypical PG (PG Bb ) that is not recycled during cell-wall turnover. Instead, this pathogen sheds PG Bb fragments into its environment during growth. Patients with LA mount a specific immunoglobulin G response against PG Bb , which is significantly higher in the synovial fluid than in the serum of the same patient. We also detect PG Bb in 94% of synovial fluid samples (32 of 34) from patients with LA, many of whom had undergone oral and intravenous antibiotic treatment. These same synovial fluid samples contain proinflammatory cytokines, similar to those produced by human peripheral blood mononuclear cells stimulated with PG Bb In addition, systemic administration of PG Bb in BALB/c mice elicits acute arthritis. Altogether, our study identifies PG Bb as a likely contributor to inflammatory responses in LA. Persistence of this antigen in the joint may contribute to synovitis after antibiotics eradicate the pathogen. Furthermore, our finding that B. burgdorferi sheds immunogenic PG Bb fragments during growth suggests a potential role for PG Bb in the immunopathogenesis of other Lyme disease manifestations., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

The scaling of organelles with cell size is thought to be exclusive to eukaryotes. Here, we demonstrate that similar scaling relationships hold for the bacterial nucleoid. Despite the absence of a nuclear membrane, nucleoid size strongly correlates with cell size, independent of changes in DNA amount and across various nutrient conditions. This correlation is observed in diverse bacteria, revealing a near-constant ratio between nucleoid and cell size for a given species. As in eukaryotes, the nucleocytoplasmic ratio in bacteria varies greatly among species. This spectrum of nucleocytoplasmic ratios is independent of genome size, and instead it appears linked to the average population cell size. Bacteria with different nucleocytoplasmic ratios have a cytoplasm with different biophysical properties, impacting ribosome mobility and localization. Together, our findings identify new organizational principles and biophysical features of bacterial cells, implicating the nucleocytoplasmic ratio and cell size as determinants of the intracellular organization of translation., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Lyme disease is the most widely reported vector-borne disease in the United States. Its incidence is rapidly increasing, and disease symptoms can be debilitating. The need to understand the biology of the disease agent, the spirochete Borrelia burgdorferi , is thus evermore pressing. Despite important advances in B. burgdorferi genetics, the array of molecular tools available for use in this organism remains limited, especially for cell biological studies. Here, we adapt a palette of bright and mostly monomeric fluorescent proteins for versatile use and multicolor imaging in B. burgdorferi We also characterize two novel antibiotic selection markers and establish the feasibility of their use in conjunction with extant markers. Last, we describe a set of promoters of low and intermediate strengths that allow fine-tuning of gene expression levels. These molecular tools complement and expand current experimental capabilities in B. burgdorferi , which will facilitate future investigation of this important human pathogen. To showcase the usefulness of these reagents, we used them to investigate the subcellular localization of BB0323, a B. burgdorferi lipoprotein essential for survival in the host and vector environments. We show that BB0323 accumulates at the cell poles and future division sites of B. burgdorferi cells, highlighting the complex subcellular organization of this spirochete. IMPORTANCE Genetic manipulation of the Lyme disease spirochete B. burgdorferi remains cumbersome, despite significant progress in the field. The scarcity of molecular reagents available for use in this pathogen has slowed research efforts to study its unusual biology. Of interest, B. burgdorferi displays complex cellular organization features that have yet to be understood. These include an unusual morphology and a highly fragmented genome, both of which are likely to play important roles in the bacterium's transmission, infectivity, and persistence. Here, we complement and expand the array of molecular tools available for use in B. burgdorferi by generating and characterizing multiple fluorescent proteins, antibiotic selection markers, and promoters of varied strengths. These tools will facilitate investigations in this important human pathogen, as exemplified by the polar and midcell localization of the cell envelope regulator BB0323, which we uncovered using these reagents., (Copyright © 2018 Takacs et al.)

We describe a general computational approach to designing self-assembling helical filaments from monomeric proteins and use this approach to design proteins that assemble into micrometer-scale filaments with a wide range of geometries in vivo and in vitro. Cryo-electron microscopy structures of six designs are close to the computational design models. The filament building blocks are idealized repeat proteins, and thus the diameter of the filaments can be systematically tuned by varying the number of repeat units. The assembly and disassembly of the filaments can be controlled by engineered anchor and capping units built from monomers lacking one of the interaction surfaces. The ability to generate dynamic, highly ordered structures that span micrometers from protein monomers opens up possibilities for the fabrication of new multiscale metamaterials., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Cell size, cell growth, and cell cycle events are necessarily intertwined to achieve robust bacterial replication. Yet, a comprehensive and integrated view of these fundamental processes is lacking. Here, we describe an image-based quantitative screen of the single-gene knockout collection of Escherichia coli and identify many new genes involved in cell morphogenesis, population growth, nucleoid (bulk chromosome) dynamics, and cell division. Functional analyses, together with high-dimensional classification, unveil new associations of morphological and cell cycle phenotypes with specific functions and pathways. Additionally, correlation analysis across ~4,000 genetic perturbations shows that growth rate is surprisingly not predictive of cell size. Growth rate was also uncorrelated with the relative timings of nucleoid separation and cell constriction. Rather, our analysis identifies scaling relationships between cell size and nucleoid size and between nucleoid size and the relative timings of nucleoid separation and cell division. These connections suggest that the nucleoid links cell morphogenesis to the cell cycle., (© 2018 The Authors. Published under the terms of the CC BY 4.0 license.)

Genetically identical cells exhibit diverse phenotypes even when experiencing the same environment. This phenomenon in part originates from cell-to-cell variability (noise) in protein expression. Although various kinetic schemes of stochastic transcription initiation are known to affect gene expression noise, how posttranscription initiation events contribute to noise at the protein level remains incompletely understood. To address this question, we developed a stochastic simulation-based model of bacterial gene expression that integrates well-known dependencies between transcription initiation, transcription elongation dynamics, mRNA degradation, and translation. We identified realistic conditions under which mRNA lifetime and transcriptional pauses modulate the protein expression noise initially introduced by the promoter architecture. For instance, we found that the short lifetime of bacterial mRNAs facilitates the production of protein bursts. Conversely, RNA polymerase (RNAP) pausing at specific sites during transcription elongation can attenuate protein bursts by fluidizing the RNAP traffic to the point of erasing the effect of a bursty promoter. Pause-prone sites, if located close to the promoter, can also affect noise indirectly by reducing both transcription and translation initiation due to RNAP and ribosome congestion. Our findings highlight how the interplay between transcription initiation, transcription elongation, translation, and mRNA degradation shapes the distribution in protein numbers. They also have implications for our understanding of gene evolution and suggest combinatorial strategies for modulating phenotypic variability by genetic engineering., (Copyright © 2018 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Spatial organization is a hallmark of all living systems. Even bacteria, the smallest forms of cellular life, display defined shapes and complex internal organization, showcasing a highly structured genome, cytoskeletal filaments, localized scaffolding structures, dynamic spatial patterns, active transport, and occasionally, intracellular organelles. Spatial order is required for faithful and efficient cellular replication and offers a powerful means for the development of unique biological properties. Here, we discuss organizational features of bacterial cells and highlight how bacteria have evolved diverse spatial mechanisms to overcome challenges cells face as self-replicating entities., (Copyright © 2018 Elsevier Inc. All rights reserved.)

To achieve robust replication, bacteria must integrate cellular metabolism and cell wall growth. While these two processes have been well characterized, the nature and extent of cross-regulation between them is not well understood. Here, using classical genetics, CRISPRi, metabolomics, transcriptomics and chemical complementation approaches, we show that a loss of the master regulator Hfq in Caulobacter crescentus alters central metabolism and results in cell shape defects in a nutrient-dependent manner. We demonstrate that the cell morphology phenotype in the hfq deletion mutant is attributable to a disruption of α-ketoglutarate (KG) homeostasis. In addition to serving as a key intermediate of the tricarboxylic acid (TCA) cycle, KG is a by-product of an enzymatic reaction required for the synthesis of peptidoglycan, a major component of the bacterial cell wall. Accumulation of KG in the hfq deletion mutant interferes with peptidoglycan synthesis, resulting in cell morphology defects and increased susceptibility to peptidoglycan-targeting antibiotics. This work thus reveals a direct crosstalk between the TCA cycle and cell wall morphogenesis. This crosstalk highlights the importance of metabolic homeostasis in not only ensuring adequate availability of biosynthetic precursors, but also in preventing interference with cellular processes in which these intermediates arise as by-products.

New strategies are needed to combat antibiotic resistance, especially against pathogens such as methicillin-resistant Staphylococcus aureus A tick antifreeze glycoprotein, IAFGP, possesses potent antibiofilm properties against a variety of clinical pathogens, including S. aureus Synergy between IAFGP, or a peptide (P1) representative of a repeat region of the protein, with different antibiotics was assessed in vitro Antibiotics that synergized with either IAFPG or P1 were further evaluated in vivo using vertebrate and invertebrate infection models. IAFGP readily enhanced the efficacy of antibiotics against S. aureus Synergy with daptomycin, an antibiotic used to treat methicillin-resistant S. aureus , was observed in vitro and in vivo using iafgp -transgenic mice and flies. Furthermore, synergy with ciprofloxacin or gentamicin, antibiotics not generally used to treat S. aureus , was also perceived. The combined effect of the antibiotic and IAFGP was associated with improved permeation of the antibiotic into the cell. Our results highlight that synergy of IAFGP with antibiotics traditionally used to treat this pathogen, and enhancement of the potency of antibiotics not commonly used against this microbe, can provide novel alternative therapeutic strategies to combat bacterial infections., (Copyright © 2017 American Society for Microbiology.)

In bacteria, chromosome dynamics and gene expression are modulated by nucleoid-associated proteins (NAPs), but little is known about how NAP activity is coupled to cell cycle progression. Using genomic techniques, quantitative cell imaging, and mathematical modeling, our study in Caulobacter crescentus identifies a novel NAP (GapR) whose activity over the cell cycle is shaped by DNA replication. GapR activity is critical for cellular function, as loss of GapR causes severe, pleiotropic defects in growth, cell division, DNA replication, and chromosome segregation. GapR also affects global gene expression with a chromosomal bias from origin to terminus, which is associated with a similar general bias in GapR binding activity along the chromosome. Strikingly, this asymmetric localization cannot be explained by the distribution of GapR binding sites on the chromosome. Instead, we present a mechanistic model in which the spatiotemporal dynamics of GapR are primarily driven by the progression of the replication forks. This model represents a simple mechanism of cell cycle regulation, in which DNA-binding activity is intimately linked to the action of DNA replication., (© 2016 The Authors.)

Arthropods transmit diverse infectious agents; however, the ways microbes influence their vector to enhance colonization are poorly understood. Ixodes scapularis ticks harbor numerous human pathogens, including Anaplasma phagocytophilum, the agent of human granulocytic anaplasmosis. We now demonstrate that A. phagocytophilum modifies the I. scapularis microbiota to more efficiently infect the tick. A. phagocytophilum induces ticks to express Ixodes scapularis antifreeze glycoprotein (iafgp), which encodes a protein with several properties, including the ability to alter bacterial biofilm formation. IAFGP thereby perturbs the tick gut microbiota, which influences the integrity of the peritrophic matrix and gut barrier-critical obstacles for Anaplasma colonization. Mechanistically, IAFGP binds the terminal d-alanine residue of the pentapeptide chain of bacterial peptidoglycan, resulting in altered permeability and the capacity of bacteria to form biofilms. These data elucidate the molecular mechanisms by which a human pathogen appropriates an arthropod antibacterial protein to alter the gut microbiota and more effectively colonize the vector., Competing Interests: The authors declare no conflict of interest.

Spatial ordering of macromolecular components inside cells is important for cellular physiology and replication. In bacteria, ParA/B systems are known to generate various intracellular patterns that underlie the transport and partitioning of low-copy-number cargos such as plasmids. ParA/B systems consist of ParA, an ATPase that dimerizes and binds DNA upon ATP binding, and ParB, a protein that binds the cargo and stimulates ParA ATPase activity. Inside cells, ParA is asymmetrically distributed, forming a propagating wave that is followed by the ParB-rich cargo. These correlated dynamics lead to cargo oscillation or equidistant spacing over the nucleoid depending on whether the cargo is in single or multiple copies. Currently, there is no model that explains how these different spatial patterns arise and relate to each other. Here, we test a simple DNA-relay model that has no imposed asymmetry and that only considers the ParA/ParB biochemistry and the known fluctuating and elastic dynamics of chromosomal loci. Stochastic simulations with experimentally derived parameters demonstrate that this model is sufficient to reproduce the signature patterns of ParA/B systems: the propagating ParA gradient correlated with the cargo dynamics, the single-cargo oscillatory motion, and the multicargo equidistant patterning. Stochasticity of ATP hydrolysis breaks the initial symmetry in ParA distribution, resulting in imbalance of elastic force acting on the cargo. Our results may apply beyond ParA/B systems as they reveal how a minimal system of two players, one binding to DNA and the other modulating this binding, can transform directionally random DNA fluctuations into directed motion and intracellular patterning., Competing Interests: The authors declare no conflict of interest.

Agents that cause Lyme disease, relapsing fever, leptospirosis, and syphilis belong to the phylum Spirochaetae-a unique lineage of bacteria most known for their long, spiral morphology. Despite the relevance to human health, little is known about the most fundamental aspects of spirochete growth. Here, using quantitative microscopy to track peptidoglycan cell-wall synthesis, we found that the Lyme disease spirochete Borrelia burgdorferi displays a complex pattern of growth. B. burgdorferi elongates from discrete zones that are both spatially and temporally regulated. In addition, some peptidoglycan incorporation occurs along the cell body, with the notable exception of a large region at the poles. Newborn cells inherit a highly active zone of peptidoglycan synthesis at midcell that contributes to elongation for most of the cell cycle. Concomitant with the initiation of nucleoid separation and cell constriction, second and third zones of elongation are established at the 1/4 and 3/4 cellular positions, marking future sites of division for the subsequent generation. Positioning of elongation zones along the cell is robust to cell length variations and is relatively precise over long distances (>30 µm), suggesting that cells ‟sense" relative, as opposed to absolute, cell length to establish zones of peptidoglycan synthesis. The transition from one to three zones of peptidoglycan growth during the cell cycle is also observed in relapsing fever Borrelia. However, this mode of growth does not extend to representative species from other spirochetal genera, suggesting that this distinctive growth mode represents an evolutionary divide in the spirochete phylum., Competing Interests: The authors declare no conflict of interest. See QnAs on page 9129.

Fluorescence nanoscopy, or super-resolution microscopy, has become an important tool in cell biological research. However, because of its usually inferior resolution in the depth direction (50-80 nm) and rapidly deteriorating resolution in thick samples, its practical biological application has been effectively limited to two dimensions and thin samples. Here, we present the development of whole-cell 4Pi single-molecule switching nanoscopy (W-4PiSMSN), an optical nanoscope that allows imaging of three-dimensional (3D) structures at 10- to 20-nm resolution throughout entire mammalian cells. We demonstrate the wide applicability of W-4PiSMSN across diverse research fields by imaging complex molecular architectures ranging from bacteriophages to nuclear pores, cilia, and synaptonemal complexes in large 3D cellular volumes., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

In bacteria, ParABS systems mediate intracellular transport of various cargos, including chromosomal regions in Caulobacter crescentus. Transport of the ParB/parS partition complex requires the DNA-binding activity of ParA, which transiently tethers the partition complex during translocation. In C. crescentus, the directionality of the transport is set up by a gradient of ParA whose concentration gradually increases from one end of the cell (old pole) to the other (new pole). Importantly, this ParA gradient is already observed before DNA replication and segregation are initiated when the partition complex is anchored at the old pole. How such micron-scale ParA pattern is established and maintained before the initiation of chromosome segregation has not been experimentally established. Although the stimulation of ParA ATPase activity by the localized ParB/parS partition complex is thought to be involved, this activity alone cannot quantitatively describe the ParA pattern observed inside cells. Instead, our experimental and theoretical study shows that the missing key component for achieving the experimentally observed steady-state ParA patterning is the slow mobility of ParA dimers (D ∼10(-3)μm(2)/s) due to intermittent DNA binding. Our model recapitulates the entire steady-state ParA distribution observed experimentally, including the shape of the gradient as well as ParA accumulation at the location of the partition complex. Stochastic simulations suggest that cell-to-cell variability in ParA pattern is due to the low ParA copy number in C. crescentus cells. The model also accounts for an apparent exclusion of ParA from regions with small spacing between partition complexes observed in filamentous cells. Collectively, our work demonstrates that in addition to its function in mediating transport, the conserved DNA-binding property of ParA has a critical function before DNA segregation by setting up a ParA pattern required for transport directionality., (Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

With the realization that bacteria display phenotypic variability among cells and exhibit complex subcellular organization critical for cellular function and behavior, microscopy has re-emerged as a primary tool in bacterial research during the last decade. However, the bottleneck in today's single-cell studies is quantitative image analysis of cells and fluorescent signals. Here, we address current limitations through the development of Oufti, a stand-alone, open-source software package for automated measurements of microbial cells and fluorescence signals from microscopy images. Oufti provides computational solutions for tracking touching cells in confluent samples, handles various cell morphologies, offers algorithms for quantitative analysis of both diffraction and non-diffraction-limited fluorescence signals and is scalable for high-throughput analysis of massive datasets, all with subpixel precision. All functionalities are integrated in a single package. The graphical user interface, which includes interactive modules for segmentation, image analysis and post-processing analysis, makes the software broadly accessible to users irrespective of their computational skills., (© 2015 John Wiley & Sons Ltd.)

Over 3 billion years ago, cellular life began anaerobically. A new study now establishes a key link between oxidative stress and proliferation of wall-less bacteria known as L-forms. The finding provides insights into both the origin of life and the potential threat posed by pathogenic L-forms., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Horizontal gene transfer allows organisms to rapidly acquire adaptive traits. Although documented instances of horizontal gene transfer from bacteria to eukaryotes remain rare, bacteria represent a rich source of new functions potentially available for co-option. One benefit that genes of bacterial origin could provide to eukaryotes is the capacity to produce antibacterials, which have evolved in prokaryotes as the result of eons of interbacterial competition. The type VI secretion amidase effector (Tae) proteins are potent bacteriocidal enzymes that degrade the cell wall when delivered into competing bacterial cells by the type VI secretion system. Here we show that tae genes have been transferred to eukaryotes on at least six occasions, and that the resulting domesticated amidase effector (dae) genes have been preserved for hundreds of millions of years through purifying selection. We show that the dae genes acquired eukaryotic secretion signals, are expressed within recipient organisms, and encode active antibacterial toxins that possess substrate specificity matching extant Tae proteins of the same lineage. Finally, we show that a dae gene in the deer tick Ixodes scapularis limits proliferation of Borrelia burgdorferi, the aetiologic agent of Lyme disease. Our work demonstrates that a family of horizontally acquired toxins honed to mediate interbacterial antagonism confers previously undescribed antibacterial capacity to eukaryotes. We speculate that the selective pressure imposed by competition between bacteria has produced a reservoir of genes encoding diverse antimicrobial functions that are tailored for co-option by eukaryotic innate immune systems.

Muscodor albus belongs to a genus of endophytic fungi that inhibit and kill other fungi, bacteria, and insects through production of a complex mixture of volatile organic compounds (VOCs). This process of mycofumigation has found commercial application for control of human and plant pathogens, but the mechanism of the VOC toxicity is unknown. Here, the mode of action of these volatiles was investigated through a series of genetic screens and biochemical assays. A single-gene knockout screen revealed high sensitivity for Escherichia coli lacking enzymes in the pathways of DNA repair, DNA metabolic process, and response to stress when exposed to the VOCs of M. albus. Furthermore, the sensitivity of knockouts involved in the repair of specific DNA alkyl adducts suggests that the VOCs may induce alkylation. Evidence of DNA damage suggests that these adducts lead to breaks during DNA replication or transcription if not properly repaired. Additional cytotoxicity profiling indicated that during VOC exposure, E. coli became filamentous and demonstrated an increase in cellular membrane fluidity. The volatile nature of the toxic compounds produced by M. albus and their broad range of inhibition make this fungus an attractive biological agent. Understanding the antimicrobial effects and the VOC mode of action will inform the utility and safety of potential mycofumigation applications for M. albus., (Copyright © 2015, American Society for Microbiology. All Rights Reserved.)

Cell size control is an intrinsic feature of the cell cycle. In bacteria, cell growth and division are thought to be coupled through a cell size threshold. Here, we provide direct experimental evidence disproving the critical size paradigm. Instead, we show through single-cell microscopy and modeling that the evolutionarily distant bacteria Escherichia coli and Caulobacter crescentus achieve cell size homeostasis by growing, on average, the same amount between divisions, irrespective of cell length at birth. This simple mechanism provides a remarkably robust cell size control without the need of being precise, abating size deviations exponentially within a few generations. This size homeostasis mechanism is broadly applicable for symmetric and asymmetric divisions, as well as for different growth rates. Furthermore, our data suggest that constant size extension is implemented at or close to division. Altogether, our findings provide fundamentally distinct governing principles for cell size and cell-cycle control in bacteria., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Several intracellular pathogens, such as Brucella abortus, display a biphasic infection process starting with a non-proliferative stage of unclear nature. Here, we study the cell cycle of B. abortus at the single-cell level, in culture and during infection of HeLa cells and macrophages. The localization of segregation and replication loci of the two bacterial chromosomes indicates that, immediately after being engulfed by host-cell endocytic vacuoles, most bacterial cells are newborn. These bacterial cells do not initiate DNA replication for the next 4 to 6 h, indicating a G1 arrest. Moreover, growth is completely stopped during that time, reflecting a global cell cycle block. Growth and DNA replication resume later, although bacteria still reside within endosomal-like compartments. We hypothesize that the predominance of G1-arrested bacteria in the infectious population, and the bacterial cell cycle arrest following internalization, may constitute a widespread strategy among intracellular pathogens to colonize new proliferation niches.

The widely conserved ParABS system plays a major role in bacterial chromosome segregation. How the components of this system work together to generate translocation force and directional motion remains uncertain. Here, we combine biochemical approaches, quantitative imaging and mathematical modeling to examine the mechanism by which ParA drives the translocation of the ParB/parS partition complex in Caulobacter crescentus. Our experiments, together with simulations grounded on experimentally-determined biochemical and cellular parameters, suggest a novel 'DNA-relay' mechanism in which the chromosome plays a mechanical function. In this model, DNA-bound ParA-ATP dimers serve as transient tethers that harness the elastic dynamics of the chromosome to relay the partition complex from one DNA region to another across a ParA-ATP dimer gradient. Since ParA-like proteins are implicated in the partitioning of various cytoplasmic cargos, the conservation of their DNA-binding activity suggests that the DNA-relay mechanism may be a general form of intracellular transport in bacteria.DOI: http://dx.doi.org/10.7554/eLife.02758.001., (Copyright © 2014, Lim et al.)

The physical nature of the bacterial cytoplasm is poorly understood even though it determines cytoplasmic dynamics and hence cellular physiology and behavior. Through single-particle tracking of protein filaments, plasmids, storage granules, and foreign particles of different sizes, we find that the bacterial cytoplasm displays properties that are characteristic of glass-forming liquids and changes from liquid-like to solid-like in a component size-dependent fashion. As a result, the motion of cytoplasmic components becomes disproportionally constrained with increasing size. Remarkably, cellular metabolism fluidizes the cytoplasm, allowing larger components to escape their local environment and explore larger regions of the cytoplasm. Consequently, cytoplasmic fluidity and dynamics dramatically change as cells shift between metabolically active and dormant states in response to fluctuating environments. Our findings provide insight into bacterial dormancy and have broad implications to our understanding of bacterial physiology, as the glassy behavior of the cytoplasm impacts all intracellular processes involving large components., (Copyright © 2014 Elsevier Inc. All rights reserved.)

It is now well appreciated that bacterial cells are highly organized, which is far from the initial concept that they are merely bags of randomly distributed macromolecules and chemicals. Central to their spatial organization is the precise positioning of certain proteins in subcellular domains of the cell. In particular, the cell poles - the ends of rod-shaped cells - constitute important platforms for cellular regulation that underlie processes as essential as cell cycle progression, cellular differentiation, virulence, chemotaxis and growth of appendages. Thus, understanding how the polar localization of specific proteins is achieved and regulated is a crucial question in bacterial cell biology. Often, polarly localized proteins are recruited to the poles through their interaction with other proteins or protein complexes that were already located there, in a so-called diffusion-and-capture mechanism. Bacteria are also starting to reveal their secrets on how the initial pole 'recognition' can occur and how this event can be regulated to generate dynamic, reproducible patterns in time (for example, during the cell cycle) and space (for example, at a specific cell pole). Here, we review the major mechanisms that have been described in the literature, with an emphasis on the self-organizing principles. We also present regulation strategies adopted by bacterial cells to obtain complex spatiotemporal patterns of protein localization.

While translational read-through of stop codons by suppressor tRNAs is common in many bacteria, archaea and eukaryotes, this phenomenon has not yet been observed in the α-proteobacterium Caulobacter crescentus. Based on a previous report that C. crescentus and Escherichia coli tRNA(His) have distinctive identity elements, we constructed E. coli tRNA(His) CUA, a UAG suppressor tRNA for C. crescentus. By examining the expression of three UAG codon- containing reporter genes (encoding a β-lactamase, the fluorescent mCherry protein, or the C. crescentus xylonate dehydratase), we demonstrated that the E. coli histidyl-tRNA synthetase/tRNA(His) CUA pair enables in vivo UAG suppression in C. crescentus. E. coli histidyl-tRNA synthetase (HisRS) or tRNA(His) CUA alone did not achieve suppression; this indicates that the E. coli HisRS/tRNA(His) CUA pair is orthogonal in C. crescentus. These results illustrate that UAG suppression can be achieved in C. crescentus with an orthogonal aminoacyl-tRNA synthetase/suppressor tRNA pair.

Background: The genetic network involved in the bacterial cell cycle is poorly understood even though it underpins the remarkable ability of bacteria to proliferate. How such network evolves is even less clear. The major aims of this work were to identify and examine the genes and pathways that are differentially expressed during the Caulobacter crescentus cell cycle, and to analyze the evolutionary features of the cell cycle network., Results: We used deep RNA sequencing to obtain high coverage RNA-Seq data of five C. crescentus cell cycle stages, each with three biological replicates. We found that 1,586 genes (over a third of the genome) display significant differential expression between stages. This gene list, which contains many genes previously unknown for their cell cycle regulation, includes almost half of the genes involved in primary metabolism, suggesting that these "house-keeping" genes are not constitutively transcribed during the cell cycle, as often assumed. Gene and module co-expression clustering reveal co-regulated pathways and suggest functionally coupled genes. In addition, an evolutionary analysis of the cell cycle network shows a high correlation between co-expression and co-evolution. Most co-expression modules have strong phylogenetic signals, with broadly conserved genes and clade-specific genes predominating different substructures of the cell cycle co-expression network. We also found that conserved genes tend to determine the expression profile of their module., Conclusion: We describe the first phylogenetic and single-nucleotide-resolution transcriptomic analysis of a bacterial cell cycle network. In addition, the study suggests how evolution has shaped this network and provides direct biological network support that selective pressure is not on individual genes but rather on the relationship between genes, which highlights the importance of integrating phylogenetic analysis into biological network studies.