AEM
Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental material
Right arrow Other Versions of this Article:
AEM.01254-07v1
73/19/6224    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by de Cárcer, D. A.
Right arrow Articles by Rivilla, R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by de Cárcer, D. A.
Right arrow Articles by Rivilla, R.
Agricola
Right arrow Articles by de Cárcer, D. A.
Right arrow Articles by Rivilla, R.

 Previous Article  |  Next Article 

Applied and Environmental Microbiology, October 2007, p. 6224-6232, Vol. 73, No. 19
0099-2240/07/$08.00+0     doi:10.1128/AEM.01254-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Changes in Bacterial Populations and in Biphenyl Dioxygenase Gene Diversity in a Polychlorinated Biphenyl-Polluted Soil after Introduction of Willow Trees for Rhizoremediation{triangledown} ,{dagger}

Daniel Aguirre de Cárcer,1 Marta Martín,1 Ulrich Karlson,2 and Rafael Rivilla1*

Departamento de Biología, Universidad Autónoma de Madrid, 28049 Madrid, Spain,1 National Environmental Research Institute, University of Aarhus, Roskilde, Denmark2

Received 6 June 2007/ Accepted 2 August 2007


   ABSTRACT
Top
 ABSTRACT
INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The aim of this study was to analyze the structural and functional changes occurring in a polychlorinated-biphenyl (PCB)-contaminated soil ecosystem after the introduction of a suitable host plant for rhizoremediation (Salix viminalis). We have studied the populations and phylogenetic distribution of key bacterial groups (Alpha- and Betaproteobacteria, Acidobacteria, and Actinobacteria) and the genes encoding iron-sulfur protein {alpha} (ISP{alpha}) subunits of the toluene/biphenyl dioxygenases in soil and rhizosphere by screening gene libraries using temperature gradient gel electrophoresis. The results, based on the analysis of 415 clones grouped into 133 operational taxonomic units that were sequence analyzed (>128 kbp), show that the rhizospheric bacterial community which evolved from the native soil community during the development of the root system was distinct from the soil community for all groups studied except for the Actinobacteria. Proteobacteria were enriched in the rhizosphere and dominated both in rhizosphere and soil. There was a higher than expected abundance of Betaproteobacteria in the native and in the planted PCB-polluted soil. The ISP{alpha} sequences retrieved indicate a high degree of catabolic and phylogenetic diversity. Many sequences clustered with biphenyl dioxygenase sequences from gram-negative bacteria. A distinct cluster that was composed of sequences from this study, some previously described environmental sequences, and a putative ISP{alpha} from Sphingomonas wittichii RW1 seems to contain greater diversity than the presently recognized toluene/biphenyl dioxygenase subfamily. Moreover, the rhizosphere selected for two ISP{alpha} sequences that accounted for almost 60% of the gene library and were very similar to sequences harbored by Pseudomonas species.


   INTRODUCTION
Top
 ABSTRACT
 INTRODUCTION
MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Soil contamination by organic compounds is a serious problem in most industrialized countries. Among pollutants, polychlorinated biphenyls (PCBs) are of special importance, and it is estimated that more than 300 million kg has been released to the environment (20). Remediation of PCB-contaminated soil has traditionally been carried out by incineration or burial of the soil in secure landfills. As a less-expensive removal strategy, rhizoremediation (53), the use of dual plant-microorganism systems in which the plant provides nutrients, support, and a greater availability of the substrate and the microorganisms drive the enzymatic remediation, has been proposed (26, 47).

The soil compartment directly under the influence of living roots is designated the rhizosphere (25). It is a niche where complex microbial communities are supported by nutrients released by root exudates, mucilage, and decaying root cells (7), which induces a characteristic change in the distribution of the microorganisms associated with plant roots in comparison with their distribution in the bulk soil, designated the rhizosphere effect. Thus, in order to exploit the potential of the emerging rhizoremediation technology, a better understanding of the diverse microbial populations and catabolic gene polymorphisms that develop in the rhizosphere is required (46). Salix sp., which produces salicylic acid and related compounds that induce the degradation of many xenobiotic molecules (such as polycyclic aromatic hydrocarbons) (41), is a good candidate for rhizoremediation of PCB-contaminated soils and has already been used in PCB rhizoremediation studies (1).

The aerobic degradation of PCBs is initiated by enzymes of the toluene/biphenyl subfamily of the Rieske nonheme iron oxygenases (13). These enzymes are multicomponent complexes composed of a terminal oxygenase iron-sulfur protein (ISP) and various electron transport proteins (4). The ISPs are heteromultimers, composed of a large ({alpha}) and a small (ß) subunit. The {alpha} subunit contains the substrate-binding site (4, 11) and thus is responsible for substrate specificity (13). Recently, various culture-independent studies have shown that the ISP{alpha} sequences obtained from cultured strains (from which most available information on ISP{alpha} structure and mechanism derives) are not likely to represent the functional gene diversity in the environment, as distinct and numerically dominant (putative) functional genes were obtained from screening various different environments (43, 44, 51). The use of PCR-based genetic-profiling techniques for the study of ISP{alpha} diversity in the environment (23, 43, 52) can be employed for monitoring the predominant ISP{alpha}'s structure and diversity and for detecting enzymes selected under the existing environmental conditions, which should prove highly valuable in future bioremediation strategies.

There are few studies analyzing the bacterial community in PCB-polluted soils, either by culture-dependent (27) or culture-independent methods (31, 32), and to our knowledge, no culture-independent studies of the rhizosphere of plants growing in PCB-polluted soil have been reported. In a previous study (1), we assessed the impact caused by the introduction of genetically engineered bacteria designed for rhizoremediation on the native microflora of a naturally PCB-polluted soil. The aim of this study was to evaluate the structural and functional changes occurring in the bacterial populations of a PCB-contaminated soil ecosystem after the introduction of a suitable host plant for rhizoremediation by screening environmentally derived gene libraries by temperature gradient gel electrophoresis (TGGE) and sequencing.


   MATERIALS AND METHODS
Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
RESULTS
 DISCUSSION
 REFERENCES
 
Experimental setup.
Willows (Salix viminalis x schwerinii, variety Björn), were pregrown from 20-cm-long cuttings in tap water for three weeks. Three willow plants were planted per galvanized iron pot in 1,900 g of previously homogenized and sieved (2-mm mesh) PCB-contaminated soil. Nonvegetated pots were included as controls. The soil, obtained from Lhenice (Czech Republic), contained Delor 103 (similar to Aroclor 1242), resulting from historic industrial activity, at a concentration of 182 ± 22 mg PCB/kg dry soil (1). All pots were incubated under controlled conditions (1). After 6 months, the plants were removed from three replicate pots, and the roots from each pot were shaken carefully to remove the bulk soil. The soil still adhering to the roots was defined as "rhizosphere" soil from the corresponding replicate. "Soil" samples were taken from three replicate control pots by coring the potted soil using polyvinyl chloride tubes (28-mm inside diameter). All samples were stored at –20 C until DNA extraction.

DNA extraction and clone library construction.
Total DNA was extracted from 0.5-g (wet weight) samples of each of the three rhizosphere and bulk soil fractions following a protocol described by Porteus et al. (35) and then purified using a Microcon YM-100 filtering device (Millipore, Bedford, MA). One µl of pooled DNA (ca. 15 ng) was used to amplify the 16S rRNA genes with group-specific primers targeting Alphaproteobacteria ({alpha}-U203 [15] and L1494 [50]), Betaproteobacteria (F948ß [15] and L1494), Acidobacteria (31F [3] and L1494), and Actinobacteria (F243HGC and R1378 [19]), and part of the gene encoding the ISP{alpha} subunit of the Rieske nonheme oxygenases of the toluene/biphenyl subfamily was targeted by employing primer pair bphAf668-3/bphAr1153-2 (51). The annealing temperatures employed for each primer pair were 56°C, 61°C, 42°C, 63°C, and 58°C, respectively.

The PCR mixture consisted of 2.5 µl of reaction buffer, 5 pmol of primers, 3.75 mM MgCl2, 0.2 mM deoxynucleoside triphosphates, 5 mg bovine serum albumin, 4% dimethyl sulfoxide, and 1 U Tth DNA polymerase (Biotools, Madrid, Spain) in a final volume of 25 µl. The PCR program was as follows: 5 min of denaturation at 94°C, followed by 25 cycles of 1 min at 94°C, 1 min for primer annealing, 2 min at 72°C for primer extension, and a final cycle at 72°C for 10 min. The products of two consecutive PCRs were then pooled and purified through extraction from agarose gels prior to cloning on pDRIVE vectors (QIAGEN, Germany). The resulting plasmids were transformed in competent Escherichia coli DH5{alpha} cells and checked for inserts of the correct size.

In order to infer the percentage of the total bacteria included in the study, the 16S rRNA genes from the Eubacteria present were amplified as well, using universal primers (27f [50] and L1494), and then cloned on pDRIVE and transformed in E. coli DH5{alpha} cells. The clones presenting an insert of the correct size were used as templates for group-specific PCRs.

TGGE.
One-microliter amounts of 1:10 dilutions of the amplified products were used for a second round of nested PCR using high-pressure liquid chromatography-purified primers (48). For 16S rRNA gene amplification, primer pair F984GC (19) and R1378 was used, resulting in ~475-bp fragments suitable for TGGE analysis. For the nested amplification of the ISP{alpha} genes, primer bphAf668-3 was modified with a GC clamp (5'-CGC CCG GGG CGC GCC CCG GGC GGG GCG GGG GCA CGG GGG G-3') at its 5' end.

For sequence-dependent separation of the PCR products derived from the gene libraries, a TGGE system (Biometra, Göttingen, Germany) was used as specified by the manufacturer. Six-percent polyacrylamide gels were polymerized in 1x Tris-acetate-EDTA buffer, 8 M urea, 20% deionized formamide, and 2% glycerol. Amounts of 5 µl of the PCR samples were resolved by electrophoresis at 130 V for 16 h through a temperature gradient of 44.5°C to 56°C for the 16S rRNA genes and 41°C to 56°C for the ISP{alpha} genes. Reference patterns consisting of various group-specific sequences were included in the gels to facilitate the analysis of the gels. The DNA bands were visualized by a routine silver-staining protocol (18). The gels were then digitally documented and analyzed (Quantity One; Bio-Rad, Hercules, CA).

Sequencing, phylogenetic analysis, and statistics.
Clones were grouped according to their corrected electrophoretic mobilities, reported as operational taxonomic units (OTUs), and named with the prefix Lhap, Lhbe, Lhad, Lhac, or LhISP depending on the sequence type (Alphaproteobacteria, Betaproteobacteria, Acidobacteria, Actinobacteria, and ISP{alpha} sequences, respectively; Lh stands for Lhenice, the origin of the soil). Their frequencies were used as an indicator of abundance. Representative clones were grown, their plasmidic DNA was extracted (RapidPURE kit; QBIOgene, Irvine, CA), and the insert was fully sequenced (Parque Científico de Madrid, Madrid, Spain). The 16S sequences obtained were submitted to the SEQUENCE MATCH program of the Ribosomal Database Project (RDP-II) (6), and the ISP{alpha} sequences were submitted to BLAST to find closely related sequences. The sequence identities between clones and the next related sequences were calculated using BIOEDIT (17). The raw sequences obtained were edited and aligned with ClustalW (45). The online programs CHECK_CHIMERA (7) and Bellerophon (21) were used to rule out the presence of chimeric sequences (hybrid artifact sequences that may arise from the PCR amplification of very similar sequences) (10). The raw data (sequences and frequencies) were used to infer possible differences, at group level, between bulk soil and rhizospheric communities. The diversities of clone libraries were calculated using the Shannon index (40). Richness and diversity analysis was performed using the DOTUR software (38). Maximum-likelihood phylogenetic trees were constructed (for the 16S sequences) using a quartet-puzzling algorithm implemented in the Tree-Puzzle software (nucleotide substitution model HKY, including gamma correction and with the molecular clock assumption enforced) (42) and were visualized with TreeView (33). For each group analyzed, a set of 16S rRNA gene sequences from the closest cultured bacteria were downloaded from the RDP-II database (6) in order to obtain a comprehensive view of the phylogeny of the sequences. The distance matrix and the {alpha} value for the gamma correction obtained were used in the subsequent analysis. MEGA software (http://www.megasoftware.net/) was used for the phylogenetic reconstruction of the ISP{alpha} sequences by using the neighbor-joining method with the JTT distance model and pairwise deletion of gaps/missing data. A consensus tree was inferred from a total of 1,000 bootstrap trees. The FST test (29) was used to compare the genetic diversity within each community to the total genetic diversity of the communities combined, using the equation FST = ({theta}t-{theta}w)/{theta}t, where {theta}t is the genetic diversity for all samples (based on the mean number of pairwise differences) and {theta}w is the genetic diversity in each community. The population differentiation using FST was calculated using the ARLEQUIN program (39); statistical significance was evaluated by randomly assigning sequences to populations and calculating the FST for 3,034 permutations. Additional estimators were also obtained with this software. Each pair of libraries was compared first in an X/Y fashion and then in the reverse Y/X using the {int}-LIBSHUFF computer program (37) which calculates the integral form of the Cramér-von Mises statistic of the formula for coverage (16) and compares it by a Monte Carlo test procedure. The Library compare program of the RDP-II was used as well. It utilizes the RDP naïve Bayesian classifier (confidence threshold, 95%) to provide classification of library sequences into the bacterial taxonomy and estimates the probability of observing the difference in a given taxon using a statistical test.

Community-level catabolic profiles.
The bacterial fraction from each replicate soil and rhizosphere sample was harvested as previously described (1). Each suspension of bacterial cells was diluted 1:10 and used to inoculate ECOlog plates (5) at 135 µl per well; the plates were then incubated at 28°C in darkness. The optical density (OD) at 595 nm was measured with a microtiter plate reader (Tecan, Zurich, Switzerland) when a maximum of wells showed activity (66 h) (14). The raw OD values were corrected for the background color in the control well without a carbon source. To account for slightly differing initial cell densities, the corrected ODs were divided by the average well color of the plate, giving the standardized OD; thus, standardized patterns rather than absolute values were compared (12). Values smaller than 0.006 were counted as 0 (no catabolic activity) to reduce noise. The resulting 31 variables (corresponding to each carbon source) were reduced to three new principal components (SPSS 12.0) which explained more than 80% of the variance (total eigenvalue sum) using the covariance of the variables. In order to test for differences between the communities, a multivariate analysis was performed (general linear model function; SPSS 12.0) using the new principal components as variables (14).

Nucleotide sequence accession numbers.
The sequence data have been submitted to the GenBank database under accession numbers DQ648899 to DQ649002 and EF565826 to EF565853.


   RESULTS
Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
DISCUSSION
 REFERENCES
 
Composition of the eubacterial populations in soil and rhizosphere.
Soil- and rhizosphere-derived eubacterial 16S rRNA gene libraries were constructed and screened for clones belonging to Alpha- and Betaproteobacteria, Acidobacteria, and Actinobacteria by amplifying random clones with each of the group-specific primer pairs. Forty-eight clones from the soil library and 43 clones from the rhizosphere library were screened. Table 1 shows the percentage of clones from each library that was amplified with each of the primer pairs. No clones were amplified with more than one primer pair, indicating that the group specific-primers were targeting nonoverlapping groups. The results also show that the targeted groups represented a significant fraction of the eubacterial populations studied (Table 1). The soil and rhizosphere eubacterial communities appear to be dominated by Proteobacteria, representing 44% and 52% of the soil and rhizosphere amplified clones, respectively. Both Alpha- and Betaproteobacteria were enriched in the rhizosphere compared to their presence in soil.


View this table:
[in this window]
[in a new window]

  		TABLE 1. Percentage of the total eubacterial clones screened that were positively amplified with the different group-specific primer sets

 
Analysis of the group-specific 16S rRNA gene libraries.
Group-specific 16S rRNA gene libraries were constructed using soil and rhizosphere DNA samples. Random clones were subjected to nested PCR and TGGE, and OTUs were defined by gel mobility within each library (see Materials and Methods). The number of times that an OTU appeared was scored, and one clone corresponding to each OTU was fully sequenced. A total of 320 clones representing 105 OTUs were analyzed, corresponding to a total of more than 115,000 bp sequenced. Table 2 shows the observed diversity for each library based on OTU abundances. The Shannon index of diversity was higher in the rhizosphere for Alpha- and Betaproteobacteria and higher in the soil for the Actinobacteria, although these differences were not statistically significant. However, a significant difference was found for the Acidobacteria community, which was more diverse in the rhizosphere. Table 2 also shows the genetic diversity, inferred by the nucleotide diversity and {theta}({pi}) (average sequence divergence) values, which were very similar in the soil and rhizosphere for the Acidobacteria and the Actinobacteria and slightly higher in the rhizosphere for the Alpha- and Betaproteobacteria. The closest relative (most similar sequence) found in the databases for each OTU, as well as the number of times it was scored, is shown in Table 3.


View this table:
[in this window]
[in a new window]

  		TABLE 2. Numbers and diversity indices for gene libraries derived from rhizosphere or soil

 

View this table:
[in this window]
[in a new window]

  		TABLE 3. Bacterial sequences from the RDP-II database with the highest similarity to each OTU (closest relative) and distribution of OTUs in the rhizosphere or soil samples

 
Alphaproteobacteria.
In the case of the Alphaproteobacteria, their closest relatives (similarities of between 83.6 and 99.4%) were usually from uncultured bacteria, although close similarities to several cultured species from the Rhizobiales order were obtained. A wide phylogenetic distribution was observed for this set of sequences (see Fig. S1A in the supplemental material). Many sequences clustered as Rhizobiales, but other orders, including Rhodobacterales, Caulobacterales, and Rhodospiralles, were also represented in the set. Several OTUs showed deep branching in the phylogenetic tree. The group formed by OTUs Lhap1 and Lhap2 shows a close relationship to bacteria that form a novel unnamed order of Alphaproteobacteria described previously (36) and probably belongs to this order. The groups formed by Lhap25 and Lhap8, as well as Lhap15 and Lhap16, do not show a close relationship to any described taxon and may represent novel lineages within Alphaproteobacteria. Sequences from the Rhizobiaceae family were only found in the rhizosphere library, and a slight rise in the population of Rhodobacteraceae was detected in the rhizosphere (P < 0.0588) by using the RDP-II library compare test.

Betaproteobacteria.
Most betaproteobacterial sequences obtained were found to be very similar (94.8 to 100%) to sequences obtained from uncultured organisms already present in the databases. The phylogenetic distribution of the betaproteobacterial sequences obtained (see Fig. S1B in the supplemental material) was clearly dominated by the Burkholderiales order, including members from the Oxalobacteraceae, Alcaligenaceae, and Comamonadaceae families, in both the soil and the rhizosphere libraries. A significant rise in the number of sequences belonging to the Comamonadaceae family (P < 0.032) was observed in the rhizosphere by using the RDP-II library compare test. This result is also clearly observable in the phylogenetic tree (see Fig. S1B in the supplemental material).

Acidobacteria.
Almost all Acidobacteria sequences obtained were most similar to uncultured species (93 to 98.8%). Most sequences, according to the grouping by Barns et al. (3), were from groups III and V (see Fig. S2A in the supplemental material). Group II sequences were detected only in the rhizosphere. No sequences were detected for the other groups. One of the OTUs (Lhad13), which was found both in the soil and rhizosphere libraries, did not belong to any of these groups, as judged from its deep branching. This OTU, together with the previously reported sequence from a sponge symbiont (GenBank accession no. AY703463), may represent a novel group within Acidobacteria. No classification or comparison was possible by using the RDP-II for the Acidobacteria sequences since there is yet no robust phylogeny for this group.

Actinobacteria.
The Actinobacteria sequences were found to be very similar (95.4 to 99.6%) to their closest relatives, which were often cultured species. Surprisingly, two main clusters were found, according to their phylogenetic distribution (see Fig. S2B in the supplemental material), the Verrucomicrobiales and the Actinobacteria clusters, which indicates a lack of true specificity for the Actinobacteria-specific primer pair under the PCR parameters employed. Most actinobacterial sequences were from the Actimomycetales order, including representatives of various families such as Sporichtyaceae, Microbacteriaceae, Mycobacteriaceae, and Microsphaeraceae, many sequences from the genus Streptomyces, and one cluster of unclassified Actinobacteria. No differences between soil- and rhizosphere-derived libraries were observed in the phylogenetic distribution of the sequences. Moreover, no significant difference between libraries was detected according to the RDP-II library compare test.

Statistical comparisons of soil- versus rhizosphere-derived 16S rRNA gene libraries.
A statistical test of differences between soil- and rhizosphere-derived group-specific libraries (Table 4) was performed by comparing the genetic diversity within each community to the total genetic diversity of the communities combined (FST test). The test results indicate significant differences between gene libraries for all groups except for the Actinobacteria. The same test results were obtained when analyzing differences between homologous and heterologous coverage curves by a Cramér-von Mises-type statistic and comparing them by a Monte Carlo test procedure ({int}-LIBSHUFF).


View this table:
[in this window]
[in a new window]

  		TABLE 4. Statistical analyses of the difference between rhizosphere and soil bacterial communities

 
Analysis of community-level catabolic profiles.
No statistical difference was found between the Shannon diversity indices of the soil and rhizospheric communities based on the use of 31 ecologically relevant carbon sources (5) (3.00 ± 0.05 [mean ± standard deviation] for soil versus 2.92 ± 0.07 for rhizosphere). Nevertheless, a significant statistical difference between the soil and rhizospheric catabolic potentials was found by analyzing the catabolic profiles (P < 0.02).

ISP{alpha} gene diversity.
A PCR-TGGE approach was employed for the detection and discrimination of the most abundant ISP{alpha} genes of the toluene/biphenyl subfamily of the aromatic ring-hydroxylating dioxygenases present in soil and rhizosphere samples from a PCB-polluted soil. Forty-eight soil-derived clones and 47 rhizosphere-derived clones were analyzed, yielding 28 different OTUs. As expected, all sequences obtained were designated by BLAST searches as Rieske nonheme iron dioxygenases. The inferred peptide sequence of the ISP{alpha} sequences obtained ranged between 37% (LhISP18) and 100% (LhISP1) identity to previously described ISP{alpha} sequences (Table 3). Most of them contained residues considered characteristic of these enzymes (P. fluorescens IP01 cumene dioxygenase [CumDO] numbering) (8): His234 and His240, involved in the coordination of mononuclear iron at the active site (8, 11); Asp231, with a role in the electron transfer bridge (8, 11); and residues Glu226, Gln227, and Tyr233 (or the equivalent Glu/Asp, Gln/Asn, and Tyr/His), which are thought likely to be involved in electron transfer by connecting the catalytic center to the Rieske center of an adjacent subunit (13, 34).

On the other hand, great variation was observed among the 14 residues proposed to define the inner surface of the substrate-binding pocket in IP01 CumDO (51) (Fig. 1). The phylogenetic analysis of the peptide sequences deduced from the clone sequences (see Fig. S3 in the supplemental material) shows that all of them were members of the toluene/biphenyl subfamily. Many sequences clustered with the gram-negative biphenyl dioxygenase group, while none of them grouped with members of the isopropyl/benzene, toluene/benzene, and gram-positive biphenyl dioxygenase clusters (grouping according to Witzig et al. [51]). Others (LhISP11, LhISP12, LhISP13, LhISP14, and LhISP15) could not be unambiguously assigned to any of the previously described clusters. Surprisingly, many sequences were grouped in a cluster represented previously only by inferred sequences from environmental clones and one putative dioxygenase from Sphingomonas wittichii RW1. As expected, no sequences clustered with members of the other Rieske nonheme oxygenase subfamilies.


Figure 1
View larger version (62K):
[in this window]
[in a new window]

  		FIG. 1. Residues contributing to the inner surface of the substrate-binding pocket of CumDO from IP01 (CumA1, gray box) and their counterparts in the deduced peptide sequences found in this study. (A) Main cluster sequences. (B) S. wittichii (S.w.) RW1 cluster sequences. The numbers above the panels correspond to the positions in the CumA1 IP01 sequence. Conserved residues are shown in gray, and highly conserved residues are shown in black. Representative sequences from each group are also shown (bold).

 
Members of the main cluster (including sequences grouped with the gram-negative biphenyl dioxygenases plus LhISP11, LhISP12, LhISP13, LhISP14, and LhISP15) had an average of 5 residues differing with respect to the 14 residues responsible for the substrate-binding pocket in IP01 (Fig. 1). Notably, all of them had fixed Gly321 and Phe384 (instead of Ala and Tyr at such positions in IP01). Surprisingly, some sequences lacked putatively important residues; LhISP4 had apparently no active site, LhISP6 lacked the residue His234, and LhISP9 lacked the residue Asp231. LhISP17 and LhISP16 (grouped in the same subcluster) had the residue His240 at position 238 and Asn227 instead of Gln. Members of the S. wittichii RW1 cluster have all Asp-Asn instead of Glu226-Gln227, and LhISP27 and LhISP28 lack the Asp231 at the active site. These last groups have an average of 11 differences in the pocket residues in comparison with those of IP01.

Comparison of soil- versus rhizosphere-derived ISP{alpha} clone libraries.
Comparison of the ISP{alpha} clone libraries derived from the soil and rhizosphere revealed differences between their catabolic gene structures; both clone libraries were statistically different according to the FST and {int}-LIBSHUFF tests (Table 4), and the Shannon, nucleotide, and {theta}({pi}) diversity indexes were higher for the soil library (Table 2). Moreover, there was a strong decrease in the abundance of clones from the S. wittichii RW1 cluster in the rhizosphere library (from 45.76 to 21.2%), especially in the group formed by sequences LhISP19, LhISP20, and LhISP21, which decreased from 27% to 4.2%. On the other hand, there was a rise in the abundance of sequences from the main cluster in the rhizospheric library (54 to 78.6%), driven by the increase in the abundance of sequences LhISP1 (0 to 8.5%), LhISP11 (0 to 19.1%), and, particularly, LhISP8 (2 to 40.4%). Nevertheless, there was a notable decrease in the abundance of other groups, such as those formed by LhISP14-LhISP13 (14.6 to 0%) and LhISP9-LhISP10 (16.6 to 2.1%).


   DISCUSSION
Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
REFERENCES
 
The aim of this study was to evaluate the structural and functional changes occurring in the bacterial communities of a PCB-contaminated soil after the introduction of a suitable host plant for rhizoremediation. This was achieved by screening soil- and rhizosphere-derived 16S rRNA and ISP{alpha} gene libraries via TGGE (since amplicons with identical length migrate dependent upon their respective primary sequences and base compositions [18]). The clones with identical mobility were defined as OTUs. Felske et al. (9) and Bano and Hollibaugh (2) had previously shown that clone screening by denaturing gradient gel electrophoresis or TGGE is a convenient and efficient way to detect and retrieve the predominant sequences in a sample. Nevertheless, to confirm that the grouping of different OTUs according to TGGE mobilities was reproducible and reliable for discriminating between clones, replicates of five different OTUs (one from each pair of gene libraries) were randomly chosen and sequenced, and no differences in sequence composition were obtained, suggesting that such grouping was efficient for this study.

The selection of the bacterial groups to be studied, which was done depending on the group-specific primers available in the literature and the published information on the groups most commonly present in soils, has allowed the study of approximately 67% to 70% of the bacterial diversity from soil and rhizosphere (Table 1). Furthermore, the long sequences obtained with these primers have given us a great deal of phylogenetic information, allowing a fine clustering of the sequences. Except for the Actinobacteria-specific primers, which also amplified sequences clearly belonging to the unrelated Verrucomicrobia, all the primer pairs seem to be truly group specific.

We have used different approaches to analyze the information obtained, each having particular advantages and limitations (29). In the present study, statistically significant differences were found between soil and rhizospheric communities, using tests based on library coverage ({int}-LIBSHUFF) and genetic diversity (FST), for all groups studied except for Actinobacteria. In addition, the catabolic profiles of both communities were found to be significantly different.

The relative abundances of the phylogenetic groups studied fall within the ranges previously reported for soil (for a review, see Janssen [22]) and rhizosphere environments (30, 24). Interestingly, the higher-than-expected abundance of Betaproteobacteria found in this PCB-polluted soil was also observed at another PCB-polluted site in Wittenberg, Germany (31, 32), while in most soils studied, this group has been found in relatively small numbers (22). Moreover, an unusual feature of the Wittenberg soil was the abundance of sequence types within the Betaproteobacteria that were related to Variovorax and Burkholderia species, which are highly represented in our Betaproteobacteria libraries as well. Nogales et al. (32) speculated that such high abundance could be due to their ability to thrive in acidic soils, but in our case the soil pH (7.2) cannot be taken as a cause. The enrichment of such species in PCB-contaminated soils of different origins is noteworthy, as it may represent selection due to the presence of PCBs. Nogales et al. (31) postulated that if the PCB concentration was high enough to account for an important part of the soil's carbon and energy reservoir (as is in our case), some of the sequence types detected should indicate the organisms likely to be involved, directly or indirectly, in the utilization of PCBs as a carbon and/or energy source. Following this rationale, it could be that some Burkholderiales species have an important functional role in the bacterial degradation of PCBs in soil.

In a recent study, based on the analysis of the culturable PCB-degrading bacteria associated with the rhizosphere of different tree species established at a PCB-contaminated site (27), Salix sp. trees were found to positively affect the size of the PCB-degrading populations compared to the results for three other tree species and non-root-containing soil. Thus, if we assume that this could also be true for nonculturable PCB-degrading bacteria, which should be the majority of PCB degraders, some of the sequence types enriched in the rhizosphere (mainly Proteobacteria species) could represent PCB degraders.

All putative ISP{alpha} sequences obtained in this study seem to be members of the toluene/biphenyl subfamily of the Rieske nonheme dioxygenases, as evidenced by BLAST searches and phylogenetic analysis. The sequences obtained show great diversity with respect to previously described Rieske nonheme oxygenases, either from well-studied cultured species or from environmental clones, with identities ranging from 37% to 100%. Many of them contained all the amino acid residues thought to be involved in the active-site functioning and those thought to play a role in electron transfer, and thus, they probably code for functional enzymes. Those sequences lacking some of the important active-site residues (LhISP4, LhISP6, LhISP9, LhISP27, and LhISP28) may represent inactive (vestigial) enzymes.

Since the samples came from a highly PCB-polluted environment (without significant amounts of other aromatic hydrocarbons), it is not surprising to find that many sequences clustered with the gram-negative biphenyl dioxygenase branch instead of the isopropylbenzene and toluene/benzene groups (since a broad correlation between this grouping and substrate range as been previously described [51]). However, the lack of sequences closely related to described gram-positive biphenyl dioxygenase sequences is noteworthy, as more than 8% of the eubacterial soil clones analyzed were Actinobacteria sequences.

We were surprised to find a previously undescribed cluster composed of many sequences from this study but only a few previously described environmental sequences and a putative ISP{alpha} from Sphingomonas wittichii RW1. Such a grouping is clearly separated from the other subfamilies, and yet it seems to contain far greater diversity than the main cluster of the toluene/biphenyl subfamily as judged from its deeper branching (see Fig. S3 in the supplemental material). Moreover, it provides sequence data exploitable for the search of full-length novel dioxygenase genes from a future metagenomic library.

The samples analyzed contained a high degree of catabolic genotype diversity, as evidenced by the low degree of homology of many of the sequences obtained (compared to that of previously described ISP{alpha} sequences), the variability found in the residues responsible for active-site functioning, and the great variation observed in the amino acid residues homologous to those defining the substrate-binding pocket in IP01 and other Rieske nonheme dioxygenases (Fig. 1). This suggests that previously described dioxygenases may not play a significant role in the PCB catabolism in the samples investigated.

Comparison of the ISP{alpha} clone libraries derived from soil and rhizosphere revealed differences between their catabolic gene structures. They were found to be statistically different, and many sequences were clearly differentially distributed in the two environments. Two ISP{alpha} sequences dominated the rhizosphere, LhISP8 and LhISP11, accounting for almost 60% of the gene library. These sequences were most similar to a putative ISP{alpha} from Pseudomonas sp. strain IC (identity 94.4%) and an ISP{alpha} from Pseudomonas putida (identity 75.5%), respectively. It would not be surprising to find out that these genes belong to Pseudomonas species, as they are known to be very good rhizosphere colonizers (28) and their broad catabolic spectrum is known to include ISP{alpha} genes (34). Thus, it would be interesting to identify the bacteria possessing these genes, as well as the substrate spectrum of the putative enzymes, in order to study their applicability for Salix sp.-mediated rhizoremediation.

The different analyses carried out in the present study support the idea that the rhizosphere bacterial community which evolved from the native PCB-polluted soil community during the development of the Salix sp. root system was distinct from its parental community at functional and structural levels. The rhizosphere was enriched in Proteobacteria species and ISP{alpha} genes closely related to the gram-negative cluster. It selected and sustained certain ISP{alpha}-possessing bacteria (putative PCB degraders), plausibly for their competence in root colonization, which could be used for future rhizoremediation strategies. Furthermore, the discovery of a new cluster containing a high degree of sequence variation opens the road for the detection of novel and useful enzymes selected under existing environmental conditions.


   ACKNOWLEDGMENTS
 
We thank Thomas Macek and Martina Mackova for supplying the Lhenice soil.

D. Aguirre de Cárcer was the recipient of a Comunidad de Madrid FPI scholarship. This work was funded by Comunidad de Madrid grant GR/AMB/0084/2004, by the research program MICROAMBIENTE-CM from the Comunidad de Madrid, by the Spanish Ministerio de Educación y Ciencia project BIO2006-08596, and by EU project QLK3-CT-2001-00101.


   FOOTNOTES
 
* Corresponding author. Mailing address: Departamento de Biología, c/Darwin 2, Universidad Autónoma de Madrid, 28049 Madrid, Spain. Phone: 34 914978188. Fax: 34 914978344. E-mail: rafael.rivilla{at}uam.es Back

{triangledown} Published ahead of print on 10 August 2007. Back

{dagger} Supplemental material for this article may be found at http://aem.asm.org/. Back


   REFERENCES
Top
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 

  1. Aguirre de Cárcer, D., M. Martín, M. Mackova, T. Macek, U. Karlson, and R. Rivilla. 2007. The introduction of genetically modified microorganisms designed for rhizoremediation induces changes on native bacteria in the rhizosphere but not in the surrounding soil. ISME J. 1:215-223.[CrossRef]
  2. Bano, N., and J. T. Hollibaugh. 2000. Diversity and distribution of DNA sequences with affinity to ammonia-oxidizing bacteria of the ß subdivision of the class Proteobacteria in the Arctic Ocean. Appl. Environ. Microbiol. 66:1960-1969.[Abstract/Free Full Text]
  3. Barns, S. M., S. L. Takala, and C. R. Kuske. 1999. Wide distribution and diversity of members of the bacterial kingdom Acidobacterium in the environment. Appl. Environ. Microbiol. 65:1731-1737.[Abstract/Free Full Text]
  4. Butler, C. S., and J. R. Mason. 1997. Structure-function analysis of the bacterial aromatic ring-hydroxylating dioxygenases. Adv. Microb. Physiol. 38:47-84.[Medline]
  5. Choi, K., and F. C. Dobbs. 1999. Comparison of two kinds of Biolog microplates (GN and ECO) in their ability to distinguish among aquatic microbial communities. J. Microbiol. Methods 36:203-213.[CrossRef][Medline]
  6. Cole, J. R., B. Chai, R. J. Farris, Q. Wang, S. A. Kulam, D. M. McGarrell, G. M. Garrity, and J. M. Tiedje. 2005. The Ribosomal Database Project (RDP-II): sequences and tools for high-throughput rRNA analysis. Nucleic Acids Res. 33:294-296.[CrossRef]
  7. Curl, E. A., and B. Truelove. 1986. The rhizosphere, vol. 15. Springer-Verlag, Giessen, Germany.
  8. Dong, X. S., S. Fushinobu, E. Fukuda, T. Terada, S. Nakamura, K. Shimizu, H. Nojiri, T. Omori, H. Shoun, and T. Wakagi. 2005. Crystal structure of the terminal oxygenase component of cumene dioxygenase from Pseudomonas fluorescens IP01. J. Bacteriol. 187:2483-2490.[Abstract/Free Full Text]
  9. Felske, A., A. Wolterink, R. Van Lis, and A. D. L. Akkermans. 1998. Phylogeny of the main bacterial 16S rRNA sequences in Drentse A grassland soils (The Netherlands). Appl. Environ. Microbiol. 64:871-879.[Abstract/Free Full Text]
  10. Friedrich, V., A. Wintzingerode, B. C. Ulf, A. Goibel, and E. Stackebrandt. 1997. Determination of microbial diversity in environmental samples: pitfalls of PCR-based rRNA analysis. FEMS Microbiol. Rev. 21:213-229.[CrossRef][Medline]
  11. Furusawa, Y., V. Nagarajan, M. Tanokura, E. Masai, M. Fukuda, and T. Senda. 2004. Crystal structure of the terminal oxygenase component of biphenyl dioxygenase derived from Rhodococcus sp. strain RHA1. J. Mol. Biol. 342:1041-1052.[CrossRef][Medline]
  12. Garland, J. L., and A. L. Mills. 1991. Classification and characterization of heterotrophic microbial communities on the basis of patterns of community-level sole-carbon-source utilization. Appl. Environ. Microbiol. 57:2351-2359.[Abstract/Free Full Text]
  13. Gibson, D. T., and R. E. Parales. 2000. Aromatic hydrocarbon dioxygenases in environmental biotechnology. Curr. Opin. Biotechnol. 11:236-243.[CrossRef][Medline]
  14. Glimm, E., H. Heuer, B. Engelen, K. Smalla, and H. Backhaus. 1997. Statistical comparisons of community catabolic profiles. J. Microbiol. Methods 30:71-80.[CrossRef]
  15. Gomes, N. C. M., H. Heuer, J. Schonfeld, R. Costa, L. Mendonca-Hagler, and K. Smalla. 2001. Bacterial diversity of the rhizosphere of maize (Zea mays) grown in tropical soil studied by temperature gradient gel electrophoresis. Plant Soil 232:167-180.[CrossRef]
  16. Good, I. J. 1953. The population frequencies of species and the estimation of population parameters. Biometrika 40:237-264.[Abstract/Free Full Text]
  17. Hall, T. A. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. 41:95-98.
  18. Heuer, H., B. Engelen, and K. Smalla. 1997. Application of denaturing gradient gel electrophoresis (DGGE) and temperature gradient gel electrophoresis (TGGE) for studying soil microbial communities, p. 353-373. In J. D. van Elsas, E. M. Hwellington, and J. T. Trevors (ed.), Modern soil microbiology. Marcel Dekker, New York, NY.
  19. Heuer, H., M. Krsek, P. Baker, K. Smalla, and E. M. H. Wellington. 1997. Analysis of Actinomycete communities by specific amplification of genes encoding 16S rRNA and gel-electrophoretic separation in denaturing gradients. Appl. Environ. Microbiol. 63:3233-3241.[Abstract]
  20. Holoubek, I. 2001. Polychlorinated biphenyl (PCB) contaminated sites worldwide, p. 17-26. In L. W. Robertson and L. G. Hansen (ed.), Recent advances in environmental toxicology and health effects. The University Press of Kentucky, Lexington, KY.
  21. Huber, T., G. Faulkner, and P. Hugenholtz. 2004. Bellerophon; a program to detect chimeric sequences in multiple sequence alignments. Bioinformatics 20:2317-2319.[Abstract/Free Full Text]
  22. Janssen, P. H. 2006. Identifying the dominant soil bacterial taxa in libraries of 16S rRNA and 16S rRNA genes. Appl. Environ. Microbiol. 72:1719-1728.[Free Full Text]
  23. Kahl, S., and B. Hofer. 2003. A genetic system for the rapid isolation of aromatic-ring-hydroxylating dioxygenase activities. Microbiology 149:1475-1481.[Abstract/Free Full Text]
  24. Kaiser, O., A. Puhler, and W. Selbitschka. 2001. Phylogenetic analysis of microbial diversity in the rhizoplane of oilseed rape (Brassica napus cv. Westar) employing cultivation-dependent and cultivation-independent approaches. Microb. Ecol. 42:136-149.[Medline]
  25. Kent, A. D., and E. W. Triplett. 2002. Microbial communities and their interactions in soil and rhizosphere ecosystems. Annu. Rev. Microbiol. 56:211-236.[CrossRef][Medline]
  26. Kuiper, I., E. L. Lagendijk, G. V. Bloemberg, and B. J. Lugtenberg. 2004. Rhizoremediation: a beneficial plant-microbe interaction. Mol. Plant-Microbe Interact. 17:6-15.[Medline]
  27. Leigh, M. B., P. Prouzová, M. Macková, T. Macek, D. P. Nagle, and J. S. Fletcher. 2006. Polychlorinated biphenyl (PCB)-degrading bacteria associated with trees in a PCB-contaminated site. Appl. Environ. Microbiol. 72:2331-2342.[Abstract/Free Full Text]
  28. Lugtenberg, B., L. Dekkers, and V. Bloemberg. 2001. Molecular determinants of rhizosphere colonization by Pseudomonas. Annu. Rev. Phytopathol. 39:461-490.[CrossRef][Medline]
  29. Martin, P. 2002. Phylogenetic approaches for describing and comparing the diversity of microbial communities. Appl. Environ. Microbiol. 68:3673-3682.[Free Full Text]
  30. McCaig, A. E., L. A. Glover, and J. I. Prosser. 1999. Molecular analysis of bacterial community structure and diversity in unimproved and improved upland grass pastures. Appl. Environ. Microbiol. 65:1721-1730.[Abstract/Free Full Text]
  31. Nogales, B., E. R. B. Moore, W. R. Abraham, and K. N. Timmis. 1999. Identification of the metabolically-active members of a bacterial community in a polychlorinated biphenyl-polluted moorland soil. Environ. Microbiol. 1:199-212.[CrossRef][Medline]
  32. Nogales, B., E. R. B. Moore, E. Llobet-Brossa, R. Rossello-Mora, R. Amann, and K. N. Timmis. 2001. Combined use of 16S ribosomal DNA and 16S rRNA to study the bacterial community of polychlorinated biphenyl-polluted soil. Appl. Environ. Microbiol. 67:1874-1884.[Abstract/Free Full Text]
  33. Page, R. D. M. 1996. TREEVIEW: an application to display phylogenetic trees on personal computers. Comput. Appl. Biosci. 12:357-358.[Free Full Text]
  34. Parales, R. E., and S. M. Resnick. 2006. Aromatic ring hydroxilating dioxygenases, p. 287-341. In J. L. Ramos and R. C. Levesque (ed.), Pseudomonas, vol. 4. Springer, New York, NY.
  35. Porteus, L. A., R. J. Seidler, and L. S. Watrud. 1997. An improved method for purifying DNA from soil for polymerase chain reaction amplification and molecular ecology applications. Mol. Ecol. 6:787-791.[CrossRef]
  36. Sait, M., P. Hugenholtz, and P. H. Janssen. 2002. Cultivation of globally distributed soil bacteria from phylogenetic lineages previously only detected in cultivation-independent surveys. Environ. Microbiol. 4:654-666.[CrossRef][Medline]
  37. Schloss, P. D., B. R. Larget, and J. Handelsman. 2004. Integration of microbial ecology and statistics: a test to compare gene libraries. Appl. Environ. Microbiol. 70:5485-5492.[Abstract/Free Full Text]
  38. Schloss, P. D., and J. Handelsman. 2005. Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness. Appl. Environ. Microbiol. 71:1501-1506.[Abstract/Free Full Text]
  39. Schneider, S., J. Kueffer, D. Roessli, and L. Excoffier. 1997. Arlequin ver. 1.1: a software for population genetic data analysis. Genetics and Biometry Laboratory, University of Geneva, Geneva, Switzerland.
  40. Shannon, C. E., and W. Weaver. 1963. The mathematical theory of communication. University of Illinois Press, Urbana, IL.
  41. Singer, A. C., D. E. Crowley, and I. P. Thompson. 2003. Secondary plant metabolites in phytoremediation and biotransformation. Trends Biotechnol. 21:123-130.[CrossRef][Medline]
  42. Strimmer, K., and A. V. Haeseler. 1996. Quartet puzzling: a quartet maximum likelihood method for reconstructing tree topologies. Mol. Biol. Evol. 13:964-969.
  43. Taylor, P. M., J. M. Medd, L. Schoenborn, B. Hodgson, and P. H. Janssen. 2002. Detection of known and novel genes encoding aromatic ring-hydroxylating dioxygenases in soils and in aromatic hydrocarbon-degrading bacteria. FEMS Microbiol. Lett. 216:61-66.[CrossRef][Medline]
  44. Taylor, P. M., and P. H. Janssen. 2005. Variations in the abundance and identity of class II aromatic ring-hydroxylating dioxygenase genes in groundwater at an aromatic hydrocarbon-contaminated site. Environ. Microbiol. 7:140-146.[CrossRef][Medline]
  45. Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. Clustal W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:4673-4680.[Abstract/Free Full Text]
  46. Van Veen, J. A., L. S. van Overbeek, and J. D. van Elsas. 1997. Fate and activity of microorganisms introduced into soil. Microbiol. Mol. Biol. Rev. 61:121-135.[Abstract]
  47. Villacieros, M., C. Whelan, M. Macková, J. Molgaard, M. Sánchez-Contreras, J. Lloret, D. Aguirre de Cárcer, R. I. Oruezábal, L. Bolaños, T. Macek, U. Karlson, D. N. Dowling, M. Martin, and R. Rivilla. 2005. Polychlorinated biphenyl rhizoremediation by Pseudomonas fluorescens F113 derivatives, using a Sinorhizobium meliloti nod system to drive bph gene expression. Appl. Environ. Microbiol. 71:2687-2694.[Abstract/Free Full Text]
  48. Villadas, P. J., F. Martinez-Abarca, and N. Toro. 2002. Polymerase chain reaction-temperature gradient gel electrophoresis requires the use of high-performance liquid chromatography-purified oligonucleotides. Anal. Biochem. 300:101-103.[CrossRef][Medline]
  49. Reference deleted.
  50. Weisburg, W. G., S. M. Barns, D. A. Pelletier, and D. J. Lane. 1991. 16S ribosomal DNA amplification for phylogenetic study. J. Bacteriol. 173:697-703.[Abstract/Free Full Text]
  51. Witzig, R., J. Howard, H. Hans-Jurgen, and D. H. Pieper. 2006. Assessment of toluene/biphenyl dioxygenase gene diversity in benzene-polluted soils: links between benzene biodegradation and genes similar to those encoding isopropylbenzene dioxygenases. Appl. Environ. Microbiol. 72:3504-3514.[Abstract/Free Full Text]
  52. Yeates, C., A. J. Holmes, and M. R. Gillings. 2000. Novel forms of ring hydroxylating dioxygenases are widespread in pristine and contaminated soils. Environ. Microbiol. 2:644-653.[CrossRef][Medline]
  53. Yee, D. C., J. A. Maynard, and T. K. Wood. 1998. Rhizoremediation of trichloroethylene by a recombinant, root-colonizing Pseudomonas fluorescens strain expressing toluene ortho-monooxygenase constitutively. Appl. Environ. Microbiol. 64:112-118.[Abstract/Free Full Text]

Applied and Environmental Microbiology, October 2007, p. 6224-6232, Vol. 73, No. 19
0099-2240/07/$08.00+0     doi:10.1128/AEM.01254-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental material
Right arrow Other Versions of this Article:
AEM.01254-07v1
73/19/6224    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by de Cárcer, D. A.
Right arrow Articles by Rivilla, R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by de Cárcer, D. A.
Right arrow Articles by Rivilla, R.
Agricola
Right arrow Articles by de Cárcer, D. A.
Right arrow Articles by Rivilla, R.

Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
J. Bacteriol. Microbiol. Mol. Biol. Rev. Eukaryot. Cell All ASM Journals