Where will conflicts between alien and rare species occur after climate and land-use change? A test with a novel combined modelling approach

Protecting native biodiversity against alien invasive species requires powerful methods to anticipate these invasions and to protect native species assumed to be at risk. Here, we describe how species distribution models (SDMs) can be used to identify areas predicted as suitable for rare native species and also predicted as highly susceptible to invasion by alien species, at present and under future climate and land-use scenarios. To assess the condition and dynamics of such conflicts, we developed a combined predictive modelling (CPM) approach, which predicts species distributions by combining two SDMs fitted using subsets of predictors classified as acting at either regional or local scales. We illustrate the CPM approach for an alien invader and a rare species associated to similar habitats in northwest Portugal. Combined models predict a wider variety of potential species responses, providing more informative projections of species distributions and future dynamics than traditional, non-combined models. They also provide more informative insight regarding current and future rare-invasive conflict areas. For our studied species, conflict areas of highest conservation relevance are predicted to decrease over the next decade, supporting previous reports that some invasive species may contract their geographic range and impact due to climate change. More generally, our results highlight the more informative character of the combined approach to address practical issues in conservation and management programs, especially those aimed at mitigating the impact of invasive plants, land-use and climate changes in sensitive regions

Regional variation in sex ratios and sex allocation in androdioecious Mercurialis annua.

In androdioecious metapopulations, where males co-occur with hermaphrodites, the absence of males from certain populations or regions may be explained by locally high selfing rates, high hermaphrodite outcross siring success (e.g. due to high pollen production by hermaphrodites), or to stochastic processes (e.g. the failure of males to invade populations or regions following colonization or range expansion by hermaphrodites). In the Iberian Peninsula and Morocco, the presence of males with hermaphrodites in the wind-pollinated androdioecious plant Mercurialis annua (Euphorbiaceae) varies both among populations within relatively small regions and among regions, with some regions lacking males from all populations. The species is known to have expanded its range into the Iberian Peninsula from a southern refugium. To account for variation in male presence in M. annua, we test the following hypotheses: (1) that males are absent in areas where plant densities are lower, because selfing rates should be correspondingly higher; (2) that males are absent in areas where hermaphrodites produce more pollen; and (3) that males are absent in areas where there is an elevated proportion of populations in which plant density and hermaphrodite pollen production disfavour their invasion. We found support for predictions two and three in Morocco (the putative Pleistocene refugium for M. annua) but no support for any hypothesis in Iberia (the expanded range). Our results are partially consistent with a hypothesis of sex-allocation equilibrium for populations in Morocco; in Iberia, the absence of males from large geographical regions is more consistent with a model of sex-ratio evolution in a metapopulation with recurrent population turnover. Our study points to the role of both frequency-dependent selection and contingencies imposed by colonization during range expansions and in metapopulations.

Will climate change drive alien invasive plants into areas of high protection value? An improved model-based regional assessment to prioritise the management of invasions.

Species distribution models (SDMs) studies suggest that, without control measures, the distribution of many alien invasive plant species (AIS) will increase under climate and land-use changes. Due to limited resources and large areas colonised by invaders, management and monitoring resources must be prioritised. Choices depend on the conservation value of the invaded areas and can be guided by SDM predictions. Here, we use a hierarchical SDM framework, complemented by connectivity analysis of AIS distributions, to evaluate current and future conflicts between AIS and high conservation value areas. We illustrate the framework with three Australian wattle (Acacia) species and patterns of conservation value in Northern Portugal. Results show that protected areas will likely suffer higher pressure from all three Acacia species under future climatic conditions. Due to this higher predicted conflict in protected areas, management might be prioritised for Acacia dealbata and Acacia melanoxylon. Connectivity of AIS suitable areas inside protected areas is currently lower than across the full study area, but this would change under future environmental conditions. Coupled SDM and connectivity analysis can support resource prioritisation for anticipation and monitoring of AIS impacts. However, further tests of this framework over a wide range of regions and organisms are still required before wide application.

Soil fungal communities of grasslands are environmentally structured at a regional scale in the Alps.

Studying patterns of species distributions along elevation gradients is frequently used to identify the primary factors that determine the distribution, diversity and assembly of species. However, despite their crucial role in ecosystem functioning, our understanding of the distribution of below-ground fungi is still limited, calling for more comprehensive studies of fungal biogeography along environmental gradients at various scales (from regional to global). Here, we investigated the richness of taxa of soil fungi and their phylogenetic diversity across a wide range of grassland types along a 2800 m elevation gradient at a large number of sites (213), stratified across a region of the Western Swiss Alps (700 km(2)). We used 454 pyrosequencing to obtain fungal sequences that were clustered into operational taxonomic units (OTUs). The OTU diversity-area relationship revealed uneven distribution of fungal taxa across the study area (i.e. not all taxa are everywhere) and fine-scale spatial clustering. Fungal richness and phylogenetic diversity were found to be higher in lower temperatures and higher moisture conditions. Climatic and soil characteristics as well as plant community composition were related to OTU alpha, beta and phylogenetic diversity, with distinct fungal lineages suggesting distinct ecological tolerances. Soil fungi, thus, show lineage-specific biogeographic patterns, even at a regional scale, and follow environmental determinism, mediated by interactions with plants.

Regional pattern of the molecular types of Cryptococcus neoformans and Cryptococcus gattii in Brazil

The molecular types of 443 Brazilian isolates of Cryptococcus neoformans and Cryptococcus gattii were analyzed to determine their geographic distribution within Brazil and their underlying host conditions. The following data, imported from previous epidemiological studies as well as two culture collections, were analyzed for: place of isolation, source (clinical or environmental), host risk factors, species, serotype, mating type, and molecular type. Molecular typing by PCR-fingerprinting using primers for the minisatellite-specific core sequence of the wild-type phage M13 or microsatellites [(GACA)4, (GTG)5], restriction fragment length polymorphism of URA5 gene analysis, and/or amplified fragment length polymorphism (AFLP) identified eight major genotypes: VNI/AFLP1, VNII/AFLP1A, VNIII/AFLP2, and VNIV/AFLP3 for C. neoformans, and VGI/AFLP4, VGII/AFLP6, VGIII/AFLP5, and VGIV/AFLP7 for C. gattii. The most common molecular type found in Brazil was VNI (64%), followed by VGII (21%), VNII (5%), VGIII (4%), VGI and VNIV (3% each), and VNIII (< 1%). Primary cryptococcosis caused by the molecular type VGII (serotype B, MAT) prevails in immunocompetent hosts in the North and Northeast regions, disclosing an endemic regional pattern for this specific molecular type in the Northern Brazil.

The molecular detection of different Leishmania species within sand flies from a cutaneous and visceral leishmaniasis sympatric area in Southeastern Brazil

Over the last 20 years, there has been an increase in the number of leishmaniasis cases in Brazil. Belo Horizonte (BH) is one of the most highly populated Brazilian cities that is affected by visceral leishmaniasis (VL). The health services in BH are coordinated by a central nucleus that is subdivided into nine sanitary districts. Historically, the highest level of human VL cases was found in the northeast sanitary district (NSD). The objective of our study was to detect Leishmania infection in the phlebotomine sand flies collected in the NSD by dissection and molecular approaches. Following the occurrence of human VL cases in 2005, entomological captures were performed from July 2006-June 2007. Out of the 245 sand flies dissected, only three Lutzomyia longipalpis spp contained flagellates. The female sand flies were grouped into 120 pools according to date, collection site and species, with approximately 10 individual sand flies in each pool. Subsquently, the DNA was extracted and Leishmania spp and other parasites were detected and identified by polymerase chain reaction (PCR) and PCR-restriction fragment length polymorfism. Leishmania infantum was present in at least 19% of the Lu. longipalpis collected, in 3.8% of the Nyssomiya whitmani collected, in 33.3% of the Evandromiya termitophila collected and in 14.3% of the Nyssomiya intermedia collected. When the females of the cortelezzii complex were compared with each other, 3.2% of the females were infected with Leishmania braziliensis, whereas 3.2% of the females were infected with trypanosomatids.

Malaria vector species in Colombia: a review

Here we present a comprehensive review of the literature on the vectorial importance of the major Anopheles malaria vectors in Colombia. We provide basic information on the geographical distribution, altitudinal range, immature habitats, adult behaviour, feeding preferences and anthropophily, endophily and infectivity rates. We additionally review information on the life cycle, longevity and population fluctuation of Colombian Anopheles species. Emphasis was placed on the primary vectors that have been epidemiologically incriminated in malaria transmission: Anopheles darlingi, Anopheles albimanus and Anopheles nuneztovari. The role of a selection of local, regional or secondary vectors (e.g., Anopheles pseudopunctipennis and Anopheles neivai) is also discussed. We highlight the importance of combining biological, morphological and molecular data for the correct taxonomical determination of a given species, particularly for members of the species complexes. We likewise emphasise the importance of studying the bionomics of primary and secondary vectors along with an examination of the local conditions affecting the transmission of malaria. The presence and spread of the major vectors and the emergence of secondary species capable of transmitting human Plasmodia are of great interest. When selecting control measures, the anopheline diversity in the region must be considered. Variation in macroclimate conditions over a species' geographical range must be well understood and targeted to plan effective control measures based on the population dynamics of the local Anopheles species.

First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species, Pinus halepensis Mill.

Aleppo pine (Pinus halepensis Mill.) is a relevant conifer species for studying adaptive responses to drought and fire regimes in the Mediterranean region. In this study, we performed Illumina next-generation sequencing of two phenotypically divergent Aleppo pine accessions with the aims of (i) characterizing the transcriptome through Illumina RNA-Seq on trees phenotypically divergent for adaptive traits linked to fire adaptation and drought, (ii) performing a functional annotation of the assembled transcriptome, (iii) identifying genes with accelerated evolutionary rates, (iv) studying the expression levels of the annotated genes and (v) developing gene-based markers for population genomic and association genetic studies. The assembled transcriptome consisted of 48,629 contigs and covered about 54.6 Mbp. The comparison of Aleppo pine transcripts to Picea sitchensis protein-coding sequences resulted in the detection of 34,014 SNPs across species, with a Ka /Ks average value of 0.216, suggesting that the majority of the assembled genes are under negative selection. Several genes were differentially expressed across the two pine accessions with contrasted phenotypes, including a glutathione-s-transferase, a cellulose synthase and a cobra-like protein. A large number of new markers (3334 amplifiable SSRs and 28,236 SNPs) have been identified which should facilitate future population genomics and association genetics in this species. A 384-SNP Oligo Pool Assay for genotyping with the Illumina VeraCode technology has been designed which showed an high overall SNP conversion rate (76.6%). Our results showed that Illumina next-generation sequencing is a valuable technology to obtain an extensive overview on whole transcriptomes of nonmodel species with large genomes.

Overcoming the rare species modelling paradox: test of a novel hierarchical framework with an Iberian endemic plant

Rare species have restricted geographic ranges, habitat specialization, and/or small population sizes. Datasets on rare species distribution usually have few observations, limited spatial accuracy and lack of valid absences; conversely they provide comprehensive views of species distributions allowing to realistically capture most of their realized environmental niche. Rare species are the most in need of predictive distribution modelling but also the most difficult to model. We refer to this contrast as the "rare species modelling paradox" and propose as a solution developing modelling approaches that deal with a sufficiently large set of predictors, ensuring that statistical models aren't overfitted. Our novel approach fulfils this condition by fitting a large number of bivariate models and averaging them with a weighted ensemble approach. We further propose that this ensemble forecasting is conducted within a hierarchic multi-scale framework. We present two ensemble models for a test species, one at regional and one at local scale, each based on the combination of 630 models. In both cases, we obtained excellent spatial projections, unusual when modelling rare species. Model results highlight, from a statistically sound approach, the effects of multiple drivers in a same modelling framework and at two distinct scales. From this added information, regional models can support accurate forecasts of range dynamics under climate change scenarios, whereas local models allow the assessment of isolated or synergistic impacts of changes in multiple predictors. This novel framework provides a baseline for adaptive conservation, management and monitoring of rare species at distinct spatial and temporal scales.

Loss of interactions with ants under cold climate in a regional myrmecophilous butterfly fauna

Aim Specialized mutualistic clades may revert and thus increase their autonomy and generalist characteristics. However, our understanding of the drivers that trigger reductions in mutualistic traits and of the consequences for the tolerance of these species to various environmental conditions remains limited. This study investigates the relationship between the environmental niche and the degree of myrmecophily (i.e. the ability to interact with ants) among members of the Lycaenidae. Location The western Swiss Alps. Methods We measured the tolerance of Lycaenidae species to low temperatures by comparing observations from a random stratified field sampling with climatic maps. We then compared the species-specific degree of myrmecophily with the species range limits at colder temperatures while controlling for phylogenetic dependence. We further evaluated whether the community-averaged degree of myrmecophily increases with temperature, as would be expected in the case of environmental filters acting on myrmecophilous species. Results Twenty-nine Lycaenidae species were found during sampling. Ancestral state reconstruction indicated that the 24 species of Polyommatinae displayed both strong myrmecophily and secondary loss of mutualism; these species were used in the subsequent statistical analyses. Species with a higher degree of ant interaction were, on average, more likely to inhabit warmer sites. Species inhabiting the coldest environments displayed little or no interaction with ants. Main conclusions Colder climates at high elevations filter out species with a high degree of myrmecophily and may have been the direct evolutionary force that promoted the loss of mutualism. A larger taxon sampling across the Holarctic may help to distinguish between the ecological and evolutionary effects of climate.

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.

Genetic differentiation in Silene dioica metapopulations: estimation of spatiotemporal effects in a successional plant species

Silene dioica is a diploid, dioecious, perennial, insect-pollinated herb and part of the deciduous phase of primary succession in Skeppsvik Archipelago, Gulf of Bothnia, Sweden. These islands are composed of material deposited and left underwater by melting ice at the end of the last ice age. A rapid and relatively constant rate of land uplift of 0.9 cm per year continually creates new islands available for colonization by plants. Because the higher deposits appear first, islands differ in age. Because it is possible to estimate the ages of islands and populations of plant species belonging to early stages of succession, the genetic dynamics occurring within an age-structured metapopulation can be investigated in this archipelago. Fifty-two island populations of S. dioica of known ages, sizes, and distances from each other were studied through electrophoretic data. A number of factors increase the degree of genetic differentiation among these island populations relative to an island model at equilibrium. Newly founded populations were more differentiated than those of intermediate age, which suggests that colonization dynamics increase genetic variance among populations. The very old populations, which decrease in size as they approach extinction, were more differentiated than intermediate-aged populations. Isolation by distance occurs in this system. Colonizers are likely to come from more than one source, and the migrant pool model best explains colonization events in the archipelago. Degree of environmental exposure also affects population differentiation.

Complementarity between human capital and trade in regional technological progress

The effect of openness and trade orientation on economic growth remains a highly contentious issue in the literature. Trade facilitates the spread of knowledge and the adoption of more advanced and efficient technologies, which hastens total factor productivity (TFP) growth and, hence, per capita income. New technologies that spread through trade require a sufficiently skilled labour force to adapt them to the domestic productive environment. Thus, openness and human capital accumulation will lead to TFP growth and the greater the complementarity between both variables, the higher the TFP growth. This paper discusses the implications of these assumptions and tests their empirical validity, using a pool of data for manufacturing industry in Spanish regions in a period in which both the stock of human capital and openness experienced a notable increase.

Plant species richness and environmental heterogeneity in a mountain landscape: effects of variability and spatial configuration

The loss of biodiversity has become a matter of urgent concern and a better understanding of local drivers is crucial for conservation. Although environmental heterogeneity is recognized as an important determinant of biodiversity, this has rarely been tested using field data at management scale. We propose and provide evidence for the simple hypothesis that local species diversity is related to spatial environmental heterogeneity. Species partition the environment into habitats. Biodiversity is therefore expected to be influenced by two aspects of spatial heterogeneity: 1) the variability of environmental conditions, which will affect the number of types of habitat, and 2) the spatial configuration of habitats, which will affect the rates of ecological processes, such as dispersal or competition. Earlier, simulation experiments predicted that both aspects of heterogeneity will influence plant species richness at a particular site. For the first time, these predictions were tested for plant communities using field data, which we collected in a wooded pasture in the Swiss Jura mountains using a four-level hierarchical sampling design. Richness generally increased with increasing environmental variability and "roughness" (i.e. decreasing spatial aggregation). Effects occurred at all scales, but the nature of the effect changed with scale, suggesting a change in the underlying mechanisms, which will need to be taken into account if scaling up to larger landscapes. Although we found significant effects of environmental heterogeneity, other factors such as history could also be important determinants. If a relationship between environmental heterogeneity and species richness can be shown to be general, recently available high-resolution environmental data can be used to complement the assessment of patterns of local richness and improve the prediction of the effects of land use change based on mean site conditions or land use history.

Complementarity between human capital and trade in regional technological progress

The effect of openness and trade orientation on economic growth remains a highly contentious issue in the literature. Trade facilitates the spread of knowledge and the adoption of more advanced and efficient technologies, which hastens total factor productivity (TFP) growth and, hence, per capita income. New technologies that spread through trade require a sufficiently skilled labour force to adapt them to the domestic productive environment. Thus, openness and human capital accumulation will lead to TFP growth and the greater the complementarity between both variables, the higher the TFP growth. This paper discusses the implications of these assumptions and tests their empirical validity, using a pool of data for manufacturing industry in Spanish regions in a period in which both the stock of human capital and openness experienced a notable increase.