Hard-Wired Control of Bacterial Processes by Chromosomal Gene Location

On the evolution of chromosomal regions with high gene strand bias in bacteria

On circular bacterial chromosomes, the majority of genes are coded on the leading strand. This gene strand bias (GSB) ranges from up to 85% in some Bacillota to a little more than 50% in other phyla. The factors determining the extent of the strand bias remain to be found. Here, we report that species in the phylum Gemmatimonadota share a unique chromosome architecture, distinct from neighboring phyla: in a conserved 600-kb region around the terminus of replication, almost all genes were located on the leading strands, while on the remaining part of the chromosome, the strand preference was more balanced. The high strand bias (HSB) region harbors the rRNA clusters, core, and highly expressed genes. Selective pressure for reduction of collisions with DNA replication to minimize detrimental mutations can explain the conservation of essential genes in this region. Repetitive and mobile elements are underrepresented, suggesting reduced recombination frequency by structural isolation from other parts of the chromosome. We propose that the HSB region forms a distinct chromosomal domain. Gemmatimonadota chromosomes evolved mainly by expansion through horizontal gene transfer and duplications outside of the ancient high strand bias region. In support of our hypothesis, we could further identify two Spiroplasma strains on a similar evolutionary path.IMPORTANCEOn bacterial chromosomes, a preferred location of genes on the leading strand has evolved to reduce conflicts between replication and transcription. Despite a vast body of research, the question why bacteria show large differences in their gene strand bias is still not solved. The discovery of "hybrid" chromosomes in different phyla, including Gemmatimonadota, in which a conserved high strand bias is found exclusively in a region at ter, points toward a role of nucleoid structure, additional to replication, in the evolution of strand preferences. A fine-grained structural analysis of the ever-increasing number of available bacterial genomes could help to better understand the forces that shape the sequential and spatial organization of the cell's information content.

Regulation of biofilm gene expression by DNA replication in Bacillus subtilis

Bacillus subtilis relies on biofilms for survival in harsh environments. Extracellular polymeric substance (EPS) is a crucial component of biofilms, yet the dynamics of EPS production in single cells remain elusive. To unveil the modulation of EPS synthesis, we built a minimal network model comprising the SinI-SinR-SlrR module, Spo0A, and EPS. Stochastic simulations revealed that antagonistic interplay between SinI and SinR enables EPS production in bursts. SlrR widens these bursts and increases their frequency by stabilizing SinR-SlrR complexes and depleting free SinR. DNA replication and chromosomal positioning of key genes dictate pulsatile changes in the slrR:sinR gene dosage ratio (gr) and Spo0A-P levels, each promoting EPS production in distinct phases of the cell cycle. As the cell cycle lengthens with nutrient stress, the duty cycle of gr pulsing decreases, whereas the amplitude of Spo0A-P pulses elevates. This coordinated response facilitates keeping a constant proportion of EPS-secreting cells within colonies across diverse nutrient conditions. Our results suggest that bacteria may 'encode' eps expression through strategic chromosomal organization. This work illuminates how stochastic protein interactions, gene copy number imbalance, and cell-cycle dynamics orchestrate EPS synthesis, offering a deeper understanding of biofilm formation.

A method to correct for local alterations in DNA copy number that bias functional genomics assays applied to antibiotic-treated bacteria

Functional genomics techniques, such as transposon insertion sequencing and RNA-sequencing, are key to studying relative differences in bacterial mutant fitness or gene expression under selective conditions. However, certain stress conditions, mutations, or antibiotics can directly interfere with DNA synthesis, resulting in systematic changes in local DNA copy numbers along the chromosome. This can lead to artifacts in sequencing-based functional genomics data when comparing antibiotic treatment to an unstressed control. Further, relative differences in gene-wise read counts may result from alterations in chromosomal replication dynamics, rather than selection or direct gene regulation. We term this artifact "chromosomal location bias" and implement a principled statistical approach to correct it by calculating local normalization factors along the chromosome. These normalization factors are then directly incorporated into statistical analyses using standard RNA-sequencing analysis methods without modifying the read counts themselves, preserving important information about the mean-variance relationship in the data. We illustrate the utility of this approach by generating and analyzing a ciprofloxacin-treated transposon insertion sequencing data set in Escherichia coli as a case study. We show that ciprofloxacin treatment generates chromosomal location bias in the resulting data, and we further demonstrate that failing to correct for this bias leads to false predictions of mutant drug sensitivity as measured by minimum inhibitory concentrations. We have developed an R package and user-friendly graphical Shiny application, ChromoCorrect, that detects and corrects for chromosomal bias in read count data, enabling the application of functional genomics technologies to the study of antibiotic stress.IMPORTANCEAltered gene dosage due to changes in DNA replication has been observed under a variety of stresses with a variety of experimental techniques. However, the implications of changes in gene dosage for sequencing-based functional genomics assays are rarely considered. We present a statistically principled approach to correcting for the effect of changes in gene dosage, enabling testing for differences in the fitness effects or regulation of individual genes in the presence of confounding differences in DNA copy number. We show that failing to correct for these effects can lead to incorrect predictions of resistance phenotype when applying functional genomics assays to investigate antibiotic stress, and we provide a user-friendly application to detect and correct for changes in DNA copy number.

Transcription-replication interactions reveal bacterial genome regulation

Organisms determine the transcription rates of thousands of genes through a few modes of regulation that recur across the genome1. In bacteria, the relationship between the regulatory architecture of a gene and its expression is well understood for individual model gene circuits2,3. However, a broader perspective of these dynamics at the genome scale is lacking, in part because bacterial transcriptomics has hitherto captured only a static snapshot of expression averaged across millions of cells4. As a result, the full diversity of gene expression dynamics and their relation to regulatory architecture remains unknown. Here we present a novel genome-wide classification of regulatory modes based on the transcriptional response of each gene to its own replication, which we term the transcription-replication interaction profile (TRIP). Analysing single-bacterium RNA-sequencing data, we found that the response to the universal perturbation of chromosomal replication integrates biological regulatory factors with biophysical molecular events on the chromosome to reveal the local regulatory context of a gene. Whereas the TRIPs of many genes conform to a gene dosage-dependent pattern, others diverge in distinct ways, and this is shaped by factors such as intra-operon position and repression state. By revealing the underlying mechanistic drivers of gene expression heterogeneity, this work provides a quantitative, biophysical framework for modelling replication-dependent expression dynamics.

EVRC: reconstruction of chromosome 3D structure models using error-vector resultant algorithm with clustering coefficient

MOTIVATION

Reconstruction of 3D structure models is of great importance for the study of chromosome function. Software tools for this task are highly needed.

RESULTS

We present a novel reconstruction algorithm, called EVRC, which utilizes co-clustering coefficients and error-vector resultant for chromosome 3D structure reconstruction. As an update of our previous EVR algorithm, EVRC now can deal with both single and multiple chromosomes in structure modeling. To evaluate the effectiveness and accuracy of the EVRC algorithm, we applied it to simulation datasets and real Hi-C datasets. The results show that the reconstructed structures have high similarity to the original/real structures, indicating the effectiveness and robustness of the EVRC algorithm. Furthermore, we applied the algorithm to the 3D conformation reconstruction of the wild-type and mutant Arabidopsis thaliana chromosomes and demonstrated the differences in structural characteristics between different chromosomes. We also accurately showed the conformational change in the centromere region of the mutant compared with the wild-type of Arabidopsis chromosome 1. Our EVRC algorithm is a valuable software tool for the field of chromatin structure reconstruction, and holds great promise for advancing our understanding on the chromosome functions.

AVAILABILITY AND IMPLEMENTATION

The software is available at https://github.com/mbglab/EVRC.

Ploidy in Vibrio natriegens: Very Dynamic and Rapidly Changing Copy Numbers of Both Chromosomes

Vibrio natriegens is the fastest-growing bacterium, with a doubling time of approximately 12-14 min. It has a high potential for basic research and biotechnological applications, e.g., it can be used for the cell-free production of (labeled) heterologous proteins, for synthetic biological applications, and for the production of various compounds. However, the ploidy level in V. natriegens remains unknown. At nine time points throughout the growth curve, we analyzed the numbers of origins and termini of both chromosomes with qPCR and the relative abundances of all genomic sites with marker frequency analyses. During the lag phase until early exponential growth, the origin copy number and origin/terminus ratio of chromosome 1 increased severalfold, but the increase was lower for chromosome 2. This increase was paralleled by an increase in cell volume. During the exponential phase, the origin/terminus ratio and cell volume decreased again. This highly dynamic and fast regulation has not yet been described for any other species. In this study, the gene dosage increase in origin-adjacent genes during the lag phase is discussed together with the nonrandom distribution of genes on the chromosomes of V. natriegens. Taken together, the results of this study provide the first comprehensive overview of the chromosome dynamics in V. natriegens and will guide the optimization of molecular biological characterization and biotechnological applications.

Chromosomal Position of Ribosomal Protein Genes Affects Long-Term Evolution of Vibrio cholerae

It is unclear how gene order within the chromosome influences genome evolution. Bacteria cluster transcription and translation genes close to the replication origin (oriC). In Vibrio cholerae, relocation of s10-spc-α locus (S10), the major locus of ribosomal protein genes, to ectopic genomic positions shows that its relative distance to the oriC correlates to a reduction in growth rate, fitness, and infectivity. To test the long-term impact of this trait, we evolved 12 populations of V. cholerae strains bearing S10 at an oriC-proximal or an oriC-distal location for 1,000 generations. During the first 250 generations, positive selection was the main force driving mutation. After 1,000 generations, we observed more nonadaptative mutations and hypermutator genotypes. Populations fixed inactivating mutations at many genes linked to virulence: flagellum, chemotaxis, biofilm, and quorum sensing. Throughout the experiment, all populations increased their growth rates. However, those bearing S10 close to oriC remained the fittest, indicating that suppressor mutations cannot compensate for the genomic position of the main ribosomal protein locus. Selection and sequencing of the fastest-growing clones allowed us to characterize mutations inactivating, among other sites, flagellum master regulators. Reintroduction of these mutations into the wild-type context led to a ≈10% growth improvement. In conclusion, the genomic location of ribosomal protein genes conditions the evolutionary trajectory of V. cholerae. While genomic content is highly plastic in prokaryotes, gene order is an underestimated factor that conditions cellular physiology and evolution. A lack of suppression enables artificial gene relocation as a tool for genetic circuit reprogramming. IMPORTANCE The bacterial chromosome harbors several entangled processes such as replication, transcription, DNA repair, and segregation. Replication begins bidirectionally at the replication origin (oriC) until the terminal region (ter) organizing the genome along the ori-ter axis gene order along this axis could link genome structure to cell physiology. Fast-growing bacteria cluster translation genes near oriC. In Vibrio cholerae, moving them away was feasible but at the cost of losing fitness and infectivity. Here, we evolved strains harboring ribosomal genes close or far from oriC. Growth rate differences persisted after 1,000 generations. No mutation was able to compensate for the growth defect, showing that ribosomal gene location conditions their evolutionary trajectory. Despite the high plasticity of bacterial genomes, evolution has sculpted gene order to optimize the ecological strategy of the microorganism. We observed growth rate improvement throughout the evolution experiment that occurred at expense of energetically costly processes such the flagellum biosynthesis and virulence-related functions. From the biotechnological point of view, manipulation of gene order enables altering bacterial growth with no escape events.

The Slowdown of Growth Rate Controls the Single-Cell Distribution of Biofilm Matrix Production via an SinI-SinR-SlrR Network

In Bacillus subtilis, master regulator Spo0A controls several cell-differentiation pathways. Under moderate starvation, phosphorylated Spo0A (Spo0A~P) induces biofilm formation by indirectly activating genes controlling matrix production in a subpopulation of cells via an SinI-SinR-SlrR network. Under severe starvation, Spo0A~P induces sporulation by directly and indirectly regulating sporulation gene expression. However, what determines the heterogeneity of individual cell fates is not fully understood. In particular, it is still unclear why, despite being controlled by a single master regulator, biofilm matrix production and sporulation seem mutually exclusive on a single-cell level. In this work, with mathematical modeling, we showed that the fluctuations in the growth rate and the intrinsic noise amplified by the bistability in the SinI-SinR-SlrR network could explain the single-cell distribution of matrix production. Moreover, we predicted an incoherent feed-forward loop; the decrease in the cellular growth rate first activates matrix production by increasing in Spo0A phosphorylation level but then represses it via changing the relative concentrations of SinR and SlrR. Experimental data provide evidence to support model predictions. In particular, we demonstrate how the degree to which matrix production and sporulation appear mutually exclusive is affected by genetic perturbations. IMPORTANCE The mechanisms of cell-fate decisions are fundamental to our understanding of multicellular organisms and bacterial communities. However, even for the best-studied model systems we still lack a complete picture of how phenotypic heterogeneity of genetically identical cells is controlled. Here, using B. subtilis as a model system, we employ a combination of mathematical modeling and experiments to explain the population-level dynamics and single-cell level heterogeneity of matrix gene expression. The results demonstrate how the two cell fates, biofilm matrix production and sporulation, can appear mutually exclusive without explicitly inhibiting one another. Such a mechanism could be used in a wide range of other biological systems.

Transcriptome Dynamics of Pseudomonas aeruginosa during Transition from Overlapping To Non-Overlapping Cell Cycles

Bacteria either duplicate their chromosome once per cell division or a new round of replication is initiated before the cells divide, thus cell cycles overlap. Here, we show that the opportunistic pathogen Pseudomonas aeruginosa switches from fast growth with overlapping cell cycles to sustained slow growth with only one replication round per cell division when cultivated under standard laboratory conditions. The transition was characterized by fast-paced, sequential changes in transcriptional activity along the ori-ter axis of the chromosome reflecting adaptation to the metabolic needs during both growth phases. Quorum sensing (QS) activity was highest at the onset of the slow growth phase with non-overlapping cell cycles. RNA sequencing of subpopulations of these cultures sorted based on their DNA content, revealed a strong gene dosage effect as well as specific expression patterns for replicating and nonreplicating cells. Expression of flagella and mexE, involved in multidrug efflux was restricted to cells that did not replicate, while those that did showed a high activity of the cell division locus and recombination genes. A possible role of QS in the formation of these subpopulations upon switching to non-overlapping cell cycles could be a subject of further research. IMPORTANCE The coordination of gene expression with the cell cycle has so far been studied only in a few bacteria, the bottleneck being the need for synchronized cultures. Here, we determined replication-associated effects on transcription by comparing Pseudomonas aeruginosa cultures that differ in their growth mode and number of replicating chromosomes. We further show that cell cycle-specific gene regulation can be principally identified by RNA sequencing of subpopulations from cultures that replicate only once per cell division and that are sorted according to their DNA content. Our approach opens the possibility to study asynchronously growing bacteria from a wide phylogenetic range and thereby enhance our understanding of the evolution of cell cycle control on the transcriptional level.

Carriage of three plasmids in a single human clinical isolate of Clostridioides difficile

A subset of clinical isolates of Clostridioides difficile contains one or more plasmids and these plasmids can harbor virulence and antimicrobial resistance determinants. Despite their potential importance, C. difficile plasmids remain poorly characterized. Here, we provide the complete genome sequence of a human clinical isolate that carries three high-copy number plasmids from three different plasmid families that are therefore compatible. For two of these, we identify a region capable of sustaining plasmid replication in C. difficile that is also compatible with the plasmid pCD630 that is found in many laboratory strains. Together, our data advance our understanding of C. difficile plasmid biology.

The role of replication-induced chromosomal copy numbers in spatio-temporal gene regulation and evolutionary chromosome plasticity

For a coherent response to environmental changes, bacterial evolution has formed a complex transcriptional regulatory system comprising classical DNA binding proteins sigma factors and modulation of DNA topology. In this study, we investigate replication-induced gene copy numbers - a regulatory concept that is unlike the others not based on modulation of promoter activity but on replication dynamics. We show that a large fraction of genes are predominantly affected by transient copy numbers and identify cellular functions and central pathways governed by this mechanism in Escherichia coli. Furthermore, we show quantitatively that the previously observed spatio-temporal expression pattern between different growth phases mainly emerges from transient chromosomal copy numbers. We extend the analysis to the plant pathogen Dickeya dadantii and the biotechnologically relevant organism Vibrio natriegens. The analysis reveals a connection between growth phase dependent gene expression and evolutionary gene migration in these species. A further extension to the bacterial kingdom indicates that chromosome evolution is governed by growth rate related transient copy numbers.

The modular pYT vector series employed for chromosomal gene integration and expression to produce carbazoles and glycolipids in P. putida

The expression of biosynthetic genes in bacterial hosts can enable access to high-value compounds, for which appropriate molecular genetic tools are essential. Therefore, we developed a toolbox of modular vectors, which facilitate chromosomal gene integration and expression in Pseudomonas putida KT2440. To this end, we designed an integrative sequence, allowing customisation regarding the modes of integration (random, at attTn7, or into the 16S rRNA gene), promoters, antibiotic resistance markers as well as fluorescent proteins and enzymes as transcription reporters. We thus established a toolbox of vectors carrying integrative sequences, designated as pYT series, of which we present 27 ready-to-use variants along with a set of strains equipped with unique 'landing pads' for directing a pYT interposon into one specific copy of the 16S rRNA gene. We used genes of the well-described violacein biosynthesis as reporter to showcase random Tn5-based chromosomal integration leading to constitutive expression and production of violacein and deoxyviolacein. Deoxyviolacein was likewise produced after gene integration into the 16S rRNA gene of rrn operons. Integration in the attTn7 site was used to characterise the suitability of different inducible promoters and successive strain development for the metabolically challenging production of mono-rhamnolipids. Finally, to establish arcyriaflavin A production in P. putida for the first time, we compared different integration and expression modes, revealing integration at attTn7 and expression with NagR/PnagAa to be most suitable. In summary, the new toolbox can be utilised for the rapid generation of various types of P. putida expression and production strains.

Patterns of abundance, chromosomal localization, and domain organization among c-di-GMP-metabolizing genes revealed by comparative genomics of five alphaproteobacterial orders

BACKGROUND

Bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP) is a bacterial second messenger that affects diverse processes in different bacteria, including the cell cycle, motility, and biofilm formation. Its cellular levels are controlled by the opposing activities of two types of enzymes, with synthesis by diguanylate cyclases containing a GGDEF domain and degradation by phosphodiesterases containing either an HD-GYP or an EAL domain. These enzymes are ubiquitous in bacteria with up to 50 encoded in some genomes, the specific functions of which are mostly unknown.

RESULTS

We used comparative analyses to identify genomic patterns among genes encoding proteins with GGDEF, EAL, and HD-GYP domains in five orders of the class Alphaproteobacteria. GGDEF-containing sequences and GGDEF-EAL hybrids were the most abundant and had the highest diversity of co-occurring auxiliary domains while EAL and HD-GYP containing sequences were less abundant and less diverse with respect to auxiliary domains. There were striking patterns in the chromosomal localizations of the genes found in two of the orders. The Rhodobacterales' EAL-encoding genes and Rhizobiales' GGDEF-EAL-encoding genes showed opposing patterns of distribution compared to the GGDEF-encoding genes. In the Rhodobacterales, the GGDEF-encoding genes showed a tri-modal distribution with peaks mid-way between the origin (ori) and terminus (ter) of replication and at ter while the EAL-encoding genes peaked near ori. The patterns were more complex in the Rhizobiales, but the GGDEF-encoding genes were biased for localization near ter.

CONCLUSIONS

The observed patterns in the chromosomal localizations of these genes suggest a coupling of synthesis and hydrolysis of c-di-GMP with the cell cycle. Moreover, the higher proportions and diversities of auxiliary domains associated with GGDEF domains and GGDEF-EAL hybrids compared to EAL or HD-GYP domains could indicate that more stimuli affect synthesis compared to hydrolysis of c-di-GMP.

E. coli chromosomal-driven expression of NADK2 from A. thaliana: A preferable alternative to plasmid-driven expression for challenging proteins

The expression and purification of large recombinant proteins or protein complexes is problematic for some biotechnology laboratories. Indeed, it is often difficult to obtain enough active proteins to perform biological characterization or reach commercialization, when large proteins or protein complexes are expressed in E. coli via the popular T7-based plasmid-driven expression system. There is also an industrial demand to decrease our dependence on plasmid-driven expression, because of its drawbacks, such as: i) the common use of antibiotics to maintain the plasmid, ii) the issue of plasmid copy number, and iii) the risk of overloading the expression system. Despite all these issues, alternative solutions, such as gene integration in the bacterial chromosome, are rarely employed and their advantages are still a matter of debate. Plant plastidial NAD kinases (NADK; ATP:NAD 2'-phosphotransferase, EC 2.7.1.23) are a classic example of proteins with high molecular weight, that are difficult to express and purify with traditional T7-based technology. We therefore compared plasmid-driven and chromosomal-driven expression of the Arabidopsis thaliana NADK2 protein, using a proprietary counter-selection tool, COLIBELT®, that allows scar-free and marker-free chromosomal modifications. Here we show that chromosomal-driven expression allowed recovery of more active NADK2 protein than classic T7 expression systems, as well as better production, thus confirming that expression from one single chromosomal copy is preferable to plasmid-driven expression and might be appealing for both basic and applied research.

Bilateral symmetry of linear streptomycete chromosomes

Here, we characterize an uncommon set of telomeres from Streptomyces rimosus ATCC 10970, the parental strain of a lineage of one of the earliest-discovered antibiotic producers. Following the closure of its genome sequence, we compared unusual telomeres from this organism with the other five classes of replicon ends found amongst streptomycetes. Closed replicons of streptomycete chromosomes were organized with respect to their phylogeny and physical orientation, which demonstrated that different telomeres were not associated with particular clades and are likely shared amongst different strains by plasmid-driven horizontal gene transfer. Furthermore, we identified a ~50 kb origin island with conserved synteny that is located at the core of all streptomycete chromosomes and forms an axis around which symmetrical chromosome inversions can take place. Despite this chromosomal bilateral symmetry, a bias in parS sites to the right of oriC is maintained across the family Streptomycetaceae and suggests that the formation of ParB/parS nucleoprotein complexes on the right replichore is a conserved feature in streptomycetes. Consequently, our studies reveal novel features of linear bacterial replicons that, through their manipulation, may lead to improvements in growth and productivity of this important industrial group of bacteria.

The Pseudoalteromonas multipartite genome: distribution and expression of pangene categories, and a hypothesis for the origin and evolution of the chromid

Bacterial genomes typically consist of one large chromosome, but can also include secondary replicons. These so-called multipartite genomes are scattered on the bacterial tree of life with the majority of cases belonging to Proteobacteria. Within the class gamma-proteobacteria, multipartite genomes are restricted to the two families Vibrionaceae and Pseudoalteromonadaceae. Whereas the genome of vibrios is well studied, information on the Pseudoalteromonadaceae genome is much scarcer. We have studied Pseudoalteromonadaceae with respect to the origin of the chromid, how pangene categories are distributed, how genes are expressed relative to their genomic location, and identified chromid hallmark genes. We calculated the Pseudoalteromonadaceae pangenome based on 25 complete genomes and found that core/softcore are significantly overrepresented in late replicating sectors of the chromid, regardless of how the chromid is replicated. On the chromosome, core/softcore and shell/cloud genes are only weakly overrepresented at the chromosomal replication origin and termination sequences, respectively. Gene expression is trending downwards with increasing distance from the chromosomal oriC, whereas the chromidal expression pattern is more complex. Moreover, we identified 78 chromid hallmark genes, and BLASTp searches suggest that the majority of them were acquired from the ancestral gene pool of Alteromonadales. Finally, our data strongly suggest that the chromid originates from a plasmid that was acquired in a relatively recent event. In summary, this study extends our knowledge on multipartite genomes, and helps us understand how and why secondary replicons are acquired, why they are maintained, and how they are shaped by evolution.

The pIT5 Plasmid Series, an Improved Toolkit for Repeated Genome Integration in E. coli

We describe a new set of tools for inserting DNA into the bacterial chromosome. The system uses site-specific recombination reactions carried out by bacteriophage integrases to integrate plasmids at up to eight phage attachment sites in E. coli MG1655. The introduction of mutant loxP sites in the integrating plasmids allows repeated removal of antibiotic resistance genes and other plasmid sequences without danger of inducing chromosomal rearrangements. The protocol for Cre-mediated antibiotic resistance gene removal is greatly simplified by introducing the Cre plasmid by phage infection. Finally, we have also developed a set of four independently inducible expression modules with tight control and high dynamic range which can be inserted at specific chromosomal locations.

Connection Between Chromosomal Location and Function of CtrA Phosphorelay Genes in Alphaproteobacteria

Most bacterial chromosomes are circular, with replication starting at one origin (ori) and proceeding on both replichores toward the terminus (ter). Several studies have shown that the location of genes relative to ori and ter can have profound effects on regulatory networks and physiological processes. The CtrA phosphorelay is a gene regulatory system conserved in most alphaproteobacteria. It was first discovered in Caulobacter crescentus where it controls replication and division into a stalked and a motile cell in coordination with other factors. The locations of the ctrA gene and targets of this response regulator on the chromosome affect their expression through replication-induced DNA hemi-methylation and specific positioning along a CtrA activity gradient in the dividing cell, respectively. Here we asked to what extent the location of CtrA regulatory network genes might be conserved in the alphaproteobacteria. We determined the locations of the CtrA phosphorelay and associated genes in closed genomes with unambiguously identifiable ori from members of five alphaproteobacterial orders. The location of the phosphorelay genes was the least conserved in the Rhodospirillales followed by the Sphingomonadales. In the Rhizobiales a trend toward certain chromosomal positions could be observed. Compared to the other orders, the CtrA phosphorelay genes were conserved closer to ori in the Caulobacterales. In contrast, the genes were highly conserved closer to ter in the Rhodobacterales. Our data suggest selection pressure results in differential positioning of CtrA phosphorelay and associated genes in alphaproteobacteria, particularly in the orders Rhodobacterales, Caulobacterales and Rhizobiales that is worth deeper investigation.

Pangenome analyses of LuxS-coding genes and enzymatic repertoires in cocoa-related lactic acid bacteria

Lactobacillaceae presents potential for interspecific Quorum Sensing (QS) in spontaneous cocoa fermentation, correlated with high abundance of luxS. Three Brazilian isolates from cocoa fermentation were characterized by Whole Genome Sequencing and luxS gene was surveyed in their genomes, in comparison with public databases. They were classified as Lactiplantibacillus plantarum, Limosilactobacillus fermentum and Pediococcus acidilactici. LuxS genes were conserved in core genomes of the novel isolates, but in some non-cocoa related Lactic Acid Bacteria (LAB) it was accessory and plasmid-borne. The conservation and horizontal acquisition of luxS reinforces that QS is determinant for bacterial adaptation in several environments, especially taking into account the luxS has been correlated with modulation of bacteriocin production, stress tolerance and biofilm formation. Therefore, in this paper, new clade and species-specific primers were designed for future application for screening of luxS gene in LAB to evaluate the adaptive potential to diverse food fermentations.

Transcription Factor Binding to Replicated DNA

Genome replication perturbs the DNA regulatory environment by displacing DNA-bound proteins, replacing nucleosomes, and introducing dosage imbalance between regions replicating at different S-phase stages. Recently, we showed that these effects are integrated to maintain transcription homeostasis: replicated genes increase in dosage, but their expression remains stable due to replication-dependent epigenetic changes that suppress transcription. Here, we examine whether reduced transcription from replicated DNA results from limited accessibility to regulatory factors by measuring the time-resolved binding of RNA polymerase II (Pol II) and specific transcription factors (TFs) to DNA during S phase in budding yeast. We show that the Pol II binding pattern is largely insensitive to DNA dosage, indicating limited binding to replicated DNA. In contrast, binding of three TFs (Reb1, Abf1, and Rap1) to DNA increases with the increasing DNA dosage. We conclude that the replication-specific chromatin environment remains accessible to regulatory factors but suppresses RNA polymerase recruitment.

Virulence and antimicrobial resistance genes are enriched in the plasmidome of clinical Escherichia coli isolates compared with wastewater isolates from western Kenya

Many low-middle income countries in Africa have poorly-developed infectious disease monitoring systems. Here, we employed whole genome sequencing (WGS) to investigate the presence/absence of antimicrobial resistance (AMR) and virulence-associated (VA) genes in a collection of clinical and municipal wastewater Escherichia coli isolates from Kakamega, west Kenya. We were particularly interested to see whether, given the association between infection and water quality, the isolates from these geographically-linked environments might display similar genomic signatures. Phylogenetic analysis based on the core genes common to all of the isolates revealed two broad divisions, corresponding to the commensal/enterotoxigenic E. coli on the one hand, and uropathogenic E. coli on the other. Although the clinical and wastewater isolates each contained a very similar mean number of antibiotic resistance-encoding genes, the clinical isolates were enriched in genes required for in-host survival. Furthermore, and although the chromosomally encoded repertoire of these genes was similar in all sequenced isolates, the genetic composition of the plasmids from clinical and wastewater E. coli was more habitat-specific, with the clinical isolate plasmidome enriched in AMR and VA genes. Intriguingly, the plasmid-borne VA genes were often duplicates of genes already present on the chromosome, whereas the plasmid-borne AMR determinants were more specific. This reinforces the notion that plasmids are a primary means by which infection-related AMR and VA-associated genes are acquired and disseminated among these strains.

Tos4 mediates gene expression homeostasis through interaction with HDAC complexes independently of H3K56 acetylation

Saccharomyces cerevisiae exhibits gene expression homeostasis, which is defined as the buffering of transcription levels against changes in DNA copy number during the S phase of the cell cycle. It has been suggested that S. cerevisiae employs an active mechanism to maintain gene expression homeostasis through Rtt109-Asf1-dependent acetylation of histone H3 on lysine 56 (H3K56). Here, we show that gene expression homeostasis can be achieved independently of H3K56 acetylation by Tos4 (Target of Swi6-4). Using Nanostring technology, we establish that Tos4-dependent gene expression homeostasis depends on its forkhead-associated (FHA) domain, which is a phosphopeptide recognition domain required to bind histone deacetylases (HDACs). We demonstrate that the mechanism of Tos4-dependent gene expression homeostasis requires its interaction with the Rpd3L HDAC complex. However, this is independent of Rpd3's well-established roles in both histone deacetylation and controlling the DNA replication timing program, as established by deep sequencing of Fluorescence-Activated Cell Sorted (FACS) S and G2 phase populations. Overall, our data reveals that Tos4 mediates gene expression homeostasis through its FHA domain-dependent interaction with the Rpd3L complex, which is independent of H3K56ac.

Genome Complexity Browser: Visualization and quantification of genome variability

Comparative genomics studies may be used to acquire new knowledge regarding genome architecture, which defines the rules for combining sets of genes in the genome of living organisms. Hundreds of thousands of prokaryotic genomes have been sequenced and assembled. However, computational tools capable of simultaneously comparing large numbers of genomes are lacking. We developed the Genome Complexity Browser, a tool that allows the visualization of gene contexts, in a graph-based format, and the quantification of variability for different segments of a genome. The graph-based visualization allows the inspection of changes in gene contents and neighborhoods across hundreds of genomes, simultaneously, which may facilitate the identification of conserved and variable segments of operons or the estimation of the overall variability associated with a particular genome locus. We introduced a measure called complexity, to quantify genome variability. Intraspecies and interspecies comparisons revealed that regions with high complexity values tended to be located in areas that are conserved across different strains and species.

The comparison of genomes among different bacteria and archaea species has revealed that many species frequently exchange genes. Occasionally, such horizontal gene transfer events result in the acquisition of pathogenic properties or antibiotic resistance in the recipient organism. Previously, the probabilities of gene insertions were found to vary, with unequal distributions along a chromosome. At some loci, referred to as hotspots, changes occur with much higher frequencies compared with other regions of the chromosome. We developed a computational method and a software tool, called Genome Complexity Browser, that allows the identification of genome variability hotspots and the visualization of changes. We compared the localization of various hotspots and revealed that some demonstrate conserved localizations, even across species, whereas others are transient. Our tool allows users to visually inspect the patterns of gene changes in graph-based format, which presents the visualization in a format that is both compact and informative.

Vibrionaceae core, shell and cloud genes are non-randomly distributed on Chr 1: An hypothesis that links the genomic location of genes with their intracellular placement

Background

The genome of Vibrionaceae bacteria, which consists of two circular chromosomes, is replicated in a highly ordered fashion. In fast-growing bacteria, multifork replication results in higher gene copy numbers and increased expression of genes located close to the origin of replication of Chr 1 (ori1). This is believed to be a growth optimization strategy to satisfy the high demand of essential growth factors during fast growth. The relationship between ori1-proximate growth-related genes and gene expression during fast growth has been investigated by many researchers. However, it remains unclear which other gene categories that are present close to ori1 and if expression of all ori1-proximate genes is increased during fast growth, or if expression is selectively elevated for certain gene categories.

Results

We calculated the pangenome of all complete genomes from the Vibrionaceae family and mapped the four pangene categories, core, softcore, shell and cloud, to their chromosomal positions. This revealed that core and softcore genes were found heavily biased towards ori1, while shell genes were overrepresented at the opposite part of Chr 1 (i.e., close to ter1). RNA-seq of Aliivibrio salmonicida and Vibrio natriegens showed global gene expression patterns that consistently correlated with chromosomal distance to ori1. Despite a biased gene distribution pattern, all pangene categories contributed to a skewed expression pattern at fast-growing conditions, whereas at slow-growing conditions, softcore, shell and cloud genes were responsible for elevated expression.

Conclusion

The pangene categories were non-randomly organized on Chr 1, with an overrepresentation of core and softcore genes around ori1, and overrepresentation of shell and cloud genes around ter1. Furthermore, we mapped our gene distribution data on to the intracellular positioning of chromatin described for V. cholerae, and found that core/softcore and shell/cloud genes appear enriched at two spatially separated intracellular regions. Based on these observations, we hypothesize that there is a link between the genomic location of genes and their cellular placement.

Correlated chromosomal periodicities according to the growth rate and gene expression

Linking genetic information to population fitness is crucial to understanding living organisms. Despite the abundant knowledge of the genetic contribution to growth, the overall patterns/features connecting genes, their expression, and growth remain unclear. To reveal the quantitative and direct connections, systematic growth assays of single-gene knockout Escherichia coli strains under both rich and poor nutritional conditions were performed; subsequently, the resultant growth rates were associated with the original expression levels of the knockout genes in the parental genome. Comparative analysis of growth and the transcriptome identified not only the nutritionally differentiated fitness cost genes but also a significant correlation between the growth rates of the single-gene knockout strains and the original expression levels of these knockout genes in the parental strain, regardless of the nutritional variation. In addition, the coordinated chromosomal periodicities of the wild-type transcriptome and the growth rates of the strains lacking the corresponding genes were observed. The common six-period periodicity was somehow attributed to the essential genes, although the underlying mechanism remains to be addressed. The correlated chromosomal periodicities associated with the gene expression-growth dataset were highly valuable for bacterial growth prediction and discovering the working principles governing minimal genetic information.

Succession of Bifidobacterium longum Strains in Response to a Changing Early Life Nutritional Environment Reveals Dietary Substrate Adaptations

Diet-microbe interactions play a crucial role in modulation of the early life microbiota and infant health. Bifidobacterium dominates the breast-fed infant gut and may persist in individuals during transition from a milk-based to a more diversified diet. Here, we investigated adaptation of Bifidobacterium longum to the changing nutritional environment. Genomic characterization of 75 strains isolated from nine either exclusively breast- or formula-fed (pre-weaning) infants in their first 18 months revealed subspecies- and strain-specific intra-individual genomic diversity with respect to carbohydrate metabolism, which corresponded to different dietary stages. Complementary phenotypic studies indicated strain-specific differences in utilization of human milk oligosaccharides and plant carbohydrates, whereas proteomic profiling identified gene clusters involved in metabolism of selected carbohydrates. Our results indicate a strong link between infant diet and B. longum diversity and provide additional insights into possible competitive advantage mechanisms of this Bifidobacterium species and its persistence in a single host.

Precision design of stable genetic circuits carried in highly-insulated E. coli genomic landing pads

Genetic circuits have many applications, from guiding living therapeutics to ordering process in a bioreactor, but to be useful they have to be genetically stable and not hinder the host. Encoding circuits in the genome reduces burden, but this decreases performance and can interfere with native transcription. We have designed genomic landing pads in Escherichia coli at high-expression sites, flanked by ultrastrong double terminators. DNA payloads >8 kb are targeted to the landing pads using phage integrases. One landing pad is dedicated to carrying a sensor array, and two are used to carry genetic circuits. NOT/NOR gates based on repressors are optimized for the genome and characterized in the landing pads. These data are used, in conjunction with design automation software (Cello 2.0), to design circuits that perform quantitatively as predicted. These circuits require fourfold less RNA polymerase than when carried on a plasmid and are stable for weeks in a recA+ strain without selection. This approach enables the design of synthetic regulatory networks to guide cells in environments or for applications where plasmid use is infeasible.

Evaluation of chromosomal insertion loci in the Pseudomonas putida KT2440 genome for predictable biosystems design

The development of Pseudomonas strains for industrial production of fuels and chemicals will require the integration of heterologous genes and pathways into the chromosome. Finding the most appropriate integration site to maximize strain performance is an essential part of the strain design process. We characterized seven chromosomal loci in Pseudomonas putida KT2440 for integration of a fluorescent protein expression construct. Insertion in five of the loci did not affect growth rate, but fluorescence varied by up to 27-fold. Three sites displaying a diversity of phenotypes with the fluorescent reporter were also chosen for the integration of a gene encoding a muconate importer. Depending on the integration locus, expression of the importer varied by approximately 3-fold and produced significant phenotypic differences. This work demonstrates the impact of the integration location on host viability, gene expression, and overall strain performance.

•

Pseudomonas putida KT2440 chromosomal loci were characterized as potential insertion targets for heterologous genes.

•

Integration location had a significant effect on heterologous protein expression and host phenotype.

•

The identification of an appropriate chromosomal insertion location is essential to optimize genetic engineering design.

Macromolecular crowding links ribosomal protein gene dosage to growth rate in Vibrio cholerae

BACKGROUND

In fast-growing bacteria, the genomic location of ribosomal protein (RP) genes is biased towards the replication origin (oriC). This trait allows optimizing their expression during exponential phase since oriC neighboring regions are in higher dose due to multifork replication. Relocation of s10-spc-α locus (S10), which codes for most of the RP, to ectopic genomic positions shows that its relative distance to the oriC correlates to a reduction on its dosage, its expression, and bacterial growth rate. However, a mechanism linking S10 dosage to cell physiology has still not been determined.

RESULTS

We hypothesized that S10 dosage perturbations impact protein synthesis capacity. Strikingly, we observed that in Vibrio cholerae, protein production capacity was independent of S10 position. Deep sequencing revealed that S10 relocation altered chromosomal replication dynamics and genome-wide transcription. Such changes increased as a function of oriC-S10 distance. Since RP constitutes a large proportion of cell mass, lower S10 dosage could lead to changes in macromolecular crowding, impacting cell physiology. Accordingly, cytoplasm fluidity was higher in mutants where S10 is most distant from oriC. In hyperosmotic conditions, when crowding differences are minimized, the growth rate and replication dynamics were highly alleviated in these strains.

CONCLUSIONS

The genomic location of RP genes ensures its optimal dosage. However, besides of its essential function in translation, their genomic position sustains an optimal macromolecular crowding essential for maximizing growth. Hence, this could be another mechanism coordinating DNA replication to bacterial growth.

Topology-dependent interference of synthetic gene circuit function by growth feedback

Growth-mediated feedback between synthetic gene circuits and host organisms leads to diverse emerged behaviors, including growth bistability and enhanced ultrasensitivity. However, the range of possible impacts of growth feedback on gene circuits remains underexplored. Here, we mathematically and experimentally demonstrated that growth feedback affects the functions of memory circuits in a network topology-dependent way. Specifically, the memory of the self-activation switch is quickly lost due to the growth-mediated dilution of the circuit products. Decoupling of growth feedback reveals its memory, manifested by its hysteresis property across a broad range of inducer concentration. On the contrary, the toggle switch is more refractory to growth-mediated dilution and can retrieve its memory after the fast-growth phase. The underlying principle lies in the different dependence of active and repressive regulations in these circuits on the growth-mediated dilution. Our results unveil the topology-dependent mechanism on how growth-mediated feedback influences the behaviors of gene circuits.

The Impact of Global Transcriptional Regulation on Bacterial Gene Order

Bacterial gene expression depends on the allocation of limited transcriptional resources provided a particular growth rate and growth condition. Early studies in a few genes suggested this global regulation to generate a unifying hyperbolic expression pattern. Here, we developed a large-scale method that generalizes these experiments to quantify the response to growth of over 700 genes that a priori do not exhibit any specific control. We distinguish a core subset following a promoter-specific hyperbolic response. Within this group, we sort genes with regard to their responsiveness to the global regulatory program to show that those with a particularly sensitive linear response are located near the origin of replication. We then find evidence that this genomic architecture is biologically significant by examining position conservation of E. coli genes in 100 bacteria. The response to the transcriptional resources of the cell results in an additional feature contributing to bacterial genome organization.

•

Cell physiology determines a global transcriptional regulatory program

•

Constitutive genes show a differential response to this global regulation

•

The most responsive constitutive genes are located near the origin of replication

•

Global transcriptional regulation acts as a gene position conservation force

Conserved Patterns of Symmetric Inversion in the Genome Evolution of Bordetella Respiratory Pathogens

Whooping cough (pertussis), primarily caused by Bordetella pertussis, has resurged in the United States, and circulating strains exhibit considerable chromosome structural fluidity in the form of rearrangement and deletion. The genus Bordetella includes additional pathogenic species infecting various animals, some even causing pertussis-like respiratory disease in humans; however, investigation of their genome evolution has been limited. We studied chromosome structure in complete genome sequences from 167 Bordetella species isolates, as well as 469 B. pertussis isolates, to gain a generalized understanding of rearrangement patterns among these related pathogens. Observed changes in gene order primarily resulted from large inversions and were only detected in species with genomes harboring multicopy insertion sequence (IS) elements, most notably B. holmesii and B. parapertussis While genomes of B. pertussis contain >240 copies of IS481, IS elements appear less numerous in other species and yield less chromosome structural diversity through rearrangement. These data were further used to predict all possible rearrangements between IS element copies present in Bordetella genomes, revealing that only a subset is observed among circulating strains. Therefore, while it appears that rearrangement occurs less frequently in other species than in B. pertussis, these clinically relevant respiratory pathogens likely experience similar mutation of gene order. The resulting chromosome structural fluidity presents both challenges and opportunity for the study of Bordetella respiratory pathogens.IMPORTANCE Bordetella pertussis is the primary agent of whooping cough (pertussis). The Bordetella genus includes additional pathogens of animals and humans, including some that cause pertussis-like respiratory illness. The chromosome of B. pertussis has previously been shown to exhibit considerable structural rearrangement, but insufficient data have prevented comparable investigation in related species. In this study, we analyze chromosome structure variation in several Bordetella species to gain a generalized understanding of rearrangement patterns in this genus. Just as in B. pertussis, we observed inversions in other species that likely result from common mutational processes. We used these data to further predict additional, unobserved inversions, suggesting that specific genome structures may be preferred in each species.

EVR: reconstruction of bacterial chromosome 3D structure models using error-vector resultant algorithm

BACKGROUND

More and more 3C/Hi-C experiments on prokaryotes have been published. However, most of the published modeling tools for chromosome 3D structures are targeting at eukaryotes. How to transform prokaryotic experimental chromosome interaction data into spatial structure models is an important task and in great need.

RESULTS

We have developed a new reconstruction program for bacterial chromosome 3D structure models called EVR that exploits a simple Error-Vector Resultant (EVR) algorithm. This software tool is particularly optimized for the closed-loop structural features of prokaryotic chromosomes. The parallel implementation of the program can utilize the computing power of both multi-core CPUs and GPUs.

CONCLUSIONS

EVR can be used to reconstruct the bacterial 3D chromosome structure based on the contact frequency matrix derived from 3C/Hi-C experimental data quickly and precisely.

A Quantitative Model Explains Single-Cell Dynamics of the Adaptive Response in Escherichia coli

DNA damage caused by alkylating chemicals induces an adaptive response in Escherichia coli that increases the tolerance of cells to further damage. Signaling of the response occurs through irreversible methylation of the Ada protein, which acts as a DNA repair protein and damage sensor. Methylated Ada induces its own gene expression through a positive feedback loop. However, random fluctuations in the abundance of Ada jeopardize the reliability of the induction signal. I developed a quantitative model to test how gene expression noise and feedback amplification affect the fidelity of the adaptive response. A remarkably simple model accurately reproduced experimental observations from single-cell measurements of gene expression dynamics in a microfluidic device. Stochastic simulations showed that delays in the adaptive response are a direct consequence of the very low number of Ada molecules present to signal DNA damage. For cells that have zero copies of Ada, response activation becomes a memoryless process that is dictated by an exponential waiting time distribution between basal Ada expression events. Experiments also confirmed the model prediction that the strength of the adaptive response drops with an increasing growth rate of cells.

High-resolution analysis of the pneumococcal transcriptome under a wide range of infection-relevant conditions

Streptococcus pneumoniae is an opportunistic human pathogen that typically colonizes the nasopharyngeal passage and causes lethal disease in other host niches, such as the lung or the meninges. The expression and regulation of pneumococcal genes at different life-cycle stages, such as commensal or pathogenic, are not entirely understood. To chart the transcriptional responses of S. pneumoniae, we used RNA-seq to quantify the relative abundance of the transcriptome under 22 different infection-relevant conditions. The data demonstrated a high level of dynamic expression and, strikingly, all annotated pneumococcal genomic features were expressed in at least one of the studied conditions. By computing the correlation values of every pair of genes across all studied conditions, we created a co-expression matrix that provides valuable information on both operon structure and regulatory processes. The co-expression data are highly consistent with well-characterized operons and regulons, such as the PyrR, ComE and ComX regulons, and have allowed us to identify a new member of the competence regulon. Lastly, we created an interactive data center named PneumoExpress (https://veeninglab.com/pneumoexpress) that enables users to access the expression data as well as the co-expression matrix in an intuitive and efficient manner, providing a valuable resource to the pneumococcal research community.

Pseudomonas putida rDNA is a favored site for the expression of biosynthetic genes

Since high-value bacterial secondary metabolites, including antibiotics, are often naturally produced in only low amounts, their efficient biosynthesis typically requires the transfer of entire metabolic pathways into suitable bacterial hosts like Pseudomonas putida. Stable maintenance and sufficient expression of heterologous pathway-encoding genes in host microbes, however, still remain key challenges. In this study, the 21 kb prodigiosin gene cluster from Serratia marcescens was used as a reporter to identify genomic sites in P. putida KT2440 especially suitable for maintenance and expression of pathway genes. After generation of a strain library by random Tn5 transposon-based chromosomal integration of the cluster, 50 strains exhibited strong prodigiosin production. Remarkably, chromosomal integration sites were exclusively identified in the seven rRNA-encoding rrn operons of P. putida. We could further demonstrate that prodigiosin production was mainly dependent on (i) the individual rrn operon where the gene cluster was inserted as well as (ii) the distance between the rrn promoter and the inserted prodigiosin biosynthetic genes. In addition, the recombinant strains showed high stability upon subculturing for many generations. Consequently, our findings demonstrate the general applicability of rDNA loci as chromosomal integration sites for gene cluster expression and recombinant pathway implementation in P. putida KT2440.

Synthetic biology approaches for chromosomal integration of genes and pathways in industrial microbial systems

Industrial biotechnology is reliant on native pathway engineering or foreign pathway introduction for efficient biosynthesis of target products. Chromosomal integration, with intrinsic genetic stability, is an indispensable step for reliable expression of homologous or heterologous genes and pathways in large-scale and long-term fermentation. With advances in synthetic biology and CRISPR-based genome editing approaches, a wide variety of novel enabling technologies have been developed for single-step, markerless, multi-locus genomic integration of large biochemical pathways, which significantly facilitate microbial overproduction of chemicals, pharmaceuticals and other value-added biomolecules. Notably, the newly discovered homology-mediated end joining strategy could be widely applicable for high-efficiency genomic integration in a number of homologous recombination-deficient microbes. In this review, we explore the fundamental principles and characteristics of genomic integration, and highlight the development and applications of targeted integration approaches in the three representative industrial microbial systems, including Escherichia coli, actinomycetes and yeasts.

Genome Location Dictates the Transcriptional Response to PolC Inhibition in Clostridium difficile

Clostridium difficile is a potentially lethal gut pathogen that causes nosocomial and community-acquired infections. Limited treatment options and reports of reduced susceptibility to current treatment emphasize the necessity for novel antimicrobials. The DNA polymerase of Gram-positive organisms is an attractive target for the development of antimicrobials. ACX-362E [N ²-(3,4-dichlorobenzyl)-7-(2-[1-morpholinyl]ethyl)guanine; MorE-DCBG] is a DNA polymerase inhibitor in preclinical development as a novel therapeutic against C. difficile infection. This synthetic purine shows preferential activity against C. difficile PolC over those of other organisms in vitro and is effective in an animal model of C. difficile infection. In this study, we have determined its efficacy against a large collection of clinical isolates. At concentrations below the MIC, the presumed slowing (or stalling) of replication forks due to ACX-362E leads to a growth defect. We have determined the transcriptional response of C. difficile to replication inhibition and observed an overrepresentation of upregulated genes near the origin of replication in the presence of PolC inhibitors, but not when cells were subjected to subinhibitory concentrations of other antibiotics. This phenomenon can be explained by a gene dosage shift, as we observed a concomitant increase in the ratio between origin-proximal and terminus-proximal gene copy number upon exposure to PolC inhibitors. Moreover, we show that certain genes differentially regulated under PolC inhibition are controlled by the origin-proximal general stress response regulator sigma factor B. Together, these data suggest that genome location both directly and indirectly determines the transcriptional response to replication inhibition in C. difficile.

Systematic Dissection of Sequence Elements Controlling σ70 Promoters Using a Genomically Encoded Multiplexed Reporter Assay in Escherichia coli

Promoters are the key drivers of gene expression and are largely responsible for the regulation of cellular responses to time and environment. In Escherichia coli, decades of studies have revealed most, if not all, of the sequence elements necessary to encode promoter function. Despite our knowledge of these motifs, it is still not possible to predict the strength and regulation of a promoter from primary sequence alone. Here we develop a novel multiplexed assay to study promoter function in E. coli by building a site-specific genomic recombination-mediated cassette exchange system that allows for the facile construction and testing of large libraries of genetic designs integrated into precise genomic locations. We build and test a library of 10898 σ70 promoter variants consisting of all combinations of a set of eight -35 elements, eight -10 elements, three UP elements, eight spacers, and eight backgrounds. We find that the -35 and -10 sequence elements can explain approximately 74% of the variance in promoter strength within our data set using a simple log-linear statistical model. Simple neural network models explain >95% of the variance in our data set by capturing nonlinear interactions with the spacer, background, and UP elements.

YocM a small heat shock protein can protect Bacillus subtilis cells during salt stress

Small heat shock proteins (sHsp) occur in all domains of life. By interacting with misfolded or aggregated proteins these chaperones fulfill a protective role in cellular protein homeostasis. Here, we demonstrate that the sHsp YocM of the Gram-positive model organism Bacillus subtilis is part of the cellular protein quality control system with a specific role in salt stress response. In the absence of YocM the survival of salt shocked cells is impaired, and increased levels of YocM protect B. subtilis exposed to heat or salt. We observed a salt and heat stress-induced localization of YocM to intracellular protein aggregates. Interestingly, purified YocM appears to accelerate protein aggregation of different model substrates in vitro. In addition, the combined presence of YocM and chemical chaperones, which accumulate in salt stressed cells, can facilitate in vitro a synergistic protective effect on protein misfolding. Therefore, the beneficial role of YocM during salt stress could be related to a mutual functional relationship with chemical chaperones and adds a new possible functional aspect to sHsp chaperone activities.

Chromosome 1 licenses chromosome 2 replication in Vibrio cholerae by doubling the crtS gene dosage

Initiation of chromosome replication in bacteria is precisely timed in the cell cycle. Bacteria that harbor multiple chromosomes face the additional challenge of orchestrating replication initiation of different chromosomes. In Vibrio cholerae, the smaller of its two chromosomes, Chr2, initiates replication after Chr1 such that both chromosomes terminate replication synchronously. The delay is due to the dependence of Chr2 initiation on the replication of a site, crtS, on Chr1. The mechanism by which replication of crtS allows Chr2 replication remains unclear. Here, we show that blocking Chr1 replication indeed blocks Chr2 replication, but providing an extra crtS copy in replication-blocked Chr1 permitted Chr2 replication. This demonstrates that unreplicated crtS copies have significant activity, and suggests that a role of replication is to double the copy number of the site that sufficiently increases its activity for licensing Chr2 replication. We further show that crtS activity promotes the Chr2-specific initiator function and that this activity is required in every cell cycle, as would be expected of a cell-cycle regulator. This study reveals how increase of gene dosage through replication can be utilized in a critical regulatory switch.

Dynamics of Bacterial Gene Regulatory Networks

The ability of bacterial cells to adjust their gene expression program in response to environmental perturbation is often critical for their survival. Recent experimental advances allowing us to quantitatively record gene expression dynamics in single cells and in populations coupled with mathematical modeling enable mechanistic understanding on how these responses are shaped by the underlying regulatory networks. Here, we review how the combination of local and global factors affect dynamical responses of gene regulatory networks. Our goal is to discuss the general principles that allow extrapolation from a few model bacteria to less understood microbes. We emphasize that, in addition to well-studied effects of network architecture, network dynamics are shaped by global pleiotropic effects and cell physiology.

Stable Regulation of Cell Cycle Events in Mycobacteria: Insights From Inherently Heterogeneous Bacterial Populations

Model bacteria, such as E. coli and B. subtilis, tightly regulate cell cycle progression to achieve consistent cell size distributions and replication dynamics. Many of the hallmark features of these model bacteria, including lateral cell wall elongation and symmetric growth and division, do not occur in mycobacteria. Instead, mycobacterial growth is characterized by asymmetric polar growth and division. This innate asymmetry creates unequal birth sizes and growth rates for daughter cells with each division, generating a phenotypically heterogeneous population. Although the asymmetric growth patterns of mycobacteria lead to a larger variation in birth size than typically seen in model bacterial populations, the cell size distribution is stable over time. Here, we review the cellular mechanisms of growth, division, and cell cycle progression in mycobacteria in the face of asymmetry and inherent heterogeneity. These processes coalesce to control cell size. Although Mycobacterium smegmatis and Mycobacterium bovis Bacillus Calmette-Guérin (BCG) utilize a novel model of cell size control, they are similar to previously studied bacteria in that initiation of DNA replication is a key checkpoint for cell division. We compare the regulation of DNA replication initiation and strategies used for cell size homeostasis in mycobacteria and model bacteria. Finally, we review the importance of cellular organization and chromosome segregation relating to the physiology of mycobacteria and consider how new frameworks could be applied across the wide spectrum of bacterial diversity.

The Transcriptome of Streptococcus pneumoniae Induced by Local and Global Changes in Supercoiling

The bacterial chromosome is compacted in a manner optimal for DNA transactions to occur. The degree of compaction results from the level of DNA-supercoiling and the presence of nucleoid-binding proteins. DNA-supercoiling is homeostatically maintained by the opposing activities of relaxing DNA topoisomerases and negative supercoil-inducing DNA gyrase. DNA-supercoiling acts as a general cis regulator of transcription, which can be superimposed upon other types of more specific trans regulatory mechanism. Transcriptomic studies on the human pathogen Streptococcus pneumoniae, which has a relatively small genome (∼2 Mb) and few nucleoid-binding proteins, have been performed under conditions of local and global changes in supercoiling. The response to local changes induced by fluoroquinolone antibiotics, which target DNA gyrase subunit A and/or topoisomerase IV, involves an increase in oxygen radicals which reduces cell viability, while the induction of global supercoiling changes by novobiocin (a DNA gyrase subunit B inhibitor), or by seconeolitsine (a topoisomerase I inhibitor), has revealed the existence of topological domains that specifically respond to such changes. The control of DNA-supercoiling in S. pneumoniae occurs mainly via the regulation of topoisomerase gene transcription: relaxation triggers the up-regulation of gyrase and the down-regulation of topoisomerases I and IV, while hypernegative supercoiling down-regulates the expression of topoisomerase I. Relaxation affects 13% of the genome, with the majority of the genes affected located in 15 domains. Hypernegative supercoiling affects 10% of the genome, with one quarter of the genes affected located in 12 domains. However, all the above domains overlap, suggesting that the chromosome is organized into topological domains with fixed locations. Based on its response to relaxation, the pneumococcal chromosome can be said to be organized into five types of domain: up-regulated, down-regulated, position-conserved non-regulated, position-variable non-regulated, and AT-rich. The AT content is higher in the up-regulated than in the down-regulated domains. Genes within the different domains share structural and functional characteristics. It would seem that a topology-driven selection pressure has defined the chromosomal location of the metabolism, virulence and competence genes, which suggests the existence of topological rules that aim to improve bacterial fitness.

Transcriptome analysis of Haloquadratum walsbyi: vanity is but the surface

Background

Haloquadratum walsbyi dominates saturated thalassic lakes worldwide where they can constitute up to 80-90% of the total prokaryotic community. Despite the abundance of the enigmatic square-flattened cells, only 7 isolates are currently known with 2 genomes fully sequenced and annotated due to difficulties to grow them under laboratory conditions. We have performed a transcriptomic analysis of one of these isolates, the Spanish strain HBSQ001 in order to investigate gene transcription under light and dark conditions.

Results

Despite a potential advantage for light as additional source of energy, no significant differences were found between light and dark expressed genes. Constitutive high gene expression was observed in genes encoding surface glycoproteins, light mediated proton pumping by bacteriorhodopsin, several nutrient uptake systems, buoyancy and storage of excess carbon. Two low expressed regions of the genome were characterized by a lower codon adaptation index, low GC content and high incidence of hypothetical genes.

Conclusions

Under the extant cultivation conditions, the square hyperhalophile devoted most of its transcriptome towards processes maintaining cell integrity and exploiting solar energy. Surface glycoproteins are essential for maintaining the large surface to volume ratio that facilitates light and organic nutrient harvesting whereas constitutive expression of bacteriorhodopsin warrants an immediate source of energy when light becomes available.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3892-2) contains supplementary material, which is available to authorized users.

Fluorescent reporters for markerless genomic integration in Staphylococcus aureus

We present integration vectors for Staphylococcus aureus encoding the fluorescent reporters mAmetrine, CFP, sGFP, YFP, mCherry and mKate. The expression is driven either from the sarA-P1 promoter or from any other promoter of choice. The reporter can be inserted markerless in the chromosome of a wide range of S. aureus strains. The integration site chosen does not disrupt any open reading frame, provides good expression, and has no detectable effect on the strains physiology. As an intermediate construct, we present a set of replicating plasmids containing the same fluorescent reporters. Also in these reporter plasmids the sarA-P1 promoter can be replaced by any other promoter of interest for expression studies. Cassettes from the replication plasmids can be readily swapped with the integration vector. With these constructs it becomes possible to monitor reporters of separate fluorescent wavelengths simultaneously.

The Proximity of Ribosomal Protein Genes to oriC Enhances Vibrio cholerae Fitness in the Absence of Multifork Replication

Recent works suggest that bacterial gene order links chromosome structure to cell homeostasis. Comparative genomics showed that, in fast-growing bacteria, ribosomal protein genes (RP) locate near the replication origin (oriC). We recently showed that Vibrio cholerae employs this positional bias as a growth optimization strategy: under fast-growth conditions, multifork replication increases RP dosage and expression. However, RP location may provide advantages in a dosage-independent manner: for example, the physical proximity of the many ribosomal components, in the context of a crowded cytoplasm, may favor ribosome biogenesis. To uncover putative dosage-independent effects, we studied isogenic V. cholerae derivatives in which the major RP locus, S10-spc-α (S10), was relocated to alternative genomic positions. When bacteria grew fast, bacterial fitness was reduced according to the S10 relative distance to oriC The growth of wild-type V. cholerae could not be improved by additional copies of the locus, suggesting a physiologically optimized genomic location. Slow growth is expected to uncouple RP position from dosage, since multifork replication does not occur. Under these conditions, we detected a fitness impairment when S10 was far from oriC Deep sequencing followed by marker frequency analysis in the absence of multifork replication revealed an up to 30% S10 dosage reduction associated with its relocation that closely correlated with fitness alterations. Hence, the impact of S10 location goes beyond a growth optimization strategy during feast periods. RP location may be important during the whole life cycle of this pathogen.IMPORTANCE The role of gene order within the bacterial chromosome is poorly understood. In fast growers, the location of genes linked with the expression of genetic information (i.e., transcription and translation) is biased toward oriC It was proposed that the location of these genes helps to maximize their expression by recruiting multifork replication during fast growth. Our results show that such genomic positioning impacts cell fitness beyond fast-growth conditions, probably across the whole life cycle of fast growers. Thus, the genomic position of key highly expressed genes, such as RP, was finely tuned during the evolution of fast-growing bacteria and may also be important in slow growers. In the near future, many more genes whose genomic position impacts bacterial phenotype will be described. These studies will contribute to discovery the rules of genome organization and application of them for the design of synthetic chromosomes and the creation of artificial life forms.

Recent Advances in Helicobacter pylori Replication: Possible Implications in Adaptation to a Pathogenic Lifestyle and Perspectives for Drug Design

DNA replication is an important step in the life cycle of every cell that ensures the continuous flow of genetic information from one generation to the next. In all organisms, chromosome replication must be coordinated with overall cell growth. Helicobacter pylori growth strongly depends on its interaction with the host, particularly with the gastric epithelium. Moreover, H. pylori actively searches for an optimal microniche within a stomach, and it has been shown that not every microniche equally supports growth of this bacterium. We postulate that besides nutrients, H. pylori senses different, unknown signals, which presumably also affect chromosome replication to maintain H. pylori propagation at optimal ratio allowing H. pylori to establish a chronic, lifelong infection. Thus, H. pylori chromosome replication and particularly the regulation of this process might be considered important for bacterial pathogenesis. Here, we summarize our current knowledge of chromosome and plasmid replication in H. pylori and discuss the mechanisms responsible for regulating this key cellular process. The results of extensive studies conducted thus far allow us to propose common and unique traits in H. pylori chromosome replication. Interestingly, the repertoire of proteins involved in replication in H. pylori is significantly different to that in E. coli, strongly suggesting that novel factors are engaged in H. pylori chromosome replication and could represent attractive drug targets.