Global transcriptional response of solvent‐sensitive and solvent‐tolerant Pseudomonas putida strains exposed to toluene

Pseudomonas putida strains are generally recognized as solvent tolerant, exhibiting varied sensitivity to organic solvents. Pan‐genome analysis has revealed that 30% of genes belong to the core‐genome of Pseudomonas. Accessory and unique genes confer high degree of adaptability and capabilities for the degradation and synthesis of a wide range of chemicals. For the use of these microbes in bioremediation and biocatalysis, it is critical to understand the mechanisms underlying these phenotypic differences. In this study, RNA‐seq analysis compared the short‐ and long‐term responses of the toluene‐sensitive KT2440 strain and the highly tolerant DOT‐T1E strain. The sensitive strain activates a larger number of genes in a higher magnitude than DOT‐T1E. This is expected because KT2440 bears one toluene tolerant pump, while DOT‐T1E encodes three of these pumps. Both strains activate membrane modifications to reduce toluene membrane permeability. The KT2440 strain activates the TCA cycle to generate energy, while avoiding energy‐intensive processes such as flagellar biosynthesis. This suggests that KT2440 responds to toluene by focusing on survival mechanisms. The DOT‐T1E strain activates toluene degradation pathways, using toluene as source of energy. Among the unique genes encoded by DOT‐T1E is a 70 kb island composed of genes of unknown function induced in response to toluene.     

CRISPR interference–mediated metabolic engineering of Corynebacterium glutamicum for homo‐butyrate production

Combinatorial metabolic engineering enabled the development of efficient microbial cell factories for modulating gene expression to produce desired products. Here, we report the combinatorial metabolic engineering of Corynebacterium glutamicum to produce butyrate by introducing a synthetic butyrate pathway including phosphotransferase and butyrate kinase reactions and repressing the essential acn gene‐encoding aconitase, which has been targeted for downregulation in a genome‐scale model. An all‐in‐one clustered regularly interspaced short palindromic repeats interference system for C. glutamicum was used for tunable downregulation of acn in an engineered strain, where by‐product‐forming reactions were deleted and the synthetic butyrate pathway was inserted, resulting in butyrate production (0.52 ± 0.02 g/L). Subsequently, biotin limitation enabled the engineered strain to produce butyrate (0.58 ± 0.01 g/L) without acetate formation for the entire duration of the culture. These results demonstrate the potential homo‐production of butyrate using engineered C. glutamicum. This method can also be applied to other industrial microorganisms.     

Genome‐scale metabolic reconstruction and constraint‐based modelling of the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125

The Antarctic strain Pseudoalteromonas haloplanktis TAC125 is one of the model organisms of cold‐adapted bacteria and is currently exploited as a new alternative expression host for numerous biotechnological applications. Here, we investigated several metabolic features of this strain through in silico modelling and functional integration of –omics data. A genome‐scale metabolic model of P. haloplanktis TAC125 was reconstructed, encompassing information on 721 genes, 1133 metabolites and 1322 reactions. The predictive potential of this model was validated against a set of experimentally determined growth rates and a large dataset of growth phenotypic data. Furthermore, evidence synthesis from proteomics, phenomics, physiology and metabolic modelling data revealed possible drawbacks of cold‐dependent changes in gene expression on the overall metabolic network of P. haloplanktis TAC125. These included, for example, variations in its central metabolism, amino acid degradation and fatty acid biosynthesis. The genome‐scale metabolic model described here is the first one reconstructed so far for an Antarctic microbial strain. It allowed a system‐level investigation of variations in cellular metabolic fluxes following a temperature downshift. It represents a valuable platform for further investigations on P. haloplanktis TAC125 cellular functional states and for the design of more focused strategies for its possible biotechnological exploitation.     

Efflux pump‐deficient mutants as a platform to search for microbes that produce antibiotics

Pseudomonas putida DOT‐T1E‐18 is a strain deficient in the major antibiotic efflux pump (TtgABC) that exhibits an overall increased susceptibility to a wide range of drugs when compared with the wild‐type strain. We used this strain as a platform to search for microbes able to produce antibiotics that inhibit growth. A collection of 2400 isolates from soil, sediments and water was generated and a drop assay developed to identify, via growth inhibition halos, strains that prevent the growth of DOT‐T1E‐18 on solid Luria–Bertani plates. In this study, 35 different isolates that produced known and unknown antibiotics were identified. The most potent inhibitor of DOT‐T1E‐18 growth was an isolate named 250J that, through multi‐locus sequence analysis, was identified as a Pseudomonas sp. strain. Culture supernatants of 250J contain four different xantholysins that prevent growth of Gram‐positive bacteria, Gram‐negative and fungi. Two of the xantholysins were produced in higher concentrations and purified. Xantholysin A was effective against Bacillus, Lysinibacillus and Rhodococcus strains, and the effect against these microbes was enhanced when used in combination with other antibiotics such as ampicillin, gentamicin and kanamycin. Xantholysin C was also efficient against Gram‐positive bacteria and showed an interesting antimicrobial effect against Pseudomonas strains, and a synergistic inhibitory effect with ampicillin, chloramphenicol and gentamicin.     

Characterization of a group of pyrroloquinoline quinone‐dependent dehydrogenases that are involved in the conversion of L‐sorbose to 2‐Keto‐L‐gulonic acid in Ketogulonicigenium vulgare WSH‐001

Ketogulonicigenium vulgare WSH‐001 is an industrial strain used for vitamin C production. Based on genome sequencing and pathway analysis of the bacterium, some of its potential pyrroloquinoline quinone (PQQ)‐dependent dehydrogenases were predicted, including KVU_pmdA_0245, KVU_2142, KVU_2159, KVU_1366, KVU_0203, KVU_0095, and KVU_pmdB_0115. BLAST and function domain searches showed that enzymes encoded by these genes may act as putative PQQ‐dependent L ‐sorbose dehydrogenases (SDH) or L ‐sorbosone dehydrogenases (SNDH). To validate whether these dehydrogenases are PQQ‐dependent or not, these seven putative dehyrogenases were overexpressed in Escherichia coli BL21 (DE3) and purified for characterization. Biochemical and kinetic characterization of the purified proteins have led to the identification of seven enzymes that possess the ability to oxidize L ‐sorbose or L ‐sorbosone to varying degrees. In addition, the dehydrogenation of sorbose in K. vulgare is validated to be PQQ dependent, identification of these PQQ‐dependent dehydrogenases expanded the PQQ‐dependent dehydrogenase family. Besides, the optimal combination of enzymes that could more efficiently catalyze the conversion of sorbose to gulonic acid was proposed. These are important in supporting the development of metabolic engineering strategies and engineering of efficient strains for one‐step production of vitamin C in the future. © 2013 American Institute of Chemical Engineers Biotechnol. Prog., 29:1398–1404, 2013     

Agrobacterium‐mediated transformation of reed (Phragmites communis Trinius) using mature seed‐derived calli

Reed (Phragmites communis) is a potential bioenergy plant. We report on its first Agrobacterium‐mediated transformation using mature seed‐derived calli. The Agrobacterium strains LBA4404, EHA105, and GV3101, each harboring the binary vector pIG121Hm, were used to optimize T‐DNA delivery into the reed genome. Bacterial strain, cocultivation period and acetosyringone concentration significantly influenced the T‐DNA transfer. About 48% transient expression and 3.5% stable transformation were achieved when calli were infected with strain EHA105 for 10 min under 800 mbar negative pressure and cocultivated for 3 days in 200 μm acetosyringone containing medium. Putative transformants were selected in 25 mg l^?1 hygromycin B. PCR, and Southern blot analysis confirmed the presence of the transgenes and their stable integration. Independent transgenic lines contained one to three copies of the transgene. Transgene expression was validated by RT‐PCR and GUS staining of stems and leaves.     

MoNSTR‐seq,a restriction site‐associated DNA sequencing technique to characterize Agrobacterium‐mediated transfer‐DNA insertions in Phomopsis longicolla

Phomopsis longicolla (Hobbs) causes Phomopsis seed decay and stem lesions in soybean (Glycine max). In this study, a novel, high‐throughput adaptation of RAD‐seq termed MoNSTR‐seq (Mutation analysis via Next‐generation DNA Sequencing of T‐DNA Regions) was developed to determine the genomic location of T‐DNA insertions in P. longicolla mutants. Insertional mutants were created via Agrobacterium tumefaciens‐mediated transformation, and one mutant, strain PL343, was further investigated due to impaired stem lesion formation. Mutation analysis via Next‐generation DNA Sequencing of T‐DNA Regions, in which DNA libraries are created with two distinct restriction enzymes and customized adapters to simultaneously enrich both T‐DNA insertion borders, was developed to characterize the genomic lesion in strain PL343. MoNSTR‐seq successfully identified a T‐DNA insertion in the predicted promoter region of a gene encoding a cellobiose dehydrogenase (CDH1), and the position of the T‐DNA insertion in strain PL343 was confirmed by Sanger sequencing. Thus, MoNSTR‐seq represents an effective tool for molecular genetics in P. longicolla, and is readily adaptable for use in diverse fungal species.

Significance and Impact of the Study

This study describes MoNSTR‐seq (Mutation analysis via Next‐generation DNA Sequencing of T‐DNA Regions), an adaptation of restriction site‐associated DNA sequencing (RAD‐seq) to identify the position of transfer‐DNA (T‐DNA) insertions in the genome of Phomopsis longicolla, an important pathogen of soybean. The technique enables high‐throughput characterization of mutants generated via Agrobacterium tumefaciens‐mediated transformation (ATMT), thus accelerating gene discovery via forward genetics. This technique represents a significant advancement over existing approaches to characterize T‐DNA insertions in fungal genomes. With minor modifications, this technique could be easily adapted to taxonomically diverse fungal pathogens and additional mutagenesis cassettes.     

Biosynthesis of the anti‐diabetic metabolite montbretin A: glucosylation of the central intermediate mini‐MbA

Type 2 diabetes (T2D) affects over 320 million people worldwide. Healthy lifestyles, improved drugs and effective nutraceuticals are different components of a response against the growing T2D epidemic. The specialized metabolite montbretin A (MbA) is being developed for treatment of T2D and obesity due to its unique pharmacological activity as a highly effective and selective inhibitor of the human pancreatic α‐amylase. MbA is an acylated flavonol glycoside found in small amounts in montbretia (Crocosmia × crocosmiiflora) corms. MbA cannot be obtained in sufficient quantities for drug development from its natural source or by chemical synthesis. To overcome these limitations through metabolic engineering, we are investigating the genes and enzymes of MbA biosynthesis. We previously reported the first three steps of MbA biosynthesis from myricetin to myricetin 3‐O‐(6′‐O‐caffeoyl)‐glucosyl rhamnoside (mini‐MbA). Here, we describe the sequence of reactions from mini‐MbA to MbA, and the discovery and characterization of the gene and enzyme responsible for the glucosylation of mini‐MbA. The UDP‐dependent glucosyltransferase CcUGT3 (UGT703E1) catalyzes the 1,2‐glucosylation of mini‐MbA to produce myricetin 3‐O‐(glucosyl‐6′‐O‐caffeoyl)‐glucosyl rhamnoside. Co‐expression of CcUGT3 with genes for myricetin and mini‐MbA biosynthesis in Nicotiana benthamiana validated its biological function and expanded the set of genes available for metabolic engineering of MbA.     

The genome‐scale metabolic network of Ectocarpus siliculosus (EctoGEM): a resource to study brown algal physiology and beyond

Brown algae (stramenopiles) are key players in intertidal ecosystems, and represent a source of biomass with several industrial applications. Ectocarpus siliculosus is a model to study the biology of these organisms. Its genome has been sequenced and a number of post‐genomic tools have been implemented. Based on this knowledge, we report the reconstruction and analysis of a genome‐scale metabolic network for E. siliculosus, EctoGEM ( http://ectogem.irisa.fr ). This atlas of metabolic pathways consists of 1866 reactions and 2020 metabolites, and its construction was performed by means of an integrative computational approach for identifying metabolic pathways, gap filling and manual refinement. The capability of the network to produce biomass was validated by flux balance analysis. EctoGEM enabled the reannotation of 56 genes within the E. siliculosus genome, and shed light on the evolution of metabolic processes. For example, E. siliculosus has the potential to produce phenylalanine and tyrosine from prephenate and arogenate, but does not possess a phenylalanine hydroxylase, as is found in other stramenopiles. It also possesses the complete eukaryote molybdenum co‐factor biosynthesis pathway, as well as a second molybdopterin synthase that was most likely acquired via horizontal gene transfer from cyanobacteria by a common ancestor of stramenopiles. EctoGEM represents an evolving community resource to gain deeper understanding of the biology of brown algae and the diversification of physiological processes. The integrative computational method applied for its reconstruction will be valuable to set up similar approaches for other organisms distant from biological benchmark models.     

A comprehensive genome‐scale reconstruction of Escherichia coli metabolism—2011

The initial genome‐scale reconstruction of the metabolic network of Escherichia coli K‐12 MG1655 was assembled in 2000. It has been updated and periodically released since then based on new and curated genomic and biochemical knowledge. An update has now been built, named iJO1366, which accounts for 1366 genes, 2251 metabolic reactions, and 1136 unique metabolites. iJO1366 was (1) updated in part using a new experimental screen of 1075 gene knockout strains, illuminating cases where alternative pathways and isozymes are yet to be discovered, (2) compared with its predecessor and to experimental data sets to confirm that it continues to make accurate phenotypic predictions of growth on different substrates and for gene knockout strains, and (3) mapped to the genomes of all available sequenced E. coli strains, including pathogens, leading to the identification of hundreds of unannotated genes in these organisms. Like its predecessors, the iJO1366 reconstruction is expected to be widely deployed for studying the systems biology of E. coli and for metabolic engineering applications.     

RNA‐guided Cas9 as an in vivo desired‐target mutator in maize

The RNA‐guided Cas9 system is a versatile tool for genome editing. Here, we established a RNA‐guided endonuclease (RGEN) system as an in vivo desired‐target mutator (DTM) in maize to reduce the linkage drag during breeding procedure, using the LIGULELESS1 (LG1) locus as a proof‐of‐concept. Our system showed 51.5%–91.2% mutation frequency in T0 transgenic plants. We then crossed the T1 plants stably expressing DTM with six diverse recipient maize lines and found that 11.79%–28.71% of the plants tested were mutants induced by the DTM effect. Analysis of successive F2 plants indicated that the mutations induced by the DTM effect were largely heritable. Moreover, DTM‐generated hybrids had significantly smaller leaf angles that were reduced more than 50% when compared with that of the wild type. Planting experiments showed that DTM‐generated maize plants can be grown with significantly higher density and hence greater yield potential. Our work demonstrate that stably expressed RGEN could be implemented as an in vivoDTM to rapidly generate and spread desired mutations in maize through hybridization and subsequent backcrossing, and hence bypassing the linkage drag effect in convention introgression methodology. This proof‐of‐concept experiment can be a potentially much more efficient breeding strategy in crops employing the RNA‐guided Cas9 genome editing.     

Chemical and genomic evolution of enzyme-catalyzed reaction networks

There is a tendency that a unit of enzyme genes in an operon-like structure in the prokaryotic genome encodes enzymes that catalyze a series of consecutive reactions in a metabolic pathway. Our recent analysis shows that this and other genomic units correspond to chemical units reflecting chemical logic of organic reactions. From all known metabolic pathways in the KEGG database we identified chemical units, called reaction modules, as the conserved sequences of chemical structure transformation patterns of small molecules. The extracted patterns suggest co-evolution of genomic units and chemical units. While the core of the metabolic network may have evolved with mechanisms involving individual enzymes and reactions, its extension may have been driven by modular units of enzymes and reactions.     

Metabolic reconstruction identifies strain‐specific regulation of virulence in Toxoplasma gondii

Increasingly, metabolic potential is proving to be a critical determinant governing a pathogen's virulence as well as its capacity to expand its host range. To understand the potential contribution of metabolism to strain‐specific infectivity differences, we present a constraint‐based metabolic model of the opportunistic parasite, Toxoplasma gondii. Dominated by three clonal strains (Type I, II, and III demonstrating distinct virulence profiles), T. gondii exhibits a remarkably broad host range. Integrating functional genomic data, our model (which we term as iCS382) reveals that observed strain‐specific differences in growth rates are driven by altered capacities for energy production. We further predict strain‐specific differences in drug susceptibilities and validate one of these predictions in a drug‐based assay, with a Type I strain demonstrating resistance to inhibitors that are effective against a Type II strain. We propose that these observed differences reflect an evolutionary strategy that allows the parasite to extend its host range, as well as result in a subsequent partitioning into discrete strains that display altered virulence profiles across different hosts, different organs, and even cell types.     

A genome‐scale metabolic network reconstruction of tomato (Solanum lycopersicum L.) and its application to photorespiratory metabolism

Tomato (Solanum lycopersicum L.) has been studied extensively due to its high economic value in the market, and high content in health‐promoting antioxidant compounds. Tomato is also considered as an excellent model organism for studying the development and metabolism of fleshy fruits. However, the growth, yield and fruit quality of tomatoes can be affected by drought stress, a common abiotic stress for tomato. To investigate the potential metabolic response of tomato plants to drought, we reconstructed iHY3410, a genome‐scale metabolic model of tomato leaf, and used this metabolic network to simulate tomato leaf metabolism. The resulting model includes 3410 genes and 2143 biochemical and transport reactions distributed across five intracellular organelles including cytosol, plastid, mitochondrion, peroxisome and vacuole. The model successfully described the known metabolic behaviour of tomato leaf under heterotrophic and phototrophic conditions. The in silico investigation of the metabolic characteristics for photorespiration and other relevant metabolic processes under drought stress suggested that: (i) the flux distributions through the mevalonate (MVA) pathway under drought were distinct from that under normal conditions; and (ii) the changes in fluxes through core metabolic pathways with varying flux ratio of RubisCO carboxylase to oxygenase may contribute to the adaptive stress response of plants. In addition, we improved on previous studies of reaction essentiality analysis for leaf metabolism by including potential alternative routes for compensating reaction knockouts. Altogether, the genome‐scale model provides a sound framework for investigating tomato metabolism and gives valuable insights into the functional consequences of abiotic stresses.     

New approaches in bioprocess‐control: Consortium guidance by synthetic cell‐cell communication based on fungal pheromones

Bioconversions in industrial processes are currently dominated by single‐strain approaches. With the growing complexity of tasks to be carried out, microbial consortia become increasingly advantageous and eventually may outperform single‐strain fermentations. Consortium approaches benefit from the combined metabolic capabilities of highly specialized strains and species, and the inherent division of labor reduces the metabolic burden for each strain while increasing product yields and reaction specificities. However, consortium‐based designs still suffer from a lack of available tools to control the behavior and performance of the individual subpopulations and of the entire consortium. Here, we propose to implement novel control elements for microbial consortia based on artificial cell–cell communication via fungal mating pheromones. Coupling to the desired output is mediated by pheromone‐responsive gene expression, thereby creating pheromone‐dependent communication channels between different subpopulations of the consortia. We highlight the benefits of artificial communication to specifically target individual subpopulations of microbial consortia and to control e.g. their metabolic profile or proliferation rate in a predefined and customized manner. Due to the steadily increasing knowledge of sexual cycles of industrially relevant fungi, a growing number of strains and species can be integrated into pheromone‐controlled sensor‐actor systems, exploiting their unique metabolic properties for microbial consortia approaches.     

Comparative assessment of the Euglena gracilis var. saccharophila variant strain as a producer of the β‐1,3‐glucan paramylon under varying light conditions

Euglena gracilis Z and a “sugar loving” variant strain E. gracilis var. saccharophila were investigated as producers of paramylon, a β‐1,3‐glucan polysaccharide with potential medicinal and industrial applications. The strains were grown under diurnal or dark growth conditions on a glucose–yeast extract medium supporting high‐level paramylon production. Both strains produced the highest paramylon yields (7.4–8 g · L^?1, respectively) while grown in the dark, but the maximum yield was achieved faster by E. gracilis var. saccharophila (48 h vs. 72 h). The glucose‐to‐paramylon yield coefficient Y_par/glu = 0.46 ± 0.03 in the E. gracilis var. saccharophila cultivation, obtained in this study, is the highest reported to date. Proteomic analysis of the metabolic pathways provided molecular clues for the strain behavior observed during cultivation. For example, overexpression of enzymes in the gluconeogenesis/glycolysis pathways including fructokinase‐1 and chloroplastic fructose‐1,6‐bisphosphatase (FBP ) may have contributed to the faster rate of paramylon accumulation in E. gracilis var. saccharophila . Differentially expressed proteins in the early steps of chloroplastogenesis pathway including plastid uroporphyrinogen decarboxylases, photoreceptors, and a highly abundant (68‐fold increase) plastid transketolase may have provided the E. gracilis var. saccharophila strain an advantage in paramylon production during diurnal cultivations. In conclusion, the variant strain E. gracilis var. saccharophila seems to be well suited for producing large amounts of paramylon. This work has also resulted in the identification of molecular targets for future improvement of paramylon production in E. gracilis , including the FBP and phosophofructokinase 1, the latter being a key regulator of glycolysis.     

MapMaker and PathTracer for tracking carbon in genome‐scale metabolic models

Constraint‐based reconstruction and analysis (COBRA) modeling results can be difficult to interpret given the large numbers of reactions in genome‐scale models. While paths in metabolic networks can be found, existing methods are not easily combined with constraint‐based approaches. To address this limitation, two tools (MapMaker and PathTracer) were developed to find paths (including cycles) between metabolites, where each step transfers carbon from reactant to product. MapMaker predicts carbon transfer maps (CTMs) between metabolites using only information on molecular formulae and reaction stoichiometry, effectively determining which reactants and products share carbon atoms. MapMaker correctly assigned CTMs for over 97% of the 2,251 reactions in an Escherichia coli metabolic model (iJO1366). Using CTMs as inputs, PathTracer finds paths between two metabolites. PathTracer was applied to iJO1366 to investigate the importance of using CTMs and COBRA constraints when enumerating paths, to find active and high flux paths in flux balance analysis (FBA) solutions, to identify paths for putrescine utilization, and to elucidate a potential CO₂ fixation pathway in E. coli. These results illustrate how MapMaker and PathTracer can be used in combination with constraint‐based models to identify feasible, active, and high flux paths between metabolites.