Comparative multi-omics systems analysis of Escherichia coli strains B and K-12
© Yoon et al.; licensee BioMed Central Ltd. 2012
Received: 9 March 2012
Accepted: 25 May 2012
Published: 25 May 2012
Elucidation of a genotype-phenotype relationship is critical to understand an organism at the whole-system level. Here, we demonstrate that comparative analyses of multi-omics data combined with a computational modeling approach provide a framework for elucidating the phenotypic characteristics of organisms whose genomes are sequenced.
We present a comprehensive analysis of genome-wide measurements incorporating multifaceted holistic data - genome, transcriptome, proteome, and phenome - to determine the differences between Escherichia coli B and K-12 strains. A genome-scale metabolic network of E. coli B was reconstructed and used to identify genetic bases of the phenotypes unique to B compared with K-12 through in silico complementation testing. This systems analysis revealed that E. coli B is well-suited for production of recombinant proteins due to a greater capacity for amino acid biosynthesis, fewer proteases, and lack of flagella. Furthermore, E. coli B has an additional type II secretion system and a different cell wall and outer membrane composition predicted to be more favorable for protein secretion. In contrast, E. coli K-12 showed a higher expression of heat shock genes and was less susceptible to certain stress conditions.
This integrative systems approach provides a high-resolution system-wide view and insights into why two closely related strains of E. coli, B and K-12, manifest distinct phenotypes. Therefore, systematic understanding of cellular physiology and metabolism of the strains is essential not only to determine culture conditions but also to design recombinant hosts.
Escherichia coli is one of the most intensively studied organisms and has been widely employed in scientific studies and industrial applications. The most widely used E. coli have been those derived from strains B and K-12, the result of pioneering work using K-12 for genetic and biochemical studies and B for studying virulent bacteriophages, restriction systems, mutagenic assays, and bacterial evolution [1, 2]. The first whole-genome sequence of strain K-12, MG1655, was determined  and compared in detail with another K-12 strain, W3110 . We have determined genome sequences of B strains [5, 6] - REL606, which has been applied to the study of long-term experimental evolution, and BL21(DE3), which has been used as a cell-factory for overproducing recombinant proteins, biofuels, and a variety of bioproducts on an industrial scale.
Previously, we showed that comparison of the genome sequences of strains B and K-12 could provide plausible explanations for some long-known differences between them . However, genome sequence alone provides limited information about the genotype-phenotype relationship . Starting with a genome sequence, comprehensive analyses of transcriptomes, proteomes, and phenomes can complement each other and be integrated to provide deeper insights into biological systems [9–12]. Furthermore, the multidimensional omics data can be integrated to reconstruct genome-wide computational models that can generate testable hypotheses concerning cellular function [13, 14].
Here, we systematically combined the results of comparative analyses of the genomes, transcriptomes, proteomes, and phenomes of E. coli B REL606 and K-12 MG1655 to decipher the whole-organism characteristics that differentiate these two strains. An in silico genome-scale metabolic model of E. coli B REL606 was reconstructed and used to determine the genetic basis of the phenotypic differences.
Bacterial cell growth
Two E. coli B strains , REL606 and BL21(DE3), and two K-12 strains , MG1655 and W3110, were cultured in Luria-Bertani (LB) medium or R/2 medium (Figure S1 in Additional file 1). Both groups of strains grew similarly when cultured in complex medium; however, the B strains grew faster than the K12 strains in minimal medium. The difference between the strain types in the accumulation of by-products such as acetic, formic, and lactic acids was negligible when cultured in complex medium compared with growth in minimal media.
Many genes showed highly distinct expression levels in REL606 and MG1655 irrespective of growth conditions (upper right for REL606 and lower left for MG1655 in Figure 2; also see Table S2 in Additional file 1). Highly expressed genes in REL606 included those encoding enzymes for biosynthesis of L-arginine (argAGDECBHI) and branched-chain amino acids (ilvGMEDA), a subunit of the L-arginine transporter (artJ), cytochrome b562 (cybC), subunits of the histidine ABC transporter (hisPJ), cytotoxins (hokED), outer membrane porin (ompF), L-arginine decarboxylase (speA), and a cell division inhibitor (sulA). In MG1655, highly expressed genes included those involved in chemotaxis (cheZYRWA, tap, trg, tsr), Lon protease (lon), C4-dicarboxylate-sensing histidine kinase (dcuS), chaperones (clpB, dnaK, groES, htpG, ibpA), the major subunit of type 1 fimbriae (fimA), a regulator of flagellar biosynthesis (flhC), glycerol-3-phosphate-dehydrogenase (glpABCD), glycerophosphoryl diester phosphodiesterase (glpQ), glycerol-3-phosphate transporter (glpT), hydrogenase 2 (hybCBO), outer membrane porins (nmpC, ompA, ompC), and galactitol transport and metabolism (gatYZC).
Proteome analyses of the total intracellular, outer membrane, and extracellular proteomes of the four E. coli strains cultured in LB medium (Figure S2 in Additional file 1) revealed 18 protein spots of B strains and 42 spots of K-12 strains that showed a more than two-fold difference in intensity and the corresponding proteins were further characterized (Table S3 in Additional file 1). Intracellular proteins that were more abundant in the B strains were enzymes required for the biosynthesis of some amino acids (AspC, ArgCDI, SerC) and maltose metabolism (MalPQ). Proteins that were more abundant in K-12 strains included enzymes for the transport and metabolism of galactitol (GatYZAB), enzymes for the degradation of amino acids (TdcE, TnaA), and stress response proteins (ClpP, CspE, YfiD). The proteome profiles of the two B strains were very similar, as were those of the two K-12 strains (Figure S2A in Additional file 1). Proteins involved in chemotaxis (CheY) and a degradation pathway for amino acids (TdcF) were synthesized at a significant level only by MG1655.
SDS-PAGE analysis of outer membrane proteins (OMPs) revealed that the B strains expressed large amounts of OmpF but not OmpC, whereas K-12 strains expressed both OmpF and OmpC (Figure S2B in Additional file 1). The expression level of OmpA was higher in K-12 strains than in B strains.
Phenotype microarray (PM) tests on REL606 and MG1655 revealed dramatic phenotypic differences between the two strains (Figure S3 in Additional file 1). The B strain was more susceptible than the K-12 strain to a variety of stressful conditions caused by osmolarity, pH, or exposure to inhibitory compounds such as salicylate and β-lactam antibiotics. MG1655 did not grow on valine dipeptides.
Reconstruction of the metabolic network
We reconstructed a genome-scale metabolic network model of REL606 by incorporating genetic differences between REL606 and MG1655 (Table S4 in Additional file 1). Starting with the metabolic network model for MG1655 (iAF1260) , 29 REL606-specific reactions and 11 REL606-specific compounds were added, 12 REL606-specific regulations were introduced, and 43 MG1655-specific reactions were excluded. The resulting metabolic model of REL606 included 1,369 metabolic reactions and 1,051 metabolites. To validate the model, flux balance analysis (FBA) was used to predict cell growth on various carbon sources, and the results were compared with the results from PM tests (Table S5 in Additional file 1). The model qualitatively agreed with 83.2% of experimental cell growth results on PM plates testing carbon source utilization.
Lower acetate accumulation of the B strain might be attributed to different regulation of the glyoxylate pathway
Acetate accumulation is one of the major problems encountered during high cell density cultivation of E. coli because it inhibits cell growth and production of foreign proteins . Low acetate accumulation by B strains is commonly attributed to the active glyoxylate shunt encoded by the acetate operon (aceBAK); however, in the exponential growth stage, up-regulation of genes involved in the glyoxylate shunt, tricarboxylic acid (TCA) cycle cycle, and gluconeogenesis was observed in K-12 strains (Figure 4). These expression patterns are similar to results from previous microarray analysis comparing transcriptomes of K-12 JM109 and BL21 during low glucose-fed batch fermentation . The need for high activity of the TCA cycle in K-12 strains cultured in low glucose media is not immediately obvious, but enhanced catabolic activity may be required for the energy-consuming motility and stress responses (see below).
Genes involved in the central metabolism and the acetate simulation/dissimilation pathways (glycolysis, TCA cycle, pentose phosphate pathway, glyoxylate pathway, gluconeogenesis, and acetate production) were virtually the same at the nucleotide level in B and K-12 strains. This was also the case for the transcriptional regulator gene fruR, which was suspected to have altered expression or regulatory strength in BL21, leading to derepression of acetate metabolism even under high glucose culture conditions . Interestingly, arpA, which has an unknown function but has been proposed to play a role in regulation of the glyoxylate pathway , is disrupted in B strains.
Many of the phenotypic differences between B and K-12 strains in utilization of carbon sources, as shown in PM tests (Figure S3 in Additional file 1), could be correlated with genotypic differences (Table S5 in Additional file 1). In particular, genes for galactitol transport and metabolism were highly expressed at both the mRNA level (gatYZC) and the protein level (GatYZAB) in MG1655, possibly because gatR, a repressor of the gat genes, is truncated in MG1655 but intact in REL606.
Enhanced amino acid biosynthesis and reduced degradation can make the B strain favorable for efficient production of recombinant proteins
Lack of flagellar biosynthetic genes and lower expression of motility genes are important properties of strain B as a protein factory
The cell envelope is the principal stress-bearing and shape maintenance element in E. coli, and its integrity is critical for cell viability. B strains have been widely used for mutagenic assays and toxicological studies because they show higher membrane permeability than K-12 strains . Structural studies on the LPS core oligosaccharide have revealed that K-12 strains are devoid of the O antigen while B strains lack the O antigen and the distal part of the core . Sequence comparison revealed that the outer membrane structure of B strains is quite different from that of K-12 strains . In both groups, IS elements are inserted at the gene clusters for O antigen biosynthesis: at wbbL for K-12 strains and between manC and wbbD for B strains. In the genome of B strains, the core part of LPS was further disrupted by the insertion of IS1 at the waaT gene encoding UDP-galatose:(glucosyl) LPS α1,2-galactosyltransferase. The B strain showed high expression of glmS, encoding a protein involved in both cell wall and LPS biosynthesis , at both the mRNA and protein level (Figure 4).
Unexpectedly, flagellar biosynthesis genes (fliY-R) were absent from the B genome, rendering B strains non-motile . This agrees with microarray data that indicated that expression levels of chemotaxis genes (cheZYRWA, tap, trg, and tsr) and of flhC, encoding the master regulator of flagellar biosynthesis, are much lower in B than in K-12 (Figure 4). This finding is also concordant with the proteomic data showing CheY, FimA, FlgF, FlgL, FliC, and FliD only in K-12 strains. Lack of flagellar biosynthetic genes and low expression of motility-related genes might be important properties when B strains are used as protein factories, because flagellar biosynthesis is energy-intensive and is not necessary under industrial conditions of constant agitation and generous nutrient supply .
Different cellular responses to environmental changes were observed for B and K-12 strains
The susceptibility of the E. coli B strain to a variety of stressful conditions and antibiotics revealed by PM tests (Figure S3 in Additional file 1) can be explained by several observations (Figure 5). First, differences in the composition of the LPS core and expression of outer membrane proteins may influence the permeability and integrity of the cell envelope. B strains produce more OmpF porin than K-12 strains because the B genome lacks micF, which post-transcriptionally prevents production of OmpF . This is further supported by the transcriptome data showing high levels of ompF expression in the B strain and high expression of ompC and ompA in the K-12 strain (Figure 4). These observations were also consistent with results of proteome analysis of the outer membrane fractions (Figure S2B in Additional file 1). Noxious agents such as antibiotics and bile acids diffuse more easily through OmpF than OmpC because the former produces a channel with a larger pore size . Second, synthesis of capsular colanic acid and biofilm formation appear to be repressed in B strains as a result of the absence of the Rcs phosphorelay system (rcsADB), the major mechanism regulating these processes . Proteome data also showed that FimA and Flu involved in biofilm formation are secreted only by K-12 strains. Third, genome sequence comparisons revealed at least five toxin-antitoxin (TA) modules in K-12, three relBE-class loci (relBE, dinJ-yafQ and yefM-yoeB) and two mazEF-homologous loci (chpRA and chpSB), whereas B strains possess only two TA loci (dinJ-yafQ, chpRA). TA modules play an important role in bacterial stress physiology and might contribute to multidrug resistance . Finally, the B strain showed lower stress responses. Heat shock genes encoding molecular chaperones (dnaK, groES, htpG, ibpA, and clpB) and ATP-dependent proteases (clpA, clpP, and hslV) were expressed at lower levels in the B strain (Figure 4). Furthermore, B strains are naturally deficient in the major protease Lon as a result of an IS186 insertion in the promoter region . Expression levels of DnaK and ClpB, which correlate with cellular stress levels, are lower in BL21 strains than in K-12 strains (JM105, HB101, and TOP10) . K-12 also showed activation of genes for cold shock protein (cspE), carbon starvation protein (cstA), osmotically inducible peroxidase (osmC), phage shock protein (pspE), and universal stress protein (uspA) (Figure 4), which may make K-12 strains more resistant than B strains to various stress conditions.
The rpoE gene, encoding heat-inducible sigma factor σE, and some genes in the operon (rseABC) regulated by this sigma factor, were highly expressed in REL606 (Figure 4). σE induces an extracytoplasmic stress response to misfolded outer membrane proteins and damage to the cell envelope . Variations in expression of rpoE regulon genes might result from differences in outer membrane composition between the two strain types.
Expression of the additional T2S genes of the B strain might be responsible for its larger amounts of secreted proteins
Interestingly, the B genome has an additional gene cluster for T2S that consists of the homologs gspOCDEFGHIJKLM. From a phylogenetic analysis of the concatenated sequences of GspDEFG, clear distinctions could be made between homologs of T2S that were common to B and K-12 and the additional T2S genes present in the B strain (Figure S4 in Additional file 1). It is generally assumed that laboratory strains of E. coli are unable to secrete many proteins because the operon encoding the T2S pathway (gsp) is silenced by the H-NS protein . The expression levels of gspB-O genes, except for gspE, were relatively low in both E. coli B and K-12 strains throughout the culture period. Expression of the additional T2S genes present in the B strain was detected by transcriptome analysis. Regarding secretion capacity, both REL606 and BL21(DE3) secreted larger amounts of protein than the K-12 strains did (Figure 3). All E. coli strains tested secreted greater amounts of protein when grown in minimal media, and the release of proteins was optimal in the stationary phase. Thus, the existence of the second T2S secretion system and enhanced capability for protein release suggest that B strains might be better suited for extracellular production of recombinant proteins.
In silico complementation test revealed genetic bases of the observed phenotypes
To identify the reactions responsible for the phenotypes unique to REL606 compared with MG1655, in silico complementation experiments were performed by reverting each of the REL606-specific modifications back to MG1655-specific reactions. Each combination of additions, deletions, and regulations unique to REL606 was reverted and growth under each condition where REL606 and MG1655 show differential substrate utilization phenotypes was predicted using the metabolic network models. Combinations of modifications whose reversion resulted in loss of REL606-specific phenotypes were regarded as reactions responsible for phenotypes unique to REL606.
Different phenotypes of E. coli B REL606 and E. coli K-12 MG1655 and modifications of the corresponding reactions
p-Hydroxyphenyl acetic acid
4H3M, 34DH23OR, 5C2HMSO, 5C2HMDI, 5O3E125TC, HPAI, HPAP, HPAtex, EX_4hpa(e)
ARABDI, RBK_D1, DARBtex, EX_arab-D(e)
Despite advances in various high-throughput technologies, extraction of discernable biological meaning from multiple omics data sets remains a challenging task [9, 10, 13, 33]. It has been emphasized that integration of multiple omics results is important and work often focuses on its usage in constructing and validating 'in silico models'. But how the different and nonlinearly correlated omics data can be combined in the real world is hardly commented on. In this study, we identified many important differences in cellular metabolism and physiology between E. coli B and K-12 strains through a combined and integrated analysis of genomic, transcriptomic, proteomic, and phenomic data (Figure 5). Differences in utilization of a variety of nutrient substrates were identified and in many cases could be correlated with genotypic and functional differences. Moreover, transcriptome and proteome data cross-validated information on differentially expressed genes and proteins, demonstrating the utility of integrating omics data collected at various levels from the genome to phenome. Genome-scale in silico models can generate and validate hypotheses concerning cellular processes.
Derivatives of E. coli B have been the major workhorse for production of recombinant foreign proteins and various biomaterials, including biofuels in the labs and in industry [1, 2]. Until now, the bioprocess optimization using B has been carried out in a rather trial-and-error manner. This inefficient approach has been unavoidable because of insufficient information about the metabolism and physiology of B. Therefore, the combined omics approach applied here to the cellular metabolism and physiology of E. coli B and K-12 strains should be pivotal in better understanding the underlying biological networks, and invaluable for designing strains with customized genomes and establishing rational fermentation strategies. Moreover, this omics-based systems approach will provide a framework for elucidating the phenetic characteristics of other organisms whose genomes are sequenced.
Materials and methods
Strains and culture conditions
E. coli B REL606 (F- mal + l s T6 r str r rm111 ara - ) was provided by Richard E Lenski, Michigan State University, and E. coli BL21(DE3) was obtained from F William Studier, Brookhaven National Laboratory. E. coli K-12 MG1655 was purchased from American Type Culture Collection (ATCC) and E. coli W3110 from the Korean Collection for Type Cultures (KCTC). Cells were grown in 250-ml flasks containing 100 ml LB medium or minimal R/2 medium (pH 6.8) supplemented with 10 g/L glucose in a shaking incubator at 37°C and 200 rpm. R/2 medium contains 2 g/L (NH4)2HPO4, 6.75 g/L KH2PO4, 0.85 g/L citric acid, 0.7 g/L MgSO4·7H2O, and 5 ml/L of a trace metal solution . Trace metal solution contains 5 M HCl, 10 g/L FeSO4·7H2O, 2.25 g/L ZnSO4·7H2O, 1 g/L CuSO4·5H2O, 0.5 g/L MnSO4·5H2O, 0.23 g/L Na2B4O7·10H2O, 2 g/L CaCl2·2H2O, and 0.1 g/L (NH4)6MO7O24.
Cell growth was monitored by measuring absorbance at 600 nm using an Ultrospec 3000 spectrophotometer (Amersham Biosciences, Uppsala, Sweden). The concentration of glucose and organic acids was determined by high-performance liquid chromatography (ProStar 210, Agilent Technologies, Santa Clara, CA, USA) with UV/visible-light (ProStar 320) and refractive index detectors (Shodex RI-71, Yokohama, Japan). A MetaCarb 87H column (300 × 7.8 mm; Varian) was eluted isocratically with 0.01 M H2SO4 at 60°C at a flow rate of 0.6 ml/minute.
Detailed procedures for genome sequencing and annotation of REL606 have been described previously . RefSeq genome data for MG1655 as of 7 June 2006 [GenBank:NC_000913] were used for comparative genomic analysis. Alignments were generated by the MUMmer 3.19 package. Genomic segments specific to each genome were identified by selecting regions >2 kb that were not aligned to each other at the equivalent position. Pseudogenes in each genome caused by IS insertion, frameshift, or in-frame stop codons were checked using the homologous gene from the other genome by BLAST or CROSS_MATCH. Genomic islands in REL606 and MG1655 were identified as described previously . Briefly, a gene was considered foreign if its G+C content (>1.5 σ) and codon usage (P < 0.05) were both aberrant. A genome scan was then performed using a ten-gene window, and regions containing four or more foreign genes were identified as genomic islands, with neighboring islands merged into larger ones.
Construction of oligonucleotide microarray
Oligonucleotide microarrays were constructed in which both REL606 and MG1655 could be tested for the synthesis of both common and unique mRNAs. Probes (70-mers) were designed by Array Designer (Premier Biosoft International, Palo Alto, CA, USA) and cross-hybridization was checked by BLASTN. A total of 4,819 probes were synthesized by Illumina. The probes represented 3,762 genes that shared over 80 nucleotides of sequence between the two E. coli strains, 382 genes found uniquely in REL606, and 534 genes unique to MG1655. The remaining probes were positive and negative controls. The probes were spotted in duplicate onto a GAPSII glass slide (Corning Inc., Corning, NY, USA) by Digital Genomics in Seoul, Korea.
The transcriptome of REL606 cells at the exponential growth stage was compared with that of MG1655 cells at the same stage, and with the transcriptome of both REL606 and MG1655 cells in stationary phase. Total RNA was isolated from cells grown in LB medium using the Qiagen RNeasy column (Valencia, CA, USA). Fluorescently labeled cDNA was synthesized by reverse transcription of total RNA with aminoallyl-dUTP and random primers followed by the coupling of Cy3 or Cy5 dyes (Amersham Pharmacia, Uppsala, Sweden). Hybridization and washes were performed according to the manufacturer's instructions (Digital Genomics, Seoul, Korea). The hybridized slides were scanned with GenePix 4000B (Molecular Devices, Sunnyvale, CA, USA) and the scanned images were analyzed by GenePix Pro 6.0. Following background subtraction, spot signal intensities were calculated as percentages of total signal intensities of Cy3 or Cy5 signals on each microarray as a means of normalization. Spots with low intensities were excluded. All microarray experiments were performed in duplicate with dye-swapping. The log2 of the intensity ratio with the two dyes for each spot was calculated from up to four values from the duplicated spot on the microarray and the dye-swap experiment. The microarray experiments were highly reproducible: the correlation coefficients between the same spots in each microarray experiment were about 0.98 and correlations between dye-swap paired microarrays were over 0.91. Genes showing growth stage-specific differences in expression levels of ≥2- or ≤0.5-fold were considered to be differentially expressed genes. The microarray data have been deposited at the Gene Expression Omnibus (GEO)  under accession number GSE13011.
Whole cellular, outer membrane, and extracellular proteins were extracted and analyzed by two-dimensional electrophoresis using the IPGphor isoelectric focusing system (GE Healthcare (Waukesha, WI, USA) and a Protean II xi electrophoresis cell (Bio-Rad, Hercules, CA, USA) as described previously  (see detailed procedures in Supplementary methods in Additional file 1). For protein identification, samples were prepared as described previously . Briefly, gel pieces were proteolysed with 0.02 μg/μl modified trypsin (Promega, Madison, WI, USA) in 40 mM ammonium bicarbonate overnight at 37°C. Tryptic peptides (10 μl aliquots) were analyzed with a nano-liquid chromatography/mass spectrometry system consisting of an Ultimate HPLC system (LC Packings, Sunnyvale, CA, USA) and a quadrupole-time-of-flight mass spectrometry (Micromass, Milford, MA, USA) equipped with a nano-electrospray ionization source, as described previously . The resulting tryptic peptide fragment data were analyzed using the MASCOT search server  for identification of protein spots. Reference databases for the identification of target proteins were an in-house E. coli B database, UniProt Knowledgebase , and NCBI .
PM tests on REL606 and MG1655 strains were carried out essentially as described previously . The PM plates (Biolog Inc., Hayward, CA, USA) consisted of 20 96-well microplates containing different sources of carbon (PM1 and PM2), nitrogen (PM3), phosphorus and sulfur (PM4), auxotrophic supplements (PM5 to PM8), or salt (PM9). The PM10 plate tests a pH stress, and plates PM11 to PM20 contain inhibitory compounds such as antibiotics, antimetabolites, and other inhibitors. Cells were grown on BUG+B agar overnight at 37°C. Colonies were picked from the agar surface and suspended in inoculating fluid containing the indicator dye tetrazolium violet. IF-0 inoculating fluid was used for plates PM1 to PM8 and IF-10 fluid for plates PM9 to PM20. Disodium succinate and ferric citrate were added to the inoculation solutions of plates PM3 to PM8. All PM plates were inoculated with cell suspensions at 100 μl/well and incubated at 37°C for 48 h in an OmniLog incubator (Biolog Inc.). For each well, the unitless areas beneath the growth-time curves for REL606 and MG1655 were determined by OmniLog-PM software and were averaged for four PM tests. These values were then subtracted from that of the negative control. Four independent tests were performed on each strain.
Metabolic network reconstruction and flux balance analysis
Based on the metabolic network model of MG1655 (iAF1260) and genome annotation of REL606 , reactions that depend on MG1655-specific genes and pseudogenes were excluded from the metabolic network model for REL606. In addition, reactions that should be turned off by REL606-specific regulation were determined using the regulatory network model for MG1655 constructed by Covert et al. . Reactions unique to REL606 were identified by genome annotation and added to the metabolic network model for it.
FBA  was used to predict cell growth on various substrates using the reconstructed metabolic network for REL606 and iAF1260 for MG1655. The predicted cell growth was compared with the results from the PM tests . Pre-specified uptake rates of each carbon source (8 mmol/gDCW/h) and oxygen (18.5 mmol/gDCW/h) were used for prediction of cell growth on PM plates testing carbon source utilization (PM1 and PM2). To predict cell growth on PM plates testing utilization of nitrogen (PM3), phosphate (PM4), and sulfur (PM4) sources, the uptake rate of each source (8 mmol/gDCW/h), succinate uptake rate (8 mmol/gDCW/h), and oxygen uptake rate (18.5 mmol/gDCW/h) were prespecified for FBA using the metabolic network models of REL606 and MG1655. The biomass objective function core  was used as an objective function. Cell growth on various media conditions was determined by the GNU linear programming kit .
flux balance analysis
Clusters of Orthologous Group of proteins
type II secretion
The authors are grateful to Jeong Im Lee for contributing to the phenotype microarray experiment and other GEM members for technical assistance; Jong Shin Yoo for helping with the mass spectrometry analysis; F William Studier for critically reading the manuscript; and Sun Chang Kim, Ohsuk Kwon, Jae-Gu Pan, and Seung-Hwan Park for helpful comments and heartful support. Special gratitude also goes to the 21C Frontier Microbial Genomics and Applications Center for financial support for fabrication of DNA microarrays. This work was supported by the National Research Foundation of Korea (2011-0017670 to JFK), the Global Frontier Intelligent Synthetic Biology Center (2011-0031952 to SHY and 2011-0031963 to SYL), the 21C Frontier Microbial Genomics and Applications Center (to JFK), and the World Class University Program (to SYL) of the Ministry of Education, Science and Technology.
- Swartz JR: Escherichia coli recombinant DNA technology. Escherichia coli and Salmonella typhimurium: Cellular and Molecular Biology. Edited by: Neidhardt FC, Curtiss III R, Ingraham JL, Lin ECC, Low KB, Magasanik B, Reznikoff WS, Riley M, Shaechter M, Umbarger HE. 1996, Washington, DC: ASM Press, 1693-1712. 2Google Scholar
- Yoon SH, Jeong H, Kwon S-K, Kim JF: Genomics, biological features, and biotechnological applications of Escherichia coli B: "Is B for better?!". Systems Biology and Biotechnology of Escherichia coli. 2009, Berlin, Germany: Springer, 1-17.View ArticleGoogle Scholar
- Blattner FR, Plunkett G 3rd, Bloch CA, Perna NT, Burland V, Riley M, Collado-Vides J, Glasner JD, Rode CK, Mayhew GF, Gregor J, Davis NW, Kirkpatrick HA, Goeden MA, Rose DJ, Mau B, Shao Y: The complete genome sequence of Escherichia coli K-12. Science. 1997, 277: 1453-1474. 10.1126/science.277.5331.1453.PubMedView ArticleGoogle Scholar
- Hayashi K, Morooka N, Yamamoto Y, Fujita K, Isono K, Choi S, Ohtsubo E, Baba T, Wanner BL, Mori H, Horiuchi T: Highly accurate genome sequences of Escherichia coli K-12 strains MG1655 and W3110. Mol Syst Biol. 2006, 2: 2006 0007Google Scholar
- Barrick JE, Yu DS, Yoon SH, Jeong H, Oh TK, Schneider D, Lenski RE, Kim JF: Genome evolution and adaptation in a long-term experiment with Escherichia coli. Nature. 2009, 461: 1243-1247. 10.1038/nature08480.PubMedView ArticleGoogle Scholar
- Jeong H, Barbe V, Lee CH, Vallenet D, Yu DS, Choi S-H, Couloux A, Lee S-W, Yoon SH, Cattolico L, Hur CG, Park H-S, Segurens B, Kim SC, Oh TK, Lenski RE, Studier FW, Daegelen P, Kim JF: Genome sequences of Escherichia coli B strains REL606 and BL21(DE3). J Mol Biol. 2009, 394: 644-652. 10.1016/j.jmb.2009.09.052.PubMedView ArticleGoogle Scholar
- Studier FW, Daegelen P, Lenski RE, Maslov S, Kim JF: Understanding the differences between genome sequences of Escherichia coli B strains REL606 and BL21(DE3) and comparison of the E. coli B and K-12 genomes. J Mol Biol. 2009, 394: 653-680. 10.1016/j.jmb.2009.09.021.PubMedView ArticleGoogle Scholar
- National Research Council: Sequence-based Classification of Select Agents: a Brighter Line. 2010, Washington, DC: The National Academies Press, [http://www.nap.edu/openbook.php?record_id=12970&page=1]Google Scholar
- Ideker T, Thorsson V, Ranish JA, Christmas R, Buhler J, Eng JK, Bumgarner R, Goodlett DR, Aebersold R, Hood L: Integrated genomic and proteomic analyses of a systematically perturbed metabolic network. Science. 2001, 292: 929-934. 10.1126/science.292.5518.929.PubMedView ArticleGoogle Scholar
- Yoon SH, Han MJ, Lee SY, Jeong KJ, Yoo JS: Combined transcriptome and proteome analysis of Escherichia coli during high cell density culture. Biotechnol Bioeng. 2003, 81: 753-767. 10.1002/bit.10626.PubMedView ArticleGoogle Scholar
- Bochner BR: New technologies to assess genotype-phenotype relationships. Nat Rev Genet. 2003, 4: 309-314.PubMedView ArticleGoogle Scholar
- Schmid AK, Reiss DJ, Kaur A, Pan M, King N, Van PT, Hohmann L, Martin DB, Baliga NS: The anatomy of microbial cell state transitions in response to oxygen. Genome Res. 2007, 17: 1399-1413. 10.1101/gr.6728007.PubMedPubMed CentralView ArticleGoogle Scholar
- Covert MW, Knight EM, Reed JL, Herrgard MJ, Palsson BO: Integrating high-throughput and computational data elucidates bacterial networks. Nature. 2004, 429: 92-96. 10.1038/nature02456.PubMedView ArticleGoogle Scholar
- Bonneau R, Facciotti MT, Reiss DJ, Schmid AK, Pan M, Kaur A, Thorsson V, Shannon P, Johnson MH, Bare JC, Longabaugh W, Vuthoori M, Whitehead K, Madar A, Suzuki L, Mori T, Chang DE, Diruggiero J, Johnson CH, Hood L, Baliga NS: A predictive model for transcriptional control of physiology in a free living cell. Cell. 2007, 131: 1354-1365. 10.1016/j.cell.2007.10.053.PubMedView ArticleGoogle Scholar
- Feist AM, Henry CS, Reed JL, Krummenacker M, Joyce AR, Karp PD, Broadbelt LJ, Hatzimanikatis V, Palsson BO: A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information. Mol Syst Biol. 2007, 3: 121-PubMedPubMed CentralView ArticleGoogle Scholar
- Eiteman MA, Altman E: Overcoming acetate in Escherichia coli recombinant protein fermentations. Trends Biotechnol. 2006, 24: 530-536. 10.1016/j.tibtech.2006.09.001.PubMedView ArticleGoogle Scholar
- Phue JN, Noronha SB, Hattacharyya R, Wolfe AJ, Shiloach J: Glucose metabolism at high density growth of E. coli B and E. coli K: differences in metabolic pathways are responsible for efficient glucose utilization in E. coli B as determined by microarrays and Northern blot analyses. Biotechnol Bioeng. 2005, 90: 805-820. 10.1002/bit.20478.PubMedView ArticleGoogle Scholar
- Diaz-Guerra M, Esteban M, Martinez JL: Growth of Escherichia coli in acetate as a sole carbon source is inhibited by ankyrin-like repeats present in the 2',5'-linked oligoadenylate-dependent human RNase L enzyme. FEMS Microbiol Lett. 1997, 149: 107-113. 10.1016/S0378-1097(97)00064-5.PubMedView ArticleGoogle Scholar
- Umbarger HE: Biosynthesis of the branched-chain amino acids. Escherichia coli and Salmonella typhimurium: Cellular and Molecular Biology. Edited by: Neidhardt FC, Curtiss III R, Ingraham JL, Lin ECC, Low KB, Magasanik B, Reznikoff WS, Riley M, Schneider D, Umbarger HE. 1996, Washington, DC: ASM Press, 442-457. 2Google Scholar
- Herrera G, Martinez A, Blanco M, O'Connor JE: Assessment of Escherichia coli B with enhanced permeability to fluorochromes for flow cytometric assays of bacterial cell function. Cytometry. 2002, 49: 62-69. 10.1002/cyto.10148.PubMedView ArticleGoogle Scholar
- Jansson PE, Lindberg AA, Lindberg B, Wollin R: Structural studies on the hexose region of the core in lipopolysaccharides from Enterobacteriaceae. Eur J Biochem. 1981, 115: 571-577.PubMedView ArticleGoogle Scholar
- Plumbridge JA, Cochet O, Souza JM, Altamirano MM, Calcagno ML, Badet B: Coordinated regulation of amino sugar-synthesizing and -degrading enzymes in Escherichia coli K-12. J Bacteriol. 1993, 175: 4951-4956.PubMedPubMed CentralGoogle Scholar
- MacNab RM: Flagella and motility. Escherichia coli and Salmonella typhimurium: Cellular and Molecular Biology. Edited by: Neidhardt FC, Curtiss III R, Ingraham JL, Lin ECC, Low KB, Magasanik B, Reznikoff WS, Riley M, Schneider D, Umbarger HE. 1996, Washington, DC: ASM Press, 123-145. 2Google Scholar
- Schneider D, Duperchy E, Depeyrot J, Coursange E, Lenski R, Blot M: Genomic comparisons among Escherichia coli strains B, K-12, and O157:H7 using IS elements as molecular markers. BMC Microbiol. 2002, 2: 18-10.1186/1471-2180-2-18.PubMedPubMed CentralView ArticleGoogle Scholar
- Nikaido H: Molecular basis of bacterial outer membrane permeability revisited. Microbiol Mol Biol Rev. 2003, 67: 593-656. 10.1128/MMBR.67.4.593-656.2003.PubMedPubMed CentralView ArticleGoogle Scholar
- Majdalani N, Gottesman S: The Rcs phosphorelay: a complex signal transduction system. Annu Rev Microbiol. 2005, 59: 379-405. 10.1146/annurev.micro.59.050405.101230.PubMedView ArticleGoogle Scholar
- Gerdes K, Christensen SK, Lobner-Olesen A: Prokaryotic toxin-antitoxin stress response loci. Nat Rev Microbiol. 2005, 3: 371-382. 10.1038/nrmicro1147.PubMedView ArticleGoogle Scholar
- saiSree L, Reddy M, Gowrishankar J: IS186 insertion at a hot spot in the lon promoter as a basis for Lon protease deficiency of Escherichia coli B: identification of a consensus target sequence for IS186 transposition. J Bacteriol. 2001, 183: 6943-6946. 10.1128/JB.183.23.6943-6946.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Seo JH, Kang DG, Cha HJ: Comparison of cellular stress levels and green-fluorescent-protein expression in several Escherichia coli strains. Biotechnol Appl Biochem. 2003, 37: 103-107. 10.1042/BA20020041.PubMedView ArticleGoogle Scholar
- Ruiz N, Silhavy TJ: Sensing external stress: watchdogs of the Escherichia coli cell envelope. Curr Opin Microbiol. 2005, 8: 122-126. 10.1016/j.mib.2005.02.013.PubMedView ArticleGoogle Scholar
- Francetic O, Belin D, Badaut C, Pugsley AP: Expression of the endogenous type II secretion pathway in Escherichia coli leads to chitinase secretion. EMBO J. 2000, 19: 6697-6703. 10.1093/emboj/19.24.6697.PubMedPubMed CentralView ArticleGoogle Scholar
- Jensen KF: The Escherichia coli K-12 "wild types" W3110 and MG1655 have an rph frameshift mutation that leads to pyrimidine starvation due to low pyrE expression levels. J Bacteriol. 1993, 175: 3401-3407.PubMedPubMed CentralGoogle Scholar
- Joyce AR, Palsson BO: The model organism as a system: integrating 'omics' data sets. Nat Rev Mol Cell Biol. 2006, 7: 198-210. 10.1038/nrm1857.PubMedView ArticleGoogle Scholar
- Lee SY, Chang HN: High cell density cultivation of Escherichia coli W using sucrose as a carbon source. Biotechnol Lett. 1993, 15: 971-974. 10.1007/BF00131766.View ArticleGoogle Scholar
- Yoon SH, Hur CG, Kang HY, Kim YH, Oh TK, Kim JF: A computational approach for identifying pathogenicity islands in prokaryotic genomes. BMC Bioinformatics. 2005, 6: 184-10.1186/1471-2105-6-184.PubMedPubMed CentralView ArticleGoogle Scholar
- Gene Expression Omnibus (GEO). [http://www.ncbi.nlm.nih.gov/geo/]
- Lee JW, Lee SY, Song H, Yoo JS: The proteome of Mannheimia succiniciproducens, a capnophilic rumen bacterium. Proteomics. 2006, 6: 3550-3566. 10.1002/pmic.200500837.PubMedView ArticleGoogle Scholar
- Han MJ, Yoon SS, Lee SY: Proteome analysis of metabolically engineered Escherichia coli producing Poly(3-hydroxybutyrate). J Bacteriol. 2001, 183: 301-308. 10.1128/JB.183.1.301-308.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- MASCOT version 1.8. [http://www.matrixscience.com/]
- Swiss-Prot and TrEMBL. [http://kr.expasy.org/]
- National Center for Biotechnology Information. [http://www.ncbi.nlm.nih.gov/]
- Zhou L, Lei XH, Bochner BR, Wanner BL: Phenotype microarray analysis of Escherichia coli K-12 mutants with deletions of all two-component systems. J Bacteriol. 2003, 185: 4956-4972. 10.1128/JB.185.16.4956-4972.2003.PubMedPubMed CentralView ArticleGoogle Scholar
- Price ND, Reed JL, Palsson BO: Genome-scale models of microbial cells: evaluating the consequences of constraints. Nat Rev Microbiol. 2004, 2: 886-897. 10.1038/nrmicro1023.PubMedView ArticleGoogle Scholar
- GNU Linear Programming Kit (GLPK). [http://www.gnu.org/software/glpk/]
- Clusters of Orthologous Groups of proteins (COGs). [http://www.ncbi.nlm.nih.gov/COG/]
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.