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.

Simultaneous identification of multiple mammalian species from mixed forensic samples based on mtDNA control region length polymorphism

Molecular species identification in mixed or contaminated biological material has always been problematic. We developed a simple and accurate method for mammal DNA identification in mixtures, based on interspecific mitochondrial DNA control region length polymorphism. Contrary to other published methods dealing with species mixtures, our protocol requires a single universal primer pair and amplification step, and is not based on a pre-defined panel of species. This protocol has been routinely employed by our laboratory for species identification in dozens of human and animal forensic caseworks. Six representative forensic caseworks involving the specific identification of mixed animal samples are reported in this paper, in order to demonstrate the applicability and usefulness of the method.

Diversity and efficiency of bradyrhizobium strains isolated from soil samples collected from around sesbania virgata roots using cowpea as trap species

The genetic diversity of ten Bradyrhizobium strains was evaluated for tolerance to high temperatures, to different salinity levels and for the efficiency of symbiosis with cowpea plants (Vigna unguiculata (L.) Walp.). Eight of these strains were isolated from nodules that appeared on cowpea after inoculation with suspensions of soil sampled from around the root system of Sesbania virgata (wand riverhemp) in ecosystems of South Minas Gerais. The other two strains used in our analyses as references, were from the Amazon and are currently recommended as cowpea inoculants. Genetic diversity was analyzed by amplifying repetitive DNA elements with the BOX primer, revealing high genetic diversity with each strain presenting a unique band profile. Leonard jar assays showed that the strains UFLA 03-30 and UFLA 03-38 had the highest N2-fixing potentials in symbiosis with cowpea. These strains had more shoot and nodule dry matter, more shoot N accumulation, and a higher relative efficiency than the strains recommended as inoculants. All strains grew in media of pH levels ranging from 4.0 to 9.0. The strains with the highest N2-fixing efficiencies in symbiosis with cowpea were also tolerant to the greatest number of antibiotics. However, these strains also had the lowest tolerance to high salt concentrations. All strains, with the exceptions of UFLA 03-84 and UFLA 03-37, tolerated temperatures of up to 40 ºC. The genetic and phenotypic characteristics of the eight strains isolated from soils of the same region were highly variable, as well as their symbiotic efficiencies, despite their common origin. This variability highlights the importance of including these tests in the selection of cowpea inoculant strains.

Evidence of climatic niche shift during biological invasion

Niche-based models calibrated in the native range by relating species observations to climatic variables are commonly used to predict the potential spatial extent of species' invasion. This climate matching approach relies on the assumption that invasive species conserve their climatic niche in the invaded ranges. We test this assumption by analysing the climatic niche spaces of Spotted Knapweed in western North America and Europe. We show with robust cross-continental data that a shift of the observed climatic niche occurred between native and non-native ranges, providing the first empirical evidence that an invasive species can occupy climatically distinct niche spaces following its introduction into a new area. The models fail to predict the current invaded distribution, but correctly predict areas of introduction. Climate matching is thus a useful approach to identify areas at risk of introduction and establishment of newly or not-yet-introduced neophytes, but may not predict the full extent of invasions.

Unifying niche shift studies: insights from biological invasions.

Assessing whether the climatic niche of a species may change between different geographic areas or time periods has become increasingly important in the context of ongoing global change. However, approaches and findings have remained largely controversial so far, calling for a unification of methods. Here, we build on a review of empirical studies of invasion to formalize a unifying framework that decomposes niche change into unfilling, stability, and expansion situations, taking both a pooled range and range-specific perspective on the niche, while accounting for climatic availability and climatic analogy. This framework provides new insights into the nature of climate niche shifts and our ability to anticipate invasions, and may help in guiding the design of experiments for assessing causes of niche changes.

Placentation in species of phylogenetic importance: the Afrotheria.

Afrotheria, one of four mammalian superorders, comprises elephants, sea cows, hyraxes, aardvark, elephant shrews, tenrecs and golden moles. Their placentas either form an equatorial band or are discoid in shape. The interhemal region, separating fetal and maternal blood, is endotheliochorial in elephants, aardvark and possibly the sea cows, but hemochorial in the remaining orders. There is a secondary epitheliochorial placenta in elephant shrews while a similar structure in tenrecs erodes maternal tissues. Specialized hemophagous regions are a striking characteristic of some of these placentas yet absent in hyraxes, elephant shrews, and golden moles. It is possible that the common ancestor of the Afrotheria had an endotheliochorial placenta. Establishment of a hemochorial condition, as seen in rock hyraxes, elephant shrews, tenrecs, and golden moles, would be a more recent development. The elephant, manatee, and aardvark all have circumferential placentas. Thus the formation of a discoid placenta with a more or less extensive secondary placenta in elephant shrews and tenrecs would also be a derived state.

Characterization of Exopolysaccharides Produced by Rhizobia Species

ABSTRACT Increasing attention has been given, over the past decades, to the production of exopolysaccharides (EPS) from rhizobia, due to their various biotechnological applications. Overall characterization of biopolymers involves evaluation of their chemical, physical, and biological properties; this evaluation is a key factor in understanding their behavior in different environments, which enables researchers to foresee their potential applications. Our focus was to study the EPS produced by Mesorhizobium huakuii LMG14107, M. loti LMG6125, M. plurifarium LMG11892,Rhizobium giardini bv. giardiniH152T, R. mongolense LMG19141, andSinorhizobium (= Ensifer)kostiense LMG19227 in a RDM medium with glycerol as a carbon source. These biopolymers were isolated and characterized by reversed-phase high-performance liquid chromatography (RP-HPLC), Fourier transform infrared (FTIR), and nuclear magnetic resonance (NMR) spectroscopies. Maximum exopolysaccharide production was 3.10, 2.72, and 2.50 g L-1for the strains LMG6125, LMG19227, and LMG19141, respectively. The purified EPS revealed prominent functional reactive groups, such as hydroxyl and carboxylic, which correspond to a typical heteropolysaccharide. The EPS are composed primarily of galactose and glucose. Minor components found were rhamnose, glucuronic acid, and galacturonic acid. Indeed, from the results of techniques applied in this study, it can be noted that the EPS are species-specific heteropolysaccharide polymers composed of common sugars that are substituted by non-carbohydrate moieties. In addition, analysis of these results indicates that rhizobial EPS can be classified into five groups based on ester type, as determined from the 13C NMR spectra. Knowledge of the EPS composition now facilitates further investigations relating polysaccharide structure and dynamics to rheological properties.

PPAR alpha structure-function relationships derived from species-specific differences in responsiveness to hypolipidemic agents.

The nuclear receptor PPAR alpha is a key regulatory transcription factor in lipid homeostasis, some liver detoxification processes and the control of inflammation. Recent findings suggest that many hypolipidemic drugs and anti-inflammatory agents can potentially act by binding to PPAR alpha and inducing its activity. Here, we identify some structure-function relationships in PPAR alpha, by using the species-specific responsiveness to the two hypolipidemic agents, Wy 14,643 and 5,8,11,14-eicosatetraynoic acid (ETYA). We first show that the species-specific differences are mediated primarily via the ligand binding domain of the receptor and that these two drugs are indeed ligands of PPAR alpha. By mutagenesis analyses we identify amino acid residues in the ligand binding domains of Xenopus, mouse and human PPAR alpha, that confer preferential responsiveness to ETYA and Wy 14,643. These findings will aid in the development of new synthetic PPAR alpha ligands as effective therapeutics for lipid-related diseases and inflammatory disorders.

What drives invasibility? A multi-model inference test and spatial modelling of alien plant species richness patterns in northern Portugal

Understanding and anticipating biological invasions can focus either on traits that favour species invasiveness or on features of the receiving communities, habitats or landscapes that promote their invasibility. Here, we address invasibility at the regional scale, testing whether some habitats and landscapes are more invasible than others by fitting models that relate alien plant species richness to various environmental predictors. We use a multi-model information-theoretic approach to assess invasibility by modelling spatial and ecological patterns of alien invasion in landscape mosaics and testing competing hypotheses of environmental factors that may control invasibility. Because invasibility may be mediated by particular characteristics of invasiveness, we classified alien species according to their C-S-R plant strategies. We illustrate this approach with a set of 86 alien species in Northern Portugal. We first focus on predictors influencing species richness and expressing invasibility and then evaluate whether distinct plant strategies respond to the same or different groups of environmental predictors. We confirmed climate as a primary determinant of alien invasions and as a primary environmental gradient determining landscape invasibility. The effects of secondary gradients were detected only when the area was sub-sampled according to predictions based on the primary gradient. Then, multiple predictor types influenced patterns of alien species richness, with some types (landscape composition, topography and fire regime) prevailing over others. Alien species richness responded most strongly to extreme land management regimes, suggesting that intermediate disturbance induces biotic resistance by favouring native species richness. Land-use intensification facilitated alien invasion, whereas conservation areas hosted few invaders, highlighting the importance of ecosystem stability in preventing invasions. Plants with different strategies exhibited different responses to environmental gradients, particularly when the variations of the primary gradient were narrowed by sub-sampling. Such differential responses of plant strategies suggest using distinct control and eradication approaches for different areas and alien plant groups.

Parallel evolution of facial stripe patterns in the Neolamprologus brichardi/pulcher species complex endemic to Lake Tanganyika.

Colour pattern diversity can be due to random processes or to natural or sexual selection. Consequently, similarities in colour patterns are not always correlated with common ancestry, but may result from convergent evolution under shared selection pressures or drift. Neolamprologus brichardi and Neolamprologus pulcher have been described as two distinct species based on differences in the arrangement of two dark bars on the operculum. Our study uses DNA sequences of the mitochondrial control region to show that relatedness of haplotypes disagrees with species assignment based on head colour pattern. This suggests repeated parallel evolution of particular stripe patterns. The complete lack of shared haplotypes between populations of the same or different phenotypes reflects strong philopatric behaviour, possibly induced by the cooperative breeding mode in which offspring remain in their natal territory and serve as helpers until they disperse to nearby territories or take over a breeding position. Concordant phylogeographic patterns between N. brichardi/N. pulcher populations and other rock-dwelling cichlids suggest that the same colonization routes have been taken by sympatric species and that these routes were affected by lake level fluctuations in the past.

On the success of generalized food deception in orchids : the role of sympatric rewarding species and pollinators

Résumé : La production de nectar assure aux plantes entomophiles un important succès reproducteur. Malgré cela, de nombreuses espèces d'orchidées ne produisent pas de nectar. La majorité de ces orchidées dites trompeuses exploitent simplement l'instinct des pollinisateurs généralistes, qui les pousse à chercher du nectar dans les fleurs. Afin d'optimiser la récolte de nectar, les pollinisateurs apprennent à différencier les fleurs trompeuses des nectarifères, et à concentrer leurs visites sur ces dernières, au détriment des plantes trompeuses. Chez les orchidées non autogames, la reproduction est assurée uniquement par les pollinisateurs. L'apprentissage des pollinisateurs a donc un impact négatif sur la reproduction des orchidées trompeuses. Cependant, les caractéristiques d'une espèce trompeuse et des espèces nectarifères au sein d'une communauté végétale peuvent affecter l'apprentissage et le taux de visite des pollinisateurs aux plantes trompeuses. J'ai réalisé des expériences en milieu naturel et en milieu contrôlé, pour déterminer si les caractéristiques florales, spatiales et temporelles des communautés affectent le taux de visite et le succès reproducteur de plantes trompeuses. Une agrégation spatiale élevée des plantes trompeuses et des plantes nectarifères diminue le succès reproducteur des plantes trompeuses. De plus, les pollinisateurs visitent plus souvent l'espèce trompeuse Iorsque ses fleurs sont de couleur similaire à celles de l'espèce nectarifère. Cet effet bénéfique de la similarité pour la couleur des fleurs s'accentue si les deux espèces sont mélangées et proches spatialement, ou si l'espèce trompeuse fleurit après l'espèce nectarifère. Enfin, le comportement des pollinisateurs n'est pas tout de suite affecté lorsque les caractéristiques de la communauté changent. Les caractéristiques des communautés végétales affectent donc la reproduction des espèces trompeuses. Bien que L'absence de coûts associés à la production de nectar, l'exportation efficace de pollen et la production de graines de qualité dont bénéficient les orchidées trompeuses favorisent Ieur maintien, les caractéristiques de la communauté peuvent aussi y contribuer. Mon étude fournit donc une explication alternative et complémentaire au maintien des orchidées trompeuses. Je conclus par une discussion des implications possibles de ces résultats sur le maintien et l'évolution des orchidées trompeuses, en tenant compte de la dynamique des caractéristiques des communautés végétales naturelles. Abstract : Despite the importance of producing food to ensure a high reproductive success, many orchid species lack such rewards. The majority of deceptive orchids simply exploit the instinctive food-foraging behaviour of generalist pollinators. This strategy is termed generalized food deception. To optimize their foraging efficiency, pollinators can learn to discriminate deceptive from rewarding flowers and to focus their visits to the rewarding plants, to the disadvantage of the deceptive plants. Because the reproductive success of non-autogamous orchids entirely relies on pollinator visitation rate, pollinator learning decreases the reproductive success of deceptive orchids. However, the characteristics of deceptive and rewarding plants within a community may affect pollinator learning and visitation rate to a deceptive orchid. Therefore, the biological characteristics of natural plant communities may be crucial to the maintenance of generalized food deceptive orchids. My study focused on the floral, spatial and temporal characteristics of plant communities. I used both in and ex sitar experiments to investigate whether these characteristics influence pollinator visitation rates and the reproductive success of deceptive orchids. A high spatial aggregation of both deceptive and rewarding species decreased the reproductive success of the deceptive species. Also, being of similar flower colour to rewarding sympatric species increased pollinator visitation rates to a deceptive species. The beneficial effect of flower colour similarity was even more pronounced when both species were spatially closely mingled or when the deceptive species flowered after the rewarding species. Finally, pollinator behaviour was unaffected in the short term by a change in the characteristics of plant communities, indicating that pollinators need time to learn under new conditions. Thus, the characteristics of plant communities may crucially affect the reproductive success of deceptive orchids. Although the absence of costs associated with nectar production, the efficient pollen export and the high seed quality of deceptive orchids may favour their maintenance, the characteristics of plant communities may also contribute to it. Therefore, my study provides an alternative yet complementary explanation to the maintenance of generalized food deceptive orchids in natural populations. I discuss the possible implications for the maintenance and the evolution of generalized food deceptive orchids with regards to the floral and temporal dynamics of natural plant communities.

Inter and intra-guild interactions in egg parasitoid species of the soybean stink bug complex

The objective of this research was to evaluate the parasitism behavior of Telenomus podisi Ashmead, Trissolcus basalis (Wollaston) e Trissolcus urichi Crawford (Hymenoptera: Scelionidae) on eggs of Nezara viridula L., Euschistus heros F., Piezodorus guildinii Westwood and Acrosternum aseadum Rolston (Heteroptera: Pentatomidae), in no choice and multiple choice experiments. For all parasitoid species, the results demonstrated the existence of a main host species that maximizes the reproductive success. The competitive interactions among the parasitoid species were investigated in experiments of sequential and simultaneous release of different combinations of parasitoid pairs on the hosts N. viridula, E. heros and A. aseadum. Exploitative competition was observed for egg batches at the genus level (Telenomus vs. Trissolcus) and interference competition at the species level (T. basalis vs. T. urichi). Trissolcus urichi was the most aggressive species, interfering with the parasitism of T. basalis. Generally, T. basalis showed an opportunistic behavior trying to parasitise eggs after T. urichi had abandoned the egg batch. The selection of parasitoid species for use in augmentative biological control programs should take into account the diversity of pentatomids present in soybean in addition to the interactions among the different species of parasitoids.

Residence time, expansion toward the equator in the invaded range and native range size matter to climatic niche shifts in non-native species

Aim Identifying climatic niche shifts and their drivers is important to accurately predict the risk of biological invasions. The niches of non-native plants and birds have recently been assessed in large-scale multi-species studies, but such large-scale tests are lacking for non-native reptiles and amphibians (herpetofauna). Furthermore, little is known about the factors contributing to niche shifts when they occur. Based on the occurrence of 71 reptile and amphibian species, we compared native and non-native realized niches in 101 invaded ranges at a worldwide scale and identified the factors that affect niche shifts. Location The world except the Antarctic. Methods We assessed climatic niche dynamics in a gridded environmental space allowing the quantification of niche overlap and expansion into climatic conditions not colonized by the species in their native range. We analyzed the factors affecting niche shifts using a model averaging approach based on generalized linear mixed-effects models. Results Approximately 57% of the invaded ranges (51% for amphibians and 61% for reptiles) showed niche shifts (≥10% expansion in the realized climatic niche). Island endemics, species introduced to Oceania and invaded ranges outside the native biogeographic realm showed a higher proportion of niche shifts. Niche shifts were more likely for species that had smaller native range sizes, were introduced earlier into a new range or invaded areas located at lower latitudes than the native range. Main conclusions The proportion of niche shifts for non-native herpetofauna was higher than those for Holarctic non-native plants and European non-native birds. The 'climate matching hypothesis' should be used with caution for species shifting their niche because it could underestimate the risk of their establishment.

Using niche-based models to improve the sampling of rare species

Because data on rare species usually are sparse, it is important to have efficient ways to sample additional data. Traditional sampling approaches are of limited value for rare species because a very large proportion of randomly chosen sampling sites are unlikely to shelter the species. For these species, spatial predictions from niche-based distribution models can be used to stratify the sampling and increase sampling efficiency. New data sampled are then used to improve the initial model. Applying this approach repeatedly is an adaptive process that may allow increasing the number of new occurrences found. We illustrate the approach with a case study of a rare and endangered plant species in Switzerland and a simulation experiment. Our field survey confirmed that the method helps in the discovery of new populations of the target species in remote areas where the predicted habitat suitability is high. In our simulations the model-based approach provided a significant improvement (by a factor of 1.8 to 4 times, depending on the measure) over simple random sampling. In terms of cost this approach may save up to 70% of the time spent in the field.

Mites as biological tags of their hosts.

Movements and spatial distribution of host populations are expected to shape the genetic structure of their parasite populations. Comparing the genetic patterns of both interacting species may improve our understanding of their evolutionary history. Moreover, genetic analyses of parasites with horizontal transmission may serve as indicators of historical events or current demographic processes that are not apparent in the genetic signature of their hosts. Here, we compared mitochondrial variation in populations of the ectoparasitic mite Spinturnix myoti with the genetic pattern of its host, the Maghrebian bat Myotis punicus in North Africa and in the islands of Corsica and Sardinia. Mite mitochondrial differentiation among populations was correlated with both host mitochondrial and nuclear differentiation, suggesting spatial co-differentiation of the lineages of the two interacting species. Therefore our results suggest that parasite dispersal is exclusively mediated by host movements, with open water between landmasses as a main barrier for host and parasite dispersal. Surprisingly the unique presence of a continental European mite lineage in Corsica was inconsistent with host phylogeographical history and strongly suggests the former presence of European mouse-eared bats on this island. Parasites may thus act as biological tags to reveal the presence of their now locally extinct host.