Machine-learning from Pseudomonas putida KT2440 transcriptomes reveals its transcriptional regulatory network

Current models in bacterial hemicellulase-encoding gene regulation

The discovery and characterization of bacterial carbohydrate-active enzymes is a fundamental component of biotechnology innovation, particularly for renewable fuels and chemicals; however, these studies have increasingly transitioned to exploring the complex regulation required for recalcitrant polysaccharide utilization. This pivot is largely due to the current need to engineer and optimize enzymes for maximal degradation in industrial or biomedical applications. Given the structural simplicity of a single cellulose polymer, and the relatively few enzyme classes required for complete bioconversion, the regulation of cellulases in bacteria has been thoroughly discussed in the literature. However, the diversity of hemicelluloses found in plant biomass and the multitude of carbohydrate-active enzymes required for their deconstruction has resulted in a less comprehensive understanding of bacterial hemicellulase-encoding gene regulation. Here we review the mechanisms of this process and common themes found in the transcriptomic response during plant biomass utilization. By comparing regulatory systems from both Gram-negative and Gram-positive bacteria, as well as drawing parallels to cellulase regulation, our goals are to highlight the shared and distinct features of bacterial hemicellulase-encoding gene regulation and provide a set of guiding questions to improve our understanding of bacterial lignocellulose utilization. KEY POINTS: • Canonical regulatory mechanisms for bacterial hemicellulase-encoding gene expression include hybrid two-component systems (HTCS), extracytoplasmic function (ECF)-σ/anti-σ systems, and carbon catabolite repression (CCR). • Current transcriptomic approaches are increasingly being used to identify hemicellulase-encoding gene regulatory patterns coupled with computational predictions for transcriptional regulators. • Future work should emphasize genetic approaches to improve systems biology tools available for model bacterial systems and emerging microbes with biotechnology potential. Specifically, optimization of Gram-positive systems will require integration of degradative and fermentative capabilities, while optimization of Gram-negative systems will require bolstering the potency of lignocellulolytic capabilities.

Meta-analysis Driven Strain Design for Mitigating Oxidative Stresses Important in Biomanufacturing

As the availability of data sets increases, meta-analysis leveraging aggregated and interoperable data types is proving valuable. This study leveraged a meta-analysis workflow to identify mutations that could improve robustness to reactive oxygen species (ROS) stresses using an industrially important melatonin production strain as an example. ROS stresses often occur during cultivation and negatively affect strain performance. Cellular response to ROS is also linked to the SOS response and resistance to pH fluctuations, which is important to strain robustness in large-scale biomanufacturing. This work integrated more than 7000 E. coli adaptive laboratory evolution (ALE) mutations across 59 experiments to statistically associate mutated genes to 2 ROS tolerance ALE conditions from 72 unique conditions. Mutant oxyR, fur, iscR, and ygfZ were significantly associated and hypothesized to contribute fitness in ROS stress. Across these genes, 259 total mutations were inspected in conjunction with transcriptomics from 46 iModulon experiments. Ten mutations were chosen for reintroduction based on mutation clustering and coinciding transcriptional changes as evidence of fitness impact. Strains with mutations reintroduced into oxyR, fur, iscR, and ygfZ exhibited increased tolerance to H2O2 and acid stress and reduced SOS response, all of which are related to ROS. Additionally, new evidence was generated toward understanding the function of ygfZ, an uncharacterized gene. This meta-analysis approach utilized aggregated and interoperable multiomics data sets to identify mutations conferring industrially relevant phenotypes with the least drawbacks, describing an approach for data-driven strain engineering to optimize microbial cell factories.

Evolution and engineering of pathways for aromatic O-demethylation in Pseudomonas putida KT2440

Biological conversion of lignin from biomass offers a promising strategy for sustainable production of fuels and chemicals. However, aromatic compounds derived from lignin commonly contain methoxy groups, and O-demethylation of these substrates is often a rate-limiting reaction that influences catabolic efficiency. Several enzyme families catalyze aromatic O-demethylation, but they are rarely compared in vivo to determine an optimal biocatalytic strategy. Here, two pathways for aromatic O-demethylation were compared in Pseudomonas putida KT2440. The native Rieske non-heme iron monooxygenase (VanAB) and, separately, a heterologous tetrahydrofolate-dependent demethylase (LigM) were constitutively expressed in P. putida, and the strains were optimized via adaptive laboratory evolution (ALE) with vanillate as a model substrate. All evolved strains displayed improved growth phenotypes, with the evolved strains harboring the native VanAB pathway exhibiting growth rates ∼1.8x faster than those harboring the heterologous LigM pathway. Enzyme kinetics and transcriptomics studies investigated the contribution of selected mutations toward enhanced utilization of vanillate. The VanAB-overexpressing strains contained the most impactful mutations, including those in VanB, the reductase for vanillate O-demethylase, PP_3494, a global regulator of vanillate catabolism, and fghA, involved in formaldehyde detoxification. These three mutations were combined into a single strain, which exhibited approximately 5x faster vanillate consumption than the wild-type strain in the first 8 h of cultivation. Overall, this study illuminates the details of vanillate catabolism in the context of two distinct enzymatic mechanisms, yielding a platform strain for efficient O-demethylation of lignin-related aromatic compounds to value-added products.

Accelerating Genetic Sensor Development, Scale-up, and Deployment Using Synthetic Biology

Living cells are exquisitely tuned to sense and respond to changes in their environment. Repurposing these systems to create engineered biosensors has seen growing interest in the field of synthetic biology and provides a foundation for many innovative applications spanning environmental monitoring to improved biobased production. In this review, we present a detailed overview of currently available biosensors and the methods that have supported their development, scale-up, and deployment. We focus on genetic sensors in living cells whose outputs affect gene expression. We find that emerging high-throughput experimental assays and evolutionary approaches combined with advanced bioinformatics and machine learning are establishing pipelines to produce genetic sensors for virtually any small molecule, protein, or nucleic acid. However, more complex sensing tasks based on classifying compositions of many stimuli and the reliable deployment of these systems into real-world settings remain challenges. We suggest that recent advances in our ability to precisely modify nonmodel organisms and the integration of proven control engineering principles (e.g., feedback) into the broader design of genetic sensing systems will be necessary to overcome these hurdles and realize the immense potential of the field.

Proteome allocation is linked to transcriptional regulation through a modularized transcriptome

It has proved challenging to quantitatively relate the proteome to the transcriptome on a per-gene basis. Recent advances in data analytics have enabled a biologically meaningful modularization of the bacterial transcriptome. We thus investigate whether matched datasets of transcriptomes and proteomes from bacteria under diverse conditions can be modularized in the same way to reveal novel relationships between their compositions. We find that; (1) the modules of the proteome and the transcriptome are comprised of a similar list of gene products, (2) the modules in the proteome often represent combinations of modules from the transcriptome, (3) known transcriptional and post-translational regulation is reflected in differences between two sets of modules, allowing for knowledge-mapping when interpreting module functions, and (4) through statistical modeling, absolute proteome allocation can be inferred from the transcriptome alone. Quantitative and knowledge-based relationships can thus be found at the genome-scale between the proteome and transcriptome in bacteria.

Deciphering nutritional stress responses via knowledge-enriched transcriptomics for microbial engineering

Understanding diverse bacterial nutritional requirements and responses is foundational in microbial research and biotechnology. In this study, we employed knowledge-enriched transcriptomic analytics to decipher complex stress responses of Vibrio natriegens to supplied nutrients, aiming to enhance microbial engineering efforts. We computed 64 independently modulated gene sets that comprise a quantitative basis for transcriptome dynamics across a comprehensive transcriptomics dataset containing a broad array of nutrient conditions. Our approach led to the i) identification of novel transporter systems for diverse substrates, ii) a detailed understanding of how trace elements affect metabolism and growth, and iii) extensive characterization of nutrient-induced stress responses, including osmotic stress, low glycolytic flux, proteostasis, and altered protein expression. By clarifying the relationship between the acetate-associated regulon and glycolytic flux status of various nutrients, we have showcased its vital role in directing optimal carbon source selection. Our findings offer deep insights into the transcriptional landscape of bacterial nutrition and underscore its significance in tailoring strain engineering strategies, thereby facilitating the development of more efficient and robust microbial systems for biotechnological applications.

Transcription factor-based biosensor: A molecular-guided approach for advanced biofuel synthesis

As a sustainable and renewable alternative to petroleum fuels, advanced biofuels shoulder the responsibility of energy saving, emission reduction and environmental protection. Traditional engineering of cell factories for production of advanced biofuels lacks efficient high-throughput screening tools and regulating systems, impeding the improvement of cellular productivity and yield. Transcription factor-based biosensors have been widely applied to monitor and regulate microbial cell factory products due to the advantages of fast detection and in-situ screening. This review updates the design and application of transcription factor-based biosensors tailored for advanced biofuels and related intermediates. The construction and genetic parts selection principle of biosensors are discussed. Strategies to enhance the performance of biosensor, including regulating promoter strength and RBS strength, optimizing plasmid copy number, implementing genetic amplifier, and modulating the structure of transcription factor, have also been summarized. We further review the application of biosensors in high-throughput screening of new metabolic engineering targets, evolution engineering, confirmation of protein function, and dynamic regulation of metabolic flux for higher production of advanced biofuels. At last, we discuss the current limitations and future trends of transcription factor-based biosensors.

Machine learning analysis of RB-TnSeq fitness data predicts functional gene modules in Pseudomonas putida KT2440

There is growing interest in engineering Pseudomonas putida KT2440 as a microbial chassis for the conversion of renewable and waste-based feedstocks, and metabolic engineering of P. putida relies on the understanding of the functional relationships between genes. In this work, independent component analysis (ICA) was applied to a compendium of existing fitness data from randomly barcoded transposon insertion sequencing (RB-TnSeq) of P. putida KT2440 grown in 179 unique experimental conditions. ICA identified 84 independent groups of genes, which we call fModules ("functional modules"), where gene members displayed shared functional influence in a specific cellular process. This machine learning-based approach both successfully recapitulated previously characterized functional relationships and established hitherto unknown associations between genes. Selected gene members from fModules for hydroxycinnamate metabolism and stress resistance, acetyl coenzyme A assimilation, and nitrogen metabolism were validated with engineered mutants of P. putida. Additionally, functional gene clusters from ICA of RB-TnSeq data sets were compared with regulatory gene clusters from prior ICA of RNAseq data sets to draw connections between gene regulation and function. Because ICA profiles the functional role of several distinct gene networks simultaneously, it can reduce the time required to annotate gene function relative to manual curation of RB-TnSeq data sets.

IMPORTANCE

This study demonstrates a rapid, automated approach for elucidating functional modules within complex genetic networks. While Pseudomonas putida randomly barcoded transposon insertion sequencing data were used as a proof of concept, this approach is applicable to any organism with existing functional genomics data sets and may serve as a useful tool for many valuable applications, such as guiding metabolic engineering efforts in other microbes or understanding functional relationships between virulence-associated genes in pathogenic microbes. Furthermore, this work demonstrates that comparison of data obtained from independent component analysis of transcriptomics and gene fitness datasets can elucidate regulatory-functional relationships between genes, which may have utility in a variety of applications, such as metabolic modeling, strain engineering, or identification of antimicrobial drug targets.

Advancing the scale of synthetic biology via cross-species transfer of cellular functions enabled by iModulon engraftment

Machine learning applied to large compendia of transcriptomic data has enabled the decomposition of bacterial transcriptomes to identify independently modulated sets of genes, such iModulons represent specific cellular functions. The identification of iModulons enables accurate identification of genes necessary and sufficient for cross-species transfer of cellular functions. We demonstrate cross-species transfer of: 1) the biotransformation of vanillate to protocatechuate, 2) a malonate catabolic pathway, 3) a catabolic pathway for 2,3-butanediol, and 4) an antimicrobial resistance to ampicillin found in multiple Pseudomonas species to Escherichia coli. iModulon-based engineering is a transformative strategy as it includes all genes comprising the transferred cellular function, including genes without functional annotation. Adaptive laboratory evolution was deployed to optimize the cellular function transferred, revealing mutations in the host. Combining big data analytics and laboratory evolution thus enhances the level of understanding of systems biology, and synthetic biology for strain design and development.

Independent component analysis reveals 49 independently modulated gene sets within the global transcriptional regulatory architecture of multidrug-resistant Acinetobacter baumannii

Acinetobacter baumannii causes severe infections in humans, resists multiple antibiotics, and survives in stressful environmental conditions due to modulations of its complex transcriptional regulatory network (TRN). Unfortunately, our global understanding of the TRN in this emerging opportunistic pathogen is limited. Here, we apply independent component analysis, an unsupervised machine learning method, to a compendium of 139 RNA-seq data sets of three multidrug-resistant A. baumannii international clonal complex I strains (AB5075, AYE, and AB0057). This analysis allows us to define 49 independently modulated gene sets, which we call iModulons. Analysis of the identified A. baumannii iModulons reveals validating parallels to previously defined biological operons/regulons and provides a framework for defining unknown regulons. By utilizing the iModulons, we uncover potential mechanisms for a RpoS-independent general stress response, define global stress-virulence trade-offs, and identify conditions that may induce plasmid-borne multidrug resistance. The iModulons provide a model of the TRN that emphasizes the importance of transcriptional regulation of virulence phenotypes in A. baumannii. Furthermore, they suggest the possibility of future interventions to guide gene expression toward diminished pathogenic potential.IMPORTANCEThe rise in hospital outbreaks of multidrug-resistant Acinetobacter baumannii infections underscores the urgent need for alternatives to traditional broad-spectrum antibiotic therapies. The success of A. baumannii as a significant nosocomial pathogen is largely attributed to its ability to resist antibiotics and survive environmental stressors. However, there is limited literature available on the global, complex regulatory circuitry that shapes these phenotypes. Computational tools that can assist in the elucidation of A. baumannii's transcriptional regulatory network architecture can provide much-needed context for a comprehensive understanding of pathogenesis and virulence, as well as for the development of targeted therapies that modulate these pathways.

Recombinant Pseudomonas growing on non-natural fluorinated substrates shows stress but overall tolerance to cytoplasmically released fluoride anion

IMPORTANCE

Society uses thousands of organofluorine compounds, sometimes denoted per- and polyfluoroalkyl substances (PFAS), in hundreds of products, but recent studies have shown some to manifest human and environmental health effects. As a class, they are recalcitrant to biodegradation, partly due to the paucity of fluorinated natural products to which microbes have been exposed. Another limit to PFAS biodegradation is the intracellular toxicity of fluoride anion generated from C-F bond cleavage. The present study identified a broader substrate specificity in an enzyme originally studied for its activity on the natural product fluoroacetate. A recombinant Pseudomonas expressing this enzyme was used here as a model system to better understand the limits and effects of a high level of intracellular fluoride generation. A fluoride stress response has evolved in bacteria and has been described in Pseudomonas spp. The present study is highly relevant to organofluorine compound degradation or engineered biosynthesis in which fluoride anion is a substrate.

Improving microbial bioproduction under low-oxygen conditions

Microbial bioconversion provides access to a wide range of sustainably produced chemicals and commodities. However, industrial-scale bioproduction process operations are preferred to be anaerobic due to the cost associated with oxygen transfer. Anaerobic bioconversion generally offers limited substrate utilization profiles, lower product yields, and reduced final product diversity compared with aerobic processes. Bioproduction under conditions of reduced oxygen can overcome the limitations of fully aerobic and anaerobic bioprocesses, but many microbial hosts are not developed for low-oxygen bioproduction. Here, we describe advances in microbial strain engineering involving the use of redox cofactor engineering, genome-scale metabolic modeling, and functional genomics to enable improved bioproduction processes under low oxygen and provide a viable path for scaling these bioproduction systems to industrial scales.

Engineering of increased L-Threonine production in bacteria by combinatorial cloning and machine learning

The goal of this study is to develop a general strategy for bacterial engineering using an integrated synthetic biology and machine learning (ML) approach. This strategy was developed in the context of increasing L-threonine production in Escherichia coli ATCC 21277. A set of 16 genes was initially selected based on metabolic pathway relevance to threonine biosynthesis and used for combinatorial cloning to construct a set of 385 strains to generate training data (i.e., a range of L-threonine titers linked to each of the specific gene combinations). Hybrid (regression/classification) deep learning (DL) models were developed and used to predict additional gene combinations in subsequent rounds of combinatorial cloning for increased L-threonine production based on the training data. As a result, E. coli strains built after just three rounds of iterative combinatorial cloning and model prediction generated higher L-threonine titers (from 2.7 g/L to 8.4 g/L) than those of patented L-threonine strains being used as controls (4-5 g/L). Interesting combinations of genes in L-threonine production included deletions of the tdh, metL, dapA, and dhaM genes as well as overexpression of the pntAB, ppc, and aspC genes. Mechanistic analysis of the metabolic system constraints for the best performing constructs offers ways to improve the models by adjusting weights for specific gene combinations. Graph theory analysis of pairwise gene modifications and corresponding levels of L-threonine production also suggests additional rules that can be incorporated into future ML models.

Independent component analysis of Corynebacterium glutamicum transcriptomes reveals its transcriptional regulatory network

Gene expression in bacteria is regulated by multiple transcription factors. Clarifying the regulation mechanism of gene expression is necessary to understand bacterial physiological activities. To further understand the structure of the transcriptional regulatory network of Corynebacterium glutamicum, we applied independent component analysis, an unsupervised machine learning algorithm, to the high-quality C. glutamicum gene expression profile which includes 263 samples from 29 independent projects. We obtained 87 robust independent regulatory modules (iModulons). These iModulons explain 76.7% of the variance in the expression profile and constitute the quantitative transcriptional regulatory network of C. glutamicum. By analyzing the constituent genes in iModulons, we identified potential targets for 20 transcription factors. We also captured the changes in iModulon activities under different growth rates and dissolved oxygen concentrations, demonstrating the ability of iModulons to comprehensively interpret transcriptional responses to environmental changes. In summary, this study provides a genome-scale quantitative transcriptional regulatory network for C. glutamicum and informs future research on complex changes in the transcriptome.

A multi-scale expression and regulation knowledge base for Escherichia coli

Transcriptomic data is accumulating rapidly; thus, scalable methods for extracting knowledge from this data are critical. Here, we assembled a top-down expression and regulation knowledge base for Escherichia coli. The expression component is a 1035-sample, high-quality RNA-seq compendium consisting of data generated in our lab using a single experimental protocol. The compendium contains diverse growth conditions, including: 9 media; 39 supplements, including antibiotics; 42 heterologous proteins; and 76 gene knockouts. Using this resource, we elucidated global expression patterns. We used machine learning to extract 201 modules that account for 86% of known regulatory interactions, creating the regulatory component. With these modules, we identified two novel regulons and quantified systems-level regulatory responses. We also integrated 1675 curated, publicly-available transcriptomes into the resource. We demonstrated workflows for analyzing new data against this knowledge base via deconstruction of regulation during aerobic transition. This resource illuminates the E. coli transcriptome at scale and provides a blueprint for top-down transcriptomic analysis of non-model organisms.

Machine learning uncovers the Pseudomonas syringae transcriptome in microbial communities and during infection

IMPORTANCE

Pseudomonas syringae pv. tomato DC3000 is a model plant pathogen that infects tomatoes and Arabidopsis thaliana. The current understanding of global transcriptional regulation in the pathogen is limited. Here, we applied iModulon analysis to a compendium of RNA-seq data to unravel its transcriptional regulatory network. We characterize each co-regulated gene set, revealing the activity of major regulators across diverse conditions. We provide new insights on the transcriptional dynamics in interactions with the plant immune system and with other bacterial species, such as AlgU-dependent regulation of flagellar genes during plant infection and downregulation of siderophore production in the presence of a siderophore cheater. This study demonstrates the novel application of iModulons in studying temporal dynamics during host-pathogen and microbe-microbe interactions, and reveals specific insights of interest.

Laboratory evolution, transcriptomics, and modeling reveal mechanisms of paraquat tolerance

Relationships between the genome, transcriptome, and metabolome underlie all evolved phenotypes. However, it has proved difficult to elucidate these relationships because of the high number of variables measured. A recently developed data analytic method for characterizing the transcriptome can simplify interpretation by grouping genes into independently modulated sets (iModulons). Here, we demonstrate how iModulons reveal deep understanding of the effects of causal mutations and metabolic rewiring. We use adaptive laboratory evolution to generate E. coli strains that tolerate high levels of the redox cycling compound paraquat, which produces reactive oxygen species (ROS). We combine resequencing, iModulons, and metabolic models to elucidate six interacting stress-tolerance mechanisms: (1) modification of transport, (2) activation of ROS stress responses, (3) use of ROS-sensitive iron regulation, (4) motility, (5) broad transcriptional reallocation toward growth, and (6) metabolic rewiring to decrease NADH production. This work thus demonstrates the power of iModulon knowledge mapping for evolution analysis.

Independent Component Analysis Reveals the Transcriptional Regulatory Modules in Bradyrhizobium diazoefficiens USDA110

The dynamic adaptation of bacteria to environmental changes is achieved through the coordinated expression of many genes, which constitutes a transcriptional regulatory network (TRN). Bradyrhizobium diazoefficiens USDA110 is an important model strain for the study of symbiotic nitrogen fixation (SNF), and its SNF ability largely depends on the TRN. In this study, independent component analysis was applied to 226 high-quality gene expression profiles of B. diazoefficiens USDA110 microarray datasets, from which 64 iModulons were identified. Using these iModulons and their condition-specific activity levels, we (1) provided new insights into the connection between the FixLJ-FixK2-FixK1 regulatory cascade and quorum sensing, (2) discovered the independence of the FixLJ-FixK2-FixK1 and NifA/RpoN regulatory cascades in response to oxygen, (3) identified the FixLJ-FixK2 cascade as a mediator connecting the FixK2-2 iModulon and the Phenylalanine iModulon, (4) described the differential activation of iModulons in B. diazoefficiens USDA110 under different environmental conditions, and (5) proposed a notion of active-TRN based on the changes in iModulon activity to better illustrate the relationship between gene regulation and environmental condition. In sum, this research offered an iModulon-based TRN for B. diazoefficiens USDA110, which formed a foundation for comprehensively understanding the intricate transcriptional regulation during SNF.

Systems biology of competency in Vibrio natriegens is revealed by applying novel data analytics to the transcriptome

Vibrio natriegens regulates natural competence through the TfoX and QstR transcription factors, which are involved in external DNA capture and transport. However, the extensive genetic and transcriptional regulatory basis for competency remains unknown. We used a machine-learning approach to decompose Vibrio natriegens's transcriptome into 45 groups of independently modulated sets of genes (iModulons). Our findings show that competency is associated with the repression of two housekeeping iModulons (iron metabolism and translation) and the activation of six iModulons; including TfoX and QstR, a novel iModulon of unknown function, and three housekeeping iModulons (representing motility, polycations, and reactive oxygen species [ROS] responses). Phenotypic screening of 83 gene deletion strains demonstrates that loss of iModulon function reduces or eliminates competency. This database-iModulon-discovery cycle unveils the transcriptomic basis for competency and its relationship to housekeeping functions. These results provide the genetic basis for systems biology of competency in this organism.

CyuR is a Dual Regulator for L-Cysteine Dependent Antimicrobial Resistance in Escherichia coli

Hydrogen sulfide (H ₂ S), mainly produced from L-cysteine (Cys), renders bacteria highly resistant to oxidative stress. This mitigation of oxidative stress was suggested to be an important survival mechanism to achieve antimicrobial resistance (AMR) in many pathogenic bacteria. CyuR (known as DecR or YbaO) is a recently characterized Cys-dependent transcription regulator, responsible for the activation of the cyuAP operon and generation of hydrogen sulfide from Cys. Despite its potential importance, the regulatory network of CyuR remains poorly understood. In this study, we investigated the roles of the CyuR regulon in a Cys-dependent AMR mechanism in E. coli strains. We found: 1) Cys metabolism has a significant role in AMR and its effect is conserved in many E. coli strains, including clinical isolates; 2) CyuR negatively controls the expression of mdlAB encoding a transporter that exports antibiotics such as cefazolin and vancomycin; 3) CyuR binds to a DNA sequence motif 'GAAwAAATTGTxGxxATTTsyCC' in the absence of Cys, confirmed by an in vitro binding assay; and 4) CyuR may regulate 25 additional genes as suggested by in silico motif scanning and transcriptome sequencing. Collectively, our findings expanded the understanding of the biological roles of CyuR relevant to antibiotic resistance associated with Cys.

RB-TnSeq identifies genetic targets for improved tolerance of Pseudomonas putida towards compounds relevant to lignin conversion

Lignin-derived mixtures intended for bioconversion commonly contain high concentrations of aromatic acids, aliphatic acids, and salts. The inherent toxicity of these chemicals places a significant bottleneck upon the effective use of microbial systems for the valorization of these mixtures. Pseudomonas putida KT2440 can tolerate stressful quantities of several lignin-related compounds, making this bacterium a promising host for converting these chemicals to valuable bioproducts. Nonetheless, further increasing P. putida tolerance to chemicals in lignin-rich substrates has the potential to improve bioprocess performance. Accordingly, we employed random barcoded transposon insertion sequencing (RB-TnSeq) to reveal genetic determinants in P. putida KT2440 that influence stress outcomes during exposure to representative constituents found in lignin-rich process streams. The fitness information obtained from the RB-TnSeq experiments informed engineering of strains via deletion or constitutive expression of several genes. Namely, ΔgacAS, ΔfleQ, ΔlapAB, ΔttgR::P_tac:ttgABC, P_tac:PP_1150:PP_1152, ΔrelA, and ΔPP_1430 mutants showed growth improvement in the presence of single compounds, and some also exhibited greater tolerance when grown using a complex chemical mixture representative of a lignin-rich chemical stream. Overall, this work demonstrates the successful implementation of a genome-scale screening tool for the identification of genes influencing stress tolerance against notable compounds within lignin-enriched chemical streams, and the genetic targets identified herein offer promising engineering targets for improving feedstock tolerance in lignin valorization strains of P. putida KT2440.

Proteome allocation is linked to transcriptional regulation through a modularized transcriptome

It has proved challenging to quantitatively relate the proteome to the transcriptome on a per-gene basis. Recent advances in data analytics have enabled a biologically meaningful modularization of the bacterial transcriptome. We thus investigated whether matched datasets of transcriptomes and proteomes from bacteria under diverse conditions could be modularized in the same way to reveal novel relationships between their compositions. We found that; 1) the modules of the proteome and the transcriptome are comprised of a similar list of gene products, 2) the modules in the proteome often represent combinations of modules from the transcriptome, 3) known transcriptional and post-translational regulation is reflected in differences between two sets of modules, allowing for knowledge-mapping when interpreting module functions, and 4) through statistical modeling, absolute proteome allocation can be inferred from the transcriptome alone. Quantitative and knowledge-based relationships can thus be found at the genome-scale between the proteome and transcriptome in bacteria.

The hitchhikers' guide to RNA sequencing and functional analysis

DNA and RNA sequencing technologies have revolutionized biology and biomedical sciences, sequencing full genomes and transcriptomes at very high speeds and reasonably low costs. RNA sequencing (RNA-Seq) enables transcript identification and quantification, but once sequencing has concluded researchers can be easily overwhelmed with questions such as how to go from raw data to differential expression (DE), pathway analysis and interpretation. Several pipelines and procedures have been developed to this effect. Even though there is no unique way to perform RNA-Seq analysis, it usually follows these steps: 1) raw reads quality check, 2) alignment of reads to a reference genome, 3) aligned reads' summarization according to an annotation file, 4) DE analysis and 5) gene set analysis and/or functional enrichment analysis. Each step requires researchers to make decisions, and the wide variety of options and resulting large volumes of data often lead to interpretation challenges. There also seems to be insufficient guidance on how best to obtain relevant information and derive actionable knowledge from transcription experiments. In this paper, we explain RNA-Seq steps in detail and outline differences and similarities of different popular options, as well as advantages and disadvantages. We also discuss non-coding RNA analysis, multi-omics, meta-transcriptomics and the use of artificial intelligence methods complementing the arsenal of tools available to researchers. Lastly, we perform a complete analysis from raw reads to DE and functional enrichment analysis, visually illustrating how results are not absolute truths and how algorithmic decisions can greatly impact results and interpretation.

Genetically encoded biosensors for microbial synthetic biology: From conceptual frameworks to practical applications

Genetically encoded biosensors are the vital components of synthetic biology and metabolic engineering, as they are regarded as powerful devices for the dynamic control of genotype metabolism and evolution/screening of desirable phenotypes. This review summarized the recent advances in the construction and applications of different genetically encoded biosensors, including fluorescent protein-based biosensors, nucleic acid-based biosensors, allosteric transcription factor-based biosensors and two-component system-based biosensors. First, the construction frameworks of these biosensors were outlined. Then, the recent progress of biosensor applications in creating versatile microbial cell factories for the bioproduction of high-value chemicals was summarized. Finally, the challenges and prospects for constructing robust and sophisticated biosensors were discussed. This review provided theoretical guidance for constructing genetically encoded biosensors to create desirable microbial cell factories for sustainable bioproduction.

Improving gene regulatory network inference and assessment: The importance of using network structure

Gene regulatory networks are graph models representing cellular transcription events. Networks are far from complete due to time and resource consumption for experimental validation and curation of the interactions. Previous assessments have shown the modest performance of the available network inference methods based on gene expression data. Here, we study several caveats on the inference of regulatory networks and methods assessment through the quality of the input data and gold standard, and the assessment approach with a focus on the global structure of the network. We used synthetic and biological data for the predictions and experimentally-validated biological networks as the gold standard (ground truth). Standard performance metrics and graph structural properties suggest that methods inferring co-expression networks should no longer be assessed equally with those inferring regulatory interactions. While methods inferring regulatory interactions perform better in global regulatory network inference than co-expression-based methods, the latter is better suited to infer function-specific regulons and co-regulation networks. When merging expression data, the size increase should outweigh the noise inclusion and graph structure should be considered when integrating the inferences. We conclude with guidelines to take advantage of inference methods and their assessment based on the applications and available expression datasets.

Genetic Code Expansion in Pseudomonas putida KT2440

Pseudomonas putida KT2440 is an emerging microbial chassis for biobased chemical production from renewable feedstocks and environmental bioremediation. However, tools for studying, engineering, and modulating protein complexes and biosynthetic enzymes in this organism are largely underdeveloped. Genetic code expansion for the incorporation of unnatural amino acids (unAAs) into proteins can advance such efforts and, furthermore, enable additional controls of biological processes of the strain. In this work, we established the orthogonality of two widely used archaeal tRNA synthetase and tRNA pairs in KT2440. Following the optimization of decoding systems, four unAAs were incorporated into proteins in response to a UAG stop codon at 34.6-78% efficiency. In addition, we demonstrated the utility of genetic code expansion through the incorporation of a photocross-linking amino acid, p-benzoyl-l-phenylalanine (pBpa), into glutathione S-transferase (GstA) and a chemosensory response regulator (CheY) for protein-protein interaction studies in KT2440. This work reported the successful genetic code expansion in KT2440 for the first time. Given the diverse structure and functions of unAAs that have been added to protein syntheses using the archaeal systems, our research lays down a solid foundation for future work to study and enhance the biological functions of KT2440.