Transcriptomic and metabolomic profiling of Zymomonas mobilis during aerobic and anaerobic fermentations

Biosensor-assisted CRISPRi high-throughput screening to identify genetic targets in Zymomonas mobilis for high d-lactate production

Lactate is an important monomer for the synthesis of poly-lactate (PLA), which is a substitute for the petrochemical plastics. To achieve the goal of high lactate titer, rate, and yield for commercial production, efficient lactate production pathway is needed as well as genetic targets that affect high lactate production and tolerance. In this study, an LldR-based d-lactate biosensor with a broad dynamic range was first applied into Zymomonas mobilis to select mutant strains with strong GFP fluorescence, which could be the mutant strains with increased d-lactate production. Then, LldR-based d-lactate biosensor was combined with a genome-wide CRISPR interference (CRISPRi) library targeting the entire genome to generate thousands of mutants with gRNA targeting different genetic targets across the whole genome. Specifically, two mutant libraries were selected containing 105 and 104 mutants with different interference sites from two rounds of fluorescence-activated cell sorting (FACS), respectively. Two genetic targets of ZMO1323 and ZMO1530 were characterized and confirmed to be associated with the increased d-lactate production, further knockout of ZMO1323 and ZMO1530 resulted in a 15% and 21% increase of d-lactate production, respectively. This work thus not only established a high-throughput approach that combines genome-scale CRISPRi and biosensor-assisted screening to identify genetic targets associated with d-lactate production in Z. mobilis, but also provided a feasible high-throughput screening approach for rapid identification of genetic targets associated with strain performance for other industrial microorganisms.

Lipid membrane remodeling and metabolic response during isobutanol and ethanol exposure in Zymomonas mobilis

BACKGROUND

Recent engineering efforts have targeted the ethanologenic bacterium Zymomonas mobilis for isobutanol production. However, significant hurdles remain due this organism's vulnerability to isobutanol toxicity, adversely affecting its growth and productivity. The limited understanding of the physiological impacts of isobutanol on Z. mobilis constrains our ability to overcome these production barriers.

RESULTS

We utilized a systems-level approach comprising LC-MS/MS-based lipidomics, metabolomics, and shotgun proteomics, to investigate how exposure to ethanol and isobutanol impact the lipid membrane composition and overall physiology of Z. mobilis. Our analysis revealed significant and distinct alterations in membrane phospholipid and fatty acid composition resulting from ethanol and isobutanol exposure. Notably, ethanol exposure increased membrane cyclopropane fatty acid content and expression of cyclopropane fatty acid (CFA) synthase. Surprisingly, isobutanol decreased cyclopropane fatty acid content despite robust upregulation of CFA synthase. Overexpression of the native Z. mobilis' CFA synthase increased cyclopropane fatty acid content in all phospholipid classes and was associated with a significant improvement in growth rates in the presence of added ethanol and isobutanol. Heterologous expression of CFA synthase from Clostridium acetobutylicum resulted in a near complete replacement of unsaturated fatty acids with cyclopropane fatty acids, affecting all lipid classes. However, this did not translate to improved growth rates under isobutanol exposure. Correlating with its greater susceptibility to isobutanol, Z. mobilis exhibited more pronounced alterations in its proteome, metabolome, and overall cell morphology-including cell swelling and formation of intracellular protein aggregates -when exposed to isobutanol compared to ethanol. Isobutanol triggered a broad stress response marked by the upregulation of heat shock proteins, efflux transporters, DNA repair systems, and the downregulation of cell motility proteins. Isobutanol also elicited widespread dysregulation of Z. mobilis' primary metabolism evidenced by increased levels of nucleotide degradation intermediates and the depletion of biosynthetic and glycolytic intermediates.

CONCLUSIONS

This study provides a comprehensive, systems-level evaluation of the impact of ethanol and isobutanol exposure on the lipid membrane composition and overall physiology of Z. mobilis. These findings will guide engineering of Z. mobilis towards the creation of isobutanol-tolerant strains that can serve as robust platforms for the industrial production of isobutanol from lignocellulosic sugars.

Better off breathing: an explanation for the seemingly detrimental impact of aerobic respiration on Zymomonas mobilis

The bacterium Zymomonas mobilis is best known for fermentatively producing more ethanol than yeast. However, Z. mobilis has also puzzled researchers for decades with the counterintuitive observation that disrupting aerobic respiration benefits aerobic growth, implying that fermentation remains favorable. Retention of detrimental respiration genes seemed to defy natural selection. New findings by Felczak et al. help clarify the importance of respiration for Z. mobilis and the factors that led to the confusing prior results (M. M. Felczak, M. P. Bernard, and M. A. TerAvest, 2023, mBio 14:e02043-23, https://doi.org/10.1128/mbio.02043-23). The team overcame redundancy from multiple genome copies to delete what turned out to be a key terminal oxidase. Unlike previous studies, wherein mutants exhibited low respiration rates and had improved aerobic growth, this mutant was incapable of respiration and had poor aerobic growth. Thus, respiration is important but surprisingly exceeds what is optimal under lab conditions. Respiration likely protects against toxic effects of oxygen, ensuring retention of respiration genes in the Z. mobilis genome.

The genetics of aerotolerant growth in an alphaproteobacterium with a naturally reduced genome

Reduced genome bacteria are genetically simplified systems that facilitate biological study and industrial use. The free-living alphaproteobacterium Zymomonas mobilis has a naturally reduced genome containing fewer than 2,000 protein-coding genes. Despite its small genome, Z. mobilis thrives in diverse conditions including the presence or absence of atmospheric oxygen. However, insufficient characterization of essential and conditionally essential genes has limited broader adoption of Z. mobilis as a model alphaproteobacterium. Here, we use genome-scale CRISPRi-seq (clustered regularly interspaced short palindromic repeats interference sequencing) to systematically identify and characterize Z. mobilis genes that are conditionally essential for aerotolerant or anaerobic growth or are generally essential across both conditions. Comparative genomics revealed that the essentiality of most "generally essential" genes was shared between Z. mobilis and other Alphaproteobacteria, validating Z. mobilis as a reduced genome model. Among conditionally essential genes, we found that the DNA repair gene, recJ, was critical only for aerobic growth but reduced the mutation rate under both conditions. Further, we show that genes encoding the F1FO ATP synthase and Rhodobacter nitrogen fixation (Rnf) respiratory complex are required for the anaerobic growth of Z. mobilis. Combining CRISPRi partial knockdowns with metabolomics and membrane potential measurements, we determined that the ATP synthase generates membrane potential that is consumed by Rnf to power downstream processes. Rnf knockdown strains accumulated isoprenoid biosynthesis intermediates, suggesting a key role for Rnf in powering essential biosynthetic reactions. Our work establishes Z. mobilis as a streamlined model for alphaproteobacterial genetics, has broad implications in bacterial energy coupling, and informs Z. mobilis genome manipulation for optimized production of valuable isoprenoid-based bioproducts. IMPORTANCE The inherent complexity of biological systems is a major barrier to our understanding of cellular physiology. Bacteria with markedly fewer genes than their close relatives, or reduced genome bacteria, are promising biological models with less complexity. Reduced genome bacteria can also have superior properties for industrial use, provided the reduction does not overly restrict strain robustness. Naturally reduced genome bacteria, such as the alphaproteobacterium Zymomonas mobilis, have fewer genes but remain environmentally robust. In this study, we show that Z. mobilis is a simplified genetic model for Alphaproteobacteria, a class with important impacts on the environment, human health, and industry. We also identify genes that are only required in the absence of atmospheric oxygen, uncovering players that maintain and utilize the cellular energy state. Our findings have broad implications for the genetics of Alphaproteobacteria and industrial use of Z. mobilis to create biofuels and bioproducts.

Current advances for omics-guided process optimization of microbial manufacturing

Currently, microbial manufacturing is widely used in various fields, such as food, medicine and energy, for its advantages of greenness and sustainable development. Process optimization is the committed step enabling the commercialization of microbial manufacturing products. However, the present optimization processes mainly rely on experience or trial-and-error method ignoring the intrinsic connection between cellular physiological requirement and production performance, so in many cases the productivity of microbial manufacturing could not been fully exploited at economically feasible cost. Recently, the rapid development of omics technologies facilitates the comprehensive analysis of microbial metabolism and fermentation performance from multi-levels of molecules, cells and microenvironment. The use of omics technologies makes the process optimization more explicit, boosting microbial manufacturing performance and bringing significant economic benefits and social value. In this paper, the traditional and omics technologies-guided process optimization of microbial manufacturing are systematically reviewed, and the future trend of process optimization is prospected.

Metabolic Engineering of Zymomonas mobilis for Acetoin Production by Carbon Redistribution and Cofactor Balance

Biorefinery to produce value-added biochemicals offers a promising alternative to meet our sustainable energy and environmental goals. Acetoin is widely used in the food and cosmetic industries as taste and fragrance enhancer. The generally regarded as safe (GRAS) bacterium Zymomonas mobilis produces acetoin as an extracellular product under aerobic conditions. In this study, metabolic engineering strategies were applied including redistributing the carbon flux to acetoin and manipulating the NADH levels. To improve the acetoin level, a heterologous acetoin pathway was first introduced into Z. mobilis, which contained genes encoding acetolactate synthase (Als) and acetolactate decarboxylase (AldC) driven by a strong native promoter Pgap. Then a gene encoding water-forming NADH oxidase (NoxE) was introduced for NADH cofactor balance. The recombinant Z. mobilis strain containing both an artificial acetoin operon and the noxE greatly enhanced acetoin production with maximum titer reaching 8.8 g/L and the productivity of 0.34 g∙L−1∙h−1. In addition, the strategies to delete ndh gene for redox balance by native I-F CRISPR-Cas system and to redirect carbon from ethanol production to acetoin biosynthesis through a dcas12a-based CRISPRi system targeting pdc gene laid a foundation to help construct an acetoin producer in the future. This study thus provides an informative strategy and method to harness the NADH levels for biorefinery and synthetic biology studies in Z. mobilis.

A plasmid-free Zymomonas mobilis mutant strain reducing reactive oxygen species for efficient bioethanol production using industrial effluent of xylose mother liquor

Genome minimization is an effective way for industrial chassis development. In this study, Zymomonas mobilis ZMNP, a plasmid-free mutant strain of Z. mobilis ZM4 with four native plasmids deleted, was constructed using native type I-F CRISPR-Cas system. Cell growth of ZMNP under different temperatures and industrial effluent of xylose mother liquor were examined to investigate the impact of native plasmid removal. Despite ZMNP grew similarly as ZM4 under different temperatures, ZMNP had better xylose mother liquor utilization than ZM4. In addition, genomic, transcriptomic, and proteomic analyses were applied to unravel the molecular changes between ZM4 and ZMNP. Whole-genome resequencing result indicated that an S267P mutation in the C-terminal of OxyR, a peroxide-sensing transcriptional regulator, probably alters the transcription initiation of antioxidant genes for stress responses. Transcriptomic and proteomic studies illustrated that the reason that ZMNP utilized the toxic xylose mother liquor better than ZM4 was probably due to the upregulation of genes in ZMNP involving in stress responses as well as cysteine biosynthesis to accelerate the intracellular ROS detoxification and nucleic acid damage repair. This was further confirmed by lower ROS levels in ZMNP compared to ZM4 in different media supplemented with furfural or ethanol. The upregulation of stress response genes due to the OxyR mutation to accelerate ROS detoxification and DNA/RNA repair not only illustrates the underlying mechanism of the robustness of ZMNP in the toxic xylose mother liquor, but also provides an idea for the rational design of synthetic inhibitor-tolerant microorganisms for economic lignocellulosic biochemical production.

The potential use of Zymomonas mobilis for the food industry

Zymomonas mobilis is a gram-negative facultative anaerobic spore, which is generally recognized as a safe. As a promising ethanologenic organism for large-scale bio-ethanol production, Z. mobilis has also shown a good application prospect in food processing and food additive synthesis for its unique physiological characteristics and excellent industrial characteristics. It not only has obvious advantages in food processing and becomes the biorefinery chassis cell for food additives, but also has a certain healthcare effect on human health. Until to now, most of the research is still in theory and laboratory scale, and further research is also needed to achieve industrial production. This review summarized the physiological characteristics and advantages of Z. mobilis in food industry for the first time and further expounds its research status in food industry from three aspects of food additive synthesis, fermentation applications, and prebiotic efficacy, it will provide a theoretical basis for its development and applications in food industry. This review also discussed the shortcomings of its practical applications in the current food industry, and explored other ways to broaden the applications of Z. mobilis in the food industry, to promote its applications in food processing.

Expression of Escherichia coli malic enzyme gene in Zymomonas mobilis for production of malic acid

Malic acid is one of the organic acids which is used in various industries including food and pharmaceuticals. Biotechnological production of malic acid by an efficient microorganism is highly desirable as the process will be eco-friendly and cost-effective. In this study, malic acid synthesis by Zymomonas mobilis was studied by expressing Escherichia coli malic enzyme gene under Pchap, Ptac and Ppdc promoters. The mae⁺ recombinants were obtained by recombineering-based genomic integration of Pchap-mae, Ptac-mae and Ppdc-mae sequences. The Ppdc promoter showed the highest expression of malic enzyme and the Pchap the lowest. However, cell growth was limited in mae⁺ recombinant containing Ppdc promoter. The metabolic analysis showed the highest level of malic acid in Ppdc-mae recombinant (2.84 g/L), which was about eight times higher than that in the wild type strain. The study showed that these three promoters can be used to produce organic acids in Z. mobilis.

Deciphering Molecular Mechanism Underlying Self-Flocculation of Zymomonas mobilis for Robust Production

Zymomonas mobilis metabolizes sugar anaerobically through the Entner-Doudoroff pathway with less ATP generated for lower biomass accumulation to direct more sugar for product formation with improved yield, making it a suitable host to be engineered as microbial cell factories for producing bulk commodities with major costs from feedstock consumption. Self-flocculation of the bacterial cells presents many advantages, such as enhanced tolerance to environmental stresses, a prerequisite for achieving high product titers by using concentrated substrates. ZM401, a self-flocculating mutant developed from ZM4, the unicellular model strain of Z. mobilis, was employed in this work to explore the molecular mechanism underlying this self-flocculating phenotype. Comparative studies between ZM401 and ZM4 indicate that a frameshift caused by a single nucleotide deletion in the poly-T tract of ZMO1082 fused the putative gene with the open reading frame of ZMO1083, encoding the catalytic subunit BcsA of the bacterial cellulose synthase to catalyze cellulose biosynthesis. Furthermore, the single nucleotide polymorphism mutation in the open reading frame of ZMO1055, encoding a bifunctional GGDEF-EAL protein with apparent diguanylate cyclase/phosphodiesterase activities, resulted in the Ala526Val substitution, which consequently compromised in vivo specific phosphodiesterase activity for the degradation of cyclic diguanylic acid, leading to intracellular accumulation of the signaling molecule to activate cellulose biosynthesis. These discoveries are significant for engineering other unicellular strains from Z. mobilis with the self-flocculating phenotype for robust production. IMPORTANCE Stress tolerance is a prerequisite for microbial cell factories to be robust in production, particularly for biorefinery of lignocellulosic biomass to produce biofuels, bioenergy, and bio-based chemicals for sustainable socioeconomic development, since various inhibitors are released during the pretreatment to destroy the recalcitrant lignin-carbohydrate complex for sugar production through enzymatic hydrolysis of the cellulose component, and their detoxification is too costly for producing bulk commodities. Although tolerance to individual stress has been intensively studied, the progress seems less significant since microbial cells are inevitably suffering from multiple stresses simultaneously under production conditions. When self-flocculating, microbial cells are more tolerant to multiple stresses through the general stress response due to enhanced quorum sensing associated with the morphological change for physiological and metabolic advantages. Therefore, elucidation of the molecular mechanism underlying such a self-flocculating phenotype is significant for engineering microbial cells with the unique multicellular morphology through rational design to boost their production performance.

Mechanism for Reducing the Horizontal Transfer Risk of the Airborne Antibiotic-Resistant Genes of Escherichia coli Species through Microwave or UV Irradiation

Antibiotic-resistant bacteria (ARBs) and antibiotic-resistant genes (ARGs) as new types of contaminants are discharged into the environment, increasing the risk of horizontal gene transfer (HGT). However, few researchers have examined the impacts of airborne ARB deactivation on HGT risk. The deactivation of airborne Escherichia coli 10667 (carrying sul genes) and the emission and removal of ARGs were mainly investigated in this study. Moreover, the potential mechanisms of HGT and transfer frequencies under microwave (MW) and ultraviolet (UV) irradiation were investigated using the nonresistant E. coli GMCC 13373 and E. coli DH5α with plasmid RP4 as the recipient and donor, respectively. E. coli CICC 10667 and E. coli DH5α with RP4 plasmid achieve log inactivation values as high as 5.5-log and 5.0-log, respectively, which were quite different from the antibiotic-sensitive strain E. coli CGMCC 13373 (3.4-log) subjected to MW irradiation. For UV disinfection, E. coli DH5α with the RP4 plasmid was reduced at 4.4-log, E. coli CGMCC 13373 was reduced at 2.3-log, and E. coli CICC 10667 was inactivated at 2.1-log. The removal rates of ARGs and HGT frequencies under MW irradiation were compared with those under UV irradiation. The ARGs removal efficiency (85.5%) obtained by MW was higher than that obtained by UV (48.2%). Consequently, the HGT frequency (0.008) of airborne ARGs released to the recipient (forward transfer) decreased and was lower than that under UV irradiation (0.014). Moreover, the plasmid RP4 was transferred from the donor to the surviving damaged E. coli 10667 as cell permeability (reverse transfer) was increased at a high HGT frequency (0.003) by MW, which was close to the value by UV (0.002). Additionally, sul1 and sul2 genes were confirmed to be more resistant to MW than the sul3 gene. These findings reveal the mechanism of HGT between damaged E. coli 10667 and surrounding environmental microbes. Microwave is a promising technology for disinfecting airborne microbes and preventing the spread of antibiotic resistance.

Zymomonas mobilis ZM4 Utilizes an NADP+-Dependent Acetaldehyde Dehydrogenase To Produce Acetate

Zymomonas mobilis is a promising bacterial host for biofuel production, but further improvement has been hindered because some aspects of its metabolism remain poorly understood. For example, one of the main by-products generated by Z. mobilis is acetate, but the pathway for acetate production is unknown. Acetaldehyde oxidation has been proposed as the major source of acetate, and an acetaldehyde dehydrogenase was previously isolated from Z. mobilis via activity guided fractionation, but the corresponding gene has never been identified. We determined that the locus ZMO1754 (also known as ZMO_RS07890) encodes an NADP⁺-dependent acetaldehyde dehydrogenase that is responsible for acetate production by Z. mobilis. Deletion of this gene from the chromosome resulted in a growth defect in oxic conditions, suggesting that acetaldehyde detoxification is an important role of acetaldehyde dehydrogenase. The deletion strain also exhibited a near complete abolition of acetate production, both in typical laboratory conditions and during lignocellulosic hydrolysate fermentation. Our results show that ZMO1754 encodes the major acetate-forming acetaldehyde dehydrogenase in Z. mobilis, and we therefore rename the gene aldB based on functional similarity to the Escherichia coli acetaldehyde dehydrogenase. IMPORTANCE Biofuel production from nonfood crops is an important strategy for reducing carbon emissions from the transportation industry, but it has not yet become commercially viable. An important avenue to improve biofuel production is to enhance the characteristics of fermentation organisms by decreasing by-product formation via genetic engineering. Here, we identified and deleted a metabolic pathway and associated gene that lead to acetate formation in Zymomonas mobilis. Acetate is one of the major by-products generated during ethanol production by Z. mobilis, so this information may be used in the future to develop better strains for commercial biofuel production.

Antimicrobial Activity of Zymomonas mobilis Is Related to Its Aerobic Catabolism and Acid Resistance

Zymomonas mobilis is an ethanologenic, facultatively anaerobic alpha-proteobacterium, known for its inhibitory effect on the growth of a wide variety of microorganisms. This property might be interesting for the design of novel antimicrobials, yet it has negative implications for biotechnology, as it hinders the use of Z. mobilis as a producer microorganism in cocultivation. So far, the chemical nature of its inhibitory compound(s) has not been established. In the present study, we demonstrate that the putative inhibitor is a low-molecular-weight (below 3 kDa), thermostable compound, resistant to protease treatment, which is synthesized under aerobic conditions in Z. mobilis strains via the active respiratory chain. It is also synthesized by aerated nongrowing, glucose-consuming cells in the presence of chloramphenicol, thus ruling out its bacteriocin-like peptide nature. The inhibitory activity is pH-dependent and strongly correlated with the accumulation of propionate and acetate in the culture medium. Although, in Z. mobilis, the synthesis pathways of these acids still need to be identified, the acid production depends on respiration, and is much less pronounced in the non-respiring mutant strain, which shows low inhibitory activity. We conclude that propionate and acetate play a central role in the antimicrobial effects of Z. mobilis, which itself is known to bear high resistance to organic acids.

Multi-Omics Integration and Network Analysis Reveal Potential Hub Genes and Genetic Mechanisms Regulating Bovine Mastitis

Mastitis, inflammation of the mammary gland, is the most prevalent disease in dairy cattle that has a potential impact on profitability and animal welfare. Specifically designed multi-omics studies can be used to prioritize candidate genes and identify biomarkers and the molecular mechanisms underlying mastitis in dairy cattle. Hence, the present study aimed to explore the genetic basis of bovine mastitis by integrating microarray and RNA-Seq data containing healthy and mastitic samples in comparative transcriptome analysis with the results of published genome-wide association studies (GWAS) using a literature mining approach. The integration of different information sources resulted in the identification of 33 common and relevant genes associated with bovine mastitis. Among these, seven genes—CXCR1, HCK, IL1RN, MMP9, S100A9, GRO1, and SOCS3—were identified as the hub genes (highly connected genes) for mastitis susceptibility and resistance, and were subjected to protein-protein interaction (PPI) network and gene regulatory network construction. Gene ontology annotation and enrichment analysis revealed 23, 7, and 4 GO terms related to mastitis in the biological process, molecular function, and cellular component categories, respectively. Moreover, the main metabolic-signalling pathways responsible for the regulation of immune or inflammatory responses were significantly enriched in cytokine–cytokine-receptor interaction, the IL-17 signaling pathway, viral protein interaction with cytokines and cytokine receptors, and the chemokine signaling pathway. Consequently, the identification of these genes, pathways, and their respective functions could contribute to a better understanding of the genetics and mechanisms regulating mastitis and can be considered a starting point for future studies on bovine mastitis.

Genome-Wide Analyses of Proteome and Acetylome in Zymomonas mobilis Under N2-Fixing Condition

Zymomonas mobilis, a promising candidate for industrial biofuel production, is capable of nitrogen fixation naturally without hindering ethanol production. However, little is known about the regulation of nitrogen fixation in Z. mobilis. We herein conducted a high throughput analysis of proteome and protein acetylation in Z. mobilis under N₂-fixing conditions and established its first acetylome. The upregulated proteins mainly belong to processes of nitrogen fixation, motility, chemotaxis, flagellar assembly, energy production, transportation, and oxidation–reduction. Whereas, downregulated proteins are mainly related to energy-consuming and biosynthetic processes. Our acetylome analyses revealed 197 uniquely acetylated proteins, belonging to major pathways such as nitrogen fixation, central carbon metabolism, ammonia assimilation pathway, protein biosynthesis, and amino acid metabolism. Further, we observed acetylation in glycolytic enzymes of central carbon metabolism, the nitrogenase complex, the master regulator NifA, and the enzyme in GS/GOGAT cycle. These findings suggest that protein acetylation may play an important role in regulating various aspects of N₂-metabolism in Z. mobilis. This study provides new knowledge of specific proteins and their associated cellular processes and pathways that may be regulated by protein acetylation in Z. mobilis.

Engineering Zymomonas mobilis for the Production of Xylonic Acid from Sugarcane Bagasse Hydrolysate

Sugarcane bagasse is an agricultural residue rich in xylose, which may be used as a feedstock for the production of high-value-added chemicals, such as xylonic acid, an organic acid listed as one of the top 30 value-added chemicals on a NREL report. Here, Zymomonas mobilis was engineered for the first time to produce xylonic acid from sugarcane bagasse hydrolysate. Seven coding genes for xylose dehydrogenase (XDH) were tested. The expression of XDH gene from Paraburkholderia xenovorans allowed the highest production of xylonic acid (26.17 ± 0.58 g L-1) from 50 g L-1 xylose in shake flasks, with a productivity of 1.85 ± 0.06 g L-1 h-1 and a yield of 1.04 ± 0.04 gAX/gX. Deletion of the xylose reductase gene further increased the production of xylonic acid to 56.44 ± 1.93 g L-1 from 54.27 ± 0.26 g L-1 xylose in a bioreactor. Strain performance was also evaluated in sugarcane bagasse hydrolysate as a cheap feedstock, which resulted in the production of 11.13 g L-1 xylonic acid from 10 g L-1 xylose. The results show that Z. mobilis may be regarded as a potential platform for the production of organic acids from cheap lignocellulosic biomass in the context of biorefineries.

Model-driven analysis of mutant fitness experiments improves genome-scale metabolic models of Zymomonas mobilis ZM4

Genome-scale metabolic models have been utilized extensively in the study and engineering of the organisms they describe. Here we present the analysis of a published dataset from pooled transposon mutant fitness experiments as an approach for improving the accuracy and gene-reaction associations of a metabolic model for Zymomonas mobilis ZM4, an industrially relevant ethanologenic organism with extremely high glycolytic flux and low biomass yield. Gene essentiality predictions made by the draft model were compared to data from individual pooled mutant experiments to identify areas of the model requiring deeper validation. Subsequent experiments showed that some of the discrepancies between the model and dataset were caused by polar effects, mis-mapped barcodes, or mutants carrying both wild-type and transposon disrupted gene copies-highlighting potential limitations inherent to data from individual mutants in these high-throughput datasets. Therefore, we analyzed correlations in fitness scores across all 492 experiments in the dataset in the context of functionally related metabolic reaction modules identified within the model via flux coupling analysis. These correlations were used to identify candidate genes for a reaction in histidine biosynthesis lacking an annotated gene and highlight metabolic modules with poorly correlated gene fitness scores. Additional genes for reactions involved in biotin, ubiquinone, and pyridoxine biosynthesis in Z. mobilis were identified and confirmed using mutant complementation experiments. These discovered genes, were incorporated into the final model, iZM4_478, which contains 747 metabolic and transport reactions (of which 612 have gene-protein-reaction associations), 478 genes, and 616 unique metabolites, making it one of the most complete models of Z. mobilis ZM4 to date. The methods of analysis that we applied here with the Z. mobilis transposon mutant dataset, could easily be utilized to improve future genome-scale metabolic reconstructions for organisms where these, or similar, high-throughput datasets are available.

Genome-Scale Transcription-Translation Mapping Reveals Features of Zymomonas mobilis Transcription Units and Promoters

Efforts to rationally engineer synthetic pathways in Zymomonas mobilis are impeded by a lack of knowledge and tools for predictable and quantitative programming of gene regulation at the transcriptional, posttranscriptional, and posttranslational levels. With the detailed functional characterization of the Z. mobilis genome presented in this work, we provide crucial knowledge for the development of synthetic genetic parts tailored to Z. mobilis. This information is vital as researchers continue to develop Z. mobilis for synthetic biology applications. Our methods and statistical analyses also provide ways to rapidly advance the understanding of poorly characterized bacteria via empirical data that enable the experimental validation of sequence-based prediction for genome characterization and annotation.

Zymomonas mobilis is an ethanologenic alphaproteobacterium with promise for the industrial conversion of renewable plant biomass into fuels and chemical bioproducts. Limited functional annotation of the Z. mobilis genome is a current barrier to both fundamental studies of Z. mobilis and its development as a synthetic biology chassis. To gain insight, we collected sample-matched multiomics data, including RNA sequencing (RNA-seq), transcription start site (TSS) sequencing (TSS-seq), termination sequencing (term-seq), ribosome profiling, and label-free shotgun proteomic mass spectrometry, across different growth conditions and used these data to improve annotation and assign functional sites in the Z. mobilis genome. Proteomics and ribosome profiling informed revisions of protein-coding genes, which included 44 start codon changes and 42 added proteins. We developed statistical methods for annotating transcript 5′ and 3′ ends, enabling the identification of 3,940 TSSs and their corresponding promoters and 2,091 transcription termination sites, which were distinguished from RNA processing sites by the lack of an adjacent RNA 5′ end. Our results revealed that Z. mobilis σ^A −35 and −10 promoter elements closely resemble canonical Escherichia coli −35 and −10 elements, with one notable exception: the Z. mobilis −10 element lacks the highly conserved −7 thymine observed in E. coli and other previously characterized σ^A promoters. The σ^A promoters of another alphaproteobacterium, Caulobacter crescentus, similarly lack the conservation of −7 thymine in their −10 elements. Our results anchor the development of Z. mobilis as a platform for synthetic biology and establish strategies for empirical genome annotation that can complement purely computational methods.

IMPORTANCE Efforts to rationally engineer synthetic pathways in Zymomonas mobilis are impeded by a lack of knowledge and tools for predictable and quantitative programming of gene regulation at the transcriptional, posttranscriptional, and posttranslational levels. With the detailed functional characterization of the Z. mobilis genome presented in this work, we provide crucial knowledge for the development of synthetic genetic parts tailored to Z. mobilis. This information is vital as researchers continue to develop Z. mobilis for synthetic biology applications. Our methods and statistical analyses also provide ways to rapidly advance the understanding of poorly characterized bacteria via empirical data that enable the experimental validation of sequence-based prediction for genome characterization and annotation.

Characterization and repurposing of the endogenous Type I-F CRISPR-Cas system of Zymomonas mobilis for genome engineering

Application of CRISPR-based technologies in non-model microorganisms is currently very limited. Here, we reported efficient genome engineering of an important industrial microorganism, Zymomonas mobilis, by repurposing the endogenous Type I-F CRISPR–Cas system upon its functional characterization. This toolkit included a series of genome engineering plasmids, each carrying an artificial self-targeting CRISPR and a donor DNA for the recovery of recombinants. Through this toolkit, various genome engineering purposes were efficiently achieved, including knockout of ZMO0038 (100% efficiency), cas2/3 (100%), and a genomic fragment of >10 kb (50%), replacement of cas2/3 with mCherry gene (100%), in situ nucleotide substitution (100%) and His-tagging of ZMO0038 (100%), and multiplex gene deletion (18.75%) upon optimal donor size determination. Additionally, the Type I-F system was further applied for CRISPRi upon Cas2/3 depletion, which has been demonstrated to successfully silence the chromosomally integrated mCherry gene with its fluorescence intensity reduced by up to 88%. Moreover, we demonstrated that genome engineering efficiency could be improved under a restriction–modification (R–M) deficient background, suggesting the perturbance of genome editing by other co-existing DNA targeting modules such as the R–M system. This study might shed light on exploiting and improving CRISPR–Cas systems in other microorganisms for genome editing and metabolic engineering practices.

A reconciliation of genome-scale metabolic network model of Zymomonas mobilis ZM4

Zymomonas mobilis ZM4 has recently been used for a variety of biotechnological purposes. To rationally enhance its metabolic performance, a reliable genome-scale metabolic network model (GEM) of this organism is required. To this end, we reconstructed a genome-scale metabolic model (iHN446) for Z. mobilis, which involves 446 genes, 859 reactions, and 894 metabolites. We started by first reconciling the existing GEMs previously constructed for Z. mobilis to obtain a draft network. Next, recent gene annotations, up-to-date literature, physiological data and biochemical databases were used to upgrade the network. Afterward, the draft network went through a curative and iterative process of gap-filling by computational tools and manual refinement. The final model was evaluated using experimental data and literature information. We next applied this model as a platform for analyzing the links between transcriptome-flux and transcriptome-metabolome. We found that experimental observations were in agreement with the predicted results from our final GEM. Taken together, this comprehensive model (iHN446) can be utilized for studying metabolism in Z. mobilis and finding rational targets for metabolic engineering applications.

Multiple Small RNAs Interact to Co-regulate Ethanol Tolerance in Zymomonas mobilis

sRNAs represent a powerful class of regulators that influences multiple mRNA targets in response to environmental changes. However, very few direct sRNA-sRNA interactions have been deeply studied in any organism. Zymomonas mobilis is a bacterium with unique ethanol-producing metabolic pathways in which multiple small RNAs (sRNAs) have recently been identified, some of which show differential expression in ethanol stress. In this study, we show that two sRNAs (Zms4 and Zms6) are upregulated under ethanol stress and have significant impacts on ethanol tolerance and production in Z. mobilis. We conducted multi-omics analysis (combining transcriptomics and sRNA-immunoprecipitation) to map gene networks under the influence of their regulation. We confirmed that Zms4 and Zms6 bind multiple RNA targets and regulate their expressions, influencing many downstream pathways important to ethanol tolerance and production. In particular, Zms4 and Zms6 interact with each other as well as many other sRNAs, forming a novel sRNA-sRNA direct interaction network. This study thus uncovers a sRNA network that co-orchestrates multiple ethanol related pathways through a diverse set of mRNA targets and a large number of sRNAs. To our knowledge, this study represents one of the largest sRNA-sRNA direct interactions uncovered so far.

Metabolic Profiling of Glucose-Fed Metabolically Active Resting Zymomonas mobilis Strains

Zymomonas mobilis is the most efficient bacterial ethanol producer and its physiology is potentially applicable to industrial-scale bioethanol production. However, compared to other industrially important microorganisms, the Z. mobilis metabolome and adaptation to various nutritional and genetic perturbations have been poorly characterized. For rational metabolic engineering, it is essential to understand how central metabolism and intracellular redox balance are maintained in Z. mobilis under various conditions. In this study, we applied quantitative mass spectrometry-based metabolomics to explore how glucose-fed non-growing Z. mobilis Zm6 cells metabolically adapt to change of oxygen availability. Mutants partially impaired in ethanol synthesis (Zm6 adhB) or oxidative stress response (Zm6 cat) were also examined. Distinct patterns of adaptation of central metabolite pools due to the change in cultivation condition and between the mutants and Zm6 reference strain were observed. Decreased NADH/NAD ratio under aerobic incubation corresponded to higher concentrations of the phosphorylated glycolytic intermediates, in accordance with predictions of the kinetic model of Entner–Doudoroff pathway. The effects on the metabolite pools of aerobic to anaerobic transition were similar in the mutants, yet less pronounced. The present data on metabolic plasticity of non-growing Z. mobilis cells will facilitate the further metabolic engineering of the respective strains and their application as biocatalysts.

Transcriptomic Profiles of Zymomonas mobilis 8b to Furfural Acute and Long-Term Stress in Both Glucose and Xylose Conditions

Zymomonas mobilis 8b is an ethanologenic bacterium engineered to utilize both glucose and xylose. The impacts of lignocellulosic hydrolyzate inhibitors on the growth of Zymomonas mobilis 8b have been investigated. However, the molecular responses of these inhibitors have not been completely elucidated yet. In this study, molecular responses to furfural were investigated using transcriptomic approaches of both chip-based microarray and a directional mRNA-Seq. Furfural acute shock time-course experiment with 3 g/L furfural supplemented when cells reached exponential phase and stress response experiment in the presence of 2 g/L furfural from the beginning of fermentation were carried out to study the physiological and transcriptional profiles of short-term and long-term effects of furfural on 8b. Furfural negatively affected 8b growth in terms of final biomass and the fermentation time. Transcriptomic studies indicated that the response of 8b to furfural was dynamic and complex, and differences existed between short-term shock and long-term stress responses. However, the gene function categories were similar with most down-regulated genes related to translation and biosynthesis, while the furfural up-regulated genes were mostly related to general stress responses. Several gene candidates have been identified and genetic studies indicated that expression of ZMO0465 and cysteine synthase operon ZMO0003-0006 driven by its native promoter in a shuttle vector enhanced the furfural tolerance of 8b. In addition, the relationship between microarray and mRNA-Seq was compared with good correlations. The directional mRNA-Seq data not only provided the gene expression profiling, but also can be applied for transcriptional architecture improvement to identify and confirm operons, novel transcripts, hypothetical gene functions, transcriptional start sites, and promoters with different strength.

Applying a New REFINE Approach in Zymomonas mobilis Identifies Novel sRNAs That Confer Improved Stress Tolerance Phenotypes

As global controllers of gene expression, small RNAs represent powerful tools for engineering complex phenotypes. However, a general challenge prevents the more widespread use of sRNA engineering strategies: mechanistic analysis of these regulators in bacteria lags far behind their high-throughput search and discovery. This makes it difficult to understand how to efficiently identify useful sRNAs to engineer a phenotype of interest. To help address this, we developed a forward systems approach to identify naturally occurring sRNAs relevant to a desired phenotype: RNA-seq Examiner for Phenotype-Informed Network Engineering (REFINE). This pipeline uses existing RNA-seq datasets under different growth conditions. It filters the total transcriptome to locate and rank regulatory-RNA-containing regions that can influence a metabolic phenotype of interest, without the need for previous mechanistic characterization. Application of this approach led to the uncovering of six novel sRNAs related to ethanol tolerance in non-model ethanol-producing bacterium Zymomonas mobilis. Furthermore, upon overexpressing multiple sRNA candidates predicted by REFINE, we demonstrate improved ethanol tolerance reflected by up to an approximately twofold increase in relative growth rate compared to controls not expressing these sRNAs in 7% ethanol (v/v) RMG-supplemented media. In this way, the REFINE approach informs strain-engineering strategies that we expect are applicable for general strain engineering.

Metabolic engineering of Zymomonas mobilis for anaerobic isobutanol production

BACKGROUND

Biofuels and value-added biochemicals derived from renewable biomass via biochemical conversion have attracted considerable attention to meet global sustainable energy and environmental goals. Isobutanol is a four-carbon alcohol with many advantages that make it attractive as a fossil-fuel alternative. Zymomonas mobilis is a highly efficient, anaerobic, ethanologenic bacterium making it a promising industrial platform for use in a biorefinery.

RESULTS

In this study, the effect of isobutanol on Z. mobilis was investigated, and various isobutanol-producing recombinant strains were constructed. The results showed that the Z. mobilis parental strain was able to grow in the presence of isobutanol below 12 g/L while concentrations greater than 16 g/L inhibited cell growth. Integration of the heterologous gene encoding 2-ketoisovalerate decarboxylase such as kdcA from Lactococcus lactis is required for isobutanol production in Z. mobilis. Moreover, isobutanol production increased from nearly zero to 100-150 mg/L in recombinant strains containing the kdcA gene driven by the tetracycline-inducible promoter Ptet. In addition, we determined that overexpression of a heterologous als gene and two native genes (ilvC and ilvD) involved in valine metabolism in a recombinant Z. mobilis strain expressing kdcA can divert pyruvate from ethanol production to isobutanol biosynthesis. This engineering improved isobutanol production to above 1 g/L. Finally, recombinant strains containing both a synthetic operon, als-ilvC-ilvD, driven by Ptet and the kdcA gene driven by the constitutive strong promoter, Pgap, were determined to greatly enhance isobutanol production with a maximum titer about 4.0 g/L. Finally, isobutanol production was negatively affected by aeration with more isobutanol being produced in more poorly aerated flasks.

CONCLUSIONS

This study demonstrated that overexpression of kdcA in combination with a synthetic heterologous operon, als-ilvC-ilvD, is crucial for diverting pyruvate from ethanol production for enhanced isobutanol biosynthesis. Moreover, this study also provides a strategy for harnessing the valine metabolic pathway for future production of other pyruvate-derived biochemicals in Z. mobilis.

Zymomonas mobilis metabolism: Novel tools and targets for its rational engineering

Zymomonas mobilis is an α-proteobacterium that interests the biofuel industry due to its perfect ethanol fermentation yields. From its first description as a bacterial isolate in fermented alcoholic beverages to date, Z. mobilis has been rigorously studied in directions basic and applied. The Z. mobilis powerful Entner-Doudoroff glycolytic pathway has been the center of rigorous biochemical studies and, aside from ethanol, it has attracted interest in terms of high-added-value chemical manufacturing. Energetic balances and the effects of respiration have been explored in fundamental directions as also in applications pursuing strain enhancement and the utilization of alternative carbon sources. Metabolic modeling has addressed the optimization of the biochemical circuitry at various conditions of growth and/or substrate utilization; it has been also critical in predicting desirable end-product yields via flux redirection. Lastly, stress tolerance has received particular attention, since it directly determines biocatalytical performance at challenging bioreactor conditions. At a genetic level, advances in the genetic engineering of the organism have brought forth beneficial manipulations in the Z. mobilis gene pool, e.g., knock-outs, knock-ins and gene stacking, aiming to broaden the metabolic repertoire and increase robustness. Recent omic and expressional studies shed light on the genomic content of the most applied strains and reveal landscapes of activity manifested at ambient or reactor-based conditions. Studies such as those reviewed in this work, contribute to the understanding of the biology of Z. mobilis, enable insightful strain development, and pave the way for the transformation of Z. mobilis into a consummate organism for biomass conversion.

Identification and characterization of ethanol-inducible promoters of Zymomonas mobilis based on omics data and dual reporter-gene system

Zymomonas mobilis is a model bacterial ethanologen and has been engineered to produce lignocellulosic biofuels and biochemicals such as 2,3-butanediol. We have previously identified promoters of different strengths using systems biology datasets and characterized them using the flow cytometry-based dual reporter-gene system. Here, we further demonstrated the capability of applying the dual reporter-gene system and omics datasets on discovering inducible promoters. Ten candidate ethanol-inducible promoters were identified through omics datasets mining and clustering. Using the dual reporter-gene system, these promoters were characterized under natural growth, ethanol stress, and ethanol-induced condition to investigate the transcriptional strength and ethanol inducibility. The results demonstrated that three promoters of P0405, P0435, and P0038 driving the expression of native genes of ZMO0405, ZMO0435, and ZMO0038, correspondingly, are potential ethanol-responsive promoters and may be growth related. This study not only identified and verified three ethanol-inducible promoters as biological parts, which can be used to synchronize the expression of heterologous pathway genes with the ethanol production process of Z. mobilis, but also demonstrated the power of combining omics datasets and dual reporter-gene system to identify biological parts for metabolic engineering and synthetic biology applications in Z. mobilis and related microorganisms.

Establishment and application of a CRISPR-Cas12a assisted genome-editing system in Zymomonas mobilis

BACKGROUND

Efficient and convenient genome-editing toolkits can expedite genomic research and strain improvement for desirable phenotypes. Zymomonas mobilis is a highly efficient ethanol-producing bacterium with a small genome size and desirable industrial characteristics, which makes it a promising chassis for biorefinery and synthetic biology studies. While classical techniques for genetic manipulation are available for Z. mobilis, efficient genetic engineering toolkits enabling rapidly systematic and high-throughput genome editing in Z. mobilis are still lacking.

RESULTS

Using Cas12a (Cpf1) from Francisella novicida, a recombinant strain with inducible cas12a expression for genome editing was constructed in Z. mobilis ZM4, which can be used to mediate RNA-guided DNA cleavage at targeted genomic loci. gRNAs were then designed targeting the replicons of native plasmids of ZM4 with about 100% curing efficiency for three native plasmids. In addition, CRISPR-Cas12a recombineering was used to promote gene deletion and insertion in one step efficiently and precisely with efficiency up to 90%. Combined with single-stranded DNA (ssDNA), CRISPR-Cas12a system was also applied to introduce minor nucleotide modification precisely into the genome with high fidelity. Furthermore, the CRISPR-Cas12a system was employed to introduce a heterologous lactate dehydrogenase into Z. mobilis with a recombinant lactate-producing strain constructed.

CONCLUSIONS

This study applied CRISPR-Cas12a in Z. mobilis and established a genome editing tool for efficient and convenient genome engineering in Z. mobilis including plasmid curing, gene deletion and insertion, as well as nucleotide substitution, which can also be employed for metabolic engineering to help divert the carbon flux from ethanol production to other products such as lactate demonstrated in this work. The CRISPR-Cas12a system established in this study thus provides a versatile and powerful genome-editing tool in Z. mobilis for functional genomic research, strain improvement, as well as synthetic microbial chassis development for economic biochemical production.

Genome Wide Phosphoproteome Analysis of Zymomonas mobilis Under Anaerobic, Aerobic, and N2-Fixing Conditions

Protein phosphorylation is a post-translational modification with widespread regulatory roles in both eukaryotes and prokaryotes. Using mass spectrometry, we performed a genome wide investigation of protein phosphorylation in the non-model organism and biofuel producer Zymomonas mobilis under anaerobic, aerobic, and N₂-fixing conditions. Our phosphoproteome analysis revealed 125 unique phosphorylated proteins, belonging to major pathways such as glycolysis, TCA cycle, electron transport, nitrogen metabolism, and protein synthesis. Quantitative analysis revealed significant and widespread changes in protein phosphorylation across growth conditions. For example, we observed increased phosphorylation of nearly all glycolytic enzymes and a large fraction of ribosomal proteins during aerobic and N₂-fixing conditions. We also observed substantial changes in the phosphorylation status of enzymes and regulatory proteins involved in nitrogen fixation and ammonia assimilation during N₂-fixing conditions, including nitrogenase, the Rnf electron transport complex, the transcription factor NifA, GS-GOGAT cycle enzymes, and the P_II regulatory protein. This suggested that protein phosphorylation may play an important role at regulating all aspects of nitrogen metabolism in Z. mobilis. This study provides new knowledge regarding the specific pathways and cellular processes that may be regulated by protein phosphorylation in this important industrial organism and provides a useful road map for future experiments that investigate the physiological role of specific phosphorylation events in Z. mobilis.

Genetic Regulation of Liver Metabolites and Transcripts Linking to Biochemical-Clinical Parameters

Given the central metabolic role of the liver, hepatic metabolites and transcripts reflect the organismal physiological state. Biochemical-clinical plasma biomarkers, hepatic metabolites, transcripts, and single nucleotide polymorphism (SNP) genotypes of some 300 pigs were integrated by weighted correlation networks and genome-wide association analyses. Network-based approaches of transcriptomic and metabolomics data revealed linked of transcripts and metabolites of the pentose phosphate pathway (PPP). This finding was evidenced by using a NADP/NADPH assay and HDAC4 and G6PD transcript quantification with the latter coding for first limiting enzyme of this pathway and by RNAi knockdown experiments of HDAC4. Other transcripts including ARG2 and SLC22A7 showed link to amino acids and biomarkers. The amino acid metabolites were linked with transcripts of immune or acute phase response signaling, whereas the carbohydrate metabolites were highly enrich in cholesterol biosynthesis transcripts. Genome-wide association analyses revealed 180 metabolic quantitative trait loci (mQTL) (p < 10^-4). Trans-4-hydroxy-L-proline (p = 6 × 10^-9), being strongly correlated with plasma creatinine (CREA), showed strongest association with SNPs on chromosome 6 that had pleiotropic effects on PRODH2 expression as revealed by multivariate analysis. Consideration of shared marker association with biomarkers, metabolites, and transcripts revealed 144 SNPs associated with 44 metabolites and 69 transcripts that are correlated with each other, representing 176 mQTL and expression quantitative trait loci (eQTL). This is the first work to report genetic variants associated with liver metabolite and transcript levels as well as blood biochemical-clinical parameters in a healthy porcine model. The identified associations provide links between variation at the genome, transcriptome, and metabolome level molecules with clinically relevant phenotypes. This approach has the potential to detect novel biomarkers displaying individual variation and promoting predictive biology in medicine and animal breeding.

New technologies provide more metabolic engineering strategies for bioethanol production in Zymomonas mobilis

Bioethanol has been considered as a potentially renewable energy source, and metabolic engineering plays an important role in the production of biofuels. As an efficient ethanol-producing bacterium, Zymomonas mobilis has garnered special attention due to its high sugar uptake, ethanol yield, and tolerance. Different metabolic engineering strategies have been used to establish new metabolic pathways for Z. mobilis to broaden its substrate range, remove competing pathways, and enhance its tolerance to ethanol and lignocellulosic hydrolysate inhibitors. Recent advances in omics technology, computational modeling and simulation, system biology, and synthetic biology contribute to the efficient re-design and manipulation of microbes via metabolic engineering at the whole-cell level. In this review, we summarize the progress of some new technologies used for metabolic engineering to improve bioethanol production and tolerance in Z. mobilis. Some successful examples of metabolic engineering used to develop strains for ethanol production are described in detail. Lastly, some important strategies for future metabolic engineering efforts are also highlighted.

Prediction and characterization of promoters and ribosomal binding sites of Zymomonas mobilis in system biology era

BACKGROUND

Zymomonas mobilis is a model bacterial ethanologen with many systems biology studies reported. Besides lignocellulosic ethanol production, Z. mobilis has been developed as a platform for biochemical production through metabolic engineering. However, identification and rigorous understanding of the genetic origins of cellular function, especially those based in non-coding region of DNA, such as promoters and ribosomal binding sites (RBSs), are still in its infancy. This knowledge is crucial for the effective application of Z. mobilis to new industrial applications of biotechnology for fuels and chemicals production.

RESULTS

In this study, we explored the possibility to systematically predict the strength of promoters based on systems biology datasets. The promoter strength was clustered based on the expression values of downstream genes (or proteins) from systems biology studies including microarray, RNA-Seq and proteomics. Candidate promoters with different strengths were selected for further characterization, which include 19 strong, nine medium, and ten weak ones. A dual reporter-gene system was developed which included appropriate reporter genes. These are the opmCherry reporter gene driven by the constitutive PlacUV5 promoter for calibration, and EGFP reporter gene driven by candidate promoters for quantification. This dual reporter-gene system was confirmed using the inducible promoter, Ptet, which was used to determine the strength of these predicted promoters with different strengths. In addition, the dual reporter-gene system was applied to determine four synthetic RBSs with different translation initiation rates based on the prediction from bioinformatics server RBS calculator. Our results showed that the correlations between the prediction and experimental results for the promoter and RBS strength are relatively high, with R ² values more than 0.7 and 0.9, respectively.

CONCLUSIONS

This study not only identified and characterized 38 promoters and four RBSs with different strengths for future metabolic engineering in Z. mobilis, but also established a flow cytometry-based dual reporter-gene system to characterize genetic elements including, but not limited to the promoters and RBSs studied in this work. This study also suggested the feasibility of predicting and selecting candidate genetic elements based on omics datasets and bioinformatics tools. Moreover, the dual reporter-gene system developed in this study can be utilized to characterize other genetic elements of Z. mobilis, which can also be applied to other microorganisms.

Systems-Level Analysis of Oxygen Exposure in Zymomonas mobilis: Implications for Isoprenoid Production

Zymomonas mobilis is an aerotolerant anaerobe and prolific ethanologen with attractive characteristics for industrial bioproduct generation. However, there is currently insufficient knowledge of the impact that environmental factors have on flux through industrially relevant biosynthetic pathways. Here, we examined the effect of oxygen exposure on metabolism and gene expression in Z. mobilis by combining targeted metabolomics, mRNA sequencing, and shotgun proteomics. We found that exposure to oxygen profoundly influenced metabolism, inducing both transient metabolic bottlenecks and long-term metabolic remodeling. In particular, oxygen induced a severe but temporary metabolic bottleneck in the methyl erythritol 4-phosphate pathway for isoprenoid biosynthesis caused by oxidative damage to the iron-sulfur cofactors of the final two enzymes in the pathway. This bottleneck was resolved with minimal changes in expression of isoprenoid biosynthetic enzymes. Instead, it was associated with pronounced upregulation of enzymes related to iron-sulfur cluster maintenance and biogenesis (i.e., flavodoxin reductase and the suf operon). We also detected major changes in glucose utilization in the presence of oxygen. Specifically, we observed increased gluconate production following exposure to oxygen, accounting for 18% of glucose uptake. Our results suggest that under aerobic conditions, electrons derived from the oxidation of glucose to gluconate are diverted to the electron transport chain, where they can minimize oxidative damage by reducing reactive oxygen species such as H2O2. This model is supported by the simultaneous upregulation of three membrane-bound dehydrogenases, cytochrome c peroxidase, and a cytochrome bd oxidase following exposure to oxygen. IMPORTANCE Microbially generated biofuels and bioproducts have the potential to provide a more environmentally sustainable alternative to fossil-fuel-derived products. In particular, isoprenoids, a diverse class of natural products, are chemically suitable for use as high-grade transport fuels and other commodity molecules. However, metabolic engineering for increased production of isoprenoids and other bioproducts is limited by an incomplete understanding of factors that control flux through biosynthetic pathways. Here, we examined the native regulation of the isoprenoid biosynthetic pathway in the biofuel producer Zymomonas mobilis. We leveraged oxygen exposure as a means to perturb carbon flux, allowing us to observe the formation and resolution of a metabolic bottleneck in the pathway. Our multi-omics analysis of this perturbation enabled us to identify key auxiliary enzymes whose expression correlates with increased production of isoprenoid precursors, which we propose as potential targets for future metabolic engineering.

Engineered Zymomonas mobilis tolerant to acetic acid and low pH via multiplex atmospheric and room temperature plasma mutagenesis

BACKGROUND

Cellulosic biofuels are sustainable compared to fossil fuels. However, inhibitors, such as acetic acid generated during lignocellulose pretreatment and hydrolysis, would significantly inhibit microbial fermentation efficiency. Microbial mutants able to tolerate high concentration of acetic acid are needed urgently to alleviate this inhibition.

RESULTS

Zymomonas mobilis mutants AQ8-1 and AC8-9 with enhanced tolerance against acetic acid were generated via a multiplex atmospheric and room temperature plasma (mARTP) mutagenesis. The growth and ethanol productivity of AQ8-1 and AC8-9 were both improved in the presence of 5.0-8.0 g/L acetic acid. Ethanol yield reached 84% of theoretical value in the presence of 8.0 g/L acetic acid (~ pH 4.0). Furthermore, a mutant tolerant to pH 3.5, named PH1-29, was generated via the third round of ARTP mutagenesis. PH1-29 showed enhanced growth and ethanol production under both sterilized/unsterilized conditions at pH 4.0 or 3.5. Intracellular NAD levels revealed that mARTP mutants could modulate NADH/NAD⁺ ratio to respond to acetic acid and low pH stresses. Moreover, genomic re-sequencing revealed that eleven single nucleic variations (SNVs) were likely related to acetic acid and low pH tolerance. Most SNVs were targeted in regions between genes ZMO0952 and ZMO0956, ZMO0152 and ZMO0153, and ZMO0373 and ZMO0374.

CONCLUSIONS

The multiplex mutagenesis strategy mARTP was efficient for enhancing the tolerance in Z. mobilis. The ARTP mutants generated in this study could serve as potential cellulosic ethanol producers.

Aerobic catabolism and respiratory lactate bypass in Ndh-negative Zymomonas mobilis

Ability to ferment in the presence of oxygen increases the robustness of bioprocesses and opens opportunity for novel industrial setups. The ethanologenic bacterium Zymomonas mobilis performs rapid and efficient anaerobic ethanol fermentation, yet its respiratory NADH dehydrogenase (Ndh)-deficient strain (ndh-) is known to produce ethanol with high yield also under oxic conditions. Compared to the wild type, it has a lower rate of oxygen consumption, and an increased expression of the respiratory lactate dehydrogenase (Ldh). Here we present a quantitative study of the product spectrum and carbon balance for aerobically growing ndh-. Ldh-deficient and Ldh-overexpressing ndh- strains were constructed and used to examine the putative role of the respiratory lactate bypass for aerobic growth and production. We show that aerobically growing ndh- strains perform fermentative metabolism with a near-maximum ethanol yield, irrespective of their Ldh expression background. Yet, Ldh activity strongly affects the aerobic product spectrum in glucose-consuming non-growing cells. Also, Ldh-deficiency hampers growth at elevated temperature (42 °C) and delays the restart of growth after 10-15 h of aerobic starvation.

Proteomic and metabolomic analysis of the cellular biomarkers related to inhibitors tolerance in Zymomonas mobilis ZM4

BACKGROUND

Toxic compounds present in both the hydrolysate and pyrolysate of lignocellulosic biomass severely hinder the further conversion of lignocellulose-derived fermentable sugars into useful chemicals by common biocatalysts like Zymomonas mobilis, which has remarkable advantages over yeast. Although the extra detoxification treatment prior to fermentation process can help biocatalysts to eliminate the inhibitory environment, it is not environment friendly and cost effective for industrial application. As also reported by previous studies, an ideal and holistic approach to solve this issue is to develop microbial strains with inhibitor tolerance. However, previously engineered strains had the limitation that they could not cope well with the synergistic effect of multiple inhibitors as they are resistant only to a single inhibitor. Hence, understanding the universal cellular responses of Z. mobilis to various inhibitors may guide the designing of rational strategies to obtain more robust engineered strains for biofuel production from lignocellulosic biomass.

RESULTS

Quantitative proteomics and metabolomics approaches were used to determine the cellular responses of Z. mobilis ZM4 to representative biomass-derived inhibitors like formic acid, acetic acid, furfural, 5-hydroxymethylfurfural, and phenol. The differentially expressed proteins identified under the challenge of single and combined inhibitors were involved in cell wall/membrane biogenesis, energy production, DNA replication, DNA recombination, DNA repair, DNA transcription, RNA translation, posttranslational modification, biosynthesis of amino acids, central carbon metabolism, etc. Metabolomics analysis showed that the up- or down-regulation pattern of metabolites was changed consistently with that of relevant proteins.

CONCLUSION

Fifteen up-regulated proteins (e.g., Isopropylmalate isomerase LeuC, transcription-repair-coupling factor Mfd, and phosphoglucose isomerase PGI) and thirteen down-regulated proteins (e.g., TonB-dependent transporter ZMO1522, transcription termination factor Rho, and S1/P1 nuclease ZMO0127) were identified as candidate proteins related to all the stress conditions, implying that these proteins are potential biomarkers for the improvement of Z. mobilis ZM4 to resist complex biomass-derived inhibitors. These data can be used to generate a database of inhibitor-tolerance biomarkers, which could provide a basis for engineering Z. mobilis that would be able to grow in the presence of multiple inhibitors and directly ferment the biomass-derived sugars into biofuels.

Comprehensive exploration of the enzymes catalysing oxygen-involved reactions and COGs relevant to bacterial oxygen utilization

To better understand the mechanisms of bacterial adaptation in oxygen environments, we explored the aerobic living-associated genes in bacteria by comparing Clusters of Orthologous Groups of proteins' (COGs) frequencies and gene expression analyses and 38 COGs were detected at significantly higher frequencies (p-value less than 1e-6) in aerobes than in anaerobes. Differential expression analyses between two conditions further narrowed the prediction to 27 aerobe-specific COGs. Then, we annotated the enzymes associated with these COGs. Literature review revealed that 14 COGs contained enzymes catalysing oxygen-involved reactions or products involved in aerobic pathways, suggesting their important roles for survival in aerobic environments. Additionally, protein-protein interaction analyses and step length comparisons of metabolic networks suggested that the other 13 COGs may function relevantly with the 14 enzymes-corresponding COGs, indicating that these genes may be highly associated with oxygen utilization. Phylogenetic and evolutionary analyses showed that the 27 COGs did not have similar trees, and all suffered purifying selection pressures. The divergent times of species containing or lacking aerobic COGs validated that the appearing time of oxygen-utilizing gene was approximately 2.80 Gyr ago. In addition to help better understand oxygen adaption, our method may be extended to identify genes relevant to other living environments.

Oleaginicity of the yeast strain Saccharomyces cerevisiae D5A

BACKGROUND

The model yeast, Saccharomyces cerevisiae, is not known to be oleaginous. However, an industrial wild-type strain, D5A, was shown to accumulate over 20% storage lipids from glucose when growth is nitrogen-limited compared to no more than 7% lipid accumulation without nitrogen stress.

METHODS AND RESULTS

To elucidate the mechanisms of S. cerevisiae D5A oleaginicity, we compared physiological and metabolic changes; as well as the transcriptional profiles of the oleaginous industrial strain, D5A, and a non-oleaginous laboratory strain, BY4741, under normal and nitrogen-limited conditions using analytic techniques and next-generation sequencing-based RNA-Seq transcriptomics. Transcriptional levels for genes associated with fatty acid biosynthesis, nitrogen metabolism, amino acid catabolism, as well as the pentose phosphate pathway and ethanol oxidation in central carbon (C) metabolism, were up-regulated in D5A during nitrogen deprivation. Despite increased carbon flux to lipids, most gene-encoding enzymes involved in triacylglycerol (TAG) assembly were expressed at similar levels regardless of the varying nitrogen concentrations in the growth media and strain backgrounds. Phospholipid turnover also contributed to TAG accumulation through increased precursor production with the down-regulation of subsequent phospholipid synthesis steps. Our results also demonstrated that nitrogen assimilation via the glutamate-glutamine pathway and amino acid metabolism, as well as the fluxes of carbon and reductants from central C metabolism, are integral to the general oleaginicity of D5A, which resulted in the enhanced lipid storage during nitrogen deprivation.

CONCLUSION

This work demonstrated the disequilibrium and rebalance of carbon and nitrogen contribution to the accumulation of lipids in the oleaginous yeast S. cerevisiae D5A. Rather than TAG assembly from acyl groups, the major switches for the enhanced lipid accumulation of D5A (i.e., fatty acid biosynthesis) are the increases of cytosolic pools of acetyl-CoA and NADPH, as well as alternative nitrogen assimilation.

Advances and prospects in metabolic engineering of Zymomonas mobilis

Biorefinery of biomass-based biofuels and biochemicals by microorganisms is a competitive alternative of traditional petroleum refineries. Zymomonas mobilis is a natural ethanologen with many desirable characteristics, which makes it an ideal industrial microbial biocatalyst for commercial production of desirable bioproducts through metabolic engineering. In this review, we summarize the metabolic engineering progress achieved in Z. mobilis to expand its substrate and product ranges as well as to enhance its robustness against stressful conditions such as inhibitory compounds within the lignocellulosic hydrolysates and slurries. We also discuss a few metabolic engineering strategies that can be applied in Z. mobilis to further develop it as a robust workhorse for economic lignocellulosic bioproducts. In addition, we briefly review the progress of metabolic engineering in Z. mobilis related to the classical synthetic biology cycle of "Design-Build-Test-Learn", as well as the progress and potential to develop Z. mobilis as a model chassis for biorefinery practices in the synthetic biology era.

Complete genome sequence and the expression pattern of plasmids of the model ethanologen Zymomonas mobilis ZM4 and its xylose-utilizing derivatives 8b and 2032

BACKGROUND

Zymomonas mobilis is a natural ethanologen being developed and deployed as an industrial biofuel producer. To date, eight Z. mobilis strains have been completely sequenced and found to contain 2-8 native plasmids. However, systematic verification of predicted Z. mobilis plasmid genes and their contribution to cell fitness has not been hitherto addressed. Moreover, the precise number and identities of plasmids in Z. mobilis model strain ZM4 have been unclear. The lack of functional information about plasmid genes in ZM4 impedes ongoing studies for this model biofuel-producing strain.

RESULTS

In this study, we determined the complete chromosome and plasmid sequences of ZM4 and its engineered xylose-utilizing derivatives 2032 and 8b. Compared to previously published and revised ZM4 chromosome sequences, the ZM4 chromosome sequence reported here contains 65 nucleotide sequence variations as well as a 2400-bp insertion. Four plasmids were identified in all three strains, with 150 plasmid genes predicted in strain ZM4 and 2032, and 153 plasmid genes predicted in strain 8b due to the insertion of heterologous DNA for expanded substrate utilization. Plasmid genes were then annotated using Blast2GO, InterProScan, and systems biology data analyses, and most genes were found to have apparent orthologs in other organisms or identifiable conserved domains. To verify plasmid gene prediction, RNA-Seq was used to map transcripts and also compare relative gene expression under various growth conditions, including anaerobic and aerobic conditions, or growth in different concentrations of biomass hydrolysates. Overall, plasmid genes were more responsive to varying hydrolysate concentrations than to oxygen availability. Additionally, our results indicated that although all plasmids were present in low copy number (about 1-2 per cell), the copy number of some plasmids varied under specific growth conditions or due to heterologous gene insertion.

CONCLUSIONS

The complete genome of ZM4 and two xylose-utilizing derivatives is reported in this study, with an emphasis on identifying and characterizing plasmid genes. Plasmid gene annotation, validation, expression levels at growth conditions of interest, and contribution to host fitness are reported for the first time.

Identification and Characterization of 5' Untranslated Regions (5'UTRs) in Zymomonas mobilis as Regulatory Biological Parts

Regulatory RNA regions within a transcript, particularly in the 5' untranslated region (5'UTR), have been shown in a variety of organisms to control the expression levels of these mRNAs in response to various metabolites or environmental conditions. Considering the unique tolerance of Zymomonas mobilis to ethanol and the growing interest in engineering microbial strains with enhanced tolerance to industrial inhibitors, we searched natural cis-regulatory regions in this microorganism using transcriptomic data and bioinformatics analysis. Potential regulatory 5'UTRs were identified and filtered based on length, gene function, relative gene counts, and conservation in other organisms. An in vivo fluorescence-based screening system was developed to confirm the responsiveness of 36 5'UTR candidates to ethanol, acetate, and xylose stresses. UTR_ZMO0347 (5'UTR of gene ZMO0347 encoding the RNA binding protein Hfq) was found to down-regulate downstream gene expression under ethanol stress. Genomic deletion of UTR_ZMO0347 led to a general decrease of hfq expression at the transcript level and increased sensitivity for observed changes in Hfq expression at the protein level. The role of UTR_ZMO0347 and other 5'UTRs gives us insight into the regulatory network of Z. mobilis in response to stress and unlocks new strategies for engineering robust industrial strains as well as for harvesting novel responsive regulatory biological parts for controllable gene expression platforms in this organism.

Generation and Characterization of Acid Tolerant Fibrobacter succinogenes S85

Microorganisms are key components for plant biomass breakdown within rumen environments. Fibrobacter succinogenes have been identified as being active and dominant cellulolytic members of the rumen. In this study, F. succinogenes type strain S85 was adapted for steady state growth in continuous culture at pH 5.75 and confirmed to grow in the range of pH 5.60–5.65, which is lower than has been reported previously. Wild type and acid tolerant strains digested corn stover with equal efficiency in batch culture at low pH. RNA-seq analysis revealed 268 and 829 genes were differentially expressed at pH 6.10 and 5.65 compared to pH 6.70, respectively. Resequencing analysis identified seven single nucleotide polymorphisms (SNPs) in the sufD, yidE, xylE, rlmM, mscL and dosC genes of acid tolerant strains. Due to the absence of a F. succinogenes genetic system, homologues in Escherichia coli were mutated and complemented and the resulting strains were assayed for acid survival. Complementation with wild-type or acid tolerant F. succinogenes sufD restored E. coli wild-type levels of acid tolerance, suggesting a possible role in acid homeostasis. Recent genetic engineering developments need to be adapted and applied in F. succinogenes to further our understanding of this bacterium.

Transgenic soybean overexpressing GmSAMT1 exhibits resistance to multiple-HG types of soybean cyst nematode Heterodera glycines

Soybean (Glycine max (L.) Merr.) salicylic acid methyl transferase (GmSAMT1) catalyses the conversion of salicylic acid to methyl salicylate. Prior results showed that when GmSAMT1 was overexpressed in transgenic soybean hairy roots, resistance is conferred against soybean cyst nematode (SCN), Heterodera glycines Ichinohe. In this study, we produced transgenic soybean overexpressing GmSAMT1 and characterized their response to various SCN races. Transgenic plants conferred a significant reduction in the development of SCN HG type 1.2.5.7 (race 2), HG type 0 (race 3) and HG type 2.5.7 (race 5). Among transgenic lines, GmSAMT1 expression in roots was positively associated with SCN resistance. In some transgenic lines, there was a significant decrease in salicylic acid titer relative to control plants. No significant seed yield differences were observed between transgenics and control soybean plants grown in one greenhouse with 22 °C day/night temperature, whereas transgenic soybean had higher yield than controls grown a warmer greenhouse (27 °C day/23 °C night) temperature. In a 1-year field experiment in Knoxville, TN, there was no significant difference in seed yield between the transgenic and nontransgenic soybean under conditions with negligible SCN infection. We hypothesize that GmSAMT1 expression affects salicylic acid biosynthesis, which, in turn, attenuates SCN development, without negative consequences to soybean yield or other morphological traits. Thus, we conclude that GmSAMT1 overexpression confers broad resistance to multiple SCN races, which would be potentially applicable to commercial production.

Zymomonas mobilis as a model system for production of biofuels and biochemicals

Zymomonas mobilis is a natural ethanologen with many desirable industrial biocatalyst characteristics. In this review, we will discuss work to develop Z. mobilis as a model system for biofuel production from the perspectives of substrate utilization, development for industrial robustness, potential product spectrum, strain evaluation and fermentation strategies. This review also encompasses perspectives related to classical genetic tools and emerging technologies in this context.

Metabolic engineering of Zymomonas mobilis for 2,3-butanediol production from lignocellulosic biomass sugars

BACKGROUND

To develop pathways for advanced biofuel production, and to understand the impact of host metabolism and environmental conditions on heterologous pathway engineering for economic advanced biofuels production from biomass, we seek to redirect the carbon flow of the model ethanologen Zymomonas mobilis to produce desirable hydrocarbon intermediate 2,3-butanediol (2,3-BDO). 2,3-BDO is a bulk chemical building block, and can be upgraded in high yields to gasoline, diesel, and jet fuel.

RESULTS

2,3-BDO biosynthesis pathways from various bacterial species were examined, which include three genes encoding acetolactate synthase, acetolactate decarboxylase, and butanediol dehydrogenase. Bioinformatics analysis was carried out to pinpoint potential bottlenecks for high 2,3-BDO production. Different combinations of 2,3-BDO biosynthesis metabolic pathways using genes from different bacterial species have been constructed. Our results demonstrated that carbon flux can be deviated from ethanol production into 2,3-BDO biosynthesis, and all three heterologous genes are essential to efficiently redirect pyruvate from ethanol production for high 2,3-BDO production in Z. mobilis. The down-selection of best gene combinations up to now enabled Z. mobilis to reach the 2,3-BDO production of more than 10 g/L from glucose and xylose, as well as mixed C6/C5 sugar streams derived from the deacetylation and mechanical refining process.

CONCLUSIONS

This study confirms the value of integrating bioinformatics analysis and systems biology data during metabolic engineering endeavors, provides guidance for value-added chemical production in Z. mobilis, and reveals the interactions between host metabolism, oxygen levels, and a heterologous 2,3-BDO biosynthesis pathway. Taken together, this work provides guidance for future metabolic engineering efforts aimed at boosting 2,3-BDO titer anaerobically.

Reconstruction of a charge balanced genome-scale metabolic model to study the energy-uncoupled growth of Zymomonas mobilis ZM1

Zymomonas mobilisis an ethanologenic bacterium and is known to be an example microorganism with energy-uncoupled growth. The reconstructed metabolic model indicate that resistance to intracellular pH reduction could be the main reason for uncoupled growth.

Genome-scale resources for Thermoanaerobacterium saccharolyticum

Background

Thermoanaerobacterium saccharolyticum is a hemicellulose-degrading thermophilic anaerobe that was previously engineered to produce ethanol at high yield. A major project was undertaken to develop this organism into an industrial biocatalyst, but the lack of genome information and resources were recognized early on as a key limitation.

Results

Here we present a set of genome-scale resources to enable the systems level investigation and development of this potentially important industrial organism. Resources include a complete genome sequence for strain JW/SL-YS485, a genome-scale reconstruction of metabolism, tiled microarray data showing transcription units, mRNA expression data from 71 different growth conditions or timepoints and GC/MS-based metabolite analysis data from 42 different conditions or timepoints. Growth conditions include hemicellulose hydrolysate, the inhibitors HMF, furfural, diamide, and ethanol, as well as high levels of cellulose, xylose, cellobiose or maltodextrin. The genome consists of a 2.7 Mbp chromosome and a 110 Kbp megaplasmid. An active prophage was also detected, and the expression levels of CRISPR genes were observed to increase in association with those of the phage. Hemicellulose hydrolysate elicited a response of carbohydrate transport and catabolism genes, as well as poorly characterized genes suggesting a redox challenge. In some conditions, a time series of combined transcription and metabolite measurements were made to allow careful study of microbial physiology under process conditions. As a demonstration of the potential utility of the metabolic reconstruction, the OptKnock algorithm was used to predict a set of gene knockouts that maximize growth-coupled ethanol production. The predictions validated intuitive strain designs and matched previous experimental results.

Conclusion

These data will be a useful asset for efforts to develop T. saccharolyticum for efficient industrial production of biofuels. The resources presented herein may also be useful on a comparative basis for development of other lignocellulose degrading microbes, such as Clostridium thermocellum.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0159-x) contains supplementary material, which is available to authorized users.