Evolutionary Models and Phylogenetic Signal Assessment via Mantel Test

High-throughput characterization of bacterial responses to complex mixtures of chemical pollutants

Our understanding of how microbes respond to micropollutants, such as pesticides, is almost wholly based on single-species responses to individual chemicals. However, in natural environments, microbes experience multiple pollutants simultaneously. Here we perform a matrix of multi-stressor experiments by assaying the growth of model and non-model strains of bacteria in all 255 combinations of 8 chemical stressors (antibiotics, herbicides, fungicides and pesticides). We found that bacterial strains responded in different ways to stressor mixtures, which could not be predicted simply from their phylogenetic relatedness. Increasingly complex chemical mixtures were both more likely to negatively impact bacterial growth in monoculture and more likely to reveal net interactive effects. A mixed co-culture of strains proved more resilient to increasingly complex mixtures and revealed fewer interactions in the growth response. These results show predictability in microbial population responses to chemical stressors and could increase the utility of next-generation eco-toxicological assays.

Benchmarking the Mantel test and derived methods for testing association between distance matrices

Testing the association between objects is central in ecology, evolution, and quantitative sciences in general. Two types of variables can describe the relationships between objects: point variables (measured on individual objects), and distance variables (measured between pairs of objects). The Mantel test and derived methods have been extensively used for distance variables. Yet, these methods have been criticized due to low statistical power and inflated type I error when spatial autocorrelation is present. Here, we assessed the statistical power between different types of tested variables and the type I error rate over a wider range of autocorrelation intensities than previously assessed, both on univariate and multivariate data. We also illustrated the performance of distance matrix statistics through computational simulations of genetic diversity. We show that the Mantel test and derived methods are not affected by inflated type I error when spatial autocorrelation affects only one variable when investigating correlations, or when either the response or the explanatory variable(s) is affected by spatial autocorrelation while investigating causal relationships. As previously noted, with autocorrelation affecting more variables, inflated type I error could be reduced by modifying the significance threshold. Additionally, the Mantel test has no problem of statistical power when the hypothesis is formulated in terms of distance variables. We highlight that transformation of variable types should be avoided because of the potential information loss and modification of the tested hypothesis. We propose a set of guidelines to help choose the appropriate method according to the type of variables and defined hypothesis.

Machine learning and statistical classification of birdsong link vocal acoustic features with phylogeny

Birdsong is a longstanding model system for studying evolution and biodiversity. Here, we collected and analyzed high quality song recordings from seven species in the family Estrildidae. We measured the acoustic features of syllables and then used dimensionality reduction and machine learning classifiers to identify features that accurately assigned syllables to species. Species differences were captured by the first 3 principal components, corresponding to basic frequency, power distribution, and spectrotemporal features. We then identified the measured features underlying classification accuracy. We found that fundamental frequency, mean frequency, spectral flatness, and syllable duration were the most informative features for species identification. Next, we tested whether specific acoustic features of species' songs predicted phylogenetic distance. We found significant phylogenetic signal in syllable frequency features, but not in power distribution or spectrotemporal features. Results suggest that frequency features are more constrained by species' genetics than are other features, and are the best signal features for identifying species from song recordings. The absence of phylogenetic signal in power distribution and spectrotemporal features suggests that these song features are labile, reflecting learning processes and individual recognition.

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.

Elevation patterns and critical environmental drivers of the taxonomic, functional, and phylogenetic diversity of small mammals in a karst mountain area

Understanding how biodiversity components are related under different environmental factors is a fundamental challenge for ecology studies, yet there is little knowledge of this interplay among the biotas, especially small mammals, in karst mountain areas. Here, we examine the elevation patterns of the taxonomic diversity (TD), phylogenetic diversity (PD), and functional diversity (FD) of small mammals in a karst mountain area, the Wuling Mountains, Southwest China, and compare these patterns between taxa (Rodentia and Eulipotyphla) and scales (broad‐ and narrow‐range species). We also disentangle the impacts of the human influence index, net primary productivity (NPP), normalized difference vegetation index (NDVI), annual precipitation (AP), and annual mean temperature (AMT) on these three facets of biodiversity by using structural equation modeling. We recorded a total of 39 small mammal species, including 26 rodents and 13 species of the order Eulipotyphla. Our study shows that the facets of biodiversity are spatially incongruent. Net primary productivity has a positive effect on the three facets for most groups, while the effect of the NDVI is negative for TD and PD in most groups. AMT temperature and AP have negative effects on FD and PD, whereas TD is dependent on the species range scale. The human influence index effect on TD and PD also depends on the species range scale. These findings provide robust evidence that the ecological drivers of biodiversity differ among different biotas and different range scales, and future research should use multifacet approach to determine biodiversity conservation strategies.

What makes a fang? Phylogenetic and ecological controls on tooth evolution in rear-fanged snakes

BACKGROUND

Fangs are a putative key innovation that revolutionized prey capture and feeding in snakes, and - along with their associated venom phenotypes - have made snakes perhaps the most medically-significant vertebrate animals. Three snake clades are known for their forward-positioned fangs, and these clades (Elapidae, Viperidae, and Atractaspidinae) contain the majority of snakes that are traditionally considered venomous. However, many other snakes are "rear-fanged": they possess potentially venom-delivering teeth situated at the rear end of the upper jaw. Quantification of fang phenotypes - and especially those of rear-fanged species - has proved challenging or impossible owing to the small size and relative rarity of many such snakes. Consequently, it has been difficult to understand the evolutionary history of both venom and prey-capture strategies across extant snakes. We quantified variation in the dentition of 145 colubriform ("advanced") snake species using microCT scanning and compared dental characters with ecological data on species' diet and prey capture method(s) to understand broader patterns in snake fang evolution.

RESULTS

Dental traits such as maxilla length, tooth number, and fang size show strong phylogenetic signal across Colubriformes. We find extreme heterogeneity and evolutionary lability in the rear-fanged phenotype in colubrid (colubrine, dipsadine, and natricine lineages) and lamprophiid snakes, in contrast to relative uniformity in the front fanged phenotypes of other groups (vipers and, to a lesser extent, elapids). Fang size and position are correlated with venom-use in vipers, elapids, and colubrid snakes, with the latter group shifting fangs anteriorly by shortening the entire maxillary bone. We find that maxilla length and tooth number may also be correlated with the evolution of dietary specialization. Finally, an ancestral state reconstruction suggests that fang loss is a recurring phenomenon in colubrid snakes, likely accompanied by shifts in diet and prey capture mode.

CONCLUSIONS

Our study provides a framework for quantifying the complex morphologies associated with venom use in snakes. Our results suggest that fang phenotypes, and particularly the rear-fanged phenotype, in snakes are both diverse and labile, facilitating a wide range of ecological strategies and contributing to spectacular radiations of these organisms in tropical and subtropical biomes worldwide.

Global patterns and determinants of lake macrophyte taxonomic, functional and phylogenetic beta diversity

Documenting the patterns of biological diversity on Earth has always been a central challenge in macroecology and biogeography. However, for the diverse group of freshwater plants, such research program is still in its infancy. Here, we examined global variation in taxonomic, functional and phylogenetic beta diversity patterns of lake macrophytes using regional data from six continents. A data set of ca. 480 lake macrophyte community observations, together with climatic, geographical and environmental variables, was compiled across 16 regions worldwide. We (a) built the very first phylogeny comprising most freshwater plant lineages; (b) exploited a wide array of functional traits that are important to macrophyte autoecology or that relate to lake ecosystem functioning; (c) assessed if different large-scale beta diversity patterns show a clear latitudinal gradient from the equator to the poles using null models; and (d) employed evolutionary and regression models to first identify the degree to which the studied functional traits show a phylogenetic signal, and then to estimate community-environment relationships at multiple spatial scales. Our results supported the notion that ecological niches evolved independently of phylogeny in macrophyte lineages worldwide. We also showed that taxonomic and phylogenetic beta diversity followed the typical global trend with higher diversity in the tropics. In addition, we were able to confirm that species, multi-trait and lineage compositions were first and foremost structured by climatic conditions at relatively broad spatial scales. Perhaps more importantly, we showed that large-scale processes along latitudinal and elevational gradients have left a strong footprint in the current diversity patterns and community-environment relationships in lake macrophytes. Overall, our results stress the need for an integrative approach to macroecology, biogeography and conservation biology, combining multiple diversity facets at different spatial scales.

Untangling the assembly of macrophyte metacommunities by means of taxonomic, functional and phylogenetic beta diversity patterns

Metacommunity ecology has broadened considerably with the recognition that measuring beta diversity beyond the purely taxonomic viewpoint may improve our understanding of the dispersal- and niche-based mechanisms across biological communities. In that perspective, we applied a novel multidimensional approach including taxonomic, functional and phylogenetic data to enhance our basic understanding of macrophyte metacommunity dynamics. For each beta diversity metric, we calculated the mean overall value and tested whether the mean value was different from that expected by chance using null models. We also employed evolutionary and spatially constrained models to first identify the degree to which the studied functional traits showed a phylogenetic signal, and then to estimate the relative importance of spatial and environmental effects on metacommunity structure. We first found that most individual ponds were inhabited by species that were merely random draws from the taxonomic and phylogenetic species pool available in the study region. Contrary to our expectations, not all measured traits were conserved along the phylogeny. We also showed that trait and phylogenetic dimensions strongly increased the amount of variation in beta diversity that can be explained by degree of environmental filtering and dispersal limitation. This suggests that accounting for functional traits and phylogeny in metacommunity ecology helps to explain idiosyncratic patterns of variation in macrophyte species distribution. Importantly, phylogenetic and functional analyses identified the influence of underlying mechanisms that would otherwise be missed in an analysis of taxonomic turnover. Together, these results let us conclude that macrophyte species have labile functional traits adapted to dispersal-based processes and some evolutionary trade-offs that drive community assembly via species sorting. Overall, our exploration of different facets of beta diversity showed how functional and phylogenetic information may be used with species-level data to test community assembly hypotheses that are more ecologically meaningful than assessments of environmental patterns based on the purely taxonomic viewpoint.

Asynchronous evolution of interdependent nest characters across the avian phylogeny

Nest building is a widespread behavior among birds that reflects their adaptation to the environment and evolutionary history. However, it remains unclear how nests evolve and how their evolution relates to the bird phylogeny. Here, by examining the evolution of three nest characters—structure, site, and attachment—across all bird families, we reveal that nest characters did not change synchronically across the avian phylogeny but had disparate evolutionary trajectories. Nest structure shows stronger phylogenetic signal than nest site, while nest attachment has little variation. Nevertheless, the three characters evolved interdependently. For example, the ability of birds to explore new nest sites might depend on the emergence of novel nest structure and/or attachment. Our results also reveal labile nest characters in passerines compared with other birds. This study provides important insights into avian nest evolution and suggests potential associations between nest diversification and the adaptive radiations that generated modern bird lineages.

Nearly all bird species build nests; however, there is extensive variation in both how and where their nests are built. Here, Fang and colleagues reconstruct the evolution of nest structure, nest site and nest attachment across all extant bird families.

Analyzing community-weighted trait means across environmental gradients: should phylogeny stay or should it go?

Functional traits mediate ecological responses of organisms to the environment, determining community structure. Community-weighted trait means (CWM) are often used to characterize communities by combining information on species traits and distribution. Relating CWM variation to environmental gradients allows for evaluating species sorting across the metacommunity, either based on correlation tests or ordinary least squares (OLS) models. Yet, it is not clear if phylogenetic signal in both traits and species distribution affect those analyses. On one hand, phylogenetic signal might indicate niche conservatism along clade evolution, reinforcing the environmental signal in trait assembly patterns. On the other hand, it might introduce phylogenetic autocorrelation to mean trait variation among communities. Under this latter scenario, phylogenetic signal might inflate type I error in analysis relating CWM variation to environmental gradients. We explore multiple ways phylogenetic history may influence analysis relating CWM to environmental gradients. We propose the concept of neutral trait diffusion, which predicts that for a functional trait x, CWM variation among local communities does not deviate from the expectation that x evolved according to a neutral evolutionary process. Based on this framework we introduce a graphical tool called neutral trait diffusion representation (NTDR) that allows for the evaluation of whether it is necessary to carry out phylogenetic correction in the trait prior to analyzing the association between CWM and environmental gradients. We illustrate the NTDR approach using simulated traits, phylogenies and metacommunities. We show that even under moderate phylogenetic signal in both the trait used to define CWM and species distribution across communities, OLS models relating CWM variation to environmental gradients lead to inflated type I error when testing the null hypothesis of no association between CWM and environmental gradient. To overcome this issue, we propose a phylogenetic correction for OLS models and evaluate its statistical performance (type I error and power). Phylogeny-corrected OLS models successfully control for type I error in analysis relating CWM variation to environmental gradients but may show decreased power. Combining the exploratory tool of NTDR and phylogenetic correction in traits, when necessary, guarantees more precise inferences about the environmental forces driving trait-mediated species sorting across metacommunities.

Morphology captures diet and locomotor types in rodents

To understand the functional meaning of morphological features, we need to relate what we know about morphology and ecology in a meaningful, quantitative framework. Closely related species usually share more phenotypic features than distant ones, but close relatives do not necessarily have the same ecologies. Rodents are the most diverse group of living mammals, with impressive ecomorphological diversification. We used museum collections and ecological literature to gather data on morphology, diet and locomotion for 208 species of rodents from different bioregions to investigate how morphological similarity and phylogenetic relatedness are associated with ecology. After considering differences in body size and shared evolutionary history, we find that unrelated species with similar ecologies can be characterized by a well-defined suite of morphological features. Our results validate the hypothesized ecological relevance of the chosen traits. These cranial, dental and external (e.g. ears) characters predicted diet and locomotion and showed consistent differences among species with different feeding and substrate use strategies. We conclude that when ecological characters do not show strong phylogenetic patterns, we cannot simply assume that close relatives are ecologically similar. Museum specimens are valuable records of species' phenotypes and with the characters proposed here, morphology can reflect functional similarity, an important component of community ecology and macroevolution.

Morphology captures diet and locomotor types in rodents

To understand the functional meaning of morphological features, we need to relate what we know about morphology and ecology in a meaningful, quantitative framework. Closely related species usually share more phenotypic features than distant ones, but close relatives do not necessarily have the same ecologies. Rodents are the most diverse group of living mammals, with impressive ecomorphological diversification. We used museum collections and ecological literature to gather data on morphology, diet and locomotion for 208 species of rodents from different bioregions to investigate how morphological similarity and phylogenetic relatedness are associated with ecology. After considering differences in body size and shared evolutionary history, we find that unrelated species with similar ecologies can be characterized by a well-defined suite of morphological features. Our results validate the hypothesized ecological relevance of the chosen traits. These cranial, dental and external (e.g. ears) characters predicted diet and locomotion and showed consistent differences among species with different feeding and substrate use strategies. We conclude that when ecological characters do not show strong phylogenetic patterns, we cannot simply assume that close relatives are ecologically similar. Museum specimens are valuable records of species' phenotypes and with the characters proposed here, morphology can reflect functional similarity, an important component of community ecology and macroevolution.