Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2006, p. 3653-3661, Vol. 72, No. 5
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.5.3653-3661.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Center for Molecular BioEngineering, Driftmier Engineering, University of Georgia, Athens, Georgia 30602,1 Biotechnology Institute,2 Department of Biochemistry, Molecular Biology and Biophysics, University of Minnesota, St. Paul, Minnesota 551083
Received 22 November 2005/ Accepted 21 February 2006
|
|
|---|
|
|
|---|
![]() |
If the rate of oxygen utilization is sufficiently high, the reduced cofactors generated by glucose consumption are reoxidized in the electron transport chain, which serves the dual purpose of maintaining an optimal redox environment and generating energy by oxidative phosphorylation. In the absence of oxygen glucose cannot be completely oxidized, and metabolic intermediates accumulate to maintain the redox balance. Even in the presence of oxygen, if the rate of glucose consumption is greater than the capacity to reoxidize the reduced equivalents generated, the response is similar to what is observed under anaerobic conditions (1, 2, 13). Since the flux from acetyl-CoA to acetate does not generate any NADH while the flux from acetyl-CoA through the TCA cycle generates 8NAD(P)H and 2FADH2, carbon flow diversion to acetate could be viewed as a means to reduce or prevent further NAD(P)H accumulation (7, 13). These inferences regarding acetate overflow have been based on physiological observations and in vitro enzyme assays, and the genetic trigger has not been identified. Details of pathways involved in acetate generation and consumption, including specific enzymes and their regulation, have recently been reviewed (39).
Aided by genomic technology, we further investigated the relationship between redox and acetate overflow by transforming wild-type E. coli with water-forming NADH oxidase (encoded by the nox gene from Streptococcus pneumoniae) in order to oxidize residual NADH (3). This enzyme decouples NADH oxidation (oxygen reduction) from respiratory energy generation (10). We report the impact of this perturbation in redox on the physiological and transcriptional activity of cells under steady-state conditions and correlate the results with acetate formation to begin to identify regulatory processes involved in overflow metabolism. Steady-state chemostat cultures (i.e., as opposed to batch cultures) are an important tool for the study of overflow metabolism and transcriptional changes because (i) a chemostat permits the system to be tuned to achieve a desired glucose consumption rate through the selection of dilution rate, in particular glucose consumption rates both above and below the threshold rate which initiates acetate formation, and (ii) steady-state conditions permit collection of transcriptional information which is itself at a steady state. These results provide insights into gene regulation in E. coli as a model facultative anaerobe during the transition from respiratory to respirofermentative metabolism.
|
|
|---|
arcA::Tet) were used in this study. QC2575 was obtained from D. Touati (l'Institut Jacques Monod, Paris, France). Growth and physiological characteristics were determined using defined media (8) composed of the following (per liter): 5 g glucose, 1.5 g NH4Cl, 0.5 g NaCl, 7.8 g Na2HPO4 · 7H2O, 3.5 g KH2PO4, 0.014 g CaCl2 · 2H2O, 0.246 g MgSO4 · 7H2O, 0.1 ml Antifoam C, 1 mg biotin, 1 mg thiamine, 100 mg ampicillin, and 10 ml trace metal solution. The trace metal solution contained the following (per liter): 16.68 g FeCl3 · 6H2O, 0.36 g ZnSO4 · 7H2O, 0.32 g CuSO4 · 5H2O, 0.2 g MnSO4 · H2O, 0.18 g CoCl2 · 6H2O, 22.4 g EDTA, and 0.1 g NaMoO4 · 2H2O.
Construction of pTrc99A-nox.
The Streptococcus pneumoniae nox gene was amplified by PCR using pPANOX7 (M.-C. Trombe, U. Paul Sabatier, Toulouse, France) as template with Pfu DNA polymerase. Primers were designed based on the published S. pneumoniae nox gene sequence (3) and contained a BamHI restriction site and a Shine-Dalgarno sequence at the beginning of the amplified fragment and a PstI restriction site at the end of the amplified fragment (forward primer, 5'-TAC TAT GGA TCC AGG AGG TAA CAG CTA TGA GTA AAA TCG TTG TAG TCG GTG C-3'; reverse primer, 5'-ATA TAG TGA TCG ATA GCA GTC TGC AGT TAT TTT TCA GCC GTA AGG GCA GC-3' [the underlined sequences are the BamHI, Shine-Dalgarno, ATG start, and PstI sites, respectively]). The resulting 1.4-kb PCR product was gel isolated, digested with BamHI and PstI, and ligated into the pTrc99A expression vector which had been digested with the same two restriction enzymes.
Chemostat cultivation.
Carbon-limited chemostat cultures of 1.5-liter working volume were grown in 2.5-liter vessels (Bioflo II; New Brunswick Scientific, NJ) at 37°C, pH 7.0, and an agitation of 500 rpm. The airflow rate was maintained at 1.5 liter/min using mass flow controllers (Unit Instruments, Orange, CA) to ensure that the dissolved oxygen concentration remained above 40% of saturation at all growth rates studied. Measurements were made after the cells attained a steady state, which required at least 7 volume changes without any perturbation. The biomass formed was quantified by washing the cells with phosphate-buffered saline (pH 7.0) and drying for 12 h at 60°C. Glucose and organic acids in the feed and effluent were measured by high-performance liquid chromatography with a detection limit of about 0.05 g/liter (6). Oxygen uptake rate and CO2 evolution rate were calculated by measuring the effluent concentrations of oxygen and CO2 (Ultramat 23 gas analyzer; Siemens, Germany). Each steady-state growth culture was freshly started from a single colony.
Quantification of NADH/NAD and glycolytic metabolites.
Metabolism was rapidly interrupted by extracting two 10-ml aliquots from a chemostat and plunging them into 40 ml methanol that had been prechilled for 4 h in a dry ice-ethanol bath. The cell pellets were resuspended in 0.2 M HCl (for extracting NAD) or 0.2 M NaOH (for extracting NADH), and the nucleotides were extracted by boiling the cell suspension. A cycling assay (4) which involves the transfer of reducing equivalents from NADH ultimately to 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) was used to measure the specific nucleotides (23). The rate of reduction of MTT as measured at 570 nm was proportional to the concentration of NADH or NAD.
The intracellular concentrations of the key metabolites glucose-6-phosphate, fructose-6-phosphate, phosphoenolpyruvate (PEP), pyruvate, and acetyl-CoA were measured enzymatically by using cell extracts prepared by the perchloric acid method (36).
Global transcription profiling.
Changes in the expression of genes at various growth rates were identified using parallel two-color hybridization to whole-genome E. coli MG1655 spotted DNA arrays corresponding to 98.8% of the annotated open reading frames. The design, printing, and probing were previously described in detail (19, 20). After attaining a steady state at a predetermined dilution rate, samples were extracted from the chemostat and placed in RNAprotect buffer (QIAGEN, Valencia, CA), and the cell pellets were frozen at 80°C. Total RNA was extracted by the hot phenol-chloroform method and treated with DNase I in the presence of RNase inhibitor for subsequent labeling by reverse transcription with Cy3-dUTP and Cy5-dUTP fluorescent dyes. Total RNA from the strain containing pTrc99A plasmid was cultured at a dilution rate of 0.1 h1 and used as the common reference (always labeled with Cy3-dUTP) against which total RNA extracted from cells cultured at five higher equally spaced dilution rates (always labeled with Cy5-dUTP) was hybridized. Similarly, for the strain containing the pTrc99A-nox plasmid, total RNA extracted from cells cultured at a dilution rate of 0.06 h1 was the reference against which total RNA from cells cultured at five higher equally spaced dilution rates was hybridized. Differential gene expression between the two references (the NOX strain cultured at a dilution rate of 0.1 h1 and labeled with Cy3-dUTP, and the NOX+ strain cultured at a dilution rate of 0.06 h1 and labeled with Cy5-dUTP) was also measured to identify transcriptional changes due to the presence of the S. pneumoniae NADH oxidase. All the hybridizations were performed at least in triplicate by using biologically independent samples as described previously (20) and incubating the labeled mixture on the arrays at 65°C overnight. The slide was subsequently washed and scanned using a GenePix 4000B microarray scanner (Axon Instruments, Union City, CA). The degree of labeling of the two dyes was quantified by measuring the intensity at a wavelength of 532 nm (for Cy3) or 635 nm (for Cy5). The relative expression of a gene was calculated as the base 2 logarithmic ratio of the background subtracted intensity from the Cy5 channel to the background subtracted intensity from the Cy3 channel, and the resulting value was referred to as the expression ratio.
Analysis of gene expression data.
Our goal for the analysis of transcription data was to identify genetic changes that corresponded to physiological observations. Specifically, we were interested in identifying those genes whose expression was sensitive to a perturbation in the redox ratio (i.e., the NADH/NAD ratio). We calculated the Pearson correlation coefficient for the expression ratio of each gene with the redox ratio for NOX+ and NOX as a function of specific glucose consumption rate. Only those genes whose expression ratios had a high correlation coefficient (R > 0.9 or R < 0.9) with the redox ratio were considered for further analysis. These highly correlated (or anticorrelated) genes were classified into 22 functional categories according to the method of Riley (35). Each functional category was tested for significant overrepresentation (P < 0.05) by using a hypergeometric distribution (15). With a priori information on the distribution of the global gene set among the 22 categories, hypergeometric distribution measures the enrichment of a functional category based on the number of genes of that particular category appearing in the cluster. The P value for each category was calculated according to the following equation:
![]() |
Only those genes from the significantly enriched functional categories were selected to study common regulatory mechanisms governing their expression. Any coregulation among these coexpressed genes was identified by searching for common transcription factor binding sites upstream of their transcription start sites. Sequences 300 bp upstream of the filtered genes were analyzed for common sequence motifs by using the hidden Markov model-based BioProspector software (24).
|
|
|---|
As shown in Fig. 1, acetate overflow is directly related to the rate at which the sole carbon source (glucose) is consumed, with acetate formation occurring only after glucose consumption surpasses some threshold rate. The presence of heterologous NADH oxidase had the effect of increasing the critical glucose consumption rate (qScrit) at which acetate first appeared and thereby delaying the entry of E. coli into respirofermentative overflow metabolism (Fig. 1). This transition between respiratory and respirofermentative metabolism occurred at a qScrit of 0.8 g/g dry cell weight (DCW) h for NOX and 1.2 g/g DCW h for NOX+. The expression of NADH oxidase therefore increased by 50% the value of qScrit. During respirofermentative metabolism, NOX+ exhibited a lower effluent acetate concentration and a lower specific acetate formation rate (qA) than NOX at any given qS. Biomass yield (YX/S) from glucose (g dry cell weight/g glucose consumed) was 0.42 to 0.48 g/g for NOX during respiratory metabolism but decreased during respirofermentative metabolism, consistent with a portion of the glucose carbon being diverted from biomass synthesis to acetate formation. For NOX+, YX/S remained 0.28 g/g at glucose consumption rates above 0.5 g/g DCW h (Fig. 1).
![]() View larger version (16K): [in a new window] |
FIG. 1. Steady-state physiological profiles of E. coli in the presence of heterologous NADH oxidase. YX/S ( , ) and qA ( , ) values are compared for NOX (open symbols and dashed lines) and NOX+ (solid symbols and lines) as functions of the specific glucose consumption rate. The highest dilution rate studied was about 80% of µmax for both strains. The arrows indicate for each strain the critical specific glucose consumption rates at which acetate formation commenced.
|
![]() View larger version (14K): [in a new window] |
FIG. 2. Steady-state respiration for NOX (open symbols and dashed lines) and NOX+ (solid symbols and lines). The steady-state qO2 ( , ) and qCO2 ( , ) values are shown as functions of qS.
|
![]() View larger version (16K): [in a new window] |
FIG. 3. In vivo molar concentration ratio of NADH/NAD for NOX ( ) and NOX+ ( ) as functions of qS. The critical value of the NADH/NAD ratio at which acetate formation commences is about 0.06 for both NOX and NOX+ (indicated by vertical lines). qA values are also shown for NOX ( ) and NOX+ () as functions of qS.
|
![]() View larger version (21K): [in a new window] |
FIG. 4. Intracellular concentrations of key glycolysis metabolites glucose-6-phosphate (red), fructose-6-phosphate (green), PEP (yellow), pyruvate (blue), and acetyl-CoA (pink) were measured under steady-state conditions in NOX (A) and NOX+ (B). qA values are also shown for NOX (open circles and dashed lines) and NOX+ (filled circles and solid lines) as functions of qS.
|
Transcriptional response to increasing glucose consumption rate.
Since most physiological events originate at the transcription level, we measured the transcriptional responses to changes in qS for NOX and NOX+ strains to establish a genetic basis for the observed physiological changes. For each strain, we used a low value for qS (corresponding to dilution rates of 0.1 h1 for NOX and 0.06 h1 for NOX+) as the reference. We first compared the transcription profiles for the two reference cultures (NOX grown at 0.1 h1 and NOX+ grown at 0.06 h1) to identify transcriptional changes only due to the presence of NADH oxidase. There were no significant transcriptional changes between these two reference cultures, suggesting limited influence of NADH oxidase at low qS. This result is consistent with the similar physiological parameters (YX/S, qS, qCO2, and qO2) and redox ratios observed for the two strains at low qS values (Fig. 1 to 3).
Next, we compared the transcriptional changes at higher values for qS relative to the appropriate reference culture for each strain. In general, we did not observe drastic changes in gene expression, but many genes exhibited a reproducible monotonically increasing or decreasing behavior relative to the reference as qS or growth rate increased. For NOX, the expression of 427 genes varied significantly with qS (P < 0.01), while only 47 genes achieved this level of significance for NOX+ and only 21 genes were common to both subsets. Among the genes whose expression varied significantly for both NOX and NOX+ were key genes in the biosynthesis of threonine, serine, and nucleotides along with acs, which encodes acetyl-CoA synthetase. Since expression of these genes changed with qS for both strains, their expression is presumably largely glucose consumption rate dependent and relatively insensitive to the redox state of the cell. The average expression ratios of all genes involved in the central metabolic pathways for NOX and NOX+ relative to their respective reference cultures are shown in Fig. 5 as functions of qS. Transcription profiles of individual genes in these central metabolic pathways for the two strains are shown in Fig. S1 in the supplemental material.
![]() View larger version (21K): [in a new window] |
FIG. 5. Transcriptional profile of central metabolic pathways for NOX (dashed lines) and NOX+ (solid lines). The mean values of the expression ratios are shown for all genes involved in glycolysis (green), the TCA cycle (red), the pentose phosphate pathway (blue), and respiration (black) as functions of qS. Vertical lines show the demarcation between respiratory and respirofermentative metabolism for NOX (dotted) and for NOX+ (solid). See Fig. S1 in the supplemental material for detailed expression profiles of individual genes.
|
The relative expression of intermediate metabolic genes involved in the biosynthesis of amino acids and nucleotides increased with qS, before either stabilizing or slightly decreasing at high qS values for both NOX and NOX+ (see Fig. S1 in the supplemental material). Moreover, we did not observe significant differences in the expression profiles of these genes between the strains at any given value of qS, except for those involved in methionine and glycine biosynthesis. The relative expression of metABCEHJL genes either remained steady or decreased with qS for NOX, while these genes were significantly upregulated for NOX+. The glycine biosynthesis genes, particularly glyA, were repressed for NOX, but the repression appeared to be reduced for NOX+. Purine and pyrimidine nucleotide biosynthesis genes monotonically increased with qS for both strains. We also observed a repression in most of the transport genes at high qS values for both NOX and NOX+. Among the genes encoding for symport or antiport proteins we found a few genes of the multifacilitator family that were upregulated (such as yhfC, yhaU, codB, uraA, and proP) while all others (including yjcG, lacY, gltS, gntT, dctA, tatC, melB, and nupG) were repressed (see Fig. S1 in the supplemental material). Genes belonging to the ATP binding cassette transporters and the PEP-dependent phosphotransferase systems for the uptake of several sugars (including glucose) were also strongly repressed for both NOX and NOX+ as qS increased (see Fig. S1 in the supplemental material). There was no particular trend observed in most unclassified genes except the gene encoding b4249 (a putative oxidoreductase), which was induced for NOX but was repressed for NOX+ as qS increased.
Identifying coregulation among coexpressed genes.
The approach governing our data analysis methodology was to first group coexpressed genes and then evaluate these gene groups for common regulatory mechanisms (see Materials and Methods). Genes involved in the biosynthesis of amino acids, cofactors, macromolecules, and nucleotides along with central and intermediary metabolic genes were positively correlated with qS (R > 0.9) for both NOX and NOX+. Among the genes that were negatively correlated with qS (R < 0.9) were those responsible for the degradation of small molecules, transport proteins, and unclassified genes. Interestingly, we observed that only the genes involved in the biosynthesis of amino acids, cofactors, and nucleotides along with the central and intermediary metabolic genes were correlated (R > 0.9) with the redox ratio for NOX (Fig. 6). There were also several partially classified genes in this subset, suggesting that the expression of a majority of these genes depends on the rate of glucose consumption and/or redox. While we identified strongly overrepresented sequences upstream of genes correlated with qS, we could not relate these sequences with any of the known promoter binding sites. However, a significantly overrepresented (P < 10170) sequence (Fig. 6) upstream of the genes correlated with the redox ratio for NOX was identified (by BioProspector) as the binding site for ArcA (25, 29, 33). The identification of an ArcA binding site upstream of genes that were correlated with the redox ratio (for NOX) is consistent with a recent discovery that cellular redox state is the signal for the activation of ArcB signal transduction (9, 28). Table S2 in the supplemental material provides a complete list of genes in NOX that showed a reduction in expression by high NADH/NAD (negatively correlated with redox ratio) and which were determined (by BioProspector) to have a binding site for ArcA. In light of these results, we speculated that the strong repression observed for several TCA cycle and respiratory genes in NOX at high qS values might be relieved by deleting arcA.
![]() View larger version (29K): [in a new window] |
FIG. 6. (Left) Hierarchical clustering of genes (rows) that are correlated (R > 0.9 or R < 0.9) with the redox ratio (NADH/NAD) in NOX as a function of increasing qS (columns). (Right) Significantly overrepresented functional categories are shown in the table, along with the number of genes in each category and the P value of its significance as calculated using a hypergeometric distribution. Several key genes involved in the TCA cycle, respiration, and biosynthesis exhibited a strong negative correlation with the redox ratio. A large portion of the genes negatively correlated to the redox ratio were partially classified, revealing redox-dependent regulation of many of these genes.
|
![]() View larger version (14K): [in a new window] |
FIG. 7. Physiological characterization of ARCA NOX (open symbols and dashed lines) and ARCA NOX+ (solid symbols and lines) in accelerostat cultures. YX/S ( , ) and qA ( , ) values are compared as functions of specific glucose consumption rate. The steady-state values of these parameters, obtained for NOX (dashed lines without symbols) and NOX+ (solid lines without symbols) by using chemostats, are also shown.
|
![]() View larger version (15K): [in a new window] |
FIG. 8. Respiration of ARCA NOX (open symbols and dashed lines) and ARCA NOX+ (solid symbols and lines) in accelerostat cultures. qO2 ( , ) and qCO2 ( , ) values are compared as functions of specific glucose consumption rate. The steady-state values of these parameters, obtained for NOX (dashed lines without symbols) and NOX+ (solid lines without symbols) by using chemostats, are also shown.
|
|
|
|---|
Although the rapid generation and subsequent oxidation of NADH in the NOX+ strain essentially introduces a futile NAD turnover, it reveals two very important metabolic events correlated with overflow metabolism in E. coli: both the redox ratio and the pyruvate/PEP ratio are correlated with the appearance of acetate. First, the redox ratio at the onset of acetate overflow was, surprisingly, identical for NOX and NOX+ (Fig. 3), indicating a relationship between the redox state of the cell and overflow metabolism. Since numerous reactions utilize or generate NADH, the redox ratio and overflow metabolism are likely to be the consequences of a complex network of metabolic events. The importance of NADH/NAD in by-product formation in E. coli has been previously demonstrated through increased reduction of NAD, which resulted not only in increased acetate but also in the appearance under aerobic conditions of typical fermentation products (5). In our study, TCA cycle genes which were generally repressed for NOX with increasing qS values commonly showed less repression upon introduction of NADH oxidase. Considering that acetate overflow has thus far been assumed to be due to rate-limiting enzymes of the TCA cycle or the electron transport chain attaining maximum reaction velocity (1, 7, 13, 27), our results provide evidence that acetate overflow occurs as a consequence of transcriptional repression of the TCA cycle and respiratory genes (see Fig. S1 in the supplemental material). The introduction of NADH oxidase appears to delay the attainment of the critical redox ratio and limit acetate formation.
Also, the pyruvate/PEP ratio appears to be related to acetate formation (Fig. 4). As pyruvate is the branch point between respiration and fermentation and a precursor to several macromolecules, its level is highly regulated. In E. coli, PEP is a cosubstrate for glucose uptake and for the principal anaplerotic pathway during growth on glucose. Since acetate is produced from pyruvate directly (via pyruvate oxidase) or indirectly (via the PDH complex and acetate pathway), a 25-fold increase in the pyruvate/PEP ratio would shift the thermodynamic equilibrium towards pyruvate utilization by these pathways. Although an increase in the pyruvate/PEP ratio may not directly cause acetate overflow, the observed shift in the control strain NOX does signal the onset of a bottleneck at the entrance to the TCA cycle, which we have shown can be modulated by redox. The introduction of NADH oxidase served to decrease the pyruvate/PEP ratio and by mass action would make acetate formation less favorable. These results provide circumstantial evidence for considering pyruvate to be one of the candidate signaling metabolites for inducing the phosphorylation of ArcB (14). NADH was previously proposed as a possible signal (14), and more recent evidence (9) indicates that the cellular redox state is the signal for the activation of Arc regulation while pyruvate is an allosteric activator (9).
The strong correlation (R > 0.9) between the redox ratio and the expression of genes involved in central and intermediary metabolism and the biosynthesis of amino acids, cofactors, and nucleotides demonstrates the important regulatory control exerted by the redox state. The identification of binding sites for the ArcA protein upstream of many of these genes suggests redox-dependent regulation of the ArcAB system and is consistent with recent studies which propose that the redox state triggers the Arc system (9). Our analysis does not rule out the possibility of secondary regulation, and therefore the relationships between redox state, ArcA, and acetate overflow could be indirect. Although ArcA-mediated repression has been reported for many individual genes (operons), their integrated effect on induction of overflow metabolism has been largely overlooked. The reduced redox ratio for NOX+ may delay significant activation of the Arc system. For the ARCA strains, several of the TCA cycle and respiratory genes were induced and qCO2 was elevated, demonstrating greater TCA cycle flux. This view is in accordance with recent 13C-labeling studies using Arc mutants of E. coli and showing increased TCA flux under both aerobic and anaerobic conditions (34). The resulting higher rate of NADH formation appears to be accommodated at least partly by the elevated qO2 which results from derepression of the respiratory chain in these strains. Importantly, although the qScrit was about 10% greater for ARCA NOX compared to NOX, acetate formation was even more pronounced for ARCA NOX at higher levels of qS. One possible explanation for this observation is that the heightened TCA cycle flux resulting from the absence of ArcA-mediated repression elevated NADH accumulation to a level beyond the capacity of the (derepressed) respiratory chain. Without a transcriptional mechanism to prevent further NADH formation in the TCA cycle, acetate formation may have occurred through some other mechanism (such as inhibition of citrate synthase). The overexpression of NADH oxidase in the arcA strain seems sufficient to provide another outlet for NADH oxidation and prevent acetate formation even at high glucose consumption rates.
In summary, our results using steady-state chemostats support a model in which an increase in the redox ratio contributes to a repression of the TCA cycle and to acetate formation, and they suggest that this overflow is due to transcriptional limitation. Providing another outlet for NADH turnover relieves TCA cycle gene repression and delays acetate formation. An arcA mutation delays the onset of acetate formation (only so far) through the maintenance of both TCA cycle flux and respiration. Moreover, a strain with an arcA mutation and heightened NADH oxidase activity appears able to both maintain an elevated TCA cycle flux and alleviate NADH buildup, thereby preventing acetate formation altogether. Considering the deleterious impact of acetate on growth (26) and recombinant protein production (37) and the wide variety of genetic and process approaches proposed to reduce acetate formation, our findings provide evidence at the level of transcription for the cause of acetate overflow as well as offer a means to overcome it.
This project was supported by the U.S. DOE (grant number DE-FG36-01ID14007), which provided a fellowship to G.N.V., and the Georgia Experiment Stations. D.P.S. was supported in part by the NSF (grant number QSB 0222636 to Friedrich Srienc and A.B.K.).
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
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»