The Multiomics Response of Bacillus subtilis to Simultaneous Genetic and Environmental Perturbations

Abstract

How bacteria respond at the systems level to both genetic and environmental perturbations imposed at the same time is one fundamental yet open question in biology. Bioengineering or synthetic biology provides an ideal system for studying such responses, as engineered strains always have genetic changes as opposed to wildtypes and are grown in conditions which often change during growth for maximal yield of desired products. So, engineered strains were used to address the outstanding question. Two Bacillus subtilis strains (MT1 and MT2) were created previously for the overproduction of N-acetylglucosamine (GlcNAc), which were grown in an environment with a carbon shift from glucose to glucose and xylose in the same culture system. We had four groups: (1) a wildtype (WT) grown with glucose at t1; (2) a WT with glucose and xylose at t2; (3) a mutant (MT1) grown with glucose at t1; and (4) MT1 with glucose and xylose at t2. By measuring transcriptomes and metabolomes, we found that GlcNAc-producing mutants, particularly MT2, had a higher yield of N-acetylglucosamine than WT but displayed a smaller maximum growth rate than the wildtype, despite MT1 reaching higher carrying capacity. Underlying the observed growth, the engineered pathways leading to N-acetylglucosamine had both higher gene expression and associated metabolite concentrations in MT1 than WT at both t1 and t2; in bioenergetics, there was higher energy supply in terms of ATP and GTP, with the energy state metric higher in MT1 than WT at both timepoints. Additionally, most top key precursor metabolites were equally abundant in MT1 and WT at either timepoints. Besides that, one prominent feature was the high consistency between transcriptomics and metabolomics in revealing the response. First, both metabolomes and transcriptomes revealed the same PCA clusters of the four groups. Second, we found that the important functions enriched both by metabolomes and transcriptomes overlapped, such as amino acid metabolism and ABC transport. Strikingly, these functions overlapped those enriched by the genes showing a high (positive or negative) correlation with metabolites. Furthermore, these functions also overlapped the enriched KEGG pathways identified using weighted gene coexpression network analysis. All these findings suggest that the responses to simultaneous genetic and environmental perturbations are well coordinated at the metabolic and transcriptional levels: they rely heavily on bioenergetics, but core metabolism does not differ much, while amino acid metabolism and ABC transport are important. This serves as a design guide for bioengineering, synthetic biology, and systems biology.