Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2008, p. 2404-2413, Vol. 74, No. 8
0099-2240/08/$08.00+0 doi:10.1128/AEM.02469-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Hyung Soo Park,1,
Jenny Been,2,
Paul Gordon,3
Christoph W. Sensen,3 and
Gerrit Voordouw1*
Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada,1 NOVA Research and Technology Corporation, 2928 16th Street NE, Calgary, Alberta, Canada,2 Sun Center of Excellence for Visual Genomics, Department of Biochemistry and Molecular Biology, University of Calgary, Calgary, Alberta, Canada3
Received 1 November 2007/ Accepted 15 February 2008
|
|
|---|
|
|
|---|
Fe2+ + 2e). In addition, oxygen is reduced under aerobic conditions, while protons are reduced under anaerobic conditions at cathodic sites drawing electrons from the metal (7). These two processes, schematically depicted in Fig. 1A, give rise to an anodic and cathodic current, iA and iC, respectively. These currents are of opposite signs and depend on the metal's potential (E), as shown in Fig. 1B. Hence, the total current measured in an external wire connected to the electrode would be iT = iA + iC, and iT = 0 at the lines' point of intersection (Fig. 1B). The identical absolute values iA and iC represent the corrosion current, icorr, and the corresponding potential is the corrosion potential, Ecorr. The metabolic activity of microorganisms can contribute to the corrosion process, and sulfate-reducing bacteria (SRB) were implicated as early as 1934 (37). SRB can increase the icorr value by using cathodically produced hydrogen to reduce sulfate to sulfide or by using produced sulfide to precipitate anodically produced ferrous ions as FeS. However, produced FeS may also form a film on the iron surface, reducing the icorr value (7). As indicated in Fig. 1B, imposing a negative potential on the iron surface boosts formation of cathodic hydrogen (iC) and lowers the corrosive flux of ferrous ions (iA). This procedure, known as cathodic protection (CP), is commonly practiced to lower the corrosion rate of steel pipelines. Most hydrogenase-containing SRB, like the well-known Desulfovibrio spp., do not effectively use electrons derived from iron as the electron donor for sulfate reduction at potentials at or above the Ecorr value (E
–0.7 V). Under these conditions, these SRB can only oxidize iron cometabolically in the presence of an organic electron donor for sulfate reduction, such as lactate (6). Desulfobacterium strains using iron efficiently as the sole electron donor for sulfate reduction (i.e., in the absence of an organic electron donor) have been isolated (12). It was suggested that these may be able to channel electrons directly from iron to the cytoplasmic sulfate reduction pathway through outer membrane electron transport complexes, as used by Fe(III)-reducing bacteria (12).
![]() View larger version (15K): [in a new window] |
FIG. 1. Anaerobic corrosion of iron. (A) Anodic and cathodic currents. (B) Dependence of log|iA| and log|iC|, where |iA| and |iC| are the absolute values of iA and iC, on electrode potential E. The intersection point has the coordinates (log|icorr|, Ecorr). (C) Tafel plot of actual experimental data in which E is plotted against log|iT|, the absolute value of the measured current: |iT| = ||iA| – |iC||.
|
|
|
|---|
Probe sequence creation and microarray manufacturing.
A total of 2,824 70-mer oligonucleotides were designed using the Osprey software program (18), one for each designated probable CDS (see Table S2 in the supplemental material). Their target melting temperature in 0.1 M NaCl was 87°C, matching that of oligonucleotides for a previously designed D. vulgaris microarray (5). The optimal design conditions for our oligonucleotides were as follows, in descending order of importance: (i) a melting temperature within 5°C of the target (i.e., from 82 to 92°C); (ii) no secondary binding of the oligonucleotide elsewhere in the genome with a melting temperature within 15°C of that for the target site; (iii) no formation of hairpins or dimers with a stability of more than 19 kcal/mol in 0.1 M NaCl; (iv) no long high-G+C stretches, to prevent "sticky" hybridization. Optimal oligonucleotide probes satisfying all four design conditions were found for 1,180 genes. A survey of design conditions violated by suboptimal probes is provided elsewhere (see Fig. S1F in the supplemental material). Additionally, eight oligonucleotides matching the optimal conditions but not targeting a CDS were included as negative hybridization controls.
The designed 70-mer oligonucleotides were custom synthesized by Invitrogen (Carlsbad, CA) on an average scale of 25 pmol/oligonucleotide. They were dissolved in 3x standard saline citrate (SSC) (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate) to a concentration of 40 pmol/µl and spotted in 0.7-nl quantities, using a ChipWriter Pro spotter (Bio-Rad, Hercules, CA), onto Corning UltraGAPS amino-silane 25-by-75-mm coated slides (Fisher Scientific). The spotting was done at the Southern Alberta Microarray Facility (SAMF). These arrays will be referred to as SAMF arrays. Spot quality was assessed by visual inspection of slides hybridized to cyanine-3-labeled random nanomers. In order to prevent localized hybridization problems from skewing the expression results, the CDS probes were laid out in a random pattern to separate probes for genes located in the same operon. All probes were spotted in duplicate with duplicate spots located in different halves of the array grid. In addition to spotting probes targeting D. vulgaris CDSs, several control probes were printed. Negative controls included blank spots, 3x SSC, and 70-mer oligonucleotide probes without homology to any D. vulgaris CDS. The positive-control spots consisted of 70-mer oligonucleotides probes to human and Arabidopsis genes for which RNA is available to add as spiked-in controls.
RNA and DNA extraction, labeling, hybridization, and data analysis.
The RNA was purified from D. vulgaris cultures with the Qiagen RNA Easy kit, and the RNA quality and quantity were assessed as described previously (3). Genomic DNA (gDNA) was isolated as described elsewhere (39). Amino-modified cDNA was generated from 10 µg of total RNA using the Invitrogen SuperScript indirect cDNA labeling system and random hexamers. Amino-modified gDNA was generated using 1 µg of gDNA with the Invitrogen BioPrime Plus Array CGH genomic labeling system. The amino-modified cDNA was then coupled with Alexa fluor 647 and the amino-modified gDNA with Alexa fluor 555. The dried cDNA and gDNA probes were mixed and diluted to 20 µl using hybridization solution (3). Following incubation, washing, and drying (3), the slides were scanned with a PerkinElmer Scanarray 5000 scanner (PerkinElmer, MA). Spot and background intensities were quantified using JCVI's Spotfinder software (32) with the gDNA channel being used to control for spot uniformity and labeling efficiency (5, 34), as described previously (3). Array analysis used RNA extracted from two independent biological samples, which were hybridized to a minimum of four arrays, each containing duplicate spots. Data were normalized (3), and differentially expressed genes were identified, using the Stanford Tools software package (32), as those with an expression change of at least twofold with a q value of less than 5%. Functional categories for global analysis were defined by clusters of orthologous genes (COGs) (35). A bootstrapped supported-hierarchical-clustering tree was created using JCVI's multiple experiment viewer from genes identified as significant by a one-way analysis of variance with adjusted Bonferroni correction at the 1% confidence level (32). The clusters were created using Euclidean distance and average linkage.
Strains, media, and culture conditions.
Desulfovibrio vulgaris subsp. vulgaris Hildenborough (NCIMB 8303), lacking the 202-kb megaplasmid (21), was the wild-type strain from which D. vulgaris Hyd100 (
hydAB Cmr), D. vulgaris Hyn100 (
hynAB), and D. vulgaris H801 (
hmc) were derived by mutagenesis as described elsewhere (13, 17, 29). These were grown in Widdel-Pfennig medium with lactate and sulfate (WP-LS) and a gas phase of 10% (vol/vol) CO2 and 90% (vol/vol) N2 or in WP-HAS medium with 5 or 50% (vol/vol) hydrogen as the electron donor in 10% CO2 and balance N2 (3, 38, 40). The latter medium had, in addition to CO2, 3 mM acetate as the carbon source. In cases where the hydrogen was supplied electrochemically, this medium will be referred to as WP-EAS medium with the value of E indicated. The experimental setup for batchwise growth with a continuous gas flow of hydrogen (100 ml/min per bioreactor) was as described elsewhere (3). For electrochemical studies, a 5% (vol/vol) inoculum of a log-phase culture, grown in WP-LS or WP-HAS, was used as the inoculum.
Electrochemical analyses.
Mixed gas (10% [vol/vol] CO2, balance N2), which was passed through a heated copper reduction tube (350 °C) to remove traces of oxygen, was bubbled into a glass biocorrosion reactor with a working volume of 400 ml at a flow rate of 50 ml/min (Fig. 2). A conventional three-electrode system was used. A low-carbon steel coupon (ASTM A366; 8 by 1 by 0.1 cm) served as the working electrode (WE). The WE was polished with ethanol, dried in air, and sterilized (160°C, 2 h). A platinum helical wire was used as the counter electrode (CE), whereas an Ag/AgCl electrode (MF-2079; Bioanalytical Systems, Inc., West Lafayette, IN) was used as the reference electrode. A diagram of the experimental setup is shown in Fig. 2. Rates of corrosion (Fig. 1B, the anodic loss of metal represented by iA) are reported for Ecorr, the potential where iA = iC = icorr. Values for icorr and Ecorr are derived from the dependence of log |iT| on imposed potential, as represented in Tafel plots (Fig. 1C), or by fitting potentiodynamic curves to equations for iA and iC (28). Data were acquired with a model 410 potentiostat (Electrosynthesis Co, NY) coupled to a multimeter (Keithley 2700), which monitored the current using Testpoint software (Fig. 2). Corrosion rates (CR) in milli-inches per year (mpy) were calculated as follows: CR = (icorr x K x EW)/(d x A), where K is 1.288 x 105, EW is the equivalent weight of iron (55.8 g per 2 mol of electrons), d is the density of iron (7.8 g cm–3), and A is the electrode area in cm2. Alternatively, icorr was derived from electrochemical impedance spectra recorded with an SI1260 or SI1287 electrochemical interference/impedance gain-phase analysis spectroscope (Solartron, United Kingdom) in the frequency range of 0.05 to 100,000 Hz. For large-scale biomass production, D. vulgaris was grown on an iron electrode with a total surface area of 17.1 cm2 in WP-EAS medium (E = –1.1 V) for approximately 2 days, at which point half of the sulfate concentration in the bioreactor was reduced. The electrode-attached biomass was scraped off with a scalpel under anaerobic conditions, after which RNA was extracted as described before (3).
![]() View larger version (33K): [in a new window] |
FIG. 2. Schematic of the biocorrosion reactors. The reactors were gassed with 10% (vol/vol) CO2 and 90% N2 at a rate of 50 ml/min per bioreactor. Each bioreactor contained an iron WE, a platinum CE kept in a separate glass compartment sealed off with an agar plug, and an Ag/AgCl reference electrode (RE). The experimental setup allows the potential (EWE) to be changed while recording the current (iT).
|
Microarray data accession number.
The microarray data presented in this paper have been entered into NCBI GEO (14) under accession number GSE10388.
|
|
|---|
Reanalysis of the D. vulgaris Hildenborough genome indicated three probable genes not previously annotated (Fig. 3). DVU0695_0696 replaces two small ORFs (Fig. 3A) (DVU0695 and DVU0696) and encodes a 211-amino-acid (aa) multiheme cytochrome in a probable operon with DVU0694 to DVU0692 for molybdopterin oxidoreductase. It is conserved in strain DP4 (Dvul_2268; 211/211) and in D. desulfuricans G20 (Dde_2932; 147/211) in what appears to be the same operon. DVU0860_new encodes a predicted glycosyltransferase (420 aa) and replaces the 33-aa ORF DVU0860 (Fig. 3B). Dvul_2122 and Dde_1117 are its homologs in the other sequenced Desulfovibrio spp. DVU2047_2048 replaces DVU2047 and DVU2048 (Fig. 3C) and encodes a predicted 684-aa DNA repair protein with homology to a DP4 gene (Dvul_2616) but no G20 homology. Probes targeting these new D. vulgaris genes were spotted on the array, and both DVU0695_0696 and DVU0860_new, but not DVU2047_2048, were subsequently found to be expressed (see Table S4 in the supplemental material).
![]() View larger version (27K): [in a new window] |
FIG. 3. Survey of three newly annotated genes found by MAGPIE in the D. vulgaris genome, as explained in the text. DVU0695_696 replaces DVU0695 and DVU0696 (A), DVU0860_new replaces DVU0860 (B), and DVU2047_2048 replaces DVU2047 and DVU2048 (C). The new locus tags (ID), their location in the genome (bp), the length of the encoded protein (aa), and their closest homolog, as identified by BLAST search, are indicated.
|
The quality of the SAMF microarrays was further tested by checking the correlation of differential expression ratios (RC,i) for genes in the same operon. This was computed as the Pearson product-moment correlation between the percentile rank of raw spot intensities (foreground minus background). Pearson's method is unaffected by downstream effects (reduced expression of 3' genes compared to that of 5' genes in many Desulfovibrio operons). Using the 639 predicted D. vulgaris operons (1), around 350 genes per replicate were considered "expressed" using a 1-mRNA/gDNA intensity as the threshold level. Random pairs of genes from these operons were tested for product-moment correlation (R, the tendency of genes to go up or down together linearly), and the average correlation for all operons was calculated. The random pair selection process was repeated 1,000 times to generate a confidence interval of R = 0.711 to 0.717, using a data set comprised of 5 hybridizations. The correlation implied that approximately 0.51% (R2 = 0.7142) of the expression change of a gene can be accounted for by other genes in the same operon. Performing the same analysis for data obtained with the VIMSS microarrays yielded data with a confidence interval of R = 0.714-0.718, suggesting data from both microarray platforms to be equally trustworthy and comparable to those published to date. These correlations are better than the similarly derived Spearman's confidence interval (R = 0.65 to 0.67) for the best analysis method for the Affymetrix Escherichia coli gene chips (19).
Global gene expression in electrode-grown D. vulgaris.
Growth of D. vulgaris on an iron electrode suspended in WP-EAS medium using the culturing system outlined in Fig. 2 required a low working electrode potential (EWE) of –1.1 V. Hydrogen evolves from the electrode, the anodic current (iA) is negligible (Fig. 1B), and the electrode is protected from corrosion (Fe0
Fe2+ + 2e) under these CP conditions. Growth was not observed at higher electrode potentials, e.g., at an EWE value of –0.7 V, corresponding to Ecorr (Fig. 1C), when the bioreactor was filled with WP-EAS medium. RNA was extracted at the point where half of the sulfate of the medium was used. For wild-type cells grown under CP conditions, 85% of total RNA was extracted from the electrode, whereas for cells grown in WP-LS medium at an E value of –0.3 V, this was only 9% (see Table S5 in the supplemental material). Hence, under CP conditions, most cells grew as an electrode-associated biofilm.
When the gene expression pattern of cells grown under CP conditions was compared with that of cells grown with gaseous hydrogen, a smaller shift in the gene expression pattern was observed (56 genes up-regulated and 40 genes down-regulated) than when that of cells grown under CP conditions was compared with that of cells grown with lactate (193 genes up- and 229 genes down-regulated) (Table 1). The bootstrap-supported hierarchical clustering tree (Fig. 4) also illustrates that the gene expression pattern in CP-grown cells is more similar to that of hydrogen-grown cells than to that of lactate-grown cells. CP-grown cells showed the greatest gene expression changes in the following COG categories: viii, cell envelope biogenesis, outer membrane; ix, cell motility and secretion; xiii, energy production and conversion; and xv, amino acid transport and metabolism.
|
View this table: [in a new window] |
TABLE 1. Numbers of genes up- or downregulated in D. vulgaris wild-type cells grown under different conditionsa
|
![]() View larger version (10K): [in a new window] |
FIG. 4. Bootstrap-supported hierarchical clustering tree for gene expression patterns for wild-type D. vulgaris grown in WP-LS (lactate), WP-HAS (50% [vol/vol] or 5% [vol/vol] hydrogen), or WP-EAS (EWE = –1.1 V; cathodic protection conditions, cells grown on cathodic hydrogen).
|
|
View this table: [in a new window] |
TABLE 2. Selected genes that are up- or down-regulated for wild-type D. vulgaris cells grown under different conditionsa
|
Determination of corrosion rates.
The corrosion current (icorr) was determined without bacteria present under both anaerobic and aerobic conditions (Fig. 1C). The derived values, converted to corrosion rates (CRs), were 4.6 mpy under anaerobic conditions and 10.3 mpy under aerobic conditions. Under the latter conditions, oxygen reduction is part of the cathodic reaction (2H+ + O2 + 4e [electrode]
2 OH–). The anaerobic abiotic rates were constant with time, as indicated in Fig. 5C and 6B.
![]() View larger version (13K): [in a new window] |
FIG. 5. (A to C) Anaerobic corrosion rates in biocorrosion reactors filled with WP-LS medium as a function of time. The reactors were inoculated with D. vulgaris wild-type (WT) or mutant strains, as indicated. The iron coupon working electrodes had a measured EWE value of –0.3 V. Corrosion rates are reported for an EWE (= Ecorr) value of –0.7 V. (D to F) Concentrations of sulfate and sulfide as a function of time.
|
![]() View larger version (16K): [in a new window] |
FIG. 6. (A and B) Anaerobic corrosion rates in biocorrosion reactors filled with WP-EAS medium as a function of time. The reactors were inoculated with D. vulgaris wild-type (WT) or mutant strains, as indicated. The iron coupon working electrodes had a set EWE value of –1.1 V. Corrosion rates are reported for an EWE (Ecorr) value of –0.7 V. (C and D) Concentrations of sulfate and sulfide as a function of time.
|
The results of three experiments conducted in WP-LS medium, one conducted with low concentrations of lactate and sulfate (5 and 3 mM, respectively; Fig. 5A) and two conducted with high concentrations of lactate and sulfate (20 and 10 mM, respectively; Fig. 5B and C), are summarized in Fig. 5. In the case of all three experiments, the initial CR was close to the anaerobic abiotic CR of 4.6 mpy. The presence of wild-type or mutant D. vulgaris strains actively reducing sulfate transiently increased the CR (Fig. 5A to C). When sulfate reduction stopped, because the limiting concentration of lactate had been used (data not shown), the CRs settled to values at or below that of the abiotic control. Hence, the presence of actively growing Desulfovibrio strains in the bioreactors, which were in part present as an electrode-associated biofilm (see Table S5 in the supplemental material; 9% for the wild type), increased the CR. The reasons for the larger increase in the CR for the wild type at low lactate and sulfate concentration (Fig. 5A;
CR = 20 mpy) than at high lactate and sulfate concentrations (Fig. 5B and C;
CR = 3 to 5 mpy) are not known. The hyn1, hyd, and hmc mutant strains increased the CR less than the wild-type strain, suggesting that [NiFe] hydrogenase 1, [Fe] hydrogenase, and the Hmc complex contribute to the corrosion-enhancing metabolism of D. vulgaris.
The results of experiments conducted in WP-EAS medium under CP conditions are presented in Fig. 6. The wild type had a significantly increased CR, which was not observed for the hyn1 mutant (Fig. 6A), whereas the hmc mutant and the wild type increased the CR similarly (Fig. 6B). Sulfate reduction with associated growth could not be established for the hyd mutant, lacking [Fe] hydrogenase, suggesting this enzyme to be essential for growth under CP conditions.
|
|
|---|
Because the hydrogen electron donor is exclusively generated at the electrode under CP conditions, the majority of the biomass is electrode associated (see Table S5 in the supplemental material). Bubbles emerging from the electrode surface did not effectively feed the planktonic population, since these were removed by the continuous gassing with N2-CO2, which also effectively stripped H2S (Fig. 5 and 6). Hence, gene expression changes of hydrogen-grown versus CP-grown cells will in part reflect the transition from the planktonic lifestyle to the biofilm lifestyle. This transition involves cells that are planktonic and motile (state A), surface associated with gliding motility through additional flagella (state B), and in a surface-attached, nonmotile, mature biofilm state (state C) (4). States A, B, and C can have increasing concentrations of the nucleotide signal cyclic di-GMP, as demonstrated for different bacteria (4, 23, 36). This signal is formed by GGDEF-domain proteins and removed by EAL-domain proteins, which synthesize and hydrolyze cyclic di-GMP, respectively. Bifunctional GGDEF/EAL-domain proteins also occur commonly. The fact that proteins in this class are among the most up-regulated (Table 1; DVU1959 and DVU3106) and the most down-regulated (see Table S4 in the supplemental material; DVU0763) of the entire genome, when CP cells and planktonically grown cells are compared, may reflect this lifestyle change.
The up-regulation of many genes for flagellar proteins in CP-grown cells (Table 2; see also Table S4 in the supplemental material), compared to cells grown in 50% (vol/vol) H2, indicates that the former are in the surface-associated, motile state (state B), not in the nonmotile, mature biofilm state (state C). The increased flagellar gene expression may allow effective gliding motility in the much more viscous substratum-associated biofilm (4, 23, 36). Relative to mature biofilm cells, cells in state B overexpress few genes for carbohydrate synthesis. In contrast, cells in state C have reduced expression of flagellum genes and increased expression of genes for enzymes involved in exopolysaccharide biosynthesis. Zhang et al. (42) have reported the gene expression pattern of D. vulgaris biofilms grown on iron coupons suspended in a continuous-culture bioreactor fed with lactate- and sulfate-containing medium for a period of 20 days. The biofilm that formed on these coupons was characterized by reduced expression of flagellum genes and increased expression of many genes involved in polysaccharide biosynthesis, indicating that these represent the mature biofilm state (state C). The striking difference between gene expression in the two biofilm types is presented in Table S7 in the supplemental material.
CP-grown biofilm and 50% (vol/vol) hydrogen-grown planktonic cells engage in similar energy metabolism with CP-grown cells experiencing a cell density that is 2 orders of magnitude higher. Perhaps due to the increased competition experienced by CP-grown cells, they overexpress virtually all known components of this energy metabolism pathway (27) relative to cells grown planktonically with 50% (vol/vol) H2, including genes for three of four hydrogenases, for the Hmc complex but not the other membrane-bound electron transfer complexes, for ATP synthase, and for the main enzymes of the sulfate reduction pathway (Table 2; see also Table S4 in the supplemental material). Of the three nickel-containing hydrogenases, genes for the large and small subunits of Hyn-1 are among the most up-regulated of the entire genome and those for Hys are also up-regulated, whereas those for Hyn-2 are down-regulated. The hyd genes for iron-only hydrogenase are also up-regulated under CP conditions, and the special enzymatic properties of this enzyme (high turnover, low affinity) may be important to explain why the hyd mutant could not be grown under CP conditions.
SRB may increase the cathodic and anodic currents (Fig. 1A). The lower CRs observed in D. vulgaris strains missing either Hyn1, Hyd, or the Hmc complex (Fig. 5 and 6) illustrate that these proteins facilitate metabolism that enhances iA, iC, or both. The gene expression results further demonstrate this connection, since the hyn1, hyd, and hmc genes are upregulated under CP conditions. The lower CRs observed when SRB metabolism has ceased (Fig. 5C) may be explained by formation of an iron sulfide film on the iron surface area, reducing the rate of the corrosion reaction (7, 25).
D. vulgaris does not appear to have an outer membrane electron transport complex (21). The Hmc complex and other membrane-bound electron transfer complexes are located in the inner membrane. Differently from Fe(III)-reducing Geobacter and Shewanella spp., which have outer membrane electron transport complexes, D. vulgaris cannot reduce extracellular Fe2O3 (26). The results obtained in this study similarly indicate that the organism does not grow by accepting electrons directly from metallic iron. Growth requires either an imposed low potential (EWE = –1.1 V), conditions under which hydrogen evolves from the metal surface, or, at higher potential (e.g., EWE = –0.3 V), the presence of the organic electron donor (lactate) to drive reduction of sulfate to sulfide. Corrosion (iron oxidation) is cometabolic with lactate oxidation under the latter conditions. SRB that use iron efficiently as an electron donor for sulfate reduction in the absence of an organic electron donor, like the Desulfobacterium strains identified by Dinh et al. (12), may be able to directly extract electrons from the iron surface, and it would be most interesting to elucidate the mechanism used.
We thank Viola Birss, Department of Chemistry, University of Calgary, for electrochemistry support. Mayi Arcellana-Panlilio and Xiuling Wang of the Southern Alberta Microarray Facility are thanked for array printing and hybridization advice. Hong Tran is thanked for helping with the hand curation of the D. vulgaris genome.
Published ahead of print on 29 February 2008. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
These two authors contributed equally to this article. ![]()
Present address: Corrosion Engineering, Advanced Materials, Alberta Research Council, 3608 33rd Street NW, Calgary, Alberta T2L 2A6, Canada. ![]()
|
|
|---|
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»