Previous Article | Next Article ![]()
Applied and Environmental Microbiology, June 2006, p. 4078-4087, Vol. 72, No. 6
0099-2240/06/$08.00+0 doi:10.1128/AEM.02969-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
R. Sean Norman,1,
Karen V. Pesce,1
Jerome J. Kukor,1,2,3 and
Gerben J. Zylstra1,4*
Biotechnology Center for Agriculture and the Environment, Rutgers University, 59 Dudley Road, New Brunswick, New Jersey 08901-8520,1 Department of Environmental Sciences, Rutgers University, 14 College Farm Road, New Brunswick, New Jersey 08901-8551,2 New Jersey Agricultural Experiment Station, 88 Lipman Drive, New Brunswick, New Jersey 08901,3 Department of Biochemistry and Microbiology, Rutgers University, 76 Lipman Drive, New Brunswick, New Jersey 08901-85254
Received 15 December 2005/ Accepted 3 April 2006
|
|
|---|
|
|
|---|
Aerobic degradation of lower molecular weight PAHs such as naphthalene (NP) by cultured microorganisms has been studied extensively (9). Bacterial degradation of PAHs under aerobic conditions begins with the addition of both molecules of molecular oxygen to the aromatic ring by a dioxygenase system. Aromatic ring dioxygenases are multicomponent enzymes which consist of an electron transport chain containing a ferredoxin and a reductase and a terminal dioxygenase (21). The dioxygenase is composed of two subunits. The alpha subunit is the catalytic component and contains two conserved regions: the [Fe2-S2] Rieske center and the mononuclear iron-containing catalytic domain. The Rieske cluster accepts electrons from the ferredoxin and passes them on to the mononuclear iron for catalysis (19, 42). The majority of information on PAH degradation pathways has come from studies on gram-negative bacteria, particularly the pseudomonads (7, 49). The best studied PAH dioxygenase is naphthalene dioxygenase from Pseudomonas putida NCIB 9816-4 (28, 35, 43), encoded by the NAH plasmid pDTG1 (14). These nah genes have been found in a wide variety of bacteria and geographic locations (1, 20, 36, 54, 55). Other more distantly related PAH degradation genes have also been described. Burkholderia sp. strain RP007, which was isolated from a PAH-contaminated site in New Zealand based on its ability to degrade phenanthrene (PH), contains a suite of PAH catabolic genes, the phn genes, which, while possessing activity similar to that of the nah genes, are only distantly related on the DNA and amino acid level (33). Competitive PCR studies by Laurie and Lloyd-Jones also showed that the Pseudomonas-type nah genes are not always dominant in the environment and that the phn-type genes can have a greater ecological significance than the nah-like genotype (34). Less information is available on PAH degradation by gram-positive bacteria, although recent reports have documented genetic and biochemical analysis of PAH degradation by Rhodococcus, Mycobacterium, Terrabacter, and Nocardioides (2, 8, 29, 41). Indeed, as more PAH-degrading bacteria have been isolated and characterized, it has become apparent that pseudomonads and the nah-like genes represent only a fraction of the PAH degradation picture. However, surveys of PAH degradation potential frequently rely on nah-based primers or probes to assess biodegradation potential in the environment (1, 52, 56).
Recently there have been a number of reports that describe primers or probes that can be used to measure PAH degradation potential in the environment (4, 17, 50, 55). Widada et al. (55) used a combination of isolation of PAH-degrading bacteria and PCR amplification of dioxygenase genes with a suite of degenerate primers to assess the diversity of PAH-degrading bacteria from different geographic locations. Even using a variety of primers, they failed to identify PAH dioxygenases genes from 7 of the 19 PAH-degrading isolates tested. Baldwin et al. (4) described a suite of aromatic oxygenase PCR primer pairs which amplify naphthalene, toluene, and biphenyl dioxygenase genes and their use in real-time PCR assessment of PAH biodegradation. While their study had the advantage of bypassing the culturing and isolation of the PAH degraders, which probably misses many of the important players in PAH biodegradation, they require the use of multiple primer sets to assess bioremediation at a site. Here we describe primers which target the conserved Rieske center of PAH dioxygenases and their use in monitoring PAH dioxygenase population shifts during degradation of naphthalene, phenanthrene, and pyrene (PY) in a series of enrichment experiments. The Rieske primers allow us to target only the dioxygenases which oxidize neutral aromatic hydrocarbons and represent a useful tool for assessing biodegradation potential in contaminated environments without cultivation and isolation of bacteria.
|
|
|---|
2 mm and consisted of 37% sand, 43% silt, 20% clay. It contained 2.77% organic carbon and had a pH of 4.95 in water. PAH residue analysis of background soil performed by Accutest Laboratories (Dayton, NJ) showed the following PAH content: 220 mg/kg phenanthrene; 457 mg/kg pyrene; 7 mg/kg acenaphthylene; 8 mg/kg dibenzo[a, h]anthracene; 11 mg/kg acenaphthene; 17 mg/kg benzo[g, h, i]perylene; 20 mg/kg indeno[1,2,3-cd]pyrene; 21 mg/kg fluorene; 31 mg/kg anthracene; 45 mg/kg benzo[a]pyrene; 68 mg/kg benzo[k]fluoranthene; 83 mg/kg benzo[b]fluoranthene; 108 mg/kg benzo[a]anthracene; 166 mg/kg chrysene; and 513 mg/kg fluoranthene.
Enrichment cultures.
Triplicate cultures were initiated in 150-ml serum bottles using 1 g of PAH-impacted soil and 25 ml of minimal medium (51) spiked with 200 mg/liter of NP, PH, or PY. Bottles were sealed using Teflon caps and shaken at 200 rpm for up to 250 h at 30°C. At each time point, samples were removed and separate triplicate cultures were used for cell enumeration and DNA extraction or PAH analysis.
Cell enumeration.
One milliliter of the enrichment culture was removed and fixed overnight at 4°C with 100 µl of 37% formaldehyde. Samples were vortexed to ensure even sample distribution, diluted 1:200 using sterile minimal medium, and placed with ice in an ultrasonic water bath for 5 min. Cells were stained with 5 µg/ml 4',6'-diamidino-2-phenylindole (DAPI [Sigma Aldrich, St. Louis, MO]) for 1 h in the dark at 4°C. Resulting solutions were filtered through a Millipore 0.2-µm prestained black 22-mm diameter filter, and a total of 9 grid fields per slide were counted using a Zeiss Axiovert 200 M epifluorescent microscope (45).
PAH analysis.
Dichloromethane (25 ml) was added to the enrichment cultures, and PAHs were extracted overnight at 200 rpm. Following extraction, the solvent layers were removed and dried with 4 g of anhydrous sodium sulfate. Extracts (1 ml) were analyzed on a Varian CP-3800 gas chromatograph with flame ionization detection using an RTX-5 column. The gas chromatograph program consisted of 6 min at 40°C followed by a 10°C/min increase to 300°C. The concentration of each PAH was calculated by comparison against individual PAH standard curves.
DNA extraction, PCR, and DGGE analysis.
At each time point, enrichment cultures were transferred to 50-ml Teflon tubes and centrifuged for 10 min at 27,000 x g. Total community DNA was extracted from the soil pellet using the UltraClean Soil DNA kit (MoBio Laboratories, Solana Beach, CA) and further purified using the Geneclean Spin kit (Q-BIOgene, Irvine, CA). For denaturing gradient gel electrophoresis (DGGE) analysis, a 193-bp sequence of the V3 region of the 16S rRNA genes was amplified using the primer set 341FGC and 534R as described by Muyzer et al. (39). PCR products were purified using QIAquick PCR purification columns (QIAGEN, Valencia, CA) and analyzed by denaturing gradient gel electrophoresis using a Dcode Universal Mutation Detection system (Bio-Rad Laboratories, Hercules, CA). Briefly, samples were run on 10% polyacrylamide gels with a denaturant gradient from 40% to 60%. Electrophoresis was carried out for 16 h at 70 V and 60°C. The gels were stained for 1 h with SYBR Green I and imaged using a gel documentation station (Kodak, Rochester, NY). Variability in the DGGE profiles for each treatment was determined by principal component analysis (PCA) based on the number of bands shared between profiles.
ARDRA.
For amplified ribosomal DNA restriction analysis (ARDRA), a 585-bp sequence of the 16S rRNA genes was amplified using the primer set 341F and 907R as described by Muyzer et al. (39). PCR products were gel purified, ligated into the pCR4-TOPO vector, and transformed into One Shot TOP10 chemically competent Escherichia coli following the manufacturer's protocol. One hundred randomly picked colonies per sample were each grown overnight at 37°C in Luria broth containing antibiotics. Plasmid DNA extraction was performed using the Qiaprep Spin Miniprep kit (QIAGEN, Valencia, CA), and the insert was amplified using primers 341F and 907R. Ten microliters of each amplicon was digested for 3 h at 37°C with 2.5 U of HaeIII, 2.5 U of RsaI, and 2.5 U of HinfI in 15-µl reactions. Restriction patterns were separated on 2% Metaphor agarose gels and imaged with SYBR Green I DNA stain. Each different restriction pattern was defined as an operational taxonomical unit (OTU). Amplicons from unique OTUs were sequenced, edited, and aligned using Lasergene sequence analysis software (DNAstar, Madison, WI). The distribution of OTUs in each treatment was determined and used to calculate the Shannon-Weaver index of diversity {H =
[ni · log(ni)], where ni represents the relative contribution of each OTU to the entire library}. Amplicons from unique OTUs were sequenced, confirmed, and hand aligned using Lasergene sequence analysis software (DNAStar, Madison, WI). Corrected sequences were screened against those in the GenBank database using Blastn.
Rieske primer design.
Amino acid sequences of the large subunit of dioxygenases targeting aromatic hydrocarbons were aligned using the ClustalW function in MegAlign (DNAStar, Madison, WI). The ProSite motif database describes the bacterial ring hydroxylating dioxygenase alpha-subunit signature as C-x-H-R-[GAR]-x (7, 8)-[GEKVI]-[NERAQ]-x (4, 5)-C-x-[FY]-H (40). However, close inspection of the aligned aromatic hydrocarbon dioxygenase amino acid sequences revealed differences in the Rieske sequences of the dioxygenases which attack nonpolar aromatic compounds, such as PAHs, polychlorinated biphenyls, benzene, toluene, and xylene, and those which attack polar aromatic compounds, such as benzoate, toluate, and phthalate. Primers were then designed based on this Rieske motif (Rieske_f, CRHRG; Rieske_r, CSYHGW) which allowed us to distinguish between dioxygenases targeting polar and nonpolar aromatic hydrocarbons. The sequences of the primers are the following: Rieske_f, TGYMGNCAYMGNGG; Rieske_r, CCANCCRTGRTANSWRCA. A second set of primers was designed substituting inosine residues for "N" residues to reduce the degeneracy of the primers. The Rieske primers are similar to those described by Cigolini et al. (12) and Kasuga et al. (27) but incorporate more degeneracies in order to amplify a wider diversity of PAH dioxygenase genes. The primers amplify a 78-bp PCR product. PCRs were prepared in either 25 or 50 µl containing 2 ng µl1 DNA, 1x PCR buffer, 2.5 mM MgCl2, 1 µM each forward and reverse primers, and 1 U Taq polymerase (Sigma, St. Louis, MO) and amplified with the following program: an initial denaturation step of 94°C for 5 min; 35 cycles of 94°C for 30 s, 48°C for 30 s, and 72°C for 30 s; and a final extension of 72°C for 5 min. The primers were tested on a series of bacterial strains which grow on benzene, toluene, naphthalene, phenanthrene, and biphenyl with successful results.
Amplification of dioxygenase genes from enrichment cultures.
DNA isolated from various time points of the PAH enrichments was screened for PAH dioxygenase genes using the Rieske primers. PAH dioxygenase genes were amplified as described above and cloned into pCR2.1 TOPO (Invitrogen, Carlsbad, CA). Colonies were picked, and plasmid DNA was isolated using the Concert 96 kit (Invitrogen). Rieske clones were sequenced with an ABI Prism 3100 Genetic Analyzer using Big Dye v3.1 chemistry (Applied Biosystems, Foster City, CA).
Phylogenetic analysis.
The nucleotide sequences were assembled, and the primer sequence was removed using SeqMan (DNAStar, Madison, WI). Rieske sequences which were 95% identical were assigned to the same clone family and aligned with reference sequences from GenBank using the ClustalW function in MegAlign (DNAStar). The phylogenetic tree was constructed by neighbor-joining analysis using Molecular Evolutionary Genetics Analysis software (32) with 1,000 bootstrap replicates.
Rieske clone library diversity analysis.
Accumulation curves were plotted for each library. The data were rarefied using the internet-based Rarefaction Calculator of Krebs and Brzustowski (31). The nonparametric richness estimates Chao1 and ACE (10, 11) and Shannon-Weaver diversity indices (H') (44) were calculated for each clone library.
Nucleotide sequence accession numbers.
The nucleotide sequences reported in this study were deposited in the GenBank database with the following accession numbers: DQ270422 to DQ270455 (16S rRNA clones) and DQ325359 to DQ325435 (Rieske clones).
|
|
|---|
![]() View larger version (10K): [in a new window] |
FIG. 1. Microbial growth ( ) and biodegradation of PAHs (O) over time. (A) Enrichment cultures with NP; (B) enrichment cultures with PH; (C) enrichment cultures with PY. Monitoring of NP cultures ended at 84 h due to rapid degradation of the compound in comparison to PH and PY. Data represent the averages and standard errors of triplicate data.
|
![]() View larger version (93K): [in a new window] |
FIG. 2. Denaturing gradient gel electrophoresis (DGGE) of PCR-amplified 16S rRNA gene fragments from PAH-amended enrichment cultures over time. The time of sampling (hours) is listed above the lanes. DGGE lanes are representative of triplicate data.
|
![]() View larger version (11K): [in a new window] |
FIG. 3. Ordination plots, by treatments and time, of bacterial communities generated by principle component analysis of bacterial species occurrence from 16S rRNA gene DGGE profiles. , unamended (UN); , naphthalene amended (NP); , phenanthrene amended (PH); , pyrene amended (PY).
|
|
View this table: [in a new window] |
TABLE 1. GenBank database sequences with the highest identity match to dominant OTUs and distribution of OTUs in PAH enrichment cultures
|
|
View this table: [in a new window] |
TABLE 2. Number of clones from unamended and PAH-amended enrichment culture 16S rRNA gene libraries representing different bacterial taxaa
|
Amplification of Rieske gene fragments from environmental samples.
The Rieske primers were utilized to assay the PAH dioxygenase gene complement from creosote-contaminated soil. Using total community DNA as the template, the Rieske primers amplified a single 78-bp PCR product (data not shown). The unamended soil library contained 25 Rieske clone families and was dominated by one family, S1, which represented 30% of the clones sequenced (Fig. 4). The closest match for this family is naphthalene dioxygenase from Comamonas testosteroni GZ42 (Blastn and BlastP identity of 100%) (24). Three other clone families, S2, S3, and S4, represented 36% of the clone library. Clone family S2 is most closely related to DxnA1, a dioxin dioxygenase from Sphingomonas sp. strain RW1 (3) (Table 3). Clone family S3 is most closely related to naphthalene dioxygenase from Polaromonas naphthalenivorans and P. putida G7 (Table 3) (49). Clone family S4, however, which represents 9% of the soil Rieske library, does not exhibit significant similarity to any dioxygenase sequences currently in the GenBank database. Indeed, 18% of the soil clone library sequences do not exhibit significant similarity to any dioxygenase sequences currently in GenBank, and a further 13% have only low matches (<50% identity) to uncharacterized dioxygenase sequences from GenBank. Some of these Rieske sequences form clusters with clone sequences from the enrichment libraries, while others are unique to the soil library (Fig. 4). Analysis of the relationship of the Rieske sequences to one another shows that clone families S5, S13, and S25 form a unique group (labeled A) and comprise 8% of the soil library. A second cluster of sequences with no known relatives in GenBank (labeled B) includes S6, S7, S21, and S23 in addition to PH6 and PH16. The remaining soil library sequences are distributed among 14 clone families and represent a wide diversity of Rieske sequences (Fig. 4).
![]() View larger version (13K): [in a new window] |
FIG. 4. Phylogenetic distribution of the Rieske gene fragment families from the soil (S) and enrichment libraries. The dendrogram was constructed from a ClustalW alignment of the Rieske sequences by neighbor-joining analysis using Mega 3.0. Nodes supported by bootstrap values greater than 50 are indicated with a filled black circle. Major clone families are in boldface and have the number of clones observed indicated in parentheses. Reference sequences from GenBank include the accession number. The scale bar represents substitutions per nucleotide.
|
|
View this table: [in a new window] |
TABLE 3. Best matches in the GenBank database for the dominant Rieske clone families in the naphthalene, phenanthrene, and pyrene enrichments and the unamended creosote-contaminated soil libraries
|
|
View this table: [in a new window] |
TABLE 4. Diversity analysis of the Rieske gene fragment librariesa
|
The PY enrichment culture clone library contained 15 clone families. The library was dominated by clone family PY1, accounting for 38% of the clones sequenced (Fig. 4). Clone family PY1 is most closely related to an alkylbenzene dioxygenase from Rhodococcus sp. strain DK17 (Table 3) (30). A second clone family, PY3, which accounts for 22% of the enrichment library, is also most closely related to alkylbenzene dioxygenase from Rhodococcus sp. strain DK17. Clone family PY2 is most closely related to an uncharacterized aromatic dioxygenase from Sphingomonas aromaticovorans (Table 3). The remaining 12 PY enrichment Rieske families occurred rarely (Fig. 4). There was significant overlap between the PH and PY libraries, with the dominant families in both libraries being most closely related to alkylbenzene dioxygenase from Rhodococcus sp. strain DK17. Indeed, clone families PY1 and PH1 have identical nucleotide sequences, as do PY2 and PH4 as well as PY3 and PH2.
Diversity and species richness of the soil and enrichment libraries.
The Shannon-Weaver diversity index and the Chao1 nonparametric species richness estimators were calculated for each Rieske library. The Shannon-Weaver index indicates that the soil enrichment library was the most diverse, followed, in order, by the PH, PY, and NP libraries (Table 4). This order is in agreement with the 16S-based bacterial community diversity reported above. The Chao1 estimate of species richness also indicates a difference in the estimated species richness of the soil and enrichment libraries. Enrichment on the different model PAHs led to selection of different dioxygenase gene populations. Species accumulation curves were plotted for each library using the internet-based rarefaction calculator (31) in order to assess how well our sequencing effort had sampled our libraries. Even though more than 100 clones were sequenced from each library, none of the accumulation curves reached an asymptote (Fig. 5). This is an indication that further sequencing would likely yield new Rieske gene fragment families. However, even our limited sampling effort indicates clear differences in the richness of the four libraries, with the NP and PY libraries being significantly less diverse than the soil and PH libraries. Furthermore, our sampling clearly identified dominant Rieske gene fragment families in each enrichment culture.
![]() View larger version (13K): [in a new window] |
FIG. 5. Rieske clone family accumulation curves showing the number of clones sequenced versus the number of Rieske families observed. Soil ( ), naphthalene ( ), phenanthrene ( ), and pyrene ( ) samples are shown.
|
|
|
|---|
In this study, we isolated DNA from a PAH-contaminated site and used the Rieske PCR primers to screen for the presence of PAH dioxygenase genes. The site was contaminated with a variety of high-molecular-weight PAHs. Enrichment of the bacterial population in the presence of model two, three, and four ring PAHs led to changes in both the microbial community structure and the dioxygenase gene profile. Our Rieske gene fragment survey revealed the presence of a wide diversity of dioxygenase genes, ranging from the well characterized naphthalene dioxygenase family to less well studied genes such as those of the dioxin dioxygenase family. A significant number of the clones sequenced had no relatives in the GenBank database.
While the unamended soil library was quite diverse, enrichment on NP led to the domination of the Rieske library by one nahAc-like Rieske sequence. It is interesting that the dominant nahAc clone from the soil library is not the dominant nucleotide sequence observed in the NP enrichment library. Sequences identical to P. putida G7 nahAc and C. testosteroni nagAc accounted for 41% of the soil library but were barely detected in the NP enrichment library. Instead, the NP enrichment library was dominated by a Rieske fragment most closely related to P. putida NCIB 9816-4 nahAc. The NP-enriched cultures exhibited a 24-h lag time before the onset of NP degradation followed by near-complete degradation by 46 h. Analysis of the initial PAH levels in the contaminated soil used in our study did not show the presence of NP, so the lag time was likely due to the absence of chronic exposure to NP, thus requiring time for an increase in the number or activity of NP-degrading populations. This is supported by molecular analysis of the NP-amended cultures showing that the microbial community remained similar to unamended cultures until after 24 h, at which point there was a shift in community structure corresponding to increased NP degradation and the development of an almost monoculture-like microbial community consisting predominantly of a single Pseudomonas fluorescens-like phylotype. It is interesting that the NP enrichment was dominated so quickly by one organism and one dioxygenase gene. Stach and Burns (50) reported a similar decrease in diversity in batch cultures compared to biofilm cultures in a study that enriched with NP and PH.
A frequent criticism of enrichment cultures is that they select only for the fastest growing organisms under the conditions utilized and do not reflect what occurs in nature (18, 23). While enrichment on NP led to a profound decrease in diversity, the PH enrichment library was more complex and more closely resembled the unamended soil library, possibly because the original soil contained relatively high PH concentrations. The prior adaptation of the microbial community to the presence of PH is supported by molecular analysis of the community structure. DGGE fingerprints for PH-amended cultures did show a slight reduction in band number; however, the overall community fingerprint changed in a way similar to that of the unamended cultures. This is not surprising, given that the three ring and higher PAHs in the original soil materials most likely caused a shift in the unamended samples. Clone libraries and phylogenetic analysis show similar enrichment for OTU 1 and OTU 2 (phylotypes closely affiliated with Burkholderia and Sphingomonas) in the PH-amended and unamended cultures. Moreover, while small numbers of OTUs 22, 31, and 32 (phylotypes affiliated with Nitrosococcus, Stenotrophomonas, and Massilia) were observed in the unamended cultures, they were enriched to high numbers in PH-amended cultures. Members of these genera have been demonstrated to degrade a range of PAHs, suggesting they may play a role in PH degradation in our cultures (5, 6, 22, 55, 57).
nahAc-like genes were rare in the PH enrichment. In contrast to the NP enrichment library, which was dominated by one Rieske sequence, the PH enrichment library contained a number of equally dominant Rieske sequences and remained quite diverse. This contrasts with the findings of Stach and Burns, who reported a decrease in dioxygenase gene diversity in both PH and NP enrichment and biofilm communities (50). The Rieske data also contrast with the microbial community profiling results in that the unamended and PH enrichments share some dominant bacterial phylotypes but have no dominant Rieske families in common. It is possible that the bacterial community analysis detects not only bacteria growing on the added PH but also on by-products from PH degradation. The use of the Rieske primers has the advantage of allowing us to target bacteria possessing the ability to catalyze the first step in PH biodegradation.
Molecular analysis of the adaptation of the microbial community to growth on PY shows that after PY degradation, the community was enriched in three different phylotypes. As in the unamended and PH-amended cultures, OTUs 1 and 2 (Burkholderia and Sphingomonas) were observed in PY-enriched cultures. However, we observed a strong selection for OTU 31 (Stenotrophomonas), resulting in a low Shannon-Weaver diversity index. Previous studies have documented the presence of members of the genus Stenotrophomonas in PAH-contaminated soils, and individual bacteria isolated from these sites have been shown to mineralize PY (6, 25, 58). The PH and PY enrichment libraries had a number of Rieske clone families in common, suggesting that these enzymes may be involved in both PH and PY degradation. Both the PH and PY libraries contained alkylbenzene dioxygenase-like sequences and an uncharacterized aromatic dioxygenase-like sequence as two of their most dominant clone families. Both libraries also contained a number of sequences that have no close relatives in the GenBank database, indicating that there remains a wealth of untapped dioxygenase sequence diversity and potentially functional diversity in contaminated environments.
The goal of our current studies is to design a better method of predicting the overall biodegradative and bioremediative potential of contaminated sites. Our approach used traditional laboratory biodegradation studies combined with molecular ecological techniques to characterize shifts in both microbial community and functional gene profiles during the biodegradation of a range of model PAHs. Overall, our study suggests that contaminated environments can harbor a wide diversity of dioxygenase genes, allowing for succession of different dioxygenase gene populations in response to exposure to different PAHs. The combination of the Rieske primers and bacterial community profiling represents a powerful tool for both assessing bioremediation potential in the environment and for the discovery of novel dioxygenase genes.
We thank Laurie Seliger for excellent DNA sequencing assistance.
These authors contributed equally to this paper. ![]()
|
|
|---|
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»