Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2005, p. 2145-2153, Vol. 71, No. 4
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.4.2145-2153.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
| SHORT REPORT |
CNRS-UMR 8030, Genoscope and Université d'Evry Val d'Essonne, Evry,1 Génopole de lInstitut Pasteur, Paris, France2
Received 21 June 2004/ Accepted 26 October 2004
|
|
|---|
|
|
|---|
In the present paper, using both qualitative and quantitative molecular approaches, we describe the occurrence of a novel high-order bacterial phylogenetic lineage never described before within an anaerobic sludge digester. Moreover we confirm the occurrence and extend the diversity of candidate divisions such as OD1, WS3, OP3, NBL-UPA2, and BRC1 previously described on the basis of a few phylotypes only (5, 11).
Anaerobic sludge samples were obtained from an anaerobic mesophilic digester at the Evry wastewater treatment plant (WWTP; 250,000 population equivalents), located about 35 km south of Paris, France. The digester temperature was 33°C, the pH was 7.2, and the digester had 37.5 days of retention time. Samples were centrifuged for 15 min at 20,000 x g. Sludge pellets were washed three times with phosphate-buffered saline and stored at 20°C. Genomic DNA extraction was performed as described by Chouari et al. (3). 16S rRNA genes were amplified from genomic DNA with a Planctomycetales-specific forward primer, Pla46F, and a universal reverse primer, 1390R (16). Molecular cloning and sequencing were done as described by Chouari et al. (3). The 16S rRNA gene sequences were treated as described by Ewing et al. (8) and compared with the complete EMBL nucleotide sequence databases. Sequences from EMBL with the best BLAST score were imported into the ARB data set when necessary (http://www.arb-home.de). Chimeric sequences were searched by the procedure described by Juretschko et al. (14) prior to phylogenetic analysis. All sequences with more than 1,200 nucleotides were imported into the ARB database and automatically aligned with the existing 16S rRNA gene sequences. The resulting alignments were manually checked and corrected when necessary, and 1,091 unambiguously aligned nucleotide positions were used for phylogenetic analyses. Phylogenetic placement was done in comparison with reference sequences representing the main lines of descent in the domain Bacteria, using the ARB program and database package. Tree topology was evaluated by applying neighbor-joining (NJ) analyses with Jukes and Cantor corrections and a 50% invariance criterion for inclusion of individual nucleotide sequence positions in the treeing analyses. We generated several trees which differed in the reference sequences, the set of alignment positions, and the outgroup sequence used. Maximum-parsimony (MP) and maximum-likelihood (ML) methods were also used. The statistical significance levels of interior nodes were determined by performing bootstrap analyses based on 100 resamplings by the NJ and MP methods. Overall, 16S rRNA gene sequence similarities were determined by using the distance matrix tool of the ARB program package (Fig. 1).
![]() View larger version (37K): [in a new window] |
FIG.1. Evolutionary distance dendrogram showing the affiliation of the environmental 16S rRNA gene sequences recovered from the anaerobic digester clone library to representative members of the divisions in the bacterial domain. The dendrogram was reconstructed by the NJ method in combination with a 50% consensus filter for the domain Bacteria with the ARB software package. The numbers at the nodes indicate the percentages of recovery of relevant branch points in 100 bootstrap resamplings (values above nodes by NJ, values below nodes by parsimony). The root was determined by using the archaeal 16S rRNA sequence of Methanobacterium formicicum (M36508) as the outgroup reference. EMBL, GenBank, and DDBJ accession numbers are given in parentheses. SR1, OP11, OD1, BRC1, OP3, and OP8 represent proposed candidate divisions. The scale bar represents the 10% estimated difference in nucleotide sequence positions.
|
![]() View larger version (46K): [in a new window] |
FIG. 2. Evolutionary distance dendrogram constructed by the NJ method showing the affiliation of the environmental WWE1 16S rRNA gene sequences recovered from the anaerobic digester clone library to representative members of the closest divisions in the bacterial domain. A large data set of spirochete 16S rRNA gene sequences was used. Aligned sequences were analyzed by three methods (NJ, MP, and ML) provided by PAUP 4.0b10 as described in the text. The numbers at the nodes indicate the percentages of recovery of relevant branch points in 100 bootstrap resamplings (values above nodes by NJ, values below nodes by parsimony). The E. coli 16S rRNA gene sequence was used as the outgroup to define the root of the tree. The scale bar represent the 5% estimated difference in nucleotide sequence positions.
|
For dot blot hybridization (DBH) and quantification, total RNA was extracted from approximately 200 mg of six different digester samples (from April to December 2001) by bead beating, phenol extraction, and ethanol-sodium acetate precipitation as described previously (22, 23). The quality of the extracted RNA was checked by formaldehyde gel electrophoresis. Total rRNA was quantitatively normalized with a universal probe as described by Sghir et al. (22). Total rRNA of Escherichia coli (Boehringer, Mannheim, Germany) and a transcribed WWE1 rRNA gene were used as an rRNA standard. A total of 200 ng of RNA was blotted onto nylon membranes (Nytran Super Charge; Schleicher and Schuell) in triplicate and hybridized with 32P-labeled oligonucleotide probes (purchased from Eurogentec, Köln, Germany). A mixture of probes (EUB338, II, and III) was used to quantify specific bacterial rRNA as a percentage of the signal intensity obtained with a universal probe (Beta Imager; Packard Instruments). For the EUB probe mixture, we used 1% SDS and 1% SSC washing buffer under a Tw of 54°C. We determined the experimental Tw for the WWE1-specific probe, and we validated the probe using the same washing buffer composition. WWE1-specific rRNA was quantified and expressed as a percentage of total rRNA measured by EUB I, II, and III as described previously (22).
Fluorescent in situ hybridization (FISH) analyses of sludge samples were done as described by Chouari et al. (3). A stringent wash step was performed for 10 min at 48°C. The newly developed WWE1-specific probe was 5'-end labeled with Cy3 (MWG Biotech, Ebersberg, Germany) and used simultaneously with the fluorescein isothiocyanate-labeled EUB338 (a mixture of probes EUB I, II, and III). We have tested the specificity of the novel probe targeting the novel WWE1 lineage, using pure cultures of Bacillus firmus and Acinetobacter calcoaceticus as negative controls; no signal was obtained with these microorganisms, whereas they showed a positive signal when the mixture of EUB338 probes was used. The systematic use of unlabeled controls in our experiments confirmed that the bright cocci (Fig. 3A) corresponded with autofluorescence. Slides were visualized with a Zeiss axioplan epifluorescence microscope.
![]() View larger version (25K): [in a new window] |
FIG. 3. In situ identification of the WWE1 candidate division in an anaerobic municipal sludge digester, using Cy3-labeled WWE1-specific probe S-*-WWE1-1181-a-A-18. (A) Cells identified, marked with WWE1-specific Cy3-labeled probe. (B) EUB338, II, and III are green. (C) DAPI (4',6'-diamidino-2-phenylindole) staining. The bright cocci in panel A represent autofluorescence. Microscopic visualization was achieved with a Zeiss axioplan epifluorescence microscope.
|
|
View this table: [in a new window] |
TABLE 1. Distribution of clone sequences and OTUs analyzed in the present studya
|
97%). WWE2 seems to form a novel high-order phylogenetic cluster within the Bacteria domain comparable to Planctomycetes or Verrucomicrobia. They can be considered as a Planctomycetes, Verrucomicrobia, or Chlamydia sister group. Our results are corroborated by recent results of Cho et al. (2) showing that Victivallis vadensis could be rather placed within the new phylogenetic phylum of Lentisphaerae (Fig. 1). Thus, WWE2 16S rRNA gene sequences are strongly associated with this new phylum. |
View this table: [in a new window] |
TABLE 2. Affiliation of the bacterial 16S rRNA gene sequences analyzed in this study
|
80% sequence homology. The authors' phylogenetic analysis was based on a small selection of 16S rRNA gene sequences. WWE1 also encompasses two 16S rRNA gene sequences (AJ390461 and AJ390485) retrieved from anoxic bulk soil of a flooded rice microcosm (5). In their phylogenetic analyses, Derakshani and coworkers showed that the two sequences cluster together in a monophyletic cluster they called SP lineage (with reference to spirochetes) showing distinct affiliation with the spirochetes with an overall level of 16S rRNA gene sequence dissimilarity ranging between 23 to 27%. If Leptospira spp., which form a deeply branching sublineage of spirochetes, are included in the phylogenetic analyses, the bootstrap support for assignment of SP lineage to the spirochete drops from 93% to 40 to 60%. In our study, we enriched this group of 16S rRNA gene sequences by 322 sequences. We used a comprehensive number of cultivated spirochete 16S rRNA gene sequences and included sequences from reports by Derakshani et al. (5) and Gu et al. (10): using the three treeing methods, phylogenetic analyses indicate that WWE1 and Spirochaetes are clearly making two monophyletic and distinct clusters (Fig. 2). We show that sequences from the reports by Derakshani or Gu et al., along with our WWE1 rRNA gene sequences, cluster in a monophyletic lineage with a bootstrap value of >93%. All treeing analyses using a number of different bacterial division sequence representatives show that WWE1 is a monophyletic lineage branching distinctly from Spirochaetes and other main lines of bacterial descent. Thus, WWE1 appears as a Spirochaetes division sister group, clearly distinguishable and forming a completely independent lineage (Fig. 2). 16S rRNA gene sequence dissimilarity between WWE1 and adjacent phyla (Lentisphaerae, Verrucomicrobia, and Planctomycetes; BRC1) is more than 20%. Sequence dissimilarity with the Spirochaetes phylum ranged between 20.9 and 28.7%. Hugenholtz et al. (12) suggested that 85% sequence identity be used as a cutoff for distinguishing new phyla; thus, by this criterion together with phylogenetic position, this clade was considered to be a novel phylum-level lineage in the domain Bacteria. Comparison of 16S rRNA sequence signatures which distinguish among Bacteria, Archaea, and Eukarya was performed as described by Woese (25). Sequence signatures of WWE1 were clearly related to Bacteria and might be considered as Spirochaetes and Verrucomicrobia sister groups, respectively (data not shown). The percentage of intradivergence within WWE1 lineage is 29%, which is comparable to that of other well-established divisions (6).
Representative sequences of these lineages have also been retrieved from two other anaerobic digesters, such as the Corbeil (2 km upstream of Evry) and Creil (100 km downstream of Evry) WWTPs. They represented 14.7 and 19.2% of the generated clone libraries (data not shown). This means that microorganisms harboring such 16S rRNA gene sequences are environmentally widely distributed. They might be important components of the trophic chain of degradation of the organic matter under anaerobic conditions.
Quantitative DBH and FISH techniques are providing quantitative information on the abundance of a microorganism or a group of microorganisms, their shape and their spatial localization within the original environmental setting. We developed a WWE1-specific probe, S-*-WWE1-1181-a-A-18 (Table 3). The washing temperatures for the WWE1 probe was determined to be 50°C, using a washing buffer containing 1x SSC and 1% SDS as described by Sghir et al. (21, 22). Application of the probe using DBH on six different total RNA extracts from anaerobic digester samples showed that WWE1 rRNA represents 11.9% ± 3.1% of the total bacterial rRNA, which is comparable to the relative rRNA index of a bacterial phylum like Proteobacteria, represented by 14% of total rRNA (unpublished data). The same probe was Cy3-5'-end labeled and used to track representatives of WWE1 within the digester samples by the FISH technique. Results show that the probe binds to rods and filamentous shaped microorganisms but not spiral shaped organisms, which is a characteristic of the Spirochaetes phylum (Fig. 3). Similar results were obtained by Rossetti et al. (19) using the SPL998 16S rRNA probe, which targets a set of sequences affiliated with the WWE1 lineage (data not shown). This provides further support for considering this major lineage as a candidate division, which may be recovered from methanogenic environments. The novel probe described in this study will enable us to determine the distribution and population dynamics of the WWE1 lineage in environmental samples such as sludge digesters and termite, animal, or human gut.
|
View this table: [in a new window] |
TABLE 3. Novel candidate division-specific probe S-*-WWE1-1181-a-A-18 based on the 5' 3' 16S rRNA sequence alignment
|
Twelve 16S rRNA gene sequences were grouped within a recently described candidate division, OD1 (11). Sequence similarity of OD1 rRNA gene sequences to known environmental sequences ranges between 77.4 and 79.6%. The BRC1 candidate division proposed on the basis of 11 sequences by Derakshani et al. (5) is represented by 12 sequences defining four operational taxonomic units (OTUs), showing between 77.9 and 85.6% similarity with the most closely affiliated uncultured soil bacterium AF390481 (Fig. 1). Thirteen more BRC1 sequences were recovered from aerobic and anoxic basins (data not shown). Other sequences were found to be affiliated with the candidate division NBL-UPA2 (5), represented by three sequences defining one OTU. Five more sequences were retrieved from the aerobic and anoxic basins (data not shown); they cluster together and show 90.6% similarity to the uncultured soil bacterium AJ390463. They extend this recently described candidate division and thus confirm the existence of this novel bacterial lineage (Fig. 1).
Due to the limitation imposed by the so-called universal primers, we explored a novel fruitful strategy based on the use of a specific division-level primer set for assessment of the diversity of quantitatively important but less known bacterial divisions. Many unique lineages were sampled within the digester ecosystem by using a specific primer set targeting the Planctomycetales division. These groups may constitute novel divisions or subphyla, greatly expanding the known higher-order diversity of Bacteria. The discovery of division or subphylum phylogenetic lineages is a step toward exploring the relationship between bacterial diversity and biogeochemical function within the digester ecosystems. However, further elucidation of their role in such complex ecosystem may require new culturing approaches or metagenomic studies that allow the linkage of functional and rRNA genes.
|
|
|---|
We are very grateful to S. Cure for reading the manuscript, the excellent technical assistance of the Genoscope sequencing team, and M. Trouvé and D. Dehon for providing samples from the WWTP.
|
|
|---|
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»