Using Phylogenetic, Functional and Trait Diversity to Understand Patterns of Plant Community Productivity

Abstract

Background

Two decades of research showing that increasing plant diversity results in greater community productivity has been predicated on greater functional diversity allowing access to more of the total available resources. Thus, understanding phenotypic attributes that allow species to partition resources is fundamentally important to explaining diversity-productivity relationships.

Methodology/Principal Findings

Here we use data from a long-term experiment (Cedar Creek, MN) and compare the extent to which productivity is explained by seven types of community metrics of functional variation: 1) species richness, 2) variation in 10 individual traits, 3) functional group richness, 4) a distance-based measure of functional diversity, 5) a hierarchical multivariate clustering method, 6) a nonmetric multidimensional scaling approach, and 7) a phylogenetic diversity measure, summing phylogenetic branch lengths connecting community members together and may be a surrogate for ecological differences. Although most of these diversity measures provided significant explanations of variation in productivity, the presence of a nitrogen fixer and phylogenetic diversity were the two best explanatory variables. Further, a statistical model that included the presence of a nitrogen fixer, seed weight and phylogenetic diversity was a better explanation of community productivity than other models.

Conclusions

Evolutionary relationships among species appear to explain patterns of grassland productivity. Further, these results reveal that functional differences among species involve a complex suite of traits and that perhaps phylogenetic relationships provide a better measure of the diversity among species that contributes to productivity than individual or small groups of traits.

About the Authors:

Marc W. Cadotte

* E-mail: [email protected]

Affiliations National Center for Ecological Analysis and Synthesis, University of California Santa Barbara, Santa Barbara, California, United States of America, Department of Biological Sciences, University of Toronto - Scarborough, Scarborough, Ontario, Canada

Jeannine Cavender-Bares

Affiliation: Department of Ecology, Evolution and Behavior, University of Minnesota, St. Paul, Minnesota, United States of America

David Tilman

Affiliation: Department of Ecology, Evolution and Behavior, University of Minnesota, St. Paul, Minnesota, United States of America

Todd H. Oakley

Affiliation: Department of Ecology, Evolution and Marine Biology, University of California Santa Barbara, Santa Barbara, California, United States of America

Introduction

For nearly two decades, researchers have tested the prediction that community productivity is positively related to plant diversity [1]–[6]. One guiding assumption has been that greater diversity in functional traits allows species to access more of the total resources [7]–[9], whether they be nutrients, water, pollinators or fungal symbionts, and allows multiple competing species to coexist [10]. Researchers have advocated measuring diversity in functional attributes relevant to those critical limiting resources and assumed that this should be the best predictor of community productivity and ecosystem functioning [11]–[13]. However, identifying the critical resources over pertinent temporal and spatial scales as well as the most relevant functional traits can be challenging.

The first approximation to classifying critical functional differences has been to group species into functional groups, which is often based on broad morphological and physiological similarities (e.g., C4, C3, legumes, etc.). The richness of functional groups has been used to potentially explain variation in community productivity [7], [13], [14]. But functional group richness is a problematic measure for two reasons. First, the removal or addition of “functionally redundant” species may have effects on community dynamics and processes [15]–[17], indicating that there are important functional differences not captured by broad groupings. Competition theory suggests there must be some niche differentiation, no matter how similar species are to stably coexist. The second reason is that functional group richness tends to predict only a limited amount of variation in productivity [18] and may even explain less variation than having randomly assigned groups [19].

Given that functional groups may be inadequate representations of critical functional diversity, ecologists have sought other ways of measuring functional diversity by measuring specific a priori selected traits [e.g.], [20], [21,22]. In contrast to measuring specific traits, other measures of functional diversity use multivariate techniques to evaluate trait differences/similarities among species without reliance on a small number of traits selected a priori. The first multivariate strategy, functional attribute diversity (FAD), introduced by Walker and colleagues [23], sums species distances in trait space as a measure of total trait or functional dissimilarity. The second strategy, functional diversity (FD) by Petchey and Gaston [24] essentially uses the FAD trait distance matrix to create a functional dendrogram from a clustering routine. FD then corresponds to the total dendrogram branch lengths connecting community members together. In a comparison of FAD and FD, Petchey et al. [25] show that FD better explains variation in community biomass accumulation. The final multivariate strategy we introduce here is a variation of the FD strategy. Since FD and FAD do not account for groups of correlated traits, we use nonmetric multidimensional scaling (NMDS) to create the distance matrix that accounts for correlated traits. Then we perform the FD clustering to produce branch lengths.

All these multivariate techniques, while relaxed from strict a priori trait decisions, may still be sensitive to which traits are included in the analyses [25]. Further, multivariate distances reflect the magnitude of scale units used (e.g., cm vs. m) and the variation in scale for different traits as opposed to the actual difference in ecological function (e.g., does a 10% difference in leaf size have the same ecological consequences as a 10% difference in seed mass?). Given these potential limitations of trait-based approaches, we advocate the use of another metric, phylogenetic diversity (PD) –that is the sum of phylogenetic branches connecting species together. If phenotypic dissimilarity is correlated with evolutionary divergence times [26]–[28], then the more divergent two species are, the greater likelihood that they differ ecologically. As a diversity measure, PD has been shown for some datasets to better explain variation in community productivity than species or functional group richness [9], [18].

In this paper, we compare the efficacy of explaining variation in community productivity with seven different types of trait diversity measures: 1) species richness, 2) variation in 10 individual traits, 3) functional group richness, 4) Walker and colleagues' (1999) Functional Attribute Diversity (FAD), 5) Petchey and Gaston's (2002) Functional Diversity (FD), 6) our Nonmetric Multi-Dimensional Scaling (NMDS) approach, and 7) community PD. Our goal is to find which of these various diversity measures provides the best possible explanation of patterns of community productivity, moving from simple, single-variable models to metrics that represent full community trait differences.

Methods

Study site

In 1993, the vegetation and seed bank in a post-abandonment agricultural field located at Cedar Creek Natural History Area, Minnesota, USA, were removed via herbicide, burn and bulldozing treatments. The following year 13×13 m plots were seeded with 1, 2, 4, 8 or 16 grassland savanna species (experiment 120). Plot composition was randomly chosen from a pool of 18 species that included four C3 grasses, C4 grasses, legumes, non-legume herbaceous forbs and two woody species. At each level of diversity 28–35 replicates were established, and plot composition was maintained by manually weeding and annual burns [for more details see: http://www.cedarcreek.umn.edu/research/exper/e120/, [2], [13], [29]]. In 1995 three more species were added to substitute for poorly germinating species from the original 18, but subsequent weeding did not target the poor germinators meaning that 21 species were actually included in this experiment (species are identified in Fig. 1).

[Figure omitted. See PDF.]

Figure 1. Three dendrograms representing relationships among species.

The first is based on maximum likelihood analyses of genetic sequences from four genes. The second uses the functional diversity methodology of Petchey and Gaston (2002) on all measured traits. The third dendrogram also uses the functional diversity method on three orthogonal dimensions from nonmetric multi-dimensional scaling.

https://doi.org/10.1371/journal.pone.0005695.g001

As an estimate of productivity, peak aboveground biomass was measured by clipping, drying and weighing four 0.1 m by 3.0 m strips per plot. Biomass was sampled annually from 1996 to 2007 and we here use the average in biomass production as our dependent variable. The long-term average biomass was used because inter-annual variation was due largely to subsampling variance and regional climatic variation [30], [31]. In plots that carried fire the oaks (Quercus ellipsoidalis and Q. macrocarpa) performed poorly –contributing little to productivity, and so they were excluded from the calculation of the trait, functional and phylogenetic diversity metrics.

Phylogeny construction

We constructed a phylogeny for the species used in this experiment and a second biodiversity experiment at Cedar Creek (i.e., Experiments 120 and 123), which included a total of 31 species (see Appendix S1 for a list of all species). In February, 2008, for each of the 31 species, we searched GenBank [32] for four gene sequences commonly used in published angiosperm phylogenies : matK, rbcl, ITS1 and 5.8s. Of the 31 species, 14 had at least one gene represented in Genbank and for a further 16 species, we used gene sequences from a congeneric relative not included in these experiments. Collectively, the species used in this experiment represent many of the deep historical angiosperm bifurcations, relative to the number of branches connecting close relatives. Therefore, the effect on branch length estimates from using congeneric species is likely minimal, so long as congenerics are monophyletic with the species used in the experiment. We also included two representatives of early diverging angiosperm lineages as outgroup species, Amborella trichopoda and Magnolia grandiflora, and we added 4 other species, each represented by multiple genes, that were included in clades whose members did not have overlapping sequences (see Appendix S1). For these 36 species we aligned sequences using MUSCLE [33]. We then selected best-fit maximum likelihood models of nucleotide substitution for each gene using the Akaike Information Criterion, as implemented in Modeltest [34], [35].

Using the aligned sequences and the best-fit models of nucleotide substitution, we estimated a maximum likelihood phylogeny using the PHYML algorithm with a BIONJ starting tree [36], [37]. To assess nodal support on maximum likelihood phylogenies, we report Approximate Likelihood Ratio Test (aLRT) scores, which have been shown to correlate with ML bootstrap scores, but require much less computational time [37]. The maximum likelihood tree is available in Appendix S1. A single species that lacked any genetic data, Rudbeckia hirta, was added as a polytomy with Liatris aspera and Coreopsis palmata because they are all considered members of the Asteroideae subfamily (see Appendix S1). For the analyses in the present paper, we pruned out the 10 species not used in experiment 120 (Fig. 1 shows the pruned phylogeny and Figure 1 in Appendix S1 shows the full phylogeny for both experiments). We did not rerun PHYML on the subtree members due to the sparseness of the gene matrix, especially for the Asteraceae species.

From the phylogeny, we calculated phylogenetic diversity (PD) for each experimental plot as the total phylogenetic branch lengths connecting only the community members together not including the root of the larger phylogeny [18]. Here we are using a single method of phylogenetic construction, but there are other methods that may alter PD estimates. However, recent analyses have shown that the method of phylogenetic construction does not appear to alter qualitative results [18], [38].

Trait data

In the summers of 2007 and 2008, we measured leaf traits in the Cedar Creek biodiversity experiment [E120]. We sampled three fully mature leaves from ten individuals of each species collected within the maintained experimental plots as well as from the unmaintained experimental plots. Each individual was identified from a randomly chosen plot to cover the range of diversity treatments. We scanned fresh leaves on a flatbed scanner on the same day as collection with petioles and sheaths removed. Leaf area, perimeter and Feret's diameter (i.e., parallel lines touching opposites ends of the leaf) were calculated from the scanned leaves using ImageJ software [39]. These measures allowed calculation of perimeter per area (P/A, cm⋅cm−2), which is empirically correlated with leaf hydraulic conductance across a wide range of taxa [40]. Perimeter per leaf area⋅Feret's diameter is a unitless measure of leaf lobedness [e.g., 41] that influences the leaf radiation balance [42]. After scanning, we dried the leaves at 65°C for three days and weighed to calculate specific leaf area (cm2⋅g−1).

We determined seed mass by collecting seed heads for ten plants per species after seeds were fully mature, air drying the seeds, and then weighing together ten seeds (and dividing by ten) to calculate a mean seed mass per plant. For five species that were not seeding, seed mass was taken from online commercial databases, including the Native Seed Network (www.nativeseednetwork.org), Wildflower Farm Inc. (www.wildflowerfarm.com), and Prairie Moon Nursery (www.prairiemoon.com). Plant height was taken from the USDA Plants Database (plants.usda.gov). See table 1 for trait codes.

[Figure omitted. See PDF.]

Table 1. Diversity variables calculated.

https://doi.org/10.1371/journal.pone.0005695.t001

Trait variation at the plot scale was estimated by the coefficient of variation in trait values for leaf area, leaf perimeter area ratio, leaf lobedness, specific leaf area, seed weight and height, while plot presence/absence was recorded for the C3, C4, forb and nitrogen fixer functional groups. We calculated four trait diversity metrics: the number of functional groups, FG; functional diversity, FD [24]; functional attribute diversity, FAD [23]; and Nonmetric Multi Dimensional Scaling, NMDS [e.g., 43]. FD estimates net species similarity or differences as branch lengths from a functional dendrogram based on a multivariate distance matrix [24]. To calculate FD, we scaled the traits to have a mean of zero and variance of one. We then calculated a Euclidean distance matrix and performed hierarchical clustering on this matrix and calculated FD as the total branch lengths connecting community members together [24]. For FAD, we again used the traits scaled to mean of zero and variance of 1 and calculated a Euclidean distance matrix. We then summed the distances for all species in a community [23]. Finally, since we were including multiple traits in these analyses, we wanted to account for groups of correlated traits. We performed NMDS on the trait matrix including functional groups (Fig. 2). We ran the analysis using two to five dimensions and choose the number of dimensions that reduced stress and which had no deviations on a dissimilarity-distance plot [44]. We selected a three-dimension model due to the stress reduction associate with this model (i.e., 8.71 versus 16.68 for the two-dimension model). We then calculated FD using the dimensions as independent traits. All variables are defined in Table 1. We have included the script, written in R 2.7.1 (www.R-project.org), to calculate PD, FD and FAD from a community membership list, a phylogeny and a trait matrix (see Appendix S2).

[Figure omitted. See PDF.]

Figure 2. The ordination plot produced by nonmetric multi-dimensional scaling.

Symbols refer to functional group membership.

https://doi.org/10.1371/journal.pone.0005695.g002

Statistical analysis

There were a total of 16 different diversity variables used in this analysis (Table 1), and the goal was to construct the most informative model explaining patterns of biomass production. To narrow down the number of potential explanatory variables, we searched for the best single-variable models where average annual productivity was regressed against each diversity metric and compared to the explanatory ability of these models using Akaike weights, which can be interpreted as the probability that a particular model is the best fit to the observed data among a set of candidate models [45]. We further used mallow's Cp to rank these single variable models. We checked diagnostic plots (e.g., residual versus fitted plots) for potential outliers and data trends.

We were not only interested in the best single variable explaining patterns of productivity, but also combining PD and trait metrics in multi-variable models. We performed a stratified all subsets model approach where PD was included in models with either single variable trait diversity metrics or multivariate functional diversity metrics, but not both since the multivariate metrics are based on the single variables. Mallow's Cp was used to select potential models up to 5 variables and Akaike weights to compare them. Alternatively, we selected the best single variables and constructed models around these, comparing them using Akaike weights.

We also asked whether any of the continuous individual traits themselves had significant phylogenetic signals. To do this we generated phylogenetically independent contrasts [28] for the various traits. We compared the summed absolute node contrast values to that expected from 1000 randomizations (mean and 95% confidence intervals) [46]. If the observed summed contrast is significantly lower than the randomly generated values, sister clades tend to be more similar to one another than random, and there is a detectable phylogenetic signal. All analyses were done using R 2.7.1 (www.R-project.org)

Results

Interestingly, the species in the FD dendrogram did not produce the same sister pairings as in the molecular phylogeny (Fig. 1), while the NMDS dendrogram shows species clustering more similar to the molecular phylogeny (Fig. 1). As would be expected because the number of branches, and thus the sum of all branch lengths, increases with the number of plant species (N), PD and N were highly correlated (Fig. 3). However, PD seems to better explain variation in productivity compared with N, FG, FD, FAD or NMDS (Table 2). Productivity was positively related to each of these measures (Fig. 4). The better fit provided by PD compared to the traditional variable, N, suggests that PD, which combines the effect of differences (via branch lengths) with the effect of the number of species (as N-1 branches) is a superior single measure of diversity than N.

[Figure omitted. See PDF.]

Figure 3. The relationship between species number and phylogenetic diversity.

https://doi.org/10.1371/journal.pone.0005695.g003

[Figure omitted. See PDF.]

Figure 4. The relationship between average annual plot productivity and six diversity metrics.

Of these six metrics, PD is the best single explanatory variable, second only to the presence of a nitrogen fixer (see Table 2).

https://doi.org/10.1371/journal.pone.0005695.g004

[Figure omitted. See PDF.]

Table 2. Results of univariate models.

https://doi.org/10.1371/journal.pone.0005695.t002

However, the single best explanatory variable was the presence of a nitrogen fixer (Table 2). Only within the plots with a nitrogen fixer present, which includes most plots, PD is still significantly related to productivity (F1,103 = 19.98, P<0.0001) and a better model (Akaike weight, AW = 0.634) than N, FAD, NMDS, FD or FG (AW = 0.142, 0.089, 0.084, 0.047 and 0.037, respectively).

While the above analyses informs our understanding of how individual models compare to one another, we were also interested in combining variables in models to better explain variation in productivity. We did this in two ways. First, we used a stratified all subsets approach (Table 3). The best models contained PD along with SW, PD and PA –while the multivariate functional diversity metrics were not selected as variables in the best models. Secondly, we took an informed approach by building various models from the top five predictors from Table 2. The best model of 26 potential models (AW≈1) is one that includes PD, SW and Nfix (Table S1).

[Figure omitted. See PDF.]

Table 3. Comparison of multivariate predictor models from the stratified all subsets routine.

https://doi.org/10.1371/journal.pone.0005695.t003

We also ran the model explicitly testing Westoby's (1998) LHS plant strategy scheme by regressing productivity against the additive effects of SLA, SW and H. While this model was significant (P<0.0001), only seed weight was significant as an individual model term (P<0.0001) while SLA and H were not (P = 0.261 and 0.190, respectively). Further, this model did not have much of an explanatory advantage over the model with seed weight alone (AIC = 1755.5 vs. 1756.7; and R2 = 0.41 vs. 0.40).

For each of the six continuous traits, we assessed whether traits covaried with relatedness using phylogenetically independent contrasts. SW had a significant phylogenetic signal (summed observed contrasts = 0.131, null expectation [95% CI] = 0.185 [0.136–0.253]), while the rest of the traits lacked any significant signal (LA: obs = 192.5, null [95% CI] = 190.4 [149.5–235.2]; LPA: obs = 306.4, null [95% CI] = 301.8 [229.3–387.9]; LL: obs = 483.1, null [95% CI] = 293.5 [197.4–489.2]; SLA: obs = 1406.4, null [95% CI] = 1374.9 [1121.8–1693.4]; and H: obs = 14.8, null [95% CI] = 14.8 [11.9–18.1]). The significant signal in seed weight is influenced by the fact that species in the Fabaceae tended to have larger seed weights (Fig. 5). We confirmed that the non-Fabaceae species lacked a significant signal by removing Fabaceae species and re-running the analysis (obs = 0.042, null [95% CI] = 0.053 [0.037–0.067]), and species within the Fabaceae also lacked a signal (obs = 0.087, null [95% CI] = 0.104 [0.066–0.135]).

[Figure omitted. See PDF.]

Figure 5. The distribution of seed size across the phylogeny.

https://doi.org/10.1371/journal.pone.0005695.g005

Discussion

Phylogenetic diversity (PD) within a plot was an important factor explaining community productivity patterns. This result is not surprising; we fully expected PD to be a significant predictor of productivity, given that previous studies have shown this [e.g.], [9,18]. On the other hand, we expected multivariate functional diversity indices to better account for productivity patterns, since trait differences should drive ecological differences –regardless of patterns of shared ancestry on the traits. Even though multivariate functional diversity metrics were significant predictors of productivity patterns [see also 25], PD was by far a better predictor. Further, individual traits (besides nitrogen fixation and seed weight) did not show detectable phylogenetic signals, meaning that PD's influence on productivity is likely driven by other, unmeasured traits. This lack of signal could be due to the limited size and coverage of our phylogeny. Further, the seed weight signal appears driven by the fact that the species in the Fabaceae clade tend to have larger seeds than other species, meaning the seed weight results are at least partially driven by the covariance between seed size and being a nitrogen fixer. Although the alternative interpretation would be that the seed weight result likely represents the effect of some other traits involved in a tradeoff with seed size.

Beyond PD, the presence of a nitrogen fixer, was an important explanatory variable. Given that grassland communities at Cedar Creek are generally nitrogen-limited, it is not surprising that the presence of a nitrogen fixer was the single most important factor explaining patterns of community productivity [7]. Yet most polycultures contained at least one nitrogen fixer, meaning that the presence of nitrogen fixers does not inform our understanding of productivity patterns among the most speciose plots. Further, controlling for the presence of N-fixers allows us to overcome a potential phylogenetic “selection” effect. Thus the analyses examining the explanatory value of PD and the various trait measures in plots containing an N-fixer are particularly illuminating. Here again PD comes out as the best single measure, but we cannot eliminate the role of selecting other highly productive species or clades.

Limitations and future directions

While the results presented here strongly support using PD to understand productivity patterns, there are at least two caveats that should be the focus of future research. First, our conclusion that PD is the best predictor of biomass production stems from a single study system at one spatial scale, and PD may not prove to be such a good predictor in other systems or at other scales. In a metaanalysis of PD effects on productivity across 29 experiments, Cadotte et al. (2008) found much variation in the explanatory value of PD (R2 ranged from 0.01 to 0.69). Further, in some systems or under some environmental conditions the critical functional traits may have much more important roles and thus more explanatory value than PD. Such critical traits may be very labile with a low phylogenetic signal, which a phylogeny fails to capture, or else these traits exhibit considerable phenotypic plasticity. Future analyses should include belowground traits, which may be critical in this system for understanding mechanisms of niche partitioning that enhance productivity. A mechanistic understanding of the productivity-diversity relationship may provide greater insight and is often more relevant for practical, management objectives than finding relationships that lack a clear mechanistic basis, as with PD.

Secondly, future research is also required to link PD to ecological differentiation. In this paper, we assume that change in PD is proportional to change in niche space or ecological function, meaning that the magnitude of variation in niche differentiation is dependent on the time since divergence (i.e., Brownian trait evolution). Under such a model, taxa covary in their niche trait based on the amount of time they are represented by a single ancestor in a phylogeny, then differentiate in that niche trait, according to a Markov process of probabilistic trait change diverging from some mean ancestral trait [27]. Other models relating phylogenetic distances to niche or phenotypic divergence could be compared, including those that do not assume that trait divergence is proportional to genetic change [47]. Further, alternative models have been proposed that incorporate explicit models of stabilizing or strong directional selection [e.g.], [48], [49,50]. Addressing alternative models of how phylogenetic distances correspond to trait changes would offer insight into what extent niches are conserved [51] and what this means for ecosystem function. Thus, future work should center on finding probable models that relate phylogenetic differentiation with actual niche and functional differences.

Given these caveats, it is still remarkable that PD is such a strong predictor of productivity patterns. PD might be further preferable because it is becoming increasingly simpler to measure PD compared with functional diversity metrics. Online databases, such as the GenBank (www.ncbi.nlm.nih.gov), the genetic sequence repository, or TreeBASE (www.treebase.org) that stores phylogenetic data from publications, or Phylomatic (www.phylodiversity.net/phylomatic/), which constructs phylogenies for species lists from available Angiosperm supertrees, provide straightforward paths to developing phylogenies. While trait repositories are being developed (e.g., TraitNet; www.columbia.edu/cu/traitnet/) there is no widely available, comprehensive source for multiple traits, and it is unlikely it would cover all the physical and biochemical traits that matter in different systems.

Relating experimental results to natural communities

While phylogenetic diversity may be useful in understanding the consequences of diversity change in ecosystems not amenable to experimental manipulation (e.g., the open ocean, communities with endangered species, communities of long-lived organisms, etc.), future research needs to relate our results to mechanisms of coexistence and productivity in natural communities [10]. In a recent study by Flombaum and Sala [52], primary productivity in natural ecosystems showed a greater response to gradients in species number compared to those from artificially assembled experimental plots, such as with the data used here. Is the stronger diversity-productivity relationship in naturally assembled communities driven more by critical functional traits, or would PD still better explain productivity? Additionally, how communities are assembled could have important ramifications for patterns of PD. Natural community assembly does not select species randomly from a regional species pool, meaning that phylogenetic relationships or ecological traits may inform basic coexistence patterns. Phylogenetic relatedness likely depends on the spatial and taxonomic scale being considered [41], where at broad scales, phylogenetic clustering corresponds to trait selection, while at narrow scales, very close relatives are less likely to co-occur [41], [53]. Certain ecological traits should be selected for in the process of community assembly, and phylogenetic patterns likely reflect the phylogenetic signal of selected traits [54], [55]. Depending on the exact nature of trait evolution and community assembly, communities could show greater phylogenetic diversity with low trait variation and vice versa [55]. The relative contribution of trait variation vs. PD on ecosystem productivity may then also vary with community assembly processes; for instance if community assembly favors high PD and lower critical trait variance [55], then PD may still represent other ecological traits, and thus be significantly related to productivity. Phylogenetic information may be a way to scale from organismal physiology to ecosystem-level processes [56], [57].

Conclusions and implications

We explicitly tested the explanatory power of PD on community productivity against a suite of individual and multivariate traits for experimentally created grassland communities. The trait metrics did not explain productivity patterns as well as PD and one major reason appears to be that we did not measure those functional traits, such as root types, rooting depth or resource requirements, which the phylogeny is representing. We assume that since PD is such a good explanatory variable, it provides a measure of diversity that must be associated with critical functional differences among species within a community that contribute to maximization of productivity.

Even without a clear mechanistic basis, these results have important implications for habitat restoration and biofuel production. If one of the goals of habitat restoration activities is to maximize community functionality, then restoration biologists should attempt to maximize the evolutionary diversity among members of the planned community. For biofuel production, especially in Midwestern prairies, in addition to choosing species with fast growth rates, biofuel-producing communities should consist of phylogenetically distinct mixtures while paying attention to critical functional groups such as nitrogen fixers.

Supporting Information

[Figure omitted. See PDF.]

Appendix S1.

Table of species, genes and GenBank accession numbers and the phylogeny.

https://doi.org/10.1371/journal.pone.0005695.s001

(0.81 MB DOC)

Appendix S2.

R scripts to calculate PD and various measures of trait diversity.

https://doi.org/10.1371/journal.pone.0005695.s002

(0.04 MB DOC)

Table S1.

Multiple model comparison.

https://doi.org/10.1371/journal.pone.0005695.s003

(0.06 MB DOC)

Acknowledgments

MWC wishes to thank Elsa Cleland and Jonathan Davies for inspiring discussions and for advice on the analyses presented here. We also thank participants of the NCEAS Ecophylogenetics working group for valuable discussions. Brian Anacker, Jason Fridley, Andy Hector and Michael Stastny kindly provided comments on an earlier manuscript. JCB wishes to thank Paul Belinski and Kari Kohler for assistance in collecting the trait data.

Author Contributions

Conceived and designed the experiments: MC JCB DT TO. Performed the experiments: JCB DT. Analyzed the data: MC TO. Contributed reagents/materials/analysis tools: MC TO. Wrote the paper: MC JCB DT TO.

Citation: Cadotte MW, Cavender-Bares J, Tilman D, Oakley TH (2009) Using Phylogenetic, Functional and Trait Diversity to Understand Patterns of Plant Community Productivity. PLoS ONE4(5): e5695. https://doi.org/10.1371/journal.pone.0005695

References

1. Tilman D, Reich PB, Knops J, Wedin D, Mielke T, et al. (2001) Diversity and productivity in a long-term grassland experiment. Science 294: 843–845.D. TilmanPB ReichJ. KnopsD. WedinT. Mielke2001Diversity and productivity in a long-term grassland experiment.Science294843845

2. Tilman D, Wedin D, Knops J (1996) Productivity and sustainability influenced by biodiversity in grassland ecosystems. Nature 379: 718–720.D. TilmanD. WedinJ. Knops1996Productivity and sustainability influenced by biodiversity in grassland ecosystems.Nature379718720

3. Naeem S, Thompson LJ, Lawler SP, Lawton JH, Woodfin RM (1994) Declining biodiversity can alter the performance of ecosystems. Nature 368: 734–737.S. NaeemLJ ThompsonSP LawlerJH LawtonRM Woodfin1994Declining biodiversity can alter the performance of ecosystems.Nature368734737

4. Hector A, Schmid B, Beierkuhnlein C, Caldeira MC, Diemer M, et al. (1999) Plant diversity and productivity experiments in European grasslands. Science 286: 1123–1127.A. HectorB. SchmidC. BeierkuhnleinMC CaldeiraM. Diemer1999Plant diversity and productivity experiments in European grasslands.Science28611231127

5. Cardinale BJ, Srivastava DS, Duffy JE, Wright JP, Downing AL, et al. (2006) Effects of biodiversity on the functioning of trophic groups and ecosystems. Nature 443: 989–992.BJ CardinaleDS SrivastavaJE DuffyJP WrightAL Downing2006Effects of biodiversity on the functioning of trophic groups and ecosystems.Nature443989992

6. Cardinale BJ, Wright JP, Cadotte MW, Carroll IT, Hector A, et al. (2007) Impacts of plant diversity on biomass production increase through time due to complimentary resource use: A meta-analysis of 45 experiments. Proceedings of the National Academy of Science, USA 104: 18123–18128.BJ CardinaleJP WrightMW CadotteIT CarrollA. Hector2007Impacts of plant diversity on biomass production increase through time due to complimentary resource use: A meta-analysis of 45 experiments.Proceedings of the National Academy of Science, USA1041812318128

7. Fornara DA, Tilman D (2008) Plant functional composition influences rates of soil carbon and nitrogen accumulation. Journal of Ecology 96: 314–322.DA FornaraD. Tilman2008Plant functional composition influences rates of soil carbon and nitrogen accumulation.Journal of Ecology96314322

8. Harpole WS, Tilman D (2007) Grassland species loss resulting from reduced niche dimension. Nature 446: 791–793.WS HarpoleD. Tilman2007Grassland species loss resulting from reduced niche dimension.Nature446791793

9. Maherali H, Klironomos JN (2007) Influence of phylogeny on fungal community assembly and ecosystem functioning. Science 316: 1746–1748.H. MaheraliJN Klironomos2007Influence of phylogeny on fungal community assembly and ecosystem functioning.Science31617461748

10. Tilman D (1999) The ecological consequences of changes in biodiversity: A search for general principles. Ecology 80: 1455–1474.D. Tilman1999The ecological consequences of changes in biodiversity: A search for general principles.Ecology8014551474

11. Heemsbergen DA, Berg MP, Loreau M, van Hal JR, Faber JH, et al. (2004) Biodiversity Effects on Soil Processes Explained by Interspecific Functional Dissimilarity. Science 306: 1019–1020.DA HeemsbergenMP BergM. LoreauJR van HalJH Faber2004Biodiversity Effects on Soil Processes Explained by Interspecific Functional Dissimilarity.Science30610191020

12. Reich PB, Tilman D, Naeem S, Ellsworth DS, Knops J, et al. (2004) Species and functional group diversity independently influence biomass accumulation and its response to CO2 and N. Proceedings of the National Academy of Sciences of the United States of America 101: 10101–10106.PB ReichD. TilmanS. NaeemDS EllsworthJ. Knops2004Species and functional group diversity independently influence biomass accumulation and its response to CO2 and N.Proceedings of the National Academy of Sciences of the United States of America1011010110106

13. Tilman D, Knops J, Wedin D, Reich P, Ritchie M, et al. (1997) The influence of functional diversity and composition on ecosystem processes. Science 277: 1300–1302.D. TilmanJ. KnopsD. WedinP. ReichM. Ritchie1997The influence of functional diversity and composition on ecosystem processes.Science27713001302

14. Lanta V, Leps J (2006) Effect of functional group richness and species richness in manipulated productivity-diversity studies: a glasshouse pot experiment. Acta Oecologica-International Journal of Ecology 29: 85–96.V. LantaJ. Leps2006Effect of functional group richness and species richness in manipulated productivity-diversity studies: a glasshouse pot experiment.Acta Oecologica-International Journal of Ecology298596

15. Mokany K, Ash J, Roxburgh S (2008) Functional identity is more important than diversity in influencing ecosystem processes in a temperate native grassland. Journal of Ecology 96: 884–893.K. MokanyJ. AshS. Roxburgh2008Functional identity is more important than diversity in influencing ecosystem processes in a temperate native grassland.Journal of Ecology96884893

16. Suding KN, Miller AE, Bechtold H, Bowman WD (2006) The consequence of species loss on ecosystem nitrogen cycling depends on community compensation. Oecologia 149: 141–149.KN SudingAE MillerH. BechtoldWD Bowman2006The consequence of species loss on ecosystem nitrogen cycling depends on community compensation.Oecologia149141149

17. Wohl DL, Arora S, Gladstone JR (2004) Functional redundancy supports biodiversity and ecosystem function in a closed and constant environment. Ecology 85: 1534–1540.DL WohlS. AroraJR Gladstone2004Functional redundancy supports biodiversity and ecosystem function in a closed and constant environment.Ecology8515341540

18. Cadotte MW, Cardinale BJ, Oakley TH (2008) Evolutionary history and the effect of biodiversity on plant productivity. Proceedings of the National Academy of Science, USA 105: 17012–17017.MW CadotteBJ CardinaleTH Oakley2008Evolutionary history and the effect of biodiversity on plant productivity.Proceedings of the National Academy of Science, USA1051701217017

19. Wright JP, Naeem S, Hector A, Lehman C, Reich PB, et al. (2006) Conventional functional classification schemes underestimate the relationship with ecosystem functioning. Ecology Letters 9: 111–120.JP WrightS. NaeemA. HectorC. LehmanPB Reich2006Conventional functional classification schemes underestimate the relationship with ecosystem functioning.Ecology Letters9111120

20. Westoby M (1998) A leaf-height-seed (LHS) plant ecology strategy scheme. Plant and Soil 199: 213–227.M. Westoby1998A leaf-height-seed (LHS) plant ecology strategy scheme.Plant and Soil199213227

21. Grime JP (1974) Vegetation Classification by Reference to Strategies. Nature 250: 26–31.JP Grime1974Vegetation Classification by Reference to Strategies.Nature2502631

22. Wright IJ, Reich PB, Westoby M, Ackerly DD, Baruch Z, et al. (2004) The worldwide leaf economics spectrum. Nature 428: 821–827.IJ WrightPB ReichM. WestobyDD AckerlyZ. Baruch2004The worldwide leaf economics spectrum.Nature428821827

23. Walker B, Kinzig A, Langridge J (1999) Plant attribute diversity, resilience, and ecosystem function: The nature and significance of dominant and minor species. Ecosystems 2: 95–113.B. WalkerA. KinzigJ. Langridge1999Plant attribute diversity, resilience, and ecosystem function: The nature and significance of dominant and minor species.Ecosystems295113

24. Petchey OL, Gaston KJ (2002) Functional diversity (FD), species richness and community composition. Ecology Letters 5: 402–411.OL PetcheyKJ Gaston2002Functional diversity (FD), species richness and community composition.Ecology Letters5402411

25. Petchey OL, Hector A, Gaston KJ (2004) How do different measures of functional diversity perform? Ecology 85: 847–857.OL PetcheyA. HectorKJ Gaston2004How do different measures of functional diversity perform?Ecology85847857

26. Darwin C (1859) The origin of the species by means of natural selection. London: Murray. C. Darwin1859The origin of the species by means of natural selectionLondonMurray

27. Felsenstein J (1985) Phylogenies and the comparative method. American Naturalist 125: 1–15.J. Felsenstein1985Phylogenies and the comparative method.American Naturalist125115

28. Harvey PH, Pagel M (1991) The comparative method in evolutionary biology. Oxford: Oxford University Press. PH HarveyM. Pagel1991The comparative method in evolutionary biologyOxfordOxford University Press

29. Knops JMH, Tilman D, Haddad NM, Naeem S, Mitchell CS, et al. (1999) Effects of plant species richness on invasion dynamics, disease outbreaks, insect abundances and diversity. Ecology Letters 2: 286–293.JMH KnopsD. TilmanNM HaddadS. NaeemCS Mitchell1999Effects of plant species richness on invasion dynamics, disease outbreaks, insect abundances and diversity.Ecology Letters2286293

30. Haddad NM, Tilman D, Knops JMH (2002) Long-term oscillations in grassland productivity induced by drought. Ecology Letters 5: 110–120.NM HaddadD. TilmanJMH Knops2002Long-term oscillations in grassland productivity induced by drought.Ecology Letters5110120

31. Tilman D, Reich PB, Knops JMH (2006) Biodiversity and ecosystem stability in a decade-long grassland experiment. Nature 441: 629–632.D. TilmanPB ReichJMH Knops2006Biodiversity and ecosystem stability in a decade-long grassland experiment.Nature441629632

32. Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL (2005) GenBank. Nucleic Acids Research 33: D34–D38.DA BensonI. Karsch-MizrachiDJ LipmanJ. OstellDL Wheeler2005GenBank.Nucleic Acids Research33D34D38

33. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research 32: 1792–1797.RC Edgar2004MUSCLE: multiple sequence alignment with high accuracy and high throughput.Nucleic Acids Research3217921797

34. Posada D, Crandall KA (1998) MODELTEST: testing the model of DNA substitution. Bioinformatics 14: 817–818.D. PosadaKA Crandall1998MODELTEST: testing the model of DNA substitution.Bioinformatics14817818

35. Posada D, Crandall KA (2001) Selecting the Best-Fit Model of Nucleotide Substitution. Systematic Biology 50: 580–601.D. PosadaKA Crandall2001Selecting the Best-Fit Model of Nucleotide Substitution.Systematic Biology50580601

36. Anisimova M, Gascuel O (2006) Approximate likelihood ratio test for branches: A fast accurate and powerful alternative. Systematic Biology 55: 539–552.M. AnisimovaO. Gascuel2006Approximate likelihood ratio test for branches: A fast accurate and powerful alternative.Systematic Biology55539552

37. Guindon S, Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic Biology 52: 696–704.S. GuindonO. Gascuel2003A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood.Systematic Biology52696704

38. Cadotte MW, Hamilton MA, Murray BR (2009) Phylogenetic relatedness and plant invader success across two spatial scales. Diversity and Distributions 15: 481–488.MW CadotteMA HamiltonBR Murray2009Phylogenetic relatedness and plant invader success across two spatial scales.Diversity and Distributions15481488

39. Rasband WS (1997–2004) ImageJ. WS Rasband1997–2004ImageJ.http://rsb.info.nih.gov/ij/ ed. Bethesda, Maryland: National Institutes of Health. http://rsb.info.nih.gov/ij/ ed. Bethesda, Maryland: National Institutes of Health.

40. Sack L, Cowan PD, Jaikumar N, Holbrook NM (2003) The ‘hydrology’ of leaves: co-ordination of structure and function in temperate woody species. Plant Cell and Environment 26: 1343–1356.L. SackPD CowanN. JaikumarNM Holbrook2003The ‘hydrology’ of leaves: co-ordination of structure and function in temperate woody species.Plant Cell and Environment2613431356

41. Cavender-Bares J, Keen A, Miles B (2006) Phylogenetic structure of floridian plant communities depends on taxonomic and spatial scale. Ecology 87: S109–S122.J. Cavender-BaresA. KeenB. Miles2006Phylogenetic structure of floridian plant communities depends on taxonomic and spatial scale.Ecology87S109S122

42. Givnish T (1976) Sizes and shapes of liane leaves. American Naturalist 110: 743–778.T. Givnish1976Sizes and shapes of liane leaves.American Naturalist110743778

43. Legendre P, Legendre L (1998) Numerical Ecology. Amsterdam: Elsevier Science B. V. P. LegendreL. Legendre1998Numerical EcologyAmsterdamElsevier Science B. V

44. Kent M, Coker P (1992) Vegetation description and analysis. New York: John Wiley & Sons. M. KentP. Coker1992Vegetation description and analysisNew YorkJohn Wiley & Sons

45. Johnson JB, Omland KS (2004) Model selection in ecology and evolution. Trends in Ecology & Evolution 19: 101–108.JB JohnsonKS Omland2004Model selection in ecology and evolution.Trends in Ecology & Evolution19101108

46. Davies TJ, Fritz SA, Grenyer R, Orme CDL, Bielby J, et al. (2008) Phylogenetic trees and the future of mammalian biodiversity. Proceedings of the National Academy of Sciences of the United States of America 105: 11556–11563.TJ DaviesSA FritzR. GrenyerCDL OrmeJ. Bielby2008Phylogenetic trees and the future of mammalian biodiversity.Proceedings of the National Academy of Sciences of the United States of America1051155611563

47. Oakley TH, Gu ZL, Abouheif E, Patel NH, Li WH (2005) Comparative methods for the analysis of gene expression evolution: An example using yeast functional genomic data. Molecular Biology and Evolution 22: 40–50.TH OakleyZL GuE. AbouheifNH PatelWH Li2005Comparative methods for the analysis of gene expression evolution: An example using yeast functional genomic data.Molecular Biology and Evolution224050

48. Butler MA, King AA (2004) Phylogenetic comparative analysis: A modeling approach for adaptive evolution. American Naturalist 164: 683–695.MA ButlerAA King2004Phylogenetic comparative analysis: A modeling approach for adaptive evolution.American Naturalist164683695

49. Hansen TF (1997) Stabilizing selection and the comparative analysis of adaptation. Evolution 51: 1341–1351.TF Hansen1997Stabilizing selection and the comparative analysis of adaptation.Evolution5113411351

50. Hansen TF, Martins EP (1996) Translating between microevolutionary process and macroevolutionary patterns: The correlation structure of interspecific data. Evolution 50: 1404–1417.TF HansenEP Martins1996Translating between microevolutionary process and macroevolutionary patterns: The correlation structure of interspecific data.Evolution5014041417

51. Wiens JJ, Graham CH (2005) Niche conservatism: Integrating evolution, ecology, and conservation biology. Annual Review of Ecology Evolution and Systematics 36: 519–539.JJ WiensCH Graham2005Niche conservatism: Integrating evolution, ecology, and conservation biology.Annual Review of Ecology Evolution and Systematics36519539

52. Flombaum P, Sala OE (2008) From the Cover: Higher effect of plant species diversity on productivity in natural than artificial ecosystems. Proceedings of the National Academy of Sciences 105: 6087–6090.P. FlombaumOE Sala2008From the Cover: Higher effect of plant species diversity on productivity in natural than artificial ecosystems.Proceedings of the National Academy of Sciences10560876090

53. Webb CO, Ackerly DD, McPeek MA, Donoghue MJ (2002) Phylogenies and community ecology. Annual Review of Ecology and Systematics 33: 475–505.CO WebbDD AckerlyMA McPeekMJ Donoghue2002Phylogenies and community ecology.Annual Review of Ecology and Systematics33475505

54. Cavender-Bares J, Wilczek A (2003) Integrating micro- and macroevolutionary processes in community ecology. Ecology 84: 592–597.J. Cavender-BaresA. Wilczek2003Integrating micro- and macroevolutionary processes in community ecology.Ecology84592597

55. Prinzing A, Reiffers R, Braakhekke WG, Hennekens SM, Tackenberg O, et al. (2008) Less lineages - more trait variation: phylogenetically clustered plant communities are functionally more diverse. Ecology Letters 11: 809–819.A. PrinzingR. ReiffersWG BraakhekkeSM HennekensO. Tackenberg2008Less lineages - more trait variation: phylogenetically clustered plant communities are functionally more diverse.Ecology Letters11809819

56. Edwards EJ, Still CJ, Donoghue MJ (2007) The relevance of phylogeny to studies of global change. Trends in Ecology & Evolution 22: 243–249.EJ EdwardsCJ StillMJ Donoghue2007The relevance of phylogeny to studies of global change.Trends in Ecology & Evolution22243249

57. Cavender-Bares J, Kozak K, Fine PVA, Kembel SW (2009) The merging of community ecology and phylogenetic biology. Ecology Letters. J. Cavender-BaresK. KozakPVA FineSW Kembel2009The merging of community ecology and phylogenetic biology.Ecology LettersIn press. In press.