Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2004, p. 483-489, Vol. 70, No. 1
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.1.483-489.2004
Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, Colorado,1 Department of Biology, San Diego State University, San Diego, California,2 Environmental Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee3
Received 23 July 2003/ Accepted 14 October 2003
|
|
|---|
|
|
|---|
Bacteria are renowned for their rapid evolution in response to novel selection pressure, and any environment subject to varying selection, either spatially or temporally, may harbor suites of bacteria that are capable of rapid change. The emergence and spread of antibiotic resistance (1) are perhaps the best known examples. In addition, the bioremediation literature is full of references to bacteria that possess unique genes that metabolize toxic chemicals (e.g., 11, 38). Many more examples of rapid evolution of metabolic characters have been described for a diverse range of bacterial species, and most cases involve the emergence of novel genes and their spread in environments that are subject to marked human impact (10, 26). Although important information about the metabolic versatility of bacteria in human-impacted environments can be gleaned from the literature, whether such versatility is a general property of natural microbial communities is less well known.
Pseudomonas, an enormously diverse genus of the
-Proteobacteria, is an important member of soil microbial communities (27). Members of the genus have been isolated from essentially all environments studied (28), including alpine soil, where it was identified as the most prevalent culturable genus in Kobresia alpine meadows (24). The genus exhibits remarkable metabolic variation (31), and a large number of different plasmids have been described for it, including enormous plasmids containing many genes (e.g., the IncP-9 TOL plasmid pWW0 in Pseudomonas putida is over 110 kb and contains 148 open reading frames [7]). The great metabolic flexibility of Pseudomonas species may allow them to inhabit variable environments. One strategy might be the evolution of strains that are capable of utilizing a large number of different carbon sources for growth. Alternatively, because the alpine environment is highly heterogeneous with pockets of specific carbon compounds, a large number of different strains that have recently gained or lost the ability to grow on particular sources of carbon may exist.
In this study we report the prevalence of fluorescent pseudomonads and describe the isolation and characterization of 17 cold-tolerant strains of Pseudomonas from high-alpine soil in Colorado. Strains were grown on 20 different carbon sources to assess their metabolic characteristics, and their 16S ribosomal DNA (rDNA) sequences were determined for phylogenetic analysis. Evaluation of the metabolic properties relative to the phylogeny revealed that Pseudomonas bacteria from high-alpine soil lack phylogenetic diversity but exhibit great metabolic versatility. The lack of concordance between the metabolic data and the inferred phylogeny may reflect rapid gains and losses of genes.
|
|
|---|
Isolation.
To obtain pure Pseudomonas isolates that were significantly represented in the soils, we used two related isolation methods. The first method, a limiting dilution culture, isolated organisms that might otherwise be outcompeted in the laboratory (4, 20). In this method, the isolates came from the highest dilution showing growth. The second method, enrichment cultures, selected fast-growing organisms by using the 10-3 dilutions. The culture medium for the isolation, enrichment, dilution, and growth studies was the same mineral salt solution with sterile soil extract as that used for the MPN counts, but with a C source (0.2 g liter-1). Table 1 outlines the dates of soil collection and the methods, C sources, and temperatures of isolation.
|
View this table: [in a new window] |
TABLE 1. Details of isolation conditions of study strains
|
Metabolic studies using different carbon sources.
We determined whether each isolate grew on each of 20 different C sources. Table 2 lists the C sources that were tested. The growth experiments were done at 22°C. A standard 96-well plate with mineral salt-glutamate medium served as the inoculation master. Duplicates for each isolate grew in adjacent wells. To maintain the separation of the isolates, we grew our isolates only in every other column of the plates; the uninoculated columns served as control blanks and also as a method to guard against cross-contamination of the wells. Five replicate control plates contained only mineral salt medium and no C compound. Using a replicator (Boekel Scientific, Inc., Feasterville, Pa.), we inoculated cells from the master plate into each control and C source plate. A Spectra Max 340PC plate reader driven by Soft Max Pro 2.6.1 software (Molecular Devices Corporation, Milpitas, Calif.) measured growth at an optical density (OD) of 595 nm. To evenly suspend the cells in the medium, we manually agitated each plate and set the plate reader to agitate the plate just before each growth reading. Growth was measured daily for 7 days. Growth was defined as present if, over the duration of each growth study, any mean OD of the two replicates was recorded as being greater than 1 standard deviation above the mean control OD.
|
View this table: [in a new window] |
TABLE 2. Growth on carbon sources for each isolate
|
|
View this table: [in a new window] |
TABLE 3. Primer dataa
|
Posterior probabilities were determined by Bayesian Markov chain Monte Carlo methods implemented by using MrBayes (13). A general time-reversible model with gamma rate heterogeneity was adopted; 500,000 generations were run, and the trees and model parameters were sampled every 100 generations. The posterior probability distribution stabilized after 36,000 generations, and so this number was adopted as the burn-in value (meaning that all parameter estimates prior to generation 36,100 were omitted). Branch lengths were estimated by using maximum likelihood and the modal parameter values estimated for the substitution model from the Bayesian analysis. The likelihood analysis was performed with PAUP (version 4.0b8a; Sinauer Associates, Inc., Sunderland, Mass.).
The concordance between the metabolic data and 16S rRNA evolution was examined in several ways. First, we constructed a tree based on the metabolic matrix by using parsimony with unordered characters, and we compared this tree with the Bayesian 16S rRNA tree by using the Shimodaira-Hasegawa likelihood-based test implemented in PAUP. Second, we explored the number of metabolic state changes by optimizing the data regarding the presence or absence of growth on the set of Bayesian 16S rDNA trees with the aid of MacClade (D. Maddison and W. Maddison, 1993). Before we did this, we rearranged the branch order on the tree for the nodes that were not resolved (i.e., nodes defined an ancestor for more than two lineages) such that the total number of metabolic character state changes was minimized. In this way, we avoided including polytomies when tracing character evolution. Significance was assessed by comparing the observed number of changes with the number of changes with assumption of no correlation between carbon source growth data and phylogeny. Significance was established by optimizing the growth data on 1,000 random trees and comparing the observed number of changes with the distribution of changes for the random trees. Only phylogenetically informative characters were subjected to this analysis. Finally, we plotted the branch lengths that were estimated for the DNA and metabolic data separately by using maximum likelihood and parsimony, respectively. Significance was assessed by linear regression.
Nucleotide sequence accession numbers. Table 4 lists the GenBank accession numbers of our newly identified isolates.
|
View this table: [in a new window] |
TABLE 4. GenBank accession numbers of newly identified isolates
|
|
|
|---|
![]() View larger version (23K): [in a new window] |
FIG. 1. MPNs of cells per gram (dry weight) of soil from a Kobresia alpine meadow on three dates (month/day/year) in autumn and winter. Yellow cells were those showing visible yellow fluorescence. Error bars represent 95% confidence intervals.
|
![]() View larger version (22K): [in a new window] |
FIG. 2. Bayesian phylogenetic tree of 16S rDNA sequences of our 17 isolates and selected Pseudomonas (sensu stricto) sequences obtained from GenBank. Posterior probability values of 95% or greater are noted. The tree was rooted with P. flavescens, which gave the same topology as Escherichia coli. GenBank accession numbers: P. aeruginosa, Z76651; P. balearica, U26418; "P. borealis," AJ012712; P. cichorii, Z76658; P. flavescens, U01916; P. fluorescens, Z76662; P. mendocina, Z76664; P. putida, Z76667; and P. syringae, Z76669.
|
![]() View larger version (40K): [in a new window] |
FIG. 3. Frequency of growth on carbon sources. Seventeen cold-tolerant alpine isolates were tested for their ability to grow on 20 C sources. All isolates grew on at least one C source tested. Frequencies of growth on C sources ranged from 25 to 90%.
|
![]() View larger version (27K): [in a new window] |
FIG. 4. Frequency of utilization of carbon sources. Frequencies ranged from 30 to 100%. The phenolic category supported growth the least, and the amino acid and protein categories supported growth the most.
|
A comparison of the number of changes of carbon source preference on the inferred phylogeny with the number of changes required for random trees suggested that phylogeny provides a poor explanation for the evolution of carbon preference. For all of the carbon sources, the observed numbers of changes fell within the distribution that was expected if the character changes occurred randomly on the tree.
A bivariate plot of the branch lengths estimated for the DNA and carbon source data revealed a lack of correlation (r2 = 0; P = 0.97). Several lineages were noteworthy. Isolate SE22°1a was inferred to have undergone seven unique changes of carbon source preference since it last shared an ancestor with another strain about 0.003 substitution ago (0.3% change in sequence). At the other extreme, isolate WR7°2 underwent an approximately 1.2% change in sequence yet did not appear to change its preference for carbon from the inferred ancestral condition.
|
|
|---|
Bacterial genomes are known to undergo high rates of evolution by a combination of point mutations, deletions, gene duplications, and acquisitions of foreign DNA. The fate of mutations is governed by selection. In the absence of a particular carbon source, the residence time for the necessary gene is probably short. Our evidence for a dynamic and versatile metabolic repertoire suggests that alpine soil environments may be tremendously heterogeneous with respect to the availability of different carbon sources, and such differential selection due to the various concentrations of different carbon pools may drive the gains and losses of metabolic genes.
Our findings of a lack of congruence between phenotypic and genotypic data are in contrast to many of the classic nutritional studies of Pseudomonas. These nutritional studies, combined with rDNA hybridization data, have been seen as giving good agreement in intrageneric clustering (27-29). However, many of the same and similar reports state that gene exchanges between clusters may also explain some of their results (9, 27, 29, 30), implying less than total congruence between genotypic (rDNA) and phenotypic (nutritional) data. Many studies of Pseudomonas have presented genotypic and phenotypic data and commented on their congruence (15, 18, 34, 39), but few have rigorously tested such congruence. Our results, rigorously showing no congruence, suggest that statements of congruence should be based on careful tests of this issue and that the famous metabolic versatility of the genus might well be founded in large measure on the easy ability to acquire opportunistic metabolic capabilities.
Unique cold-soil clade.
Our phylogenetic analysis reveals that 7 of our 17 isolates (41%) fall into a novel, well-supported clade (posterior probability, 99%) that is closely related to the P. syringae lineage (25) (Fig. 2). This clade includes "P. borealis," an organism that was also isolated from cold soils. Although unpublished, the 16S rDNA sequence of "P. borealis" was deposited in GenBank after being isolated from Swedish tundra soil north of the arctic circle (M. Hokeberg, personal communication). An organism with 99.3% sequence identity has also been isolated from soil from Signy Island off Antarctica (B. Stallwood, personal communication). Since 41% of our isolates fall into this unique cold-soil clade, our data suggest that a significant proportion of the Pseudomonas isolates in our alpine soil are specifically adapted to cold soils. Also, the finding of very similar 16S isolates in widely separated, persistently cold climates suggests that this strain is not endemic but perhaps ubiquitous in extremely cold environments. This suggestion is in contrast to the results of a four-continent study of 38 mesophyllic soil strains of Pseudomonas sensu stricto which showed endemicity at distances of less than 197 km (5).
Significant role in the Kobresia ecosystem.
Our growth data alone strongly suggest that pseudomonads play a significant biogeochemical role in the Kobresia dry meadow community. All 17 of our isolates grow on levan (Fig. 4), a 2,6-linked polymer of fructose. This polysaccharide, thought to be important in the frost resistance of Kobresia, is specifically accumulated in significant quantities by Kobresia (T. Rosenstiel, personal communication). Since Kobresia dominates this community and is a significant source of litter and root-derived organic matter, we believe that levan is a major C source for the soil microbes here.
Several other interesting properties of alpine Pseudomonas emerged from our data. The isolates demonstrated a high use of maltose but a low use of starch (Fig. 4). Since alpine soil has significant amylase activity in both summer and winter (21), our results may indicate that Pseudomonas strains do not expend resources on excreting amylase but rather rely on the amylases of other organisms. More than 85% of the isolates grew on casein, a milk protein (Fig. 4), which supports other data suggesting that amino acids from the degradation of peptides account for most of the N cycled in the alpine (19). Also, in contrast to the low utilization of glycine by alpine microbes (20), our isolates have a high utilization of glycine (Fig. 4), suggesting that Pseudomonas strains may be important competitors with plants for amino acids, an important N source for plants in this ecosystem (33).
We found a high absolute and relative prevalence of fluorescent pseudomonads in our alpine soils, with the highest percentages present in autumn and the lowest percentage present in winter (Fig. 1). These data are in contrast to those of Mancinelli (24), who found Pseudomonas MPN counts to be 2 orders of magnitude lower and the lowest proportions to be present in autumn (24). These differences in results may represent year-to-year fluctuations, inadequate spatial sampling, or differences in technique. In any case, our results demonstrate large numbers of fluorescent pseudomonads in our study soils, suggesting that they play a prominent role in biogeochemical processes there, perhaps more actively in the autumn, when more substrate is available (6).
The interesting growth pattern of our isolates confirms earlier laboratory work on bacterial aromatic-compound degradation pathways. Many aromatics are converted to either protocatechuate or to catechol before ring cleavage (37). Of the phenolics we tested, earlier studies show that vanillate is converted to protocatechuate while salicylate is converted to catechol (37). Interestingly, the growth patterns of our isolates on vanillate and salicylate were mutually exclusive: those that grew on one did not grow on the other (Table 2), which is consistent with the laboratory work on vanillate and salicylate degradation.
Caveats.
Although we suspect that the observed variation in metabolic capacities among the strains reflects the characteristics of natural populations, the conditions used for the isolation of the strains may have enhanced the observed metabolic variation. Strains were isolated at five different times on seven different carbon sources and at three different temperatures. It is certainly true that different isolation conditions selected for different metabolic characteristics; however, it is also possible that specific metabolic pathways were lost due to mutation during isolation.
This work was supported by NSF grants MCB-0084223 and IBN-9817164.
|
|
|---|
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»