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

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.

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.

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.

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.

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.

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.

Environmental and topographic drivers of amphibian phylogenetic diversity and endemism in the Iberian Peninsula

Understanding the ecological and evolutionary processes driving biodiversity patterns and allowing their persistence is of utmost importance. Many hypotheses have been proposed to explain spatial diversity patterns, including water-energy availability, habitat heterogeneity, and historical climatic refugia. The main goal of this study is to identify if general spatial drivers of species diversity patterns of phylogenetic diversity (PD) and phylogenetic endemism (PE) at the global scale are also predictive of PD and PE at regional scales, using Iberian amphibians as a case study. Our main hypothesis assumes that topography along with contemporary and historical climate are drivers of phylogenetic diversity and endemism, but that the strength of these predictors may be weaker at the regional scale than it tends to be at the global scale. We mapped spatial patterns of Iberian amphibians' phylogenetic diversity and endemism, using previously published phylogenetic and distribution data. Furthermore, we compiled spatial data on topographic and climatic variables related to the water-energy availability, topography, and historical climatic instability hypotheses. To test our hypotheses, we used Spatial Autoregressive Models and selected the best model to explain diversity patterns based on Akaike Information Criterion. Our results show that, out of the variables tested in our study, water-energy availability and historical climate instability are the most important drivers of amphibian diversity in Iberia. However, as predicted, the strength of these predictors in our case study is weaker than it tends to be at global scales. Thus, additional drivers should also be investigated and we suggest caution when interpreting these predictors as surrogates for different components of diversity.

Decomposing niche components reveals simultaneous effects of opposite deterministic processes structuring alpine small mammal assembly

Introduction

Species distribution in alpine areas is constrained by multiple abiotic and biotic stressors. This leads to discrepant assembly patterns between different locations and study objects as opposite niche-based processes—limiting similarity and habitat filtering—simultaneously structure communities, masking overall patterns. We aimed to address how these processes structure small mammal communities in the alpine tree line transition zone, one of the most distinct vegetation transitions between alpine and montane habitats.

Methods

We compiled a dataset of species checklist, phylogeny, and functional traits from field collection and published sources spanning 18 mountains in southwest China. We first examined hypothetical niche-based processes with frequently used phylogenetic and trait approaches using this dataset. The species traits were decomposed into different niche components to explore the respective effects of specific stressors. Indices representing evolutionary history, trait space, and pairwise species distance were estimated and compared with null model expectations. Linear mixed-effect models were used to assess the association patterns between diversity indices and elevation.

Results

The results indicated that phylogenetic and functional richness were positively correlated with species richness. In contrast, distance-based indices were either negatively or weakly positively correlated with species richness. Null model analyses suggested no evidence of non-random phylogenetic or overall trait patterns. However, the resource acquisition niche tended to be more overdispersed (positive slopes), while the habitat affinity niche tended to be more clustered (negative slopes) beyond the high elevation tree line.

Discussion

These findings show that opposite niche-based processes simultaneously structure small mammal communities in alpine areas. Overall, the present study provides vital insights into the complexity of assembly processes in these habitats. It also highlights the importance of relating relevant traits to distinguish the influences of specific abiotic and biotic stressors.

A simplified protocol for DNA extraction from FTA cards for faecal microbiome studies

As metagenomic studies continue to increase in size and complexity, they are often required to incorporate data from geographically isolated locations or longitudinal time samples. This represents a technical challenge, given that many of the commonly used methods used for sample collection, storage, and DNA extraction are sensitive to differences related to the time, storage and chemistry involved. FTA cards have been previously proposed as a simple, reliable and cost-efficient method for the preservation of animal faecal microbiomes. In this study, we report a simplified extraction methodology for recovering microbiome DNA from faeces stored on FTA cards and compare its performance to a common alternative means of characterising such microbiomes; namely, immediate freezing of the faeces followed by DNA extraction using the Qiagen PowerSoil DNA isolation kit. Our results show that overall the application of our simplified DNA extraction methodology yields microbial community results that have higher diversity and an expanded core microbiome than that found using the PowerSoil methodology. This suggests that the FTA card extraction method presented here is a viable alternative for metagenomic studies using faecal material when traditional freeze-based storage methods are not feasible.

Bat diversity is driven by elevation and distance to the nearest watercourse in a terra firme forest in the northeastern Brazilian Amazon

Variations in environmental conditions along gradients play an important role in species distribution through environmental filtering of morphological and physiological traits; however, their effects on bat diversity remain poorly understood. Here, we investigate the effect of the distance to the nearest watercourse, terrain elevation, vegetation clutter, basal area and canopy height on taxonomic, functional and phylogenetic diversity and on the predominance of some functional traits (body mass, wing morphology and trophic level) of bat assemblages (phyllostomid and mormoopid bats) in a terra firme forest, in the northeastern Brazilian Amazon. We captured bats using mist nets in 15 permanent plots over a 25 km2 area of continuous forest. We captured 279 individuals belonging to 28 species with a total of 77.760 m2.h of sampling effort. Our results showed that bat richness increases as a function of distance to the nearest watercourse and that the assemblage also changes, with more diverse taxonomic and functional groups in areas further from the watercourse. Furthermore, elevation positively affects species richness, and the basal area of the forest positively influences the average body mass of bats. Taken together, our results demonstrate that subtle variations in the environmental conditions along a local scale gradient impact on the main dimensions of bat diversity in primary forests.

Changes in Phylogenetic and Functional Diversity of Ciliates along the Course of a Mediterranean Karstic River

Ciliates are a group of phagotrophic protists found in a wide variety of ecosystems. This study builds on recent studies of ciliates in the Krka river and investigates changes in the phylogenetic and functional diversity of ciliates in biofilm to predict the phylogenetic and functional structure of ciliates in other karstic rivers. Biofilm samples were collected from four representative locations: upstream (Krka spring), midstream (Marasovine), and downstream (Roški slap, Skradinski buk) of the Krka river to test for differences in phylogenetic and functional diversity of ciliates in relation to location and positioning on tufa stones (light/dark-exposed side of tufa stone). Our results showed that Krka spring had higher phylogenetic species variability, lower phylogenetic diversity, and lower functional richness than Skradinski buk, suggesting phylogenetic overdispersal at Krka spring. This could be due to environmental filtering, competitive exclusion, or a combination of these factors. As the first study of its kind in the Mediterranean, our results shed light on the phylogenetic and functional diversity of ciliates in karst ecosystems and provide a basis for future ecological and conservation efforts.

The community ecology perspective of omics data

The measurement of uncharacterized pools of biological molecules through techniques such as metabarcoding, metagenomics, metatranscriptomics, metabolomics, and metaproteomics produces large, multivariate datasets. Analyses of these datasets have successfully been borrowed from community ecology to characterize the molecular diversity of samples (ɑ-diversity) and to assess how these profiles change in response to experimental treatments or across gradients (β-diversity). However, sample preparation and data collection methods generate biases and noise which confound molecular diversity estimates and require special attention. Here, we examine how technical biases and noise that are introduced into multivariate molecular data affect the estimation of the components of diversity (i.e., total number of different molecular species, or entities; total number of molecules; and the abundance distribution of molecular entities). We then explore under which conditions these biases affect the measurement of ɑ- and β-diversity and highlight how novel methods commonly used in community ecology can be adopted to improve the interpretation and integration of multivariate molecular data. Video Abstract.

Hedgerows increase the diversity and modify the composition of arbuscular mycorrhizal fungi in Mediterranean agricultural landscapes

Sustainable agriculture is essential to address global challenges such as climate change and biodiversity loss. Hedgerows enhance aboveground biodiversity and provide ecosystem services, but little is known about their impact on soil biota. Arbuscular mycorrhizal (AM) fungi are one of the key components of belowground biodiversity. We compared the diversity and composition of AM fungal communities at four farmland sites located in Central Spain, where 132 soil samples in total were collected to assess soil physical and chemical properties and the AM fungal communities. We compared the richness (number of AM fungal taxa), taxonomic, functional, and phylogenetic diversity, and structure of the AM fungal communities across three farmland habitat types, namely hedgerows, woody crops (olive groves and vineyard), and herbaceous crops (barley, sunflower, and wheat). Our results showed positive effects of hedgerows on most diversity metrics. Almost 60% of the AM fungal taxa were shared among the three farmland habitat types. Hedgerows increased AM fungal taxonomic richness (31%) and alpha diversity (25%), and especially so compared to herbaceous crops (45% and 28%, respectively). Hedgerows harbored elevated proportions of AM fungi with non-ruderal life-history strategies. AM fungal communities were more similar between hedgerows and woody crops than between hedgerows and adjacent herbaceous crops, possibly because of differences in tillage and fertilization. Unexpectedly, hedgerows reduced phylogenetic diversity, which might be related to more selective associations of AM fungi with woody plants than with herbaceous crops. Overall, the results suggest that planting hedgerows contributes to maintain belowground diversity. Thus, European farmers should plant more hedgerows to attain the goals of the EU Biodiversity Strategy for 2030.