Rubén G. Mateo

Do stacked species distribution models reflect altitudinal diversity patterns

The objective of this study was to evaluate the performance of stacked species distribution models in predicting the alpha and gamma species diversity patterns of two important plant clades along elevation in the Andes. We modelled the distribution of the species in the Anthurium genus (53 species) and the Bromeliaceae family (89 species) using six modelling techniques. We combined all of the predictions for the same species in ensemble models based on two different criteria: the average of the rescaled predictions by all techniques and the average of the best techniques. The rescaled predictions were then reclassified into binary predictions (presence/absence). By stacking either the original predictions or binary predictions for both ensemble procedures, we obtained four different species richness models per taxa. The gamma and alpha diversity per elevation band (500 m) was also computed. To evaluate the prediction abilities for the four predictions of species richness and gamma diversity, the models were compared with the real data along an elevation gradient that was independently compiled by specialists. Finally, we also tested whether our richness models performed better than a null model of altitudinal changes of diversity based on the literature. Stacking of the ensemble prediction of the individual species models generated richness models that proved to be well correlated with the observed alpha diversity richness patterns along elevation and with the gamma diversity derived from the literature. Overall, these models tend to overpredict species richness. The use of the ensemble predictions from the species models built with different techniques seems very promising for modelling of species assemblages. Stacking of the binary models reduced the over-prediction, although more research is needed. The randomisation test proved to be a promising method for testing the performance of the stacked models, but other implementations may still be developed.

Approximate Bayesian Computation Reveals the Crucial Role of Oceanic Islands for the Assembly of Continental Biodiversity.

The perceived low levels of genetic diversity, poor interspecific competitive and defensive ability, and loss of dispersal capacities of insular lineages have driven the view that oceanic islands are evolutionary dead ends. Focusing on the Atlantic bryophyte flora distributed across the archipelagos of the Azores, Madeira, the Canary Islands, Western Europe, and northwestern Africa, we used an integrative approach with species distribution modeling and population genetic analyses based on approximate Bayesian computation to determine whether this view applies to organisms with inherent high dispersal capacities. Genetic diversity was found to be higher in island than in continental populations, contributing to mounting evidence that, contrary to theoretical expectations, island populations are not necessarily genetically depauperate. Patterns of genetic variation among island and continental populations consistently fitted those simulated under a scenario of de novo foundation of continental populations from insular ancestors better than those expected if islands would represent a sink or a refugium of continental biodiversity. We, suggest that the northeastern Atlantic archipelagos have played a key role as a stepping stone for transoceanic migrants. Our results challenge the traditional notion that oceanic islands are the end of the colonization road and illustrate the significant role of oceanic islands as reservoirs of novel biodiversity for the assembly of continental floras.

How do temperate bryophytes face the challenge of a changing environment? Lessons from the past and predictions for the future

Bryophytes are a group of early land plants, whose specific ecophysiological and biological features, including poikilohydry, sensitivity to moderately high temperature and high dispersal ability, make them ideal candidates for investigating the impact of climate changes. Employing a combined approach of species distribution modelling (SDM) and molecular phylogeography in the temperate moss Homalothecium sericeum, we explore the significance of the Mediterranean refugia, contrasting the southern and northern refugia hypotheses, determine the extent to which recolonization of previously glaciated areas has been facilitated by the high dispersal ability of the species and make predictions on the extent to which it will be impacted by ongoing climate change. The Mediterranean areas exhibit the highest nucleotidic diversities and host a mixture of ancestral, endemic and more recently derived haplotypes. Extra-Mediterranean areas exhibit low genetic diversities and Euro-Siberian populations display a significant signal of expansion that is identified to be of Euro-Siberian origin, pointing to the northern refugia hypothesis. The SDMs predict a global net increase in range size owing to ongoing climate change, but substantial range reductions in southern areas. Presence of a significant phylogeographical signal at different spatial scales suggests, however, that dispersal limitations might constitute, as opposed to the traditional view of spore-producing plants as efficient dispersers, a constraint for migration. This casts doubts about the ability of the species to face the massive extinctions predicted in the southern areas, threatening their status of reservoir of genetic diversity.

Modeling Species Distributions from Heterogeneous Data for the Biogeographic Regionalization of the European Bryophyte Flora

The definition of biogeographic regions provides a fundamental framework for a range of basic and applied questions in biogeography, evolutionary biology, systematics and conservation. Previous research suggested that environmental forcing results in highly congruent regionalization patterns across taxa, but that the size and number of regions depends on the dispersal ability of the taxa considered. We produced a biogeographic regionalization of European bryophytes and hypothesized that (1) regions defined for bryophytes would differ from those defined for other taxa due to the highly specific eco-physiology of the group and (2) their high dispersal ability would result in the resolution of few, large regions. Species distributions were recorded using 10,000 km2 MGRS pixels. Because of the lack of data across large portions of the area, species distribution models employing macroclimatic variables as predictors were used to determine the potential composition of empty pixels. K-means clustering analyses of the pixels based on their potential species composition were employed to define biogeographic regions. The optimal number of regions was determined by v-fold cross-validation and Moran’s I statistic. The spatial congruence of the regions identified from their potential bryophyte assemblages with large-scale vegetation patterns is at odds with our primary hypothesis. This reinforces the notion that post-glacial migration patterns might have been much more similar in bryophytes and vascular plants than previously thought. The substantially lower optimal number of clusters and the absence of nested patterns within the main biogeographic regions, as compared to identical analyses in vascular plants, support our second hypothesis. The modelling approach implemented here is, however, based on many assumptions that are discussed but can only be tested when additional data on species distributions become available, highlighting the substantial importance of developing integrated mapping projects for all taxa in key biogeographically areas of Europe, and the Mediterranean peninsulas in particular.

The mossy north: an inverse latitudinal diversity gradient in European bryophytes

It remains hotly debated whether latitudinal diversity gradients are common across taxonomic groups and whether a single mechanism can explain such gradients. Investigating species richness (SR) patterns of European land plants, we determine whether SR increases with decreasing latitude, as predicted by theory, and whether the assembly mechanisms differ among taxonomic groups. SR increases towards the south in spermatophytes, but towards the north in ferns and bryophytes. SR patterns in spermatophytes are consistent with their patterns of beta diversity, with high levels of nestedness and turnover in the north and in the south, respectively, indicating species exclusion towards the north and increased opportunities for speciation in the south. Liverworts exhibit the highest levels of nestedness, suggesting that they represent the most sensitive group to the impact of past climate change. Nevertheless, although the extent of liverwort species turnover in the south is substantially and significantly lower than in spermatophytes, liverworts share with the latter a higher nestedness in the north and a higher turn-over in the south, in contrast to mosses and ferns. The extent to which the similarity in the patterns displayed by spermatophytes and liverworts reflects a similar assembly mechanism remains, however, to be demonstrated.

Legume diversity patterns in West Central Africa: influence of species biology on distribution models.

Objectives Species Distribution Models (SDMs) are used to produce predictions of potential Leguminosae diversity in West Central Africa. Those predictions are evaluated subsequently using expert opinion. The established methodology of combining all SDMs is refined to assess species diversity within five defined vegetation types. Potential species diversity is thus predicted for each vegetation type respectively. The primary aim of the new methodology is to define, in more detail, areas of species richness for conservation planning. Methodology Using Maxent, SDMs based on a suite of 14 environmental predictors were generated for 185 West Central African Leguminosae species, each categorised according to one of five vegetation types: Afromontane, coastal, non-flooded forest, open formations, or riverine forest. The relative contribution of each environmental variable was compared between different vegetation types using a nonparametric Kruskal-Wallis analysis followed by a post-hoc Kruskal-Wallis Paired Comparison contrast. Legume species diversity patterns were explored initially using the typical method of stacking all SDMs. Subsequently, five different ensemble models were generated by partitioning SDMs according to vegetation category. Ecological modelers worked with legume specialists to improve data integrity and integrate expert opinion in the interpretation of individual species models and potential species richness predictions for different vegetation types. Results/Conclusions Of the 14 environmental predictors used, five showed no difference in their relative contribution to the different vegetation models. Of the nine discriminating variables, the majority were related to temperature variation. The set of variables that played a major role in the Afromontane species diversity model differed significantly from the sets of variables of greatest relative important in other vegetation categories. The traditional approach of stacking all SDMs indicated overall centers of diversity in the region but the maps indicating potential species richness by vegetation type offered more detailed information on which conservation efforts can be focused.

Climate threat on the Macaronesian endemic bryophyte flora

Oceanic islands are of fundamental importance for the conservation of biodiversity because they exhibit high endemism rates coupled with fast extinction rates. Nowhere in Europe is this pattern more conspicuous than in the Macaronesian biogeographic region. A large network of protected areas within the region has been developed, but the question of whether these areas will still be climatically suitable for the globally threatened endemic element in the coming decades remains open. Here, we make predictions on the fate of the Macaronesian endemic bryophyte flora in the context of ongoing climate change. The potential distribution of 35 Macaronesian endemic bryophyte species was assessed under present and future climate conditions using an ensemble modelling approach. Projections of the models under different climate change scenarios predicted an average decrease of suitable areas of 62–87% per species and a significant elevational increase by 2070, so that even the commonest species were predicted to fit either the Vulnerable or Endangered IUCN categories. Complete extinctions were foreseen for six of the studied Macaronesian endemic species. Given the uncertainty regarding the capacity of endemic species to track areas of suitable climate within and outside the islands, active management associated to an effective monitoring program is suggested.

Biodiversity Models: What If Unsaturation Is the Rule?

Improving biodiversity predictions is essential if we are to meet the challenges posed by global change. As knowledge is key to feed models, we need to evaluate how debated theory can affect models. An important ongoing debate is whether environmental constraints limit the number of species that can coexist in a community (saturation), with recent findings suggesting that species richness in many communities might be unsaturated. Here, we propose that biodiversity models could address this issue by accounting for a duality: considering communities as unsaturated but where species composition is constrained by different scale-dependent biodiversity drivers. We identify a variety of promising advances for incorporating this duality into commonly applied biodiversity modelling approaches and improving their spatial predictions.

Improving spatial predictions of taxonomic, functional and phylogenetic diversity

In this study, we compare two community modelling approaches to determine their ability to predict the taxonomic, functional and phylogenetic properties of plant assemblages along a broad elevation gradient and at a fine resolution. The first method is the standard stacking individual species distribution modelling (SSDM) approach, which applies a simple environmental filter to predict species assemblages. The second method couples the SSDM and macroecological modelling (MEMSSDM-MEM) approaches to impose a limit on the number of species co-occurring at each site. Because the detection of diversity patterns can be influenced by different levels of phylogenetic or functional trees, we also examine whether performing our analyses from broad to more exact structures in the trees influences the performance of the two modelling approaches when calculating diversity indices. We found that coupling the SSDM with the MEM improves the overall predictions for the three diversity facets compared with those of the SSDM alone. The accuracy of the SSDM predictions for the diversity indices varied greatly along the elevation gradient, and when considering broad to more exact structure in the functional and phylogenetic trees, the SSDM-MEM predictions were more stable. SSDM-MEM moderately but significantly improved the prediction of taxonomic diversity, which was mainly driven by the corrected number of predicted species. The performance of both modelling frameworks increased when predicting the functional and phylogenetic diversity indices. In particular, fair predictions of the taxonomic composition by SSDM-MEM led to increasingly accurate predictions of the functional and phylogenetic indices, suggesting that the compositional errors were associated with species that were functionally or phylogenetically close to the correct ones; however, this did not always hold for the SSDM predictions.Synthesis. In this study, we tested the use of a recently published approach that couples species distribution and macroecological models to provide the first predictions of the distribution of multiple facets of plant diversity: taxonomic, functional and phylogenetic. Moderate but significant improvements were obtained; thus, our results open promising avenues for improving our ability to predict the different facets of biodiversity in space and time across broad environmental gradients when functional and phylogenetic information is available.

How to describe species richness patterns for bryophyte conservation

Abstract A large amount of data for inconspicuous taxa is stored in natural history collections; however, this information is often neglected for biodiversity patterns studies. Here, we evaluate the performance of direct interpolation of museum collections data, equivalent to the traditional approach used in bryophyte conservation planning, and stacked species distribution models (S‐SDMs) to produce reliable reconstructions of species richness patterns, given that differences between these methods have been insufficiently evaluated for inconspicuous taxa. Our objective was to contrast if species distribution models produce better inferences of diversity richness than simply selecting areas with the higher species numbers. As model species, we selected Iberian species of the genus Grimmia (Bryophyta), and we used four well‐collected areas to compare and validate the following models: 1) four Maxent richness models, each generated without the data from one of the four areas, and a reference model created using all of the data and 2) four richness models obtained through direct spatial interpolation, each generated without the data from one area, and a reference model created with all of the data. The correlations between the partial and reference Maxent models were higher in all cases (0.45 to 0.99), whereas the correlations between the spatial interpolation models were negative and weak (−0.3 to −0.06). Our results demonstrate for the first time that S‐SDMs offer a useful tool for identifying detailed richness patterns for inconspicuous taxa such as bryophytes and improving incomplete distributions by assessing the potential richness of under‐surveyed areas, filling major gaps in the available data. In addition, the proposed strategy would enhance the value of the vast number of specimens housed in biological collections.