A guide to phylogenetic metrics for conservation, community ecology and macroecology

Non-native plants tend to be phylogenetically distant but functionally similar to native plants under intense disturbance at the Three Gorges Reservoir Area

Darwin's two opposing hypotheses, proposing that non-native species closely or distantly related to native species are more likely to succeed, are known as 'Darwin's Naturalization Conundrum'. Recently, invasion ecologists have sought to unravel these hypotheses. Studies that incorporate rich observational data in disturbed ecosystems that integrate phylogenetic and functional perspectives have potential to shed light on the conundrum. Using 313 invaded plant communities including 46 invasive plant species and 531 native plant species across the Three Gorges Reservoir Area in China, we aim to evaluate the coexistence mechanisms of invasive and native plants by integrating phylogenetic and functional dimensions at spatial and temporal scales. Our findings revealed that invasive plants tended to co-occur more frequently with native plant species that were phylogenetically distant but functionally similar in the reservoir riparian zone. Furthermore, our study demonstrated that the filtering of flood-dry-flood cycles played a significant role in deepening functional similarities of native communities and invasive-native species over time. Our study highlights the contrasting effects of phylogenetic relatedness and functional similarity between invasive and native species in highly flood-disturbed habitats, providing new sights into Darwin's Naturalization Conundrum.

Phylogenetic tree statistics: A systematic overview using the new R package 'treestats'

Phylogenetic trees are believed to contain a wealth of information on diversification processes. However, comparing phylogenetic trees is not straightforward due to their high dimensionality. Researchers have therefore defined a wide range of low-dimensional summary statistics. Currently, it remains unexplored to what extent these summary statistics cover the same underlying information and what summary statistics best explain observed variation across phylogenies. Furthermore, a large subset of available summary statistics focusses on measuring the topological features of a phylogenetic tree, but are often only explored at the extreme edge cases of the fully balanced or imbalanced tree and not for trees of intermediate balance. Here, we introduce a new R package called 'treestats', that provides speed optimized code to compute 70 summary statistics. We study correlations between summary statistics on empirical trees and on trees simulated using several diversification models. Furthermore, we introduce an algorithm to create intermediately balanced trees in a well-defined manner, in order to explore variation in summary statistics across a balance gradient. We find that almost all summary statistics are correlated with tree size, and find that it is difficult, if not impossible, to correct for tree size, unless the tree generating model is known. Furthermore, we find that across empirical and simulated trees, at least three large clusters of correlated summary statistics can be found, where statistics group together based on information used (topology or branching times). However, the finer grained correlation structure appears to depend strongly on either the taxonomic group studied (in empirical studies) or the tree generating model (in simulation studies). Amongst statistics describing the (im)balance of a tree, we find that almost all statistics vary non-linearly, and sometimes even non-monotonically, with our generated balance gradient. This indicates that balance is perhaps a more complex property of a tree than previously thought. Furthermore, using our new imbalancing algorithm, we devise a numerical test to identify balance statistics, and identify several statistics as balance statistics that were not previously considered as such. Lastly, our results lead to several recommendations on which statistics to select when analyzing and comparing phylogenetic trees.

The global loss of avian functional and phylogenetic diversity from anthropogenic extinctions

Humans have been driving a global erosion of species richness for millennia, but the consequences of past extinctions for other dimensions of biodiversity-functional and phylogenetic diversity-are poorly understood. In this work, we show that, since the Late Pleistocene, the extinction of 610 bird species has caused a disproportionate loss of the global avian functional space along with ~3 billion years of unique evolutionary history. For island endemics, proportional losses have been even greater. Projected future extinctions of more than 1000 species over the next two centuries will incur further substantial reductions in functional and phylogenetic diversity. These results highlight the severe consequences of the ongoing biodiversity crisis and the urgent need to identify the ecological functions being lost through extinction.

Relationship Between Evolutionary Diversity and Aboveground Biomass During 150 Years of Natural Vegetation Regeneration in Temperate China

While the link between plant species diversity and biomass has been well-studied, the impact of evolutionary diversity on community biomass across long timescales or ongoing change remains a subject of debate. We elucidated the association between evolutionary diversity and community aboveground biomass (AGB) using an ideal experimental system with over 150-year history of natural vegetation regeneration. Higher phylogenetic diversity facilitated the sampling effect under the influence of environmental filtering, and caused an increase in AGB. Phylogenetic structure varied from aggregation to dispersion during the later period of vegetation recovery (70-150 years), which was correlated with increases in niche complementarity and increasing AGB. Woody plant evolutionary diversity was used as a key to predict the relationship between vegetation recovery and AGB, with a total explanatory power of ~84.7%. Mixed forests composed of evergreen conifers and deciduous broadleaf forests had higher carbon sequestration potential than that of pure forests, which is advantageous for increasing top-stage AGB. This research expands our knowledge of the causes and effects of biodiversity and ecosystem function dynamics over time and space, which is important for accurately predicting future climate change effects.

Inferring the Evolutionary Model of Community-Structuring Traits with Convolutional Kitchen Sinks

When communities are assembled through processes such as filtering or limiting similarity acting on phylogenetically conserved traits, the evolutionary signature of those traits may be reflected in patterns of community membership. We show how the model of trait evolution underlying community-structuring traits can be inferred from community membership data using both a variation of a traditional eco-phylogenetic metric-the mean pairwise phylogenetic distance (MPD) between taxa-and a recent machine learning tool, Convolutional Kitchen Sinks (CKS). Both methods perform well across a range of phylogenetically informative evolutionary models, but CKS outperforms MPD as tree size increases. We demonstrate CKS by inferring the evolutionary history of freeze tolerance in angiosperms. Our analysis is consistent with a late burst model, suggesting freeze tolerance evolved recently. We suggest that multiple data types that are ordered on phylogenies, such as trait values, species interactions, or community presence/absence, are good candidates for CKS modeling because the generative models produce structured differences between neighboring points that CKS is well-suited for. We introduce the R package kitchen to perform CKS for generic application of the technique.

Strong and consistent effects of waterbird composition on HPAI H5 occurrences across Europe

Since 2014, highly pathogenic avian influenza (HPAI) H5 viruses of clade 2.3.4.4 have been dominating the outbreaks across Europe, causing massive deaths among poultry and wild birds. However, the factors shaping these broad-scale outbreak patterns, especially those related to waterbird community composition, remain unclear. In particular, we do not know whether these risk factors differ from those of other H5 clades. Addressing this knowledge gap is important for predicting and preventing future HPAI outbreaks. Using extensive waterbird survey datasets from about 6883 sites, we here explored the effect of waterbird community composition on HPAI H5Nx (clade 2.3.4.4) spatial patterns in the 2016/2017 and 2020/2021 epidemics in Europe, and compared it with the 2005/2006 HPAI H5N1 (clade 2.2) epidemic. We showed that HPAI H5 occurrences in wild birds in the three epidemics were strongly associated with very similar waterbird community attributes, which suggested that, in nature, similar interspecific transmission processes operate between the HPAI H5 subtypes or clades. Importantly, community phylogenetic diversity consistently showed a negative association with H5 occurrence in all three epidemics, suggesting a dilution effect of phylogenetic diversity. In contrast, waterbird community variables showed much weaker associations with HPAI H5Nx occurrence in poultry. Our results demonstrate that models based on previous epidemics can predict future HPAI H5 patterns in wild birds, implying that it is important to include waterbird community factors in future HPAI studies to predict outbreaks and improve surveillance activities.

Time since fire as a driver of taxonomic and phylogenetic patterns of grassland plant communities

Fire plays a key role in grasslands, determining the distribution and evolution of species and boundaries with neighboring ecosystems. Evidence of community-wide responses to fire is largely based on taxonomic and functional descriptors, while the phylogenetic dimension is overlooked. Here we evaluated how the taxonomic and phylogenetic structure of grassland plant communities responded to a time since fire (TSF) gradient. We sampled 12 communities in Southern Brazil under varying TSF and calculated taxonomic species richness (S) and dominance (D), phylogenetic diversity (PD), and mean phylogenetic distances (MPD). We used Structural Equation Models to test the relationships between the environmental gradient and community descriptors. Communities with longer TSF presented higher PD and MPD but lower species richness and increased taxonomic dominance. These sites were dominated by monocots, specifically C4 grasses, but also presented exclusive clades, whereas recently-burned sites presented lower taxonomic dominance and more species distributed in a wider variety of clades. Our results indicate that these scenarios are interchangeable and dependent on fire management. Fire adaptation was not constrained by phylogenetic relatedness, contrasting with previous findings for tropical savannahs and indicating that temperate and tropical non-forest ecosystems from South America respond differently to fire, possibly due to different evolutionary histories.

Urbanization reduces diversity, simplifies community and filter bird species based on their functional traits in a tropical city

Understanding how organisms are coping with major changes imposed by urban intensification is a complex task. In fact, our understanding of the impacts of urbanization on biodiversity is scarce in the global south compared to the north. In this study, we evaluated how bird communities are affected by impact of urban intensification in a tropical city. Thus, we assessed whether increased urban intensification 1) jeopardizes bird diversity (taking into account taxonomic-TD, phylogenetic-PD, and functional-FD dimensions), 2) drives changes in bird community composition and enables the detection of indicator species of such impact, and 3) leads to changes in bird functional traits linked to reproduction, resource acquisition, and survival. We found that urban intensification has a direct impact on the bird community, reducing all three types of diversity. Communities in areas of greater urban intensity are represented by fewer species, and these species are PD and FD less distinct. In addition, we detected at least ten species of areas of lower urban intensity that proved to be more sensitive to urban intensification. With regard to bird traits, we found no significant responses from reproductive, habitat use and feeding variables. Body weight and tail length were the only variables with significant results, with higher urbanization intensity areas selecting for species with lower weights and longer tails. Given the global biodiversity loss we are observing, this information can guide urban managers and planners in designing urban landscapes to maintain biodiversity in cities.

Phylogenetic Diversity Indices from an Affine and Projective Viewpoint

Phylogenetic diversity indices are commonly used to rank the elements in a collection of species or populations for conservation purposes. The derivation of these indices is typically based on some quantitative description of the evolutionary history of the species in question, which is often given in terms of a phylogenetic tree. Both rooted and unrooted phylogenetic trees can be employed, and there are close connections between the indices that are derived in these two different ways. In this paper, we introduce more general phylogenetic diversity indices that can be derived from collections of subsets (clusters) and collections of bipartitions (splits) of the given set of species. Such indices could be useful, for example, in case there is some uncertainty in the topology of the tree being used to derive a phylogenetic diversity index. As well as characterizing some of the indices that we introduce in terms of their special properties, we provide a link between cluster-based and split-based phylogenetic diversity indices that uses a discrete analogue of the classical link between affine and projective geometry. This provides a unified framework for many of the various phylogenetic diversity indices used in the literature based on rooted and unrooted phylogenetic trees, generalizations and new proofs for previous results concerning tree-based indices, and a way to define some new phylogenetic diversity indices that naturally arise as affine or projective variants of each other or as generalizations of tree-based indices.

Disentangling the drivers of urban bird diversity in the non-breeding season: A general synthesis

Current knowledge about the impacts of urbanisation on bird assemblages is based on evidence from studies partly or wholly undertaken in the breeding season. In comparison, the non-breeding season remains little studied, despite the fact that winter conditions at higher latitudes are changing more rapidly than other seasons. During the non-breeding season, cities may attract or retain bird species because they offer milder conditions or better feeding opportunities than surrounding habitats. However, the range of climatic, ecological and anthropogenic mechanisms shaping different facets of urban bird diversity in the non-breeding season are poorly understood. We explored these mechanisms using structural equation modelling to assess how urbanisation affects the taxonomic, phylogenetic and functional diversity of avian assemblages sampled worldwide in the non-breeding season. We found that minimum temperature, elevation, urban area and city age played a critical role in determining taxonomic diversity while a range of factors-including productivity, precipitation, elevation, distance to coasts and rivers, socio-economic (as a proxy of human facilitation) and road density-each contributed to patterns of phylogenetic and functional diversity. The structure and function of urban bird assemblages appear to be predominantly shaped by temperature, productivity and city age, with effects of these factors differing across seasons. Our results underline the importance of considering multiple hypotheses, including seasonal effects, when evaluating the impacts of urbanisation on biodiversity.

When indices disagree: facing conceptual and practical challenges

Hypothesis testing requires meaningful ways to quantify biological phenomena and account for alternative mechanisms that could explain the same pattern. Researchers combine experiments, statistics, and indices to account for these confounding mechanisms. Key concepts in ecology and evolution, such as niche breadth (NB) or fitness, can be represented by several indices, which often provide uncorrelated estimates. Is this because the indices use different types of noisy data or because the targeted phenomenon is complex and multidimensional? We discuss implications of these scenarios and propose five steps to aid researchers in identifying and combining indices, experiments, and statistics. Building on prior efforts to construct databases of hypotheses and indices and document assumptions, these steps help provide a formal strategy to reduce self-confirmatory bias.

Discriminating woody species assemblages from National Forest Inventory data based on phylogeny in Georgia

Classifications of forest vegetation types and characterization of related species assemblages are important analytical tools for mapping and diversity monitoring of forest communities. The discrimination of forest communities is often based on β-diversity, which can be quantified via numerous indices to derive compositional dissimilarity between samples. This study aims to evaluate the applicability of unsupervised classification for National Forest Inventory data from Georgia by comparing two cluster hierarchies. We calculated the mean basal area per hectare for each woody species across 1059 plot observations and quantified interspecies distances for all 87 species. Following an unspuervised cluster analysis, we compared the results derived from the species-neutral dissimilarity (Bray-Curtis) with those based on the Discriminating Avalanche dissimilarity, which incorporates interspecies phylogenetic variation. Incorporating genetic variation in the dissimilarity quantification resulted in a more nuanced discrimination of woody species assemblages and increased cluster coherence. Favorable statistics include the total number of clusters (23 vs. 20), mean distance within clusters (0.773 vs. 0.343), and within sum of squares (344.13 vs. 112.92). Clusters derived from dissimilarities that account for genetic variation showed a more robust alignment with biogeographical units, such as elevation and known habitats. We demonstrate that the applicability of unsupervised classification of species assemblages to large-scale forest inventory data strongly depends on the underlying quantification of dissimilarity. Our results indicate that by incorporating phylogenetic variation, a more precise classification aligned with biogeographic units is attained. This supports the concept that the genetic signal of species assemblages reflects biogeographical patterns and facilitates more precise analyses for mapping, monitoring, and management of forest diversity.

Ecological drivers of taxonomic, functional, and phylogenetic diversity of bryophytes in an oceanic island

Montane oceanic islands possess unique geographic and ecological attributes, rendering them valuable for assessing patterns and drivers of alpha and beta taxonomic, functional, and phylogenetic diversity along elevational gradients. Such comparisons of diversity facets can provide insights into the mechanisms governing community assembly on islands. Herein, we aimed to characterize taxonomic, functional, and phylogenetic bryophyte diversity on Madeira Island within and across areas at varying elevations. We also assessed how these diversity facets for the alpha and beta components relate to ecological and anthropogenic factors. We estimated and compared alpha and beta taxonomic, functional, and phylogenetic diversity using 80 plots of 0.5 m × 0.5 m across the whole elevational gradient of the island. We compiled trait databases and supplemented them with our own observations. Phylogenetic information was sourced from the Moss and Liverwort Tree of Life. To assess the impact of ecological and anthropogenic factors on the three facets, we applied linear mixed-effects models and generalized dissimilarity models to alpha- and beta-diversity matrices, respectively. All facets of diversity exhibited strong correlations within both mosses and liverworts, indicating a substantial congruence when alpha and beta are analyzed separately. The bryophyte groups categorized by the growth form demonstrated contrasting patterns, aligning with their distinctive ecological requirements. While a mid-elevation peak emerged as a common pattern across the three facets of alpha diversity, beta diversity often displayed the opposite trend. Although the relative influence of environmental factors varied depending on the diversity facet and bryophyte grouping considered, we found that alpha and beta diversity of bryophytes are more influenced by climatic factors and the predominant type of vegetation than by anthropogenic factors. In the current context of global change, these results should be interpreted with caution, but they point to the resilience of bryophytes to survive in relatively well-preserved natural microhabitats within anthropogenic landscapes. In this study on Madeira Island, we investigated patterns and drivers of alpha and beta taxonomic, functional, and phylogenetic diversity along elevational gradients. We found that alpha and beta diversity of bryophytes are more strongly influenced by climatic factors and the predominant type of vegetation than by anthropogenic factors.

Wildfire Impacts on Molecular Changes of Dissolved Organic Matter during Its Passage through Soil

The frequency and intensity of global wildfires are escalating, leading to an increase in derived pyrogenic dissolved organic matter (pyDOM), which potentially influences the riverine carbon reservoir and poses risks to drinking water safety. However, changes in pyDOM properties as it traverses through soil to water bodies are highly understudied due to the challenges of simulating such processes under laboratory conditions. In this study, we extracted soil DOM along hillslope gradients and soil depths in both burned and unburned catchments post wildfire. Using high-resolution mass spectrometry and a substrate-explicit model, we observed significant increases in the relative abundance of condensed aromatics (ConAC) and tannins in wildfire-affected soil DOM. Wildfire-affected soil DOM also displayed a broader spectrum of molecular and thermodynamic properties, indicative of its diverse composition and reactivity. Furthermore, as the fire-induced weakening of topsoil microbial reprocessing abilities hindered the transformation of plant-derived DOM, the relative abundance of lignin-like compounds increased with soil depth in the fire regions. Meanwhile, the distribution of shared molecular formulas along the hillslope gradient (from shoulder to toeslope) exhibited analogous patterns in both burned and unburned catchments. Although there was an increased prevalence of ConAC and tannin in the burned catchments, the relative abundance of these fractions diminished along the hillslope in all three catchments. Based on the substrate-explicit model, the biodegradability exhibited by wildfire-affected DOM fractions offers the possibility of its conversion along hillslopes. Our findings reveal the spatial distribution of DOM properties after a wildfire, facilitating accurate evaluation of dissolved organic carbon composition involved in the watershed-scale carbon cycle.

Identifying global conservation priorities for terrestrial vertebrates based on multiple dimensions of biodiversity

The Kunming-Montreal Global Biodiversity Framework of the Convention on Biological Diversity calls for an expansion of the current protected areas (PAs) to cover at least 30% of global land and water areas by 2030 (i.e., the 30×30 target). Efficient spatial planning for PA expansion is an urgent need for global conservation practice. A spatial prioritization framework considering multiple dimensions of biodiversity is critical for improving the efficiency of the spatial planning of PAs, yet it remains a challenge. We developed an index for the identification of priority areas based on functionally rare, evolutionarily distinct, and globally endangered species (FREDGE) and applied it to 21,536 terrestrial vertebrates. We determined species distributions, conservation status (global endangerment), molecular phylogenies (evolutionary distinctiveness), and life-history traits (functional rarity). Madagascar, Central America, and the Andes were of high priority for the conservation of multiple dimensions of terrestrial vertebrate biodiversity. However, 68.8% of grid cells in these priority areas had <17% of their area covered by PAs, and these priority areas were under intense anthropogenic and climate change threats. These results highlight the difficulties of conserving multiple dimensions of biodiversity. Our global analyses of the geographical patterns of multiple dimensions of terrestrial vertebrate biodiversity demonstrate the insufficiency of the conservation of different biodiversity dimensions, and our index, based on multiple dimensions of biodiversity, provides a useful tool for guiding future spatial prioritization of PA expansion to achieve the 30×30 target under serious pressures.

Phylogenetic Biodiversity Metrics Should Account for Both Accumulation and Attrition of Evolutionary Heritage

Phylogenetic metrics are essential tools used in the study of ecology, evolution and conservation. Phylogenetic diversity (PD) in particular is one of the most prominent measures of biodiversity and is based on the idea that biological features accumulate along the edges of phylogenetic trees that are summed. We argue that PD and many other phylogenetic biodiversity metrics fail to capture an essential process that we term attrition. Attrition is the gradual loss of features through causes other than extinction. Here we introduce "EvoHeritage", a generalization of PD that is founded on the joint processes of accumulation and attrition of features. We argue that while PD measures evolutionary history, EvoHeritage is required to capture a more pertinent subset of evolutionary history including only components that have survived attrition. We show that EvoHeritage is not the same as PD on a tree with scaled edges; instead, accumulation and attrition interact in a more complex non-monophyletic way that cannot be captured by edge lengths alone. This leads us to speculate that the one-dimensional edge lengths of classic trees may be insufficiently flexible to capture the nuances of evolutionary processes. We derive a measure of EvoHeritage and show that it elegantly reproduces species richness and PD at opposite ends of a continuum based on the intensity of attrition. We demonstrate the utility of EvoHeritage in ecology as a predictor of community productivity compared with species richness and PD. We also show how EvoHeritage can quantify living fossils and resolve their associated controversy. We suggest how the existing calculus of PD-based metrics and other phylogenetic biodiversity metrics can and should be recast in terms of EvoHeritage accumulation and attrition.

Gut metagenomic analysis of gastric cancer patients reveals Akkermansia, Gammaproteobacteria, and Veillonella microbiota as potential non-invasive biomarkers

The goal of the study was to investigate the changes in the gut microbiota during the advancement of gastric cancer (GC) and identify pertinent taxa associated with the disease. We used a public fecal amplicon gastric cancer dataset from the Sequence Retrieval Archive (SRA), of patients with GC, gastritis, and healthy individuals. We did sequence pre-processing, including quality filtering of the sequences. Then, we performed a diversity analysis, evaluating α- and β-diversity. Next, taxonomic composition analysis was performed and the relative abundances of different taxa at the phylum and genus levels were compared between GC, gastritis, and healthy controls. The obtained results were subsequently subjected to statistical validation. To conclude, metagenomic function prediction was carried out, followed by correlation analysis between the microbiota and KEGG pathways. α analysis revealed a significant difference between male and female categories, while β analysis demonstrated significant distinctions between GC, gastritis, and healthy controls, as well as between sexes within the GC and gastritis groups. The statistically confirmed taxonomic composition analysis highlighted the presence of the microbes Bacteroides and Veillonella. Furthermore, through metagenomic prediction analysis and correlation analysis with pathways, three taxa, namely Akkermansia, Gammaproteobacteria, and Veillonella, were identified as potential biomarkers for GC. Additionally, this study reports, for the first time, the presence of two bacteria, Desulfobacteriota and Synergistota, in GC, necessitating further investigation. Overall, this research sheds light on the potential involvement of gut microbiota in GC pathophysiology; however, additional studies are warranted to explore its functional significance.

Taxonomic, phylogenetic, and functional diversity of mollusk death assemblages in coral reef and seagrass sediments from two shallow gulfs in Western Cuban Archipelago

Mollusk death assemblages are formed by shell remnants deposited in the surficial mixed layer of the seabed. Diversity patterns in tropical marine habitats still are understudied; therefore, we aimed to investigate the taxonomic, phylogenetic, and functional diversity of mollusk death assemblages at regional and local scales in coral reef sands and seagrass meadows. We collected sediment samples at 11 sites within two shallow gulfs in the Northwestern Caribbean Sea and Southeastern Gulf of Mexico. All the shells were counted and identified to species level and classified into biological traits. We identified 7113 individuals belonging to 393 species (290 gastropods, 94 bivalves, and nine scaphopods). Diversity and assemblage structure showed many similarities between gulfs given their geological and biogeographical commonalities. Reef sands had higher richness than seagrasses likely because of a more favorable balance productivity-disturbance. Reef sands were dominated by epifaunal herbivores likely feeding on microphytobenthos and bysally attached bivalves adapted to intense hydrodynamic regime. In seagrass meadows, suspension feeders dominated in exposed sites and chemosynthetic infaunal bivalves dominated where oxygen replenishment was limited. Time averaging of death assemblages was likely in the order of 100 years, with stronger effects in reef sands compared to seagrass meadows. Our research provides evidence of the high taxonomic, phylogenetic, and functional diversity of mollusk death assemblages in tropical coastal sediments as result of the influence of scale-related processes and habitat type. Our study highlights the convenience of including phylogenetic and functional traits, as well as dead shells, for a more complete assessment of mollusk biodiversity.

Relationships between functional alpha and beta diversities of flea parasites and their small mammalian hosts

We studied the relationships between functional alpha and beta diversities of fleas and their small mammalian hosts in 4 biogeographic realms (the Afrotropics, the Nearctic, the Neotropics and the Palearctic), considering 3 components of alpha diversity (functional richness, divergence and regularity). We asked whether (a) flea alpha and beta diversities are driven by host alpha and beta diversities; (b) the variation in the off-host environment affects variation in flea alpha and beta diversities; and (c) the pattern of the relationship between flea and host alpha or beta diversities differs between geographic realms. We analysed alpha diversity using modified phylogenetic generalized least squares and beta diversity using modified phylogenetic generalized dissimilarity modelling. In all realms, flea functional richness and regularity increased with an increase in host functional richness and regularity, respectively, whereas flea functional divergence correlated positively with host functional divergence in the Nearctic only. Environmental effects on the components of flea alpha diversity were found only in the Holarctic realms. Host functional beta diversity was invariantly the best predictor of flea functional beta diversity in all realms, whereas the effects of environmental variables on flea functional beta diversity were much weaker and differed between realms. We conclude that flea functional diversity is mostly driven by host functional diversity, whereas the environmental effects on flea functional diversity vary (a) geographically and (b) between components of functional alpha diversity.

Global analysis of Poales diversification - parallel evolution in space and time into open and closed habitats

Poales are one of the most species-rich, ecologically and economically important orders of plants and often characterise open habitats, enabled by unique suites of traits. We test six hypotheses regarding the evolution and assembly of Poales in open and closed habitats throughout the world, and examine whether diversification patterns demonstrate parallel evolution. We sampled 42% of Poales species and obtained taxonomic and biogeographic data from the World Checklist of Vascular Plants database, which was combined with open/closed habitat data scored by taxonomic experts. A dated supertree of Poales was constructed. We integrated spatial phylogenetics with regionalisation analyses, historical biogeography and ancestral state estimations. Diversification in Poales and assembly of open and closed habitats result from dynamic evolutionary processes that vary across lineages, time and space, most prominently in tropical and southern latitudes. Our results reveal parallel and recurrent patterns of habitat and trait transitions in the species-rich families Poaceae and Cyperaceae. Smaller families display unique and often divergent evolutionary trajectories. The Poales have achieved global dominance via parallel evolution in open habitats, with notable, spatially and phylogenetically restricted divergences into strictly closed habitats.

Spatial heterogeneity of neo- and paleo-endemism for plants in Madagascar

Madagascar is a biogeographically unique island with a remarkably high level of endemism. However, endemic taxa in Madagascar are massively threatened due to unprecedented pressures from anthropogenic habitat modification and climate change. A comprehensive phylogeny-based biodiversity evaluation of the island remains lacking. Here, we identify hotspots of taxonomic and phylogenetic plant diversity and neo- and paleo-endemism by generating a novel dated tree of life for the island. The tree is based on unprecedented sampling of 3,950 species (33% of the total known species) and 1,621 genera (93% of the total known genera and 69% of endemic genera) of Malagasy vascular plants. We find that island-endemic genera are concentrated in multiple lineages combining high taxonomic and phylogenetic diversity. Integrating phylogenetic and geographic distribution data, our results reveal that taxon richness and endemism are concentrated in the northern, eastern, and southeastern humid forests. Paleo-endemism centers are concentrated in humid eastern and central regions, whereas neo-endemism centers are concentrated in the dry and spiny forests in western and southern Madagascar. Our statistical analysis of endemic genera in each vegetation region supports a higher proportion of ancient endemic genera in the east but a higher proportion of recent endemic genera in the south and west. Overlaying centers of phylogenetic endemism with protected areas, we identify conservation gaps concentrated in western and southern Madagascar. These gaps should be incorporated into conservation strategies to aid the protection of multiple facets of biodiversity and their benefits to the Malagasy people.

Diversity of an Odonata assemblage from a tropical dry forest in San Buenaventura, Jalisco, Mexico (Insecta, Odonata)

Background

The patterns of richness, diversity, and abundance of an odonate assemblage from San Buenaventura, Jalisco are presented here. A total of 1087 specimens from seven families, 35 genera and 66 species were obtained through monthly samplings of five days each during a period of one year. Libellulidae was the most diverse family (28 species), followed by Coenagrionidae (21), Gomphidae (7), Aeshnidae (6), Calopterygidae (2), Lestidae (1) and Platystictidae (1). Argia was the most speciose genus. The highest species richness and Shannon diversity were found during August and September, whereas the highest abundance was observed in June and the highest Simpson diversity was recorded in September - all of which were associated with the rainy season. The highest values of phylogenetic diversity were found from June to October. The different diversity facets of this assemblage were positively correlated with precipitation and minimum temperature, whereas maximum temperature showed no influence. In addition, we found that this odonate diversity was higher than most Mexican localities with tropical dry forest (TDF) studied.

New information

We continue our efforts to describe the patterns of richness, diversity and abundance of some insect groups associated with the tropical dry forest ecosystem in Mexico, following a latitudinal gradient of the distribution of this ecosystem in the country. Our emphasis here was to evaluate the spatial and temporal patterns of richness and diversity of an Odonata assemblage from Jalisco, Mexico.

Dense city centers support less evolutionary unique bird communities than sparser urban areas

Urbanization alters avian communities, generally lowering the number of species and contemporaneously increasing their functional relatedness, leading to biotic homogenization. Urbanization can also negatively affect the phylogenetic diversity of species assemblages, potentially decreasing their evolutionary distinctiveness. We compare species assemblages in a gradient of building density in seventeen European cities to test whether the evolutionary distinctiveness of communities is shaped by the degree of urbanization. We found a significant decline in the evolutionary uniqueness of avian communities in highly dense urban areas, compared to low and medium-dense areas. Overall, communities from dense city centers supported one million years of evolutionary history less than communities from low-dense urban areas. Such evolutionary homogenization was due to a filtering process of the most evolutionarily unique birds. Metrics related to evolutionary uniqueness have to play a role when assessing the effects of urbanization and can be used to identify local conservation priorities.

Dental complexity and diet in amniotes: A meta-analysis

Tooth morphology is among the most well-studied indicators of ecology. For decades, researchers have examined the gross morphology and wear patterns of teeth as indicators of diet, and recent advances in scanning and computer analysis have allowed the development of new and more quantitative measures of tooth morphology. One of the most popular of these new methods is orientation patch count (OPC). OPC, a measure of surface complexity, was originally developed to distinguish the more complex tooth crowns of herbivores from the less complex tooth crowns of faunivores. OPC and a similar method derived from it, orientation patch count rotated (OPCR), have become commonplace in analyses of both modern and fossil amniote dietary ecology. The widespread use of these techniques makes it possible to now re-assess the utility of OPC and OPCR. Here, we undertake a comprehensive review of OPC(R) and diet and perform a meta-analysis to determine the overall difference in complexity between herbivores and faunivores. We find that the relationship between faunivore and herbivore OPC or OPCR values differs substantially across studies, and although some support the initial assessment of greater complexity in herbivores, others do not. Our meta-analysis does not support an overall pattern of greater complexity in herbivores than faunivores across terrestrial amniotes. It appears that the relationship of OPC or OPCR to diet is taxon-specific and dependent on the type of faunivory of the group in question, with insectivores often having values similar to herbivores. We suggest extreme caution in comparing OPC and OPCR values across studies and offer suggestions for how OPCR can constructively be used in future research.

Optimal trade-off control in machine learning-based library design, with application to adeno-associated virus (AAV) for gene therapy

Adeno-associated viruses (AAVs) hold tremendous promise as delivery vectors for gene therapies. AAVs have been successfully engineered-for instance, for more efficient and/or cell-specific delivery to numerous tissues-by creating large, diverse starting libraries and selecting for desired properties. However, these starting libraries often contain a high proportion of variants unable to assemble or package their genomes, a prerequisite for any gene delivery goal. Here, we present and showcase a machine learning (ML) method for designing AAV peptide insertion libraries that achieve fivefold higher packaging fitness than the standard NNK library with negligible reduction in diversity. To demonstrate our ML-designed library's utility for downstream engineering goals, we show that it yields approximately 10-fold more successful variants than the NNK library after selection for infection of human brain tissue, leading to a promising glial-specific variant. Moreover, our design approach can be applied to other types of libraries for AAV and beyond.

Phylogenetic diversity and community wide-trait means offer different insights into mechanisms regulating aboveground carbon storage

Both attributes of functional traits and phylogenetic diversity influence ecosystem functions, but which of these factors is most important is still poorly understood in natural systems. Using data from West African forests and tree savannas, we analyse how (i) phylogenetic diversity complements attributes of functional traits in explaining aboveground carbon (AGC); (ii) phylogenetic diversity relates with attributes of functional traits along gradients of phylogenetic signal; and (iii) pathways between phylogenetic diversity and attributes of functional traits relate AGC to soil and climate. Phylogenetic diversity was measured as standardised effect size of Mean Pairwise Distance (sesMPD) and Mean Nearest Taxon Distance (sesMNTD). Functional dispersion (FDis) and community weighted mean (CWM) were calculated for four traits related to leaf economics spectrum and plant life-history. Functional traits-based models explained 11 % of AGC variability. With two out of the four traits being phylogenetically conserved, incorporating phylogenetic diversity in the models increased the explained variance in AGC by 15 %. The slope of phylogenetic diversity-trait relationship was more responsive to trait conservatism for FDis than CWM. AGC was positively influenced by sesMPD and CWM of plant maximum height. In turn, CWM of plant maximum height increased with higher soil nitrogen and climate moisture, whereas sesMPD was negatively related with climate moisture. Although FDis was positively associated with sesMPD, it was not as important as sesMPD and CWM of plant maximum height in influencing and relating AGC to soil nitrogen and climate moisture. Our results suggest that phylogenetic diversity is important for AGC but does not fully reflect the functional mechanisms pertaining to community-wide trait means. The study also demonstrates the role of environment in regulating AGC, which operates through differences in community fitness driven by tall plant stature, and evolutionary processes whereby closely related species are maintained in less arid environments.

Diversity across organisational scale emerges through dispersal ability and speciation dynamics in tropical fish

BACKGROUND

Biodiversity exists at different levels of organisation: e.g. genetic, individual, population, species, and community. These levels of organisation all exist within the same system, with diversity patterns emerging across organisational scales through several key processes. Despite this inherent interconnectivity, observational studies reveal that diversity patterns across levels are not consistent and the underlying mechanisms for variable continuity in diversity across levels remain elusive. To investigate these mechanisms, we apply a spatially explicit simulation model to simulate the global diversification of tropical reef fishes at both the population and species levels through emergent population-level processes.

RESULTS

We find significant relationships between the population and species levels of diversity which vary depending on both the measure of diversity and the spatial partitioning considered. In turn, these population-species relationships are driven by modelled biological trait parameters, especially the divergence threshold at which populations speciate.

CONCLUSIONS

To explain variation in multi-level diversity patterns, we propose a simple, yet novel, population-to-species diversity partitioning mechanism through speciation which disrupts continuous diversity patterns across organisational levels. We expect that in real-world systems this mechanism is driven by the molecular dynamics that determine genetic incompatibility, and therefore reproductive isolation between individuals. We put forward a framework in which the mechanisms underlying patterns of diversity across organisational levels are universal, and through this show how variable patterns of diversity can emerge through organisational scale.

A new universal system of tree shape indices

The comparison and categorization of tree diagrams is fundamental to large parts of biology, linguistics, computer science, and other fields, yet the indices currently applied to describing tree shape have important flaws that complicate their interpretation and limit their scope. Here we introduce a new system of indices with no such shortcomings. Our indices account for node sizes and branch lengths and are robust to small changes in either attribute. Unlike currently popular phylogenetic diversity, phylogenetic entropy, and tree balance indices, our definitions assign interpretable values to all rooted trees and enable meaningful comparison of any pair of trees. Our self-consistent definitions further unite measures of diversity, richness, balance, symmetry, effective height, effective outdegree, and effective branch count in a coherent system, and we derive numerous simple relationships between these indices. The main practical advantages of our indices are in 1) quantifying diversity in non-ultrametric trees; 2) assessing the balance of trees that have non-uniform branch lengths or node sizes; 3) comparing the balance of trees with different leaf counts or outdegrees; 4) obtaining a coherent, generic, multidimensional quantification of tree shape that is robust to sampling error and inferential error. We illustrate these features by comparing the shapes of trees representing the evolution of HIV and of Uralic languages, and trees generated by computational models of tumour evolution. Given the ubiquity of tree structures, we identify a wide range of applications across diverse domains.

Functional diversity of sharks and rays is highly vulnerable and supported by unique species and locations worldwide

Elasmobranchs (sharks, rays and skates) are among the most threatened marine vertebrates, yet their global functional diversity remains largely unknown. Here, we use a trait dataset of >1000 species to assess elasmobranch functional diversity and compare it against other previously studied biodiversity facets (taxonomic and phylogenetic), to identify species- and spatial- conservation priorities. We show that threatened species encompass the full extent of functional space and disproportionately include functionally distinct species. Applying the conservation metric FUSE (Functionally Unique, Specialised, and Endangered) reveals that most top-ranking species differ from the top Evolutionarily Distinct and Globally Endangered (EDGE) list. Spatial analyses further show that elasmobranch functional richness is concentrated along continental shelves and around oceanic islands, with 18 distinguishable hotspots. These hotspots only marginally overlap with those of other biodiversity facets, reflecting a distinct spatial fingerprint of functional diversity. Elasmobranch biodiversity facets converge with fishing pressure along the coast of China, which emerges as a critical frontier in conservation. Meanwhile, several components of elasmobranch functional diversity fall in high seas and/or outside the global network of marine protected areas. Overall, our results highlight acute vulnerability of the world's elasmobranchs' functional diversity and reveal global priorities for elasmobranch functional biodiversity previously overlooked.

Climate change and land use threaten global hotspots of phylogenetic endemism for trees

Across the globe, tree species are under high anthropogenic pressure. Risks of extinction are notably more severe for species with restricted ranges and distinct evolutionary histories. Here, we use a global dataset covering 41,835 species (65.1% of known tree species) to assess the spatial pattern of tree species' phylogenetic endemism, its macroecological drivers, and how future pressures may affect the conservation status of the identified hotspots. We found that low-to-mid latitudes host most endemism hotspots, with current climate being the strongest driver, and climatic stability across thousands to millions of years back in time as a major co-determinant. These hotspots are mostly located outside of protected areas and face relatively high land-use change and future climate change pressure. Our study highlights the risk from climate change for tree diversity and the necessity to strengthen conservation and restoration actions in global hotspots of phylogenetic endemism for trees to avoid major future losses of tree diversity.

A latitudinal gradient in Darwin's naturalization conundrum at the global scale for flowering plants

Darwin's naturalization conundrum describes two seemingly contradictory hypotheses regarding whether alien species closely or distantly related to native species should be more likely to naturalize in regional floras. Both expectations have accumulated empirical support, and whether such apparent inconsistency can be reconciled at the global scale is unclear. Here, using 219,520 native and 9,531 naturalized alien plant species across 487 globally distributed regions, we found a latitudinal gradient in Darwin's naturalization conundrum. Naturalized alien plant species are more closely related to native species at higher latitudes than they are at lower latitudes, indicating a greater influence of preadaptation in harsher climates. Human landscape modification resulted in even steeper latitudinal clines by selecting aliens distantly related to natives in warmer and drier regions. Our results demonstrate that joint consideration of climatic and anthropogenic conditions is critical to reconciling Darwin's naturalization conundrum.

Developing a New Phylogeny-Driven Random Forest Model for Functional Metagenomics

Metagenomics is an unobtrusive science linking microbial genes to biological functions or environmental states. Classifying microbial genes into their functional repertoire is an important task in the downstream analysis of Metagenomic studies. The task involves Machine Learning (ML) based supervised methods to achieve good classification performance. Random Forest (RF) has been applied rigorously to microbial gene abundance profiles, mapping them to functional phenotypes. The current research targets tuning RF by the evolutionary ancestry of microbial phylogeny, developing a Phylogeny-RF model for functional classification of metagenomes. This method facilitates capturing the effects of phylogenetic relatedness in an ML classifier itself rather than just applying a supervised classifier over the raw abundances of microbial genes. The idea is rooted in the fact that closely related microbes by phylogeny are highly correlated and tend to have similar genetic and phenotypic traits. Such microbes behave similarly; and hence tend to be selected together, or one of these could be dropped from the analysis, to improve the ML process. The proposed Phylogeny-RF algorithm has been compared with state-of-the-art classification methods including RF and the phylogeny-aware methods of MetaPhyl and PhILR, using three real-world 16S rRNA metagenomic datasets. It has been observed that the proposed method not only achieved significantly better performance than the traditional RF model but also performed better than the other phylogeny-driven benchmarks (p < 0.05). For example, Phylogeny-RF attained a highest AUC of 0.949 and Kappa of 0.891 over soil microbiomes in comparison to other benchmarks.

Components of functional diversity revisited: A new classification and its theoretical and practical implications

Functional diversity is regarded as a key concept for understanding the link between ecosystem function and biodiversity. The different and ecologically well-defined aspects of the concept are reflected by the so-called functional components, for example, functional richness and divergence. Many authors proposed that components be distinguished according to the multivariate technique on which they rely, but more recent studies suggest that several multivariate techniques, providing different functional representations (such as dendrograms and ordinations) of the community can in fact express the same functional component. Here, we review the relevant literature and find that (1) general ecological acceptance of the field is hampered by ambiguous terminology and (2) our understanding of the role of multivariate techniques in defining components is unclear. To address these issues, we provide new definitions for the three basic functional diversity components namely functional richness, functional divergence and functional regularity. In addition, we present a classification of presence-/absence-based approaches suitable for quantifying these components. We focus exclusively on the binary case for its relative simplicity. We find illogical, as well as logical but unused combinations of components and representations; and reveal that components can be quantified almost independently from the functional representation of the community. Finally, theoretical and practical implications of the new classification are discussed.

Native diversity buffers against severity of non-native tree invasions

Determining the drivers of non-native plant invasions is critical for managing native ecosystems and limiting the spread of invasive species1,2. Tree invasions in particular have been relatively overlooked, even though they have the potential to transform ecosystems and economies3,4. Here, leveraging global tree databases5-7, we explore how the phylogenetic and functional diversity of native tree communities, human pressure and the environment influence the establishment of non-native tree species and the subsequent invasion severity. We find that anthropogenic factors are key to predicting whether a location is invaded, but that invasion severity is underpinned by native diversity, with higher diversity predicting lower invasion severity. Temperature and precipitation emerge as strong predictors of invasion strategy, with non-native species invading successfully when they are similar to the native community in cold or dry extremes. Yet, despite the influence of these ecological forces in determining invasion strategy, we find evidence that these patterns can be obscured by human activity, with lower ecological signal in areas with higher proximity to shipping ports. Our global perspective of non-native tree invasion highlights that human drivers influence non-native tree presence, and that native phylogenetic and functional diversity have a critical role in the establishment and spread of subsequent invasions.

Evolutionary history shapes grassland productivity through opposing effects on complementarity and selection

Phylogenetic diversity (PD), the evolutionary history of the organisms comprising a community, is increasingly recognized as an important driver of ecosystem function. However, biodiversity-ecosystem function experiments have rarely included PD as an a priori treatment. Thus, PD's effects in existing experiments are often confounded by covarying differences in species richness and functional trait diversity (FD). Here we report an experimental demonstration of strong PD effects on grassland primary productivity that are independent of FD, which was separately manipulated, and species richness, which was planted uniformly high to mimic diverse natural grasslands. Partitioning diversity effects demonstrated that higher PD increased complementarity (niche partitioning and/or facilitation) but lowered selection effects (probability of sampling highly productive species). Specifically, for every 5% increase in PD, complementarity increased by 26% on average (±8% SE), while selection effects decreased more modestly (8 ± 16%). PD also shaped productivity through clade-level effects on functional traits, that is, trait values associated with particular plant families. This clade effect was most pronounced in the Asteraceae (sunflower family), which, in tallgrass prairies, generally comprises tall, high-biomass species with low phylogenetic distinctiveness. FD also reduced selection effects but did not alter complementarity. Our results show that PD, independent of richness and FD, mediates ecosystem function through contrasting effects on complementarity and selection. This adds to growing evidence that consideration of phylogenetic dimensions of biodiversity can advance ecological understanding and inform conservation and restoration.

Phylogenetically-conserved candidate genes unify biodiversity-ecosystem function relationships and eco-evolutionary dynamics across biological scales

The intra- and interspecific facets of biodiversity have traditionally been analysed separately, limiting our understanding of how evolution has shaped biodiversity, how biodiversity (as a whole) alters ecological dynamics and hence eco-evolutionary feedbacks at the community scale. Here, we propose using candidate genes phylogenetically-conserved across species and sustaining functional traits as an inclusive biodiversity unit transcending the intra- and interspecific boundaries. This framework merges knowledge from functional genomics and functional ecology, and we first provide guidelines and a concrete example for identifying phylogenetically-conserved candidate genes (PCCGs) within communities and for measuring biodiversity from PCCGs. We then explain how biodiversity measured at PCCGs can be linked to ecosystem functions, which unifies recent observations that both intra- and interspecific biodiversity are important for ecosystem functions. We then highlight the eco-evolutionary processes shaping PCCG diversity patterns and argue that their respective role can be inferred from concepts derived from population genetics. Finally, we explain how PCCGs may shift the field of eco-evolutionary dynamics from a focal-species approach to a more realistic focal-community approach. This framework provides a novel perspective to investigate the global ecosystem consequences of diversity loss across biological scales, and how these ecological changes further alter biodiversity evolution.

Resolving the Pathogenesis of Anaplastic Wilms Tumors through Spatial Mapping of Cancer Cell Evolution

PURPOSE

While patients with intermediate-risk (IR) Wilms tumors now have an overall survival (OS) rate of almost 90%, those affected by high-stage tumors with diffuse anaplasia have an OS of only around 50%. We here identify key events in the pathogenesis of diffuse anaplasia by mapping cancer cell evolution over anatomic space in Wilms tumors.

EXPERIMENTAL DESIGN

We spatially mapped subclonal landscapes in a retrospective cohort of 20 Wilms tumors using high-resolution copy-number profiling and TP53 mutation analysis followed by clonal deconvolution and phylogenetic reconstruction. Tumor whole-mount sections (WMS) were utilized to characterize the distribution of subclones across anatomically distinct tumor compartments.

RESULTS

Compared with non-diffuse anaplasia Wilms tumors, tumors with diffuse anaplasia showed a significantly higher number of genetically distinct tumor cell subpopulations and more complex phylogenetic trees, including high levels of phylogenetic species richness, divergence, and irregularity. All regions with classical anaplasia showed TP53 alterations. TP53 mutations were frequently followed by saltatory evolution and parallel loss of the remaining wild-type (WT) allele in different regions. Morphologic features of anaplasia increased with copy-number aberration (CNA) burden and regressive features. Compartments demarcated by fibrous septae or necrosis/regression were frequently (73%) associated with the emergence of new clonal CNAs, although clonal sweeps were rare within these compartments.

CONCLUSIONS

Wilms tumors with diffuse anaplasia display significantly more complex phylogenies compared with non-diffuse anaplasia Wilms tumors, including features of saltatory and parallel evolution. The subclonal landscape of individual tumors was constrained by anatomic compartments, which should be considered when sampling tissue for precision diagnostics.

Spaces of Phylogenetic Diversity Indices: Combinatorial and Geometric Properties

Biodiversity is a concept most naturally quantified and measured across sets of species. However, for some applications, such as prioritising species for conservation efforts, a species-by-species approach is desirable. Phylogenetic diversity indices are functions that apportion the total biodiversity value of a set of species across its constituent members. As such, they aim to measure each species' individual contribution to, and embodiment of, the diversity present in that set. However, no clear definition exists that encompasses the diversity indices in current use. This paper presents conditions that define diversity indices arising from the phylogenetic diversity measure on rooted phylogenetic trees. In this context, the diversity index 'score' given to a species represents a measure of its unique and shared evolutionary history as displayed in the underlying phylogenetic tree. Our definition generalises the diversity index notion beyond the popular Fair Proportion and Equal-Splits indices. These particular indices may now be seen as two points in a convex space of possible diversity indices, for which the boundary conditions are determined by the underlying shape of each phylogenetic tree. We calculated the dimension of the convex space associated with each tree shape and described the extremal points.

Observed and projected functional reorganization of riverine fish assemblages from global change

Climate and land-use/land-cover change ("global change") are restructuring biodiversity, globally. Broadly, environmental conditions are expected to become warmer, potentially drier (particularly in arid regions), and more anthropogenically developed in the future, with spatiotemporally complex effects on ecological communities. We used functional traits to inform Chesapeake Bay Watershed fish responses to future climate and land-use scenarios (2030, 2060, and 2090). We modeled the future habitat suitability of focal species representative of key trait axes (substrate, flow, temperature, reproduction, and trophic) and used functional and phylogenetic metrics to assess variable assemblage responses across physiographic regions and habitat sizes (headwaters through large rivers). Our focal species analysis projected future habitat suitability gains for carnivorous species with preferences for warm water, pool habitats, and fine or vegetated substrates. At the assemblage level, models projected decreasing habitat suitability for cold-water, rheophilic, and lithophilic individuals but increasing suitability for carnivores in the future across all regions. Projected responses of functional and phylogenetic diversity and redundancy differed among regions. Lowland regions were projected to become less functionally and phylogenetically diverse and more redundant while upland regions (and smaller habitat sizes) were projected to become more diverse and less redundant. Next, we assessed how these model-projected assemblage changes 2005-2030 related to observed time-series trends (1999-2016). Halfway through the initial projecting period (2005-2030), we found observed trends broadly followed modeled patterns of increasing proportions of carnivorous and lithophilic individuals in lowland regions but showed opposing patterns for functional and phylogenetic metrics. Leveraging observed and predicted analyses simultaneously helps elucidate the instances and causes of discrepancies between model predictions and ongoing observed changes. Collectively, results highlight the complexity of global change impacts across broad landscapes that likely relate to differences in assemblages' intrinsic sensitivities and external exposure to stressors.

Range restricted old and young lineages show the southern Western Ghats to be both a museum and a cradle of diversity for woody plants

The Western Ghats (WG) mountain chain is a global biodiversity hotspot with high diversity and endemicity of woody plants. The latitudinal breadth of the WG offers an opportunity to determine the evolutionary drivers of latitudinal diversity patterns. We examined the spatial patterns of evolutionary diversity using complementary phylogenetic diversity and endemism measures. To examine if different regions of the WG serve as a museum or cradle of evolutionary diversity, we examined the distribution of 470 species based on distribution modelling and occurrence locations across the entire region. In accordance with the expectation, we found that the southern WG is both a museum and cradle of woody plant evolutionary diversity, as a higher proportion of both old and young evolutionary lineages are restricted to the southern WG. The diversity gradient is likely driven by high geo-climatic stability in the south and phylogenetic niche conservatism for moist and aseasonal sites. This is corroborated by persistent lineage nestedness at almost all evolutionary depths (10-135 million years), and a strong correlation of evolutionary diversity with drought seasonality, precipitation and topographic heterogeneity. Our results highlight the global value of the WG, demonstrating, in particular, the importance of protecting the southern WG-an engine of plant diversification and persistence.

EVI and NDVI as proxies for multifaceted avian diversity in urban areas

Most ecological studies use remote sensing to analyze broad-scale biodiversity patterns, focusing mainly on taxonomic diversity in natural landscapes. One of the most important effects of high levels of urbanization is species loss (i.e., biotic homogenization). Therefore, cost-effective and more efficient methods to monitor biological communities' distribution are essential. This study explores whether the Enhanced Vegetation Index (EVI) and the Normalized Difference Vegetation Index (NDVI) can predict multifaceted avian diversity, urban tolerance, and specialization in urban landscapes. We sampled bird communities among 15 European cities and extracted Landsat 30-meter resolution EVI and NDVI values of the pixels within a 50-m buffer of bird sample points using Google Earth Engine (32-day Landsat 8 Collection Tier 1). Mixed models were used to find the best associations of EVI and NDVI, predicting multiple avian diversity facets: Taxonomic diversity, functional diversity, phylogenetic diversity, specialization levels, and urban tolerance. A total of 113 bird species across 15 cities from 10 different European countries were detected. EVI mean was the best predictor for foraging substrate specialization. NDVI mean was the best predictor for most avian diversity facets: taxonomic diversity, functional richness and evenness, phylogenetic diversity, phylogenetic species variability, community evolutionary distinctiveness, urban tolerance, diet foraging behavior, and habitat richness specialists. Finally, EVI and NDVI standard deviation were not the best predictors for any avian diversity facets studied. Our findings expand previous knowledge about EVI and NDVI as surrogates of avian diversity at a continental scale. Considering the European Commission's proposal for a Nature Restoration Law calling for expanding green urban space areas by 2050, we propose NDVI as a proxy of multiple facets of avian diversity to efficiently monitor bird community responses to land use changes in the cities.

Plant community phylogeny responses to protections and its main drivers in boreal forests, China: General pattern and implications

Patterns of the phylogenetic structure have been broadly applied to predict community assembly processes. However, the distribution pattern of evolutionary diversity and its drivers under nature conservation are still poorly understood in boreal forests. Here, we investigated 1738 sampling plots and subplots from distinct protection intensities (PIs) zones in five representative National Nature Reserves (NNRs). Multiple comparisons, redundancy analysis, and linear mixed model were performed to identify the changes in community phylogeny across different PIs and NNRs and the drivers for these variations. Our results showed considerable plant community phylogeny variations in different NNRs. As indicated by SesMPD (standardized mean pairwise distance) and SesMNTD (standardized the mean nearest taxon distance), trees, shrubs, and herbs presented overdispersed, clustered, and random distribution patterns, respectively, in different PIs. Protection resulted in the phylogenetic structure between the nearest species of trees showing a more overdispersed pattern (p < 0.05). Protection decreased the phylogenetically clustered degree between the nearest species of shrubs (p > 0.05), while the herbs still maintained a random pattern. Community traits explained the most to phylogeny variation of different communities (24 %-71 %, p < 0.01), followed by geoclimatic factors (2 %-24 %) and conservation processes (1 %-21 %). The higher mean annual precipitation and under branch height at the lower latitude area accompanied the higher SesMPD and SesMNTD. The higher PIs attended with higher tree SesMPD, and the longer protection time resulted in higher shrub PSR (phylogenetic species richness) and PSV (phylogenetic species variability). Including the location of NNRs, community traits, and years of protection, rather than only emphasizing PI itself, could optimize community phylogenetic structure and preserve the evolutionary potential of biodiversity.

Applying image clustering to phylogenetic analysis: A trial

•Molecular phylogenetic analysis can be supplemented by image clustering analysis that uses pretrained machine learning tools.

Perennial herb diversity contributes more than annual herb diversity to multifunctionality in dryland ecosystems of North-western China

BACKGROUND

Considerable attention has been given to how different aspects of biodiversity sustain ecosystem functions. Herbs are a critical component of the plant community of dryland ecosystems, but the importance of different life form groups of herbs is often overlooked in experiments on biodiversity-ecosystem multifunctionality. Hence, little is known about how the multiple attributes of diversity of different life form groups of herbs affect changes to the multifunctionality of ecosystems.

METHODS

We investigated geographic patterns of herb diversity and ecosystem multifunctionality along a precipitation gradient of 2100 km in Northwest China, and assessed the taxonomic, phylogenetic and functional attributes of different life form groups of herbs on the multifunctionality.

RESULTS

We found that subordinate (richness effect) species of annual herbs and dominant (mass ratio effect) species of perennial herbs were crucial for driving multifunctionality. Most importantly, the multiple attributes (taxonomic, phylogenetic and functional) of herb diversity enhanced the multifunctionality. The functional diversity of herbs provided greater explanatory power than did taxonomic and phylogenetic diversity. In addition, the multiple attribute diversity of perennial herbs contributed more than annual herbs to multifunctionality.

CONCLUSIONS

Our findings provide insights into previously neglected mechanisms by which the diversity of different life form groups of herbs affect ecosystem multifunctionality. These results provide a comprehensive understanding of the relationship between biodiversity and multifunctionality, and will ultimately contribute to multifunctional conservation and restoration programs in dryland ecosystems.

Stratification and recovery time jointly shape ant functional reassembly in a neotropical forest

Microhabitat differentiation of species communities such as vertical stratification in tropical forests contributes to species coexistence and thus biodiversity. However, little is known about how the extent of stratification changes during forest recovery and influences community reassembly. Environmental filtering determines community reassembly in time (succession) and in space (stratification), hence functional and phylogenetic composition of species communities are highly dynamic. It is poorly understood if and how these two concurrent filters-forest recovery and stratification-interact. In a tropical forest chronosequence in Ecuador spanning 34 years of natural recovery, we investigated the recovery trajectory of ant communities in three overlapping strata (ground, leaf litter, lower tree trunk) by quantifying 13 traits, as well as the functional and phylogenetic diversity of the ants. We expected that functional and phylogenetic diversity would increase with recovery time and that each ant community within each stratum would show a distinct functional reassembly. We predicted that traits related to ant diet would show divergent trajectories reflecting an increase in niche differentiation with recovery time. On the other hand, traits related to the abiotic environment were predicted to show convergent trajectories due to a more similar microclimate across strata with increasing recovery age. Most of the functional traits and the phylogenetic diversity of the ants were clearly stratified, confirming previous findings. However, neither functional nor phylogenetic diversity increased with recovery time. Community-weighted trait means had complex relationships to recovery time and the majority were shaped by a statistical interaction between recovery time and stratum, confirming our expectations. However, most trait trajectories converged among strata with increasing recovery time regardless of whether they were related to ant diet or environmental conditions. We confirm the hypothesized interaction among environmental filters during the functional reassembly in tropical forests. Communities in individual strata respond differently to recovery, and possible filter mechanisms likely arise from both abiotic (e.g. microclimate) and biotic (e.g. diet) conditions. Since vertical stratification is prevalent across animal and plant taxa, our results highlight the importance of stratum-specific analysis in dynamic ecosystems and may generalize beyond ants.

Recent fire in a Mediterranean ecosystem strengthens hoverfly populations and their interaction networks with plants

Fire affects many critical ecological processes, including pollination, and effects of climate change on fire regimes may have profound consequences that are difficult to predict. Considerable work has examined effects of fire on pollinator diversity, but relatively few studies have examined these effects on interaction networks including those of pollinators other than bees. We examined the effects of a severe wildfire on hoverfly pollinators in a Mediterranean island system. Using data collected over 3 consecutive years at burnt and unburnt sites, we documented differences in species diversity, abundance, and functional traits, as well as hoverfly interactions with flowering plants. Hoverfly abundance and species richness peaked during the first post-fire flowering season (year 1), which coincided with the presence of many opportunistic species. Also in year 1, hoverfly pollination networks were larger, less specialized, more nested, and less modular at burnt (vs. unburnt) sites; furthermore, these networks exhibited higher phylogenetic host-plant diversity. These effects declined over the next 2 years, with burnt and unburnt sites converging in similarity to hoverfly communities and interaction networks. While data obtained over 3 years provide a clear timeline of initial post-fire recovery, we emphasize the importance of longer-term monitoring for understanding the responses of natural communities to wildfires, which are projected to become more frequent and more destructive in the future.

Global models and predictions of plant diversity based on advanced machine learning techniques

Despite the paramount role of plant diversity for ecosystem functioning, biogeochemical cycles, and human welfare, knowledge of its global distribution is still incomplete, hampering basic research and biodiversity conservation. Here, we used machine learning (random forests, extreme gradient boosting, and neural networks) and conventional statistical methods (generalized linear models and generalized additive models) to test environment-related hypotheses of broad-scale vascular plant diversity gradients and to model and predict species richness and phylogenetic richness worldwide. To this end, we used 830 regional plant inventories including c. 300 000 species and predictors of past and present environmental conditions. Machine learning showed a superior performance, explaining up to 80.9% of species richness and 83.3% of phylogenetic richness, illustrating the great potential of such techniques for disentangling complex and interacting associations between the environment and plant diversity. Current climate and environmental heterogeneity emerged as the primary drivers, while past environmental conditions left only small but detectable imprints on plant diversity. Finally, we combined predictions from multiple modeling techniques (ensemble predictions) to reveal global patterns and centers of plant diversity at multiple resolutions down to 7774 km² . Our predictive maps provide accurate estimates of global plant diversity available at grain sizes relevant for conservation and macroecology.

Application of Phytosociological Information in the Evaluation of the Management of Protected Areas

The classification system of plant communities using phytosociological methods can be applied to their conservation in protected areas, as well as in establishing adequate protections and granting legal status to such areas. A new integrative index is developed to classify plant communities for the evaluation of the conservation status of protected areas, obtained from the product of three statistical indices of diversity: Syntaxonomic Distinctness, Rarefaction and Areas Prioritisation, which has been named DRA (acronym of the three indices used). The DRA is used to assess whether the status granted to Protected Areas matches the values provided by the plant communities within them and which were the basis for the identification and description of the Habitats of Community Interest (Habitats Directive-92/43/CEE). The proposed method was applied to the network of protected natural areas on the Andalusian coast, including 14 areas with different protection status, where, once the plant communities they contain were identified, the DRA index was applied to each of them and compared with the Legal Protection Index, i.e., the current protection regime; it becomes clear, objectively, that not all the statuses assigned, whether the IUCN criteria or those of the Andalusian government, correspond to the real levels of protection they should have on the basis of their plant communities.