Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2008, p. 7585-7595, Vol. 74, No. 24
0099-2240/08/$08.00+0 doi:10.1128/AEM.01221-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Peter G. Green,2 and
Patricia A. Holden1*
Donald Bren School of Environmental Science and Management, University of California, Santa Barbara, California 93106-5131,1 Department of Civil and Environmental Engineering, University of California, Davis, California 956162
Received 2 June 2008/ Accepted 18 September 2008
|
|
|---|
|
|
|---|
Though denitrifier communities in soil (5, 49, 51, 64, 77) and in marine sediments (8, 10, 45, 53, 68), where they may be influenced by season, nutrient concentration, soil types, redox gradients, and fertilizer types, have been widely studied, those in salt marsh systems have received little attention. Salt marshes differ fundamentally from soils and marine sediments in that marsh sediments are tidally influenced with associated redox potential, salinity, and inundation fluctuations (50). Studies of total bacterial communities in salt marsh sediments indicate that while they are highly spatially variable, they are significantly influenced by geographic location, elevation, and selected pollutants, such as metals (12, 18). Yet, rhizosphere diazotroph communities were shown to be relatively stable and persistent (46) in the face of short-term carbon and nutrient variability (61, 62). How denitrifier communities and their functions may vary with environmental and pollutant factors in salt marshes remains largely unknown. Therefore, the objective of this study was to examine, for a salt marsh system with high nitrate loading (57), the spatial gradients of microbial community composition and denitrification enzyme activities (DEAs), the relationships of community composition and function to environmental factors as well as to each other, and the potential effects of metals considered here to be sedimentary pollutants that could affect denitrifiers and denitrification (12, 18).
|
|
|---|
![]() View larger version (81K): [in a new window] |
FIG. 1. Study location including 12 sampling stations. The stations are divided into four groups: stations along the main channel of West Creek (A, B, C, D, E, and I), stations on branches of West Creek (G and J), stations along the central road (H, F, and K), and a station further east of the marsh (L).
|
Sediment characterization.
Moisture content and total organic carbon (TOC) were measured for all 36 samples. Moisture content was determined for triplicate sediment samples (1.5 g each) according to standard methods (23). Dried (105°C; O/N) sediments were combusted in a muffle oven at 450°C to determine TOC (52). The remaining sediments were preserved at –20°C until microbial community and sediment eluate analyses (as described below).
Dissolved nutrient analysis.
Bulk water samples collected from surfaces immediately above the sediment at the low elevations were analyzed for concentrations of nitrate (38), phosphate (as soluble reactive phosphorus [39] and ammonium [37]), using a model 2300-000 ion analyzer (QuikChem II; Lachat, Milwaukee, WI). Before nutrient analysis, the water samples were filtered through a filter sandwich (0.7-µm Whatman GF/F filter overlaying a 0.2-µm Millipore filter) to eliminate the suspended particles and microbes in the samples. Filtered water samples were stored at –20°C until analysis (57).
Sediment eluate analyses.
Bioavailable metals in sediment eluates were quantified. To recover eluate, 15 g of thawed and homogenized sediment was weighed into 50-ml polypropylene centrifugation tubes and 15 ml chilled (4°C) sterile artificial seawater (ASW) was added. After an initial brief vigorous shaking to ensure thorough mixing, samples were horizontally shaken (2 h, 4°C) and then centrifuged (2500 x g, 30 min, 4°C). Eluted water samples were collected and analyzed for the concentrations of metals and dissolved organic carbon (DOC). Eluates were similarly collected and analyzed for additional fresh sediment samples to determine the effect of the sediment freeze-thaw process on DOC and metal concentrations. For all samples, analysis controls included either ASW or ultrapure water (carbon and ion free) in either centrifuge tubes or combusted (450°C; O/N) glass vials.
An aliquot of each eluate (250 µl) was diluted 50-fold in 1% nitric acid (trace metal grade; Fisher Scientific, Pittsburgh, PA) in a new polypropylene centrifuge tube. A blank tube was measured as a control, confirming negligible contamination. Metal quantification was performed by direct aspiration using an inductively coupled plasma mass spectrometer (model 7500i; Agilent, Palo Alto, CA) and compared to external standards, covering a range from 0.01 µg l–1 to 200 µg l–1 for minor elements and up to 100 mg l–1 for major elements. Analysis was performed as prescribed by EPA method 6020 (74). The instrument was operated with very low-charge molecular ions, a robust plasma characterized by 0.4% CeO/Ce. For DOC, eluates were diluted 100-fold and analyzed via high-temperature combustion, using a Shimadzu TOC-V instrument (Shimadzu Scientific Instruments, Columbia, MD), as described previously (13). Three to five replicate 100-µl samples were injected into the combustion tube for analysis. The injections continued until at least three injections met the specified levels (a standard deviation of 0.1 area counts and a coefficient of variation of
2%); otherwise, the best three of five injections were considered. Both DOC and metal concentrations were corrected for analysis controls and normalized to sediment dry weight.
DEA.
DEA was measured using a standard method (72) for triplicates of each sediment sample within 36 h of sample collection. Briefly, ca. 3 g of sediment was placed in a 15-ml sterile, airtight, Hungate-type anaerobic culture tube (Bellco Glass, Inc., Vineland, NJ) containing 1 mM glucose, 1 mM KNO3, and 1 g l–1 chloramphenicol. The headspace was flushed with 7 volumes of filter-sterilized (0.2 µm) helium. Acetylene (10% in He) was added with a syringe through the septum, and the sediment slurry was incubated on a rotary shaker (30°C, 2.5 h). During the incubation, the headspace was sampled (3 ml) and replaced three times with 10% acetylene at 0.5, 1.5, and 2.5 h. The headspace concentration of N2O was determined by gas chromatography with an electron capture detector (GC-14A; Shimadzu Corporation, Columbia, MD), using a chromatography laboratory automated software system (CLASS-VP; Shimadzu Corporation, Columbia, MD) from a calibration curve generated with N2O standards (Scott Specialty Gases, Inc., Plumsteadville, PA). The above-mentioned sampling and incubation times were determined experimentally to achieve the best linearity of DEA for the sediment samples (72).
Abundance of bacterial cells.
Enumeration of total bacterial cells was performed in triplicate, using a method similar to that of Holden et al. (25). Briefly, 150 µl of the sediment suspension was stained with a nucleic acid-specific stain (Sybr gold; Molecular Probes, Eugene, OR) and then suction filtered onto a 13-mm-diameter, 0.2-µm-pore-size black Isopore membrane (Millipore, Bedford, MA). Total cells were directly counted using a Nikon E800 epifluorescence microscope with a green fluorescent protein filter set (Chroma, Brattleboro, VT); counts were normalized to sediment dry weight.
Microbial community analysis.
The total bacterial and denitrifier communities were profiled via terminal restriction fragment length polymorphisms (TRFLPs) as described previously (40). Briefly, 2 g of sediment was extracted using an UltraClean soil DNA kit (MoBio, Solana Beach, CA), following the manufacturer's protocol. Extracted DNA was purified (by size exclusion with Sepharose spin columns) and checked for quality using A260/230 spectrophotometry and for quantity using Picogreen (Molecular Probes, Eugene, OR) fluorometry. For profiling of the total bacterial community, genes encoding 16S rRNA were PCR amplified from purified community DNA by using universal eubacterial primers 8F hex (fluorescently labeled forward primer; 5'AGAGTTTGATCCTGGCTCAG) (44) and 1389R (5'ACGGGCGGTGTGTACAAG) (56). For the denitrifying community, genes encoding the nitrite reductase NirS were PCR amplified using the primers nirS1F hex (fluorescently labeled forward primer; 5'CCTAYTGGCCGCCRCART) (9) and nirS6R (5'CGTTGAACTTRCCGGT) (9) as described previously (10). PCR products were purified and digested using either HhaI or MspI in two separate digestions (21). The restriction enzymes were inactivated (65°C, 20 min). The lengths of fluorescently labeled terminal restriction fragments (TRFs) were determined with an Applied Biosystems Instruments (ABI; Foster City, CA) model 373A automated sequencer. Data from HhaI and MspI digestions were separately aligned and normalized using a method similar to that of Dunbar et al. (20) and then concatenated to form a collective data set. All TRFLP data were analyzed in binary form (i.e., presence and absence of TRFs) unless otherwise indicated. Data analyses using relative abundances of TRFs (i.e., normalized to the total fluorescence peak height) were conducted for comparison. Numbers of TRFs were tabulated for a comparative measure of community richness.
Data analyses.
All data analyses were conducted using either S-plus version 6.1 (Insightful Corp., Seattle, WA), Primer version 6 (PRIMER-E Ltd, Plymouth, United Kingdom), or CANOCO version 4.53 (Microcomputer Power, Inc., Ithaca, NY). Continuous variables (organic content and metal concentrations) were standardized to zero mean and unit variance before either regression or multivariate correlation analyses in order to account for differences in the values of these variables. The correlation structure among metal concentrations was determined by calculating pairwise Pearson's correlation coefficients with S-plus. Effects of elevation on DOC, TOC, and metal concentrations were evaluated via analysis of variance. The relationship between any two continuous variables was investigated by simple regression analysis. To investigate the effects of metals on a single biological variable (e.g., DEA and total cell counts, as opposed to multivariate community profiles) in the presence of other environmental variables (such as DOC, TOC, and elevation), the residual variability in either DEA or total cell counts (after adjustment for carbon and/or elevation effects) was correlated with concentrations of specific metals (36) with S-plus. Bray-Curtis and Euclidean distances were calculated for microbial community TRFLP profiles and for abiotic characteristics (i.e., environmental variables and metal concentrations), respectively, to summarize their relative spatial distributions in the marsh. Analyses of correlation between these distance measures were then conducted through the RELATE program in Primer version 6 (16) to compare, for example, the spatial distributions of microbial communities in the 16S TRFLP and the nirS TRFLP profiles for all samples or at each elevation. The Spearman coefficient (32) was used as a measure of correlation, and P values were calculated via permutation.
Three multivariate techniques based on correspondence analysis were employed to analyze the TRFLP profiles (58, 69, 70). Detrended correspondence analysis (DCA) was used to characterize the overall variation in the TRFLP profiles. Canonical correspondence analysis (CCA) was used to directly assess the relationships between microbial community profiles and environmental variables such as elevation, TOC, and DOC, with dummy variables included to account for sample elevation. CCA was also used to examine relationships between TRFLP profiles and other biotic measures, such as DEA and total cell counts. A third technique, partial CCA (pCCA), was used to investigate the effects of metals on microbial communities after adjustment for natural variability due to elevation and carbon content. Briefly, variables representing elevation and carbon content were designated covariables, and the residual variabilities in the TRFLP profiles were correlated with metal concentrations, as previously described (12).
|
|
|---|
![]() View larger version (20K): [in a new window] |
FIG. 2. Dissolved nutrient concentrations in bulk water samples overlying low-elevation sampling sites at each station: nitrate (A), ammonium (B), and phosphate (C). Concentrations were based on single-sample measurements. Inset graphs exclude station A to emphasize differences across the remaining stations. Vertical dashed lines separate the groups of stations described in the legend to Fig. 1.
|
![]() View larger version (28K): [in a new window] |
FIG. 3. TOC (A) and DOC (B) for three elevations (high, middle, and low) at each station. Vertical dashed lines separate the groups of stations described in the legend to Fig. 1. Asterisks in panel A denote stations where TOC differed significantly across elevations. DOC differed significantly across elevations in all 12 stations. The bars represent the mean values, and the error bars represent the standard deviations (n = 3).
|
0.85) existed between the concentrations of the representative metals and the remaining metals: Al-Fe, 0.94; Al-Rb, 0.89; Al-As, 0.88; Ca-Mg, 0.95; Ca-Sr, 0.94; Ca-K, 0.92; Ca-Rb, 0.86; Ca-As, 0.85; Ca-Fe, 0.88; Cr-Mg, 0.90; Cr-K, 0.89; Cr-Rb, 0.89; Cr-As, 0.88; Li-Mg, 0.94; Li-K, 0.88; Ni-As, 0.96; Ni-Rb, 0.86; Ni-Sr, 0.85; Se-Mg, 0.91; Se-K, 0.91; and Se-Rb, 0.87. Also, many of the concentrations of representative metals were well correlated: Ca-Se, 0.86; Ca-Ni, 0.87; Ca-Cr, 0.85; Ca-Li, 0.85; Cr-Li, 0.87; Cr-Ni, 0.87; Cr-Se, 0.86; and Cr-Ca, 0.85. Cu and Ni were also correlated, albeit comparatively weakly (0.80). Among the 14 representative metals, Zn was significantly positively related to both DOC (R2 = 0.13; P = 0.030) and TOC (R2 = 0.18; P = 0.010), and V was significantly positively related to DOC (R2 = 0.11; P = 0.047). |
View this table: [in a new window] |
TABLE 1. Metal concentrations in sediment eluatesa
|
s = 0.21) (Table 3) and along the 12 stations at the middle elevation (
s = 0.35) (Table 3), suggest that the average trends in spatial distribution throughout the marsh were similar for the two communities. |
View this table: [in a new window] |
TABLE 2. Numbers of TRFs in the 16S TRFLP and nirS TRFLP profiles
|
|
View this table: [in a new window] |
TABLE 3. Correlations between 16S and nirS TRFLP community profiles for all samples at each elevation
|
Because significant natural variability (i.e., related to elevation and carbon content) was present in the microbial community profiles, three variables (high elevation, low elevation, and DOC) were included as covariables for correlating metal concentrations with the community profiles via pCCA. The relationships of metals to community profiles were determined in three ways: by forward selection, by manual selection, and by correlation with metals that were selected by the previous two methods. More emphasis is placed on relationships determined by the first two methods because those metals were directly correlated with community composition while those in the last category were inferred. Via automated forward selection in pCCA, three metals explained the residual variation in the 16S TRFLP profiles: chromium (Cr), zinc (Zn), and copper (Cu) (Fig. 4). Each of the 14 metal pollutant variables was also tested manually for a correlation (i.e., ability to explain the residual variation) with 16S TRFLP profiles, using the Monte Carlo permutation test, and Cu (P = 0.075), Cr (P = 0.081), Ni (P = 0.066), and Li (P = 0.083) were marginally correlated in the presence/absence of 16S TRFLP data but significantly (P < 0.05) correlated in a relative abundance of 16S TRFLP data. Se, which was highly correlated with Cr (0.86), was marginally correlated in a relative abundance of 16S TRFLP data (P = 0.07). Mg, Rb, As, and K, which were not included in pCCA as representative metals, may also be related to 16S TRFLP profiles due to their high correlation with Cr (>0.85; see above).
![]() View larger version (17K): [in a new window] |
FIG. 4. Ordination diagram from pCCA of 16S TRFLP data, where elevations and DOC were used as covariables. Metals are represented by arrows, the lengths of which indicate the strength of correlation between the metal and the community data. Samples are labeled using their station names and denoted by triangles and circles. Open triangles pointing upward indicate high-elevation samples, filled triangles pointing downward indicate low-elevation samples, and filled circles indicate middle-elevation samples.
|
![]() View larger version (17K): [in a new window] |
FIG. 5. Ordination diagram from pCCA of nirS TRFLP data, where elevations and DOC were the covariables. Metals are represented by arrows, the lengths of which indicate the strength of correlation between the metal and the community data. Samples are labeled using their station names and denoted by triangles and circles. Open triangles pointing upward indicate high-elevation samples, filled triangles pointing downward indicate low-elevation samples, and filled circles indicate middle-elevation samples.
|
![]() View larger version (31K): [in a new window] |
FIG. 6. Total cell counts. Vertical dashed lines separate the groups of stations described in the legend to Fig. 1. Asterisks denote stations where total cell counts differed significantly across elevations. The bars represent the mean values, and the error bars represent the standard deviations (n = 3).
|
![]() View larger version (25K): [in a new window] |
FIG. 7. DEA measurements for three elevations (high, middle, and low) at each station. The inset graph excludes stations A and B to facilitate display. Vertical dashed lines separate the groups of stations described in the legend to Fig. 1. The data points represent the mean values, and the error bars represent the standard deviations (n = 3).
|
![]() View larger version (22K): [in a new window] |
FIG. 8. Relationship between DEA and DOC. Overall regression analysis across all three elevations yielded R2 values of 0.092 and P values of 0.072. The data points represent the mean values, and the error bars represent the standard deviations (n = 3).
|
When the variables of elevation and carbon were included in the statistical analyses, DEA was not significantly correlated with either the sodium concentrations or the metal levels. Therefore, to specifically evaluate the effects of metals on DEA, the variability residuals remaining after adjustment for the elevation-dependent DOC effects were correlated with metal concentrations. Four metals then appeared to be significantly negatively correlated with DEA: Pb (P = 0.0060), Al (P = 0.0068), Ni (P = 0.038), and Ba (P = 0.047). As, Fe, Ca, Rb, and Sr may also be negatively related to DEA due to their high positive correlations with Ni and/or Al (correlation coefficients of >0.85; see above).
|
|
|---|
The nirS communities were highly variable, more so than the 16S communities, and appeared to be subjected to multiple gradients or factors in the salt marsh. Similarly, a previous study also showed that functional genes with faster evolutionary rates, such as nirS, had higher levels of diversity than the more conserved genes encoding 16S rRNA (8). In this study, elevation appeared to account for the extensive variations in microbial community composition. This was similar to what was found in one of our previous studies (12) but different from what was found in another study regarding TRFLP profiles across several marshes (18). The contrast occurs because prevailing explanatory variables will change depending on the spatial and temporal context of a study. For example, in Córdova-Kreylos et al. (18), six marshes at different latitudes were analyzed in multiple years, and marsh, latitude, and year were the strongest factors influencing TRFLP profiles. On the other hand, Cao et al. (12) found that elevation was a significant factor influencing TRFLP profiles when only two marshes from one single year were studied. In a retrospective analysis of the latter study, performed here for comparison, the effect of elevation became even stronger when the differences between the two marshes (i.e., the "marsh" gradient) were statistically removed, i.e., elevation then explained 27% of the variation in 16S TRFLP profiles for one marsh (Fig. 1, stations A, B, C, and D), compared to 15% of the variation explained when both marshes were considered (12). Even when only one marsh is of concern (in this study), the differences between the subsystems (West Creek, the central road, and the east marsh) inside the single marsh are also potentially influential.
DEA was also highly variable within the marsh and was significantly correlated with DOC, but the correlations differed among elevations. The DEAs in this study were of magnitudes similar to those in prior reports, i.e., for a riparian zone affected by swine wastewater (3 to 1,660 µg N2O-N kg–1 dry soil h–1) (29), a rocky estuarine biofilm (140 to 840 µg N2-N kg–1 wet weight h–1) (48), and for salt marsh sediments amended with nitrate (1,500 µg N2O-N kg–1 dry soil h–1) (1). These, in turn, were all slightly higher than the levels for intertidal sediments (14 to 126 µg N2-N kg–1 wet weight h–1) (48). DEA is a measure of denitrifier abundance, denitrification pathway gene expression levels, and denitrification enzyme efficiencies. While the significant correlation between DEA and DOC in this study was consistent with the general knowledge (72), as most denitrifiers are heterotrophs (33) and often C limited (11, 15, 59), the correlations varied within the salt marsh, where multiple gradients or patchiness exist. For example, when samples from the three elevations, with their differing characteristics, were analyzed together, DEA was only marginally related to DOC (P = 0.072) and only 9.2% of DEA variability could be explained (Fig. 8). In contrast, when each elevation was analyzed separately, not only could a larger percentage of variability (on average, 50%) in DEA be explained by DOC, but also DEA in high elevations was much less responsive to increased C. This suggests an interaction between the elevation gradient and the DOC gradient. The interaction may result from factors affecting denitrification, such as competition of marsh plants for nitrogen (61), quality (61) and abundance (15) of DOC, and rhizosphere influences that can dampen the effects of environmental variations on microbes (6). These factors are part of the elevation gradient but perhaps vary inconsistently across the 12 stations. Also, because other factors, such as nitrate concentration, salinity, temperature, pH, oxygen supply, and even particle size (43), likely affect denitrification (17, 73), only a few or even none would appear dominant in a complex environment where the multiple gradients may simultaneously increase and/or decrease DEA. For example, nitrate concentration, sometimes reported as an important or even primary factor affecting sediment denitrification (54, 63), may exert minimal control on DEA in marsh sediments (this study and reference 48) because of the coexisting salinity effect (48).
Besides elevation and carbon availability, selected metals were also found to be significantly related to community composition and to DEA. While some metals affected several microbial measurements (e.g., Cu for both profiles and cell counts and Ni for both profiles and DEA), others (e.g., V) were specific to denitrifier community composition. Although effects of metals on total microbial community composition have been widely studied (for reviews, see references 24 and 35), and similar metals were observed influencing 16S TRFLP profiles in two previous studies involving the CSM (12, 18), few assessments have been made for denitrification (47). Further, similar metal amendments caused site-specific responses for denitrification rate (47), suggesting that variations in the "transducers" (76) of denitrification, i.e., denitrifier communities, could modulate effects of metals on the denitrification process. In studies of denitrifying fractions of whole communities, using extracted cells from a sandy loam, significant reductions in denitrification rates and in denitrifier growth rates were caused by heavy metals (Cd, Cu, and Zn), even when communities were adapted to the metals (26). In the present study, Cu, Ni, V, Se, and Cr were present in the eluates at concentrations that caused potential adverse effects on denitrifiers (Table 1) and affected denitrifier community composition. Although the concentration range of Cu (0 to 277 ng g–1 dry sediment) (Table 1) was 1,000-fold lower than that recently reported for estuarine sediments (47) and 2 orders of magnitude lower than the dose used in the previous study (80 mg kg–1 dry soil), only a tiny fraction (0.70 to 0.95 mg kg–1) of the latter was extractable and thus bioavailable (26). The concentration range of Ni (0.03 to 0.3 mg l–1) was similar to a level reportedly negatively affecting denitrification rates and diminishing nirS genes to nondetectable levels in river biofilms (0.5 mg l–1) (42). In this study, however, the metal concentrations in eluates were measured and thus represented not only the bioavailable portion of the sediment metals but also the fraction affecting long-term exposure to microbes (12). Indeed, long-term exposures to microbes at metal concentrations in ranges comparable to those used in this study have been shown to be toxic when free metal ion concentrations were considered (24).
Indirectly, the observed similar gradients (elevation, DOC, and metals) of DEA and denitrifying community composition could imply that community structure and function are coupled. However, a lack of direct correlation between denitrifier activity and community composition indicates an uncoupling, as in a previous study (66). Yet, denitrifier community composition has been shown to affect denitrification rate (14, 30, 65) and its response to environmental conditions, such as oxygen level and pH (14). The apparent inconsistency may be explained in several ways. First, the presence of functional genes (e.g., nirS in the current study and nosZ in reference 66) does not ensure transcription and hence DEA (60). This might be referred to as "methodological uncoupling." Second, since the ability to denitrify is widespread among diverse phylogenies, denitrifiers are functionally redundant (78), but meanwhile, enzymatic efficiency may vary greatly among species (22). Thus, shifts in community composition may not necessarily lead to changes in DEA of similar magnitudes. A plausible scenario would be that an environmental factor exerts effects on a highly abundant and/or efficient denitrifying group, and such effects are detected on both denitrifier community composition and DEA. Here, the heavy metal Ni (plus As and Rb) was significantly related to nirS TRFLP and also significantly affected DEA, which was similar to a prior report that Ni negatively impacted both nirS gene abundance and denitrification rate (42). Third, denitrifiers could be functionally complementary, as different subcommunities become active under different environmental conditions (60). For example, vanadium (V), although correlated with nirS community composition, did not correlate with DEA, perhaps because correlation between V and community composition may reflect the presence of vanadate-reducing denitrifiers (3, 4) causing a complementary shift (60) in community composition rather than inhibition of denitrification activity. Also, a positive correlation between V and DOC might have masked a correlation between V and DEA. Fourth, although both microbial community composition and environmental factors control denitrification activity (65), one or the other may dominate in different cases. Finally, although the cytochrome cd1-containing nitrite reductase (NirS) is generally more common (22), the copper-containing NirK protein may also be present and contribute to DEA in the current study. Perhaps Cu may induce a shift to NirK, thus influencing nirS TRFLP profiles while not significantly altering the overall DEA. Studies that include nirK TRFLP could be used to assess this possibility.
We thank Tom Young for the availability of the inductively coupled plasma mass spectrometer and Scott Olson, Laurie Van De Werfhorst, and Daniel Brooks for their assistance. This work was performed at the University of California Natural Reserve System Carpinteria Salt Marsh Reserve.
Published ahead of print on 31 October 2008. ![]()
Present address: Southern California Coastal Water Research Project, 3535 Harbor Blvd., Suite 110, Costa Mesa, CA 92626. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»