Genetic changes during a laboratory adaptive evolution process that allowed fast growth in glucose to an Escherichia coli strain lacking the major glucose transport system

BackgroundEscherichia coli strains lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system (PTS), which is the major bacterial component involved in glucose transport and its phosphorylation, accumulate high amounts of phosphoenolpyruvate that can be diverted to the synthesis of commercially relevant products. However, these strains grow slowly in glucose as sole carbon source due to its inefficient transport and metabolism. Strain PB12, with 400% increased growth rate, was isolated after a 120 hours adaptive laboratory evolution process for the selection of faster growing derivatives in glucose. Analysis of the genetic changes that occurred in the PB12 strain that lacks PTS will allow a better understanding of the basis of its growth adaptation and, therefore, in the design of improved metabolic engineering strategies for enhancing carbon diversion into the aromatic pathways.ResultsWhole genome analyses using two different sequencing methodologies: the Roche NimbleGen Inc. comparative genome sequencing technique, and high throughput sequencing with Illumina Inc. GAIIx, allowed the identification of the genetic changes that occurred in the PB12 strain. Both methods detected 23 non-synonymous and 22 synonymous point mutations. Several non-synonymous mutations mapped in regulatory genes (arcB, barA, rpoD, rna) and in other putative regulatory loci (yjjU, rssA and ypdA). In addition, a chromosomal deletion of 10,328 bp was detected that removed 12 genes, among them, the rppH, mutH and galR genes. Characterization of some of these mutated and deleted genes with their functions and possible functions, are presented.ConclusionsThe deletion of the contiguous rppH, mutH and galR genes that occurred simultaneously, is apparently the main reason for the faster growth of the evolved PB12 strain. In support of this interpretation is the fact that inactivation of the rppH gene in the parental PB11 strain substantially increased its growth rate, very likely by increasing glycolytic mRNA genes stability. Furthermore, galR inactivation allowed glucose transport by GalP into the cell. The deletion of mutH in an already stressed strain that lacks PTS is apparently responsible for the very high mutation rate observed.

[1]  Thomas D. Schmittgen,et al.  Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. , 2001, Methods.

[2]  B. Palsson,et al.  Scalable method to determine mutations that occur during adaptive evolution of Escherichia coli , 2003, Biotechnology Letters.

[3]  Sangeeta Banerji,et al.  Patatin-like proteins: a new family of lipolytic enzymes present in bacteria? , 2004, Microbiology.

[4]  F. Bolivar,et al.  New Insights into the Role of Sigma Factor RpoS as Revealed in Escherichia coli Strains Lacking the Phosphoenolpyruvate:Carbohydrate Phosphotransferase System , 2007, Journal of Molecular Microbiology and Biotechnology.

[5]  L. Chao,et al.  Understanding the Evolutionary Fate of Finite Populations: The Dynamics of Mutational Effects , 2007, PLoS biology.

[6]  D. Georgellis,et al.  Identification of a quinone-sensitive redox switch in the ArcB sensor kinase. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[7]  A. Danchin,et al.  Organised Genome Dynamics in the Escherichia coli Species Results in Highly Diverse Adaptive Paths , 2009, PLoS genetics.

[8]  M. Schofield,et al.  DNA mismatch repair: molecular mechanisms and biological function. , 2003, Annual review of microbiology.

[9]  Xueqiao Liu,et al.  Probing the ArcA-P Modulon of Escherichia coli by Whole Genome Transcriptional Analysis and Sequence Recognition Profiling* , 2004, Journal of Biological Chemistry.

[10]  T. Richmond,et al.  Mutation discovery in bacterial genomes: metronidazole resistance in Helicobacter pylori , 2005, Nature Methods.

[11]  I. Matic,et al.  Evolutionary significance of stress-induced mutagenesis in bacteria. , 2004, Trends in microbiology.

[12]  R. Kaplan,et al.  Accumulation of nucleotides by starved Escherichia coli cells as a probe for the involvement of ribonucleases in ribonucleic acid degradation , 1977, Journal of bacteriology.

[13]  S. Darst,et al.  Structure of the Bacterial RNA Polymerase Promoter Specificity σ Subunit , 2002 .

[14]  B. Wanner,et al.  One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[15]  Wing-Kin Sung,et al.  PE-Assembler: de novo assembler using short paired-end reads , 2011, Bioinform..

[16]  D. Georgellis,et al.  Identification of UvrY as the Cognate Response Regulator for the BarA Sensor Kinase in Escherichia coli * , 2001, The Journal of Biological Chemistry.

[17]  S A Bustin,et al.  Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. , 2002, Journal of molecular endocrinology.

[18]  Marie Nguyen,et al.  A practical approach to RT-qPCR-Publishing data that conform to the MIQE guidelines. , 2010, Methods.

[19]  V. Beneš,et al.  The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. , 2009, Clinical chemistry.

[20]  Francisco Bolívar,et al.  Nutrient-Scavenging Stress Response in an Escherichia coli Strain Lacking the Phosphoenolpyruvate:Carbohydrate Phosphotransferase System, as Explored by Gene Expression Profile Analysis , 2006, Journal of Molecular Microbiology and Biotechnology.

[21]  F. Bolivar,et al.  Transcription Analysis of Central Metabolism Genes in Escherichia coli. Possible Roles of σ38 in Their Expression, as a Response to Carbon Limitation , 2009, PloS one.

[22]  Tom M. Conrad,et al.  Whole-genome resequencing of Escherichia coli K-12 MG1655 undergoing short-term laboratory evolution in lactate minimal media reveals flexible selection of adaptive mutations , 2009, Genome Biology.

[23]  Monica Riley,et al.  A functional update of the Escherichia coli K-12 genome , 2001, Genome Biology.

[24]  T. Silhavy,et al.  RpoS-Dependent Transcriptional Control ofsprE: Regulatory Feedback Loop , 2001, Journal of bacteriology.

[25]  Francisco Bolívar,et al.  Metabolic engineering for the production of shikimic acid in an evolved Escherichia coli strain lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system , 2010, Microbial cell factories.

[26]  H. Čelešnik,et al.  The bacterial enzyme RppH triggers messenger RNA degradation by 5′ pyrophosphate removal , 2008, Nature.

[27]  L. Wyns,et al.  Nonspecific base recognition mediated by water bridges and hydrophobic stacking in ribonuclease I from Escherichia coli , 2008, Protein science : a publication of the Protein Society.

[28]  D. Apirion,et al.  Decay of ribosomal ribonucleic acid in Escherichia coli cells starved for various nutrients. , 1975, The Journal of biological chemistry.

[29]  F. Bolivar,et al.  Metabolic engineering and protein directed evolution increase the yield of L‐phenylalanine synthesized from glucose in Escherichia coli , 2004, Biotechnology and bioengineering.

[30]  Jeffrey H. Miller,et al.  In vitro and in vivo studies of MutS, MutL and MutH mutants: correlation of mismatch repair and DNA recombination. , 2003, DNA repair.

[31]  Valeria Souza,et al.  Stress-Induced Mutagenesis in Bacteria , 2003, Science.

[32]  M. Berriman,et al.  Improving draft assemblies by iterative mapping and assembly of short reads to eliminate gaps , 2010, Genome Biology.

[33]  Alfredo Martínez,et al.  Growth Recovery on Glucose under Aerobic Conditions of an Escherichia coli Strain Carrying a Phosphoenolpyruvate:Carbohydrate Phosphotransferase System Deletion by Inactivating arcA and Overexpressing the Genes Coding for Glucokinase and Galactose Permease , 2007, Journal of Molecular Microbiology and Biotechnology.

[34]  S. Adhya,et al.  Functional characterization of roles of GalR and GalS as regulators of the gal regulon , 1997, Journal of bacteriology.

[35]  Hans G. Schlegel,et al.  Biology of the prokaryotes , 1999 .

[36]  Alfredo Martínez,et al.  Role of Pyruvate Oxidase in Escherichia coli Strains Lacking the Phosphoenolpyruvate:Carbohydrate Phosphotransferase System , 2005, Journal of Molecular Microbiology and Biotechnology.

[37]  F. Neidhardt,et al.  Phosphoenolpyruvate:carbohydrate phosphotransferase systems , 1996 .

[38]  B. Palsson,et al.  An evaluation of Comparative Genome Sequencing (CGS) by comparing two previously-sequenced bacterial genomes , 2007, BMC Genomics.

[39]  S. Kanaya,et al.  Large-scale identification of protein-protein interaction of Escherichia coli K-12. , 2006, Genome research.

[40]  Kris Popendorf,et al.  Murasaki: A Fast, Parallelizable Algorithm to Find Anchors from Multiple Genomes , 2010, PloS one.

[41]  N. W. Davis,et al.  The complete genome sequence of Escherichia coli K-12. , 1997, Science.

[42]  Aaron E. Darling,et al.  Reordering contigs of draft genomes using the Mauve Aligner , 2009, Bioinform..

[43]  Francisco Bolívar,et al.  Adaptation for fast growth on glucose by differential expression of central carbon metabolism and gal regulon genes in an Escherichia coli strain lacking the phosphoenolpyruvate:carbohydrate phosphotransferase system. , 2005, Metabolic engineering.

[44]  F. Bolivar,et al.  Pathway engineering for the production of aromatic compounds in Escherichia coli , 1996, Nature Biotechnology.

[45]  Chen Yang,et al.  Analysis of Gene Expression in Escherichia coli in Response to Changes of Growth-Limiting Nutrient in Chemostat Cultures , 2004, Applied and Environmental Microbiology.

[46]  E. Lin,et al.  Mutational analysis of signal transduction by ArcB, a membrane sensor protein responsible for anaerobic repression of operons involved in the central aerobic pathways in Escherichia coli , 1992, Journal of bacteriology.

[47]  J. W. Campbell,et al.  Experimental Determination and System Level Analysis of Essential Genes in Escherichia coli MG1655 , 2003, Journal of bacteriology.

[48]  S. Normark,et al.  The Escherichia coli BarA-UvrY Two-Component System Is Needed for Efficient Switching between Glycolytic and Gluconeogenic Carbon Sources , 2003, Journal of bacteriology.

[49]  F. Bolivar,et al.  Determination of 3-deoxy-D-arabino-heptulosonate 7-phosphate productivity and yield from glucose in Escherichia coli devoid of the glucose phosphotransferase transport system. , 2001, Biotechnology and bioengineering.

[50]  A. D. de Graaf,et al.  Analysis of carbon metabolism in Escherichia coli strains with an inactive phosphotransferase system by (13)C labeling and NMR spectroscopy. , 2002, Metabolic Engineering.

[51]  Howard Ochman,et al.  Neutral mutations and neutral substitutions in bacterial genomes. , 2003, Molecular biology and evolution.

[52]  B. Bachmann,et al.  Derivations and genotypes of some mutant derivatives of Escherichia coli K12 , 1987 .

[53]  F. Blattner,et al.  Mauve: multiple alignment of conserved genomic sequence with rearrangements. , 2004, Genome research.

[54]  Anu Raghunathan,et al.  Comparative genome sequencing of Escherichia coli allows observation of bacterial evolution on a laboratory timescale , 2006, Nature Genetics.

[55]  W. J. Kent,et al.  BLAT--the BLAST-like alignment tool. , 2002, Genome research.

[56]  Bernhard Ø. Palsson,et al.  Genetic Basis of Growth Adaptation of Escherichia coli after Deletion of pgi, a Major Metabolic Gene , 2010, PLoS genetics.

[57]  Guillermo Gosset,et al.  Coutilization of glucose and glycerol enhances the production of aromatic compounds in an Escherichia coli strain lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system , 2008, Microbial cell factories.

[58]  Ken Chen,et al.  VarScan: variant detection in massively parallel sequencing of individual and pooled samples , 2009, Bioinform..

[59]  Robert Gentleman,et al.  ShortRead: a bioconductor package for input, quality assessment and exploration of high-throughput sequence data , 2009, Bioinform..

[60]  Cole Trapnell,et al.  Ultrafast and memory-efficient alignment of short DNA sequences to the human genome , 2009, Genome Biology.