Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2007, p. 4881-4891, Vol. 73, No. 15
0099-2240/07/$08.00+0 doi:10.1128/AEM.00177-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Biotechnology, Delft University of Technology, Julianalaan 67, 2628 BC Delft,1 Bird Engineering B.V., Westfrankelandsedijk 1, 3115 HG Schiedam, The Netherlands2
Received 24 January 2007/ Accepted 26 May 2007
|
|
|---|
|
|
|---|
Saccharomyces cerevisiae is presently the organism of choice for industrial ethanol production. Although wild-type S. cerevisiae strains rapidly ferment hexoses with high efficiency, they cannot grow on or use pentoses, such as D-xylose and L-arabinose (3). In addition to the development of pentose-consuming bacteria such as Zymomonas mobilis, Escherichia coli, and Klebsiella oxytoca as alternative biocatalysts for ethanol production (5), this situation has inspired various studies to expand the substrate range of S. cerevisiae. The combination of metabolic and evolutionary engineering with the heterologous expression of either yeast xylose reductase and xylitol dehydrogenase (14, 32, 34, 35, 41) or a fungal xylose isomerase (19-22) has already enabled the anaerobic fermentation of D-xylose by S. cerevisiae. The next challenge is the fermentation of other pentoses, such as L-arabinose. Although several yeasts and fungi can utilize L-arabinose as a carbon and energy source, most of them are unable to ferment it into ethanol. At best, very low ethanol yields and production rates are accompanied by the formation of L-arabinitol under microaerophilic or oxygen-limited conditions (6, 18, 26). The rarity of ethanolic arabinose fermentation may be due to a redox imbalance in the fungal arabinose pathway as a result of the use of NADPH for the reductive reactions and the production of NADH in the oxidation reactions (6, 36) (Fig. 1). Consistent with this, the overexpression of all the structural genes of the fungal L-arabinose pathway (XYL1, lad1, lxr1, XYL2, and XKS1) in S. cerevisiae does not result in fast and efficient fermentation of L-arabinose into ethanol (29). Although the engineered strain produced only 0.35 mg of ethanol g–1 h–1 under anaerobic conditions, it provided the first example of ethanolic arabinose fermentation by S. cerevisiae.
![]() View larger version (26K): [in a new window] |
FIG. 1. Schematic representation of D-xylose and L-arabinose catabolism in metabolically engineered S. cerevisiae strains described in the literature. The engineering steps involved in this work are indicated by the underlined gene names. Components of the catabolism are as follows: aldose/xylose reductase (GRE3/XYL1), xylitol dehydrogenase (XYL2), xylulokinase (XKS1), D-xylose isomerase (xylA), arabinitol 4-dehydrogenase (lad1), L-xylulose reductase (lxr1), L-arabinose isomerase (araA), L-ribulokinase (araB), L-ribulose-5-phosphate 4-epimerase (araD), transaldolase (TAL1), transketolase (TKL1), D-ribulose-5-phosphate 3-epimerase (RPE1), and ribose-5-phosphate ketol-isomerase (RKI1). PEP, phosphoenolpyruvate.
|
The present study aimed for fast and efficient fermentation of arabinose by S. cerevisiae under anaerobic conditions. To this end, a genetically engineered S. cerevisiae strain was constructed and subjected to extensive evolutionary engineering for improved anaerobic arabinose utilization.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. S. cerevisiae strains and plasmids used in this study
|
In order to get a high level of expression, the L. plantarum araA and araD genes were ligated into the xylA-bearing plasmid pAKX002. The araA expression cassette was constructed by amplifying the S. cerevisiae TDH3 promoter (PTDH3) with oligonucleotides SpeI5'Ptdh3 and 5'AraAPtdh3, the araA gene with oligonucleotides Ptdh5'AraA and Tadh3'AraA, and the ADH1 terminator (TADH1) with oligonucleotides 3'AraATadh1 and 3'Tadh1SpeI. The three fragments were extracted from gel and mixed in roughly equimolar amounts. With this mixture, a PCR using the SpeI5'Ptdh3 and 3'Tadh1SpeI oligonucleotides was performed. The resulting PTDH3-araA-TADH1 cassette was gel purified, cut at the 5' and 3' SpeI sites, and then ligated into NheI-digested pAKX002, resulting in plasmid pRW230.
The araD construct was made by first amplifying a truncated version of the HXT7 promoter (PHXT7) with oligonucleotides SalI5'Phxt7 and 5'AraDPhxt, the araD gene with Phxt5'AraD and Tpgi3'AraD, and the PGI1 terminator region (TPGI) with the 3'AraDTpgi and 3'TpgiSalI oligonucleotides. The resulting fragments were extracted from gel and mixed in roughly equimolar amounts, after which a PCR using oligonucleotides SalI5'Phxt7 and 3'Tpgi1SalI was performed. The resulting PHXT7-araD-TPGI1 cassette was gel purified, cut at the 5' and 3' SalI sites, and then ligated into XhoI-digested pRW230, resulting in plasmid pRW231 (Fig. 2).
![]() View larger version (13K): [in a new window] |
FIG. 2. Maps of plasmids pRW231 and pRW243 containing the genes encoding L-arabinose isomerase (araA), L-ribulokinase (araB), and L-ribulose-5-phosphate epimerase (araD) from L. plantarum, the Piromyces sp. strain E2 xylose isomerase gene (xylA), and the yeast xylulokinase gene (XKS1).
|
Plasmids pRW243 and pRW231 were used to transform S. cerevisiae strain RWB220, a strain that is derived from strain RWB219. Strain RWB219 was obtained by the selection of RWB217 colonies for improved growth on D-xylose, similar to the way in which strain RWB218 was developed (21), with the difference that RWB218 was obtained by plating and restreaking onto MYD plates while RWB219 was selected from solid synthetic medium supplemented with D-xylose. Strain RWB219 was grown nonselectively on YPD in order to facilitate the loss of plasmids pAKX002 and p415ADHXKS1 (20), harboring the URA3 and LEU2 selective markers, respectively. After plating onto YPD, single colonies were screened for plasmid loss by testing for uracil and leucine auxotrophy. In order to remove a KANMX cassette, still present after the integration of the RKI1 overexpression construct (20), a strain that had lost both plasmids was transformed with pSH47, containing the Cre recombinase gene (10). Transformants containing pSH47 were resuspended in complex medium with 1% D-galactose and incubated for 1 h at 30°C. Cells were plated onto YPD, and colonies were screened for the loss of the KANMX marker (G418 resistance) and pSH47 (URA3). A strain that had lost both the KANMX marker and the pSH47 plasmid was designated RWB220. Strain RWB220 was transformed with pRW231 and pRW243 (Table 1; Fig. 2), resulting in strain IMS0001. The selection of strain IMS0001 for (anaerobic) growth on L-arabinose, as described in Results, yielded strain IMS0002.
Cultivation.
Shake flask and fermenter cultivations were performed at 30°C in MY containing 5 g liter–1 (NH4)2SO4, 3 g liter–1 KH2PO4, 0.5 g liter–1 MgSO4·7H2O, 0.05 ml liter–1 silicon antifoam, and trace elements (39). For the cultivation in shake flasks, the pH of the medium was adjusted to 6.0 with 2 M KOH prior to sterilization. After heat sterilization (121°C for 20 min), a filter-sterilized vitamin solution (39) and an appropriate carbon and energy source were added. Shake flask cultures and precultures for anaerobic batch cultivations were prepared by inoculating 100 ml of medium containing the appropriate sugar into a 500-ml shake flask with a frozen stock culture. After incubation at 30°C in an orbital shaker (200 rpm), precultures were used to inoculate either another 500-ml shake flask containing MY with an appropriate carbon source or an anaerobic fermenter. For anaerobic cultivation in fermenters, the synthetic medium was supplemented with 0.01 g liter–1 ergosterol and 0.42 g liter–1 Tween 80 dissolved in ethanol (1, 2). Anaerobic (sequencing) batch cultivation was carried out at 30°C in 2-liter laboratory fermenters (Applikon, Schiedam, The Netherlands) with a working volume of 1 liter. The culture pH was maintained at 5.0 by the automatic addition of 2 M KOH. Cultures were stirred at 800 rpm and sparged with 0.5 liter min–1 nitrogen gas (<10 ppm oxygen). To minimize the diffusion of oxygen, fermenters were equipped with Norprene tubing (Cole-Palmer Instrument Company, Vernon Hills, IL). Dissolved oxygen was monitored with an oxygen electrode (Applisens, Schiedam, The Netherlands). Oxygen-limited conditions in the same experimental setup were achieved by headspace aeration at approximately 0.05 liter min–1.
Determination of biomass dry weight.
Culture samples (10.0 ml) were filtered over preweighed nitrocellulose filters (pore size, 0.45 µm; Gelman Laboratory, Ann Arbor, MI). After filtration of the broth, the biomass was washed with demineralized water, dried in a microwave oven for 20 min at 360 W, and weighed. Duplicate determinations varied by less than 1%.
Gas analysis.
Exhaust gas was cooled in a condenser (2°C) and dried with a type MD-110-48P-4 dryer (Perma Pure, Toms River, NJ). Oxygen and carbon dioxide concentrations were determined with an NGA 2000 analyzer (Rosemount Analytical, Orrville, OH). The exhaust gas flow rate and specific carbon dioxide production rates were determined as described previously (37, 42). In calculating the cumulative carbon dioxide production, volume changes caused by withdrawing culture samples were taken into account.
Metabolite analysis.
Glucose, xylose, arabinose, xylitol, organic acids, glycerol, and ethanol were analyzed by high-performance liquid chromatography (HPLC) using an Alliance 2690 HPLC system (Waters, Milford, MA) supplied with an HPX 87H column (Bio-Rad, Hercules, CA), a Waters 2410 refractive-index detector, and a Waters 2487 UV detector. The column was eluted at 60°C with 0.5 g liter–1 sulfuric acid at a flow rate of 0.6 ml min–1.
Rate calculations.
For calculation of the specific rates of arabinose consumption and ethanol production, the time-dependent arabinose and ethanol data were fitted with Boltzmann sigmoidal equations. For each time point, the specific arabinose consumption rate and the ethanol production rate were calculated by dividing the derivative of the fitted curve by the dry weight.
Carbon recovery.
Carbon recoveries were calculated as the amount of carbon in the products formed divided by the total amount of sugar carbon consumed and were based on a carbon content of the biomass of 48% (38). To correct for ethanol evaporation during the fermentations, the amount of ethanol produced was assumed to be equal to the measured cumulative amount of CO2 produced minus the amounts of the CO2 produced due to biomass synthesis (5.85 mmol of CO2 per g of biomass) (38) and the CO2 associated with acetate formation.
Total RNA extraction and cDNA synthesis.
Strains RWB219, IMS0001, and IMS0002 were cultivated in shake flasks containing MYD. At an optical density at 660 nm (OD660) of approximately 2, when glucose was not yet depleted, cells were harvested by centrifuging 40 ml of the culture at 0°C. Total RNA extraction was performed as described previously (28). Isolated RNA was DNase I treated according to the RNeasy mini kit (QIAGEN) protocol for RNA cleanup. First-strand cDNA synthesis was carried out using 2 µg of DNase-treated RNA, random hexamer primers, and Moloney murine leukemia virus reverse transcriptase (Invitrogen, Carlsbad, CA).
Expression profiling using real-time quantitative PCR.
Quantitative real-time PCR was run on a DNA engine Opticon I system (Bio-Rad, Hercules, CA) with the following settings: 94°C for 2 min, 94°C for 10 s, 66°C for 10 s, 72°C for 10 s, and plate reading. The denaturation, annealing, elongation, and reading steps were repeated for 39 cycles. A melting curve from 50 to 94°C was generated at the end of the reaction. The 20-µl reaction mixture consisted of 10 µl of SYBR green Taq ready mix (catalog no. S1816; Sigma), 0.2 µM forward primer, 0.2 µM reverse primer, and cDNA (20, 2, 0.2, or 0.02 ng per PCR). The primers used are listed in Table 2. The cycle threshold value was calculated with the Opticon Monitor software version 1.08 (Bio-Rad, Hercules, CA) by setting the threshold for significant detection levels to 10 times the standard deviation over the cycle range from 3 to 7. Each PCR was carried out in triplicate. The relative expression of xylA, XKS1, araA, araB, and araD was quantified by using the comparative 2–
CT method (24) with ACT1 as the reference gene.
|
View this table: [in a new window] |
TABLE 2. Oligonucleotides used in this work
|
|
|
|---|
Since several relevant substrates for enzymes encoded by araA, araB, and araD are not commercially available (16), thus precluding enzyme activity determination, the expression of these genes in strain IMS0001 at the mRNA level was analyzed by real-time quantitative PCR. Figure 3 displays the levels of expression of the genes relative to the expression of ACT1. This highly expressed housekeeping gene encoding actin is commonly used as an internal standard for quantitative mRNA analysis (24). The highest levels of expression were found to be those of xylA and araA, encoding the D-xylose isomerase and L-arabinose isomerase, respectively. Both genes were transcribed at levels approximately 8.5-fold higher than the reference gene ACT1. However, the levels of expression of araB, araD, and XKS1 were lower than that of ACT1, at 0.35-, 0.02-, and 0.04-fold the level of ACT1 expression, respectively. Although these quantitative PCR results confirm the transcription of araA, araB, and araD, strain IMS0001 was not capable of growing on solid MYA.
![]() View larger version (22K): [in a new window] |
FIG. 3. Expression of xylA, XKS1, araA, araB, and araD in strains RWB219 (black bars), IMS0001 (gray bars), and IMS0002 (white bars). The levels of expression displayed are relative to the level of expression of the reference gene ACT1 and are the averages calculated from cycle threshold values measured in triplicate reactions with mixtures containing 0.02, 0.2, 2, and 20 ng of cDNA. The error bars indicate the standard deviations calculated by using standard propagation of error methods.
|
![]() View larger version (18K): [in a new window] |
FIG. 4. (A) Specific growth rate (µ) of S. cerevisiae IMS0001 in shake flask cultures containing synthetic medium with 2% (wt/vol) L-arabinose during serial transfers. Each point represents the growth rate estimated from the OD660 measured during (exponential) growth. The closed and open circles represent results from duplicate serial transfer experiments. (B) Growth rate and CO2 profile during an anaerobic SBR fermentation of S. cerevisiae IMS0001 in synthetic medium with 2% (wt/vol) L-arabinose. Each point (closed circle) represents the growth rate estimated from the CO2 profile (lines) during exponential growth.
|
mol of ATP per mol of arabinose that is fermented into ethanol, which is much lower than the 13
mol of ATP per mol of arabinose produced under respiratory conditions, assuming a ratio of ATP produced per oxygen atom of 1.0 for all redox equivalents in S. cerevisiae (40). Therefore, under anaerobic (fermentative) conditions, an 8-fold-higher arabinose influx is required to achieve a specific ATP production rate equal to that achieved under aerobic (respiratory) conditions. In batch cultures, cells are growing at the maximum growth rate, and the sequential transfer of these batch cultures should select for cells with an increasingly higher specific growth rate and thus a higher arabinose consumption rate. In contrast, chemostat cultures usually select for cells with a lower saturation constant (Ks) for the specific substrate (27). To automate the sequential transfers of the culture, a sequencing batch reactor (SBR) setup was used (21, 22). The first cycle of anaerobic batch fermentation was initiated by inoculating synthetic medium supplemented with 2% L-arabinose with aerobically evolved L-arabinose-grown S. cerevisiae cells. During the aerobic serial transfer protocol, the maximum aerobic growth rate of these cells on arabinose had increased to approximately 0.12 h–1. Initially, growth under anaerobic conditions was not observed. To allow for a more gradual transfer to anaerobic conditions, the first cycle of SBR fermentation was performed under oxygen-limited conditions. When growth was observed, the culture was switched to anaerobic conditions in the next batch cycle. Upon arabinose depletion, indicated by the decrease of the CO2 percentage to below 0.05% after the CO2 production peak (Fig. 4B), a new cycle was initiated by either manual or automated replacement of approximately 90% of the culture with fresh synthetic medium containing 20 g liter–1 L-arabinose. For each cycle, the maximum specific growth rate was estimated from the CO2 profile (Fig. 4B). In 13 cycles, the anaerobic specific growth rate increased from 0.025 to 0.08 h–1. In the next seven cycles, the growth rate did not increase noticeably. After a total of 20 cycles, single colonies were isolated on solid MYA. A culture originating from a single colony grown in a shake flask containing MYA was designated strain IMS0002 and was used for further characterization.
The analysis of the expression of araA, araB, araD, xylA, and XKS1 in strain IMS0002 by real-time quantitative PCR revealed that all five genes were expressed at higher levels than they were in the parental strain IMS0001 (Fig. 3). In particular, araB and araD were expressed at much higher levels: 52- and 90-fold higher than in strain IMS0001. The levels of expression of the genes displaying the highest relative expression levels, araA and xylA, increased only six- and twofold, respectively. Surprisingly, the expression of XKS1, encoding xylulokinase, also increased by 34-fold in strain IMS0002.
Characterization of an S. cerevisiae strain engineered and selected for anaerobic growth on L-arabinose.
To characterize the growth of and product formation by the evolved strain IMS0002, anaerobic batch cultivations of strain IMS0002 in synthetic medium with 20 g liter–1 of L-arabinose as the sole carbon source were performed. Precultures for these anaerobic batch fermentations were prepared in aerobic shake flasks containing 100 ml of MYA. With L-arabinose as the sole carbon source, strain IMS0002 was capable of growing anaerobically at a specific growth rate of 0.051 ± 0.001 h–1 (Table 3). After 52 h, no further increase of the biomass was observed (Fig. 5C). With an initial biomass of 0.28 g (dry weight) liter–1, strain IMS0002 completely consumed 20 g liter–1 L-arabinose within 70 h (Fig. 5A) with a maximum specific arabinose consumption rate of 0.70 g h–1 g (dry weight)–1 (Table 4). The ethanol production was deduced from the CO2 production to correct for ethanol evaporation during the batch fermentations (22). From 138 ± 1 mmol liter–1 L-arabinose, 194 ± 2 mmol liter–1 ethanol was produced, corresponding to a high ethanol yield of 0.43 g g–1 of arabinose. The maximum specific ethanol production rate was determined to be 0.29 g h–1 g (dry weight)–1 (Table 4). By-products such as glycerol, lactate, and succinate were formed only in small amounts, similar to quantities observed previously in the fermentation of xylose by engineered S. cerevisiae strains (20-22). The formation of arabinitol, which in other studies was found to be produced in substantial amounts by (engineered) yeasts fermenting arabinose (6, 16, 33), was not observed during the anaerobic growth of strain IMS0002 on L-arabinose.
|
View this table: [in a new window] |
TABLE 3. Growth parameters, final substrate and product concentrations, and product yields determined for anaerobic batch cultivations of strain IMS0002a
|
![]() View larger version (31K): [in a new window] |
FIG. 5. Sugar consumption and product formation by strain IMS0002 during anaerobic batch fermentations in synthetic medium supplemented with 20 g liter–1 L-arabinose (A and C) and 20 g liter–1 glucose and 20 g liter–1 L-arabinose (B and D). Data points represent the averages of results from two independent anaerobic batch fermentations. (A and B) Levels of glucose (), L-arabinose ( ), ethanol (calculated from cumulative CO2 production) ( ), glycerol ( ), and cumulative CO2 (—) produced. l, liter. (C and D) Levels of acetate () and lactate ( ) produced and biomass dry weight ( ).
|
|
View this table: [in a new window] |
TABLE 4. Arabinose consumption rates, ethanol production rates, and ethanol and arabinitol production yields for engineered arabinose-utilizing S. cerevisiae strainsa
|
To test the capability for cofermentation of D-glucose, D-xylose, and L-arabinose, strain IMS0002 was grown anaerobically in a mixture of 30 g liter–1 D-glucose, 15 g liter–1 D-xylose, and 15 g liter–1 L-arabinose. Despite high levels of expression of both the xylose isomerase and xylulokinase genes (Fig. 3), strain IMS0002 was not able to ferment xylose into ethanol during an anaerobic batch fermentation with mixed substrates. Instead, a considerable amount of xylitol was produced (data not shown).
|
|
|---|
Long-term selection in synthetic medium with L-arabinose as the sole carbon source resulted in an efficiently L-arabinose fermenting strain (IMS0002) that was able to grow on L-arabinose as the sole carbon source under anaerobic conditions with a specific growth rate of 0.05 h–1 (Table 3), deduced from measurements of dry weight. The increase of 52- and 90-fold in (relative) levels of araB and araD transcripts, respectively, after the extensive selection on arabinose strongly supports the idea that an increase in the activity of L-ribulokinase and L-ribulose-5-phosphate 4-epimerase contributed to the improved (anaerobic) growth on arabinose and that these enzymes may have been limiting the arabinose utilization in the parental strain, IMS0001. In addition, the expression of araA increased, indicating that increased isomerase activity also contributed to the improved utilization of arabinose. The increased level of araB expression, however, is not consistent with the data in an earlier report of Becker and Boles, who described an engineered S. cerevisiae strain selected for growth on L-arabinose (4). Their evolved strain had acquired a mutation in the araB gene, resulting in decreased L-ribulokinase activity compared to that in the original (unevolved) strain. This finding was interpreted as an effect of a high (uncontrolled) level of expression of sugar kinases such as L-ribulokinase and D-xylulokinase on yeast physiology (e.g., rapid ATP depletion and/or toxic levels of phosphorylated sugars) (15). It seems plausible that the levels of activity of the enzymes of the pentose phosphate pathway in strain IMS0002, which were much higher than those in the strain constructed by Becker and Boles, may have prevented the accumulation of the phosphorylated sugars L-ribulose-5-phosphate and/or D-xylulose-5-phosphate. In addition, a higher level of flux through the nonoxidative pentose phosphate pathway, resulting in an increased flux through glycolysis and thereby increased ATP generation, may allow for higher in vivo activity of the initial phosphorylating enzyme L-ribulokinase.
The evolved strain IMS0002 displayed promising kinetics with respect to arabinose consumption and ethanol production during anaerobic cultivation on arabinose or mixed sugars. To our knowledge, the maximum specific arabinose consumption rate of 0.70 g h–1 g (dry weight)–1, the ethanol production rate of 0.29 g h–1 g (dry weight)–1, and the yield of 0.43 g g–1 are the highest rates and yields reported for engineered S. cerevisiae strains (Table 4). In contrast to previously described yeast strains or S. cerevisiae strains engineered for growth on arabinose, strain IMS0002 did not form any detectable arabinitol (Table 3), and this lack of arabinitol formation contributed to high ethanol yields.
Although the arabinose in the anaerobic cultures was completely consumed, a decreasing arabinose consumption rate over time at arabinose concentrations below approximately 50 mmol liter–1 was observed (Fig. 4 and 5). This decrease may indicate a low affinity for arabinose of the arabinose transporter. Alternatively, the relatively high Km of the L. plantarum arabinose isomerase for arabinose (28 mM) (11) may have contributed to the decreasing arabinose consumption rates. In addition, strain IMS0002 exhibited a reduced arabinose consumption rate when grown in mixtures of arabinose and glucose compared to that in the single-sugar fermentation with arabinose as the sole carbon source. This observation is similar to that for xylose consumption by the engineered S. cerevisiae strain RWB217 in mixtures of glucose and xylose (20) and may indicate that the expression of the transporter(s) responsible for arabinose transport (possibly Gal2p [4, 17]) is hampered by the preceding growth on glucose.
The xylose-fermenting S. cerevisiae strain RWB218, which was used as a platform for the present study on L-arabinose fermentation, was shown to retain its xylose fermentation characteristics during several cycles of nonselective batch cultivation on glucose (M. Kuyper et al., unpublished results). Nevertheless, the ability to ferment xylose in anaerobic cultures was lost during the long-term selection for improved L-arabinose fermentation. This finding should not be interpreted to indicate a specific, inherent instability of this xylose-fermenting strain background. During long-term (200-generation) glucose-limited, aerobic chemostat cultivation of wild-type S. cerevisiae, even the high glycolytic capacity of wild-type S. cerevisiae is strongly reduced (13, 25). The latter observation has been interpreted as the result of selection for spontaneous mutants that synthesize fewer "excess," energetically expensive glycolytic proteins. Prolonged cultivation (>1,000 generations) on a carbon source other than xylose may well exert a similar selective pressure on the expression of genes involved in xylose metabolism. Future work will focus on functionally expressing the two pentose pathways in a single strain and on further improvement of the kinetics of mixed-sugar consumption and ethanol production.
We thank Marinka Almering for her contribution to the quantitative real-time PCR experiments.
Published ahead of print on 1 June 2007. ![]()
|
|
|---|

CT method. Methods 25:402-408.[CrossRef][Medline]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»