Demogenetics in metapopulations

SUMMARY : Human-induced habitat fragmentation constitutes a major threat to biodiversity. Small and isolated populations suffer from increased stochasticity and from limited rescue effects. These two factors may be sufficient to cause local extinctions but fragmentation induces some genetic consequences that can also contribute significantly to extinction risks. Increased genetic drift reduces the effectiveness of selection against deleterious mutations, leading to their progressive accumulation. Drift also decreases both the standing genetic variation and the rate of fixation of beneficial mutations, limiting the evolutionary potential of isolated populations. Demography and genetics further interact and feed back on each other, progressively driving fragmented populations into "extinction vortices". The aim of the thesis was to better understand the processes occurring in fragmented populations. For this, I combined simulation studies and empirical data from three species that live in structured habitats. Chapter 1 and 2 investigate the demography of two shrew species in fragmented habitats. I showed that connectivity and habitat quality strongly affect the demography of the greater white-tooted shrew, although demographic stochasticity was extremely high. I also demonstrated that habitat fragmentation is one of the leading factors allowing the local coexistence of two competing shrew species. Chapter 3 and 4 focus on measuring connectivity in fragmented populations based on genetic data. In particular, I showed that genetic data can be used to detect the landscape elements impeding dispersal. In Chapter 5 that deals with the accumulation of deleterious mutations in fragmented populations, I demonstrated that mutation accumulation, as well a time to extinction, can be predicted from simple demographic and genetic measures. In the last two chapters, I monitored individual reproductive success in an isolated tree frogs population. These data allowed quantifying the effective population size, a measure closely linked to population evolutionary potential. To conclude, this thesis brings some new insights into the processes occurring in fragmented populations, and I hope it will contribute to the improvement of the management and conservation of fragmented populations.

How a haemosporidian parasite of bats gets around: the genetic structure of a parasite, vector and host compared.

Parasite population structure is often thought to be largely shaped by that of its host. In the case of a parasite with a complex life cycle, two host species, each with their own patterns of demography and migration, spread the parasite. However, the population structure of the parasite is predicted to resemble only that of the most vagile host species. In this study, we tested this prediction in the context of a vector-transmitted parasite. We sampled the haemosporidian parasite Polychromophilus melanipherus across its European range, together with its bat fly vector Nycteribia schmidlii and its host, the bent-winged bat Miniopterus schreibersii. Based on microsatellite analyses, the wingless vector, and not the bat host, was identified as the least structured population and should therefore be considered the most vagile host. Genetic distance matrices were compared for all three species based on a mitochondrial DNA fragment. Both host and vector populations followed an isolation-by-distance pattern across the Mediterranean, but not the parasite. Mantel tests found no correlation between the parasite and either the host or vector populations. We therefore found no support for our hypothesis; the parasite population structure matched neither vector nor host. Instead, we propose a model where the parasite's gene flow is represented by the added effects of host and vector dispersal patterns.

Linking landscape- and population ecology for large population management modelling: the case of Ibex (capra ibex) in Switzerland

La biologie de la conservation est communément associée à la protection de petites populations menacées d?extinction. Pourtant, il peut également être nécessaire de soumettre à gestion des populations surabondantes ou susceptibles d?une trop grande expansion, dans le but de prévenir les effets néfastes de la surpopulation. Du fait des différences tant quantitatives que qualitatives entre protection des petites populations et contrôle des grandes, il est nécessaire de disposer de modèles et de méthodes distinctes. L?objectif de ce travail a été de développer des modèles prédictifs de la dynamique des grandes populations, ainsi que des logiciels permettant de calculer les paramètres de ces modèles et de tester des scénarios de gestion. Le cas du Bouquetin des Alpes (Capra ibex ibex) - en forte expansion en Suisse depuis sa réintroduction au début du XXème siècle - servit d?exemple. Cette tâche fut accomplie en trois étapes : En premier lieu, un modèle de dynamique locale, spécifique au Bouquetin, fut développé : le modèle sous-jacent - structuré en classes d?âge et de sexe - est basé sur une matrice de Leslie à laquelle ont été ajoutées la densité-dépendance, la stochasticité environnementale et la chasse de régulation. Ce modèle fut implémenté dans un logiciel d?aide à la gestion - nommé SIM-Ibex - permettant la maintenance de données de recensements, l?estimation automatisée des paramètres, ainsi que l?ajustement et la simulation de stratégies de régulation. Mais la dynamique d?une population est influencée non seulement par des facteurs démographiques, mais aussi par la dispersion et la colonisation de nouveaux espaces. Il est donc nécessaire de pouvoir modéliser tant la qualité de l?habitat que les obstacles à la dispersion. Une collection de logiciels - nommée Biomapper - fut donc développée. Son module central est basé sur l?Analyse Factorielle de la Niche Ecologique (ENFA) dont le principe est de calculer des facteurs de marginalité et de spécialisation de la niche écologique à partir de prédicteurs environnementaux et de données d?observation de l?espèce. Tous les modules de Biomapper sont liés aux Systèmes d?Information Géographiques (SIG) ; ils couvrent toutes les opérations d?importation des données, préparation des prédicteurs, ENFA et calcul de la carte de qualité d?habitat, validation et traitement des résultats ; un module permet également de cartographier les barrières et les corridors de dispersion. Le domaine d?application de l?ENFA fut exploré par le biais d?une distribution d?espèce virtuelle. La comparaison à une méthode couramment utilisée pour construire des cartes de qualité d?habitat, le Modèle Linéaire Généralisé (GLM), montra qu?elle était particulièrement adaptée pour les espèces cryptiques ou en cours d?expansion. Les informations sur la démographie et le paysage furent finalement fusionnées en un modèle global. Une approche basée sur un automate cellulaire fut choisie, tant pour satisfaire aux contraintes du réalisme de la modélisation du paysage qu?à celles imposées par les grandes populations : la zone d?étude est modélisée par un pavage de cellules hexagonales, chacune caractérisée par des propriétés - une capacité de soutien et six taux d?imperméabilité quantifiant les échanges entre cellules adjacentes - et une variable, la densité de la population. Cette dernière varie en fonction de la reproduction et de la survie locale, ainsi que de la dispersion, sous l?influence de la densité-dépendance et de la stochasticité. Un logiciel - nommé HexaSpace - fut développé pour accomplir deux fonctions : 1° Calibrer l?automate sur la base de modèles de dynamique (par ex. calculés par SIM-Ibex) et d?une carte de qualité d?habitat (par ex. calculée par Biomapper). 2° Faire tourner des simulations. Il permet d?étudier l?expansion d?une espèce envahisseuse dans un paysage complexe composé de zones de qualité diverses et comportant des obstacles à la dispersion. Ce modèle fut appliqué à l?histoire de la réintroduction du Bouquetin dans les Alpes bernoises (Suisse). SIM-Ibex est actuellement utilisé par les gestionnaires de la faune et par les inspecteurs du gouvernement pour préparer et contrôler les plans de tir. Biomapper a été appliqué à plusieurs espèces (tant végétales qu?animales) à travers le Monde. De même, même si HexaSpace fut initialement conçu pour des espèces animales terrestres, il pourrait aisément être étndu à la propagation de plantes ou à la dispersion d?animaux volants. Ces logiciels étant conçus pour, à partir de données brutes, construire un modèle réaliste complexe, et du fait qu?ils sont dotés d?une interface d?utilisation intuitive, ils sont susceptibles de nombreuses applications en biologie de la conservation. En outre, ces approches peuvent également s?appliquer à des questions théoriques dans les domaines de l?écologie des populations et du paysage.<br/><br/>Conservation biology is commonly associated to small and endangered population protection. Nevertheless, large or potentially large populations may also need human management to prevent negative effects of overpopulation. As there are both qualitative and quantitative differences between small population protection and large population controlling, distinct methods and models are needed. The aim of this work was to develop theoretical models to predict large population dynamics, as well as computer tools to assess the parameters of these models and to test management scenarios. The alpine Ibex (Capra ibex ibex) - which experienced a spectacular increase since its reintroduction in Switzerland at the beginning of the 20th century - was used as paradigm species. This task was achieved in three steps: A local population dynamics model was first developed specifically for Ibex: the underlying age- and sex-structured model is based on a Leslie matrix approach with addition of density-dependence, environmental stochasticity and culling. This model was implemented into a management-support software - named SIM-Ibex - allowing census data maintenance, parameter automated assessment and culling strategies tuning and simulating. However population dynamics is driven not only by demographic factors, but also by dispersal and colonisation of new areas. Habitat suitability and obstacles modelling had therefore to be addressed. Thus, a software package - named Biomapper - was developed. Its central module is based on the Ecological Niche Factor Analysis (ENFA) whose principle is to compute niche marginality and specialisation factors from a set of environmental predictors and species presence data. All Biomapper modules are linked to Geographic Information Systems (GIS); they cover all operations of data importation, predictor preparation, ENFA and habitat suitability map computation, results validation and further processing; a module also allows mapping of dispersal barriers and corridors. ENFA application domain was then explored by means of a simulated species distribution. It was compared to a common habitat suitability assessing method, the Generalised Linear Model (GLM), and was proven better suited for spreading or cryptic species. Demography and landscape informations were finally merged into a global model. To cope with landscape realism and technical constraints of large population modelling, a cellular automaton approach was chosen: the study area is modelled by a lattice of hexagonal cells, each one characterised by a few fixed properties - a carrying capacity and six impermeability rates quantifying exchanges between adjacent cells - and one variable, population density. The later varies according to local reproduction/survival and dispersal dynamics, modified by density-dependence and stochasticity. A software - named HexaSpace - was developed, which achieves two functions: 1° Calibrating the automaton on the base of local population dynamics models (e.g., computed by SIM-Ibex) and a habitat suitability map (e.g. computed by Biomapper). 2° Running simulations. It allows studying the spreading of an invading species across a complex landscape made of variously suitable areas and dispersal barriers. This model was applied to the history of Ibex reintroduction in Bernese Alps (Switzerland). SIM-Ibex is now used by governmental wildlife managers to prepare and verify culling plans. Biomapper has been applied to several species (both plants and animals) all around the World. In the same way, whilst HexaSpace was originally designed for terrestrial animal species, it could be easily extended to model plant propagation or flying animals dispersal. As these softwares were designed to proceed from low-level data to build a complex realistic model and as they benefit from an intuitive user-interface, they may have many conservation applications. Moreover, theoretical questions in the fields of population and landscape ecology might also be addressed by these approaches.

Interviewer and Respondent Survey Quality Effects in a CATI Panel

Especially in panel surveys, respondent attrition, respondent learning, and interviewer experience effects play a crucial role with respect to data quality. We examine three interview survey quality indicators in the same survey in a cross sectional as well as in a longitudinal way. In the cross sectional analysis we compare data quality in the mature original sample with that in a refreshment sample, surveyed in the same wave. Because in the same wave an interviewer survey was conducted, collecting attitudes on their socio demography, survey attitudes and burden measures, we are able to consider interviewer fixed effects as well. The longitudinal analysis gives more insight in the respondent learning effects with respect to the quality indicators considered by considering the very same respondents across waves. The Swiss Household Panel, a CATI survey representative of the Swiss residential population, forms an ideal modelling database: the interviewer - respondent assignment is random, both within and across waves. This design avoids possible confusion with other effects stemming from a non-random assignment of interviewers, e.g. area effects or effects from assigning the best interviewers to the hard cases. In order to separate interviewer, respondent and wave effects, we build cross-classified multilevel models.

Sex ratio bias, male aggression, and population collapse in lizards.

The adult sex ratio (ASR) is a key parameter of the demography of human and other animal populations, yet the causes of variation in ASR, how individuals respond to this variation, and how their response feeds back into population dynamics remain poorly understood. A prevalent hypothesis is that ASR is regulated by intrasexual competition, which would cause more mortality or emigration in the sex of increasing frequency. Our experimental manipulation of populations of the common lizard (Lacerta vivipara) shows the opposite effect. Male mortality and emigration are not higher under male-biased ASR. Rather, an excess of adult males begets aggression toward adult females, whose survival and fecundity drop, along with their emigration rate. The ensuing prediction that adult male skew should be amplified and total population size should decline is supported by long-term data. Numerical projections show that this amplifying effect causes a major risk of population extinction. In general, such an "evolutionary trap" toward extinction threatens populations in which there is a substantial mating cost for females, and environmental changes or management practices skew the ASR toward males.

High connectivity in a long-lived High-Arctic seabird, the ivory gull Pagophila eburnea

Species may cope with rapid habitat changes by distribution shifts or adaptation to new conditions. A common feature of these responses is that they depend on how the process of dispersal connects populations, both demographically and genetically. We analyzed the genetic structure of a near-threatened high-Arctic seabird, the ivory gull (Pagophila eburnea) in order to infer the connectivity among gull colonies. We analyzed 343 individuals sampled from 16 localities across the circumpolar breeding range of ivory gulls, from northern Russia to the Canadian Arctic. To explore the roles of natal and breeding dispersal, we developed a population genetic model to relate dispersal behavior to the observed genetic structure of worldwide ivory gull populations. Our key finding is the striking genetic homogeneity of ivory gulls across their entire distribution range. The lack of population genetic structure found among colonies, in tandem with independent evidence of movement among colonies, suggests that ongoing effective dispersal is occurring across the Arctic Region. Our results contradict the dispersal patterns generally observed in seabirds where species movement capabilities are often not indicative of dispersal patterns. Model predictions show how natal and breeding dispersal may combine to shape the genetic homogeneity among ivory gull colonies separated by up to 2800 km. Although field data will be key to determine the role of dispersal for the demography of local colonies and refine the respective impacts of natal versus breeding dispersal, conservation planning needs to consider ivory gulls as a genetically homogeneous, Arctic-wide metapopulation effectively connected through dispersal.