From modules to networks: A systems-­level analysis of the bacitracin stress response in <i>Bacillus subtilis</i>

Bacterial resistance against antibiotics often involves multiple mechanisms that are interconnected to ensure robust protection. So far, the knowledge about underlying regulatory features of those resistance networks is sparse, since they can hardly be determined by experimentation alone. Here, we present the first computational approach to elucidate the interplay between multiple resistance modules against a single antibiotic and how regulatory network structure allows the cell to respond to and compensate for perturbations of resistance. Based on the response of Bacillus subtilis toward the cell wall synthesis-inhibiting antibiotic bacitracin, we developed a mathematical model that comprehensively describes the protective effect of two well-studied resistance modules (BceAB and BcrC) on the progression of the lipid II cycle. By integrating experimental measurements of expression levels, the model accurately predicts the efficacy of bacitracin against the B. subtilis wild type as well as mutant strains lacking one or both of the resistance modules. Our study reveals that bacitracin-induced changes in the properties of the lipid II cycle itself control the interplay between the two resistance modules. In particular, variations in the concentrations of UPP, the lipid II cycle intermediate that is targeted by bacitracin, connect the effect of the BceAB transporter and the homeostatic response via BcrC to an overall resistance response. We propose that monitoring changes in pathway properties caused by a stressor allows the cell to fine-tune deployment of multiple resistance systems and may serve as a cost-beneficial strategy to control the overall response toward this stressor

Epistasis, core-genome disharmony, and adaptation in recombining bacteria

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome’s co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background—facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche

Regulation of heterologous subtilin production in Bacillus subtilis W168

Background: Subtilin is a peptide antibiotic (lantibiotic) natively produced by Bacillus subtilis ATCC6633. It is encoded in a gene cluster spaBTCSIFEGRK (spa-locus) consisting of four transcriptional units: spaS (subtilin pre-peptide), spaBTC (modification and export), spaIFEG (immunity) and spaRK (regulation). Despite the pioneer understanding on subtilin biosynthesis, a robust platform to facilitate subtilin research and improve subtilin production is still a poorly explored spot. Results: In this work, the intact spa-locus was successfully integrated into the chromosome of Bacillus subtilis W168, which is the by far best-characterized Gram-positive model organism with powerful genetics and many advantages in industrial use. Through systematic analysis of spa-promoter activities in B. subtilis W168 wild type and mutant strains, our work demonstrates that subtilin is basally expressed in B. subtilis W168, and the transition state regulator AbrB strongly represses subtilin biosynthesis in a growth phase-dependent manner. The deletion of AbrB remarkably enhanced subtilin gene expression, resulting in comparable yield of bioactive subtilin production as for B. subtilis ATCC6633. However, while in B. subtilis ATCC6633 AbrB regulates subtilin gene expression via SigH, which in turn activates spaRK, AbrB of B. subtilis W168 controls subtilin gene expression in SigH-independent manner, except for the regulation of spaBTC. Furthermore, the work shows that subtilin biosynthesis in B. subtilis W168 is regulated by the two-component regulatory system SpaRK and strictly relies on subtilin itself as inducer to fulfill the autoregulatory circuit. In addition, by incorporating the subtilin-producing system (spa-locus) and subtilin-reporting system (PpsdA-lux) together, we developed “online” reporter strains to efficiently monitor the dynamics of subtilin biosynthesis. Conclusions: Within this study, the model organism B. subtilis W168 was successfully established as a novel platform for subtilin biosynthesis and the underlying regulatory mechanism was comprehensively characterized. This work will not only facilitate genetic (engineering) studies on subtilin, but also pave the way for its industrial production. More broadly, this work will shed new light on the heterologous production of other lantibiotics