Ecological traits influence the phylogenetic structure of bird species co-occurrences worldwide

Species-level drivers of avian centrality within seed-dispersal networks across different levels of organisation

Bird-plant seed-dispersal networks are structural components of ecosystems. The role of bird species in seed-dispersal networks (from less [peripheral] to more connected [central]), determines the interaction patterns and their ecosystem services. These roles may be driven by morphological and functional traits as well as evolutionary, geographical and environmental properties acting at different spatial extents. It is still unknown if such drivers are equally important in determining species centrality at different network levels, from individual local networks to the global meta-network representing interactions across all local networks. Using 308 networks covering five continents and 11 biogeographical regions, we show that at the global meta-network level species' range size was the most important driver of species centrality, with more central species having larger range sizes, which would facilitate the interaction with a higher number of plants and thus the maintenance of seed-dispersal interactions. At the local network level, body mass was the only driver with a significant effect, implying that local factors related to resource availability are more important at this level of network organisation than those related to broad spatial factors such as range sizes. This could also be related to the mismatch between species-level traits, which do not consider intraspecific variation, and the local networks that can depend on such variation. Taken together, our results show that the drivers determining species centrality are relative to the levels of network organisation, suggesting that prediction of species functional roles in seed-dispersal interactions requires combined local and global approaches.

Global impacts of climate change on avian functional diversity

Climate change is predicted to drive geographical range shifts, leading to fluctuations in species richness (SR) worldwide. However, the effect of these changes on functional diversity (FD) remains unclear, in part because comprehensive species-level trait data are generally lacking at global scales. Here, we use morphometric and ecological traits for 8268 bird species to estimate the impact of climate change on avian FD. We show that future bird assemblages are likely to undergo substantial shifts in trait structure, with a magnitude of change greater than predicted from SR alone, and a direction of change varying according to geographical location and trophic guild. For example, our models predict that FD of insect predators will increase at higher latitudes with concurrent losses at mid-latitudes, whereas FD of seed dispersing birds will fluctuate across the tropics. Our findings highlight the potential for climate change to drive continental-scale shifts in avian FD with implications for ecosystem function and resilience.

How deregulation, drought and increasing fire impact Amazonian biodiversity

Biodiversity contributes to the ecological and climatic stability of the Amazon Basin^1,2, but is increasingly threatened by deforestation and fire^3,4. Here we quantify these impacts over the past two decades using remote-sensing estimates of fire and deforestation and comprehensive range estimates of 11,514 plant species and 3,079 vertebrate species in the Amazon. Deforestation has led to large amounts of habitat loss, and fires further exacerbate this already substantial impact on Amazonian biodiversity. Since 2001, 103,079-189,755 km² of Amazon rainforest has been impacted by fires, potentially impacting the ranges of 77.3-85.2% of species that are listed as threatened in this region⁵. The impacts of fire on the ranges of species in Amazonia could be as high as 64%, and greater impacts are typically associated with species that have restricted ranges. We find close associations between forest policy, fire-impacted forest area and their potential impacts on biodiversity. In Brazil, forest policies that were initiated in the mid-2000s corresponded to reduced rates of burning. However, relaxed enforcement of these policies in 2019 has seemingly begun to reverse this trend: approximately 4,253-10,343 km² of forest has been impacted by fire, leading to some of the most severe potential impacts on biodiversity since 2009. These results highlight the critical role of policy enforcement in the preservation of biodiversity in the Amazon.

Functional diversity of avian communities increases with canopy height: From individual behavior to continental-scale patterns

Vegetation complexity is an important predictor of animal species diversity. Specifically, taller vegetation should provide more potential ecological niches and thus harbor communities with higher species richness and functional diversity (FD). Resource use behavior is an especially important functional trait because it links species to their resource base with direct relevance to niche partitioning. However, it is unclear how exactly the diversity of resource use behavior changes with vegetation complexity. To address this question, we studied avian FD in relation to vegetation complexity along a continental-scale vegetation gradient. We quantified foraging behavior of passerine birds in terms of foraging method and substrate use at 21 sites (63 transects) spanning 3,000 km of woodlands and forests in Australia. We also quantified vegetation structure on 630 sampling points at the same sites. Additionally, we measured morphological traits for all 111 observed species in museum collections. We calculated individual-based, abundance-weighted FD in morphology and foraging behavior and related it to species richness and vegetation complexity (indexed by canopy height) using structural equation modeling, rarefaction analyses, and distance-based metrics. FD of morphology and foraging methods was best predicted by species richness. However, FD of substrate use was best predicted by canopy height (ranging 10-30 m), but only when substrates were categorized with fine resolution (17 categories), not when categorized coarsely (8 categories). These results suggest that, first, FD might increase with vegetation complexity independently of species richness, but whether it does so depends on the studied functional trait. Second, patterns found might be shaped by how finely we categorize functional traits. More complex vegetation provided larger "ecological space" with more resources, allowing the coexistence of more species with disproportionately more diverse foraging substrate use. We suggest that the latter pattern was driven by nonrandom accumulation of functionally distinct species with increasing canopy height.

Habitat Quality and Social Behavioral Association Network in a Wintering Waterbirds Community

Migratory waterbirds concentrated in freshwater ecosystems in mosaic environments rely on quality habitats for overwintering. At West Dongting Lake National Nature Reserve (WDLNNR), China, land-use change and hydrology alternation are compounding factors that have affected important wintering areas for migratory waterbirds. Presently, changes in the hydrology and landscape have reshaped natural wintering habitats and their availability, though the impact of hydrological management on habitat selection of wintering waterbirds is largely unknown. In this study, we classified differentially managed habitats and calculated their area using the normalized difference vegetation index (NDVI) to evaluate suitable habitat availability over the study period (2016–2017 and 2017–2018 wintering periods). We then used social behavioral association network (SBAN) model to compare habitat quality through species-species social interactions and species-habitat associations in lakes with different hydrological management. The results indicated that social interactions between and within species structured wintering waterbirds communities, which could be dominated by one or more species, while dominant species control the activities of other co-existing species. Analysis of variance (ANOVA) tests indicated significant differences in SBAN metrics between lakes (p = 0.0237) and habitat (p < 0.0001) levels. Specifically, lakes with managed hydrology were preferred by more species. The managed lakes had better habitat quality in terms of significantly higher habitat areas (p < 0.0001) and lower habitat transitions (p = 0.0113). Collectively, our findings suggest that proper hydrological management can provide continuous availability of quality habitats, especially mudflats and shallow waters, for a stable SBAN to ensure a wintering waterbirds community with more sympatric species in a dynamic environment.

Targeting Conservation Actions at Species Threat Response Thresholds

Given the failure of the world's governments to improve the status of biodiversity by 2020, a new strategic plan for 2030 is being developed. In order to be successful, a step-change is needed to not just simply halt biodiversity loss, but to bend the curve of biodiversity loss to stable or increasing species' populations. Here, we propose a framework that quantifies species' responses across gradients of threat intensity to implement more efficient and better targeted conservation actions. Our framework acknowledges the variation in threat intensities as well as the differences among species in their capacity to respond, and is implemented at a relevant scale for national and international policy-making.

Global functional and phylogenetic structure of avian assemblages across elevation and latitude

Mountain systems are exceptionally species rich, yet the associated elevational gradients in functional and phylogenetic diversity and their consistency across latitude remain little understood. Here, we document how avian functional and phylogenetic diversity and structure vary along all major elevational gradients worldwide and uncover strong latitudinal differences. Assemblages in warm tropical lowlands and cold temperate highlands are marked by high functional overdispersion and distinctiveness, whereas tropical highlands and temperate lowlands appear strongly functionally clustered and redundant. We additionally find strong geographic variation in the interplay of phylogenetic and functional structure, with strongest deviations between the two in temperate highlands. This latitudinal and elevational variation in assemblage functional structure is underpinned by nuanced shifts in the position, shape and composition of multivariate trait space. We find that, independent of latitude, high-elevation assemblages emerge as exceptionally susceptible to functional change.

Trait-environment relationships differ between mixed-species flocking and nonflocking bird assemblages

Hypotheses about the mechanisms of community assembly suggest that biotic and abiotic filters constrain species establishment through selection on their functional traits. It is unclear how differences in traits influence the niche dimensions of closely related bird species when they coexist in spatiotemporally heterogeneous environments. Further, it is necessary to take into account their participation in mixed-species flocks, social systems that can include both competition and facilitation. For 6 yr, we conducted counts of forest bird species and took measurements of environmental variables along an elevational gradient in the Nanling Mountains, China. To disentangle different deterministic and historical/stochastic processes between flocking and nonflocking bird assemblages, we first compared phylogenetic and functional structure, and community-weighted mean trait values (CWM). We further assessed elevational variations in trait-environment relationships. We found that the flocking and nonflocking bird assemblages were structured by environmental gradients in contrasting ways. The nonflocking assemblage showed a strong change from over-dispersed to clustered community structure with increasing elevations, consistent with the strong selective pressures of a harsh environment (i.e., environmental filtering). The nonflocking assemblage also displayed significant trait-environment relationships in bivariate correlations and multivariate ordination space, including specific morphological and foraging traits that are linked to vegetation characteristics (e.g., short trees at high elevations). By contrast, flocking birds were more resilient to habitat change with elevation, with relatively consistent community membership, and showed fewer trait-environment associations. CWM of traits that are known to be associated with species' propensity to join mixed-species flocks, including small body size and broad habitat specificity, were linked to the flocking assemblage consistently across the elevational gradient. Collectively, our trait-based analyses provide strong evidence that trait-environment relationships differ between flocking and nonflocking bird assemblages. Besides serving as bellwethers of changing environments, emergent properties of flock systems may increase the resilience of animal communities undergoing environmental change. Mixed-species flocks present an ideal model with which to explore cooccurrence of closely related species, because habitat filtering may be buffered, and the patterns observed are therefore the outcomes of species interactions including both competition and facilitation.

Solitary ecology as a phenomenon extending beyond insular systems: exaptive evolution in Anolis lizards

The mechanisms driving phenotypic evolution have been of interest to biologists since Darwin. Ecological release-wherein adaptive evolution occurs following relaxation of constraining selective pressures-and environmental filtering-wherein exaptive traits allow colonization of a new area-have been studied in several insular cases. Anolis lizards, which may exist in solitude or sympatry with multiple congeners, are an excellent system for evaluating whether ecological release and environmental filtering are associated with phenotypic shifts across phylogenetic and geographical scales. Insular solitary Anolis exhibit phenotypic differentiation in body size and sexual size dimorphism-SSD-through exaptive and adaptive evolution, respectively. But, the generality of these effects has not yet been addressed. Here, we analyse the evolution of body size and SSD relative to sympatry in mainland Anolis. We found that mainland species co-occurring with few congeners exhibit uniform body size and greater SSD relative to other random mainland assemblages, consistent with the insular solitary pattern. The locations of evolutionary shifts for both traits do not coincide with evolutionary transitions to decreased levels of sympatry. These results are consistent with exaptive environmental filtering but not adaptive ecological release. Future studies should be conducted at local scales to evaluate the role of these factors in the evolution of solitary existence in mainland and island species.

Drivers of Phylogenetic Assemblage Structure of the Furnariides, a Widespread Clade of Lowland Neotropical Birds

Species co-occurrence in local assemblages is shaped by distinct processes at different spatial and temporal scales. Here we focus on historical explanations and examine the phylogenetic structure of local assemblages of the Furnariides clade (Aves: Passeriformes), assessing the influence of diversification rates on the assembly and species co-occurrence within those assemblages. Using 120 local assemblages across Bolivia and Argentina and a nearly complete phylogeny for the clade, we analyzed assemblage phylogenetic structure, applying a recently developed model (DAMOCLES, or dynamic assembly model of colonization, local extinction, and speciation) accounting for the historical processes of speciation, colonization, and local extinction. We also evaluated how diversification rates determine species co-occurrence. We found that the assembly of Furnariides assemblages can be explained largely by speciation, colonization, and local extinction without invoking current local species interactions. Phylogenetic structure of open habitat assemblages mainly showed clustering, characterized by faster rates of colonization and local extinction than in forest habitats, whereas forest habitat assemblages were congruent with the model's equal rates expectation, thus highlighting the influence of habitat preferences on assembly and co-occurrence patterns. Our results suggest that historical processes are sufficient to explain local assemblage phylogenetic structure, while there is little evidence for species ecological interactions in avian assemblage diversity and composition.

Body size information in large-scale acoustic bat databases

Background

Citizen monitoring programs using acoustic data have been useful for detecting population and community patterns. However, they have rarely been used to study broad scale patterns of species traits. We assessed the potential of acoustic data to detect broad scale patterns in body size. We compared geographical patterns in body size with acoustic signals in the bat species Pipistrellus pipistrellus. Given the correlation between body size and acoustic characteristics, we expected to see similar results when analyzing the relationships of body size and acoustic signals with climatic variables.

Methods

We assessed body size using forearm length measurements of 1,359 bats, captured by mist nets in France. For acoustic analyses, we used an extensive dataset collected through the French citizen bat survey. We isolated each bat echolocation call (n = 4,783) and performed automatic measures of signals, including the frequency of the flattest part of the calls (characteristic frequency). We then examined the relationship between forearm length, characteristic frequencies, and two components resulting from principal component analysis for geographic (latitude, longitude) and climatic variables.

Results

Forearm length was positively correlated with higher precipitation, lower seasonality, and lower temperatures. Lower characteristic frequencies (i.e., larger body size) were mostly related to lower temperatures and northern latitudes. While conducted on different datasets, the two analyses provided congruent results.

Discussion

Acoustic data from citizen science programs can thus be useful for the detection of large-scale patterns in body size. This first analysis offers a new perspective for the use of large acoustic databases to explore biological patterns and to address both theoretical and applied questions.

A taxonomic, functional, and phylogenetic perspective on the community assembly of passerine birds along an elevational gradient in southwest China

Integrating multiple facets of biodiversity to describe spatial and temporal distribution patterns is one way of revealing the mechanisms driving community assembly. We assessed the species, functional, and phylogenetic composition and structure of passerine bird communities along an elevational gradient both in wintering and breeding seasons in the Ailao Mountains, southwest China, in order to identify the dominant ecological processes structuring the communities and how these processes change with elevation and season. Our research confirms that the highest taxonomic diversity, and distinct community composition, was found in the moist evergreen broadleaf forest at high elevation in both seasons. Environmental filtering was the dominant force at high elevations with relatively cold and wet climatic conditions, while the observed value of mean pairwise functional and phylogenetic distances of low elevation was constantly higher than expectation in two seasons, suggested interspecific competition could play the key role at low elevations, perhaps because of relative rich resource result from complex vegetation structure and human-induced disturbance. Across all elevations, there was a trend of decreasing intensity of environmental filtering whereas increasing interspecific competition from wintering season to breeding season. This was likely due to the increased resource availability but reproduction-associated competition in the summer months. In general, there is a clear justification for conservation efforts to protect entire elevational gradients in the Ailao Mountains, given the distinct taxonomic, functional, and phylogenetic compositions and also elevational migration pattern in passerine bird communities.

Phylogenetic and Functional Structure of Wintering Waterbird Communities Associated with Ecological Differences

Ecological differences may be related to community component divisions between Oriental (west) and Sino-Japanese (east) realms, and such differences may result in weak geographical breaks in migratory species that are highly mobile. Here, we conducted comparative phylogenetic and functional structure analyses of wintering waterbird communities in southern China across two realms and subsequently examined possible climate drivers of the observed patterns. An analysis based on such highly migratory species is particularly telling because migration is bound to reduce or completely eliminate any divergence between communities. Phylogenetic and functional structure of eastern communities showed over-dispersion while western communities were clustered. Basal phylogenetic and functional turnover of western communities was significant lower than that of eastern communities. The break between eastern and western communities was masked by these two realms. Geographic patterns were related to mean temperature changes and temperature fluctuations, suggesting that temperature may filter waterbird lineages and traits, thus underlying geographical community divisions. These results suggest phylogenetic and functional divisions in southern China, coinciding with biogeography. This study shows that temperature fluctuations constitute an essential mechanism shaping geographical divisions that have largely gone undetected previously, even under climate change.

The role of environment, dispersal and competition in explaining reduced co-occurrence among related species

The composition of ecological assemblages depends on a variety of factors including environmental filtering, biotic interactions and dispersal limitation. By evaluating the phylogenetic pattern of assemblages, we gain insight into the relative contribution of these mechanisms to generating observed assemblages. We address some limitations in the field of community phylogenetics by using simulations, biologically relevant null models, and cost distance analysis to evaluate simultaneous mechanisms leading to observed patterns of co-occurrence. Building from past studies of phylogenetic community structure, we applied our approach to hummingbird assemblages in the Northern Andes. We compared the relationship between relatedness and co-occurrence among predicted assemblages, based on estimates of suitable habitat and dispersal limitation, and observed assemblages. Hummingbird co-occurrence peaked at intermediate relatedness and decreased when a closely-related species was present. This result was most similar to simulations that included simultaneous effects of phylogenetic conservatism and repulsion. In addition, we found older sister taxa were only weakly more separated by geographic barriers, suggesting that time since dispersal is unlikely to be the sole factor influencing co-occurrence of closely related species. Our analysis highlights the role of multiple mechanisms acting simultaneously, and provides a hypothesis for the potential importance of competition at regional scales.

Climate-driven changes in functional biogeography of Arctic marine fish communities

Climate change triggers poleward shifts in species distribution leading to changes in biogeography. In the marine environment, fish respond quickly to warming, causing community-wide reorganizations, which result in profound changes in ecosystem functioning. Functional biogeography provides a framework to address how ecosystem functioning may be affected by climate change over large spatial scales. However, there are few studies on functional biogeography in the marine environment, and none in the Arctic, where climate-driven changes are most rapid and extensive. We investigated the impact of climate warming on the functional biogeography of the Barents Sea, which is characterized by a sharp zoogeographic divide separating boreal from Arctic species. Our unique dataset covered 52 fish species, 15 functional traits, and 3,660 stations sampled during the recent warming period. We found that the functional traits characterizing Arctic fish communities, mainly composed of small-sized bottom-dwelling benthivores, are being rapidly replaced by traits of incoming boreal species, particularly the larger, longer lived, and more piscivorous species. The changes in functional traits detected in the Arctic can be predicted based on the characteristics of species expected to undergo quick poleward shifts in response to warming. These are the large, generalist, motile species, such as cod and haddock. We show how functional biogeography can provide important insights into the relationship between species composition, diversity, ecosystem functioning, and environmental drivers. This represents invaluable knowledge in a period when communities and ecosystems experience rapid climate-driven changes across biogeographical regions.

Phylogenetic fields through time: temporal dynamics of geographical co-occurrence and phylogenetic structure within species ranges

Species co-occur with different sets of other species across their geographical distribution, which can be either closely or distantly related. Such co-occurrence patterns and their phylogenetic structure within individual species ranges represent what we call the species phylogenetic fields (PFs). These PFs allow investigation of the role of historical processes--speciation, extinction and dispersal--in shaping species co-occurrence patterns, in both extinct and extant species. Here, we investigate PFs of large mammalian species during the last 3 Myr, and how these correlate with trends in diversification rates. Using the fossil record, we evaluate species' distributional and co-occurrence patterns along with their phylogenetic structure. We apply a novel Bayesian framework on fossil occurrences to estimate diversification rates through time. Our findings highlight the effect of evolutionary processes and past climatic changes on species' distributions and co-occurrences. From the Late Pliocene to the Recent, mammal species seem to have responded in an individualistic manner to climate changes and diversification dynamics, co-occurring with different sets of species from different lineages across their geographical ranges. These findings stress the difficulty of forecasting potential effects of future climate changes on biodiversity.

The effect of spatial scale on relative influences of assembly processes in temperate stream fish assemblages

Understanding how biotic and abiotic processes influence community assembly is a fundamental aim in ecology. Although spatial scales at which communities are studied may affect the relative importance of such assembly processes, spatial influences on community assembly have not been thoroughly addressed. We tested how spatial scale affects inferences of habitat filtering and competitive exclusion assembly processes in darter (Percidae: Etheostomatinae) assemblages across four temperate stream systems. We predicted competitive exclusion would influence assembly in fine-scale assemblages, and habitat filtering would be more influential as spatial scale increased. Moreover, we assumed that habitat heterogeneity would increase with scale, and consequently alleviate direct competitive exclusion acting at finer scales. Using a framework that incorporated genetic relatedness, morphological traits, and habitat use among co-occurring darter species, we identified ecological and evolutionary patterns of structure, which allowed us to elucidate processes of assembly. Based on phylogenetic structure, assemblages showed an increase in habitat filtering (i.e., increased phylogenetic clustering) as we scaled up from fine to intermediate assemblages; however, we found mixed signals for habitat filtering and competitive exclusion at the broadest spatial scale. While habitat filtering was found to have an overall high relative importance during assembly, we also found influence of competitive exclusion processes based on limited morphological similarity among co-occurring species. Our results generally support an increased influence of habitat filtering processes in broader scale assemblages. Moreover, we suggest that habitat filtering and competitive exclusion processes act simultaneously during assembly, although the relative influence of each process may be spatial-scale dependent.

Omnivory in birds is a macroevolutionary sink

Diet is commonly assumed to affect the evolution of species, but few studies have directly tested its effect at macroevolutionary scales. Here we use Bayesian models of trait-dependent diversification and a comprehensive dietary database of all birds worldwide to assess speciation and extinction dynamics of avian dietary guilds (carnivores, frugivores, granivores, herbivores, insectivores, nectarivores, omnivores and piscivores). Our results suggest that omnivory is associated with higher extinction rates and lower speciation rates than other guilds, and that overall net diversification is negative. Trait-dependent models, dietary similarity and network analyses show that transitions into omnivory occur at higher rates than into any other guild. We suggest that omnivory acts as macroevolutionary sink, where its ephemeral nature is retrieved through transitions from other guilds rather than from omnivore speciation. We propose that these dynamics result from competition within and among dietary guilds, influenced by the deep-time availability and predictability of food resources.

Diet is known to influence speciation, but much less is known about how this process operates at macroevolutionary scales. Using a global dietary database of birds, Burin et al. show that omnivory is associated with higher extinction and lower speciation rates compared to other guilds.

Testing projected wild bee distributions in agricultural habitats: predictive power depends on species traits and habitat type

Species distribution models (SDM) are increasingly used to understand the factors that regulate variation in biodiversity patterns and to help plan conservation strategies. However, these models are rarely validated with independently collected data and it is unclear whether SDM performance is maintained across distinct habitats and for species with different functional traits. Highly mobile species, such as bees, can be particularly challenging to model. Here, we use independent sets of occurrence data collected systematically in several agricultural habitats to test how the predictive performance of SDMs for wild bee species depends on species traits, habitat type, and sampling technique. We used a species distribution modeling approach parametrized for the Netherlands, with presence records from 1990 to 2010 for 193 Dutch wild bees. For each species, we built a Maxent model based on 13 climate and landscape variables. We tested the predictive performance of the SDMs with independent datasets collected from orchards and arable fields across the Netherlands from 2010 to 2013, using transect surveys or pan traps. Model predictive performance depended on species traits and habitat type. Occurrence of bee species specialized in habitat and diet was better predicted than generalist bees. Predictions of habitat suitability were also more precise for habitats that are temporally more stable (orchards) than for habitats that suffer regular alterations (arable), particularly for small, solitary bees. As a conservation tool, SDMs are best suited to modeling rarer, specialist species than more generalist and will work best in long-term stable habitats. The variability of complex, short-term habitats is difficult to capture in such models and historical land use generally has low thematic resolution. To improve SDMs' usefulness, models require explanatory variables and collection data that include detailed landscape characteristics, for example, variability of crops and flower availability. Additionally, testing SDMs with field surveys should involve multiple collection techniques.

A new dynamic null model for phylogenetic community structure

Phylogenies are increasingly applied to identify the mechanisms structuring ecological communities but progress has been hindered by a reliance on statistical null models that ignore the historical process of community assembly. Here, we address this, and develop a dynamic null model of assembly by allopatric speciation, colonisation and local extinction. Incorporating these processes fundamentally alters the structure of communities expected due to chance, with speciation leading to phylogenetic overdispersion compared to a classical statistical null model assuming equal probabilities of community membership. Applying this method to bird and primate communities in South America we show that patterns of phylogenetic overdispersion – often attributed to negative biotic interactions – are instead consistent with a species neutral model of allopatric speciation, colonisation and local extinction. Our findings provide a new null expectation for phylogenetic community patterns and highlight the importance of explicitly accounting for the dynamic history of assembly when testing the mechanisms governing community structure.

Establishing macroecological trait datasets: digitalization, extrapolation, and validation of diet preferences in terrestrial mammals worldwide

Ecological trait data are essential for understanding the broad-scale distribution of biodiversity and its response to global change. For animals, diet represents a fundamental aspect of species’ evolutionary adaptations, ecological and functional roles, and trophic interactions. However, the importance of diet for macroevolutionary and macroecological dynamics remains little explored, partly because of the lack of comprehensive trait datasets. We compiled and evaluated a comprehensive global dataset of diet preferences of mammals (“MammalDIET”). Diet information was digitized from two global and cladewide data sources and errors of data entry by multiple data recorders were assessed. We then developed a hierarchical extrapolation procedure to fill-in diet information for species with missing information. Missing data were extrapolated with information from other taxonomic levels (genus, other species within the same genus, or family) and this extrapolation was subsequently validated both internally (with a jack-knife approach applied to the compiled species-level diet data) and externally (using independent species-level diet information from a comprehensive continentwide data source). Finally, we grouped mammal species into trophic levels and dietary guilds, and their species richness as well as their proportion of total richness were mapped at a global scale for those diet categories with good validation results. The success rate of correctly digitizing data was 94%, indicating that the consistency in data entry among multiple recorders was high. Data sources provided species-level diet information for a total of 2033 species (38% of all 5364 terrestrial mammal species, based on the IUCN taxonomy). For the remaining 3331 species, diet information was mostly extrapolated from genus-level diet information (48% of all terrestrial mammal species), and only rarely from other species within the same genus (6%) or from family level (8%). Internal and external validation showed that: (1) extrapolations were most reliable for primary food items; (2) several diet categories (“Animal”, “Mammal”, “Invertebrate”, “Plant”, “Seed”, “Fruit”, and “Leaf”) had high proportions of correctly predicted diet ranks; and (3) the potential of correctly extrapolating specific diet categories varied both within and among clades. Global maps of species richness and proportion showed congruence among trophic levels, but also substantial discrepancies between dietary guilds. MammalDIET provides a comprehensive, unique and freely available dataset on diet preferences for all terrestrial mammals worldwide. It enables broad-scale analyses for specific trophic levels and dietary guilds, and a first assessment of trait conservatism in mammalian diet preferences at a global scale. The digitalization, extrapolation and validation procedures could be transferable to other trait data and taxa.

Computing the skewness of the phylogenetic mean pairwise distance in linear time

Background

The phylogenetic Mean Pairwise Distance (MPD) is one of the most popular measures for computing the phylogenetic distance between a given group of species. More specifically, for a phylogenetic tree and for a set of species R represented by a subset of the leaf nodes of , the MPD of R is equal to the average cost of all possible simple paths in that connect pairs of nodes in R.

Among other phylogenetic measures, the MPD is used as a tool for deciding if the species of a given group R are closely related. To do this, it is important to compute not only the value of the MPD for this group but also the expectation, the variance, and the skewness of this metric. Although efficient algorithms have been developed for computing the expectation and the variance the MPD, there has been no approach so far for computing the skewness of this measure.

Results

In the present work we describe how to compute the skewness of the MPD on a tree optimally, in Θ(n) time; here n is the size of the tree . So far this is the first result that leads to an exact, let alone efficient, computation of the skewness for any popular phylogenetic distance measure. Moreover, we show how we can compute in Θ(n) time several interesting quantities in , that can be possibly used as building blocks for computing efficiently the skewness of other phylogenetic measures.

Conclusions

The optimal computation of the skewness of the MPD that is outlined in this work provides one more tool for studying the phylogenetic relatedness of species in large phylogenetic trees. Until now this has been infeasible, given that traditional techniques for computing the skewness are inefficient and based on inexact resampling.