Below is the uncorrected machine-read text of this chapter, intended to provide our own search engines and external engines with highly rich, chapter-representative searchable text of each book. Because it is UNCORRECTED material, please consider the following text as a useful but insufficient proxy for the authoritative book pages.
7 Microbes on Mountainsides: Contrasting Elevational Patterns of Bacterial and Plant Diversity JESSICA A. BRYANT,* CHRISTINE LAMANNA,â HÃLÃNE MORLON,* ANDREW J. KERKHOFF,â¡ BRIAN J. ENQUIST,â Â§î© and JESSICA L. GREEN* The study of elevational diversity gradients dates back to the foundation of biogeography. Although elevational patterns of plant and animal diversity have been studied for centuries, such pat- terns have not been reported for microorganisms and remain poorly understood. Here, in an effort to assess the generality of elevational diversity patterns, we examined soil bacterial and plant diversity along an elevation gradient. To gain insight into the forces that structure these patterns, we adopted a multifaceted approach to incorporate information about the structure, diversity, and spatial turnover of montane communities in a phylogenetic context. We found that observed patterns of plant and bacterial diversity were fundamentally different. While bacterial taxon rich- ness and phylogenetic diversity decreased monotonically from the lowest to highest elevations, plants followed a unimodal pat- tern, with a peak in richness and phylogenetic diversity at mid- elevations. At all elevations bacterial communities had a tendency to be phylogenetically clustered, containing closely related taxa. In contrast, plant communities did not exhibit a uniform phyloge- netic structure across the gradient: they became more overdis- persed with increasing elevation, containing distantly related taxa. *Center for Ecology and Evolutionary Biology, University of Oregon, Eugene, OR 97403; â Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ 85721; â¡Departments of Biology and Mathematics, Kenyon College, Gambier, OH 43022; Â§Santa Fe Institute, Santa Fe, NM 87501; and î©Center for Applied Biodiversity, Science Conservation International, Arlington, VA 22202 127
128â /â Jessica A. Bryant et al. Finally, a metric of phylogenetic beta-diversity showed that bacte- rial lineages were not randomly distributed, but rather exhibited significant spatial structure across the gradient, whereas plant lineages did not exhibit a significant phylogenetic signal. Quantify- ing the influence of sample scale in intertaxonomic comparisons remains a challenge. Nevertheless, our findings suggest that the forces structuring microorganism and macroorganism communi- ties along elevational gradients differ. R oughly 250 years ago, Carolus Linnaeus (1781) documented how distinct plant and animal communities characterized the succes- sion of climatic zones along the slopes of mountains. Such eleva- tional gradients are characterized by dramatic changes in climate and biotic turnover over short geographic distances. The patterns observed by Linnaeus and his contemporaries played a foundational role in the development of ecology and biogeography (Briggs and Humphries, 2004). Studies of how individual taxa and community composition respond to elevational gradients have led to a search for generalized elevational pat- terns of biodiversity (Brown, 2001; Lomolino, 2001; McCain, 2005). These studies have documented elevational patterns of diversity across a wide variety of taxonomic groups, including trees, mammals, birds, reptiles, insects, and amphibians. In sum, this work has shown that taxa generally exhibit either monotonically decreasing or hump-shaped richness patterns with elevation (Stevens, 1992; Rahbek, 2005). However, despite a large number of proposed hypotheses to explain elevation patterns of diversity, their causes remain poorly understood. Improved knowledge of elevation gradients is fundamental to advancing basic ecology and predicting the potential consequences of climate change. Species in montane regions are often cited as being very sensitive to the impacts of warming (McDonald and Brown, 1992; Parmesan, 2006; Thuiller, 2007). Although elevational patterns of diversity for plants and animals are well established, we know very little about how microbial diversity varies across elevational gradients. This is a serious gap in our general under- standing of biodiversity, given that microbes are abundant and diverse, play a central role in ecosystem functioning, and will likely be an impor- tant component of ecosystem response to global warming (Rillig et al., 2002; Monson et al., 2006; Carney et al., 2007). Elevational diversity studies that consider empirical patterns of macroorganisms and microorganisms in parallel are needed to provide a more unified framework for under- standing diversity patterns in Earthâs major environmental gradients and predicting systemwide ecological responses to climatic change. Traditional elevational diversity studies have focused on how patterns of species richness, abundance, and range size change with altitude. These
Microbes on Mountainsidesâ /â 129 analyses have used a nomenclatural approach by focusing on species identities. However, the increasing availability of molecular phylogenies has renewed interest in using phylogenetic approaches to study the forces that influence patterns of biodiversity and biogeography [e.g., Chave et al. (2007)]. Because many species traits are generally conserved during the evolution of a lineage, one would expect a positive relationship between a measure of the phylogenetic relatedness of two species and a measure of their overall ecological similarity (phylogenetic niche conservatism) (Har- vey and Pagel, 1991). As a result, analysis of the degree of phylogenetic relatedness of taxa found within and across communities should provide insight into the ecological and evolutionary processes that organize these communities. Here, in an effort to assess the generality of elevational diversity patterns and the forces that structure these patterns, we quantified both plant and soil bacterial diversity patterns along an elevational gradient in the Colorado Rocky Mountains. A parsimonious hypothesis is that if the forces structuring biodiversity across the gradient are the same for bacte- ria and plants, then the resulting taxon and phylogenetic biogeographic patterns will be similar for both groups. Alternatively, if ecological and evolutionary processes along elevational gradients differ between the two groups (e.g., the taxa differ in their dispersal ability, response to envi- ronmental heterogeneity, interspecific interactions, or speciation rates), then we would expect them to be characterized by distinct patterns of diversity. To test these hypotheses, we adopted a multifaceted approach that examines diversity in the context of both ecological and evolutionary patterns. Therefore, in addition to the established convention of quantify- ing patterns of taxon richness and taxon turnover along the gradient [e.g., Whittaker (1960, 1967)], we examined several biodiversity measures that incorporate information about the phylogenetic structure, phylogenetic diversity, and phylogenetic turnover of plant and bacterial communities. ELEVATIONAL DIVERSITY IN A PHYLOGENETIC CONTEXT While the sampling methods and taxonomy used to quantify plant diversity are well established and standardized, microbial surveys vary greatly in their approach to characterizing diversity (Eisen, 2007). We determined the bacterial community composition of our soil samples by analyzing a PCR-amplified region of 16S ribosomal DNA, the most commonly used indicator of microbial biodiversity. Because bacteria are overwhelmingly diverse in soils, we chose PCR primers that narrowed our focal group to the phylum Acidobacteria. This subgroup of bacteria is diverse and ubiquitous in soils (Janssen, 2006) and thought to play an important role in biogeochemical cycling (Eichorst et al., 2007).
130â /â Jessica A. Bryant et al. We followed the classic approach to intertaxonomic diversity analy- sis by comparing patterns of species richness and phylotype richness of plants and bacteria, respectively, along the gradient. We also quantified the phylogenetic diversity of every sampled community by calculating the sum of the branch length in a phylogeny that connects all species in a community and the root (Faith, 1992b). Phylogenetic diversity is more inclusive than a simple count of species or types, in that it quantifies the evolutionary history of a group of taxa (Vane-Wright et al., 1991). Conser- vation biologists are interested in preserving phylogenetic diversity, as this is fundamental to maximizing evolutionary options for the future (Faith, 1994; Myers and Knoll, 2001; Sechrest et al., 2002; Forest et al., 2007). Phy- logenetic diversity is also believed to correspond to ââfeature diversity,ââ meaning the number of evolutionarily derived traits within a biological community (Faith, 1992b). In addition to measuring phylogenetic diversity, we quantified com- munity phylogenetic structure along the gradient by using two commonly used metrics: a mean pairwise distance metric sensitive to phylogeny- wide patterns [net relatedness index (NRI)] and a nearest-taxon-based measure sensitive to patterns at the ââtipsââ of the phylogeny [nearest taxon index (NTI)] (Webb et al., 2002). The degree of phylogenetic related- ness quantified by these metrics provides insight into drivers of commu- nity assembly. Assuming phylogenetic niche conservatism, phylogenetic clustering within a local assemblage is considered consistent with the hypothesis that selective filters (e.g., environmental conditions) cause local assemblages to comprise closely related taxa (Webb et al., 2002). Phylogenetic overdispersion, on the other hand, can be explained by two possible biotic interactions: competition (Webb et al., 2002) or facilitation (Lortie, 2007; Valiente-Banuet and Verdu, 2007). In the case of competition, more closely related species are hypothesized to compete more strongly with one another. This results in competitive exclusion, which leads to a community of distantly related species. In the case of facilitation, facilita- tor species are hypothesized to create microhabitats that permit distantly related species adapted to different environments to persist within a local assemblage. In addition to considering patterns in the diversity and phyloge- netic structure within communities along the elevation gradient (alpha- diversity), we investigated how community composition changes across a landscape (beta-diversity). Ecologists have long recognized that beta- diversity is important for understanding the biodiversity of montane ecosystems (Jaccard, 1912; Whittaker, 1960; Harte et al., 1999; Brehm et al., 2003; Mena and Vazequez-Dominguez, 2005). We examined beta- diversity in terms of compositional similarity, defined as the fraction of taxa shared between two samples (SÃ¸rensen index), and phylogenetic
Microbes on Mountainsidesâ /â 131 similarity, defined as the fraction of branch lengths shared between two samples. By analogy with the well-established distanceâdecay relation- ship, which describes the decrease in compositional similarity between two communities with increasing geographic distance (or equivalently elevational separation) between them (Soininen et al., 2007), we described the decrease in phylogenetic similarity with distance (phylogenetic dis- tanceâdecay). Our objective in exploring both measures of beta-diversity is to understand not only if there are shifts in compositional similarity with increasing elevational distance, as expected along an environmental gradient, but to quantify the phylogenetic nature of the shifts. Phylogenetic similarity reflects the combined additive influence of: (i) lineages that are shared between two communities that lead to shared taxa, and (ii) lineages that are shared but ultimately lead to unshared taxa. One can test whether the phylogenetic similarity between two communities is solely a consequence of compositional similarity, or if it is also caused by a nonrandom structure of shared and unshared lineages. A signifi- cant phylogenetic distanceâdecay pattern (i.e., one that differs from that expected by taxa turnover alone; see Materials and Methods) reflects signifi- cant spatial variability in lineage composition across a landscape. Based on the assumption of phylogenetic niche conservatism described above, changes in lineage composition should correspond to changes in the traits of species. Under this model, a significant phylogenetic distanceâdecay relationship should reflect strong variability in the ecologically relevant traits of biological communities across a landscape. RESULTS AND DISCUSSION Whereas bacterial richness decreased monotonically from the lowest to highest elevations, plant richness followed a unimodal pattern with a peak in species richness at midelevations (Fig. 7.1A). These contrasting diversity patterns emerged when richness values were calculated for bac- terial and plant samples individually and also when the samples for each respective group were pooled together at every elevational band [follow- ing the protocol suggested by Whittaker (1960)]. To our knowledge, an alti- tudinal richness pattern has never been reported for microorganisms. The patterns observed here for microbes and plants are consistent, respectively, with the classical monotonically decreasing and hump-shaped patterns observed across most macroorganism groups (Stevens, 1992; McCain, 2005; Rahbek, 2005). It has been argued that these two contrasting rich- ness patterns may emerge as a result of inconsistent sampling approaches among different studies, rather than an underlying ecological mechanism (Lomolino, 2001; Rahbek, 2005). By implementing a parallel sample design for the bacteria and plants, we controlled for two potential biases: varia-
132â /â Jessica A. Bryant et al. FIGURE 7.1â Variation in taxon richness (A) and phylogenetic diversity (B) across the elevation gradient. Data are presented as the fraction of total richness and phylogenetic diversity across the gradient. Solid symbols indicate sample rich- ness (core or quadrat), and open symbols indicate the pooled richness at each elevational site (n = 5 per site). At the sample level, Acidobacteria richness and phylogenetic diversity linearly decrease with elevation (regression analysis, r2 = 0.22, P < 0.05; r2 = 0.23, P < 0.05, respectively), whereas angiosperm richness and phylogenetic diversity patterns are hump-shaped (regression analysis, r2 = 0.53, PJBryant_Figure1_bw.tif respectively). Model choice was based on Akaike < 0.0005; r2 = 0.47, P < 0.005, information criteria. tion in sampling extent (i.e., the geographic distance between the furthest sampled elevations), which often occurs among studies, and sampling intensity (or effort) along the gradient within a taxonomic group, which often occurs within studies. As a result, the disparity in elevational rich- ness patterns observed between bacteria and plants is likely caused by differences in how ecological and evolutionary processes have operated across the gradient (although see below for a discussion on the potential influence of scaling effects). As expected, for both bacteria and plants (Forest et al., 2007) we found that the patterns of phylogenetic diversity mirrored those of taxon richness (Fig. 7.1B). However, a more detailed look at the phylogenetic structure of the bacterial and plant communities revealed another significant differ- ence. At all elevations bacterial communities had a tendency to be more phylogenetically clustered than expected by chance (Fig. 7.2). This obser- vation is consistent with results reported by Horner-Devine and Bohannan (2006) who found that bacterial communities in a wide range of environ- ments tended to be phylogenetically more closely related than expected
Microbes on Mountainsidesâ /â 133 FIGURE 7.2â Variation in community phylogenetic relatedness along the elevation gradient as measured with the NRI (A) and NTI (B). Positive index values indicate phylogenetic clustering, and negative values indicate phylogenetic overdisper- sion. Observed community phylogenetic structures unlikely to arise by chance (P < 0.05) are depicted by solid symbols. All microbial communities are clustered, zpq9990837750002.c.tif with >50% being significantly clustered. Angiosperm communities are not uni- formly clustered or dispersed across the gradient, but rather become increasingly overdispersed with increasing elevation. This trend in increased overdispersion with elevation is significant when measuring relatedness with the NRI (solid line; r2 = 0.70, P < 0.001). by chance. Given the parsimonious hypothesis that closely related taxa are more ecologically similar (i.e., phylogenetic niche conservatism), our results suggest that abiotic filtering tends to be a more prominent force in the structuring of bacterial communities along the gradient. Several studies have suggested that for most macroorganisms, ecological traits are phylogenetically conserved (Prinzing et al., 2001; Blomberg et al., 2003; Cavender-Bares et al., 2004). It is important to emphasize that although this statement may be correct for macrorganisms, the generality of niche conservatism for microorganisms and in particular bacteria is unknown. Observed phylogenetic clustering in microbial communities could also be the result of radiation events combined with dispersal limitation (Horner- Devine and Bohannan, 2006). As we discuss below, alternative explana- tions for the patterns we observe relate to the phylogenetic and spatial scale of our analyses (Swenson et al., 2006, 2007). Scaling issues are relevant to all of the biodiversity patterns we examined. In contrast to bacteria, plant communities did not show a uniform phylogenetic structure across the gradient. Plant communities tended
134â /â Jessica A. Bryant et al. to exhibit either random phylogenetic structure or phylogenetic over- dispersion. Surprisingly, our analyses indicated that plant communities also tended to become increasingly overdispersed at higher elevations (Fig. 7.2). Given niche conservatism, phylogenetic overdispersion is con- sistent with the importance of biotic forces (competitive exclusion or facili- tation) structuring community diversity. Recent experimental evidence suggests that both of these forces are important drivers in alpine plant community assembly, with a shift from competition at lower elevations, where conditions are less physically stressful, to facilitation at higher elevations where abiotic stress is high (Callaway et al., 2002). Increased overdispersion at high elevations suggests that the influence of facilitation on high-elevation communities is stronger than the influence of competi- tion at low elevations. An alternative explanation is that the evolution of traits necessary to cope with environmental conditions at high elevations has occurred independently in distantly related lineages (i.e., convergent evolution in high alpine plants) (Webb et al., 2002). This explanation goes against the assumption of phylogenetic niche conservatism. We observed that both plant and bacterial compositional similarity sig- nificantly decreased with elevational distance (Fig. 7.3). Plant and bacterial communities differed, however, in their phylogenetic distanceâdecay pat- terns. The bacterial phylogenetic distanceâdecay curve was significantly steeper than expected from the observed bacterial taxa turnover alone (Fig. 7.3A). In contrast, the plant phylogenetic distanceâdecay curve was not significantly different from expected from the observed decay in plant compositional similarity. These results are consistent with those reported above for the NRI and NTI measures of community phylogenetic structure, indicating that bacteria lineages were not randomly distributed across the elevation gradient. Rather, bacterial lineages exhibited a spatially struc- tured pattern across the gradient. Given the parsimonious hypothesis that closely related taxa are more ecologically (or functionally) similar, our observations suggest that bacterial lineages harbor increasingly disparate ecological features (or functions) at increased elevational distances as a probable consequence of abiotic filtering. These findings highlight the utility of gathering information on phylogenetic relationships between communities in montane regions as a means to quantify the potential consequences of selectively trimming evolutionary lineages under the scenario of mountaintop extinctions in response to global warming. Although our study was not designed to directly examine the environ- mental drivers of elevational diversity patterns, our results do illuminate their potential role in shaping biodiversity patterns across the gradient. The contrasting phylogenetic diversity patterns we observed in plants and microbes suggest a differing role in how abiotic forces structure communities across the gradient. Soil temperature and pH were consis-
Microbes on Mountainsidesâ /â 135 FIGURE 7.3â Compositional and phylogenetic similarity of Acidobacteria communi- ties (A) and angiosperm communities (B), as a function of the elevation separating the communities. The compositional (hollow black triangle) and phylogenetic (solid black circles) similarity for both angiosperm and Acidobacteria communi- ties significantly decrease with increasing elevational separation (Mantel test, P < 0.001). Lines represent best-fit regressions of similarity versus change in elevation (see Materials and Methods). The slope of the decay of phylogenetic similarity be- tween Acidobacteria communities is significantly steeper than predicted by a null model constrained by the decrease in taxon turnover (solid gray circles) (P < 0.05). The slope of the decay in phylogenetic similarity across the angiosperm communi- ties is not significantly different from the null prediction given species turnover. tently correlated with diversity in both plants and bacteria, and bacterial diversity was also strongly correlated with slope in a univariate analysis (Table 7.1). A multivariate analysis suggested that soil temperature was the major explanatory variable of taxon richness and phylogenetic diversity for both plants and bacteria (P < 0.001 in six of eight multivariate models). Turnover in taxon and phylogenetic composition of plant and bacterial communities was significantly correlated with changes in the majority of our measured environmental parameters (7.2); however, the combined influence of soil temperature, pH, and total nitrogen was the most impor- tant predictor for both groups. After controlling for these environmental parameters, geographic distance between samples significantly correlated with all turnover patterns (partial Mantel test, P < 0.001). These results imply that dispersal limitation could be occurring, but given the small geographic range of our gradient, they are more likely caused by effects of environmental heterogeneity that we did not characterize. The correlation of richness and turnover with temperature and pH is consistent with the
136 TABLE 7.1 Pearson Correlation Coefficients from Univariate Models Predicting Alpha Diversity Measures as a Function of Environmental Parameters % % % Diversity Type Elevation Temperature pH Slope Aspect Nitrogen Carbon Moisture Microbial richness â0.222* 0.174* 0.205* â0.239* 0.044 â0.002 0.0002 â0.148 Plant richness 0.533*** â0.502*** â0.189* 0.004 0.039 0.271 0.077 0.045 Microbial â0.228* 0.216* 0.132 â0.311* 0.030 â0.0001 0.116 â0.0009 phylogenetic diversity Plant phylogenetic 0.473** â0.443** â0.283* 0.069 0.002 0.229 0.022 0.056 diversity Microbial NRI 0.056 â0.063 â0.025 0.227 0.018 0.037 0.0001 0.002 Plant NRI â0.696*** 0.699*** 0.536*** â0.411** 0.132 0.094 0.056 0.001 Microbial NTI 0.101 â0.120 â0.054 0.048 â0.008 0.0004 â0.233 â0.026 Plant NTI â0.109 0.109 â0.015 â0.015 0.002 â0.059 â0.037 â0.004 For each environmental parameter both a linear and a quadratic model were fit. The best fit model was chosen using the lowest Akaike Informa- tion Criterion (AIC) value (see Materials and Methods). Quadratic models are indicated with italics, and a minus sign (â) is used to indicate the
TABLE 7.2â Mantel Correlations Between Community and Environmental Similarities % % % Elevation Temperature pH Slope Aspect Nitrogen Carbon Moisture Acidobacteria compositional â0.57** â0.57** â0.64*** â0.19* â0.10 â0.32** â0.17* 0.17* similaritya Acidobacteria phylogenetic â0.78 â0.76** 0.76** â0.40** â0.24* â0.18* â0.06 â0.04 similarityb Angiosperm compositional â0.69** â0.72** â0.44** â0.24* â0.26* â0.33* â0.25* â0.07 similaritya Angiosperm phylogenetic â0.65** â0.65** â0.37** â0.16 â0.19 â0.38* â0.27* â0.23* similarityb a b Quantified with SÃ¸rensen index. Quantified with PhyloSor. *P â¤ 0.05, **P â¤ 0.005. 137
138â /â Jessica A. Bryant et al. findings of other studies of plant (Allen et al., 2002; Hawkins et al., 2003; Currie et al., 2004) and microbial (Fierer and Jackson, 2006) diversity. It is well documented that the scale over which biodiversity is sampled will strongly influence observed patterns. For example, recent empirical studies have shown that decreasing the spatial grain at which organisms are sampled shifts their diversity patterns (Rahbek, 2005; Cavender- Bares et al., 2006; Slingsby and Verboom, 2006; Swenson et al., 2006, 2007). Although the spatial extent of our study was the same for bacteria and plants, the grain of our sample observations was different between these two groups. The spatial scales over which bacteria interact with each other are likely to be several orders of magnitude smaller than the scale at which they were sampled. Therefore, relative to plants, bacteria were likely sampled at a coarser grain, and thus we may have included a greater amount of environmental heterogeneity within a bacterial sample. Sampling bacteria at a spatial scale that more closely approaches the ââeco- logically equivalentââ grain of plants may result in convergent biodiversity patterns between these two groups. Taxonomic scale also influences biodiversity patterns. For example, taxonomic breadth, which defines how broadly or narrowly a target com- munity is defined from a phylogenetic perspective (e.g., bacteria versus Acidobacteria), can shift the degree of observed overdispersion or cluster- ing in that community (Swenson et al., 2006). Species are a natural taxo- nomic unit by which to measure plants (Mayr, 1942). Such an intuitive unit does not exist for prokaryotes. In this study we classified partial Acidobacteria 16S ribosomal DNA sequences into taxonomic units based on the commonly used 99% sequence similarity designation (see Materi- als and Methods). It is unknown how taxonomic resolution, defined as the threshold at which individuals are binned into taxonomic units, should influence phylogenetic patterns, although it has been shown to impact taxonomic patterns such as the taxaâarea relationship (Horner-Devine et al., 2004b). We found that binning bacteria into increasingly broader taxonomic units (i.e., 97%, 94%, and 90% sequence similarity) tended to dampen the strength of all observed elevational diversity patterns. How- ever, general trends did not qualitatively change (Figs. 7.4â7.6), suggesting that taxonomic resolution is not the cause of disparate bacterial and plant biodiversity patterns in this study. Alternative approaches to defining bac- terial taxonomic units such as ââecotypesââ (Cohan and Perry, 2007) could significantly change the results and lead to plant and microbial diversity patterns that more resemble one another. Differences in the approach to building the Acidobacteria and angio- sperm phylogenetic trees should also be considered when comparing phy- logenetic patterns between these two groups. The Acidobacteria phylogeny was estimated solely from molecular data identified in this study, whereas
Microbes on Mountainsidesâ /â 139 FIGURE 7.4â Variation in Acidobacteria (A) taxon richness and (B) phylogenetic di- versity across the elevation gradient at four different taxonomic resolutions. Taxo- nomic richness is the total number of taxa (phylotypes) and phylogenetic diversity is the minimum total branch length connecting all taxa in the community and the root. Richness and phylogenetic diversity are best modeled as linearly decreasing with elevation at all taxonomic resolutions (Akaike information criteria). Slopes of the taxon richness patterns: 99%, -0.006982*; 97%, -0.005158*; 94%, -0.001939; JBryant_SFig1_bw.tif 90%, -0.002364*. Slopes of the phylogenetic diversity pattern: 99%, -0.0007653*; 97%, -0.0004603*; 94%, -0.0002757*; 90%, -0.0001679. Asterisks indicate that the decrease in diversity with elevation is significant (P < 0.05) based on linear regres- sion analysis. the angiosperm phylogenic tree topology was constructed by using the widely accepted supertree approach (Webb and Donoghue, 2005), and branch lengths were assigned based on estimates of the minimum age of internal nodes (see Materials and Methods). Comparative analyses using molecular approaches alone for both plants and microbes would improve our confidence in such phylogenetic comparisons. Such approaches will be facilitated in the future by increased accessibility to molecular data. Microorganisms (especially prokaryotes) are very diverse in soils (Torsvik et al., 1990; Janssen, 2006). On par with most microbial diversity studies, it is likely that we sampled the most abundant taxa in each soil core along the elevational gradient. Sampling effort (i.e., the proportion of a community that is sampled) is known to significantly influence taxonomic biodiversity patterns (Plotkin et al., 2002; Woodcock et al., 2006; Green and Plotkin, 2007; Morlon et al., 2008). To our knowledge the influence of sampling effort on phylogenetic biodiversity patterns has not been explored. For example, estimators are available to predict the taxon richness (Hughes et al., 2001) and taxon similarity (Chao et al.,
140â /â Jessica A. Bryant et al. FIGURE 7.5â Microbial NRI and NTI values across the elevational gradient at four different taxonomic resolutions. Observed community phylogenetic structures chance (P < 0.05) are depicted with solid circles. unlikely to arise by JBryant_SFig2_bw.tif
Microbes on Mountainsidesâ /â 141 FIGURE 7.6â (A) Compositional and (B) phylogenetic similarity of Acidobacteria communities at four taxonomic resolutions. The compositional and phylogenetic similarity between Acidobacteria communities significantly decreases with increas- ing elevational separation (Mantel test, P < 0.001) at all taxonomic resolutions presented. The slope of the decay of phylogenetic similarity between Acidobacteria communities is significantly steeper than predicted by a null model constrained by the decrease in taxon similarity (P < 0.05). As in Figure 7.3, lines represent an exponential model fit to the observed data. JBryant_SFig3_bw.tif 2005) of ecological communities from sample data, but there are currently no estimators to predict phylogenetic richness, phylogenetic structure, or phylogenetic turnover from sample data. A new generation of estimators is needed for future comprehensive studies that examine taxonomic and phylogenetic diversity patterns in parallel. As discussed by others, a promising approach to understanding eleva- tional diversity patterns (Rahbek, 2005), and more generally biodiver- sity patterns (Green and Bohannan, 2006), is to conduct intertaxonomic comparisons to elucidate the spatial and taxonomic scales and degree of sampling effort over which microbial biodiversity relationships approach those of macroorganisms. Such an approach is ambitious, but increasingly tractable as molecular approaches advance our ability to comprehensively characterize biodiversity. Here, we have shown that across an elevation gradient, plant and microbial communities exhibit different patterns of diversity. Phylogenetic-based analyses suggest that the evolutionary and ecological processes driving the biogeographic patterns may differ signifi- cantly between these two domains of life. Further work is needed to link the phylogenetic patterns to functional differences among plant and bacte- rial taxa. Such comparative analyses are needed to provide the empirical
142â /â Jessica A. Bryant et al. foundation for a truly inclusive and predictive theory linking patterns of biodiversity to ecosystem function. MATERIALS AND METHODS Study Site and Sampling We sampled angiosperm and Acidobacteria communities at five sites along an elevational transect located near the Rocky Mountain Biologi- cal Laboratory, Gunnison County, Colorado. The sites extend from 2,460 to 3,380 m above sea level and spanned a geographic distance of 39 km. Within each study site, we placed five 1-m by 1-m quadrats in a subtran- sect running down the slope. Three soil samples, separated by 1 cm, were collected adjacent to the middle, highest, and lowest quadrats at each site (nine total). All soil samples were collected from the B-horizon by using sterile glass collection jars. After collection, the soil samples were homog- enized and stored at â80Â°C until analysis. Average soil temperature at each site was measured by placing Hobo Temperature Data Loggers (OnSet) at 10-cm depth in relatively open patches and recording soil temperature every hour for the month of July in 2007. Total carbon and nitrogen in the soil samples were measured by using a Costech ECS 4010 CHNS-O system. Soil pH was measured after shaking a soil water (1:3 wt/vol) suspension for 30 min. Soil moisture was measured gravimetrically. ArcGIS data and area photos were used to cal- culate slope and aspect of the sites. These data had a 15.24-cm resolution per pixel and a horizontal accuracy of 60â90 cm. Characterization of Acidobacteria Communities with 16S Clone Libraries At each site, the bacterial communities within the three soil samples collected adjacent to the middle quadrat and one soil sample adjacent to the lower and upper quadrats were characterized by using sequence analysis of clone libraries (five total). DNA was extracted by using Mobio Power Soil DNA Isolation kits (MoBio Laboratories). Triplicate PCRs were carried out on each soil extraction by using the Acidobacteria-specific PCR primer set Acid31/Eub518 (Barns et al., 1999). The 25-Î¼l PCR mixtures were composed of 10 Î¼l of 5 Prime MasterMix (5 Prime, Inc.), 14 Î¼l of water, 1 Î¼mol of each primer, and 1 Î¼l of DNA extract. The PCR conditions used were as follows: 3 min at 94Â°C, 25 cycles of 30 s at 94Â°C, 30 s at 50Â°C (Fierer et al., 2005b), 60 s at 72Â°C, and a final extension for 5 min at 72Â°C. Triplicate PCRs were pooled and then gel-purified by using a Min- Elute PCR Purification kit (Qiagen). Amplicons were ligated into pCR4-
Microbes on Mountainsidesâ /â 143 TOPO vectors and cloned by using a TOPO-TA cloning kit (Invitrogen). Ninety-six clones from each soil sample were selected for sequencing. Plasmid purification and sequencing of cloned PCR products was done at the Qiagen Genomic Services/Sequencing facility with an ABI 377 or 377xl sequencer (Applied Biosystems). A total of 2,239 cloned 16s sequences were aligned with the NAST alignment tool (DeSantis et al., 2006b), and the alignments were manu- ally edited based on conserved primary sequence and secondary struc- ture information in the ARB software package (Ludwig et al., 2004). The phylum affiliation of each sequence was checked by using the BLAST tool within the National Center for Biotechology Information (Altschul et al., 1990). Potentially chimeric sequences were identified by using the Bellerophon server (Huber et al., 2004). Putative chimeric sequences were manually assessed by building trees in ARB that contained a set of refer- ence sequences obtained from the Greengenes database (DeSantis et al., 2006a) and the 5â² and 3â² sides of the putative chimeras. Sequences were removed from the analysis that had 5â² and 3â² ends affiliating with different groups of reference sequences in the tree (Horner-Devine et al., 2004a). There is no standard definition of microbial species. Therefore we grouped our 2,196 nonchimeric Acidobacteria sequences into phylotypes with a <99% sequence similarity cutoff by using the programs PHYLIP (Felsenstein, 1989) and DOTUR (Schloss and Handelsman, 2005). This is a commonly used phylotype designation (Kroes et al., 1999), which pro- vides high phylogenetic resolution. One sequence was randomly chosen to represent each phylotype. The representative sequences were used to build a phylogenetic tree by maximum-likelihood methods using the program phyML (Guindon and Gascuel, 2003).We used JukesâCantor and gamma substitution models where the gamma distribution parameter was estimated from the data. Only informative base positions were used to bin sequences into phylotypes and build the microbial phylogenetic tree. All diversity analyses were later repeated by using 97%, 94%, and 90% sequence similarity cutoffs (Figs. 7.4â7.6). Characterization of Plant Communities Angiosperms within each quadrat were identified to species level and checked against Rocky Mountain Biological Laboratory (RMBL) Her- barium specimens jointly by B.J.E., A.J.K., and C.L. Vouchers are being prepared for deposition in the RMBL Herbarium and the University of Arizona Herbarium. All plants were identified in 2005, except for the plants at the lowest site, which were sampled in 2006. Plots were sampled near the peak of the growing season, and thus some individuals with later phenologies could not be identified to species. We staggered the plant
144â /â Jessica A. Bryant et al. sampling dates with the aim of sampling each community at the same rela- tive phenological time point. Any individuals that could not be identified to species or differentiated from known species were excluded from the analysis. This affected between 10% and 15% of the possible species at each site. We used version R20031202 of Phylomatic to construct a tree topology consisting of all of the angiosperms identified in all our quadrats, based on the Angiosperm Phylogeny Group (APG) II backbone (Angiosperm Phy- logeny Group II, 2003; Webb and Donoghue, 2005) and used results from recent plant cladistics studies to resolve polytomies (Worberg et al., 2007; Winkworth et al., 2008). The final tree we used for our analyses was almost completely resolved to the family level. We assigned branch lengths to the tree by using the Phylocom module BLADJ to constrain the internal nodes with available age estimates (Wikstrom et al., 2001) and interpolated the other nodes for which direct age estimates are not available. Diversity Analyses Taxon Richness and Phylogenetic Diversity We define the term community as all phylotypes originating from a single soil core (bacteria) or species identified in a single quadrat (plants). Taxon richness within each community was quantified as the total number of species or phylotypes within that community. Phylogenetic diversity within each community was quantified as the minimum total branch length connecting all species within the community to the root of the phylogenetic tree (Faith, 1992b). Phylogenetic diversity was calculated by using the pd module within Phylocom-3.40 (by C. O. Webb, D. D. Ackerly, and S. W. Kemble; available at http://phylodiversity.net/phylocom/). We used a rarefaction sampling approach to account for the unequal sample sizes of each microbial community (number of clones) by calculating the mean of the taxon richness and phylogenetic diversity of 1,000 random- ized subsamples of each community. Each community was subsampled by the number of clones in the smallest library (75 clones). Phylogenetic Community Structure Using the classical NRI and NTI (Webb, 2000; Webb et al., 2002), we measured the extent to which co-occurring species in a community are phylogenetically related compared with what is expected by chance. With both indices, the phylogenetic structure of the observed community was compared to a null expectation obtained by randomly sampling the pool of all of the species identified in the study 1,000 times, while constraining both the number of taxa in the community and species occurrence across
Microbes on Mountainsidesâ /â 145 communities (Kembel and Hubbell, 2006). Observed values smaller or larger than 975 of the randomizations were considered significantly struc- tured (P < 0.05). Compositional and Phylogenetic Similarity Compositional similarity between all pairwise comparisons of com- munities was quantified with the SÃ¸rensen Index: Sij Sorij = , 1 (Si + Sj ) 2 where Sij is the number of taxa common to both communities i and j, and Si and Sj are the total number of species found in community i and j, respectively (Krebs, 1998). By analogy, phylogenetic similarity between two communities was quantified by using an index, coined PhyloSor: BLij PhyloSorij = . 1 (BLi + BL j ) 2 Here, BLij is the branch length common to both communities i and j, and BLi and BLj are the total branch lengths of community i and j, respectively. The PhyloSor index ranges from indefinably close to 0 (two communi- ties only share a very small root) to 1 (both communities are composed of the same taxa). Similar approaches have been carried out by Lozupone and Knight (2005) and Ferrier et al. (2007) when considering the closely related Jaccard and Bray-Curtis similarity indices. Using PhyloSor, one can test whether two communities are phylogenetically more or less similar than what is expected given their taxa similarity. This is done by comparing the phylogenetic similarity of the observed communities to a null expectation obtained by randomly sampling the pool of all of the species identified in the study while constraining the number of taxa in each community and the number of taxa shared by the two communities. By analogy with the well-established distanceâdecay relationship, which describes the decrease in compositional similarity between two communities with increasing geographic distance (or equivalently eleva- tional separation) between them (Soininen et al., 2007), we used PhyloSor to quantify the decrease in phylogenetic similarity with distance (phy- logenetic distanceâdecay). We tested whether the slope in the decay of phylogenetic similarity was greater or less than what was expected given the taxonomic decay in similarity by comparing the observed slope with a distribution of distanceâdecay slopes obtained by randomizing the loca-
146â /â Jessica A. Bryant et al. tion of taxa at the tips of the community phylogenetic tree 1,000 times. This is equivalent to randomly sampling the taxa while constraining the number of taxa in each community, the number of taxa shared by any two communities, and taxa occurrence across all communities. The observed phylogenetic slope was assumed significantly different from the null if it was greater than or less than 975 of the slopes of the randomizations (two-tailed test, P < 0.05). To address the unequal sampling of microbial communities, we repeatedly calculated and tested the significance of distanceâdecay slope on subsampled communities, where communities were subsampled by the number of clones in the smallest library. The significance of results did not differ between repeated subsamples. Linking Diversity Measurements to Environmental Parameters To determine the correlation between elevation and potential environ- mental drivers with the observed diversity patterns, we used polynomial regression analysis. For each environmental variable, we fit a linear and a quadratic regression model. The best model was determined based on Akaikes Information Criterion differences (Burnham and Anderson, 2002). Using a stepwise regression to select variables and interactions, a multi- variate model was constructed for each alpha diversity measurement. Mantel tests (999 permutations) were used to determine whether compositional and phylogenetic similarity decayed significantly with elevational distance (Legendre and Legendre, 1998). Similarity values between pairwise comparisons of microbial communities were the aver- ages of 1,000 rarefaction samples, as described above. The best fit and the most homoscedastic residuals were found in models that used the log transformation of similarity against elevational distance, with the exception of angiosperm taxa similarity, which was best described by a linearâlinear distanceâdecay model. We used Mantel tests to examine cor- relations between community similarity and environmental similarity [for a discussion of these methods, see Legendre et al. (2005) and Toumisto and Ruokolainen (2006)]. We chose the combination of environmental variables that best explained the changes in angiosperm and Acidobacteria commu- nity composition between samples with BIO-ENV (Clark and Ainsworth, 1993) and tested the importance of these variables after controlling for geographic distance and vice versa by using partial Mantel tests. For all analyses, moisture, carbon, and nitrogen were arcsin(sqrt(y))-transformed and aspect was 1/y-transformed (Legendre and Legendre, 1998).
Microbes on Mountainsidesâ /â 147 ACKNOWLEDGMENTS We thank our field crew from the University of Arizona (Laura Crum- bacher, Peter Gaube, Melissa Wilson, James Stegen, Jason Pither, and Robin Sleith) for their help collecting the plant data and soil cores; the staff and researchers of the Rocky Mountain Biological Laboratory, in particular John Harte and Ian Billick for their support; the staff and fac- ulty at the University of California, Merced, for their support; Brendan Bohannan, Kathryn Docherty, Michael Donoghue, Rebecca Mueller, James OâDwyer, Nathan Swenson, Ian Wright, and an anonymous reviewer for useful comments; and John Avise, Francisco Ayala, and Stephen Hubbell for organizing the Sackler Colloquium. Logistical support was provided by National Science Foundation Grant DBI 0420919 (to Ian Billick). J.A.B., H.M., and J.L.G. were supported by National Science Foundation Grant MCB 0500124 (to J.L.G.). B.J.E., A.J.K., and C.L. were supported by a National Science Foundation CAREER award, the University of Arizona, and a Los Alamos (Department of Energy) Laboratory Directed Research and Development award (to B.J.E.).