Taxon sampling effects on the quantification and comparison of community phylogenetic diversity

Assessing digital accessible botanical knowledge and priorities for exploration and discovery of plant diversity across Mesoamerica

Digital accessible biodiversity knowledge has the potential to greatly advance botanical research and guide conservation efforts. Evaluating its shortfalls is key to understanding its limits and prioritising regions in need of renewed survey efforts. We used the Royal Botanical Gardens Kew's World Checklist of Vascular Plants to parse publicly available occurrence data downloaded from the Global Biodiversity Information Facility and quantify the spatial distribution of spatial, phylogenetic, and temporal data shortfalls across Mesoamerica. After processing 3578 777 occurrence records for 32 522 species of vascular plants across Mesoamerica, we found evidence of poor data coverage: incomplete characterisation of species diversity, old occurrence records, and low phylogenetic representation. One-third of the region showed large gaps for at least one of these dimensions (hotspots) and < 15% had adequate data coverage across dimensions. Overall, the shortfalls we identified compromise the quality of digitally available occurrence data and hamper research on spatial phylogenetics and species dynamics under anthropogenic disturbances. Our analyses identified areas of opportunity for increased efforts in data digitisation, botanical exploration, sequencing, and biodiversity monitoring. These efforts would serve to increase and rejuvenate knowledge on the geographic distribution of vascular plants in Mesoamerica.

Phylogeny-based assignment of functional traits to DNA barcodes outperforms distance-based, in a comparison of approaches

The full potential for using DNA barcodes for profiling functional trait diversity has yet to be determined in plants and animals; thus, we outline a general framework for quantifying functional trait diversity of insect community DNA and propose and assess the accuracy of three methods for achieving this. We built a novel dataset of traits and DNA barcodes for wild bees in China. An informatics framework was developed for phylogeny-based integration of these data and prediction of traits for any subject barcodes, which was compared with two distance-based methods. For Phylogenetic Assignment, we additionally conducted a species-level analysis of publically available bee trait data. Under the specimen-level dataset, the rate of trait assignment was negatively correlated with distance between the query and the nearest trait-known reference, for all methods. Phylogenetic Assignment was found to perform best under several criteria; particularly, it had the lowest false-positive rate (rarely returning a state prediction where success was unlikely; where the distance from query to the nearest reference was high). For a wider range of compiled traits, conservative life-history traits showed the highest rates of assignment; for example, sociality was predicted with confidence at 53%, parasitism at 44% and nest location at 33%. As outlined herein, automated trait assignment might be applied at scale to either barcodes or metabarcodes. With further compilation and databasing of DNA barcode and trait data, the rate and accuracy of trait assignment is expected to increase to the point of being a widely viable and informative approach.

Community‐level phylogenetic diversity does not differ between rare and common lineages across tallgrass prairies in the northern Great Plains

Niche differentiation has served as one explanation for species coexistence, and phylogenetic relatedness provides a means to approximate how ecologically similar species are to each other. To explore the contribution of rare species to community phylogenetic diversity, we sampled 21 plant communities across the Prairie Coteau ecoregion, an area of high conservation concern. We used breakpoint analysis through the iterative addition of less abundant species to the phylogenetic tree for each community to assess the contribution of rare species to community phylogenetic diversity. We also quantify the phylogenetic signal of abundance using Blomberg's K statistic and calculated the phylogenetic similarity between rare and common species using a phylogenetic beta‐diversity metric (D_nn). To estimate the phylogenetic structuring of these prairie communities, we calculated two common metrics that capture evolutionary relatedness at different scales (MPD and MNTD). Additionally, we examine the correlation between Faith's PD, MPD, and MNTD and species richness. We found rare species do not generally contribute higher levels of phylogenetic diversity than common species. Eight communities had significant breakpoints, with only four communities having an increasing trend for the rarest species. The phylogenetic signal for abundance was low but significant in only four communities, and communities had lower phylogenetic diversity than expected from the regional species pool. Finally, the strength of the correlation between species richness and phylogenetic diversity was mixed. Our results indicate niche differentiation does not explain the persistence of rare species in tallgrass prairies, as they were more closely related than expected from random, suggesting high functional redundancy between rare and common species. This is promising for the long‐term resilience of this ecosystem, but only insofar as enough species remain in the system. With ongoing biodiversity loss, it is essential that we understand the role rare species play in their communities.

Leaf out time correlates with wood anatomy across large geographic scales and within local communities

There is a long-standing idea that the timing of leaf production in seasonally cold climates is linked to xylem anatomy, specifically vessel diameter because of the hydraulic requirements of expanding leaves. We tested for a relationship between the timing of leaf out and vessel diameter in 220 plants in three common gardens accounting for species' phylogenetic relationships. We investigated how vessel diameter related to wood porosity, plant height and leaf length. We also used dye perfusion tests to determine whether plants relied on xylem produced during the previous growing season at the time of leaf out. In all three gardens, there was later leaf out in species with wider vessels. Ring-porous species had the widest vessels, exhibited latest leaf out and relied less on xylem made during the previous growing season than diffuse-porous species. Wood anatomy and leaf phenology did not exhibit a phylogenetic signal. The timing of leaf out is correlated with wood anatomy across species regardless of species' geographic origin and phylogenetic relationships. This correlation could be a result of developmental and physiological links between leaves and wood or tied to a larger safety efficiency trade-off.

Spatial phylogenetics of Japanese ferns: Patterns, processes, and implications for conservation

PREMISE

Biodiversity is often only measured with species richness; however, this metric ignores evolutionary history and is not sufficient for making conservation decisions. Here, we characterize multiple facets and drivers of biodiversity to understand how these relate to bioregions and conservation status in the ferns of Japan.

METHODS

We compiled a community data set of 1239 grid cells (20 × 20 km each) including 672 taxa based on >300,000 specimen records. We combined the community data with a phylogeny and functional traits to analyze taxonomic, phylogenetic, and functional diversity and modeled biodiversity metrics in response to environmental factors and reproductive mode. Hierarchical clustering was used to delimit bioregions. Conservation status and threats were assessed by comparing the overlap of significantly diverse grid cells with conservation zones and range maps of native Japanese deer.

RESULTS

Taxonomic richness was highest at mid-latitudes. Phylogenetic and functional diversity and phylogenetic endemism were highest in small southern islands. Relative phylogenetic and functional diversity were high at high and low latitudes, and low at mid-latitudes. Grid cells were grouped into three (phylogenetic) or four (taxonomic) major bioregions. Temperature and apomixis were identified as drivers of biodiversity patterns. Conservation status was generally high for grid cells with significantly high biodiversity, but the threat due to herbivory by deer was greater for taxonomic richness than other metrics.

CONCLUSIONS

Our integrative approach reveals previously undetected patterns and drivers of biodiversity in the ferns of Japan. Future conservation efforts should recognize that threats can vary by biodiversity metric and consider multiple metrics when establishing conservation priorities.

Phylogenetic signal of sub-arctic beetle communities

Postglacial dispersal and colonization processes have shaped community patterns in sub-Arctic regions such as Churchill, Manitoba, and Canada. This study investigates evolutionary community structure within the beetle (Coleoptera) families of Churchill and tests whether biological traits have played a role in governing colonization patterns from refugial and southerly geographic regions. This study quantifies sub-Arctic beetle phylogenetic community structure for each family using the net relatedness index (NRI) and nearest taxon index (NTI), calculated using publicly available data from the Barcode of Life Data Systems (BOLD); compares patterns across families with different traits (habitat, diet) using standard statistical analysis (ANOVA) as well as phylogenetic generalized least squares (PGLS) using a family-level beetle phylogeny obtained from the literature; and compares community structure in Churchill with a region in southern Canada (Guelph, Ontario). These analyses were also repeated at a genus level. The dominant pattern detected in our study was that aquatic families were much better represented in Churchill compared to terrestrial families, when compared against richness sampled from across Canada and Alaska. Individually, most families showed significant phylogenetic clustering in Churchill, likely due to the strong environmental filtering present in Arctic environments. There was no significant difference in phylogenetic structure between Churchill and Guelph but with a trend toward stronger clustering in the North. Fungivores were significantly more overdispersed than other feeding modes, predators were significantly more clustered, and aquatic families showed significantly stronger clustering compared to terrestrial. This study contributes to our understanding of the traits and processes structuring insect biodiversity and macroecological trends in the sub-Arctic.

Testing Phylogenetic Stability with Variable Taxon Sampling

Over the past three decades, computational capabilities have grown at such a rapid rate that they have given rise to many computationally heavy science fields such as phylogenomics. As increasingly more genomes are sequenced in the three domains of life, larger and more species-complete phylogenetic tree reconstructions are leading to a better understanding of the tree of life and the evolutionary histories in deep times. However, these large datasets pose unique challenges from a modeling and computational perspective: accurately describing the evolutionary process of thousands of species is still beyond the capability of current models, while the computational burden limits our ability to test multiple hypotheses. Thus, it is common practice to reduce the size of a dataset by selecting species to represent a clade (taxon sampling). Unfortunately, this process is subjective, and comparisons of large tree of life studies show that choice and number of species used in a dataset can alter the topology obtained. Thus, taxon sampling is, in itself, a process that needs to be fully investigated to determine its effect on phylogenetic stability. Here, we present the theory and practical application of an automated pipeline that can be easily implemented to explore the effect of taxon sampling on phylogenetic reconstructions. The application of this approach was recently discussed in a study of Terrabacteria and shows its power in investigating the accuracy of deep nodes of a phylogeny.

OpenTree: A Python Package for Accessing and Analyzing Data from the Open Tree of Life

The Open Tree of Life project constructs a comprehensive, dynamic, and digitally available tree of life by synthesizing published phylogenetic trees along with taxonomic data. Open Tree of Life provides web-service application programming interfaces (APIs) to make the tree estimate, unified taxonomy, and input phylogenetic data available to anyone. Here, we describe the Python package opentree, which provides a user friendly Python wrapper for these APIs and a set of scripts and tutorials for straightforward downstream data analyses. We demonstrate the utility of these tools by generating an estimate of the phylogenetic relationships of all bird families, and by capturing a phylogenetic estimate for all taxa observed at the University of California Merced Vernal Pools and Grassland Reserve.[Evolution; open science; phylogenetics; Python; taxonomy.].

Phylogenetic uncertainty and the inference of patterns in community ecology and comparative studies

Progress in phylogenetic community ecology is often limited by the availability of phylogenetic information and the lack of appropriate methods and solutions to deal with this problem. We estimate the effect of the lack of phylogenetic information on the relations among taxa measured by commonly used phylogenetic metrics in comparative studies and community ecology, namely: Blomberg's K phylogenetic signal, Faith's Phylogenetic Diversity (PD), Mean Phylogenetic Distance (MPD) and Mean Nearest Taxon Distance (MNTD). To overcome this problem, we tested two possible solutions: Polytomic trees and Operational trees. Our results show that the effects on K values strongly depended on the level of phylogenetic signal. In the case of the community metrics, the effects were insensitive to the patterns of species distribution in the communities. Community metrics tended to be overestimated with both Polytomic and Operational trees, but the overestimation was higher with Polytomic trees. PD and MPD metrics were less biased than MNTD metric. We show that the lack of phylogenetic resolution is not necessarily problematic for all analyses and that its effect will depend on the chosen metric and on the solutions used to deal with the problem. Based on our results, we suggest that ecologists should prefer the Operational tree solution to remove polytomies in the phylogenetic tree and take careful consideration while designing experiments, and analyzing and interpreting the results of phylogenetic metrics.

Building a Robust, Densely-Sampled Spider Tree of Life for Ecosystem Research

Phylogenetic relatedness is a key diversity measure for the analysis and understanding of how species and communities evolve across time and space. Understanding the nonrandom loss of species with respect to phylogeny is also essential for better-informed conservation decisions. However, several factors are known to influence phylogenetic reconstruction and, ultimately, phylogenetic diversity metrics. In this study, we empirically tested how some of these factors (topological constraint, taxon sampling, genetic markers and calibration) affect phylogenetic resolution and uncertainty. We built a densely sampled, species-level phylogenetic tree for spiders, combining Sanger sequencing of species from local communities of two biogeographical regions (Iberian Peninsula and Macaronesia) with a taxon-rich backbone matrix of Genbank sequences and a topological constraint derived from recent phylogenomic studies. The resulting tree constitutes the most complete spider phylogeny to date, both in terms of terminals and background information, and may serve as a standard reference for the analysis of phylogenetic diversity patterns at the community level. We then used this tree to investigate how partial data affect phylogenetic reconstruction, phylogenetic diversity estimates and their rankings, and, ultimately, the ecological processes inferred for each community. We found that the incorporation of a single slowly evolving marker (28S) to the DNA barcode sequences from local communities, had the highest impact on tree topology, closely followed by the use of a backbone matrix. The increase in missing data resulting from combining partial sequences from local communities only had a moderate impact on the resulting trees, similar to the difference observed when using topological constraints. Our study further revealed substantial differences in both the phylogenetic structure and diversity rankings of the analyzed communities estimated from the different phylogenetic treatments, especially when using non-ultrametric trees (phylograms) instead of time-stamped trees (chronograms). Finally, we provide some recommendations on reconstructing phylogenetic trees to infer phylogenetic diversity within ecological studies.

Darwin's naturalization conundrum can be explained by spatial scale

Darwin proposed two seemingly contradictory hypotheses regarding factors influencing the outcome of biological invasions. He initially posited that nonnative species closely related to native species would be more likely to successfully establish, because they might share adaptations to the local environment (preadaptation hypothesis). However, based on observations that the majority of naturalized plant species in the United States belonged to nonnative genera, he concluded that the lack of competitive exclusion would facilitate the establishment of alien invaders phylogenetically distinct from the native flora (competition-relatedness hypothesis). To date, no consensus has been reached regarding these opposing hypotheses. Here, following Darwin, we use the flora of the United States to examine patterns of taxonomic and phylogenetic relatedness between native and nonnative taxa across thousands of nested locations ranging in size and extent, from local to regional scales. We find that the probability of observing the signature of environmental filtering over that of competition increases with spatial scale. Further, native and nonnative species tended to be less related in warm, humid environments. Our work provides an empirical assessment of the role of observation scale and climate in biological invasions and demonstrates that Darwin's two opposing hypotheses need not be mutually exclusive.

Endemism patterns are scale dependent

Areas of endemism are important in biogeography because they capture facets of biodiversity not represented elsewhere. However, the scales at which they are relevant to research and conservation are poorly analysed. Here, we calculate weighted endemism (WE) and phylogenetic endemism (PE) separately for all birds and amphibians across the globe. We show that scale dependence is widespread for both indices and manifests across grain sizes, spatial extents and taxonomic treatments. Variations in taxonomic opinions-whether species are treated by systematic 'lumping' or 'splitting'-can profoundly affect the allocation of WE hotspots. Global patterns of PE can provide insights into complex evolutionary processes but this congruence is lost at the continental to country extents. These findings are explained by environmental heterogeneity at coarser grains, and to a far lesser extent at finer resolutions. Regardless of scale, we find widespread deficits of protection for endemism hotspots. Our study presents a framework for assessing areas for conservation that are robust to assumptions on taxonomy, spatial grain and extent.

Clade composition of a plant community indicates its phylogenetic diversity

Phylogenetic diversity quantification is based on indices computed from phylogenetic distances among species, which are derived from phylogenetic trees. This approach requires phylogenetic expertise and available molecular data, or a fully sampled synthesis-based phylogeny. Here, we propose and evaluate a simpler alternative approach based on taxonomic coding. We developed metrics, the clade indices, based on information about clade proportions in communities and species richness of a community or a clade, which do not require phylogenies. Using vegetation records from herbaceous plots from Central Europe and simulated vegetation plots based on a megaphylogeny of vascular plants, we examined fit accuracy of our proposed indices for all dimensions of phylogenetic diversity (richness, divergence, and regularity). For real vegetation data, the clade indices fitted phylogeny-based metrics very accurately (explanatory power was usually higher than 80% for phylogenetic richness, almost always higher than 90% for phylogenetic divergence, and often higher than 70% for phylogenetic regularity). For phylogenetic regularity, fit accuracy was habitat and species richness dependent. For phylogenetic richness and divergence, the clade indices performed consistently. In simulated datasets, fit accuracy of all clade indices increased with increasing species richness, suggesting better precision in species-rich habitats and at larger spatial scales. Fit accuracy for phylogenetic divergence and regularity was unreliable at large phylogenetic scales, suggesting inadvisability of our method in habitats including many distantly related lineages. The clade indices are promising alternative measures for all projects with a phylogenetic framework, which can trade-off a little precision for a significant speed-up and simplification, such as macroecological analyses or where phylogenetic data is incomplete.

Effects of taxon sampling and tree reconstruction methods on phylodiversity metrics

The amount and patterns of phylodiversity in a community are often used to draw inferences about the local and historical factors affecting community assembly and can be used to prioritize communities and locations for conservation. Because measures of phylodiversity are based on the topology and branch lengths of phylogenetic trees, which are affected by the number and diversity of taxa in the tree, these analyses may be sensitive to changes in taxon sampling and tree reconstruction methods.To investigate the effects of taxon sampling and tree reconstruction methods on measures of phylodiversity, we investigated the community phylogenetics of the Ordway-Swisher Biological Station (Florida), which is home to over 600 species of vascular plants. We studied the effects of (a) the number of taxa included in the regional phylogeny; (b) random versus targeted sampling of species to assemble the regional species pool; (c) including only species from specific clades rather than broad sampling; (d) using trees reconstructed directly for the taxa under study compared to trees pruned from a larger reconstructed tree; and (e) using phylograms compared to chronograms.We found that including more taxa in a study increases the likelihood of observing significantly nonrandom phylogenetic patterns. However, there were no consistent trends in the phylodiversity patterns based on random taxon sampling compared to targeted sampling, or within individual clades compared to the complete dataset. Using pruned and reconstructed phylogenies resulted in similar patterns of phylodiversity, while chronograms in some cases led to significantly different results from phylograms.The methods commonly used in community phylogenetic studies can significantly impact the results, potentially influencing both inferences of community assembly and conservation decisions. We highlight the need for both careful selection of methods in community phylogenetic studies and appropriate interpretation of results, depending on the specific questions to be addressed.

For common community phylogenetic analyses, go ahead and use synthesis phylogenies

Should we build our own phylogenetic trees based on gene sequence data, or can we simply use available synthesis phylogenies? This is a fundamental question that any study involving a phylogenetic framework must face at the beginning of the project. Building a phylogeny from gene sequence data (purpose-built phylogeny) requires more effort, expertise, and cost than subsetting an already available phylogeny (synthesis-based phylogeny). However, we still lack a comparison of how these two approaches to building phylogenetic trees influence common community phylogenetic analyses such as comparing community phylogenetic diversity and estimating trait phylogenetic signal. Here, we generated three purpose-built phylogenies and their corresponding synthesis-based trees (two from Phylomatic and one from the Open Tree of Life, OTL). We simulated 1,000 communities and 12,000 continuous traits along each purpose-built phylogeny. We then compared the effects of different trees on estimates of phylogenetic diversity (alpha and beta) and phylogenetic signal (Pagel's λ and Blomberg's K). Synthesis-based phylogenies generally yielded higher estimates of phylogenetic diversity when compared to purpose-built phylogenies. However, resulting measures of phylogenetic diversity from both types of phylogenies were highly correlated (Spearman's ρ  > 0.8 in most cases). Mean pairwise distance (both alpha and beta) is the index that is most robust to the differences in tree construction that we tested. Measures of phylogenetic diversity based on the OTL showed the highest correlation with measures based on the purpose-built phylogenies. Trait phylogenetic signal estimated with synthesis-based phylogenies, especially from the OTL, was also highly correlated with estimates of Blomberg's K or close to Pagel's λ from purpose-built phylogenies when traits were simulated under Brownian motion. For commonly employed community phylogenetic analyses, our results justify taking advantage of recently developed and continuously improving synthesis trees, especially the Open Tree of Life.

Increasing phylogenetic stochasticity at high elevations on summits across a remote North American wilderness

PREMISE

At the intersection of ecology and evolutionary biology, community phylogenetics can provide insights into overarching biodiversity patterns, particularly in remote and understudied ecosystems. To understand community assembly of the high alpine flora in the Sawtooth National Forest, USA, we analyzed phylogenetic structure within and between nine summit communities.

METHODS

We used high-throughput sequencing to supplement existing data and infer a nearly completely sampled community phylogeny of the alpine vascular flora. We calculated mean nearest taxon distance (MNTD) and mean pairwise distance (MPD) to quantify phylogenetic divergence within summits, and assessed whether maximum elevation explains phylogenetic structure. To evaluate similarities between summits, we quantified phylogenetic turnover, taking into consideration microhabitats (talus vs. meadows).

RESULTS

We found different patterns of community phylogenetic structure within the six most species-rich orders, but across all vascular plants phylogenetic structure was largely not different from random. There was a significant negative correlation between elevation and tree-wide phylogenetic diversity (MPD) within summits: overdispersion degraded as elevation increased. Between summits, we found high phylogenetic turnover driven by greater niche heterogeneity on summits with alpine meadows.

CONCLUSIONS

Our results provide further evidence that stochastic processes may also play an important role in the assembly of vascular plant communities in high alpine habitats at regional scales. However, order-specific patterns suggest that adaptations are still important for assembly of specific sectors of the plant tree of life. Further studies quantifying functional diversity will be important in disentangling the interplay of eco-evolutionary processes that likely shape broad community phylogenetic patterns in extreme environments.