Leveraging high-resolution omics data for predicting responses and adverse events to immune checkpoint inhibitors.

A long-standing goal of personalized and precision medicine is to enable accurate prediction of the outcomes of a given treatment regimen for patients harboring a disease. Currently, many clinical trials fail to meet their endpoints due to underlying factors in the patient population that contribute to either poor responses to the drug of interest or to treatment-related adverse events. Identifying these factors beforehand and correcting for them can lead to an increased success of clinical trials. Comprehensive and large-scale data gathering efforts in biomedicine by omics profiling of the healthy and diseased individuals has led to a treasure-trove of host, disease and environmental factors that contribute to the effectiveness of drugs aiming to treat disease. With increasing omics data, artificial intelligence allows an in-depth analysis of big data and offers a wide range of applications for real-world clinical use, including improved patient selection and identification of actionable targets for companion therapeutics for improved translatability across more patients. As a blueprint for complex drug-disease-host interactions, we here discuss the challenges of utilizing omics data for predicting responses and adverse events in cancer immunotherapy with immune checkpoint inhibitors (ICIs). The omics-based methodologies for improving patient outcomes as in the ICI case have also been applied across a wide-range of complex disease settings, exemplifying the use of omics for in-depth disease profiling and clinical use.

An automated workflow for multi-omics screening of microbial model organisms.

Multi-omics datasets are becoming of key importance to drive discovery in fundamental research as much as generating knowledge for applied biotechnology. However, the construction of such large datasets is usually time-consuming and expensive. Automation might enable to overcome these issues by streamlining workflows from sample generation to data analysis. Here, we describe the construction of a complex workflow for the generation of high-throughput microbial multi-omics datasets. The workflow comprises a custom-built platform for automated cultivation and sampling of microbes, sample preparation protocols, analytical methods for sample analysis and automated scripts for raw data processing. We demonstrate possibilities and limitations of such workflow in generating data for three biotechnologically relevant model organisms, namely Escherichia coli, Saccharomyces cerevisiae, and Pseudomonas putida.

Oral supplementation of nicotinamide riboside alters intestinal microbial composition in rats and mice, but not humans.

The gut microbiota impacts systemic levels of multiple metabolites including NAD+ precursors through diverse pathways. Nicotinamide riboside (NR) is an NAD+ precursor capable of regulating mammalian cellular metabolism. Some bacterial families express the NR-specific transporter, PnuC. We hypothesized that dietary NR supplementation would modify the gut microbiota across intestinal sections. We determined the effects of 12 weeks of NR supplementation on the microbiota composition of intestinal segments of high-fat diet-fed (HFD) rats. We also explored the effects of 12 weeks of NR supplementation on the gut microbiota in humans and mice. In rats, NR reduced fat mass and tended to decrease body weight. Interestingly, NR increased fat and energy absorption but only in HFD-fed rats. Moreover, 16S rRNA gene sequencing analysis of intestinal and fecal samples revealed an increased abundance of species within Erysipelotrichaceae and Ruminococcaceae families in response to NR. PnuC-positive bacterial strains within these families showed an increased growth rate when supplemented with NR. The abundance of species within the Lachnospiraceae family decreased in response to HFD irrespective of NR. Alpha and beta diversity and bacterial composition of the human fecal microbiota were unaltered by NR, but in mice, the fecal abundance of species within Lachnospiraceae increased while abundances of Parasutterella and Bacteroides dorei species decreased in response to NR. In conclusion, oral NR altered the gut microbiota in rats and mice, but not in humans. In addition, NR attenuated body fat mass gain in rats, and increased fat and energy absorption in the HFD context.

Laboratory evolution reveals general and specific tolerance mechanisms for commodity chemicals.

Although strain tolerance to high product concentrations is a barrier to the economically viable biomanufacturing of industrial chemicals, chemical tolerance mechanisms are often unknown. To reveal tolerance mechanisms, an automated platform was utilized to evolve Escherichia coli to grow optimally in the presence of 11 industrial chemicals (1,2-propanediol, 2,3-butanediol, glutarate, adipate, putrescine, hexamethylenediamine, butanol, isobutyrate, coumarate, octanoate, hexanoate), reaching tolerance at concentrations 60%-400% higher than initial toxic levels. Sequencing genomes of 223 isolates from 89 populations, reverse engineering, and cross-compound tolerance profiling were employed to uncover tolerance mechanisms. We show that: 1) cells are tolerized via frequent mutation of membrane transporters or cell wall-associated proteins (e.g., ProV, KgtP, SapB, NagA, NagC, MreB), transcription and translation machineries (e.g., RpoA, RpoB, RpoC, RpsA, RpsG, NusA, Rho), stress signaling proteins (e.g., RelA, SspA, SpoT, YobF), and for certain chemicals, regulators and enzymes in metabolism (e.g., MetJ, NadR, GudD, PurT); 2) osmotic stress plays a significant role in tolerance when chemical concentrations exceed a general threshold and mutated genes frequently overlap with those enabling chemical tolerance in membrane transporters and cell wall-associated proteins; 3) tolerization to a specific chemical generally improves tolerance to structurally similar compounds whereas a tradeoff can occur on dissimilar chemicals, and 4) using pre-tolerized starting isolates can hugely enhance the subsequent production of chemicals when a production pathway is inserted in many, but not all, evolved tolerized host strains, underpinning the need for evolving multiple parallel populations. Taken as a whole, this study provides a comprehensive genotype-phenotype map based on identified mutations and growth phenotypes for 223 chemical tolerant isolates.

Cell-free fetal DNA for genetic evaluation in Copenhagen Pregnancy Loss Study (COPL): a prospective cohort study.

BACKGROUND: One in four pregnancies end in a pregnancy loss. Although the effect on couples is well documented, evidence-based treatments and prediction models are absent. Fetal aneuploidy is associated with a higher chance of a next successful pregnancy compared with euploid pregnancy loss in which underlying maternal conditions might be causal. Ploidy diagnostics are therefore advantageous but challenging as they require collection of the pregnancy tissue. Cell-free fetal DNA (cffDNA) from maternal blood has the potential for evaluation of fetal ploidy status, but no large-scale validation of the method has been done. METHODS: In this prospective cohort study, women with a pregnancy loss were recruited as a part of the Copenhagen Pregnancy Loss (COPL) study from three gynaecological clinics at public hospitals in Denmark. Women were eligible for inclusion if older than 18 years with a pregnancy loss before gestational age 22 weeks (ie, 154 days) and with an intrauterine pregnancy confirmed by ultrasound (including anembryonic sac), and women with pregnancies of unknown location or molar pregnancies were excluded. Maternal blood was collected while pregnancy tissue was still in situ or within 24 h after pregnancy tissue had passed and was analysed by genome-wide sequencing of cffDNA. Direct sequencing of the pregnancy tissue was done as reference. FINDINGS: We included 1000 consecutive women, at the time of a pregnancy loss diagnosis, between Nov 12, 2020, and May 1, 2022. Results from the first 333 women with a pregnancy loss (recruited between Nov 12, 2020, and Aug 14, 2021) were used to evaluate the validity of cffDNA-based testing. Results from the other 667 women were included to evaluate cffDNA performance and result distribution in a larger cohort of 1000 women in total. Gestational age of fetus ranged from 35-149 days (mean of 70·5 days [SD 16·5], or 10 weeks plus 1 day). The cffDNA-based test had a sensitivity for aneuploidy detection of 85% (95% CI 79-90) and a specificity of 93% (95% CI 88-96) compared with direct sequencing of the pregnancy tissue. Among 1000 cffDNA-based test results, 446 (45%) were euploid, 405 (41%) aneuploid, 37 (4%) had multiple aneuploidies, and 112 (11%) were inconclusive. 105 (32%) of 333 women either did not manage to collect the pregnancy tissue or collected a sample classified as unknown tissue giving a high risk of being maternal. INTERPRETATION: This validation of cffDNA-based testing in pregnancy loss shows the potential and feasibility of the method to distinguish euploid and aneuploid pregnancy loss for improved clinical management and benefit of future reproductive medicine and women's health research. FUNDING: Ole Kirks Foundation, BioInnovation Institute Foundation, and the Novo Nordisk Foundation.

Identification and Engineering of Transporters for Efficient Melatonin Production in Escherichia coli.

Transporter discovery and engineering play an important role in cell factory development. Decreasing the intracellular concentration of the product reduces product inhibition and/or toxicity. Lowering intracellular concentrations is especially beneficial for achieving a robust strain at high titers. However, the identification of transporters for xenobiotic chemicals in the host strain is challenging. Here we present a high-throughput workflow to discover Escherichia coli transporters responsible for the efflux of the inhibitory xenobiotic compound melatonin. We took advantage of the Keio collection and screened about 400 transporter knockouts in the presence of a high concentration of melatonin. We found five transporters that when knocked out showed decreased tolerance to melatonin, indicating they are exporters of melatonin. We overexpressed these five genes individually in the production strain and found that one of them, yhjV, encoding a transporter with unknown substrates, resulted in a 27% titer increase in cultivation mimicking fed-batch fermentation. This study demonstrates how microbial cell factories can be improved through transporter identification and engineering. Further, these results lay the foundation for the scale-up of melatonin production in E. coli.

Membrane transporter identification and modulation via adaptive laboratory evolution.

Membrane transport proteins are potential targets for medical and biotechnological applications. However, more than 30% of reported membrane transporter families are either poorly characterized or lack adequate functional annotation. Here, adaptive laboratory evolution was leveraged to identify membrane transporters for a set of four amino acids as well as specific mutations that modulate the activities of these transporters. Specifically, Escherichia coli was adaptively evolved under increasing concentrations of L-histidine, L-phenylalanine, L-threonine, and L-methionine separately with multiple replicate evolutions. Evolved populations and isolated clones displayed growth rates comparable to the unstressed ancestral strain at elevated concentrations (four-to six-fold increases) of the targeted amino acids. Whole genome sequencing of the evolved strains revealed a diverse number of key mutations, including SNPs, small deletions, and copy number variants targeting the transporters leuE for histidine, yddG for phenylalanine, yedA for methionine, and brnQ and rhtC for threonine. Reverse engineering of the mutations in the ancestral strain established mutation causality of the specific mutations for the tolerant phenotypes. The functional roles of yedA and brnQ in the transport of methionine and threonine, respectively, are novel assignments and their functional roles were validated using a flow cytometry cellular accumulation assay. To demonstrate how the identified transporters can be leveraged for production, an L-phenylalanine overproduction strain was shown to be a superior producer when the identified yddG exporter was overexpressed. Overall, the results revealed the striking efficiency of laboratory evolution to identify transporters and specific mutational mechanisms to modulate their activities, thereby demonstrating promising applicability in transporter discovery efforts and strain engineering.

A dual-reporter system for investigating and optimizing protein translation and folding in E. coli.

Strategies for investigating and optimizing the expression and folding of proteins for biotechnological and pharmaceutical purposes are in high demand. Here, we describe a dual-reporter biosensor system that simultaneously assesses in vivo protein translation and protein folding, thereby enabling rapid screening of mutant libraries. We have validated the dual-reporter system on five different proteins and find an excellent correlation between reporter signals and the levels of protein expression and solubility of the proteins. We further demonstrate the applicability of the dual-reporter system as a screening assay for deep mutational scanning experiments. The system enables high throughput selection of protein variants with high expression levels and altered protein stability. Next generation sequencing analysis of the resulting libraries of protein variants show a good correlation between computationally predicted and experimentally determined protein stabilities. We furthermore show that the mutational experimental data obtained using this system may be useful for protein structure calculations.

Genome-scale metabolic modeling of P. thermoglucosidasius NCIMB 11955 reveals metabolic bottlenecks in anaerobic metabolism.

Parageobacillus thermoglucosidasius represents a thermophilic, facultative anaerobic bacterial chassis, with several desirable traits for metabolic engineering and industrial production. To further optimize strain productivity, a systems level understanding of its metabolism is needed, which can be facilitated by a genome-scale metabolic model. Here, we present p-thermo, the most complete, curated and validated genome-scale model (to date) of Parageobacillus thermoglucosidasius NCIMB 11955. It spans a total of 890 metabolites, 1175 reactions and 917 metabolic genes, forming an extensive knowledge base for P. thermoglucosidasius NCIMB 11955 metabolism. The model accurately predicts aerobic utilization of 22 carbon sources, and the predictive quality of internal fluxes was validated with previously published 13C-fluxomics data. In an application case, p-thermo was used to facilitate more in-depth analysis of reported metabolic engineering efforts, giving additional insight into fermentative metabolism. Finally, p-thermo was used to resolve a previously uncharacterised bottleneck in anaerobic metabolism, by identifying the minimal required supplemented nutrients (thiamin, biotin and iron(III)) needed to sustain anaerobic growth. This highlights the usefulness of p-thermo for guiding the generation of experimental hypotheses and for facilitating data-driven metabolic engineering, expanding the use of P. thermoglucosidasius as a high yield production platform.

SmartPeak Automates Targeted and Quantitative Metabolomics Data Processing.

Technological advances in high-resolution mass spectrometry (MS) vastly increased the number of samples that can be processed in a life science experiment, as well as volume and complexity of the generated data. To address the bottleneck of high-throughput data processing, we present SmartPeak (https://github.com/AutoFlowResearch/SmartPeak), an application that encapsulates advanced algorithms to enable fast, accurate, and automated processing of capillary electrophoresis-, gas chromatography-, and liquid chromatography (LC)-MS(/MS) data and high-pressure LC data for targeted and semitargeted metabolomics, lipidomics, and fluxomics experiments. The application allows for an approximate 100-fold reduction in the data processing time compared to manual processing while enhancing quality and reproducibility of the results.

Dietary Change Enables Robust Growth-Coupling of Heterologous Methyltransferase Activity in Yeast.

Genetic modifications of living organisms and proteins are made possible by a catalogue of molecular and synthetic biology tools, yet proper screening assays for genetic variants of interest continue to lag behind. Synthetic growth-coupling (GC) of enzyme activities offers a simple, inexpensive way to track such improvements. In this follow-up study we present the optimization of a recently established GC design for screening of heterologous methyltransferases (MTases) and related pathways in the yeast Saccharomyces cerevisiae. Specifically, upon testing different media compositions and genetic backgrounds, improved GC of different heterologous MTase activities is obtained. Furthermore, we demonstrate the strength of the system by screening a library of catechol O-MTase variants converting protocatechuic acid into vanillic acid. We demonstrated high correlation (R2 = 0.775) between vanillic acid and cell density as a proxy for MTase activity. We envision that the improved MTase GC can aid evolution-guided optimization of biobased production processes for methylated compounds with yeast in the future.

Elucidating aromatic acid tolerance at low pH in Saccharomyces cerevisiae using adaptive laboratory evolution.

Toxicity from the external presence or internal production of compounds can reduce the growth and viability of microbial cell factories and compromise productivity. Aromatic compounds are generally toxic for microorganisms, which makes their production in microbial hosts challenging. Here we use adaptive laboratory evolution to generate Saccharomyces cerevisiae mutants tolerant to two aromatic acids, coumaric acid and ferulic acid. The evolution experiments were performed at low pH (3.5) to reproduce conditions typical of industrial processes. Mutant strains tolerant to levels of aromatic acids near the solubility limit were then analyzed by whole genome sequencing, which revealed prevalent point mutations in a transcriptional activator (Aro80) that is responsible for regulating the use of aromatic amino acids as the nitrogen source. Among the genes regulated by Aro80, ESBP6 was found to be responsible for increasing tolerance to aromatic acids by exporting them out of the cell. Further examination of the native function of Esbp6 revealed that this transporter can excrete fusel acids (byproducts of aromatic amino acid catabolism) and this role is shared with at least one additional transporter native to S. cerevisiae (Pdr12). Besides conferring tolerance to aromatic acids, ESBP6 overexpression was also shown to significantly improve the secretion in coumaric acid production strains. Overall, we showed that regulating the activity of transporters is a major mechanism to improve tolerance to aromatic acids. These findings can be used to modulate the intracellular concentration of aromatic compounds to optimize the excretion of such products while keeping precursor molecules inside the cell.

Adaptive laboratory evolution of Pseudomonas putida KT2440 improves p-coumaric and ferulic acid catabolism and tolerance.

Pseudomonas putida KT2440 is a promising bacterial chassis for the conversion of lignin-derived aromatic compound mixtures to biofuels and bioproducts. Despite the inherent robustness of this strain, further improvements to aromatic catabolism and toxicity tolerance of P. putida will be required to achieve industrial relevance. Here, tolerance adaptive laboratory evolution (TALE) was employed with increasing concentrations of the hydroxycinnamic acids p-coumaric acid (pCA) and ferulic acid (FA) individually and in combination (pCA â€‹+ â€‹FA). The TALE experiments led to evolved P. putida strains with increased tolerance to the targeted acids as compared to wild type. Specifically, a 37 â€‹h decrease in lag phase in 20 â€‹g/L pCA and a 2.4-fold increase in growth rate in 30 â€‹g/L FA was observed. Whole genome sequencing of intermediate and endpoint evolved P. putida populations revealed several expected and non-intuitive genetic targets underlying these aromatic catabolic and toxicity tolerance enhancements. PP_3350 and ttgB were among the most frequently mutated genes, and the beneficial contributions of these mutations were verified via gene knockouts. Deletion of PP_3350, encoding a hypothetical protein, recapitulated improved toxicity tolerance to high concentrations of pCA, but not an improved growth rate in high concentrations of FA. Deletion of ttgB, part of the TtgABC efflux pump, severely inhibited growth in pCA â€‹+ â€‹FA TALE-derived strains but did not affect growth in pCA â€‹+ â€‹FA in a wild type background, suggesting epistatic interactions. Genes involved in flagellar movement and transcriptional regulation were often mutated in the TALE experiments on multiple substrates, reinforcing ideas of a minimal and deregulated cell as optimal for domesticated growth. Overall, this work demonstrates increased tolerance towards and growth rate at the expense of hydroxycinnamic acids and presents new targets for improving P. putida for microbial lignin valorization.

CRISPR interference of nucleotide biosynthesis improves production of a single-domain antibody in Escherichia coli.

Growth decoupling can be used to optimize the production of biochemicals and proteins in cell factories. Inhibition of excess biomass formation allows for carbon to be utilized efficiently for product formation instead of growth, resulting in increased product yields and titers. Here, we used CRISPR interference to increase the production of a single-domain antibody (sdAb) by inhibiting growth during production. First, we screened 21 sgRNA targets in the purine and pyrimidine biosynthesis pathways and found that the repression of 11 pathway genes led to the increased green fluorescent protein production and decreased growth. The sgRNA targets pyrF, pyrG, and cmk were selected and further used to improve the production of two versions of an expression-optimized sdAb. Proteomics analysis of the sdAb-producing pyrF, pyrG, and cmk growth decoupling strains showed significantly decreased RpoS levels and an increase of ribosome-associated proteins, indicating that the growth decoupling strains do not enter stationary phase and maintain their capacity for protein synthesis upon growth inhibition. Finally, sdAb production was scaled up to shake-flask fermentation where the product yield was improved 2.6-fold compared to the control strain with no sgRNA target sequence. An sdAb content of 14.6% was reached in the best-performing pyrG growth decoupling strain.

Combining Environmental and Economic Performance for Bioprocess Optimization.

Biochemical production faces economic and environmental challenges that need to be overcome to enable a viable and sustainable bioeconomy. We propose an assessment framework that consistently combines environmental and economic indicators to support optimized biochemical production at early development stages. We define internally consistent system boundaries and a comprehensive set of quantitative indicators from life cycle assessment (LCA) and techno-economic assessment (TEA) to combine environmental and economic performance in a single score. Our framework enables the identification of trade-offs across environmental and economic aspects over the entire biochemical life cycle. This approach provides input for the optimization of future biochemicals in terms of overall sustainability, to overcome prevailing obstacles in the development of biochemical production processes.

Fluctuations in glucose availability prevent global proteome changes and physiological transition during prolonged chemostat cultivations of Saccharomyces cerevisiae.

Chemostat cultivation mode imposes selective pressure on the cells, which may result in slow adaptation in the physiological state over time. We applied a two-compartment scale-down chemostat system imposing feast-famine conditions to characterize the long-term (100 s of hours) response of Saccharomyces cerevisiae to fluctuating glucose availability. A wild-type strain and a recombinant strain, expressing an insulin precursor, were cultured in the scale-down system, and analyzed at the physiological and proteomic level. Phenotypes of both strains were compared with those observed in a well-mixed chemostat. Our results show that S. cerevisiae subjected to long-term chemostat conditions undergoes a global reproducible shift in its cellular state and that this transition occurs faster and is larger in magnitude for the recombinant strain including a significant decrease in the expression of the insulin product. We find that the transition can be completely avoided in the presence of fluctuations in glucose availability as the strains subjected to feast-famine conditions under otherwise constant culture conditions exhibited constant levels of the measured proteome for over 250 hr. We hypothesize possible mechanisms responsible for the observed phenotypes and suggest experiments that could be used to test these mechanisms.

Publisher Correction: MEMOTE for standardized genome-scale metabolic model testing.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.