Multiomics Study of Bacterial Growth Arrest in a Synthetic Biology Application

Penalized likelihood optimization for censored missing value imputation in proteomics

Label-free bottom-up proteomics using mass spectrometry and liquid chromatography has long been established as one of the most popular high-throughput analysis workflows for proteome characterization. However, it produces data hindered by complex and heterogeneous missing values, which imputation has long remained problematic. To cope with this, we introduce Pirat, an algorithm that harnesses this challenge using an original likelihood maximization strategy. Notably, it models the instrument limit by learning a global censoring mechanism from the data available. Moreover, it estimates the covariance matrix between enzymatic cleavage products (ie peptides or precursor ions), while offering a natural way to integrate complementary transcriptomic information when multi-omic assays are available. Our benchmarking on several datasets covering a variety of experimental designs (number of samples, acquisition mode, missingness patterns, etc.) and using a variety of metrics (differential analysis ground truth or imputation errors) shows that Pirat outperforms all pre-existing imputation methods. Beyond the interest of Pirat as an imputation tool, these results pinpoint the need for a paradigm change in proteomics imputation, as most pre-existing strategies could be boosted by incorporating similar models to account for the instrument censorship or for the correlation structures, either grounded to the analytical pipeline or arising from a multi-omic approach.

Improving the genetic stability of bacterial growth control for long-term bioproduction

Using microorganisms for bioproduction requires the reorientation of metabolic fluxes from biomass synthesis to the production of compounds of interest. We previously engineered a synthetic growth switch in Escherichia coli based on inducible expression of the β- and β'-subunits of RNA polymerase. Depending on the level of induction, the cells stop growing or grow at a rate close to that of the wild-type strain. This strategy has been successful in transforming growth-arrested bacteria into biofactories with a high production yield, releasing cellular resources from growth towards biosynthesis. However, high selection pressure is placed on a growth-arrested population, favoring mutations that allow cells to escape from growth control. Accordingly, we made the design of the growth switch more robust by building in genetic redundancy. More specifically, we added the rpoA gene, encoding for the α-subunit of RNA polymerase, under the control of a copy of the same inducible promoter used for expression control of ββ'. The improved growth switch is much more stable (escape frequency <10-9), while preserving the capacity to improve production yields. Moreover, after a long period of growth inhibition the population can be regenerated within a few generations. This opens up the possibility to alternate biomass accumulation and product synthesis over a longer period of time and is an additional step towards the dynamical control of bioproduction.

Phenotypic heterogeneity follows a growth-viability tradeoff in response to amino acid identity

In their natural environments, microorganisms mainly operate at suboptimal growth conditions with fluctuations in nutrient abundance. The resulting cellular adaptation is subject to conflicting tasks: growth or survival maximisation. Here, we study this adaptation by systematically measuring the impact of a nitrogen downshift to 24 nitrogen sources on cellular metabolism at the single-cell level. Saccharomyces lineages grown in rich media and exposed to a nitrogen downshift gradually differentiate to form two subpopulations of different cell sizes where one favours growth while the other favours viability with an extended chronological lifespan. This differentiation is asymmetrical with daughter cells representing the new differentiated state with increased viability. We characterise the metabolic response of the subpopulations using RNA sequencing, metabolic biosensors and a transcription factor-tagged GFP library coupled to high-throughput microscopy, imaging more than 800,000 cells. We find that the subpopulation with increased viability is associated with a dormant quiescent state displaying differences in MAPK signalling. Depending on the identity of the nitrogen source present, differentiation into the quiescent state can be actively maintained, attenuated, or aborted. These results establish amino acids as important signalling molecules for the formation of genetically identical subpopulations, involved in chronological lifespan and growth rate determination.

Coupling and uncoupling growth and product formation for producing chemicals

Microbial fermentation employs two strategies: growth- and nongrowth-coupled productions. Stoichiometric metabolic models with flux balance analysis enable pathway engineering to couple target synthesis with growth, yielding numerous successful results. Growth-coupled engineering also contributes to improving bottleneck flux through subsequent adaptive evolution. However, because growth-coupled production inevitably shares resources between biomass and target syntheses, the cost-effective production of bulk chemicals mandates a nongrowth-coupled approach. In such processes, understanding how and when to transition the metabolic state from growth to production modes becomes crucial, as does maintaining cellular activity during the nongrowing state to achieve high productivity. In this paper, we review recent technologies for growth-coupled and nongrowth-coupled production, considering their advantages and disadvantages.

Controlling for false discoveries subsequently to large scale one-way ANOVA testing in proteomics: Practical considerations

In discovery proteomics, as well as many other "omic" approaches, the possibility to test for the differential abundance of hundreds (or of thousands) of features simultaneously is appealing, despite requiring specific statistical safeguards, among which controlling for the false discovery rate (FDR) has become standard. Moreover, when more than two biological conditions or group treatments are considered, it has become customary to rely on the one-way analysis of variance (ANOVA) framework, where a first global differential abundance landscape provided by an omnibus test can be subsequently refined using various post-hoc tests (PHTs). However, the interactions between the FDR control procedures and the PHTs are complex, because both correspond to different types of multiple test corrections (MTCs). This article surveys various ways to orchestrate them in a data processing workflow and discusses their pros and cons.

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

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.