Predicting species distribution: offering more than simple habitat models

n the last two decades, interest in species distribution models (SDMs) of plants and animals has grown dramatically. Recent advances in SDMs allow us to potentially forecast anthropogenic effects on patterns of biodiversity at different spatial scales. However, some limitations still preclude the use of SDMs in many theoretical and practical applications. Here, we provide an overview of recent advances in this field, discuss the ecological principles and assumptions underpinning SDMs, and highlight critical limitations and decisions inherent in the construction and evaluation of SDMs. Particular emphasis is given to the use of SDMs for the assessment of climate change impacts and conservation management issues. We suggest new avenues for incorporating species migration, population dynamics, biotic interactions and community ecology into SDMs at multiple spatial scales. Addressing all these issues requires a better integration of SDMs with ecological theory.

Evolutionary branching in deme-structured populations.

Adaptive dynamics shows that a continuous trait under frequency dependent selection may first converge to a singular point followed by spontaneous transition from a unimodal trait distribution into a bimodal one, which is called "evolutionary branching". Here, we study evolutionary branching in a deme-structured population by constructing a quantitative genetic model for the trait variance dynamics, which allows us to obtain an analytic condition for evolutionary branching. This is first shown to agree with previous conditions for branching expressed in terms of relatedness between interacting individuals within demes and obtained from mutant-resident systems. We then show this branching condition can be markedly simplified when the evolving trait affect fecundity and/or survival, as opposed to affecting population structure, which would occur in the case of the evolution of dispersal. As an application of our model, we evaluate the threshold migration rate below which evolutionary branching cannot occur in a pairwise interaction game. This agrees very well with the individual-based simulation results.

Unfolding nature "in silico" : overcoming practical and metodological limitations of species distribution models

Abstract : The existence of a causal relationship between the spatial distribution of living organisms and their environment, in particular climate, has been long recognized and is the central principle of biogeography. In turn, this recognition has led scientists to the idea of using the climatic, topographic, edaphic and biotic characteristics of the environment to predict its potential suitability for a given species or biological community. In this thesis, my objective is to contribute to the development of methodological improvements in the field of species distribution modeling. More precisely, the objectives are to propose solutions to overcome limitations of species distribution models when applied to conservation biology issues, or when .used as an assessment tool of the potential impacts of global change. The first objective of my thesis is to contribute to evidence the potential of species distribution models for conservation-related applications. I present a methodology to generate pseudo-absences in order to overcome the frequent lack of reliable absence data. I also demonstrate, both theoretically (simulation-based) and practically (field-based), how species distribution models can be successfully used to model and sample rare species. Overall, the results of this first part of the thesis demonstrate the strong potential of species distribution models as a tool for practical applications in conservation biology. The second objective this thesis is to contribute to improve .projections of potential climate change impacts on species distributions, and in particular for mountain flora. I develop and a dynamic model, MIGCLIM, that allows the implementation of dispersal limitations into classic species distribution models and present an application of this model to two virtual species. Given that accounting for dispersal limitations requires information on seed dispersal, distances, a general methodology to classify species into broad dispersal types is also developed. Finally, the M~GCLIM model is applied to a large number of species in a study area of the western Swiss Alps. Overall, the results indicate that while dispersal limitations can have an important impact on the outcome of future projections of species distributions under climate change scenarios, estimating species threat levels (e.g. species extinction rates) for a mountainous areas of limited size (i.e. regional scale) can also be successfully achieved when considering dispersal as unlimited (i.e. ignoring dispersal limitations, which is easier from a practical point of view). Finally, I present the largest fine scale assessment of potential climate change impacts on mountain vegetation that has been carried-out to date. This assessment involves vegetation from 12 study areas distributed across all major western and central European mountain ranges. The results highlight that some mountain ranges (the Pyrenees and the Austrian Alps) are expected to be more affected by climate change than others (Norway and the Scottish Highlands). The results I obtain in this study also indicate that the threat levels projected by fine scale models are less severe than those derived from coarse scale models. This result suggests that some species could persist in small refugias that are not detected by coarse scale models. Résumé : L'existence d'une relation causale entre la répartition des espèces animales et végétales et leur environnement, en particulier le climat, a été mis en évidence depuis longtemps et est un des principes centraux en biogéographie. Ce lien a naturellement conduit à l'idée d'utiliser les caractéristiques climatiques, topographiques, édaphiques et biotiques de l'environnement afin d'en prédire la qualité pour une espèce ou une communauté. Dans ce travail de thèse, mon objectif est de contribuer au développement d'améliorations méthodologiques dans le domaine de la modélisation de la distribution d'espèces dans le paysage. Plus précisément, les objectifs sont de proposer des solutions afin de surmonter certaines limitations des modèles de distribution d'espèces dans des applications pratiques de biologie de la conservation ou dans leur utilisation pour évaluer l'impact potentiel des changements climatiques sur l'environnement. Le premier objectif majeur de mon travail est de contribuer à démontrer le potentiel des modèles de distribution d'espèces pour des applications pratiques en biologie de la conservation. Je propose une méthode pour générer des pseudo-absences qui permet de surmonter le problème récurent du manque de données d'absences fiables. Je démontre aussi, de manière théorique (par simulation) et pratique (par échantillonnage de terrain), comment les modèles de distribution d'espèces peuvent être utilisés pour modéliser et améliorer l'échantillonnage des espèces rares. Ces résultats démontrent le potentiel des modèles de distribution d'espèces comme outils pour des applications de biologie de la conservation. Le deuxième objectif majeur de ce travail est de contribuer à améliorer les projections d'impacts potentiels des changements climatiques sur la flore, en particulier dans les zones de montagnes. Je développe un modèle dynamique de distribution appelé MigClim qui permet de tenir compte des limitations de dispersion dans les projections futures de distribution potentielle d'espèces, et teste son application sur deux espèces virtuelles. Vu que le fait de prendre en compte les limitations dues à la dispersion demande des données supplémentaires importantes (p.ex. la distance de dispersion des graines), ce travail propose aussi une méthode de classification simplifiée des espèces végétales dans de grands "types de disperseurs", ce qui permet ainsi de d'obtenir de bonnes approximations de distances de dispersions pour un grand nombre d'espèces. Finalement, j'applique aussi le modèle MIGCLIM à un grand nombre d'espèces de plantes dans une zone d'études des pré-Alpes vaudoises. Les résultats montrent que les limitations de dispersion peuvent avoir un impact considérable sur la distribution potentielle d'espèces prédites sous des scénarios de changements climatiques. Cependant, quand les modèles sont utilisés pour évaluer les taux d'extinction d'espèces dans des zones de montages de taille limitée (évaluation régionale), il est aussi possible d'obtenir de bonnes approximations en considérant la dispersion des espèces comme illimitée, ce qui est nettement plus simple d'un point dé vue pratique. Pour terminer je présente la plus grande évaluation à fine échelle d'impact potentiel des changements climatiques sur la flore des montagnes conduite à ce jour. Cette évaluation englobe 12 zones d'études réparties sur toutes les chaines de montages principales d'Europe occidentale et centrale. Les résultats montrent que certaines chaines de montagnes (les Pyrénées et les Alpes Autrichiennes) sont projetées comme plus sensibles aux changements climatiques que d'autres (les Alpes Scandinaves et les Highlands d'Ecosse). Les résultats obtenus montrent aussi que les modèles à échelle fine projettent des impacts de changement climatiques (p. ex. taux d'extinction d'espèces) moins sévères que les modèles à échelle large. Cela laisse supposer que les modèles a échelle fine sont capables de modéliser des micro-niches climatiques non-détectées par les modèles à échelle large.

The role of climatic tolerances and seed traits in reduced extinction rates of temperate polygonaceae.

The latitudinal diversity gradient (LDG) is one of the most striking and consistent biodiversity patterns across taxonomic groups. We investigate the species richness gradient in the buckwheat family, Polygonaceae, which exhibits a reverse LDG and is, thus, decoupled from dominant gradients of energy and environmental stability that increase toward the tropics and confound mechanistic interpretations. We test competing age and evolutionary diversification hypotheses, which may explain the diversification of this plant family over the past 70 million years. Our analyses show that the age hypothesis, which posits that clade richness is positively correlated with the ecological and evolutionary time since clade origin, fails to explain the richness gradient observed in Polygonaceae. However, an evolutionary diversification hypothesis is highly supported, with diversification rates being 3.5 times higher in temperate clades compared to tropical clades. We demonstrate that differences in rates of speciation, migration, and molecular evolution insufficiently explain the observed patterns of differential diversification rates. We suggest that reduced extinction rates in temperate clades may be associated with adaptive responses to selection, through which seed morphology and climatic tolerances potentially act to minimize risk in temporally variable environments. Further study is needed to understand causal pathways among these traits and factors correlated with latitude.

An evolutionary model explaining the Neolithic transition from egalitarianism to leadership and despotism.

The Neolithic was marked by a transition from small and relatively egalitarian groups to much larger groups with increased stratification. But, the dynamics of this remain poorly understood. It is hard to see how despotism can arise without coercion, yet coercion could not easily have occurred in an egalitarian setting. Using a quantitative model of evolution in a patch-structured population, we demonstrate that the interaction between demographic and ecological factors can overcome this conundrum. We model the coevolution of individual preferences for hierarchy alongside the degree of despotism of leaders, and the dispersal preferences of followers. We show that voluntary leadership without coercion can evolve in small groups, when leaders help to solve coordination problems related to resource production. An example is coordinating construction of an irrigation system. Our model predicts that the transition to larger despotic groups will then occur when: (i) surplus resources lead to demographic expansion of groups, removing the viability of an acephalous niche in the same area and so locking individuals into hierarchy; (ii) high dispersal costs limit followers' ability to escape a despot. Empirical evidence suggests that these conditions were probably met, for the first time, during the subsistence intensification of the Neolithic.

Temporal and spatial origin of Gesneriaceae in the New World inferred from plastid DNA sequences

Gesneriaceae are represented in the New World (NW) by a major clade (c. 1000 species) currently recognized as subfamily Gesnerioideae. Radiation of this group occurred in all biomes of tropical America and was accompanied by extensive phenotypic and ecological diversification. Here we performed phylogenetic analyses using DNA sequences from three plastid loci to reconstruct the evolutionary history of Gesnerioideae and to investigate its relationship with other lineages of Gesneriaceae and Lamiales. Our molecular data confirm the inclusion of the South Pacific Coronanthereae and the Old World (OW) monotypic genus Titanotrichum in Gesnerioideae and the sister-group relationship of this subfamily to the rest of the OW Gesneriaceae. Calceolariaceae and the NW genera Peltanthera and Sanango appeared successively sister to Gesneriaceae, whereas Cubitanthus, which has been previously assigned to Gesneriaceae, is shown to be related to Linderniaceae. Based on molecular dating and biogeographical reconstruction analyses, we suggest that ancestors of Gesneriaceae originated in South America during the Late Cretaceous. Distribution of Gesneriaceae in the Palaeotropics and Australasia was inferred as resulting from two independent long-distance dispersals during the Eocene and Oligocene, respectively. In a short time span starting at 34 Mya, ancestors of Gesnerioideae colonized several Neotropical regions including the tropical Andes, Brazilian Atlantic forest, cerrado, Central America and the West Indies. Subsequent diversification within these areas occurred largely in situ and was particularly extensive in the mountainous systems of the Andes, Central America and the Brazilian Atlantic forest. Only two radiations account for 90% of the diversity of Gesneriaceae in the Brazilian Atlantic forest, whereas half of the species richness in the northern Andes and Central America originated during the last 10 Myr from a single radiation.

Ecological assembly rules in plant communities--approaches, patterns and prospects.

Understanding how communities of living organisms assemble has been a central question in ecology since the early days of the discipline. Disentangling the different processes involved in community assembly is not only interesting in itself but also crucial for an understanding of how communities will behave under future environmental scenarios. The traditional concept of assembly rules reflects the notion that species do not co-occur randomly but are restricted in their co-occurrence by interspecific competition. This concept can be redefined in a more general framework where the co-occurrence of species is a product of chance, historical patterns of speciation and migration, dispersal, abiotic environmental factors, and biotic interactions, with none of these processes being mutually exclusive. Here we present a survey and meta-analyses of 59 papers that compare observed patterns in plant communities with null models simulating random patterns of species assembly. According to the type of data under study and the different methods that are applied to detect community assembly, we distinguish four main types of approach in the published literature: species co-occurrence, niche limitation, guild proportionality and limiting similarity. Results from our meta-analyses suggest that non-random co-occurrence of plant species is not a widespread phenomenon. However, whether this finding reflects the individualistic nature of plant communities or is caused by methodological shortcomings associated with the studies considered cannot be discerned from the available metadata. We advocate that more thorough surveys be conducted using a set of standardized methods to test for the existence of assembly rules in data sets spanning larger biological and geographical scales than have been considered until now. We underpin this general advice with guidelines that should be considered in future assembly rules research. This will enable us to draw more accurate and general conclusions about the non-random aspect of assembly in plant communities.

Inbreeding and sex-biased gene flow in the ant Formica exsecta.

The objective of this study was to assess breeding and dispersal patterns of both males and females in a monogyne (a single queen per colony) population of ants. Monogyny is commonly associated with extensive nuptial flights, presumably leading to considerable gene flow over large areas. Opposite to these expectations we found evidence of both inbreeding and sex-biased gene flow in a monogyne population of Formica exsecta. We found a significant degree of population subdivision at a local scale (within islands) for queens (females heading established colonies) and workers, but not for colony fathers (the males mated to the colony queens). However, we found little evidence of population subdivision at a larger scale (among islands). More conclusive support for sex-biased gene flow comes from the analysis of isolation by distance on the largest island, and from assignment tests revealing differences in female and male philopatry. The genetic similarity between pairs of queens decreased significantly when geographical distance increased, demonstrating limited dispersal and isolation by distance in queens. By contrast, we found no such pattern for colony fathers. Furthermore, a significantly greater fraction of colony queens were assigned as having originated from the population of residence, as compared to colony fathers. Inbreeding coefficients were significantly positive for workers, but not for mother queens. The queen-male relatedness coefficient of 0.23 (regression relatedness) indicates that mating occurs between fairly close relatives. These results suggest that some monogyne species of ants have complex dispersal and mating systems that can result in genetic isolation by distance over small geographical scales. More generally, this study also highlights the importance of identifying the relevant scale in analyses of population structure and dispersal.

Queen acceptance in a socially polymorphic ant

A central question in social evolution is what processes regulate the number of breeders in each social group. Here, we tested whether differences in the rate of acceptance of new queens by resident workers could be a proximate cause explaining the coexistence of single- and multiple-queen colonies in an ant population. We found that Formica selysi workers discriminated against foreign (non-nestmate) queens, which contributes to maintaining the genetic integrity of the social group essential to kin selection. All the young and newly mated foreign queens introduced into experimental groups of workers died within 48 h. In contrast, workers frequently accepted young newly mated nestmate queens. The survival of nestmate queens was not significantly lower in groups of workers originating from single- queen colonies than in groups of workers originating from multiple-queen colonies. Finally, virgin queens had significantly higher survival than mated queens. Together, these results show that the maintenance of single-queen and multiple-queen colonies in the same population is unlikely to be caused by strong differences between the two types of colonies in their rate of acceptance of new queens by workers. They also suggest that the discrimination of queens by resident workers restricts the dispersal of foreign queens among colonies, but not the acceptance of additional nestmate queens.

Potential impact of climate change on vegetation in the European Alps: a review

Based on conclusions drawn from general climatic impact assessment in mountain regions, the review synthesizes results relevant to the European Alps published mainly from 1994 onward in the fields of population genetics, ecophysiology, phenology, phytogeography, modeling, paleoecology and vegetation dynamics. Other important factors of global change interacting synergistically with climatic factors are also mentioned, such as atmospheric CO2 concentration, eutrophication, ozone or changes in land-use. Topics addressed are general species distribution and populations (persistence, acclimation, genetic variability, dispersal, fragmentation, plant/animal interaction, species richness, conservation), potential response of vegetation (ecotonal shift - area, physiography - changes in the composition, structural changes), phenology, growth and productivity, and landscape. In conclusion, the European Alps appear to have a natural inertia and thus to tolerate an increase of 1-2 K of mean air temperature as far as plant species and ecosystems are concerned in general. However, the impact of land-use is very likely to negate this buffer in many areas. For a change of the order of 3 K or more, profound changes may be expected.

The abrupt climate change at the Eocene-Oligocene boundary and the emergence of South-East Asia triggered the spread of sapindaceous lineages.

Background and Aims Paleoclimatic data indicate that an abrupt climate change occurred at the Eocene-Oligocene (E-O) boundary affecting the distribution of tropical forests on Earth. The same period has seen the emergence of South-East (SE) Asia, caused by the collision of the Eurasian and Australian plates. How the combination of these climatic and geomorphological factors affected the spatio-temporal history of angiosperms is little known. This topic is investigated by using the worldwide sapindaceous clade as a case study. Methods Analyses of divergence time inference, diversification and biogeography (constrained by paleogeography) are applied to a combined plastid and nuclear DNA sequence data set. Biogeographical and diversification analyses are performed over a set of trees to take phylogenetic and dating uncertainty into account. Results are analysed in the context of past climatic fluctuations. Key Results An increase in the number of dispersal events at the E-O boundary is recorded, which intensified during the Miocene. This pattern is associated with a higher rate in the emergence of new genera. These results are discussed in light of the geomorphological importance of SE Asia, which acted as a tropical bridge allowing multiple contacts between areas and additional speciation across landmasses derived from Laurasia and Gondwana. Conclusions This study demonstrates the importance of the combined effect of geomorphological (the emergence of most islands in SE Asia approx. 30 million years ago) and climatic (the dramatic E-O climate change that shifted the tropical belt and reduced sea levels) factors in shaping species distribution within the sapindaceous clade.

Invader immunology: invasion history alters immune system function in cane toads (Rhinella marina) in tropical Australia.

Because an individual's investment into the immune system may modify its dispersal rate, immune function may evolve rapidly in an invader. We collected cane toads (Rhinella marina) from sites spanning their 75-year invasion history in Australia, bred them, and raised their progeny in standard conditions. Evolved shifts in immune function should manifest as differences in immune responses among the progeny of parents collected in different locations. Parental location did not affect the offspring's cell-mediated immune response or stress response, but blood from the offspring of invasion-front toads had more neutrophils, and was more effective at phagocytosis and killing bacteria. These latter measures of immune function are negatively correlated with rate of dispersal in free-ranging toads. Our results suggest that the invasion of tropical Australia by cane toads has resulted in rapid genetically based compensatory shifts in the aspects of immune responses that are most compromised by the rigours of long-distance dispersal.

Phylogenetic relatedness and proboscis length contribute to structuring bumblebee communities in the extremes of abiotic and biotic gradients

Aim To improve our understanding of how biological communities assemble, we investigated changes in bumblebee communities in space along an elevation gradient. We assessed how much deterministic abiotic and biotic factors shape community assembly. We focused on proboscis length (influencing the species' dietary regime) and phylogenetic relatedness to investigate if competition and environmental filtering occur in more and less productive climates, respectively. Location Western Swiss Alps. Methods We recorded bumblebee species in 149 plots along a 1800-m wide elevation gradient. We contrasted two major clades of bumblebees, a short-tongued and a long-tongued clade. We calculated the phylogenetic and proboscis-length diversity of the bumblebee communities and compared these observed data with a random distribution to detect clustering likely to be caused by environmental filtering or overdispersion likely to be caused by competition. We compared the prevalence of clustered and overdispersed communities along the gradients of plant species richness (biotic) and temperature (abiotic). Results Under colder conditions, where plant species richness is lower and floral resources are scarcer, the clade with shorter proboscides prevails over the clade with longer proboscides, and communities are functionally and phylogenetic clustered. Under warmer conditions, we found phylogenetic but not functional overdispersion in communities. Main conclusions We show for the first time a strong correlation between phylogenetic relatedness, proboscis length and species distribution along temperature and plant richness gradients shaping bumblebee communities. The low temperatures and low levels of plant species richness limit the dispersal of the species from the long-tongued clade, which have more specialized diets, into high-elevation areas. Competition under warmer conditions may produce communities composed of less closely related species that share distinct ecological preferences. Our empirical results corroborate theoretical expectation as well as experiments on the prevalence of deterministic processes in the most severe and most productive parts of environmental gradients.

Climate oscillations and species interactions: large-scale congruence but regional differences in the phylogeographic structures of an alpine plant and its monophagous insect

Aim. To predict the fate of alpine interactions involving specialized species, using a monophagous beetle and its host-plant as a case study. Location. The Alps. Methods. We investigated genetic structuring of the herbivorous beetle Oreina gloriosa and its specific host-plant Peucedanum ostruthium. We used genome fingerprinting (in the insect and the plant) and sequence data (in the insect) to compare the distribution of the main gene pools in the two associated species and to estimate divergence time in the insect, a proxy for the temporal origin of the interaction. We quantified the similarity in spatial genetic structures by performing a Procrustes analysis, a tool from the shape theory. Finally, we simulated recolonization of an empty space analogous to the deglaciated Alps just after ice retreat by two lineages from two species showing unbalanced dependence, to examine how timing of the recolonization process, as well as dispersal capacities of associated species, could explain the observed pattern. Results. Contrasting with expectations based on their asymmetrical dependence, patterns in the beetle and plant were congruent at a large scale. Exceptions occurred at a regional scale in areas of admixture, matching known suture zones in Alpine plants. Simulations using a lattice-based model suggested these empirical patterns arose during or soon after recolonization, long after the estimated origin of the interaction c. 0.5 million years ago. Main conclusions. Species-specific interactions are scarce in alpine habitats because glacial cycles have limited opportunities for coevolution. Their fate, however, remains uncertain under climate change. Here we show that whereas most dispersal routes are paralleled at large scale, regional incongruence implies that the destinies of the species might differ under changing climate. This may be a consequence of the host-dependence of the beetle that locally limits the establishment of dispersing insects.

Ecology and evolution of social-organization: insights from fire ants and other highly eusocial insects

Social organisms exhibit conspicuous intraspecific variation in all facets of their social organization. A prominent example of such variation in the highly eusocial Hymenoptera is differences in the number of reproductive queens per colony, Differences in queen number in ants are associated with differences in a host of reproductive and social traits, including queen phenotype and breeding strategy, mode of colony reproduction, and pattern of sex allocation. We examine the causes and consequences of changes in colony queen number and associated traits using the fire ant Solenopsis invicta as a principal model. Ecological constraints on mode of colony founding may act as important selective forces causing the evolution of queen number in this and many other ants, with social organization generally perpetuated across generations by means of the social environment molding appropriate queen phenotypes and reproductive strategies. Shifts in colony queen number have profound effects on genetic structure within nests and may also influence genetic structure at higher levels (aggregations of nests or local demes) because of the association of queen number with particular mating and dispersal habits. Divergence of breeding habits between populations with different social organizations has the potential to promote genetic differentiation between these social variants. Thus, evolution of social organization can be important in generating intrinsic selective regimes that channel subsequent social evolution and in initiating the development of significant population genetic structure, including barriers to gene flow important in cladogenesis.