Previous Article | Next Article ![]()
Applied and Environmental Microbiology, July 2007, p. 4191-4197, Vol. 73, No. 13
0099-2240/07/$08.00+0 doi:10.1128/AEM.00387-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
,
Netherlands Institute of Ecology (NIOO-KNAW), Centre for Terrestrial Ecology, Boterhoeksestraat 48, 6666 GA Heteren, The Netherlands
Received 19 February 2007/ Accepted 24 April 2007
|
|
|---|
|
|
|---|
Until now, collimonads have been quantified only for the acidic dune grassland soils from which they were originally isolated (11). In these soils, the numbers of collimonads ranged from 103 to 105 CFU per g (dry weight) soil. Enumeration was based on plate counts of chitin-degrading colonies on agar plates containing colloidal chitin. On such plates, Collimonas strains can be recognized as halo-producing bacteria due to clearing of chitin, with concomitant production of translucent biomass. However, identification of collimonads on the basis of colony morphology cannot be conclusive without a more specific identification method. For this purpose, we developed and describe here a Collimonas-specific restriction fragment length polymorphism (RFLP) assay based on restriction analysis of PCR-amplified 16S rRNA.
Enumeration on chitin agar plates provides an indication of the abundance of collimonads. However, collimonads are relatively slow-growing bacteria, and, therefore, they can remain undetected when fast-growing chitinolytic bacteria are present as well. Additionally, this method does not detect potentially nonculturable collimonads. As an alternative to plate enumeration, we developed and applied a culture-independent real-time quantitative PCR (qPCR) assay (3, 21) for quantification of collimonads in soil. Real-time PCR has been used successfully to detect and quantify bacterial cells in various environmental samples (1, 3, 20, 25, 26), including soils (16, 21, 30, 34, 38, 42). In the current study, the presence of soil collimonads was examined for a wide range of soils (40 sites) using both the plate count-RFLP method and the real-time PCR assay. In order to identify the possible factors that determine collimonad soil population sizes, the presence or absence of culturable collimonads and the real-time PCR-based total numbers of collimonads were compared to fungal density based on ergosterol measurements, as well as other soil properties, vegetation composition, and management practices.
|
|
|---|
Estimates of soil fungal biomass were based on measurements of the soil ergosterol content (mg per kg soil). Ergosterol is the major sterol in the membrane of most fungi and is not common outside the fungal kingdom (46, 47). Soil ergosterol was extracted using an alkaline extraction procedure and was analyzed by high-performance liquid chromatography as described elsewhere (12).
Total DNA from soil samples was extracted using a Power Soil DNA isolation kit (MOBIO Laboratories, Solana Beach, CA) according to the manufacturer's instructions except that two 30-s bead-beating treatments with a Mixer Mill MM301 (Retsch, Haan, Germany) were used instead of vortex mixing. In addition, DNA was eluted in a final volume of 50 µl instead of 100 µl. The DNA extract was diluted 10-fold before it was used as a template for real-time PCR quantification.
Detection and identification of culturable collimonads by RFLP analysis of 16S rRNA.
Soil suspensions were prepared as described elsewhere (10). Fifty-microliter portions of 10-fold dilutions were plated on chitin-yeast extract agar with the composition described previously (11) but with an extra 0.1 g per liter of the fungal inhibitor delvocid (DSM, Delft, The Netherlands) (11). The plates were sealed with Parafilm, incubated at 20°C for 2 weeks, and inspected regularly for halo-forming (chitinolytic) bacterial colonies. Among the chitinolytic bacterial colonies, we screened for Collimonas-like isolates based on the morphology described by de Boer et al. (11). In addition, other chitinolytic colony types were sampled. In total, 205 chitinolytic isolates were identified. All isolates were streaked on 0.1x tryptone soy broth (TSB) (Oxoid) agar containing chitin (11). This TSB-chitin agar was used to see if the chitinolytic ability of the strains was repressed by the presence of TSB, as such catabolic repression of chitinase production has been reported for all Collimonas strains described so far (11). Plates were screened for chitinolytic activity after 2 weeks of incubation at 20°C.
Total DNA of the 205 chitinolytic isolates was extracted from 3-ml King's broth cultures (King's B) incubated at 27°C on a rotary shaker (200 rpm) for 48 h, using the procedure used for extraction of DNA from soil samples.
Nearly complete fragments of 16S rRNA genes were amplified from isolated genomic DNA using the universal bacterial primers pA (5'-AGAGTTTGATCCTGGCTCAG-3') and 1492r (5'-GRTACCTTGTTACGACTT-3') (15). Bacterial DNA (1 µl of 10x diluted genomic DNA) was added to a 25-µl (final volume) mixture containing 0.6 µM of each primer, 200 µM of each deoxynucleoside triphosphate, 2.5 µl of 10x buffer (Promega, Leiden, The Netherlands), and 0.056 U Taq polymerase (Promega). The PCR was performed using a touchdown program in which the annealing temperature initially decreased from 65 to 55°C by 2°C per cycle, followed by 12 cycles at 55°C for 1 min. The denaturing step was 30 s at 92°C, and the extension step was 2 min at 68°C.
The restriction enzyme BstBI (New England Biolabs, MA) was used to digest the PCR amplicons. Based on the current information from the Ribosomal Database Project (Ribosomal Database Project II; http://rdp.cme.msu.edu/), we determined that collimonads are the only members of the family Oxalobacteraceae which share a BstBI restriction site in the 16S rRNA gene at the position corresponding to nucleotides 993 to 998 of the 16S rRNA gene of Collimonas fungivorans Ter331 (accession number AJ310395). Restriction fragments were examined on 1.8% agarose gels in 0.5x TSB buffer (11), using O'GeneRuler 1-kb DNA Ladder Plus (Fermentas, St. Leon-Rot, Germany).
The specificity of the RFLP assay for collimonads was tested using five representative strains belonging to genera closely related to Collimonas, namely, three Herbaspirillum type strains (H. seropedicae DSM 6445T, H. rubrisubalbicans ATCC 19308T, and H. frisingense DSM 13128T) and two Janthinobacterium type strains (J. agaricidamnosum DSM 9268T and J. lividum DSM 1522T) (11). C. fungivorans Ter331 was used as a positive control (11).
PCR products that were identified by RFLP analysis as products from Collimonas, as well as 38 PCR products that were identified as not originating from Collimonas, were sequenced (Greenomics, Wageningen, The Netherlands) using the universal primer U1115R (5'-TCCCGCAACGAGCGCAACC-3') (4) as the sequencing primer. DNA sequences up to 650 bp long were compared with sequences available in the GenBank (http://www.ncbi.nlm.nih.gov/) and Ribosomal Database Project II (http://rdp.cme.msu.edu/) databases. The sequences were aligned and compared using Clustal-W in the Lasergene DNA and protein analysis software (DNASTAR, Madison, WI).
Development, validation, and application of a Collimonas-specific real-time qPCR assay.
In a dual-labeled probe assay, the primers used for real-time PCR quantification of collimonads were Eddy3for and Eddy3rev. The forward primer Eddy3for (5'-GTACAGAATCCCGAAGAGATTTGG-3') was based on a previously reported fluorescent in situ hybridization probe specific for Collimonas species (11). In combination with the nonspecific reverse primer Eddy3rev (5'-ACTTAACCCAACATCTCACGACA-3'), it yields a 100-bp amplicon.
To obtain assay specificity, we designed the Collimonas-specific probe Sophie (5'-FAM-CGAAAGA+AA+GC+TG+TA+ACACAGG-BHQ-3') (Fig. 1), which contains 6-carboxyfluorescein (FAM) as the fluorophore, butylhydroquinone (BHQ) as a quencher, and several locked nucleic acids (a plus sign indicates a locked nucleic acid base). A locked nucleic acid is a modified nucleic acid with increased binding affinity for cDNA sequences (22). The primers and probe were synthesized by Biolegio BV (Nijmegen, The Netherlands) and Sigma-Proligo (Paris, France), respectively.
![]() View larger version (8K): [in a new window] |
FIG. 1. Alignment of partial 16S rRNA sequences from previously described collimonads and closely related species. The numbers correspond to the numbers in the 16S rRNA sequence of C. fungivorans Ter331 (accession no. AJ310395). Dashes indicate nucleotides identical to those in the Ter331 sequence. The locations of primers Eddy3for and Eddy3rev, as well as the dual-labeled probe Sophie, are indicated. The boxed nucleotides indicate the BstBI restriction site unique to collimonads. The accession numbers for the 16S rRNA sequences of strains Ter166, Ter300, Ter228, Ter14, Ter165, Ter299, Ter90, Ter94, Ter227, Ter291, DSM 9628T, DSM 1522T, DSM 6445T, ATCC 19308T, and DSM 13128T are AY281140, AY281145, AY281148, AY281135, AY281139, AY281144, AY281136, AY281138, AJ496445, AY281143, Y08845, Y08846, Y10146, AB021424, and AJ238358, respectively.
|
The ABsolute qPCR mixture (ABgene) was used at a final concentration of 1x for the real-time reaction. Primers Eddy3for and Eddy3rev were added to final concentrations of 100 and 70 nM, respectively. The final concentration of the Sophie probe was 50 nM. The volume of the reaction mixture was adjusted to 20 µl using nuclease-free water. The reaction mixture contained 400 ng of bovine serum albumin per µl to minimize soil PCR inhibition (16, 24). After addition of 5 µl DNA extracts from soil or control DNA (see below), amplifications were performed using the following conditions: 15 min at 95°C, followed by 45 cycles of 15 s for 95°C and 45 s for 66°C. Each DNA extract was tested in duplicate. In all cases, negative controls containing nuclease-free water instead of DNA extract were included.
For quantification of collimonads, a standard curve was generated using serial dilutions of genomic DNA isolated from C. fungivorans Ter331. Standard curves were generated by plotting threshold cycles (Ct) versus genome equivalents of strain Ter331. The undiluted genomic DNA extract contained 40 µg/ml genomic DNA, as measured spectrophotometrically. The size of the genomic DNA of C. fungivorans Ter331 has been estimated to be 5.1 Mbp (unpublished data). Considering that 1 kb of double-stranded DNA is equivalent to 6.5 x 105 Da (1 Da = 1.65 x 1024 g), one genome of C. fungivorans Ter331 weighs 5.5 x 1015 g. Thus, undiluted DNA extract contained 7.3 x 106 genomic DNA equivalents per 1 µl. We considered one genome equivalent to one cell.
Ct values indicate the minimum number of PCR cycles needed to obtain fluorescence signals that are significantly greater than the background for the exponential phase of PCR amplification. The Rotor-Gene 6 software (Corbett Research, Sydney, Australia) was used to establish Ct values for each sample.
Statistics.
The chi-square test was used to compare the patterns of distribution of culturable collimonads in arable, forest, and (semi)natural grassland soils. For this test, sites cropped with maize and fertilized grassland were categorized as arable land, nonfertilized grassland and formerly arable lands that were abandoned for >10 years were categorized as (semi)natural grasslands, and recently abandoned arable sites (abandoned for <10 years) were not included in the analyses of land management since they represent a transient stage from arable lands to natural areas.
Data for real-time PCR-based (log-transformed) numbers of Collimonas cells and soil ergosterol contents were analyzed by one-way analysis of variance with different land management practices (arable, forest, and [semi]natural grassland) and different layers (mineral and organic) of forest soil as treatments. Differences between groups were tested for significance with a modified Tukey's honestly significant difference test at a P value of <0.05. Possible relationships between real-time PCR-based numbers of Collimonas cells and physical and chemical characteristics of soil were examined using correlation analyses. Principal-component analysis was performed to ordinate the soils on the basis of soil characteristics. Scaling was focused on intersample distances. Species scores (soil characteristics) were divided by the standard deviation and used for centering. The analysis was performed using CANOCO 4.5 (41). All other statistical analyses were performed using STATISTICA (Statsoft Inc., Tulsa, OK).
|
|
|---|
![]() View larger version (19K): [in a new window] |
FIG. 2. (a and b) Average ergosterol contents (mg kg1 [dry weight] soil) in the upper layer (0 to 10 cm) of the mineral soils at arable, forest, and (semi)natural grassland sites (a) and in the organic and mineral layers at five forest sites (b). (c and d) Average log-transformed numbers of Collimonas genome equivalents per gram (dry weight) of soil in the upper layer (0 to 10 cm) of the mineral soils at arable, forest, and (semi)natural grassland sites (c) and in organic and mineral layers at five forest sites (d). Different letters indicate significant differences at the 0.05% level (as determined by a modified Tukey's honestly significant difference test).
|
![]() View larger version (63K): [in a new window] |
FIG. 3. RFLP analysis of BstBI-digested 16S rRNA genes amplified almost to completion from J. agaricidamnosum DSM 9628T (lane A), J. lividum DSM 1522T (lane B), H. rubrisubalbicans ATCC 19308T (lane C), H. seropedicae DSM 6445T (lane D), H. frisingense DSM 13128T (lane E), C. fungivorans Ter331 (lane F), Collimonas sp. from soil sample 19 (lanes G and H), Collimonas sp. from soil sample 26 (lane I), Collimonas sp. from soil sample 15 (lane J), Collimonas sp. from soil sample 26 (lane K), isolate 142 from soil sample 16 (= Flavobacterium sp.) (lane L), isolate 134 from soil sample 27 (= Burkholderia sp.) (lane M), isolate 186 from soil sample 36 (= Rhodanobacter sp.) (lane N), and isolate 88 from soil sample 9 (= Pedobacter sp.) (lane O). The lanes indicated by an asterisk contained a 1-kb marker. The strains in lanes A to E were obtained from DSMZ (Braunschweig, Germany). The strains in lanes G to O were chitinolytic isolates from our inventory study. Identification of the strains in lanes G to O (see above) was based on sequence analysis of the 16S rRNA gene fragment amplified by primers pA and 1492r (15).
|
The 26 confirmed collimonads originated from nine different sampling sites, sites 4, 5, 6, 15, 18, 19, 22, 26, and 27. These sites were eight (semi)natural grassland sites and one heathland soil (Table 1). All 26 Collimonas isolates had a colony morphology that resembled the colony morphologies described for collimonads (11). However, two isolates obtained from grassland site 26 produced a purple pigment, which was previously described for species belonging to the related genus Janthinobacterium (11, 28, 33). The purple color was not clear on the original counting plates. However, when colonies were restreaked on both chitin-yeast extract agar and TSB-chitin plates for purity, the purple color became apparent. Production of purple pigments by collimonads has been reported previously (32).
|
View this table: [in a new window] |
TABLE 1. Numbers of arable sites, (semi)natural grasslands, and forests in which culturable collimonads were present or absenta
|
Development and validation of a Collimonas-specific dual-labeled probe assay.
Although the aforementioned combination of plate counts and RFLP-based identification gave us an indication of the distribution of collimonads in different soil types, there were still many limitations, including the presence of nonculturable collimonads and the apparent absence of collimonads due to the presence of fast-growing chitinolytic gliding bacteria. Therefore, we developed a 16S rRNA gene-based real-time PCR assay. A Collimonas-specific primer/probe set was designed by alignment of 16S rRNA gene sequences of 12 representative collimonads (11) and the 22 most closely related bacterial species belonging to the family Oxalobacteraceae (accession numbers AY167838, Y08845, Y08846, AF174648, AY247410, Y10146, AB021424, AJ238358, AB074524, U49757, PSJ001384, AF543312, AB008506, AB024305, AY133107, AF529336, AF154097, AY177773, AF529095, AF358019, AY214204, and AF407411).
The probe, which was designed based on a 22-bp stretch shared only by collimonads, had a single mismatch (G or A at position 10) (Fig. 1). The single mismatch of the probe did not influence the Ct value (data not shown), indicating that the probe was specific for real-time PCR analysis of collimonads.
The specificity of the dual-labeled probe assays was tested empirically using genomic DNA of representative species of genera closely related to Collimonas (i.e., DNA of H. seropedicae DSM 6445T, H. rubrisubalbicans ATCC 19308T, H. frisingense DSM 13128T, J. agaricidamnosum DSM 9268T, and J. lividum DSM 1522T). The dual-labeled probe assay did not give a signal during 40 cycles of the reaction (results not shown).
Direct detection and quantification of Collimonas in soils using the dual-labeled probe PCR assay.
Figure 4 shows a composite standard curve that was derived using Ct values from five independent runs of PCR assays. The standard curve was obtained for the range from 3.7 x 101 to 3.7 x 106 genome equivalents per reaction mixture with an R2 value of 0.96 and a coefficient of variation of 10.3% ± 2.58%. The amplification efficiency, calculated using the methods described by Pfaffl (37), was 1.54 ± 0.07. The lower limit of quantification was 3.7 x 101 genome equivalents per reaction mixture, which corresponds to a mean Ct value of 38, whereas the negative controls did not give a signal during the 40 cycles of the reaction. The lower limit, 3.7 x 101 genome equivalents per reaction mixture, equals 3.7 x 101 genome equivalents per 5 µl of 10x-diluted DNA extract or 3.7 x 103 genome equivalents per 50 µl of undiluted DNA extract. The 50 µl was extracted from 0.25 g (dry weight) soil, which implies that the quantification limit of our real-time PCR assay was 1.5 x 104 Collimonas genome equivalents (cells) g1 (dry weight) soil. The standard curve that was constructed was used to estimate the number of genome equivalents or cells of Collimonas in all soil samples.
![]() View larger version (14K): [in a new window] |
FIG. 4. Standard curve for the Collimonas-specific dual-labeled probe PCR assays. This curve was created using known amounts of DNA ranging from 3.7 x 101 to 3.7 x 106 genome equivalents of C. fungivorans Ter331 genomic DNA per reaction mixture. The values shown are representative results derived from multiple independent assays (n = 5).
|
|
|
|---|
Real-time PCR has the advantage that it can be used to quantify both culturable and nonculturable bacteria in the soil microbial community (21, 30, 42). The discrepancy between real-time PCR detection and the results of the isolation experiments (Table 1) could be explained by the inability of some Collimonas strains to grow on the enumeration plates (chitin-yeast extract agar). Another explanation could be that fast-growing antagonistic or gliding bacteria prevented or masked the (slow) development of Collimonas colonies.
In theory, the quantification limit for plate enumeration is 2.0 x 102 Collimonas cells g1 (dry weight) soil. However, in practice, the detection limit is higher due to suppression of colony development by fast-growing antagonistic and gliding chitinolytic bacteria. The quantification limit for our real-time PCR assay was 1.5 x 104 Collimonas cells g1 (dry weight) soil. Our qPCR results indicate that the size of the population of collimonads was about 105 g1 (dry weight) soil in 60% of soils that we tested. In the other soils containing Collimonas (20%), the numbers were close to the quantification limit. The presence of 105 cells g1 soil corresponds to at most 0.01% of the total soil bacteria (based on 4',6'-diamidino-2-phenylindole [DAPI] counts) (result not shown). In the nine soils where they were detected, culturable collimonads accounted for 0.1% of the total culturable chitinolytic bacteria (data not shown). Hence, it seems that the abundance of collimonads is generally low. This does not indicate a strong competitive ability in soil, even though it has been shown that collimonads are metabolically versatile (11). Collimonas was originally isolated in a search for chitinolytic bacteria. However, it is known that collimonads are weak chitin degraders compared to fungi and actinomycetes (8). Hence, their growth on chitin sources, mainly dead fungal hyphae and exoskeletons of arthropods, in soil is probably insignificant.
As another approach to get an indication of the abundance of collimonads, we used BLAST to screen for the presence of Collimonas sequences in published gene clone libraries of bacteria. Collimonads were in 16S rRNA gene clone libraries of soils in only three studies, supporting our finding that these organisms do not occur in high numbers in soil. In one of these studies, collimonads were found in the lichen-dominated surface of forest soils, and the authors suggested that they may play a role in the decay of the lichen biomass, which consists largely of fungi (32). Another study showed that collimonads were found in a forest soil-based microbial biofilm exposed to 2,4-dichlorophenoxyacetic acid (2). Opelt and Berg found collimonads in the bacterial community associated with bryophytes living in the Baltic Sea coast area (36).
Real-time PCR results indicated that collimonads were widely distributed, but their numbers were significantly higher in (semi)natural ecosystems than in arable soils. In addition, culturable Collimonas cells were detected only in natural grasslands. The relatively low numbers of collimonads in arable soils and production grasslands are in line with mycophagous growth, as these sites are fungus poor due to agricultural management practices (44). The low numbers are also compatible with a parasitic or predatory growth strategy (12). For example, Bdellovibrio species, which are bacterial predators of other bacteria, also occur in low numbers (102 to 104 CFU per g1) in soils (23).
Despite the apparent preference for the more fungus-rich natural ecosystems, the numbers of collimonads were not pairwise correlated with ergosterol-based soil fungal biomass. The lack of a statistically significant correlation between the fungal biomass and the numbers of collimonads could imply that mycophagy is not an essential growth strategy for collimonads under natural conditions. Alternatively, enumeration of collimonads by qPCR may also include as-yet-unknown strains that are not mycophagous. Furthermore, preferential feeding of Collimonas strains on specific taxonomic or functional fungal groups may interfere with a correlation between ergosterol-based soil fungal biomass and numbers of collimonads. Indeed, under lab conditions preferential growth of Collimonas strains was demonstrated. The proliferation of collimonads on the zygomycete Mucor hiemalis was much greater than that on the ascomycetes Chaetomium globosum and Fusarium culmorum (10).
Surprisingly, within the forest soil profile, collimonads were present in relatively low numbers in the fungus-rich forest organic layers. If mycophagy is an important in situ growth strategy for collimonads, this must indicate that the environmental conditions or the type of fungal species in the organic layer of forest soils is not favorable for collimonads. Recently, Lindahl et al. (29) reported spatial separation in the vertical distribution of different functional groups of fungi in the organic layer of boreal forest soils. It was observed that saprotrophic fungi were primarily confined to the top organic layer (litter), whereas ectomycorrhizal fungi dominated the lower layers (29, 39). Hence, it may be that the higher numbers of collimonads in the mineral layer are due to a preference for hyphae of ectomycorrhizal fungi.
A functional group of fungi with a short turnover time for hyphae is the arbuscular mycorrhizal (AM) fungi (40). If these fungi are an important source of food for collimonads, then this should not be apparent from ergosterol measurements as AM fungi do not contain ergosterol (35).
Obviously, other approaches are needed to indicate whether mycophagy is important for collimonads under natural growth conditions. One possible strategy is to introduce fungal hyphae as bait in soil samples that naturally contain collimonads. Application of the real-time PCR assay that we developed here would allow us to follow changes in Collimonas numbers and assess whether the mycophagous property is essential for in situ growth. Future studies also need to be aimed at increasing our knowledge of the types of fungal species that are preferably targeted by collimonads. This should enable reassessment of a possible link between the abundance of collimonads and the presence of fungi by focusing more specifically on the presence of fungal species that are actually susceptible to collimonads. Depending on their specificity for different fungi, collimonads could have different functions in soil environments and could range from being deleterious (e.g., due to disturbance of symbiosis between plants and AM fungi) to being beneficial to plants by targeting plant-pathogenic fungi.
Published ahead of print on 4 May 2007. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
Publication no. 4047 NIOO-KNAW of The Netherlands Institute of Ecology. ![]()
|
|
|---|
milauer. 2002. CANOCO reference manual and CanoDraw for Windows user's guide: software for canonical community ordination (version 4.5). Microcomputer Power, Ithaca, NY.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»