Combinatorial properties of phylogenetic diversity indices

Journal of Mathematical Biology - Phylogenetic diversity indices provide a formal way to apportion ‘evolutionary heritage’ across species. Two natural diversity indices are Fair...     

A comparative test of phylogenetic diversity indices

Traditional measures of biodiversity, such as species richness, usually treat species as being equal. As this is obviously not the case, measuring diversity in terms of features accumulated over evolutionary history provides additional value to theoretical and applied ecology. Several phylogenetic diversity indices exist, but their behaviour has not yet been tested in a comparative framework. We provide a test of ten commonly used phylogenetic diversity indices based on 40 simulated phylogenies of varying topology. We restrict our analysis to a topological fully resolved tree without information on branch lengths and species lists with presence-absence data. A total of 38,000 artificial communities varying in species richness covering 5-95% of the phylogenies were created by random resampling. The indices were evaluated based on their ability to meet a priori defined requirements. No index meets all requirements, but three indices turned out to be more suitable than others under particular conditions. Average taxonomic distinctness (AvTD) and intensive quadratic entropy (J) are calculated by averaging and are, therefore, unbiased by species richness while reflecting phylogeny per se well. However, averaging leads to the violation of set monotonicity, which requires that species extinction cannot increase the index. Total taxonomic distinctness (TTD) sums up distinctiveness values for particular species across the community. It is therefore strongly linked to species richness and reflects phylogeny per se weakly but satisfies set monotonicity. We suggest that AvTD and J are best applied to studies that compare spatially or temporally rather independent communities that potentially vary strongly in their phylogenetic composition-i.e. where set monotonicity is a more negligible issue, but independence of species richness is desired. In contrast, we suggest that TTD be used in studies that compare rather interdependent communities where changes occur more gradually by species extinction or introduction. Calculating AvTD or TTD, depending on the research question, in addition to species richness is strongly recommended.     

Scale‐dependence of phylogenetic signal in ecological traits of ectoparasites

The non‐independence of traits among closely related species is a well‐documented phenomenon underpinning modern methods for comparative analyses or prediction of trait values in new species. Surprisingly such studies have mainly focused on life‐history or morphological traits of free‐living organisms, ignoring ecological attributes of parasite species in spite of the fact that they are critical for conservation and human health. We tested for a phylogenetic signal acting on two ecological traits, abundance and host specificity, using data for 218 flea species parasitic on small mammals in 19 regions of the Palaearctic and Nearctic, and a phylogenetic tree for these species. We tested for the presence of a phylogenetic signal at both regional and continental scales using three measures (Abouheif/Moran's I, Pagel's λ, and Blomberg et al.'s K). Our results show 1) a consistent positive phylogenetic signal for flea abundance, but only a weaker and erratic signal for host specificity, and 2) a clear dependence on scale, with the signals being stronger at the continental scale and relatively weaker or inconsistent at the regional scale. Whenever values of Blomberg et al.'s K were found significant, they were <1 suggesting that the effects of phylogeny on the evolution of abundance and host specificity in fleas are weaker than expected from a Brownian motion model. The most striking finding is that, within a continental fauna, closely‐related flea species are characterized by similar levels of abundance, though this pattern is weaker within local assemblages, possibly eroded by local biotic or abiotic conditions. We discuss the link between history (represented by phylogeny) and pattern of variation among species in morphological and ecological traits, and use comparisons between the Palaearctic and Nearctic to infer a role of historical events in the probability of detecting phylogenetic signals.     

Influence of tree shape and evolutionary time‐scale on phylogenetic diversity metrics

During the last decades, describing, analysing and understanding the phylogenetic structure of species assemblages has been a central theme in both community ecology and macro‐ecology. Among the wide variety of phylogenetic structure metrics, three have been predominant in the literature: Faith's phylogenetic diversity (PD_Faith), which represents the sum of the branch lengths of the phylogenetic tree linking all species of a particular assemblage, the mean pairwise distance between all species in an assemblage (MPD) and the pairwise distance between the closest relatives in an assemblage (MNTD). Comparisons between studies using one or several of these metrics are difficult because there has been no comprehensive evaluation of the phylogenetic properties each metric captures. In particular it is unknown how PD_Faith relates to MDP and MNTD. Consequently, it is possible that apparently opposing patterns in different studies might simply reflect differences in metric properties. Here, we aim to fill this gap by comparing these metrics using simulations and empirical data. We first used simulation experiments to test the influence of community structure and size on the mismatch between metrics whilst varying the shape and size of the phylogenetic tree of the species pool. Second we investigated the mismatch between metrics for two empirical datasets (gut microbes and global carnivoran assemblages). We show that MNTD and PD_Faith provide similar information on phylogenetic structure, and respond similarly to variation in species richness and assemblage structure. However, MPD demonstrate a very different behaviour, and is highly sensitive to deep branching structure. We suggest that by combining complementary metrics that are sensitive to processes operating at different phylogenetic depths (i.e. MPD and MNTD or PD_Faith) we can obtain a better understanding of assemblage structure.     

A phylogenetic perspective on species diversity, β‐diversity and biogeography for the microbial world

There is an increasing interest to combine phylogenetic data with distributional and ecological records to assess how natural communities arrange under an evolutionary perspective. In the microbial world, there is also a need to go beyond the problematic species definition to deeply explore ecological patterns using genetic data. We explored links between evolution/phylogeny and community ecology using bacterial 16S rRNA gene information from a high‐altitude lakes district data set. We described phylogenetic community composition, spatial distribution, and β‐diversity and biogeographical patterns applying evolutionary relatedness without relying on any particular operational taxonomic unit definition. High‐altitude lakes districts usually contain a large mosaic of highly diverse small water bodies and conform a fine biogeographical model of spatially close but environmentally heterogeneous ecosystems. We sampled 18 lakes in the Pyrenees with a selection criteria focused on capturing the maximum environmental variation within the smallest geographical area. The results showed highly diverse communities nonrandomly distributed with phylogenetic β‐diversity patterns mainly shaped by the environment and not by the spatial distance. Community similarity based on both bacterial taxonomic composition and phylogenetic β‐diversity shared similar patterns and was primarily structured by similar environmental drivers. We observed a positive relationship between lake area and phylogenetic diversity with a slope consistent with highly dispersive planktonic organisms. The phylogenetic approach incorporated patterns of common ancestry into bacterial community analysis and emerged as a very convenient analytical tool for direct inter‐ and intrabiome biodiversity comparisons and sorting out microbial habitats with potential application in conservation studies.     

Breeding systems and phylogenetic diversity of seed plants along a large‐scale elevational gradient

Aim The factors affecting the distribution of dioecious species have received surprisingly little attention and, as a consequence, are poorly understood. For example, there is a well‐documented negative association between dioecy and latitude, for which we have no candidate mechanisms. Conversely, it has been suggested that the relative proportion of dioecious species should be positively correlated with changes in elevation. Location Costa Rica, Central America. Methods We made use of data on the distribution of 175 seed plant species from a series of plots along a transect in Costa Rica that ranged from 30 to 2600 m a.s.l. to test the prediction that dioecy is correlated with elevation. Specifically, we examined relationships between the proportion of dioecy, at the species and individual levels, and the elevation, species richness, number of individuals, and phylogenetic diversity (PD) of plots. For comparison, we repeated all analyses with monoecious species, which also have unisexual flowers but do not suffer from reduced mate assurance and the seed shadow effect that may be the outcomes of having spatially separated sexes. Results The relative proportions of dioecious species and individuals displayed a unimodal relationship with elevation, both peaking at 750 m a.s.l. In contrast, the relative proportions of monoecious species and individuals displayed negative associations with elevation. In addition, the proportion of dioecious species was significantly positively associated with species richness and with the PD of plots. The proportion of monoecious species was not associated with species richness and was weakly positively associated with the PD of plots. Main conclusions Our results suggest that the selection pressure of elevation on the pollinators and life‐history strategy of dioecious species is much less than expected, and is overshadowed by the as‐yet unexplained correlation of dioecy with species richness. Additional studies of the ecology of cosexual and unisexual (i.e. dioecious and monoecious) species along large environmental gradients are needed.     

Non‐random patterns of invasion and extinction reduce phylogenetic diversity in island bird assemblages

Anthropogenically driven changes in bird communities on oceanic islands exemplify the biotic upheaval experienced by island floras and faunas. While the influence of invasions and extinctions on species richness and beta‐diversity of island bird assemblages has been explored, little is known about the impact of these invasions and extinctions on phylogenetic diversity. Here we quantify phylogenetic diversity of island bird assemblages resulting from extinctions alone, invasions alone, and the combination of extinctions and invasions in the historic time period (1500 CE to the current), and compare it to the expected phylogenetic diversity that would result if these processes involved randomly selected island bird species. We assessed phylogenetic diversity and structure at the scale of the island (n = 152), the archipelago containing the islands (n = 22), and the four oceans containing the archipelagos using three measures. We found that extinction, invasion, and the combination of invasion and extinction generally resulted in lower phylogenetic diversity than expected, regardless of the spatial scale examined. We conclude that extinction and invasion of birds on islands are non‐random with respect to phylogeny and that these processes generally leave bird assemblages with lower phylogenetic diversity than we would expect under random invasion or extinction.     

Spatial patterns of phylogenetic diversity

Ecologists and conservation biologists have historically used species-area and distance-decay relationships as tools to predict the spatial distribution of biodiversity and the impact of habitat loss on biodiversity. These tools treat each species as evolutionarily equivalent, yet the importance of species' evolutionary history in their ecology and conservation is becoming increasingly evident. Here, we provide theoretical predictions for phylogenetic analogues of the species-area and distance-decay relationships. We use a random model of community assembly and a spatially explicit flora dataset collected in four Mediterranean-type regions to provide theoretical predictions for the increase in phylogenetic diversity - the total phylogenetic branch-length separating a set of species - with increasing area and the decay in phylogenetic similarity with geographic separation. These developments may ultimately provide insights into the evolution and assembly of biological communities, and guide the selection of protected areas.     

The phylogenetic diversity of metagenomes

Phylogenetic diversity--patterns of phylogenetic relatedness among organisms in ecological communities--provides important insights into the mechanisms underlying community assembly. Studies that measure phylogenetic diversity in microbial communities have primarily been limited to a single marker gene approach, using the small subunit of the rRNA gene (SSU-rRNA) to quantify phylogenetic relationships among microbial taxa. In this study, we present an approach for inferring phylogenetic relationships among microorganisms based on the random metagenomic sequencing of DNA fragments. To overcome challenges caused by the fragmentary nature of metagenomic data, we leveraged fully sequenced bacterial genomes as a scaffold to enable inference of phylogenetic relationships among metagenomic sequences from multiple phylogenetic marker gene families. The resulting metagenomic phylogeny can be used to quantify the phylogenetic diversity of microbial communities based on metagenomic data sets. We applied this method to understand patterns of microbial phylogenetic diversity and community assembly along an oceanic depth gradient, and compared our findings to previous studies of this gradient using SSU-rRNA gene and metagenomic analyses. Bacterial phylogenetic diversity was highest at intermediate depths beneath the ocean surface, whereas taxonomic diversity (diversity measured by binning sequences into taxonomically similar groups) showed no relationship with depth. Phylogenetic diversity estimates based on the SSU-rRNA gene and the multi-gene metagenomic phylogeny were broadly concordant, suggesting that our approach will be applicable to other metagenomic data sets for which corresponding SSU-rRNA gene sequences are unavailable. Our approach opens up the possibility of using metagenomic data to study microbial diversity in a phylogenetic context.     

Patterns of species diversity and phylogenetic structure of vascular plants on the Qinghai‐Tibetan Plateau

Large-scale patterns of species richness and the underlying mechanisms regulating these patterns have long been the central issues in biogeography and macroecology. Phylogenetic community structure is a result of combined effects of contemporary ecological interactions, environmental filtering, and evolutionary history, and it links community ecology with biogeography and trait evolution. The Qinghai-Tibetan Plateau provides a good opportunity to test the influence of contemporary climate on shaping species richness because of its unique geological history, cold climate, and high biodiversity. In this study, based on high-resolution distributions of ˜9000 vascular plant species, we explored how species richness and phylogenetic structure of vascular plants correlate with climates on the highest (and species rich) plateau on the Earth. The results showed that most of the vascular plants were distributed on the eastern part of the plateau; there was a strong association between species richness and climate, even after the effects of habitat heterogeneity were controlled. However, the responses of richness to climate remarkably depended on life-forms. Richness of woody plants showed stronger climatic associations than that of herbaceous plants; energy and water availability together regulated richness pattern of woody plants; whereas water availability predominantly regulated richness pattern of herbaceous plants. The phylogenetic structure of vascular species clustered in most areas of the plateau, suggesting that rapid speciation and environment filtering dominated the assembly of communities on the plateau. We further propose that biodiversity conservation in this area should better take into account ecological features for different life-forms and phylogenetic lineages.     

Restoring phylogenetic diversity through facilitation

Phylogenetic diversity enhances ecosystem functioning but restoration ecology has not taken advantage of this knowledge. We propose plant facilitation as a mechanism to promote phylogenetic diversity in restoration practices. We planted three functionally different species (Gypsophila struthium, Sedum album, and Limonium sucronicum) in a degraded gypsum ecosystem in Spain and found that after 7 years, the species with nurse traits (G. struthium) survived longer and facilitated the establishment of new species forming phylogenetically diverse neighborhoods. These facilitation‐driven phylodiverse communities may potentially produce a cascade of benefits on ecosystem functioning.