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.

Mineralogical, petrographic and geochemical characterisation of white and coloured Iberian marbles in the context of the provenancing of some artefacts from Thamusida (Kenitra, Morocco)

A multi-analytical study has been carried out on a collection of white and coloured Iberian marbles. A total of 135 marble specimens were collected in Spain and Portugal from the Betic chain (Alhaurin de la Torre, Mijas, Macael), Ossa Morena (Alconera, Almaden de la Plata and Viana do Alentejo), and the Estremoz Anticline (Bencatel, Borba and Estremoz) areas. X-ray diffractometry and carbon and oxygen stable isotope analysis were carried out on these samples; 38 samples were also investigated by optical and scanning electron microscopy. The results provide a set of diagnostic parameters that allow discriminating the sampled marble quarries. The carbonate minerals composition is distinctive for the Mijas and Alhaurin de la Torre marbles; the isotopic analysis allows discriminating also between these two dolomitic marble quarries. The Ossa Morena and Estremoz Anticline marbles share a similar stable isotope composition; the accessory mineral content, the maximum grain size (MGS) and the fabric are particularly useful in the distinction between them. In the framework of archaeometric provenance studies on Thamusida (Kenitra, Morocco) Roman marble artefacts, a specific comparison between this new Iberian database and archaeological findings has been carried out. The hypothesis of commercial exchanges between the Iberian regions and Roman Morocco is supported by the results of the provenance study, which suggested the Almaden de la Plata and Mijas quarries as possible sources of raw materials for the production of archaeological artefacts.

Conservation of the Eurasian Lynx (Lynx lynx) in a fragmented landscape - habitat models, dispersal and potential distribution

Abstract: The expansion of a recovering population - whether re-introduced or spontaneously returning - is shaped by (i) biological (intrinsic) factors such as the land tenure system or dispersal, (ii) the distribution and availability of resources (e.g. prey), (iii) habitat and landscape features, and (iv) human attitudes and activities. In order to develop efficient conservation and recovery strategies, we need to understand all these factors and to predict the potential distribution and explore ways to reach it. An increased number of lynx in the north-western Swiss Alps in the nineties lead to a new controversy about the return of this cat. When the large carnivores were given legal protection in many European countries, most organizations and individuals promoting their protection did not foresee the consequences. Management plans describing how to handle conflicts with large predators are needed to find a balance between "overabundance" and extinction. Wildlife and conservation biologists need to evaluate the various threats confronting populations so that adequate management decisions can be taken. I developed a GIS probability model for the lynx, based on habitat information and radio-telemetry data from the Swiss Jura Mountains, in order to predict the potential distribution of the lynx in this mountain range, which is presently only partly occupied by lynx. Three of the 18 variables tested for each square kilometre describing land use, vegetation, and topography, qualified to predict the probability of lynx presence. The resulting map was evaluated with data from dispersing subadult lynx. Young lynx that were not able to establish home ranges in what was identified as good lynx habitat did not survive their first year of independence, whereas the only one that died in good lynx habitat was illegally killed. Radio-telemetry fixes are often used as input data to calibrate habitat models. Radio-telemetry is the only way to gather accurate and unbiased data on habitat use of elusive larger terrestrial mammals. However, it is time consuming and expensive, and can therefore only be applied in limited areas. Habitat models extrapolated over large areas can in turn be problematic, as habitat characteristics and availability may change from one area to the other. I analysed the predictive power of Ecological Niche Factor Analysis (ENFA) in Switzerland with the lynx as focal species. According to my results, the optimal sampling strategy to predict species distribution in an Alpine area lacking available data would be to pool presence cells from contrasted regions (Jura Mountains, Alps), whereas in regions with a low ecological variance (Jura Mountains), only local presence cells should be used for the calibration of the model. Dispersal influences the dynamics and persistence of populations, the distribution and abundance of species, and gives the communities and ecosystems their characteristic texture in space and time. Between 1988 and 2001, the spatio-temporal behaviour of subadult Eurasian lynx in two re-introduced populations in Switzerland was studied, based on 39 juvenile lynx of which 24 were radio-tagged to understand the factors influencing dispersal. Subadults become independent from their mothers at the age of 8-11 months. No sex bias neither in the dispersal rate nor in the distance moved was detected. Lynx are conservative dispersers, compared to bear and wolf, and settled within or close to known lynx occurrences. Dispersal distances reached in the high lynx density population - shorter than those reported in other Eurasian lynx studies - are limited by habitat restriction hindering connections with neighbouring metapopulations. I postulated that high lynx density would lead to an expansion of the population and validated my predictions with data from the north-western Swiss Alps where about 1995 a strong increase in lynx abundance took place. The general hypothesis that high population density will foster the expansion of the population was not confirmed. This has consequences for the re-introduction and recovery of carnivores in a fragmented landscape. To establish a strong source population in one place might not be an optimal strategy. Rather, population nuclei should be founded in several neighbouring patches. Exchange between established neighbouring subpopulations will later on take place, as adult lynx show a higher propensity to cross barriers than subadults. To estimate the potential population size of the lynx in the Jura Mountains and to assess possible corridors between this population and adjacent areas, I adapted a habitat probability model for lynx distribution in the Jura Mountains with new environmental data and extrapolated it over the entire mountain range. The model predicts a breeding population ranging from 74-101 individuals and from 51-79 individuals when continuous habitat patches < 50 km2 are disregarded. The Jura Mountains could once be part of a metapopulation, as potential corridors exist to the adjoining areas (Alps, Vosges Mountains, and Black Forest). Monitoring of the population size, spatial expansion, and the genetic surveillance in the Jura Mountains must be continued, as the status of the population is still critical. ENFA was used to predict the potential distribution of lynx in the Alps. The resulting model divided the Alps into 37 suitable habitat patches ranging from 50 to 18,711 km2, covering a total area of about 93,600 km2. When using the range of lynx densities found in field studies in Switzerland, the Alps could host a population of 961 to 1,827 residents. The results of the cost-distance analysis revealed that all patches were within the reach of dispersing lynx, as the connection costs were in the range of dispersal cost of radio-tagged subadult lynx moving through unfavorable habitat. Thus, the whole Alps could once be considered as a metapopulation. But experience suggests that only few disperser will cross unsuitable areas and barriers. This low migration rate may seldom allow the spontaneous foundation of new populations in unsettled areas. As an alternative to natural dispersal, artificial transfer of individuals across the barriers should be considered. Wildlife biologists can play a crucial role in developing adaptive management experiments to help managers learning by trial. The case of the lynx in Switzerland is a good example of a fruitful cooperation between wildlife biologists, managers, decision makers and politician in an adaptive management process. This cooperation resulted in a Lynx Management Plan which was implemented in 2000 and updated in 2004 to give the cantons directives on how to handle lynx-related problems. This plan was put into practice e.g. in regard to translocation of lynx into unsettled areas. Résumé: L'expansion d'une population en phase de recolonisation, qu'elle soit issue de réintroductions ou d'un retour naturel dépend 1) de facteurs biologiques tels que le système social et le mode de dispersion, 2) de la distribution et la disponibilité des ressources (proies), 3) de l'habitat et des éléments du paysage, 4) de l'acceptation de l'espèce par la population locale et des activités humaines. Afin de pouvoir développer des stratégies efficaces de conservation et de favoriser la recolonisation, chacun de ces facteurs doit être pris en compte. En plus, la distribution potentielle de l'espèce doit pouvoir être déterminée et enfin, toutes les possibilités pour atteindre les objectifs, examinées. La phase de haute densité que la population de lynx a connue dans les années nonante dans le nord-ouest des Alpes suisses a donné lieu à une controverse assez vive. La protection du lynx dans de nombreux pays européens, promue par différentes organisations, a entraîné des conséquences inattendues; ces dernières montrent que tout plan de gestion doit impérativement indiquer des pistes quant à la manière de gérer les conflits, tout en trouvant un équilibre entre l'extinction et la surabondance de l'espèce. Les biologistes de la conservation et de la faune sauvage doivent pour cela évaluer les différents risques encourus par les populations de lynx, afin de pouvoir rapidement prendre les meilleuresmdécisions de gestion. Un modèle d'habitat pour le lynx, basé sur des caractéristiques de l'habitat et des données radio télémétriques collectées dans la chaîne du Jura, a été élaboré afin de prédire la distribution potentielle dans cette région, qui n'est que partiellement occupée par l'espèce. Trois des 18 variables testées, décrivant pour chaque kilomètre carré l'utilisation du sol, la végétation ainsi que la topographie, ont été retenues pour déterminer la probabilité de présence du lynx. La carte qui en résulte a été comparée aux données télémétriques de lynx subadultes en phase de dispersion. Les jeunes qui n'ont pas pu établir leur domaine vital dans l'habitat favorable prédit par le modèle n'ont pas survécu leur première année d'indépendance alors que le seul individu qui est mort dans l'habitat favorable a été braconné. Les données radio-télémétriques sont souvent utilisées pour l'étalonnage de modèles d'habitat. C'est un des seuls moyens à disposition qui permette de récolter des données non biaisées et précises sur l'occupation de l'habitat par des mammifères terrestres aux moeurs discrètes. Mais ces méthodes de- mandent un important investissement en moyens financiers et en temps et peuvent, de ce fait, n'être appliquées qu'à des zones limitées. Les modèles d'habitat sont ainsi souvent extrapolés à de grandes surfaces malgré le risque d'imprécision, qui résulte des variations des caractéristiques et de la disponibilité de l'habitat d'une zone à l'autre. Le pouvoir de prédiction de l'Analyse Ecologique de la Niche (AEN) dans les zones où les données de présence n'ont pas été prises en compte dans le calibrage du modèle a été analysée dans le cas du lynx en Suisse. D'après les résultats obtenus, la meilleure mé- thode pour prédire la distribution du lynx dans une zone alpine dépourvue d'indices de présence est de combiner des données provenant de régions contrastées (Alpes, Jura). Par contre, seules les données sur la présence locale de l'espèce doivent être utilisées pour les zones présentant une faible variance écologique tel que le Jura. La dispersion influence la dynamique et la stabilité des populations, la distribution et l'abondance des espèces et détermine les caractéristiques spatiales et temporelles des communautés vivantes et des écosystèmes. Entre 1988 et 2001, le comportement spatio-temporel de lynx eurasiens subadultes de deux populations réintroduites en Suisse a été étudié, basé sur le suivi de 39 individus juvéniles dont 24 étaient munis d'un collier émetteur, afin de déterminer les facteurs qui influencent la dispersion. Les subadultes se sont séparés de leur mère à l'âge de 8 à 11 mois. Le sexe n'a pas eu d'influence sur le nombre d'individus ayant dispersés et la distance parcourue au cours de la dispersion. Comparé à l'ours et au loup, le lynx reste très modéré dans ses mouvements de dispersion. Tous les individus ayant dispersés se sont établis à proximité ou dans des zones déjà occupées par des lynx. Les distances parcourues lors de la dispersion ont été plus courtes pour la population en phase de haute densité que celles relevées par les autres études de dispersion du lynx eurasien. Les zones d'habitat peu favorables et les barrières qui interrompent la connectivité entre les populations sont les principales entraves aux déplacements, lors de la dispersion. Dans un premier temps, nous avons fait l'hypothèse que les phases de haute densité favorisaient l'expansion des populations. Mais cette hypothèse a été infirmée par les résultats issus du suivi des lynx réalisé dans le nord-ouest des Alpes, où la population connaissait une phase de haute densité depuis 1995. Ce constat est important pour la conservation d'une population de carnivores dans un habitat fragmenté. Ainsi, instaurer une forte population source à un seul endroit n'est pas forcément la stratégie la plus judicieuse. Il est préférable d'établir des noyaux de populations dans des régions voisines où l'habitat est favorable. Des échanges entre des populations avoisinantes pourront avoir lieu par la suite car les lynx adultes sont plus enclins à franchir les barrières qui entravent leurs déplacements que les individus subadultes. Afin d'estimer la taille de la population de lynx dans le Jura et de déterminer les corridors potentiels entre cette région et les zones avoisinantes, un modèle d'habitat a été utilisé, basé sur un nouveau jeu de variables environnementales et extrapolé à l'ensemble du Jura. Le modèle prédit une population reproductrice de 74 à 101 individus et de 51 à 79 individus lorsque les surfaces d'habitat d'un seul tenant de moins de 50 km2 sont soustraites. Comme des corridors potentiels existent effectivement entre le Jura et les régions avoisinantes (Alpes, Vosges, et Forêt Noire), le Jura pourrait faire partie à l'avenir d'une métapopulation, lorsque les zones avoisinantes seront colonisées par l'espèce. La surveillance de la taille de la population, de son expansion spatiale et de sa structure génétique doit être maintenue car le statut de cette population est encore critique. L'AEN a également été utilisée pour prédire l'habitat favorable du lynx dans les Alpes. Le modèle qui en résulte divise les Alpes en 37 sous-unités d'habitat favorable dont la surface varie de 50 à 18'711 km2, pour une superficie totale de 93'600 km2. En utilisant le spectre des densités observées dans les études radio-télémétriques effectuées en Suisse, les Alpes pourraient accueillir une population de lynx résidents variant de 961 à 1'827 individus. Les résultats des analyses de connectivité montrent que les sous-unités d'habitat favorable se situent à des distances telles que le coût de la dispersion pour l'espèce est admissible. L'ensemble des Alpes pourrait donc un jour former une métapopulation. Mais l'expérience montre que très peu d'individus traverseront des habitats peu favorables et des barrières au cours de leur dispersion. Ce faible taux de migration rendra difficile toute nouvelle implantation de populations dans des zones inoccupées. Une solution alternative existe cependant : transférer artificiellement des individus d'une zone à l'autre. Les biologistes spécialistes de la faune sauvage peuvent jouer un rôle important et complémentaire pour les gestionnaires de la faune, en les aidant à mener des expériences de gestion par essai. Le cas du lynx en Suisse est un bel exemple d'une collaboration fructueuse entre biologistes de la faune sauvage, gestionnaires, organes décisionnaires et politiciens. Cette coopération a permis l'élaboration du Concept Lynx Suisse qui est entré en vigueur en 2000 et remis à jour en 2004. Ce plan donne des directives aux cantons pour appréhender la problématique du lynx. Il y a déjà eu des applications concrètes sur le terrain, notamment par des translocations d'individus dans des zones encore inoccupées.

Carbon-isotope stratigraphy and ammonite faunal turnover for the Middle Jurassic in the Southern Iberian palaeomargin

Variations in the stable carbon-isotope ratio of marine and continental sediments can reflect changes in sink and flux modifications of the palaeocarbon cycle. Here we report carbon-isotope compositions of Middle Jurassic marine carbonates from the Betic Cordillera (southern Spain), which represents an ideal region to link the stable carbon-isotope curves directly to ammonite zones and subzones, and thereby for the first time achieve an accurate chronostratigraphic calibration. The five sections studied represent basin and high swell deposits of the Southern Iberian palaeomargin. We find a similar delta C-13 of carbonates between different oceanic areas, suggesting a homogeneous carbon-isotope oceanic reservoir through the Middle Jurassic. The Aalenian-Bajocian transition is a critical period in ammonite evolution; hence the Early Jurassic fauna are replaced by new ammonite families which become dominant throughout the Middle and Late Jurassic. For this reason, we compared the delta C-13 values of carbonates with ammonite diversity and extinction rates at different taxonomical levels in order to explore the possible relationship between the carbon cycle and ammonite evolution. The carbon-isotope values of carbonates are not exactly linearly correlated with the extinction rate and ammonite diversity, but the main faunal turnovers follow minimum delta C-13 values, where extinct taxa are replaced by new ones. Likewise, radiation episodes are associated with increasing delta C-13 values and with transgressive sea-level rise. All these data support the idea that perturbations in the global carbon cycle reflect rapid palaeoenvironmental changes. We made detailed analyses of these faunal turnovers, using them as a proxy to identify major palaeoenvironmental crises in their ecosystems forced by modification in the carbon cycle. (c) 2006 Elsevier B.V All rights reserved.

A New Species of Sand Racer, Psammodromus (Squamata: Lacertidae), from the Western Iberian Peninsula

A new species of lacertid lizard of the genus Psammodromus is described from the Iberian Peninsula. Genetic and recently published phenotypic data support the differentiation of Psammodromus hispanicus into three, and not as previously suggested two, distinct lineages. Age estimates, lineage allopatry, the lack of mitochondrial and nuclear haplotype sharing between lineages, ecological niche divergence, and the current biogeographic distribution, indicated that the three lineages correspond to three independent species. Here, we describe a new species, Psammodromus occidentalis sp. n., which is genetically different from the other sand racers and differentiated by the number of femoral pores, number of throat scales, snout shape, head ratio, green nuptial coloration, and number of supralabial scales below the subocular scale. We also propose to upgrade the two previously recognized subspecies, Psammodromus hispanicus hispanicus Fitzinger, 1826 from central Spain and Psammodromus hispanicus edwardsianus (Dugès, 1829) from eastern Spain, to the species level: Psammodromus hispanicus stat. nov. and Psammodromus edwardsianus stat. nov. Given that the holotype of Psammodromus hispanicus was lost, we designate a neotype. We also analysed museum specimens of P. blanci, P. microdactylus and P. algirus to describe differentiation of the Psammodromus hispanicus lineages/species from their closest relatives.

Population structure of three Psammodromus species in the Iberian Peninsula.

The knowledge of a species' population structure is essential for the development of adequate conservation actions as well as for the understanding of its evolution. The population structure is unknown in all species of the Genus Psammodromus, including the Western Sand Racer (Psammodromus occidentalis; a recently described species), the Edward's Sand Racer (P. edwardsianus) and the Spanish Sand Racer (P. hispanicus). In this article, the genetic variability and population structure of Psammodromus edwardsianus, P. hispanicus, and P. occidentalis were studied in the Iberian Peninsula covering their natural geographic distribution. Mitochondrial DNA showed genetically different units in all species with higher genetic variability in their southern populations (latitudinal variation). Genetic differentiation was different among species and contrasted to those of species with similar characteristics. Our results therefore highlight the importance of species-specific studies analysing population structure.

Threats from climate change to terrestrial vertebrate hotspots in europe.

We identified hotspots of terrestrial vertebrate species diversity in Europe and adjacent islands. Moreover, we assessed the extent to which by the end of the 21(st) century such hotspots will be exposed to average monthly temperature and precipitation patterns which can be regarded as extreme if compared to the climate experienced during 1950-2000. In particular, we considered the entire European sub-continent plus Turkey and a total of 1149 species of terrestrial vertebrates. For each species, we developed species-specific expert-based distribution models (validated against field data) which we used to calculate species richness maps for mammals, breeding birds, amphibians, and reptiles. Considering four global circulation model outputs and three emission scenarios, we generated an index of risk of exposure to extreme climates, and we used a bivariate local Moran's I to identify the areas with a significant association between hotspots of diversity and high risk of exposure to extreme climates. Our results outline that the Mediterranean basin represents both an important hotspot for biodiversity and especially for threatened species for all taxa. In particular, the Iberian and Italian peninsulas host particularly high species richness as measured over all groups, while the eastern Mediterranean basin is particularly rich in amphibians and reptiles; the islands (both Macaronesian and Mediterranean) host the highest richness of threatened species for all taxa occurs. Our results suggest that the main hotspots of biodiversity for terrestrial vertebrates may be extensively influenced by the climate change projected to occur over the coming decades, especially in the Mediterranean bioregion, posing serious concerns for biodiversity conservation.

Subduction and obduction processes in the Swiss Alps

The significance of the Brianconnais domain in the Alpine orogen is reviewed in the light of data concerning its collision with the active Adriatic margin and the passive Helvetic margin. The Brianconnais which formerly belonged to the Iberian plate, was located on the northern margin of the Alpine Tethys (Liguro-Piemont ocean) since its opening in the early-Middle Jurassic. Together with the Iberian plate the Brianconnais terrane was separated from the European plate in the Late Jurassic-Early Cretaceous, following the northern Atlantic, Bay of Biscay, Valais ocean opening. This was accompanied by the onset of subduction along the northern margin of Adria and the closure of the Alpine Tethys. Stratigraphic and metamorphic data regarding this subduction and the geohistory of the Brianconnais allows the scenario of subduction-obduction processes during the Late Cretaceous-early Tertiary in the eastern and western Alps to be specified. HP-LT metamorphism record a long-lasting history of oceanic subduction-accretion, followed in the Middle Eocene by the incorporation of the Brianconnais as an exotic terrane into the accretionary prism. Middle to Late Eocene cooling ages of the Brianconnais basement and the presence of pelagic, anorogenic sedimentation lasting until the Middle Eocene on the Brianconnais preclude any sort of collision before that time between this domain and the active Adria margin or the Helvetic margin. This is confirmed by plate reconstructions constrained by magnetic anomalies in the Atlantic domain. Only a small percentage of the former Brianconnais domain was obducted, most of the crust and lithospheric roots were subducted. This applies also to domains formerly belonging to the southern Alpine Tethys margin (Austroalpine-inner Carpathian domain). It is proposed that there was a single Palaeogene subduction zone responsible for the Alpine orogen formation (from northern Spain to the East Carpathians), with the exception of a short-lived Late Cretaceous partial closure of the Valais ocean. Subduction in the western Tethyan domain originated during the closure of the Meliata ocean during the Jurassic incorporating the Austroalpine-Carpathian domain as terranes during the Cretaceous. The subduction zone propagated into the northern margin of Adria and then to the northern margin of the Iberian plate, where it gave birth to the Pyrenean-Provencal orogenic belt. This implies the absence of a separated Cretaceous subduction zone within the Austro-Carpathian Penninic ocean. Collision of Iberia with Europe forced the subduction to jump to the SE margin of Iberia in the Eocene, creating the Apenninic orogenic wedge and inverting the vergence of subduction from south- to north-directed. (C) 1998 Elsevier Science B.V. All rights reserved.

THE USE OF SIMULATIONS IN EVOLUTIONARY POPULATION GENETICS: APPLICATIONS ON HUMANS, OWLS AND VIRTUAL ORGANISMS

Summary (in English) Computer simulations provide a practical way to address scientific questions that would be otherwise intractable. In evolutionary biology, and in population genetics in particular, the investigation of evolutionary processes frequently involves the implementation of complex models, making simulations a particularly valuable tool in the area. In this thesis work, I explored three questions involving the geographical range expansion of populations, taking advantage of spatially explicit simulations coupled with approximate Bayesian computation. First, the neutral evolutionary history of the human spread around the world was investigated, leading to a surprisingly simple model: A straightforward diffusion process of migrations from east Africa throughout a world map with homogeneous landmasses replicated to very large extent the complex patterns observed in real human populations, suggesting a more continuous (as opposed to structured) view of the distribution of modern human genetic diversity, which may play a better role as a base model for further studies. Second, the postglacial evolution of the European barn owl, with the formation of a remarkable coat-color cline, was inspected with two rounds of simulations: (i) determine the demographic background history and (ii) test the probability of a phenotypic cline, like the one observed in the natural populations, to appear without natural selection. We verified that the modern barn owl population originated from a single Iberian refugium and that they formed their color cline, not due to neutral evolution, but with the necessary participation of selection. The third and last part of this thesis refers to a simulation-only study inspired by the barn owl case above. In this chapter, we showed that selection is, indeed, effective during range expansions and that it leaves a distinguished signature, which can then be used to detect and measure natural selection in range-expanding populations. Résumé (en français) Les simulations fournissent un moyen pratique pour répondre à des questions scientifiques qui seraient inabordable autrement. En génétique des populations, l'étude des processus évolutifs implique souvent la mise en oeuvre de modèles complexes, et les simulations sont un outil particulièrement précieux dans ce domaine. Dans cette thèse, j'ai exploré trois questions en utilisant des simulations spatialement explicites dans un cadre de calculs Bayésiens approximés (approximate Bayesian computation : ABC). Tout d'abord, l'histoire de la colonisation humaine mondiale et de l'évolution de parties neutres du génome a été étudiée grâce à un modèle étonnement simple. Un processus de diffusion des migrants de l'Afrique orientale à travers un monde avec des masses terrestres homogènes a reproduit, dans une très large mesure, les signatures génétiques complexes observées dans les populations humaines réelles. Un tel modèle continu (opposé à un modèle structuré en populations) pourrait être très utile comme modèle de base dans l'étude de génétique humaine à l'avenir. Deuxièmement, l'évolution postglaciaire d'un gradient de couleur chez l'Effraie des clocher (Tyto alba) Européenne, a été examiné avec deux séries de simulations pour : (i) déterminer l'histoire démographique de base et (ii) tester la probabilité qu'un gradient phénotypique, tel qu'observé dans les populations naturelles puisse apparaître sans sélection naturelle. Nous avons montré que la population actuelle des chouettes est sortie d'un unique refuge ibérique et que le gradient de couleur ne peux pas s'être formé de manière neutre (sans l'action de la sélection naturelle). La troisième partie de cette thèse se réfère à une étude par simulations inspirée par l'étude de l'Effraie. Dans ce dernier chapitre, nous avons montré que la sélection est, en effet, aussi efficace dans les cas d'expansion d'aire de distribution et qu'elle laisse une signature unique, qui peut être utilisée pour la détecter et estimer sa force.

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.

Integrative analyses of speciation and divergence in Psammodromus hispanicus (Squamata: Lacertidae).

ABSTRACT: BACKGROUND: Genetic, phenotypic and ecological divergence within a lineage is the result of past and ongoing evolutionary processes, which lead ultimately to diversification and speciation. Integrative analyses allow linking diversification to geological, climatic, and ecological events, and thus disentangling the relative importance of different evolutionary drivers in generating and maintaining current species richness. RESULTS: Here, we use phylogenetic, phenotypic, geographic, and environmental data to investigate diversification in the Spanish sand racer (Psammodromus hispanicus). Phylogenetic, molecular clock dating, and phenotypic analyses show that P. hispanicus consists of three lineages. One lineage from Western Spain diverged 8.3 (2.9-14.7) Mya from the ancestor of Psammodromus hispanicus edwardsianus and P. hispanicus hispanicus Central lineage. The latter diverged 4.8 (1.5-8.7) Mya. Molecular clock dating, together with population genetic analyses, indicate that the three lineages experienced northward range expansions from southern Iberian refugia during Pleistocene glacial periods. Ecological niche modelling shows that suitable habitat of the Western lineage and P. h. edwardsianus overlap over vast areas, but that a barrier may hinder dispersal and genetic mixing of populations of both lineages. P. h. hispanicus Central lineage inhabits an ecological niche that overlaps marginally with the other two lineages. CONCLUSIONS: Our results provide evidence for divergence in allopatry and niche conservatism between the Western lineage and the ancestor of P. h. edwardsianus and P. h. hispanicus Central lineage, whereas they suggest that niche divergence is involved in the origin of the latter two lineages. Both processes were temporally separated and may be responsible for the here documented genetic and phenotypic diversity of P. hispanicus. The temporal pattern is in line with those proposed for other animal lineages. It suggests that geographic isolation and vicariance played an important role in the early diversification of the group, and that lineage diversification was further amplified through ecological divergence.

Cryptic diversity among Western Palearctic tree frogs: postglacial range expansion, range limits, and secondary contacts of three European tree frog lineages (Hyla arborea group).

We characterize divergence times, intraspecific diversity and distributions for recently recognized lineages within the Hyla arborea species group, based on mitochondrial and nuclear sequences from 160 localities spanning its whole distribution. Lineages of H. arborea, H. orientalis, H. molleri have at least Pliocene age, supporting species level divergence. The genetically uniform Iberian H. molleri, although largely isolated by the Pyrenees, is parapatric to H. arborea, with evidence for successful hybridization in a small Aquitanian corridor (southwestern France), where the distribution also overlaps with H. meridionalis. The genetically uniform H. arborea, spread from Crete to Brittany, exhibits molecular signatures of a postglacial range expansion. It meets different mtDNA clades of H. orientalis in NE-Greece, along the Carpathians, and in Poland along the Vistula River (there including hybridization). The East-European H. orientalis is strongly structured genetically. Five geographic mitochondrial clades are recognized, with a molecular signature of postglacial range expansions for the clade that reached the most northern latitudes. Hybridization with H. savignyi is suggested in southwestern Turkey. Thus, cryptic diversity in these Pliocene Hyla lineages covers three extremes: a genetically poor, quasi-Iberian endemic (H. molleri), a more uniform species distributed from the Balkans to Western Europe (H. arborea), and a well-structured Asia Minor-Eastern European species (H. orientalis).

The Toarcian in the Subbetic basin (southern Spain): Bio-events (ammonite and calcareous nannofossils) and carbon-isotope stratigraphy

A detailed carbon-isotope stratigraphic study for the uppermost Pliensbachian lowermost Aalenian interval in the Median Subbetic palaeogeographic domain (External zones of the Betic Cordillera, southern Spain) has been carried out. During the Early Jurassic, the Median Subbetic, which represents a typical basin of the Hispanic Corridor connecting the Tethys and the Eastern Pacific, was located in the westernmost Tethys. The analyzed sections encompass the entire Toarcian stage as represented in the southern Iberian palaeomargin. Rocks are mainly rhythmic sequences of grey marls and marly limestones containing a rich ammonite fauna, nannofossils, and benthic foraminifers-all these provide an accurate biostratigraphic control. The lower and upper Toarcian boundaries are well represented in some of these sections and therefore represent optimal sites to link the carbon-isotope curves to ammonite zones, and to nannofossil events. delta C-13 values of bulk carbonates from the different localities of the Subbetic basin have similar variations from the uppermost Pliensbachian to the lowermost Aalenian, suggesting changes in the original DIC carbon isotope composition along the Hispanic corridor. The transition from Pliensbachian to Toarcian is marked by increasing delta C-13 values from similar to 12 to 2.0 parts per thousand, interrupted in the Serpentinum Zone by a negative shift concomitant with the Toarcian oceanic anoxic event (T-OAE), with the major ammonite extinction event of the Toarcian, and an important turnover of calcareous nannoplankton. The negative shift observed in the Serpentinum Zone confirms the global perturbation of the carbon cycling documented along the Tethys and the palaeo-Pacific in organic material and in marine carbonates. However, the amplitude of the negative excursion (similar to - 1.5 parts per thousand) is not compatible with an isotopic homogeneous seawater DIC and/or CO2 atmospheric reservoirs. The interval from the middle to the top of the Toarcian delta C-13 shows relatively constant values, minor ammonite turnovers, and is associated with increasing diversity of calcareous nannoplankton. (c) 2012 Elsevier B.V. All rights reserved.

Pliocene and Pleistocene diversification and multiple refugia in a Eurasian shrew (Crocidura suaveolens group).

We sequenced 998 base pairs (bp) of mitochondrial DNA cytochrome b and 799 bp of nuclear gene BRCA1 in the Lesser white-toothed shrew (Crocidura suaveolens group) over its geographic range from Portugal to Japan. The aims of the study were to identify the main clades within the group and respective refugia resulting from Pleistocene glaciations. Analyses revealed the Asian lesser white-toothed shrew (C. shantungensis) as the basal clade, followed by a major branch of C. suaveolens, subdivided sensu stricto into six clades, which split-up in the Upper Pliocene and Lower Pleistocene (1.9-0.9 Myr). The largest clade, occurring over a huge range from east Europe to Mongolia, shows evidence of population expansion after a bottleneck. West European clades originated from Iberian and Italo-Balkanic refugia. In the Near East, three clades evolved in an apparent hotspot of refugia (west Turkey, south-west and south-east of the Caucasus). Most clades include specimens of different morphotypes and the validity of many taxa in the C. suaveolens group has to be re-evaluated.

Molecular and cellular analysis of genodermatoses

Summary Skin is the essential interface between our body and its environment; not only does it prevent water loss and protect us from external insults it also plays an essential role in the central nervous system acting as a major sense organ primarily for touch and pain. The main cell type present in skin, keratinocyte, undergoes a differentiation process leading to the formation of this protecting barrier. This work is intended to contribute to the understanding of how keratinocyte differentiates and skin functions. To do this, we studied two genetic skin diseases: Erythrokeratodermia variabilis and Mal de Meleda. Our approach was to examine the expression and localization of proteins implicated in these two pathologies in normal and diseased tissues and to determine the influence of mutant proteins at the molecular and cellular levels. Connexins are major components of gap junctions, channels allowing direct communication between cells. Our laboratory has identified mutations in both connexin 30.3 (Cx30.3) and 31 (Cx31) to be causally involved in erythrokeratodermia variabilis (EKV), an autosomal dominant disorder of keratinization. In the first chapter, we show a new mutation of Cx31, L209P-Cx31, in 3 EKV patients, extending the field of EKV-causing mutations although the mechanism by which connexin mutations lead to the disease is unclear. In the second chapter, we studied the effect of F137L-Cx30.3 on expression, trafficking and localization of cotransfected Cx31 and Cx30.3 in connexin-deficient HeLa cells. The F137 amino acid, highly conserved in connexin family, is oriented towards the channel pore and F137L mutation in either Cx30.3 or Cx31 lead to EKV. As two genes can lead to EKV when mutated, our hypothesis was that Cx31 and Cx30.3 might cooperate at a molecular level. We were able to demonstrate a physical interaction between Cx31 and Cx30.3. The presence of F137L-Cx30.3 disturbed the trafficking of both connexins, less connexins were integrated into gap junctions and thus, the coupling between cell was diminished. Connexins formed in the presence of F137L-Cx30.3 are degraded at their exit from the endoplasmic reticulum. In conclusion, our results indicate that the genetic heterogeneity of EKV is due to mutations in two interacting proteins. F137L-Cx30.3 has a dominant negative effect and affects Cx31, disturbing cellular communication in epidermal cells. Mal de Meleda is an autosomal recessive inflammatory and a keratotic palmoplantar skin disorder due to mutations in SLURP1 (secreted LY6/PLAUR-related protein 1). SLURP1 belongs to the LY6/PLAUR family of proteins and has the particularity of being secreted instead of being GPI-anchored. The high degree of structural similarity between SLURP1 and the three fingers motif of snake neurotoxins and LYNX 1-C suggests that this protein could interact with the neuronal acetylcholine receptors. In the third chapter, we show that SLURP1 potentiates responses of the a7 nicotinic acetylcholine receptor (nAchR) to acetylcholine. These results identify SLURP1 as a secreted epidermal neuromodulator that is likely to be essential for palmoplantar skin. In the fourth chapter, we show that SLURP1 is expressed in the granular layer of the epidermis but is absent from skin biopsies of Mal de Meleda patients. SLURP1 is also present in secretions such as sweat, tears or saliva. An in vitro analysis on two mutant of SLURP-I demonstrates that W15R-SLURP1 is absent in cells while G86R-SLURP1 is expressed and secreted, suggesting that SLURP1 can lead to the disease by either an absent or an abnormal protein. Finally, in the fifth chapter, we analyse the expression and biological properties of other LY6/PLAUR members, clustered around SLURP] on chromosome 8. Their GPI-anchored or secreted status were analysed in vitro. SLURP1, LYNX1-A and -B are secreted while LYPDC2 and LYNX 1-C are GPI anchored. Three of these proteins are expressed in the epidermis and in cultured keratinocytes. These results suggest that these LY6/PLAUR members may have an important role in skin homeostasis. Résumé Résumé La peau est la barrière essentielle entre notre corps et l'environnement, nous protégeant des agressions extérieures, de la déshydratation et assurant aussi un rôle dans le système nerveux central en tant qu'organe du toucher et de la douleur. Le principal type de cellules présent dans la peau est le kératinocyte qui suit un processus de différenciation aboutissant à la formation de cette barrière protectrice. Ce travail est destiné à comprendre la différenciation des kératinocytes et le fonctionnement de la peau. Pour cela, nous avons étudié deux maladies génodermatoses : l'Erthrokeratodermia Variabilis (EKV) et le Mal de Meleda. Nous avons examiné l'expression et la localisation des protéines impliquées dans ces deux pathologies dans des tissus normaux et malades puis déterminé l'influence des protéines mutantes aux niveaux moléculaires et cellulaires. Les connexines (Cx) sont les composants majeurs des jonctions communicantes, canaux permettant la communication directe entre les cellules. Notre laboratoire a identifié des mutations dans les Cx30.3 et Cx31 comme responsables de l'EKV, génodermatose de transmission autosomique dominante. Dans le ler chapitre, nous décrivons une nouvelle mutation de Cx31, L209-Cx31, et contribuons à l'établissement du catalogue des mutations de Cx31 entraînant cette maladie. Cependant, le mécanisme par lequel les mutations de Cx31 et C3x0.3 provoquent l'EKV est inconnu. Dans le 2ème chapitre, nous étudions les effets de la mutation F137L-Cx30.3 sur l'expression, le trafic et la localisation des Cx31 et Cx30.3 transfectées dans des cellules HeLa, déficientes en connexines. Comme deux gènes peuvent causer une EKV quand ils sont mutés, notre hypothèse était que Cx31 et Cx30.3 pourraient coopérer au niveau moléculaire. Nous avons montré l'existence d'une interaction physique entre ces deux connexines. La présence de la mutation F137L-Cx30.3 perturbe le trafic des deux connexines, moins de connexines sont intégrées dans les jonctions communicantes et donc le couplage entre les cellules est diminué. Les connexons formés en présence de cette mutation sont dégradés à leur sortie du réticulum endoplasmique. En conclusion, nos résultats indiquent que l'hétérogénéité génétique de EKV est due à des mutations dans deux protéines qui interagissent. F137L-Cx30.3 a un effet dominant négatif et affecte Cx31, perturbant la communication entre les cellules épidermiques. Le Mal de Meleda est une maladie récessive de la peau palmoplantaire due à des mutations dans SLURP1. SLURP1 appartient à la famille des protéines contenant un domaine LY6/PLAUR et a la particularité d'être sécrétée. La grande homologie de structure existant entre SLURP1, les neurotoxines de serpent et LYNX1-C suggère que la protéine pourrait interagir avec des récepteurs à acétylcholine (Ach). Dans le 3ème chapitre, nous montrons que SLURP1 module la réponse à l'Ach du récepteur nicotinique α7. Ces résultats identifient SLURP1 comme un neuromodulateur épidermique sécrété, probablement essentiel pour la peau palmoplantaire. Dans le 4ème chapitre, nous montrons que SLURP1 est exprimé dans la couche granuleuse de l'épiderme et qu'il est absent des biopsies des patients. SLURP1 a aussi été détecté dans des sécrétions telles que la sueur, les lamies et la salive. Une analyse in vitro de deux mutants de SLURP1 a montré que W15R-SLURP1 est absent des cellules tandis que G86R-SLURP1 est exprimé et sécrété, suggérant qu'une absence ou une anomalie de SLURP1 peuvent causer la maladie. Finalement, dans le 5ème chapitre, nous analysons l'expression et les propriétés biologiques d'autres membres de la famille LY6/PLAUR localisés autour de SLURP1 sur le chromosome 8. Leur statut de protéines sécrétées ou liées à la membrane par une ancre GPI est analysé in vitro. SLURP1, LYNXI-A et -B sont sécrétées alors que LYPDC2 et LYNX1-C sont liés à la membrane. Trois de ces protéines sont exprimées dans l'épiderme et dans des kératinocytes cultivés. Ces résultats suggèrent que la famille LY6/PLAUR pourrait avoir un rôle important dans l'homéostasie de la peau.