Genetic and biocatalytic basis of formate dependent growth of Escherichia coli strains evolved in continuous culture.

The reductive glycine pathway was described as the most energetically favorable synthetic route of aerobic formate assimilation. Here we report the successful implementation of formatotrophy in Escherichia coli by means of a stepwise adaptive evolution strategy. Medium swap and turbidostat regimes of continuous culture were applied to force the channeling of carbon flux through the synthetic pathway to pyruvate establishing growth on formate and CO2 as sole carbon sources. Labeling with 13C-formate proved the assimilation of the C1 substrate via the pathway metabolites. Genetic analysis of intermediate isolates revealed a mutational path followed throughout the adaptation process. Mutations were detected affecting the copy number (gene ftfL) or the coding sequence (genes folD and lpd) of genes which specify enzymes implicated in the three steps forming glycine from formate and CO2, the central metabolite of the synthetic pathway. The mutation R191S present in methylene-tetrahydrofolate dehydrogenase/cyclohydrolase (FolD) abolishes the inhibition of cyclohydrolase activity by the substrate formyl-tetrahydrofolate. The mutation R273H in lipoamide dehydrogenase (Lpd) alters substrate affinities as well as kinetics at physiological substrate concentrations likely favoring a reactional shift towards lipoamide reduction. In addition, genetic reconstructions proved the necessity of all three mutations for formate assimilation by the adapted cells. The largely unpredictable nature of these changes demonstrates the usefulness of the evolutionary approach enabling the selection of adaptive mutations crucial for pathway engineering of biotechnological model organisms.

Activating Silent Glycolysis Bypasses in Escherichia coli.

All living organisms share similar reactions within their central metabolism to provide precursors for all essential building blocks and reducing power. To identify whether alternative metabolic routes of glycolysis can operate in E. coli, we complementarily employed in silico design, rational engineering, and adaptive laboratory evolution. First, we used a genome-scale model and identified two potential pathways within the metabolic network of this organism replacing canonical Embden-Meyerhof-Parnas (EMP) glycolysis to convert phosphosugars into organic acids. One of these glycolytic routes proceeds via methylglyoxal and the other via serine biosynthesis and degradation. Then, we implemented both pathways in E. coli strains harboring defective EMP glycolysis. Surprisingly, the pathway via methylglyoxal seemed to immediately operate in a triosephosphate isomerase deletion strain cultivated on glycerol. By contrast, in a phosphoglycerate kinase deletion strain, the overexpression of methylglyoxal synthase was necessary to restore growth of the strain. Furthermore, we engineered the "serine shunt" which converts 3-phosphoglycerate via serine biosynthesis and degradation to pyruvate, bypassing an enolase deletion. Finally, to explore which of these alternatives would emerge by natural selection, we performed an adaptive laboratory evolution study using an enolase deletion strain. Our experiments suggest that the evolved mutants use the serine shunt. Our study reveals the flexible repurposing of metabolic pathways to create new metabolite links and rewire central metabolism.

Change in Cofactor Specificity of Oxidoreductases by Adaptive Evolution of an Escherichia coli NADPH-Auxotrophic Strain.

The nicotinamide cofactor specificity of enzymes plays a key role in regulating metabolic processes and attaining cellular homeostasis. Multiple studies have used enzyme engineering tools or a directed evolution approach to switch the cofactor preference of specific oxidoreductases. However, whole-cell adaptation toward the emergence of novel cofactor regeneration routes has not been previously explored. To address this challenge, we used an Escherichia coli NADPH-auxotrophic strain. We continuously cultivated this strain under selective conditions. After 500 to 1,100 generations of adaptive evolution using different carbon sources, we isolated several strains capable of growing without an external NADPH source. Most isolated strains were found to harbor a mutated NAD+-dependent malic enzyme (MaeA). A single mutation in MaeA was found to switch cofactor specificity while lowering enzyme activity. Most mutated MaeA variants also harbored a second mutation that restored the catalytic efficiency of the enzyme. Remarkably, the best MaeA variants identified this way displayed overall superior kinetics relative to the wild-type variant with NAD+. In other evolved strains, the dihydrolipoamide dehydrogenase (Lpd) was mutated to accept NADP+, thus enabling the pyruvate dehydrogenase and 2-ketoglutarate dehydrogenase complexes to regenerate NADPH. Interestingly, no other central metabolism oxidoreductase seems to evolve toward reducing NADP+, which we attribute to several biochemical constraints, including unfavorable thermodynamics. This study demonstrates the potential and biochemical limits of evolving oxidoreductases within the cellular context toward changing cofactor specificity, further showing that long-term adaptive evolution can optimize enzyme activity beyond what is achievable via rational design or directed evolution using small libraries. IMPORTANCE In the cell, NAD(H) and NADP(H) cofactors have different functions. The former mainly accepts electrons from catabolic reactions and carries them to respiration, while the latter provides reducing power for anabolism. Correspondingly, the ratio of the reduced to the oxidized form differs for NAD+ (low) and NADP+ (high), reflecting their distinct roles. We challenged the flexibility of E. coli's central metabolism in multiple adaptive evolution experiments using an NADPH-auxotrophic strain. We found several mutations in two enzymes, changing the cofactor preference of malic enzyme and dihydrolipoamide dehydrogenase. Upon deletion of their corresponding genes we performed additional evolution experiments which did not lead to the emergence of any additional mutants. We attribute this restricted number of mutational targets to intrinsic thermodynamic barriers; the high ratio of NADPH to NADP+ limits metabolic redox reactions that can regenerate NADPH, mainly by mass action constraints.

Continuous Culture Adaptation of Methylobacterium extorquens AM1 and TK 0001 to Very High Methanol Concentrations.

The bio-economy relies on microbial strains optimized for efficient large scale production of chemicals and fuels from inexpensive and renewable feedstocks under industrial conditions. The reduced one carbon compound methanol, whose production does not involve carbohydrates needed for the feed and food sector, can be used as sole carbon and energy source by methylotrophic bacteria like Methylobacterium extorquens AM1. This strain has already been engineered to produce various commodity and high value chemicals from methanol. The toxic effect of methanol limits its concentration as feedstock to 1% v/v. We obtained M. extorquens chassis strains tolerant to high methanol via adaptive directed evolution using the GM3 technology of automated continuous culture. Turbidostat and conditional medium swap regimes were employed for the parallel evolution of the recently characterized strain TK 0001 and the reference strain AM1 and enabled the isolation of derivatives of both strains capable of stable growth with 10% methanol. The isolates produced more biomass at 1% methanol than the ancestor strains. Genome sequencing identified the gene metY coding for an O-acetyl-L-homoserine sulfhydrylase as common target of mutation. We showed that the wildtype enzyme uses methanol as substrate at elevated concentrations. This side reaction produces methoxine, a toxic homolog of methionine incorporated in polypeptides during translation. All mutated metY alleles isolated from the evolved populations coded for inactive enzymes, designating O-acetyl-L-homoserine sulfhydrylase as a major vector of methanol toxicity. A whole cell transcriptomic analysis revealed that genes coding for chaperones and proteases were upregulated in the evolved cells as compared with the wildtype, suggesting that the cells had to cope with aberrant proteins formed during the adaptation to increasing methanol exposure. In addition, the expression of ribosomal proteins and enzymes related to energy production from methanol like formate dehydrogenases and ATP synthases was boosted in the evolved cells upon a short-term methanol stress. D-lactate production from methanol by adapted cells overexpressing the native D-lactate dehydrogenase was quantified. A significant higher lactate yield was obtained compared with control cells, indicating an enhanced capacity of the cells resistant to high methanol to assimilate this one carbon feedstock more efficiently.

Evolution of a Biomass-Fermenting Bacterium To Resist Lignin Phenolics.

Increasing the resistance of plant-fermenting bacteria to lignocellulosic inhibitors is useful to understand microbial adaptation and to develop candidate strains for consolidated bioprocessing. Here, we study and improve inhibitor resistance in Clostridium phytofermentans (also called Lachnoclostridium phytofermentans), a model anaerobe that ferments lignocellulosic biomass. We survey the resistance of this bacterium to a panel of biomass inhibitors and then evolve strains that grow in increasing concentrations of the lignin phenolic, ferulic acid, by automated, long-term growth selection in an anaerobic GM3 automat. Ultimately, strains resist multiple inhibitors and grow robustly at the solubility limit of ferulate while retaining the ability to ferment cellulose. We analyze genome-wide transcription patterns during ferulate stress and genomic variants that arose along the ferulate growth selection, revealing how cells adapt to inhibitors through changes in gene dosage and regulation, membrane fatty acid structure, and the surface layer. Collectively, this study demonstrates an automated framework for in vivo directed evolution of anaerobes and gives insight into the genetic mechanisms by which bacteria survive exposure to chemical inhibitors.IMPORTANCE Fermentation of plant biomass is a key part of carbon cycling in diverse ecosystems. Further, industrial biomass fermentation may provide a renewable alternative to fossil fuels. Plants are primarily composed of lignocellulose, a matrix of polysaccharides and polyphenolic lignin. Thus, when microorganisms degrade lignocellulose to access sugars, they also release phenolic and acidic inhibitors. Here, we study how the plant-fermenting bacterium Clostridium phytofermentans resists plant inhibitors using the lignin phenolic, ferulic acid. We examine how the cell responds to abrupt ferulate stress by measuring changes in gene expression. We evolve increasingly resistant strains by automated, long-term cultivation at progressively higher ferulate concentrations and sequence their genomes to identify mutations associated with acquired ferulate resistance. Our study develops an inhibitor-resistant bacterium that ferments cellulose and provides insights into genomic evolution to resist chemical inhibitors.

Contrasted microcolinearity and gene evolution within a homoeologous region of wheat and barley species.

We study here the evolution of genes located in the same physical locus using the recently sequenced Ha locus in seven wheat genomes in diploid, tetraploid, and hexaploid species and compared them with barley and rice orthologous regions. We investigated both the conservation of microcolinearity and the molecular evolution of genes, including coding and noncoding sequences. Microcolinearity is restricted to two groups of genes (Unknown gene-2, VAMP, BGGP, Gsp-1, and Unknown gene-8 surrounded by several copies of ATPase), almost conserved in rice and barley, but in a different relative position. Highly conserved genes between wheat and rice run along with genes harboring different copy numbers and highly variable sequences between close wheat genomes. The coding sequence evolution appeared to be submitted to heterogeneous selective pressure and intronic sequences analysis revealed that the molecular clock hypothesis is violated in most cases.

Molecular basis of evolutionary events that shaped the hardness locus in diploid and polyploid wheat species (Triticum and Aegilops).

The Hardness (Ha) locus controls grain hardness in hexaploid wheat (Triticum aestivum) and its relatives (Triticum and Aegilops species) and represents a classical example of a trait whose variation arose from gene loss after polyploidization. In this study, we investigated the molecular basis of the evolutionary events observed at this locus by comparing corresponding sequences of diploid, tertraploid, and hexaploid wheat species (Triticum and Aegilops). Genomic rearrangements, such as transposable element insertions, genomic deletions, duplications, and inversions, were shown to constitute the major differences when the same genomes (i.e., the A, B, or D genomes) were compared between species of different ploidy levels. The comparative analysis allowed us to determine the extent and sequences of the rearranged regions as well as rearrangement breakpoints and sequence motifs at their boundaries, which suggest rearrangement by illegitimate recombination. Among these genomic rearrangements, the previously reported Pina and Pinb genes loss from the Ha locus of polyploid wheat species was caused by a large genomic deletion that probably occurred independently in the A and B genomes. Moreover, the Ha locus in the D genome of hexaploid wheat (T. aestivum) is 29 kb smaller than in the D genome of its diploid progenitor Ae. tauschii, principally because of transposable element insertions and two large deletions caused by illegitimate recombination. Our data suggest that illegitimate DNA recombination, leading to various genomic rearrangements, constitutes one of the major evolutionary mechanisms in wheat species.

The DNA sequence and analysis of human chromosome 14.

Chromosome 14 is one of five acrocentric chromosomes in the human genome. These chromosomes are characterized by a heterochromatic short arm that contains essentially ribosomal RNA genes, and a euchromatic long arm in which most, if not all, of the protein-coding genes are located. The finished sequence of human chromosome 14 comprises 87,410,661 base pairs, representing 100% of its euchromatic portion, in a single continuous segment covering the entire long arm with no gaps. Two loci of crucial importance for the immune system, as well as more than 60 disease genes, have been localized so far on chromosome 14. We identified 1,050 genes and gene fragments, and 393 pseudogenes. On the basis of comparisons with other vertebrate genomes, we estimate that more than 96% of the chromosome 14 genes have been annotated. From an analysis of the CpG island occurrences, we estimate that 70% of these annotated genes are complete at their 5' end.