The transcriptional profile of early to middle sporulation in Bacillus subtilis

The Transcriptomic Response of Cells of the Thermophilic Bacterium Geobacillus icigianus to Terahertz Irradiation

As areas of application of terahertz (THz) radiation expand in science and practice, evidence is accumulating that this type of radiation can affect not only biological molecules directly, but also cellular processes as a whole. In this study, the transcriptome in cells of the thermophilic bacterium Geobacillus icigianus was analyzed immediately after THz irradiation (0.23 W/cm2, 130 μm, 15 min) and at 10 min after its completion. THz irradiation does not affect the activity of heat shock protein genes and diminishes the activity of genes whose products are involved in peptidoglycan recycling, participate in redox reactions, and protect DNA and proteins from damage, including genes of chaperone protein ClpB and of DNA repair protein RadA, as well as genes of catalase and kinase McsB. Gene systems responsible for the homeostasis of transition metals (copper, iron, and zinc) proved to be the most sensitive to THz irradiation; downregulation of these systems increased significantly 10 min after the end of the irradiation. It was also hypothesized that some negative effects of THz radiation on metabolism in G. icigianus cells are related to disturbances in activities of gene systems controlled by metal-sensitive transcription factors.

Complex sporulation-specific expression of transcription termination factor Rho highlights its involvement in Bacillus subtilis cell differentiation

Termination factor Rho, responsible for the main factor-dependent pathway of transcription termination and the major inhibitor of antisense transcription, is an emerging regulator of various physiological processes in microorganisms. In Gram-positive bacterium Bacillus subtilis, Rho is involved in the control of cell adaptation to starvation and, in particular, in the control of sporulation, a complex differentiation program leading to the formation of a highly resistant dormant spore. While the initiation of sporulation requires a decrease in Rho protein levels during the transition to stationary phase, the mechanisms regulating the expression of rho gene throughout the cell cycle remain largely unknown. Here we show that a drop in the activity of the vegetative SigA-dependent rho promoter causes the inhibition of rho expression in stationary phase. However, after the initiation of sporulation, rho gene is specifically reactivated in two compartments of the sporulating cell using distinct mechanisms. In the mother cell, rho expression occurs by read-through transcription initiated at the SigH-dependent promoter of the distal spo0F gene. In the forespore, rho gene is transcribed from the intrinsic promoter recognized by the alternative sigma factor SigF. These regulatory elements ensure the activity of Rho during sporulation, which appears important for the proper formation of spores. We provide experimental evidence that disruption of the spatiotemporal expression of rho during sporulation affects the resistance properties of spores, their morphology, and the ability to return to vegetative growth under favorable growth conditions.

Altruistic feeding and cell-cell signaling during bacterial differentiation actively enhance phenotypic heterogeneity

Starvation triggers bacterial spore formation, a committed differentiation program that transforms a vegetative cell into a dormant spore. Cells in a population enter sporulation nonuniformly to secure against the possibility that favorable growth conditions, which put sporulation-committed cells at a disadvantage, may resume. This heterogeneous behavior is initiated by a passive mechanism: stochastic activation of a master transcriptional regulator. Here, we identify a cell-cell communication pathway containing the proteins ShfA (YabQ) and ShfP (YvnB) that actively promotes phenotypic heterogeneity, wherein Bacillus subtilis cells that start sporulating early use a calcineurin-like phosphoesterase to release glycerol, which simultaneously acts as a signaling molecule and a nutrient to delay nonsporulating cells from entering sporulation. This produced a more diverse population that was better poised to exploit a sudden influx of nutrients compared to those generating heterogeneity via stochastic gene expression alone. Although conflict systems are prevalent among microbes, genetically encoded cooperative behavior in unicellular organisms can evidently also boost inclusive fitness.

Transcription factor AbrB regulates ROS generation and clearance in Bacillus licheniformis

Oxidative damage caused by the accumulation of reactive oxygen species (ROS) is one of the main obstacles to the improvement of microbial cell growth and fermentation characteristics under adverse environments. And the antioxidant capacity of cells will increase with the cell growth. Here, we found that a transition state transcription factor AbrB related to changes in cell growth status could regulate the accumulation of ROS and antioxidant capacity in Bacillus licheniformis. The results showed that the accumulation of intracellular ROS was reduced by 23.91 % and the cell survival rates were increased by 1.77-fold under 0.5 mM H2O2 when AbrB was knocked out. We further mapped regulatory target genes of AbrB related to ROS generation or clearance based on our previously analyzed transcriptome sequencing. It proved that AbrB could promote ROS generation via upregulating the synthesis of oxidase and siderophores, and negatively regulating the synthesis of iron chelators (pulcherriminic acid, and H2S). Additionally, AbrB could inhibit ROS clearance by negatively regulating the synthesis of antioxidase (superoxide dismutase, catalase, peroxidase, thioredoxin, thioredoxin reductase) and cysteine. Those results illustrated that the inactivation of AbrB during the stationary phase, along with its control over ROS generation and clearance, might represent a vital self-protection mechanism during cell evolution. Overall, the systematic investigation of the multi-pathway regulation network of ROS generation and clearance highlights the important function of AbrB in maintaining intracellular redox balance.

Translation in Bacillus subtilis is spatially and temporally coordinated during sporulation

The transcriptional control of sporulation in Bacillus subtilis is reasonably well understood, but its translational control is underexplored. Here, we use RNA-seq, ribosome profiling and fluorescence microscopy to study the translational dynamics of B. subtilis sporulation. We identify two events of translation silencing and describe spatiotemporal changes in subcellular localization of ribosomes during sporulation. We investigate the potential regulatory role of ribosomes during sporulation using a strain lacking zinc-independent paralogs of three zinc-dependent ribosomal proteins (L31, L33 and S14). The mutant strain exhibits delayed sporulation, reduced germination efficiency, dysregulated translation of metabolic and sporulation-related genes, and disruptions in translation silencing, particularly in late sporulation.

Influence of PepF peptidase and sporulation on microcin J25 production in Bacillus subtilis

The lasso peptide microcin J25 (MccJ25) possesses strong antibacterial properties and is considered a potential effective component of bacterial disease treatment drugs and safe food preservatives. Although MccJ25 can be heterologously expressed in Bacillus subtilis as we have previously reported, its regulation and accumulation are yet to be understood. Here, we investigated the expression level and stability of MccJ25 in B. subtilis strains with disruption in peptidase genes pepA, pepF, and pepT. Oligoendopeptidase F (PepF) was found to be involved in reduction of the production of MccJ25 by degradation of its precursor peptide. In the pepF mutant, the MccJ25 reached a concentration of 1.68 µM after a cultivation time exceeding 60 hours, while the wild-type strain exhibited a concentration of only 0.14 µM. Moreover, the production of MccJ25 in B. subtilis downregulated the genes associated with sporulation, and this may contribute to its accumulation. Finally, this study provides a strategy to improve the stability and production of MccJ25 in B. subtilis.

IMPORTANCE

MccJ25 displays significant antibacterial activity, a well-defined mode of action, exceptional safety, and remarkable stability. Hence, it presents itself as a compelling candidate for an optimal antibacterial or anti-endotoxin medication. The successful establishment of exogenous production of MccJ25 in Bacillus subtilis provides a strategy for reducing its production cost and diversifying its utilization. In this study, we have provided evidence indicating that both peptidase PepF and sporulation are significant factors that limit the expression of MccJ25 in B. subtilis. The ΔpepF and ΔsigF mutants of B. subtilis express MccJ25 with higher production yield and enhanced stability. To sum up, this study developed several better engineered strains of B. subtilis, which greatly reduced the consumption of MccJ25 during the nutrient depletion stage of the host strain, improved its production, and elucidated factors that may be involved in reducing MccJ25 accumulation in B. subtilis.

Altruistic feeding and cell-cell signaling during bacterial differentiation actively enhance phenotypic heterogeneity

Starvation triggers bacterial spore formation, a committed differentiation program that transforms a vegetative cell into a dormant spore. Cells in a population enter sporulation non-uniformly to secure against the possibility that favorable growth conditions, which puts sporulation-committed cells at a disadvantage, may resume. This heterogeneous behavior is initiated by a passive mechanism: stochastic activation of a master transcriptional regulator. Here, we identify a cell-cell communication pathway that actively promotes phenotypic heterogeneity, wherein Bacillus subtilis cells that start sporulating early utilize a calcineurin-like phosphoesterase to release glycerol, which simultaneously acts as a signaling molecule and a nutrient to delay non-sporulating cells from entering sporulation. This produced a more diverse population that was better poised to exploit a sudden influx of nutrients compared to those generating heterogeneity via stochastic gene expression alone. Although conflict systems are prevalent among microbes, genetically encoded cooperative behavior in unicellular organisms can evidently also boost inclusive fitness.

Spore-forming properties and enhanced oxygen tolerance of butyrate-producing Anaerostipes spp

OBJECTIVES

Butyrate producing bacteria are promising candidates for next-generation probiotics. However, they are extremely sensitive to oxygen, which is a significant obstacle to their inclusion in food matrices in a viable form. The present study characterized the spore-forming properties and stress tolerance of human gut butyrate-producing Anaerostipes spp.

METHODS

Spore formation properties in six species of Anaerostipes spp. were studied by in vitro and in silico tests.

RESULTS

Spores were observed from the cells of three species using microscopic analyses, while the remaining three did not form spores under the tested conditions. Spore-forming properties were confirmed by an ethanol treatment. The spores of Anaerostipes caccae were tolerant to oxygen and survived for 15 weeks under atmospheric conditions. Spores tolerated heat stress at 70 °C, but not at 80 °C. An in silico analysis of the conservation of potential sporulation signature genes revealed that the majority of human gut butyrate-producing bacteria were classified as potential spore formers. Comparative genomics revealed that three spore-forming Anaerostipes spp. specifically possessed the spore formation-related genes of bkdR, sodA, and splB, which may be key genes for different sporulation properties in Anaerostipes spp.

CONCLUSIONS

The present study demonstrated the enhanced stress tolerance of butyrate producing Anaerostipes spp. for future probiotic application. Presence of specific gene(s) are possibly keys for sporulation in Anaerostipes spp.

Time-Course Transcriptome Analysis of Bacillus subtilis DB104 during Growth

Bacillus subtilis DB104, an extracellular protease-deficient derivative of B. subtilis 168, is widely used for recombinant protein expression. An understanding of the changes in gene expression during growth is essential for the commercial use of bacterial strains. Transcriptome and proteome analyses are ideal methods to study the genomic response of microorganisms. In this study, transcriptome analysis was performed to monitor changes in the gene expression level of B. subtilis DB104 while growing on a complete medium. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, K-mean cluster analysis, gene ontology (GO) enrichment analysis, and the function of sigma factors were used to divide 2122 differentially expressed genes (DEGs) into 10 clusters and identified gene functions according to expression patterns. The results of KEGG pathway analysis indicated that ABC transporter is down-regulated during exponential growth and metabolic changes occur at the transition point where sporulation starts. At this point, several stress response genes were also turned on. The genes involved in the lipid catabolic process were up-regulated briefly at 15 h as an outcome of the programmed cell death that postpones sporulation. The results suggest that changes in the gene expression of B. subtilis DB104 were dependent on the initiation of sporulation. However, the expression timing of the spore coat gene was only affected by the relevant sigma factor. This study can help to understand gene expression and regulatory mechanisms in B. subtilis species by providing an overall view of transcriptional changes during the growth of B. subtilis DB104.

Effects of Omeprazole on Recurrent Clostridioides difficile Infection Caused by ST81 Strains and Their Potential Mechanisms

Clostridioides difficile infection (CDI) is associated with high recurrence rates that have substantial effects on patients' quality of life. To investigate the risk factors and potential mechanisms contributing to recurrent CDI (rCDI), a total of 243 cases were enrolled in this study. The history of omeprazole (OME) medication and ST81 strain infection were considered the two independent risks with the highest odds ratios in rCDI. In the presence of OME, we detected concentration-dependent increases in the MIC values of fluoroquinolone antibiotics against ST81 strains. Mechanically, OME facilitated ST81 strain sporulation and spore germination by blocking the pathway of purine metabolism and also promoted an increase in cell motility and toxin production by turning the flagellar switch to the ON state. In conclusion, OME affects several biological processes during C difficile growth, which have fundamental impacts on the development of rCDI caused by ST81 strains. Programmed OME administration and stringent surveillance of the emerging ST81 genotype are matters of considerable urgency and significance in rCDI prevention.

The putative role of the epipeptide EpeX in Bacillus subtilis intra-species competition

Bacteria engage in competitive interactions with neighbours that can either be of the same or different species. Multiple mechanisms are deployed to ensure the desired outcome and one tactic commonly implemented is the production of specialised metabolites. The Gram-positive bacterium Bacillus subtilis uses specialized metabolites as part of its intra-species competition determinants to differentiate between kin and non-kin isolates. It is, however, unknown if the collection of specialized metabolites defines competitive fitness when the two isolates start as a close, interwoven community that grows into a densely packed colony biofilm. Moreover, the identity of specialized metabolites that have an active role in defining the outcome of an intra-species interaction has not been revealed. Here, we determine the competition outcomes that manifest when 21 environmental isolates of B. subtilis are individually co-incubated with the model isolate NCIB 3610 in a colony biofilm. We correlated these data with the suite of specialized metabolite biosynthesis clusters encoded by each isolate. We found that the epeXEPAB gene cluster was primarily present in isolates with a strong competitive phenotype. This cluster is responsible for producing the epipeptide EpeX. We demonstrated that EpeX is a competition determinant of B. subtilis in an otherwise isogenic context for NCBI 3610. However, when we competed the NCIB 3610 EpeX-deficient strain against our suite of environmental isolates we found that the impact of EpeX in competition is isolate-specific, as only one of the 21 isolates showed increased survival when EpeX was lacking. Taken together, we have shown that EpeX is a competition determinant used by B. subtilis that impacts intra-species interactions but only in an isolate-specific manner.

Engineering of global transcription factors in Bacillus, a genetic tool for increasing product yields: a bioprocess overview

Transcriptional factors are well studied in bacteria for their global interactions and the effects they produce at the phenotypic level. Particularly, Bacillus subtilis has been widely employed as a model Gram-positive microorganism used to characterize these network interactions. Bacillus species are currently used as efficient commercial microbial platforms to produce diverse metabolites such as extracellular enzymes, antibiotics, surfactants, industrial chemicals, heterologous proteins, among others. However, the pleiotropic effects caused by the genetic modification of specific genes that codify for global regulators (transcription factors) have not been implicated commonly from a bioprocess point of view. Recently, these strategies have attracted the attention in Bacillus species because they can have an application to increase production efficiency of certain commercial interest metabolites. In this review, we update the recent advances that involve this trend in the use of genetic engineering (mutations, deletion, or overexpression) performed to global regulators such as Spo0A, CcpA, CodY and AbrB, which can provide an advantage for the development or improvement of bioprocesses that involve Bacillus species as production platforms. Genetic networks, regulation pathways and their relationship to the development of growth stages are also discussed to correlate the interactions that occur between these regulators, which are important to consider for application in the improvement of commercial-interest metabolites. Reported yields from these products currently produced mostly under laboratory conditions and, in a lesser extent at bioreactor level, are also discussed to give valuable perspectives about their potential use and developmental level directed to process optimization at large-scale.

Mechanism of Action of Ribosomally Synthesized and Post-Translationally Modified Peptides

Ribosomally synthesized and post-translationally modified peptides (RiPPs) are a natural product class that has undergone significant expansion due to the rapid growth in genome sequencing data and recognition that they are made by biosynthetic pathways that share many characteristic features. Their mode of actions cover a wide range of biological processes and include binding to membranes, receptors, enzymes, lipids, RNA, and metals as well as use as cofactors and signaling molecules. This review covers the currently known modes of action (MOA) of RiPPs. In turn, the mechanisms by which these molecules interact with their natural targets provide a rich set of molecular paradigms that can be used for the design or evolution of new or improved activities given the relative ease of engineering RiPPs. In this review, coverage is limited to RiPPs originating from bacteria.

Conservation and Evolution of the Sporulation Gene Set in Diverse Members of the Firmicutes

The current classification of the phylum Firmicutes (new name, Bacillota) features eight distinct classes, six of which include known spore-forming bacteria. In Bacillus subtilis, sporulation involves up to 500 genes, many of which do not have orthologs in other bacilli and/or clostridia. Previous studies identified about 60 sporulation genes of B. subtilis that were shared by all spore-forming members of the Firmicutes. These genes are referred to as the sporulation core or signature, although many of these are also found in genomes of nonsporeformers. Using an expanded set of 180 firmicute genomes from 160 genera, including 76 spore-forming species, we investigated the conservation of the sporulation genes, in particular seeking to identify lineages that lack some of the genes from the conserved sporulation core. The results of this analysis confirmed that many small acid-soluble spore proteins (SASPs), spore coat proteins, and germination proteins, which were previously characterized in bacilli, are missing in spore-forming members of Clostridia and other classes of Firmicutes. A particularly dramatic loss of sporulation genes was observed in the spore-forming members of the families Planococcaceae and Erysipelotrichaceae. Fifteen species from diverse lineages were found to carry skin (sigK-interrupting) elements of different sizes that all encoded SpoIVCA-like recombinases but did not share any other genes. Phylogenetic trees built from concatenated alignments of sporulation proteins and ribosomal proteins showed similar topology, indicating an early origin and subsequent vertical inheritance of the sporulation genes. IMPORTANCE Many members of the phylum Firmicutes (Bacillota) are capable of producing endospores, which enhance the survival of important Gram-positive pathogens that cause such diseases as anthrax, botulism, colitis, gas gangrene, and tetanus. We show that the core set of sporulation genes, defined previously through genome comparisons of several bacilli and clostridia, is conserved in a wide variety of sporeformers from several distinct lineages of Firmicutes. We also detected widespread loss of sporulation genes in many organisms, particularly within the families Planococcaceae and Erysipelotrichaceae. Members of these families, such as Lysinibacillus sphaericus and Clostridium innocuum, could be excellent model organisms for studying sporulation mechanisms, such as engulfment, formation of the spore coat, and spore germination.

Mechanisms and Applications of Bacterial Sporulation and Germination in the Intestine

Recent studies have suggested a major role for endospore forming bacteria within the gut microbiota, not only as pathogens but also as commensal and beneficial members contributing to gut homeostasis. In this review the sporulation processes, spore properties, and germination processes will be explained within the scope of the human gut. Within the gut, spore-forming bacteria are known to interact with the host’s immune system, both in vegetative cell and spore form. Together with the resistant nature of the spore, these characteristics offer potential for spores’ use as delivery vehicles for therapeutics. In the last part of the review, the therapeutic potential of spores as probiotics, vaccine vehicles, and drug delivery systems will be discussed.

Calcium Prevents Biofilm Dispersion in Bacillus subtilis

Biofilm dispersion is the final stage of biofilm development, during which biofilm cells actively escape from biofilms in response to deteriorating conditions within the biofilm. Biofilm dispersion allows cells to spread to new locations and form new biofilms in better locations. However, dispersal mechanisms have been elucidated only in a limited number of bacteria. Here, we investigated biofilm dispersion in Bacillus subtilis. Biofilm dispersion was clearly observed when B. subtilis was grown under static conditions in modified LB medium containing glycerol and manganese. Biofilm dispersion was synergistically caused by two mechanisms: decreased expression of the epsA operon encoding exopolysaccharide synthetases and the induction of sporulation. Indeed, constitutive expression of the epsA operon in the sporulation-defective ΔsigK mutant prevented biofilm dispersion. The addition of calcium to the medium prevented biofilm dispersion without significantly affecting the expression of the epsA operon and sporulation genes. In synthetic medium, eliminating calcium did not prevent the expression of biofilm matrix genes and, thereby, biofilm formation, but it attenuated biofilm architecture. These results indicate that calcium structurally stabilizes biofilms and causes resistance to biofilm dispersion mechanisms. Sporulation-dependent biofilm dispersion required the spoVF operon, encoding dipicolinic acid (DPA) synthase. During sporulation, an enormous amount of DPA is synthesized and stored in spores as a chelate with calcium. We speculate that, during sporulation, calcium bound to biofilm matrix components may be transported to spores as a calcium-DPA complex, which weakens biofilm structure and leads to biofilm dispersion. IMPORTANCE Bacteria growing as biofilms are notoriously difficult to eradicate and sometimes pose serious threats to public health. Bacteria escape from biofilms by degrading them when biofilm conditions deteriorate. This process, called biofilm dispersion, has been studied as a promising strategy for safely controlling biofilms. However, the regulation and mechanism of biofilm dispersion has been elucidated only in a limited number of bacteria. Here, we identified two biofilm dispersion mechanisms in the Gram-positive, spore-forming bacterium Bacillus subtilis. The addition of calcium to the medium stabilized biofilms and caused resistance to dispersal mechanisms. Our findings provide new insights into biofilm dispersion and biofilm control.

Sporofaciens musculi gen. nov., sp. nov., a novel bacterium isolated from the caecum of an obese mouse

A bacterial strain, designated WCA-9-b2^T, was isolated from the caecal content of an 18-week-old obese C57BL/6NTac male mouse. According to phenotypic analyses, the isolate was rod-shaped, strictly anaerobic, spore-forming, non-motile and Gram-stain-positive, under the conditions tested. Colonies were irregular and non-pigmented. Analysis of the 16S rRNA gene sequence indicated that the isolate belonged to the order Clostridiales with Dorea longicatena ATCC 27755^T (94.9 % sequence identity), Ruminococcus gnavus ATCC 29149^T (94.8%) and Clostridium scindens ATCC 35704^T (94.3%) being the closest relatives. Whole genome sequencing showed an average nucleotide identity <74.23 %, average amino acid identity <64.52-74.67 % and percentage of conserved proteins values <50 % against the nine closest relatives (D. longicatena, Ruminococcus gnavus, C. scindens, Dorea formicigenerans, Ruminococcus lactaris, Clostridium hylemonae, Merdimonas faecis, Faecalicatena contorta and Faecalicatena fissicatena). The genome-based G+C content of genomic DNA was 44.4 mol%. The major cellular fatty acids were C_16 : 0 (24.5%), C_18 : 1 cis9 (19.8 %), C_16 : 0 DMA (11.7%), C_18 : 0 (8.4%) and C_14 : 0 (6.6%). Respiratory quinones were not detected. The predominant metabolic end products of glucose fermentation were acetate and succinate. Production of CO₂ and H₂ were detected. Based on these data, we propose that strain WCA-9-b2^T represents a novel species within a novel genus, for which the name Sporofaciens musculi gen. nov., sp. nov. is proposed. The type strain is WCA-9-b2^T (=DSM 106039^T=CECT 30156^T).

Antibacterial Activity of Bacillus inaquosorum Strain T1 against pirABVp -Bearing Vibrio parahaemolyticus: Genetic and Physiological Characterization

Acute hepatopancreatic necrosis disease (AHPND) is caused by PirAB toxin-producing Vibrio parahaemolyticus and has devastated the global shrimp aquaculture industry. One approach for preventing the growth of AHPND-producing Vibrio spp. is through the application of beneficial bacteria capable of inhibiting these pathogens. In this study, we focused on the inhibitory activity of Bacillus inaquosorum strain T1, which hinders V. parahaemolyticus growth in coculture experiments in a density-dependent manner; inhibition was also observed using cell-free supernatants from T1 stationary-phase cultures. Using mariner-based transposon mutagenesis, 17 mutants having a complete or partial loss of inhibitory activity were identified. Of those displaying a total loss of activity, 13 had insertions within a 42.6-kb DNA region comprising 15 genes whose deduced products were homologous to nonribosomal polypeptide synthetases (NRPSs), polyketide synthases (PKSs), and related activities, which were mapped as one transcriptional unit. Mutants with partial activity contained insertions in spo0A and oppA, indicating stationary-phase control. The levels of expression of NRPS and PKS lacZ transcriptional fusions were negligible during growth and were the highest during early stationary phase. Inactivation of sigH resulted in a loss of inhibitor activity, indicating a role for σ^H in transcription. Disruption of abrB resulted in NRPS and PKS gene overexpression during growth as well as enhanced growth inhibition. Our characterization of the expression and control of an NRPS-PKS gene cluster in B. inaquosorum T1 provides an understanding of the factors involved in inhibitor production, enabling this strain's development for use as a tool against AHPND-causing Vibrio pathogens in shrimp aquaculture.IMPORTANCE The shrimp aquaculture industry has been significantly impacted by acute hepatopancreatic necrosis disease (AHPND), resulting in significant financial losses annually. AHPND is caused by strains of the bacterial pathogen Vibrio parahaemolyticus, and treatment of AHPND involves the use of antibiotics, which leads to a rise in the number of antibiotic-resistant strains. Alternative treatments include the application of beneficial microorganisms having inhibitory activities against pathogens causing AHPND. In this study, we examined the ability of Bacillus inaquosorum strain T1 to inhibit the growth of an AHPND-causing Vibrio strain, and we show that this activity involves a gene cluster associated with antibacterial compound production. We found that gene expression is under stationary-phase control and that enhanced activity occurs upon inactivation of a global transition state regulator. Our approach for understanding the factors involved in producing B. inaquosorum strain T1 inhibitory activity will allow for the development of this strain as a tool for AHPND prevention and treatment.

Inference of Bacterial Small RNA Regulatory Networks and Integration with Transcription Factor-Driven Regulatory Networks

Small noncoding RNAs (sRNAs) are key regulators of bacterial gene expression. Through complementary base pairing, sRNAs affect mRNA stability and translation efficiency. Here, we describe a network inference approach designed to identify sRNA-mediated regulation of transcript levels. We use existing transcriptional data sets and prior knowledge to infer sRNA regulons using our network inference tool, the Inferelator. This approach produces genome-wide gene regulatory networks that include contributions by both transcription factors and sRNAs. We show the benefits of estimating and incorporating sRNA activities into network inference pipelines using available experimental data. We also demonstrate how these estimated sRNA regulatory activities can be mined to identify the experimental conditions where sRNAs are most active. We uncover 45 novel experimentally supported sRNA-mRNA interactions in Escherichia coli, outperforming previous network-based efforts. Additionally, our pipeline complements sequence-based sRNA-mRNA interaction prediction methods by adding a data-driven filtering step. Finally, we show the general applicability of our approach by identifying 24 novel, experimentally supported, sRNA-mRNA interactions in Pseudomonas aeruginosa, Staphylococcus aureus, and Bacillus subtilis. Overall, our strategy generates novel insights into the functional context of sRNA regulation in multiple bacterial species.

IMPORTANCE Individual bacterial genomes can have dozens of small noncoding RNAs with largely unexplored regulatory functions. Although bacterial sRNAs influence a wide range of biological processes, including antibiotic resistance and pathogenicity, our current understanding of sRNA-mediated regulation is far from complete. Most of the available information is restricted to a few well-studied bacterial species; and even in those species, only partial sets of sRNA targets have been characterized in detail. To close this information gap, we developed a computational strategy that takes advantage of available transcriptional data and knowledge about validated and putative sRNA-mRNA interactions for inferring expanded sRNA regulons. Our approach facilitates the identification of experimentally supported novel interactions while filtering out false-positive results. Due to its data-driven nature, our method prioritizes biologically relevant interactions among lists of candidate sRNA-target pairs predicted in silico from sequence analysis or derived from sRNA-mRNA binding experiments.

From Root to Tips: Sporulation Evolution and Specialization in Bacillus subtilis and the Intestinal Pathogen Clostridioides difficile

Bacteria of the Firmicutes phylum are able to enter a developmental pathway that culminates with the formation of highly resistant, dormant endospores. Endospores allow environmental persistence, dissemination and for pathogens, are also infection vehicles. In both the model Bacillus subtilis, an aerobic organism, and in the intestinal pathogen Clostridioides difficile, an obligate anaerobe, sporulation mobilizes hundreds of genes. Their expression is coordinated between the forespore and the mother cell, the two cells that participate in the process, and is kept in close register with the course of morphogenesis. The evolutionary mechanisms by which sporulation emerged and evolved in these two species, and more broadly across Firmicutes, remain largely unknown. Here, we trace the origin and evolution of sporulation using the genes known to be involved in the process in B. subtilis and C. difficile, and estimating their gain-loss dynamics in a comprehensive bacterial macroevolutionary framework. We show that sporulation evolution was driven by two major gene gain events, the first at the base of the Firmicutes and the second at the base of the B. subtilis group and within the Peptostreptococcaceae family, which includes C. difficile. We also show that early and late sporulation regulons have been coevolving and that sporulation genes entail greater innovation in B. subtilis with many Bacilli lineage-restricted genes. In contrast, C. difficile more often recruits new sporulation genes by horizontal gene transfer, which reflects both its highly mobile genome, the complexity of the gut microbiota, and an adjustment of sporulation to the gut ecosystem.

Transcriptome dynamics of the Myxococcus xanthus multicellular developmental program

The bacterium Myxococcus xanthus exhibits a complex multicellular life cycle. In the presence of nutrients, cells prey cooperatively. Upon starvation, they enter a developmental cycle wherein cells aggregate to produce macroscopic fruiting bodies filled with resistant myxospores. We used RNA-Seq technology to examine the transcriptome of the 96 hr developmental program. These data revealed that 1415 genes were sequentially expressed in 10 discrete modules, with expression peaking during aggregation, in the transition from aggregation to sporulation, or during sporulation. Analysis of genes expressed at each specific time point provided insights as to how starving cells obtain energy and precursors necessary for assembly of fruiting bodies and into developmental production of secondary metabolites. This study offers the first global view of developmental transcriptional profiles and provides important tools and resources for future studies.

Laboratory strains of Bacillus anthracis exhibit pervasive alteration in expression of proteins related to sporulation under laboratory conditions relative to genetically related wild strains

The spore forming pathogen Bacillus anthracis is the etiologic agent of anthrax in humans and animals. It cycles through infected hosts as vegetative cells and is eventually introduced into the environment where it generates an endospore resistant to many harsh conditions. The endospores are subsequently taken up by another host to begin the next cycle. Outbreaks of anthrax occur regularly worldwide in wildlife and livestock, and the potential for human infection exists whenever humans encounter infected animals. It is also possible to encounter intentional releases of anthrax spores, as was the case in October 2001. Consequently, it is important to be able to rapidly establish the provenance of infectious strains of B. anthracis. Here, we compare protein expression in seven low-passage wild isolates and four laboratory strains of B. anthracis grown under identical conditions using LC-MS/MS proteomic analysis. Of the 1,023 total identified proteins, 96 had significant abundance differences between wild and laboratory strains. Of those, 28 proteins directly related to sporulation were upregulated in wild isolates, with expression driven by Spo0A, CodY, and AbrB/ScoC. In addition, we observed evidence of changes in cell division and fatty acid biosynthesis between the two classes of strains, despite being grown under identical experimental conditions. These results suggest wild B. anthracis cells are more highly tuned to sporulate than their laboratory cousins, and this difference should be exploited as a method to differentiate between laboratory and low passage wild strains isolated during an anthrax outbreak. This knowledge should distinguish between intentional releases and exposure to strains in nature, providing a basis for the type of response by public health officials and investigators.

Non-canonical processing of DNA photodimers with Bacillus subtilis UV-endonuclease YwjD, 5'→3' exonuclease YpcP and low-fidelity DNA polymerases YqjH and YqjW

It has been shown that mutation frequency decline (Mfd) and nucleotide excision repair (NER) deficiencies promote UV-induced mutagenesis in B. subtilis sporangia. As replication is halted in sporulating B. subtilis cells, in this report, we investigated if this response may result from an error-prone repair event involving the UV-endonuclease YwjD and low fidelity (LF) DNA synthesis. Accordingly, disruption of YwjD generated B. subtilis sporangia that were more susceptible to UV-C radiation than sporangia of the WT strain and such susceptibility increased even more after the single or simultaneous inactivation of the LF DNA polymerases YqjH and YqjW. To further explore this concept, functional His₆-tagged YwjD and Y-DNA polymerases YqjH and YqjW were produced and purified to homogeneity. In vitro repair assays showed that YwjD hydrolyzed the phosphodiester bond immediately located 5´-end of a ds-DNA substrate bearing either, cyclobutane pyrimidine dimers (CPD), 6-4 photoproducts (6-4 PD) or Dewar isomers (DWI). Notably, the 6-4 PD and DWI but not the CPDs repair intermediaries of YwjD were efficiently processed by the LF polymerase YqjH suggesting that an additional 5'→3' exonuclease event was necessary to process PD. Accordingly, the LF polymerase YqjW efficiently processed the incision-excision repair products derived from YwjD and exonuclease YpcP attack over CPD-containing DNA. In summary, our results unveiled a novel non-canonical repair pathway that employs YwjD to incise PD-containing DNA and low fidelity synthesis contributing thus to mutagenesis, survival and spore morphogenesis in B. subtilis.

Random knock-in expression system for high yield production of heterologous protein in Bacillus subtilis

Chromosome-integrated recombinant protein expression in bacteria has advantages for the stable maintenance of genes without any use of antibiotics during large-scale fermentation. Even though different levels of gene expression were reported, depending upon their chromosomal position in bacterial species, only a limited number of integration sites have been used in B. subtilis. In this study, we randomly integrated the GFP and AprE expression cassettes into the B. subtilis genome to determine integration sites that can produce a high yield of heterologous protein expression. Our mariner transposon-based expression cassette integration system was able to find integration sites, which can produce up to 2.9-fold and 1.5-fold increased expression of intracellular GFP and extracellular AprE, respectively, compared to the common integration site amyE. By analyzing the location of integration sites, we observed an adjacent promoter effect, gene dosage effect, and gene knock-out effect all complexly contributing to the increased level of integrated gene expression. Besides obtaining a high yield of heterologous protein expression, our system can also provide a wide-range of expression to expand the systematic application for steady-state metabolic protein production.

A Mutation in the Bacillus subtilis rsbU Gene That Limits RNA Synthesis during Sporulation

Mutants of Bacillis subtilis that are temperature sensitive for RNA synthesis during sporulation were isolated after selection with a ³²P suicide agent. Whole-genome sequencing revealed that two of the mutants carried an identical lesion in the rsbU gene, which encodes a phosphatase that indirectly activates SigB, the stress-responsive RNA polymerase sigma factor. The mutation appeared to cause RsbU to be hyperactive, because the mutants were more resistant than the parent strain to ethanol stress. In support of this hypothesis, pseudorevertants that regained wild-type levels of sporulation at high temperature had secondary mutations that prevented expression of the mutant rsbU gene. The properties of these RsbU mutants support the idea that activation of SigB diminishes the bacterium's ability to sporulate.IMPORTANCE Most bacterial species encode multiple RNA polymerase promoter recognition subunits (sigma factors). Each sigma factor directs RNA polymerase to different sets of genes; each gene set typically encodes proteins important for responses to specific environmental conditions, such as changes in temperature, salt concentration, and nutrient availability. A selection for mutants of Bacillus subtilis that are temperature sensitive for RNA synthesis during sporulation unexpectedly yielded strains with a point mutation in rsbU, a gene that encodes a protein that normally activates sigma factor B (SigB) under conditions of salt stress. The mutation appears to cause RsbU, and therefore SigB, to be active inappropriately, thereby inhibiting, directly or indirectly, the ability of the cells to transcribe sporulation genes.

Cell Death Pathway That Monitors Spore Morphogenesis

The use of quality control mechanisms to stall developmental pathways or completely remove defective cells from a population is a widespread strategy to ensure the integrity of morphogenetic programs. Endospore formation (sporulation) is a well conserved microbial developmental strategy in the Firmicutes phylum wherein a progenitor cell that faces starvation differentiates to form a dormant spore. Despite the conservation of this strategy, it has been unclear what selective pressure maintains the fitness of this developmental program, composed of hundreds of unique genes, during multiple rounds of vegetative growth when sporulation is not required. Recently, a quality control pathway was discovered in Bacillus subtilis which monitors the assembly of the spore envelope and specifically eliminates, through cell lysis, sporulating cells that assemble the envelope incorrectly. Here, we review the use of checkpoints that govern the entry into sporulation in B. subtilis and discuss how the use of regulated cell death pathways during bacterial development may help maintain the fidelity of the sporulation program in the species.

Crystal structure of the inactive state of the receiver domain of Spo0A from Paenisporosarcina sp. TG-14, a psychrophilic bacterium isolated from an Antarctic glacier

The two-component phosphorelay system is the most prevalent mechanism for sensing and transducing environmental signals in bacteria. Spore formation, which relies on the two-component phosphorelay system, enables the long-term survival of the glacial bacterium Paenisporosarcina sp. TG-14 in the extreme cold environment. Spo0A is a key response regulator of the phosphorelay system in the early stage of spore formation. The protein is composed of a regulatory N-terminal phospho-receiver domain and a DNA-binding C-terminal activator domain. We solved the three-dimensional structure of the unphosphorylated (inactive) form of the receiver domain of Spo0A (PaSpo0A-R) from Paenisporosarcina sp. TG-14. A structural comparison with phosphorylated (active form) Spo0A from Bacillus stearothermophilus (BsSpo0A) showed minor notable differences. A molecular dynamics study of a model of the active form and the crystal structures revealed significant differences in the α4 helix and the preceding loop region where phosphorylation occurs. Although an oligomerization study of PaSpo0A-R by analytical ultracentrifugation (AUC) has shown that the protein is in a monomeric state in solution, both crosslinking and crystal-packing analyses indicate the possibility of weak dimer formation by a previously undocumented mechanism. Collectively, these observations provide insight into the mechanism of phosphorylation-dependent activation unique to Spo0A.

Metabolic flexibility of a butyrate pathway mutant of Clostridium acetobutylicum

Clostridium acetobutylicum possesses two homologous buk genes, buk (or buk1) and buk2, which encode butyrate kinases involved in the last step of butyrate formation. To investigate the contribution of buk in detail, an in-frame deletion mutant was constructed. However, in all the Δbuk mutants obtained, partial deletions of the upstream ptb gene were observed, and low phosphotransbutyrylase and butyrate kinase activities were measured. This demonstrates that i) buk (CA_C3075) is the key butyrate kinase-encoding gene and that buk2 (CA_C1660) that is poorly transcribed only plays a minor role; and ii) strongly suggests that a Δbuk mutant is not viable if the ptb gene is not also inactivated, probably due to the accumulation of butyryl-phosphate, which might be toxic for the cell. One of the ΔbukΔptb mutants was subjected to quantitative transcriptomic (mRNA molecules/cell) and fluxomic analyses in acidogenic, solventogenic and alcohologenic chemostat cultures. In addition to the low butyrate production, drastic changes in metabolic fluxes were also observed for the mutant: i) under acidogenic conditions, the primary metabolite was butanol and a new metabolite, 2-hydroxy-valerate, was produced ii) under solventogenesis, 58% increased butanol production was obtained compared to the control strain under the same conditions, and a very high yield of butanol formation (0.3gg^-1) was reached; and iii) under alcohologenesis, the major product was lactate. Furthermore, at the transcriptional level, adhE2, which encodes an aldehyde/alcohol dehydrogenase and is known to be a gene specifically expressed in alcohologenesis, was surprisingly highly expressed in all metabolic states in the mutant. The results presented here not only support the key roles of buk and ptb in butyrate formation but also highlight the metabolic flexibility of C. acetobutylicum in response to genetic alteration of its primary metabolism.

Sigma Factor Regulated Cellular Response in a Non-solvent Producing Clostridium beijerinckii Degenerated Strain: A Comparative Transcriptome Analysis

Clostridium beijerinckii DG-8052, derived from NCIMB 8052, cannot produce solvent or form spores, a phenomenon known as degeneration. To explore the mechanisms of degeneration at the gene level, transcriptomic profiles of the wild-type 8052 and DG-8052 strains were compared. Expression of 5168 genes comprising 98.6% of the genome was assessed. Interestingly, 548 and 702 genes were significantly up-regulated in the acidogenesis and solventogenesis phases of DG-8052, respectively, and mainly responsible for the phosphotransferase system, sugar metabolic pathways, and chemotaxis; meanwhile, 699 and 797 genes were significantly down-regulated, respectively, and mainly responsible for sporulation, oxidoreduction, and solventogenesis. The functions of some altered genes, including 286 and 333 at the acidogenesis and solventogenesis phases, respectively, remain unknown. Dysregulation of the fermentation machinery was accompanied by lower transcription levels of glycolysis rate-limiting enzymes (pfk and pyk), and higher transcription of cell chemotaxis genes (cheA, cheB, cheR, cheW, and cheY), controlled mainly by σ⁵⁴ at acidogenesis. Meanwhile, abnormal spore formation was associated with repressed spo0A, sigE, sigF, sigG, and sigK which are positively regulated by σ⁷⁰, and correspondingly inhibited expression of CoA-transferase at the solventogenesis phase. These findings indicated that morphological and physiological changes in the degenerated Clostridium strain may be related to altered expression of sigma factors, providing valuable targets for strain development of Clostridium species.

Role of Supercomputers in Bioinformatics

Due to the involvement of effective and client-friendly components (i.e. supercomputers), rapid data analysis is being accomplished. In Bioinformatics, it is expanding many areas of research such as genomics, proteomics, metabolomics, etc. Structure-based drug design is one of the major areas of research to cure human malady. This chapter initiates a discussion on supercomputing in sequence analysis with a detailed table summarizing the software and Web-based programs used for sequence analysis. A brief talk on the supercomputing in virtual screening is given where the databases like DOCK, ZINC, EDULISS, etc. are introduced. As the chapter transitions to the next phase, the intricacies of advanced Quantitative Structure-Activity Relationship technologies like Fragment-Based 2D QSAR, Multiple-Field 3D QSAR, and Amino Acid-Based Peptide Prediction are put forth in a manner similar to the concept of abstraction. The supercomputing in docking studies is stressed where docking software for Protein-Ligand docking, Protein-Protein docking, and Multi-Protein docking are provided. The chapter ends with the applications of supercomputing in widely used microarray data analysis.

A Novel Cell Type Enables B. subtilis to Escape from Unsuccessful Sporulation in Minimal Medium

Sporulation is the most enduring survival strategy developed by several bacterial species. However, spore development of the model organism Bacillus subtilis has mainly been studied by means of media or conditions optimized for the induction of sporogenesis. Here, I show that during prolonged growth during stationary phase in minimal medium, B. subtilis undergoes an asymmetric cell division that produces small and round-shaped, DNA containing cells. In contrast to wild-type cells, mutants harboring spo0A or spoIIIE/sftA double mutations neither sporulate nor produce this special cell type, providing evidence that the small round cells emerge from the abortion of endospore formation. In most cases observed, the small round cells arise in the presence of sigma H but absence of sigma F activity, different from cases of abortive sporulation described for rich media. These data suggest that in minimal media, many cells are able to initiate but fail to complete spore development, and therefore return to normal growth as rods. This work reveals that the continuation of asymmetric cell division, which results in the formation of the small round cells, is a way for cells to delay or escape from—unsuccessful—sporulation. Based on these findings, I suggest to name the here described cell type as “dwarf cells” to distinguish them from the well-known minicells observed in mutants defective in septum placement or proper chromosome partitioning.

Time Series Analysis of the Bacillus subtilis Sporulation Network Reveals Low Dimensional Chaotic Dynamics

Chaotic behavior refers to a behavior which, albeit irregular, is generated by an underlying deterministic process. Therefore, a chaotic behavior is potentially controllable. This possibility becomes practically amenable especially when chaos is shown to be low-dimensional, i.e., to be attributable to a small fraction of the total systems components. In this case, indeed, including the major drivers of chaos in a system into the modeling approach allows us to improve predictability of the systems dynamics. Here, we analyzed the numerical simulations of an accurate ordinary differential equation model of the gene network regulating sporulation initiation in Bacillus subtilis to explore whether the non-linearity underlying time series data is due to low-dimensional chaos. Low-dimensional chaos is expectedly common in systems with few degrees of freedom, but rare in systems with many degrees of freedom such as the B. subtilis sporulation network. The estimation of a number of indices, which reflect the chaotic nature of a system, indicates that the dynamics of this network is affected by deterministic chaos. The neat separation between the indices obtained from the time series simulated from the model and those obtained from time series generated by Gaussian white and colored noise confirmed that the B. subtilis sporulation network dynamics is affected by low dimensional chaos rather than by noise. Furthermore, our analysis identifies the principal driver of the networks chaotic dynamics to be sporulation initiation phosphotransferase B (Spo0B). We then analyzed the parameters and the phase space of the system to characterize the instability points of the network dynamics, and, in turn, to identify the ranges of values of Spo0B and of the other drivers of the chaotic dynamics, for which the whole system is highly sensitive to minimal perturbation. In summary, we described an unappreciated source of complexity in the B. subtilis sporulation network by gathering evidence for the chaotic behavior of the system, and by suggesting candidate molecules driving chaos in the system. The results of our chaos analysis can increase our understanding of the intricacies of the regulatory network under analysis, and suggest experimental work to refine our behavior of the mechanisms underlying B. subtilis sporulation initiation control.

Coordination of cell decisions and promotion of phenotypic diversity in B. subtilis via pulsed behavior of the phosphorelay

The phosphorelay of Bacillus subtilis, a kinase cascade that activates master regulator Spo0A ~ P in response to starvation signals, is the core of a large network controlling the cell's decision to differentiate into sporulation and other phenotypes. This article reviews recent advances in understanding the origins and purposes of the complex dynamical behavior of the phosphorelay, which pulses with peaks of activity coordinated with the cell cycle. The transient imbalance in the expression of two critical genes caused by their strategic placement at opposing ends of the chromosome proved to be the key for this pulsed behavior. Feedback control loops in the phosphorelay use these pulses to implement a timer mechanism, which creates several windows of opportunity for phenotypic transitions over multiple generations. This strategy allows the cell to coordinate multiple differentiation programs in a decision process that fosters phenotypic diversity and adapts to current conditions.

Salt-sensitivity of σ(H) and Spo0A prevents sporulation of Bacillus subtilis at high osmolarity avoiding death during cellular differentiation

The spore-forming bacterium Bacillus subtilis frequently experiences high osmolarity as a result of desiccation in the soil. The formation of a highly desiccation-resistant endospore might serve as a logical osmostress escape route when vegetative growth is no longer possible. However, sporulation efficiency drastically decreases concomitant with an increase in the external salinity. Fluorescence microscopy of sporulation-specific promoter fusions to gfp revealed that high salinity blocks entry into the sporulation pathway at a very early stage. Specifically, we show that both Spo0A- and SigH-dependent transcription are impaired. Furthermore, we demonstrate that the association of SigH with core RNA polymerase is reduced under these conditions. Suppressors that modestly increase sporulation efficiency at high salinity map to the coding region of sigH and in the regulatory region of kinA, encoding one the sensor kinases that activates Spo0A. These findings led us to discover that B. subtilis cells that overproduce KinA can bypass the salt-imposed block in sporulation. Importantly, these cells are impaired in the morphological process of engulfment and late forespore gene expression and frequently undergo lysis. Altogether our data indicate that B. subtilis blocks entry into sporulation in high-salinity environments preventing commitment to a developmental program that it cannot complete.

Structure of Spo0M, a sporulation-control protein from Bacillus subtilis

Spo0M is a sporulation-control protein that is thought to play an essential role in the early stage of endospore formation. While little is known about the functions of Spo0M, a recent phylogenetic study suggests that, based on its amino-acid sequence, Spo0M might belong to the arrestin clan. The crystal structure of the Spo0M protein was determined at a resolution of 2.3 Å. Ten amino acids at the end of the N-terminus were removed to improve the thermal stability of the purified Spo0M protein and the crystal structure of Spo0M was determined by SAD. Spo0M has a well conserved N-terminal domain with an arrestin-like fold, which consists of a β-strand sandwich structure. Surprisingly, the C-terminal domain of Spo0M, which has no structural homology to arrestin-clan proteins, bears significant structural similarity to the FP domain of the human PI31 protein. In addition, Spo0M harbours a potential polar-core structure connecting the N- and C-terminal domains with several salt bridges, as seen in the crystal structures of arrestin and VPS26. The structure reported here constitutes the first structural information on a bacterial protein that shares significant structural homology to members of the arrestin clan and the FP domain.

An experimentally supported model of the Bacillus subtilis global transcriptional regulatory network

Organisms from all domains of life use gene regulation networks to control cell growth, identity, function, and responses to environmental challenges. Although accurate global regulatory models would provide critical evolutionary and functional insights, they remain incomplete, even for the best studied organisms. Efforts to build comprehensive networks are confounded by challenges including network scale, degree of connectivity, complexity of organism–environment interactions, and difficulty of estimating the activity of regulatory factors. Taking advantage of the large number of known regulatory interactions in Bacillus subtilis and two transcriptomics datasets (including one with 38 separate experiments collected specifically for this study), we use a new combination of network component analysis and model selection to simultaneously estimate transcription factor activities and learn a substantially expanded transcriptional regulatory network for this bacterium. In total, we predict 2,258 novel regulatory interactions and recall 74% of the previously known interactions. We obtained experimental support for 391 (out of 635 evaluated) novel regulatory edges (62% accuracy), thus significantly increasing our understanding of various cell processes, such as spore formation.

Transcriptional Profile of Bacillus subtilis sigF-Mutant during Vegetative Growth

Sigma factor F is the first forespore specific transcription factor in Bacillus subtilis and controls genes required for the early stages of prespore development. The role of sigF is well studied under conditions that induce sporulation. Here, the impact of sigF disruption on the transcriptome of exponentially growing cultures is studied by micro-array analysis. Under these conditions that typically don’t induce sporulation, the transcriptome showed minor signs of sporulation initiation. The number of genes differentially expressed and the magnitude of expression were, as expected, quite small in comparison with sporulation conditions. The genes mildly down-regulated were mostly involved in anabolism and the genes mildly up-regulated, in particular fatty acid degradation genes, were mostly involved in catabolism. This is probably related to the arrest at sporulation stage II occurring in the sigF mutant, because continuation of growth from the formed disporic sporangia may require additional energy. The obtained knowledge is relevant for various experiments, such as industrial fermentation, prolonged experimental evolution or zero-growth studies, where sporulation is an undesirable trait that should be avoided, e.g by a sigF mutation.

Antimicrobial copper alloy surfaces are effective against vegetative but not sporulated cells of gram-positive Bacillus subtilis

This study explores the role of membrane phospholipid peroxidation in the copper alloy mediated contact killing of Bacillus subtilis, a spore-forming gram-positive bacterial species. We found that B. subtilis endospores exhibited significant resistance to copper alloy surface killing but vegetative cells were highly sensitive to copper surface exposure. Cell death and lipid peroxidation occurred in B. subtilis upon copper alloy surface exposure. In a sporulation-defective strain carrying a deletion of almost the entire SpoIIA operon, lipid peroxidation directly correlated with cell death. Moreover, killing and lipid peroxidation initiated immediately and at a constant rate upon exposure to the copper surface without the delay observed previously in E. coli. These findings support the hypothesis that membrane lipid peroxidation is the initiating event causing copper surface induced cell death of B. subtilis vegetative cells. The findings suggest that the observed differences in the kinetics of copper-induced killing compared to E. coli result from differences in cell envelop structure. As demonstrated in E. coli, DNA degradation was shown to be a secondary effect of copper exposure in a B. subtilis sporulation-defective strain.

A plasmid-born Rap-Phr system regulates surfactin production, sporulation and genetic competence in the heterologous host, Bacillus subtilis OKB105

According to the change of environment, soil-dwelling Bacillus species differentiate into distinct subpopulations, such as spores and competent cells. Rap-Phr systems have been found to be involved in this differentiation circuit by interacting with major regulatory proteins, such as Spo0A, ComA, and DegU. In this study, we report that the plasmid-born RapQ-PhrQ system found in Bacillus amyloliquefaciens B3 affects three regulatory pathways in the heterologous host Bacillus subtilis. Expression of rapQ in B. subtilis OKB105 strongly suppressed its sporulation efficiency, transformation efficiency, and surfactin production. Co-expression of phrQ or addition of synthesized PhrQ pentapeptide in vitro could compensate for the suppressive effects caused by rapQ. We also found that expression of rapQ decreased the transcriptional level of the sporulation-related gene spoIIE and surfactin synthesis-related gene srfA; meanwhile, the transcriptional levels of these genes could be rescued by co-expression of phrQ and in vitro addition of PhrQ pentapeptide. Electrophoretic mobility shift (EMSA) result also showed that RapQ could bind to ComA without interacting with ComA binding to DNA, and PhrQ pentapeptide antagonized RapQ activity in vitro. These results indicate that this new plasmid-born RapQ-PhrQ system controls sporulation, competent cell formation, and surfactin production in B. subtilis OKB105.

High-throughput analysis of gene essentiality and sporulation in Clostridium difficile

Clostridium difficile is the most common cause of antibiotic-associated intestinal infections and a significant cause of morbidity and mortality. Infection with C. difficile requires disruption of the intestinal microbiota, most commonly by antibiotic usage. Therapeutic intervention largely relies on a small number of broad-spectrum antibiotics, which further exacerbate intestinal dysbiosis and leave the patient acutely sensitive to reinfection. Development of novel targeted therapeutic interventions will require a detailed knowledge of essential cellular processes, which represent attractive targets, and species-specific processes, such as bacterial sporulation. Our knowledge of the genetic basis of C. difficile infection has been hampered by a lack of genetic tools, although recent developments have made some headway in addressing this limitation. Here we describe the development of a method for rapidly generating large numbers of transposon mutants in clinically important strains of C. difficile. We validated our transposon mutagenesis approach in a model strain of C. difficile and then generated a comprehensive transposon library in the highly virulent epidemic strain R20291 (027/BI/NAP1) containing more than 70,000 unique mutants. Using transposon-directed insertion site sequencing (TraDIS), we have identified a core set of 404 essential genes, required for growth in vitro. We then applied this technique to the process of sporulation, an absolute requirement for C. difficile transmission and pathogenesis, identifying 798 genes that are likely to impact spore production. The data generated in this study will form a valuable resource for the community and inform future research on this important human pathogen.

Clostridium difficile is a common cause of potentially fatal intestinal infections in hospital patients, particularly those who have been treated with antibiotics. Our knowledge of this bacterium has been hampered by a lack of tools for dissecting the organism. We have developed a method to study the function of every gene in the bacterium simultaneously. Using this tool, we have identified a set of 404 genes that are required for growth of the bacteria in the laboratory. C. difficile also produces a highly resistant spore that can survive in the environment for a long time and is a requirement for transmission of the bacteria between patients. We have applied our genetic tool to identify all of the genes required for production of a spore. All of these genes represent attractive targets for new drugs to treat infection.

Novel mechanisms of controlling the activities of the transcription factors Spo0A and ComA by the plasmid-encoded quorum sensing regulators Rap60-Phr60 in Bacillus subtilis

Bacillus subtilis and its closest relatives have multiple rap-phr quorum sensing gene pairs that coordinate a variety of physiological processes with population density. Extra-chromosomal rap-phr genes are also present on mobile genetic elements, yet relatively little is known about their function. In this work, we demonstrate that Rap60-Phr60 from plasmid pTA1060 coordinates a variety of biological processes with population density including sporulation, cannibalism, biofilm formation and genetic competence. Similar to other Rap proteins that control sporulation, Rap60 modulates phosphorylation of the transcription factor Spo0A by acting as a phosphatase of Spo0F∼P, an intermediate of the sporulation phosphorelay system. Additionally, Rap60 plays a noncanonical role in regulating the autophosphorylation of the sporulation-specific kinase KinA, a novel activity for Rap proteins. In contrast, Rap proteins that modulate genetic competence interfere with DNA binding by the transcription factor ComA. Rap60 regulates the activity of ComA in a unique manner by forming a Rap60-ComA-DNA ternary complex that inhibits transcription of target genes. Taken together, this work provides new insight into two novel mechanisms of regulating Spo0A and ComA by Rap60 and expands our general understanding of how plasmid-encoded quorum sensing pairs regulate important biological processes.

Transcriptome analysis of cyst formation in Rhodospirillum centenum reveals large global changes in expression during cyst development

Background

Rhodospirillum centenum is a photosynthetic member of the Gram-negative Azospirillum clade members of which exhibit a complex developmental life-cycle featuring morphologically distinct cell types. Under periods of nutrient deprivation, replicative vegetative cells differentiate into metabolically dormant cysts that survive harsh environmental stresses such as desiccation. Encystment involves a multi-stage developmental process that includes the rounding of cells, production of large intracellular storage granules of poly-hydroxybutyrate (PHB) and the excretion of a protective exopolysaccharide coating that envelops dormant cysts.

Results

To study the process of cyst development, we performed RNA-seq studies on cells that were induced to undergo cyst development. To assay for temporal changes in gene expression, RNA was extracted at 4, 24, 48, 72, 96 hours during development and subjected to deep sequence analysis. These results show that 812 genes exhibit log₂ ≥ 1.5-fold changes in expression over a 96 hour cyst induction period demonstrating large global changes in gene expression during cyst development.

Conclusions

Notable changes in expression occurred in numerous genes involved in cell wall and lipid biosynthesis, metabolic enzymes, and numerous regulatory genes such as histidine kinases and transcription factors. Many genes involved in protein synthesis and DNA replication were also significantly reduced during late stages of cyst development. Genes previously identified by genetic screens as being critical for cyst development also exhibited changes of expression during cyst induction. This study provides the first transcriptome profile of global changes in gene expression that occur during development of cysts in a Gram-negative species.

Electronic supplementary material

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

Physiological and cell morphology adaptation of Bacillus subtilis at near-zero specific growth rates: a transcriptome analysis

Nutrient scarcity is a common condition in nature, but the resulting extremely low growth rates (below 0.025 h(-1) ) are an unexplored research area in Bacillus subtilis. To understand microbial life in natural environments, studying the adaptation of B. subtilis to near-zero growth conditions is relevant. To this end, a chemostat modified for culturing an asporogenous B. subtilis sigF mutant strain at extremely low growth rates (also named a retentostat) was set up, and biomass accumulation, culture viability, metabolite production and cell morphology were analysed. During retentostat culturing, the specific growth rate decreased to a minimum of 0.00006 h(-1) , corresponding to a doubling time of 470 days. The energy distribution between growth and maintenance-related processes showed that a state of near-zero growth was reached. Remarkably, a filamentous cell morphology emerged, suggesting that cell separation is impaired under near-zero growth conditions. To evaluate the corresponding molecular adaptations to extremely low specific growth, transcriptome changes were analysed. These revealed that cellular responses to near-zero growth conditions share several similarities with those of cells during the stationary phase of batch growth. However, fundamental differences between these two non-growing states are apparent by their high viability and absence of stationary phase mutagenesis under near-zero growth conditions.

Computational modelling and analysis of the molecular network regulating sporulation initiation in Bacillus subtilis

BACKGROUND

Bacterial spores are important contaminants in food, and the spore forming bacteria are often implicated in food safety and food quality considerations. Spore formation is a complex developmental process involving the expression of more than 500 genes over the course of 6 to 8 hrs. The process culminates in the formation of resting cells capable of resisting environmental extremes and remaining dormant for long periods of time, germinating when conditions promote further vegetative growth. Experimental observations of sporulation and germination are problematic and time consuming so that reliable models are an invaluable asset in terms of prediction and risk assessment. In this report we develop a model which assists in the interpretation of sporulation dynamics.

RESULTS

This paper defines and analyses a mathematical model for the network regulating Bacillus subtilis sporulation initiation, from sensing of sporulation signals down to the activation of the early genes under control of the master regulator Spo0A. Our model summarises and extends other published modelling studies, by allowing the user to execute sporulation initiation in a scenario where Isopropyl β-D-1-thiogalactopyranoside (IPTG) is used as an artificial sporulation initiator as well as in modelling the induction of sporulation in wild-type cells. The analysis of the model results and the comparison with experimental data indicate that the model is good at predicting inducible responses to sporulation signals. However, the model is unable to reproduce experimentally observed accumulation of phosphorelay sporulation proteins in wild type B. subtilis. This model also highlights that the phosphorelay sub-component, which relays the signals detected by the sensor kinases to the master regulator Spo0A, is crucial in determining the response dynamics of the system.

CONCLUSION

We show that there is a complex connectivity between the phosphorelay features and the master regulatory Spo0A. Additional we discovered that the experimentally observed regulation of the phosphotransferase Spo0B for wild-type B. subtilis may be playing an important role in the network which suggests that modelling of sporulation initiation may require additional experimental support.

A mother cell-to-forespore channel: current understanding and future challenges

Formation of endospores allows some bacteria to survive extreme nutrient limitation. The resulting dormant cell, the spore, persists in the environment and is highly resistant to physical and chemical stresses. During spore formation, cells divide asymmetrically and the mother cell engulfs the developing spore, encasing it within a double membrane and isolating it from the medium. Communication between mother cell and isolated forespore involves a specialised connection system that allows nurturing of the forespore and continued macromolecular synthesis, required to finalise spore maturation. Here, we review current understanding of this feeding channel formed by a forespore protein, SpoIIQ, and a mother cell protein, SpoIIIAH, in the model organism Bacillus subtilis and the important human pathogen Clostridium difficile. We also analyse the presence of this channel across endospore-forming bacteria and highlight the main questions still remaining.

Transcriptomic analysis of Clostridium thermocellum Populus hydrolysate-tolerant mutant strain shows increased cellular efficiency in response to Populus hydrolysate compared to the wild type strain

Background

The thermophilic, anaerobic bacterium, Clostridium thermocellum is a model organism for consolidated processing due to its efficient fermentation of cellulose. Constituents of dilute acid pretreatment hydrolysate are known to inhibit C. thermocellum and other microorganisms. To evaluate the biological impact of this type of hydrolysate, a transcriptomic analysis of growth in hydrolysate-containing medium was conducted on 17.5% v/v Populus hydrolysate-tolerant mutant (PM) and wild type (WT) strains of C. thermocellum.

Results

In two levels of Populus hydrolysate medium (0% and 10% v/v), the PM showed both gene specific increases and decreases of gene expression compared to the wild-type strain. The PM had increased expression of genes in energy production and conversion, and amino acid transport and metabolism in both standard and 10% v/v Populus hydrolysate media. In particular, expression of the histidine metabolism increased up to 100 fold. In contrast, the PM decreased gene expression in cell division and sporulation (standard medium only), cell defense mechanisms, cell envelope, cell motility, and cellulosome in both media. The PM downregulated inorganic ion transport and metabolism in standard medium but upregulated it in the hydrolysate media when compared to the WT. The WT differentially expressed 1072 genes in response to the hydrolysate medium which included increased transcription of cell defense mechanisms, cell motility, and cellulosome, and decreased expression in cell envelope, amino acid transport and metabolism, inorganic ion transport and metabolism, and lipid metabolism, while the PM only differentially expressed 92 genes. The PM tolerates up to 17.5% v/v Populus hydrolysate and growth in it elicited 489 genes with differential expression, which included increased expression in energy production and conversion, cellulosome production, and inorganic ion transport and metabolism and decreased expression in transcription and cell defense mechanisms.

Conclusion

These results suggest the mechanisms of tolerance for the Populus hydrolysate-tolerant mutant strain of C. thermocellum are based on increased cellular efficiency caused apparently by downregulation of non-critical genes and increasing the expression of genes in energy production and conversion rather than tolerance to specific hydrolysate components. The wild type, conversely, responds to hydrolysate media by down-regulating growth genes and up-regulating stress response genes.

Transcription of the lysine-2,3-aminomutase gene in the kam locus of Bacillus thuringiensis subsp. kurstaki HD73 is controlled by both σ54 and σK factors

Lysine 2,3-aminomutase (KAM; EC 5.4.3.2) catalyzes the interconversion of l-lysine and l-β-lysine. The transcription and regulation of the kam locus, including lysine-2,3-aminomutase-encoding genes, in Bacillus thuringiensis were analyzed in this study. Reverse transcription-PCR (RT-PCR) analysis revealed that this locus forms two operons: yodT (yodT-yodS-yodR-yodQ-yodP-kamR) and kamA (kamA-yokU-yozE). The transcriptional start sites (TSSs) of the kamA gene were determined using 5' rapid amplification of cDNA ends (RACE). A typical -12/-24 σ(54) binding site was identified in the promoter PkamA, which is located upstream of the kamA gene TSS. A β-galactosidase assay showed that PkamA, which directs the transcription of the kamA operon, is controlled by the σ(54) factor and is activated through the σ(54)-dependent transcriptional regulator KamR. The kamA operon is also controlled by σ(K) and regulated by the GerE protein in the late stage of sporulation. kamR and kamA mutants were prepared by homologous recombination to examine the role of the kam locus. The results showed that the sporulation rate in B. thuringiensis HD(ΔkamR) was slightly decreased compared to that in HD73, whereas that in HD(ΔkamA) was similar to that in HD73. This means that other genes regulated by KamR are important for sporulation.

Bacillus subtilis Systems Biology: Applications of -Omics Techniques to the Study of Endospore Formation

Endospore-forming bacteria, with Bacillus subtilis being the prevalent model organism, belong to the phylum Firmicutes. Although the last common ancestor of all Firmicutes is likely to have been an endospore-forming species, not every lineage in the phylum has maintained the ability to produce endospores (hereafter, spores). In 1997, the release of the full genome sequence for B. subtilis strain 168 marked the beginning of the genomic era for the study of spore formation (sporulation). In this original genome sequence, 139 of the 4,100 protein-coding genes were annotated as sporulation genes. By the time a revised genome sequence with updated annotations was published in 2009, that number had increased significantly, especially since transcriptional profiling studies (transcriptomics) led to the identification of several genes expressed under the control of known sporulation transcription factors. Over the past decade, genome sequences for multiple spore-forming species have been released (including several strains in the Bacillus anthracis / Bacillus cereus group and many Clostridium species), and phylogenomic analyses have revealed many conserved sporulation genes. Parallel advances in transcriptomics led to the identification of small untranslated regulatory RNAs (sRNAs), including some that are expressed during sporulation. An extended array of -omics techniques, i.e., techniques designed to probe gene function on a genome-wide scale, such as proteomics, metabolomics, and high-throughput protein localization studies, have been implemented in microbiology. Combined with the use of new computational methods for predicting gene function and inferring regulatory relationships on a global scale, these -omics approaches are uncovering novel information about sporulation and a variety of other bacterial cell processes.