Previous Article | Next Article ![]()
Journal of Bacteriology, September 2003, p. 5442-5451, Vol. 185, No. 18
0021-9193/03/$08.00+0 DOI: 10.1128/JB.185.18.5442-5451.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Dong-Eun Lee,1 Bheong-Uk Lee,2 and Hak-Sung Kim1*
Department of Biological Sciences, Korea Advanced Institute of Science and Technology, 373-1, Kusung-dong, Yusung-gu, Taejon, 305-701,1 Department of Biological Science, Kosin University, Busan, Korea2
Received 5 November 2002/ Accepted 26 June 2003
|
|
|---|
|
|
|---|
Amino acids have been the prominent target metabolites from microorganisms in bioindustry due to large commercial demands for flavor enhancers, animal feed, sweeteners, and therapeutic agents. Of them, L-threonine, one of the essential amino acids, is widely used as a feed and food additive, and various industrial strains that more efficiently produce L-threonine have been successfully developed by traditional approaches, including the deregulation of enzymes, the elimination of competitive pathways, and the amplification of genes (7, 20, 26, 36). In a previous study, we developed the L-threonine-producing Escherichia coli strain TF5015 via recursive mutations (24). This strain requires both L-methionine and L-isoleucine for growth and shows resistance to various chemical analogues. It was suggested that the mechanism of L-threonine production of TF5015 probably resulted from releasing the feedback regulation and blocking the carbon flow into undesirable by-products. Regulation mechanism involved in L-threonine biosynthesis in E. coli has been relatively well characterized (30). However, to understand the overall regulatory mechanism and the physiological events in response to the accumulation of L-threonine in TF5015, detailed information, including the cellular regulations, entire metabolic fluxes, and dynamic responses of the complex metabolic networks is crucial.
In the present study, to gain further insight into the global regulatory mechanism for L-threonine biosynthesis, we carried out the comparative analyses of transcriptome, proteome, and nucleotide sequences between the prototrophic E. coli W3110 and the L-threonine-producing E. coli TF5015. Expression patterns of the genes and proteins were investigated for both strains by using DNA macroarrays containing virtually every gene of E. coli and two-dimensional gel electrophoresis. The profiles were analyzed in terms of the accumulation of L-threonine and physiological consequences in the mutant strain. The details of these analyses are presented below.
|
|
|---|
-amino-ß-hydroxyvaleric acid (AHV; threonine analogue), 3 mg of 2-aminoethyl-L-cysteine (AEC; lysine analogue), and 1 mg of L-azetidine-2-carboxylic acid (AC; proline analogue)/ml (Fig. 1) (24). This strain requires both methionine and isoleucine for growth. E. coli Hfr 3000 YA73 (thrB relA1 spoT1 thi-1) was used for cloning of the thr operon (39).
![]() View larger version (21K): [in a new window] |
FIG. 1. Genealogy of E. coli mutants producing L-threonine from E. coli W3110. Abbreviations: Met-, L-methionine auxotroph; AHVr, resistant to -amino-ß-hydroxyvaleric acid; AECr, resistant to 2-aminoethyl-L-cysteine; IleL, L-isoleucine leaky character; ACr, resistant to L-azetidine-2-carboxylic acid. The production of L-threonine was conducted in a 250-ml flask (24).
|
RNA purification. For analysis of transcription levels between the prototrophic strain and threonine-producing mutant at the threonine production stage, both strains were cultivated to an early stationary phase in the fermentation medium (Fig. 2), and culture broth was taken for extraction and purification of total RNA by using RNeasy kits (Qiagen) according to the procedures reported elsewhere (1, 17, 32, 38). Briefly, 0.2 ml of culture broth was lysed in boiling lysis buffer (1), and the resulting solution was extracted twice by using acidic hot phenol (pH 4.3) at 65°C with vigorous vortex and incubation for 5 min, followed by acid phenol-chloroform (5:1) extraction. The RNA was precipitated with ethanol, redissolved in water, treated with DNase I (Ambion), and applied to RNeasy column. The purified RNA was redissolved in RNase-free water and stored at -72°C. The absence of genomic DNA was confirmed according to the reported procedure (32). The RNA concentration was determined by using the spectrophotometer.
![]() View larger version (26K): [in a new window] |
FIG. 2. Growth curve of E. coli W3110 and TF5015 in the fermentation medium. Cells were harvested at the time indicated in arrows.
|
-33P]dCTP (
2,500 Ci/mmol; Amersham-Pharmacia Biotech) was added to the cDNA labeling mixture, followed by incubation for 2 h 30 min at 42°C. Unincorporated-radiolabeled nucleotides were removed by applying the reaction mixture to a Sephadex G-25 gel filtration spin column. Hybridizations were performed in roller bottles (3.5 by 30 cm) in a hybridization oven. After the blots were rinsed in 50 ml of 2x SSPE (1x SSPE is 0.18 M NaCl, 10 mM NaH2PO4, and 1 mM EDTA [pH 7.7]) for 5 min, the Panorama gene arrays (Sigma-GenoSys) were prehybridized for 2 h at 65°C in 5 ml of hybridization solution (17, 32, 38). The radiolabeled cDNA mixture was added to the hybridization solution and preincubated at 90 to 95°C for 10 min in a water bath. The denatured labeled cDNA in hybridization solution was added to the blots in the roller bottles and hybridized at 65°C for 18 h. After the hybridization solution was decanted, the array membranes were first washed with 50 ml of wash solution three times for 5 min each time at room temperature, followed by a three 20-min washes with 100 ml of prewarmed wash solution at 65°C. The membranes were then wrapped in clear wrap and exposed to a phosphorimager screen (BAS-SR2025 imaging plate [20 by 40 cm]; Fuji Photo Film Co., Ltd.) for 24 to 72 h. The exposed imaging screens were scanned with a pixel size of 50 µm for greater resolution of spots on an FLA 3000 phosphorimager (Fuji).
Analysis of DNA arrays. The spot intensities were analyzed by using ArrageGauge software from Fuji. Background values were automatically subtracted by using the global subtraction method. The density values for the duplicate spots were averaged and then normalized by the global normalization method described in the software. In this case, the sums of total spot densities, including data and controls, were the same in the two arrays. Identification of each gene and its functional grouping was carried out by using the manufacturer's E. coli array information file. The correlation coefficients of spot intensities were calculated from the duplicate spots on the same membrane and from the duplicate experiments with total RNA purified from TF5015, and the coefficients ranged from 0.997 to 0.999 and from 0.910 to 0.980, respectively. The expression ratio (relative intensity of genes in TF5015 to those in W3110) for genes showing confidence levels higher than 99.9% in two sets of experiments was 1.78; thus, the genes showing a level greater than this value were considered to be upregulated. The genes with an expression ratio less than 0.56 were regarded as downregulated. In addition to above criteria, genes whose signal intensities were lower than the mean background intensity plus three standard deviations in both W3110 and TF5015 samples were excluded (32).
Two-dimensional gel electrophoresis and mass spectrometric analysis. Two-dimensional gel electrophoresis was carried out according to a procedure described elsewhere (13, 14, 17) with a slight modification. For protein extraction, cells grown on fermentation medium were harvested at the same growth phase as for total RNA preparation. A 2-ml portion of culture broth was centrifuged at 2,500 x g and 4°C for 5 min and then washed four times with a wash solution (14). The resulting pellet was resuspended in 10 mM Tris buffer (pH 8.0) containing 1 mM phenylmethylsulfonyl fluoride, sonicated for 2 min on ice, and centrifuged at 10,000 x g for 30 min. The supernatant was treated with 1 U of DNase I at 37°C for 15 min, centrifuged for 30 min at 10,000 x g and 4°C, and dialyzed in 5 mM Tris buffer overnight at 4°C. After protein quantification, 30-, 50-, and 100-µg portions of protein extracts were dried by vacuum centrifugation and stored at -72°C until use.
Extracted protein dissolved in 400 µl of rehydration solution (8 M urea; CHAPS {3-[(3-cholamidopropyl)-dimethylammonio]-1-propanesulfonate}; 18 mM dithiothreitol; 0.5% IPG buffer at pH 4 to 7 L or pH 3 to 10 NL [Amersham-Pharmacia Biotech]) was loaded onto 24-cm IPG strips (pH 4 to 7 L or pH 3 to 10 NL type [Amersham-Pharmacia Biotech]). After rehydration of strips for 12 h, proteins were focused by using a series of voltage increases at 500 V for 1 h, 1,000 V for 1 h, and 8,000 V for 12 h. The second dimension was carried out in a 12.5% polyacrylamide gel (Ettan DALT gel, 255 by 196 by 1 mm [Amersham-Pharmacia Biotech]). After fixation and silver staining, the wet gel was scanned with an ImageScanner and quantified with ImageMaster 2D Elite software (Amersham Pharmacia Biotech). Spot densities were determined on three gel images from duplicate cultures and then normalized to the total spot volume of the protein spots on each gel.
Protein spots showing distinctly different expression patterns were excised and destained by a silver stain-destain protocol (Scripps Center for Mass Spectrometry [http://masspec.scripps.edu/services/protein/]) and subjected to gel digestion with 20 ng of sequencing-grade trypsin (Boehringer Mannheim)/µl in 8 µl of 2.5 mM Tris-HCl at pH 8.5 (13, 14, 35). The resulting samples were analyzed by matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) mass spectrometry (PE Biosystems Voyager System 4095 [PerSeptive Biosystems, Inc.]). Peptide fragment peaks (at 805.416 and 2,163.056 m/z) produced by autodigestion of trypsin were used as an internal calibration. Proteins were identified by using the ProteinProspector server (http://prospector.ucsf.edu/) and by comparison with the E. coli reference gel in the SWISS-2DPAGE database (http://kr.expasy.org/ch2d/).
Recombinant DNA techniques. DNA manipulations and PCR were carried out by standard procedures (34). Chromosomal DNA and PCR fragments were purified by using a Wizard genomic DNA purification kit (Promega) and a QIAquick gel extraction kit (Qiagen), respectively. Fifteen genes were amplified by PCR with Pfu DNA polymerase (Strategene) from genomic DNA of E. coli W3110 and TF5015. The nucleotide sequence was determined by the dideoxy chain termination method (34).
Cloning of thrABC operon. The wild-type and mutated thrABC operons were cloned in a pBR322 plasmid (New England Biolabs) by complementation of E. coli mutant. Genomic DNAs of E. coli W3110 and TF5015 were completely digested with HindIII and BamHI, and DNA fragments ranging from 5 to 7 kb were isolated, respectively. The digested fragments were ligated with 4.3-kb pBR322/HindIII/BamHI and transformed into thrB auxotrophic E. coli Hfr 3000 YA73 (39). The plasmids, pAW88 (10.5kb) and pATF92 (10.5kb) containing wild-type and mutated thr operons, respectively, were obtained from thr+ transformants. The presence of thrABC genes and its regulatory region in pAW88 and pATF92 was confirmed by DNA sequencing. The constructed plasmids were introduced into W3110 by electroporation (34).
Analysis. Cell growth was monitored by measuring the optical density at 562 nm by using a spectrophotometer (Beckman DU650). Amino acids were analyzed by high-pressure liquid chromatography (Waters). Acetic acid was determined by using a DX-600 ion chromatograph (Dionex).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Genes showing a significant differential expression
|
|
View this table: [in a new window] |
TABLE 2. Functional classification of differentially expressed genes
|
![]() View larger version (27K): [in a new window] |
FIG. 3. Metabolic map showing relative expression ratios (intensity of TF5015 versus that of W3110) of transcript levels in the central metabolic pathway (A) and the threonine biosynthetic pathway (B). Abbreviations: GLC, glucose; PEP, phosphoenolpyruvate; PYR, pyruvate; LAC, lactate; Ac-CoA, acetyl-CoA; AC-P, acetyl-phosphate; ACE, acetate; CIT, citrate; ICT, isocitrate; -KG, -ketoglutarate; SUC-CoA, succinyl-CoA; FUM, fumarate; MAL, malate; OAA, oxaloacetate; GOL, glyoxylate; GLU, glutamate; GLN, glutamine; PRO, proline; ASP, aspartate; LYS, lysine; MET, methionine; ß-AS-P, ß-aspartyl phosphate; ASS, aspartate semialdehyde; HOM, homoserine; HOM-P, homoserine phosphate; THR, threonine; GLY, glycine; ILE, isoleucine; LEU, leucine; VAL, valine. The shaded and boxed numbers in the figure represent up- and downregulated genes, respectively, in TF5015.
|
|
View this table: [in a new window] |
TABLE 3. Accumulation levels of amino acids and acetic acid in E. coli W3110 and TF5015 at the early stationary phase
|
4-fold lower than in W3110 at the early stationary phase (Table 3). Many ribosomal genes were downregulated in TF5015 (Table 1). Proteome analysis between W3110 and TF5015. In order to analyze the global gene expression at protein levels and to confirm whether the transcriptome analysis using the DNA arrays would correlate with the protein expression levels, two-dimensional gel electrophoresis was carried out for proteins extracted from W3110 and TF5015 at the same growth stage as for total RNA preparation. As shown in Fig. 4, ca. 500 protein spots were visualized by silver staining. Of them, 18 protein spots displayed significantly different expression patterns; 10 proteins were found to be highly expressed, whereas 8 proteins were downregulated in TF5015 (Fig. 4A, 4B). We identified 14 proteins from 18 spots by using MALDI-TOF (Fig. 4C). Four proteinsAceA, IcdA, OppA, and ThrCappeared as a doublet spot in two-dimensional gel. Proteins showing a higher expression level in TF5015 were identified as AceA, AldA, IcdA, LeuC, LeuD, ThrC, and Udp, whereas the downregulated proteins were ArgG, OppA, and YfiD. Of the identified 10 proteins, 8 proteins were in good agreement with the expression profiles of the corresponding genes, and these include AceA, AldA, IcdA, LeuC, LeuD, OppA, ThrC, and YfiD. However, the protein levels of ArgG and Udp were not directly correlated with the expression levels of each transcript.
![]() ![]() ![]() View larger version (240K): [in a new window] |
FIG. 4. Silver-stained two-dimensional gel patterns of E. coli W3110 (A) and E. coli TF5015 (B) at the same growth stage as DNA macroarray and quantification of protein spots showing differential expression patterns (C). In panels A and B, the horizontal axes represent pH, and the vertical axes represent molecular masses in kilodaltons. The transcript fold values (i.e., the relative intensity of transcript in TF5015 versus that of W3110) of corresponding genes are shown below the differently expressed proteins in panel C. 18 proteins (10 spots [ ] of higher expression in TF5015 and 8 spots [ ] of higher expression in W3110) showing differential expression patterns were selected, and 14 protein spots were identified by MALDI-TOF and shown in (C). Four proteins (IcdA, ThrC, AceA, and OppA) displayed double spots in the two-dimensional gels, respectively. The symbol ( ) indicates unidentified protein spots. Spot intensities were measured and normalized as described in Materials and Methods. Error bars represent the standard deviation of the mean intensity. Continued.
|
Comparison of nucleotide sequences between W3110 and TF5015. Analysis of nucleotide sequences between W3110 and TF5015 is expected to provide some insights into regulation mechanism of genes. The analyses of transcriptome and proteome between the two strains revealed that a number of the genes showing the significant expression changes in TF5015 are mainly related to the metabolic pathways affecting directly the biosynthesis and/or metabolism of threonine. Thus, we selected 15 genes in the central metabolic pathway and threonine biosynthetic pathway, and analyzed their nucleotide sequences. As a result, it was found that mutations occurred in the open reading frames (ORFs) of thrA and ilvA (Table 4) in TF5015. The mutation sites of thrA and ilvA genes were S345F (designated thrA345) and S97F (designated ilvA97), respectively. No mutation was observed in other 13 genes.
|
View this table: [in a new window] |
TABLE 4. Changes in nucleotide sequences of genes in E. coli TF5015 compared to strain W3110
|
|
View this table: [in a new window] |
TABLE 5. Production of L-threonine by recombinant E. coli strain
|
|
|
|---|
Comparison of gene expression profiles between W3110 and TF5015 showed that only 54 of 4,290 genes (1.3%) exhibit differential transcript expression patterns, and this is interesting since TF5015 is the mutant producing a much higher level of threonine compared to the parent strain W3110. From the viewpoint of growth rate, metabolic fluxes, regulation mechanisms, genetic modifications, and accumulation of threonine, although the physiology of TF5015 would greatly differ from that of W3110, global comparisons of the expression profiles suggest that the general trends at the transcriptional level are very similar to each other.
Both transcriptome and proteome analyses revealed the upregulation of glyoxylate shunt and TCA cycle in TF5015 compared to W3110 (Fig. 3). It was reported that the supply of oxaloacetate from PEP via PEP carboxylase and glyoxylate bypass is required for maximal yield of threonine by theoretical stoichiometric analysis (40) and experimental evidence (37). Since the biosynthesis of oxaloacetate via PEP carboxylase is main route from glucose and PEP carboxylase is inhibited by aspartate in E. coli (11), the regulation of ppc and its gene product is the pivotal step for biosynthesis of aspartate family amino acids. When we analyzed the ppc gene at genetic and transcriptional levels between W3110 and TF5015, distinctive differences were not observed (Fig. 3 and Table 4). In this regard, our analyses strongly imply that TF5015 adopts glyoxylate shunt and TCA cycle rather than PEP carboxylase for glucose metabolism at the transcript level.
Pyruvate formation is a main route from glucose metabolism according to stoichiometric analysis of central metabolism in E. coli (15); thus, it seems that the conversion of pyruvate to oxaloacetate without losing carbon as CO2 can be a major factor for improving the production of threonine. Therefore, it seems to be reasonable that upregulation of the ace operon could partly contribute to an increase of acetyl coenzyme A (acetyl-CoA) flux into oxaloacetate and increased production of threonine, even though the precise balanced flow of isocitrate at the branch point is not clear in TF5015. It was known that E. coli could excrete 10 to 30% of carbon flux from glucose as acetate under aerobic conditions (16). Accordingly, we suggest that upregulation of the ace operon can be a way of preventing the accumulation of acetate or accommodating an efficient utilization of acetate (8). This is well supported by a reduced level of acetate in TF5015 compared to W3110 (Table 3). Since glyoxylate shunt is normally repressed during the growth of E. coli on glucose, upregulation of the ace operon is not consistent with the previous observation (5). The expression of aceBAK is known to be significantly influenced by various regulation factors, including arcA, arcB, fadR, fruR, himA, himD, and iclR (5), but expression levels of these factors were almost the same between W3110 and TF5015. In addition, no modifications of DNA sequences in an ORF of aceBAK, an upstream region of aceBAK promoter, and an iclR in both strains were detected by DNA sequencing (Table 4). Thus, some other factors, such as a change of metabolites responding to threonine production or unknown mechanisms, seem to affect the transcriptional control of the ace operon.
Upregulation of the thr operon in TF5015 was confirmed by both transcriptome and proteome analyses (Fig. 3 and 4). In a previous study, we observed that the expression level of aspartokinase I (AKI) in TF427, the parent strain of TF5015, was elevated ca. three to fourfold compared to its parent strain, TF125, in the presence of threonine and isoleucine in culture medium (24). The production level of threonine by TF427 is much higher than TF125. By DNA sequencing of the thr operon, we confirmed the replacement of Ser with Phe at position 345 in the thrA product, AKI-homoserine dehydrogenase I (HDI), of TF5015 (Table 4). Also, the enzyme assay of mutated aspartokinase showed that the activity of AKI was not inhibited by threonine (24). In a previous result about mutant analysis of Serratia marcescens thr operon, threonine-mediated feedback inhibition of both AKI activity and HDI activity was released by single amino acid substitution (the exchange of Gly to Asp at position 330 or Ser to Phe at position 352) in the central region of the AKI-HDI product (29). Similarly, a change at S300Y of Corynebacterium glutamicum AK yielded a threonine- and lysine-insensitive aspartokinase product (21). Genetic analysis revealed that the mutation of thrA results in a release of feedback inhibition of AKI-HDI by threonine in TF5015 rather than a deregulation of feedback repression by threonine plus isoleucine. Meanwhile, amino acid substitution of threonine deaminase (TD; Ser-97
Phe) was detected from the analysis of the nucleotide sequence between W3110 and TF5015 (Table 4). The nutritional requirement of isoleucine and enzyme assay of TD (24) represented that the loss of enzyme activity came from the replacement of Ser by Phe that is located in N2 subdomain of TD (12). Auxotrophic mutants containing the amino acid substitutions (Lys62
Glu, Ala66
Val, or Pro156
Ser) in the N2 subdomain and its nearby loops (12) were reported through genetic experiments (10), which support our interpretation. In this respect, we suggest that the upregulation of thr operon in TF5015 would be caused by limitation of isoleucine due to the inactivation of TD. As shown in Table 5, simple overexpression of the wild-type thr operon resulted in a slight increase in the level of threonine, whereas expression of the mutated thrABC (pATF92) in W3110 gave rise to a significant incremental effect on threonine production. Thus, it is obvious that upregulation of the mutated thr operon in TF5015 led to the overproduction of threonine from the analyses of transcriptome, proteome, and DNA sequences, as well as from the enzyme assay and fermentation experiments. Furthermore, since sequence changes of metL (encodes AKII-HDII) and lysC (encodes AKIII) in contrast to thrA were not detected in ORFs and regulatory regions of both genes (30), we suggest that the biosynthesis of threonine in TF5015 is mainly catalyzed by AKI-HDI rather than by AKII-HDII and AKIII. It is necessary to determine why the expression of thrB was slightly increased in TF5015 in two sets of array experiments. One possibility is that an internal promoter at the 3' end of thrA allows the formation of thrB transcript in addition to those initiated at the major promoter of thrA upstream (30).
Considering the increased accumulation of lysine in the course of fermentation by TF5015 and feedback regulation of the biosynthetic pathway of aspartate family amino acids by lysine (3), blocking of the lysine pathway is likely to be indispensable for the prevention of lysine accumulation, the derepression of asd, and finally an increase in threonine production in TF5015. The increase of threonine production noted in a previous study was accomplished by auxotrophic mutation of
,
-diaminopimelic acid (26). The asd gene, which is essential for threonine biosynthesis, was downregulated in TF5015 (Fig. 3), and this seems to be due to the accumulation of lysine (Table 3) (30).
We observed drastic changes in transcript or protein levels in some genes in response to the accumulation of other amino acids or difference of growth rate. At some promoters of the Lrp regulon, Lrp action was greatly modified by the presence of leucine (27). The expressions of dadAX, leu operon, glnA, glyA, ompF, oppA, oppB, and oppF seem to directly or indirectly respond to an accumulation of leucine, Lrp, or leucine/Lrp (4, 27). Upregulation of aldA is likely to arise from an accumulation of glutamate (Tables 1 and 3 and Fig. 4) (33). The differential expressions of mdoB, proA, and ompF, which encode phosphoglycerol transferase I, a gamma-glutamyl phosphate reductase, and an outer membrane protein F precursor, respectively, might be involved in the response of osmotic regulation of TF5015 (Table 1) (9, 36, 41). As observed in earlier study (38), the downregulation of many ribosomal genes would result from the slow growth rate of TF5015 compared to W3110 (Fig. 2 and Table 1).
The accumulation of acetate in TF5015 was about fourfold lower than in W3110 at the early stationary phase, although acetate was simultaneously accumulated in both W3110 and TF5015 (Table 3), reaching up to 200 and 66.7 mM, respectively, with fermentation time (data not shown). Previous studies revealed that RpoS-regulated genes, periplasmic transporters for amino acids and peptides, and metabolic enzymes are induced either by acetate or at low pH (2, 3, 13, 17, 23). Of these, we suppose that downregulations of b1795, hdeAB operon, oppA, and yfiD are related to a lower accumulation of acetate in TF5015 compared to W3110. In particular, the induction of YfiD can be a strong indication for internal acidification. A low level of acetate and repression of YfiD in TF5015 imply that this mutant maintains intracellular homeostasis at the early stationary phase, even though the overproduction of threonine was expected to affect to some extent the cellular physiology.
We have demonstrated that a global analysis of expression profiles at mRNA and protein levels between prototrophic strain and threonine-producing mutant provides crucial information for understanding the mechanism of threonine overproduction and the physiological consequences in TF5015. Integrated knowledge regarding the threonine-producing mutant is expected to offer more rational strategies for developing microorganisms with greater potential.
This work was supported by the BK21 Program of the Ministry of Education and the National Research Laboratory Program of the Ministry of Science and Technology of Korea.
Present address: R&D Center of Bioproducts, Institute of Science and Technology, CJ Corp., Kyonggi-do, Korea. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»