- Open Access
Comparative multi-omics systems analysis of Escherichia coli strains B and K-12
Genome Biology volume 13, Article number: R37 (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.
The average nucleotide identity of the aligned genomic regions of REL606 [GenBank:NC_012967]  and MG1655 [GenBank:NC_000913]  was >99.1%. Only about 4% of the total genome accounts for strain-specific regions, including prophages and seemingly recently transferred genomic islands (Figure 1). Interestingly, the B genome encodes an additional set of genes for type II secretion (T2S) and D-arabinose utilization, and lacks the gene cluster for flagellar biosynthesis and the very short-patch repair system. Different sets of genes were observed for the Qin prophage, O antigen biosynthesis, catabolism of aromatic compounds, and lipopolysaccharide (LPS) oligosaccharide biosynthesis. Different gene clusters for the catabolism of aromatic compounds were also detected; the hpa cluster for degradation of 3- and 4-hydroxy phenyl acetic acid in the B strain and the paa cluster for catabolism of phenyl acetic acid in the K-12 strain. We found numerous strain-specific gene disruptions caused by deletions, frameshifts, or insertions of insertion sequence (IS) elements that might affect bacterial phenotypes (Table S1 in Additional file 1).
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.
Transcriptome profiles were analyzed using samples taken in the exponential and stationary growth phases during cultivation of REL606 and MG1655 in LB medium. During the exponential growth phase, the most highly expressed genes in the B strain were those involved in replication, translation, or nucleotide transport and metabolism, while many of the highly expressed genes in the K-12 strain were involved in cell motility, transcription, carbohydrate transport, or energy production (Figure 2; Table S2 in Additional file 1). At the stationary phase, many of the genes that were highly expressed in the B strain were involved in transport and metabolism of various amino acids and carbohydrates, whereas highly expressed genes in the K-12 strain had functions related to cell motility, ribosomal subunit protein production, or energy generation.
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.
Examination of the extracellular proteome (Figure S2C in Additional file 1) revealed that B strains released larger amounts of protein in the stationary growth phase than K-12 strains (Figure 3). Proteins commonly released from both B and K-12 strains were outer membrane proteins (OmpA, OmpF, and OmpX) and transport/binding proteins (ArgT, ArtI, ArtJ, CusF, GlnH, HisJ, and LivJ). The outer membrane proteins OmpF and YaeT, the periplasmic substrate transport and binding proteins CusF and RbsB, and the cold shock protein CspC were released at higher levels from B than from K-12 strains, whereas motility-related proteins (FimA, FlgF, FlgL, FliC, and FliD) and Flu were secreted only by K-12 strains (top of Figure 4).
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
Major differences were observed in the expression of genes involved in the biosynthesis and degradation of amino acids (Figure 4). Most of the genes involved in biosynthetic pathways of amino acids, especially those required for biosynthesis of arginine and branched-chain amino acids (leucine, isoleucine, and valine), were more highly expressed in the B strain than in K-12. Genes encoding the three isozymes of acetohydroxy acid synthase I (ilvBN), II (ilvGM) and III (ilvIH) exhibited complicated expression patterns (Figure 4), which may be attributed to transcriptional control by attenuation . A frameshift mutation was found in the ilvG gene of the K-12 strains, whereas this gene is intact in the B strains. The gene ilvG mediates valine resistance and allows growth on various valine dipeptides in PM tests (Figure S3 in Additional file 1); this explains the failure of MG1655 to grow on valine dipeptides. Proteome analysis also indicated that enzymes required for the biosynthesis of arginine (ArgDCI), serine (SerC), and aspartate (AspC) were synthesized at high levels in the B strains. Interestingly, TnaA, a protein involved in the tryptophan and cysteine degradation pathway, was highly expressed at both the mRNA and protein levels in the K-12 strains. The same was true for TdcE and TdcF, required for anaerobic degradation of L-threonine, GadB, involved in glutamate degradation and acid resistance, and AspA, needed for the reversible conversion of L-aspartate to fumarate and ammonia. These observations suggest that the B strains are, in general, disposed to enhanced amino acid biosynthesis and reduced degradation, conditions favorable for efficient production of recombinant proteins (Figure 5).
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.
In silico complementation experiments of the reconstructed model verified 23 reactions responsible for the phenotypes unique to REL606 (Table 1). The model can be further applied to test hypotheses derived from sequence and expression data; for example, slow growth of K-12 in pyrimidine-free media can be attributed to the frame-shift mutation of rph, thereby decreasing the expression of pyrE . FBA results demonstrated that the flux level associated with pyrE shows linear correlation with the specific growth rate. Indeed, the rph gene was found to be intact in B strains, which grow better than K-12 strains in minimal media (Figure S1B in Additional file 1). Reconstruction of a computational metabolic network of E. coli B and in silico complementation experiments revealed genetic bases behind the phenotypic characteristics of B.
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 T6r strr 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
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. 2
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.
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.
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 0007
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.
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.
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.
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]
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.
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.
Bochner BR: New technologies to assess genotype-phenotype relationships. Nat Rev Genet. 2003, 4: 309-314.
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.
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.
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.
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-
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.
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.
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.
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. 2
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.
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.
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.
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. 2
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.
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.
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.
Gerdes K, Christensen SK, Lobner-Olesen A: Prokaryotic toxin-antitoxin stress response loci. Nat Rev Microbiol. 2005, 3: 371-382. 10.1038/nrmicro1147.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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/]
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.
The authors declare that they have no competing interests.
JFK conceived, organized and supervised the project, and analyzed the omics information; SHY designed DNA microarray and phenotype microarray experiments, and analyzed and accessed the biological significance of genome, transcriptome, proteome, and phenome data; M-JH and X-XX designed and performed proteome experiments, and analyzed the data; HJ contributed to analyzing the genomic differences between B and K-12; CHL reconstructed the metabolic model of B and performed in silico complementation tests; JHS performed the phylogenetic analysis of type II secretion systems; SYL and TKO participated in conceiving the project and analyzing the biological implications of omics data; SHY and M-JH wrote the manuscript; and JFK, SYL, and D-HL revised the draft. All authors have read and approved the manuscript for publication.
Electronic supplementary material
Additional file 1: . Supplementary Methods: preparation of whole cellular, membrane and extracellular proteins, two-dimensional gel electrophoresis and image analysis, and metabolite abbreviations in Figure 4. Figure S1: growth curves of E. coli B strains (REL606 and BL21(DE3) and K-12 strains (MG1655 and W3110) in complex LB medium and minimal R/2 medium supplemented with 10 g/l glucose. Figure S2: total cellular, outer membrane, and extracellular proteomes of the E. coli strains. Figure S3: phenotype microarray comparison of E. coli B REL606 and K-12 MG1655. Figure S4: phylogenetic position of two T2S systems in E. coli B REL606. Table S1: pseudogene comparison between E. coli B REL606 and K-12 MG1655. Table S2: genes that were highly expressed at both the exponential and stationary growth phases during growth of E. coli B REL606 and K-12 MG1655 in LB medium. Table S3: proteins exhibiting significant quantitative differences between E. coli B and K-12 strains. Table S4: metabolic reactions modified in the metabolic network model for E. coli B REL606 compared to the model for E. coli K-12 MG1655. Table S5: phenotypic differences of E. coli B REL606 and K-12 MG1655 in PM1 and PM2 and in silico predictions of cell growth on each carbon source. Supplementary References. (PDF 2 MB)
About this article
Cite this article
Yoon, S.H., Han, MJ., Jeong, H. et al. Comparative multi-omics systems analysis of Escherichia coli strains B and K-12. Genome Biol 13, R37 (2012). https://doi.org/10.1186/gb-2012-13-5-r37
- Outer Membrane Protein
- Flux Balance Analysis
- Genomic Island
- Acetate Accumulation
- Phenotype Microarray