Phylogeography, conservation genetics and parasitology of chimpanzees (Pan troglodytes versus) in Guinea-Bissau, West Africa

Tese de doutoramento em Antropologia, especialidade em Antropologia Biológica e Etnoecologia

Conservation genetics and population monitoring of the Alpine Salamander (Salamandra Atra) :

SUMMARYIn the context of the biodiversity crisis, amphibians are experiencing the most severe worldwide decline of all vertebrates and are in urgent need of better management. Efficient conservation strategies rely on sound knowledge of the species biology and of the genetic and demographic processes that might impair their welfare. Nonetheless, these processes are poorly understood in amphibians. Delineating population boundaries remains consequently problematic for these species, while it is of critical importance to define adequate management units for conservation. In this study, our attention focused on the alpine salamander (Salamandra atra), a species that deserves much interest in terms of both conservation biology and evolution. This endemic alpine species shows peculiar life-history traits (viviparity, reduced activity period, slow maturation) and has a slow population turnover, which might be problematic for its persistence in a changing environment. Due to its elusive behaviour (individuals spend most of their time underground and are unavailable for sampling), dynamic processes of gene and individuals were poorly understood for that species. Consequently, its conservation status could hardly be reliably assessed. Similarly the fire salamander (Salamandra salamandra) also poses special challenges for conservation, as no clear demarcation of geographical populations exists and dispersal patterns are poorly known. Through a phylogeographic analysis, we first studied the evolutionary history of the alpine salamander to better document the distribution of the genetic diversity along its geographical range. This study highlighted the presence of multiple divergent lineages in Italy together with a clear genetic divergence between populations from Northern and Dinaric Alps. These signs of cryptic genetic differentiation, which are not accounted for by the current taxonomy of the species, should not be neglected for further definition of conservation units. In addition, our data supported glacial survival of the species in northern peripheral glacial réfugia and nunataks, a pattern rarely documented for long-lived species. Then, we evaluated the level of gene flow between populations at the local scale and tested for asymmetries in male versus female dispersal using both field-based (mark-recapture) and genetic approaches. This study revealed high level of gene flow between populations, which stems mainly from male dispersal. This corroborated the idea that salamanders are much better dispersers than hitherto thought and provided a well- supported example of male-biased dispersal in amphibians. In a third step, based on a mark- recapture survey, we addressed the problem of sampling unavailability in alpine salamanders and evaluated its impact on two monitoring methods. We showed that about three quarters of individuals were unavailable for sampling during sampling sessions, a proportion that can vary with climatic conditions. If not taken into account, these complexities would result in false assumptions on population trends and misdirect conservation efforts. Finally, regarding the daunting task of delineating management units, our attention was drawn on the fire salamander. We conducted a local population genetic study that revealed high levels of gene flow among sampling sites. Management units for this species should consequently be large. Interestingly, despite the presence of several landscape features often reported to act as barriers, genetic breaks occurred at unexpected places. This suggests that landscape features may rather have idiosyncratic effects on population structure. In conclusion, this work brought new insights on both genetic and demographic processes occurring in salamanders. The results suggest that some biological paradigms should be taken with caution when particular species are in focus. Species- specific studies remain thus fundamental for a better understanding of species evolution and conservation, particularly in the context of current global changes.RESUMEDans le contexte de la crise de la biodiversité actuelle, les amphibiens subissent le déclin le plus important de tous les vertébrés et ont urgemment besoin d'une meilleure protection. L'établissement de stratégies de conservation efficaces repose sur des connaissances solides de la biologie des espèces et des processus génétiques et démographiques pouvant menacer leur survie. Ces processus sont néanmoins encore peu étudiés chez les amphibiens.Dans cette étude, notre attention s'est portée sur la salamandre noire (Salamandra atra), une espèce endémique des Alpes dont les traits d'histoire de vie atypiques (viviparité, phase d'activité réduite, lent turnover des populations) pourraient la rendre très vulnérable face aux changements environnementaux. Par ailleurs, en raison de son comportement cryptique (les individus passent la plupart de leur temps sous terre) la dynamique des gènes et des individus est mal comprise chez cette espèce. Il est donc difficile d'évaluer son statut de conservation de manière fiable. La salamandre tachetée {Salamandra salamandra), pour qui il n'existe aucune démarcation géographique apparente des populations, pose également des problèmes en termes de gestion. Dans un premier temps, nous avons étudié l'histoire évolutive de la salamandre noire afin de mieux décrire la distribution de sa diversité génétique au sein de son aire géographique. Cela a permis de mettre en évidence la présence de multiples lignées en Italie, ainsi qu'une nette divergence entre les populations du nord des Alpes et des Alpes dinariques. Ces résultats seront à prendre en compte lorsqu'il s'agira de définir des unités de conservation pour cette espèce. D'autre part, nos données soutiennent l'hypothèse d'une survie glaciaire dans des refuges nordiques périglaciaires ou dans des nunataks, fait rarement documenté pour une espèce longévive. Nous avons ensuite évalué la différentiation génétique des populations à l'échelle locale, ce qui a révélé d'important flux de gènes, ainsi qu'une asymétrie de dispersion en faveur des mâles. Ces résultats corroborent l'idée que les amphibiens dispersent mieux que ce que l'on pensait, et fournissent un exemple robuste de dispersion biaisée en faveur des mâles chez les amphibiens. Nous avons ensuite abordé le problème de Γ inaccessibilité des individus à la capture. Nous avons montré qu'environ trois quarts des individus sont inaccessibles lors des échantillonnages, une proportion qui peut varier en fonction des conditions climatiques. Ignoré, ce processus pourrait entraîner une mauvaise interprétation des fluctuations de populations ainsi qu'une mauvaise allocation des efforts de conservation. Concernant la définition d'unités de gestion pour la salamandre tachetée, nous avons pu mettre en évidence un flux de gènes important entre les sites échantillonnés. Les unités de gestion pour cette espèce devraient donc être étendues. Etonnamment, malgré la présence de nombreuses barrières potentielles au flux de gènes, les démarcations génétiques sont apparues à des endroits inattendus. En conclusion, ce travail a apporté une meilleure compréhension des processus génétiques et démographiques en action chez les salamandres. Les résultats suggèrent que certains paradigmes biologiques devraient être considérés avec précaution quand il s'agit de les appliquer à des espèces particulières. Les études spécifiques demeurent donc fondamentales pour une meilleure compréhension de l'évolution des espèces et leur conservation, tout particulièrement dans le contexte des changements globaux actuels.

Population genetics of Capercaillie (Tetrao urogallus) in the Jura and the Pyrenees: a non-invasive approach to avian conservation genetics

RÉSUMÉ Le Grand tétras est un galliforme de montagne apparenté au faisan et au tétras lyre. Il est distribué de manière continue à travers la toundra et les montagnes de moyenne altitude en Europe de l'ouest. Toutefois, les populations d'Europe de l'ouest ont subi un déclin constant au cours des derniers siècles. Les causes de ce déclin sont probablement liées à l'activité humaine, telle .que l'élevage ou le tourisme, qui ont engendré une modification et une fragmentation de l'habitat de l'espèce. Malheureusement, les populations soumises à de forts déclins démographiques peuvent subir des effets génétiques (augmentation de la consanguinité et perte de diversité génétique) pouvant diminuer leur potentiel de reproduction et conduire irrémédiablement à l'extinction. Cette thèse présente les analyses conduites dans le but d'estimer l'impact du déclin démographique des populations de Grand tétras sur l'étendue et la distribution de leur variabilité génétique dans le Jura et dans les Pyrénées. Du fait de la législation locale protégeant les tétraonidés en général, mais également en raison de la biologie très cryptique du Grand tétras, l'ensemble des analyses de cette étude a été réalisé à partir de matériel génétique extrait des fientes (ou échantillonnage génétique non invasif). Dans la première partie de l'étude, je détaille les protocoles d'extraction. d'ADN et d'amplification par PCR modifiés à partir des protocoles classiques utilisant des échantillons conventionnels, riches en ADN. L'utilisation d'ADN fécal impose des contraintes dues à la mauvaise qualité et à la faible quantité du matériel génétique à disposition dans les fientes. Ces contraintes ont pu être partiellement contournées en réalisant des répétitions multiples du génotypage afin d'obtenir un degré de fiabilité suffisante. J'ai également analysé les causes de la dégradation de l'ADN dans les excréments. Parmi les causes les plus communes, telles que l'activité bactérienne, l'hydrolyse spontanée et la dégradation enzymatique par les DNases libres, c'est ce dernier facteur qui apparaît comme étant la cause majeure et la plus rapide responsable de la dégradation de la qualité des échantillons. La rapidité de l'action enzymatique suggère que les plans d'échantillonnages de excréments sur le terrain pourraient être optimisés en les réalisant dans des conditions climatiques froides et sèches, favorisant ainsi l'inhibition des DNases. La seconde partie de la thèse est une étude par simulation visant à déterminer la capacité du logiciel Structure à identifier les structures génétiques complexes et hiérarchiques fréquemment rencontrées dans les populations naturelles, et ce en utilisant différents types de marqueurs génétiques. Les troisième et quatrième parties de cette thèse décrivent le statut génétique des populations résiduelles du Jura et des Pyrénées à partir de l'analyse de 11 loci microsatellites. Nous n'avons pas pu mettre en évidence dans les deux populations des effets liés à la consanguinité ou à la réduction de la diversité génétique. De plus, la différenciation génétique entre les patches d'habitats favorables reste modérée et corrélée à la distance géographique, ce qui suggère que la dispersion d'individus entre les patches a été importante au moins pendant ces dernières générations. La comparaison des paramètres de la diversité génétique avec ceux d'autres populations de Grand tétras, ou d'autres espèces proches, indique que la population du Jura a retenu une proportion importante de sa diversité originelle. Ces résultats suggèrent que le déclin récent des populations a jusqu'ici eu un impact modéré sur les facteurs génétiques et que ces populations semblent avoir conservé le potentiel génétique nécessaire à leur survie à long terme. Finalement, en cinquième partie, l'analyse de l'apparentement entre les mâles qui participent à la parade sur les places de chant (leks) indique que ces derniers sont distribués en agrégats de manière non aléatoire, préférentiellement entre individus apparentés. De plus, la corrélation entre les distances génétique et géographique entre les leks est en accord avec les motifs d'isolement par la distance mis en évidence à d'autres niveaux hiérarchiques (entre patches d'habitat et populations), ainsi qu'avec les études menées sur d'autres espèces ayant choisi ce même système de reproduction. En conclusion, cette première étude basée uniquement sur de l'ADN nucléaire aviaire extrait à partir de fèces a fourni des informations nouvelles qui n'auraient pas pu être obtenues par une méthode d'observation sur le terrain ou d'échantillonnage génétique classique. Aucun oiseau n'a été dérangé ou capturé, et les résultats sont comparables à d'autres études concernant des espèces proches. Néanmoins, la taille de ces populations approche des niveaux au-dessous desquels la survie à long terme est fortement incertaine. La persistance de la diversité génétique pour les prochaines générations reste en conséquence liée à la survie des adultes et à une reprise du succès de la reproduction. ABSTRACT Capercaillie (Tetrao urogallus) is a large grouse that is continuously distributed across the tundra and the mid-high mountains of Western Europe. However, the populations in Western Europe have been showing a constant decline during the last decades. The causes for this decline are possibly related to human activities, such as cattle breeding and tourism that have both led to habitat modification and fragmentation. Unfortunately, populations that have undergone drastic demographic bottlenecks often go through genetic processes of inbreeding and loss of diversity that decrease their fitness and eventually lead to extinction. This thesis presents the investigations conducted to estimate the impact of the demographic decline of capercaillie populations on the extent and distribution of their genetic variability in the Jura and in the Pyrenees mountains. Because grouse are protected by wildlife legislation, and also because of the cryptic behaviour of capercaillie, all DNA material used in this study was extracted from faeces (non-invasive genetic sampling). In the first part of my thesis, I detail the protocols of DNA extraction and PCR amplification adapted from classical methods using conventional DNA-rich samples. The use of faecal DNA imposes specific constraints due to the low quantity and the highly degraded genetic material available. These constraints are partially overcome by performing multiple genotyping repetitions to obtain sufficient reliability. I also investigate the causes of DNA degradation in faeces. Among the main degraders, namely bacterial activity, spontaneous hydrolysis, and free-¬DNase activities, the latter was pointed out as the most important according to our experiments. These enzymes degrade DNA very rapidly, and, as a consequence, faeces sampling schemes must be planned preferably in cold and dry weather conditions, allowing for enzyme activity inhibition. The second part of the thesis is a simulation study aiming to assess the capacity of the software Structure to detect population structure in hierarchical models relevant to situations encountered in wild populations, using several genetic markers. The methods implemented in Structure appear efficient in detecting the highest hierarchical structure. The third and fourth parts of the thesis describe the population genetics status of the remaining Jura and Pyrenees populations using 11 microsatellite loci. In either of these populations, no inbreeding nor reduced genetic diversity was detected. Furthermore, the genetic differentiation between patches defined by habitat suitability remains moderate and correlated with geographical distance, suggesting that significant dispersion between patches was at work at least until the last generations. The comparison of diversity indicators with other species or other populations of capercaillie indicate that population in the Jura has retained a large part of its original genetic diversity. These results suggest that the recent decline has had so forth a moderate impact on• genetic factors and that these populations might have retained the potential for long term survival, if the decline is stopped. Finally, in the fifth part, the analysis of relatedness between males participating in the reproduction parade, or lek, indicate that capercaillie males, like has been shown for some other grouse species, gather on leks• among individuals that are more related than the average of the population. This pattern appears to be due to both population structure and kin-association. As a conclusion, this first study relying exclusively on nuclear DNA extracted from faeces has provided novel information that was not available through field observation or classical genetic sampling. No bird has been captured or disturbed, and the results are consistent with other studies of closely related species. However, the size of these populations is approaching thresholds below which long-term survival is unlikely. The persistence of genetic diversity for the forthcoming generations remains therefore bond to adult survival and to the increase of reproduction success.

Conservation Genetics of Exploited Finnish Salmonid Fishes

Genetic diversity is one of the levels of biodiversity that the World Conservation Union (IUCN) has recognized as being important to preserve. This is because genetic diversity is fundamental to the future evolution and to the adaptive flexibility of a species to respond to the inherently dynamic nature of the natural world. Therefore, the key to maintaining biodiversity and healthy ecosystems is to identify, monitor and maintain locally-adapted populations, along with their unique gene pools, upon which future adaptation depends. Thus, conservation genetics deals with the genetic factors that affect extinction risk and the genetic management regimes required to minimize the risk. The conservation of exploited species, such as salmonid fishes, is particularly challenging due to the conflicts between different interest groups. In this thesis, I conduct a series of conservation genetic studies on primarily Finnish populations of two salmonid fish species (European grayling, Thymallus thymallus, and lake-run brown trout, Salmo trutta) which are popular recreational game fishes in Finland. The general aim of these studies was to apply and develop population genetic approaches to assist conservation and sustainable harvest of these populations. The approaches applied included: i) the characterization of population genetic structure at national and local scales; ii) the identification of management units and the prioritization of populations for conservation based on evolutionary forces shaping indigenous gene pools; iii) the detection of population declines and the testing of the assumptions underlying these tests; and iv) the evaluation of the contribution of natural populations to a mixed stock fishery. Based on microsatellite analyses, clear genetic structuring of exploited Finnish grayling and brown trout populations was detected at both national and local scales. Finnish grayling were clustered into three genetically distinct groups, corresponding to northern, Baltic and south-eastern geographic areas of Finland. The genetic differentiation among and within population groups of grayling ranged from moderate to high levels. Such strong genetic structuring combined with low genetic diversity strongly indicates that genetic drift plays a major role in the evolution of grayling populations. Further analyses of European grayling covering the majority of the species’ distribution range indicated a strong global footprint of population decline. Using a coalescent approach the beginning of population reduction was dated back to 1 000-10 000 years ago (ca. 200-2 000 generations). Forward simulations demonstrated that the bottleneck footprints measured using the M ratio can persist within small populations much longer than previously anticipated in the face of low levels of gene flow. In contrast to the M ratio, two alternative methods for genetic bottleneck detection identified recent bottlenecks in six grayling populations that warrant future monitoring. Consistent with the predominant role of random genetic drift, the effective population size (N_e) estimates of all grayling populations were very low with the majority of N_e estimates below 50. Taken together, highly structured local populations, limited gene flow and the small N_e of grayling populations indicates that grayling populations are vulnerable to overexploitation and, hence, monitoring and careful management using the precautionary principles is required not only in Finland but throughout Europe. Population genetic analyses of lake-run brown trout populations in the Inari basin (northernmost Finland) revealed hierarchical population structure where individual populations were clustered into three population groups largely corresponding to different geographic regions of the basin. Similar to my earlier work with European grayling, the genetic differentiation among and within population groups of lake-run brown trout was relatively high. Such strong differentiation indicated that the power to determine the relative contribution of populations in mixed fisheries should be relatively high. Consistent with these expectations, high accuracy and precision in mixed stock analysis (MSA) simulations were observed. Application of MSA to indigenous fish caught in the Inari basin identified altogether twelve populations that contributed significantly to mixed stock fisheries with the Ivalojoki river system being the major contributor (70%) to the total catch. When the contribution of wild trout populations to the fisheries was evaluated regionally, geographically nearby populations were the main contributors to the local catches. MSA also revealed a clear separation between the lower and upper reaches of Ivalojoki river system – in contrast to lower reaches of the Ivalojoki river that contributed considerably to the catch, populations from the upper reaches of the Ivalojoki river system (>140 km from the river mouth) did not contribute significantly to the fishery. This could be related to the available habitat size but also associated with a resident type life history and increased cost of migration. The studies in my thesis highlight the importance of dense sampling and wide population coverage at the scale being studied and also demonstrate the importance of critical evaluation of the underlying assumptions of the population genetic models and methods used. These results have important implications for conservation and sustainable fisheries management of Finnish populations of European grayling and brown trout in the Inari basin.

Mating system in Myracrodruon urundeuva (Anarcadiaceae): implications for conservation genetics

The aims of this study were to investigate the mating system of a fragmented population of the dioecious tropical tree Myracrodruon urundeuva Allemão, using five microsatellite loci and the mixed mating and correlated mating models. The study was conducted in the Estação Ecológica de Paulo de Farias (436 ha), where the population occupies about 142 ha. The mating system was estimated using 514 open-pollinated offspring, collected from 30 seed-trees. Estimates of the multilocus outcrossing rate confirm that the species is dioecious (t m = 1.0). Low levels of mating among relatives were detected in the population (1 - t s = 0.020). The estimate of paternity correlation (r p(m)) indicated that offsprings were composed of mixtures of half-sibs and full-sibs, with the latter occurring at a low frequency (average of 0.148). The estimated coancestry coefficient within families (Θ = 0.147) was larger and the effective population size (Ne(v)) was lower (Ne(v) = 2.98) than expected in progenies from panmictic populations (Θ = 0.125, Ne(v) = 4, respectively). In terms of conservation, the results indicate that to retain an effective population size of 150, is necessary to collect seeds from at least 50 seed-trees.

Mating system in Myracrodruon urundeuva (Anarcadiaceae): implications for conservation genetics

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Conservation genetics of the endangered Pampas deer (Ozotoceros bezoarticus)

The Pampas deer (Ozotoceros bezoarticus L. 1758) is the most endangered neotropical cervid, and in the past occupied a wide range of open habitats including grassland, pampas, savanna, and cerrado (Brazil) from 5 degrees to 41 degrees S. To better understand the effect of habitat fragmentation on gene flow and genetic variation, and to uncover genetic units for conservation, we examined DNA sequences from the mitochondrial control region of 54 individuals from six localities distributed throughout the present geographical range of the Pampas deer. Our results suggest that the control region of the Pampas deer is one of the most polymorphic of any mammal. This remarkably high variability probably reflects large historic population sizes of millions of individuals in contrast to numbers of fewer than 80 000 today. Gene flow between populations is generally close to one migrant per generation and, with the exception of two populations from Argentina, all populations are significantly differentiated. The degree of gene flow was correlated with geographical distance between populations, a result consistent with limited dispersal being the primary determinant of genetic differentiation between populations. The molecular genetic results provide a mandate for habitat restoration and reintroduction of Pampas deer so that levels of genetic variation can be preserved and historic patterns of abundance can be reconstructed. However, the source of individuals for reintroduction generally should be from populations geographically closest to those now in danger of extinction.

Novel tools for conservation genetics in marine fish: population structure and evolution of the European hake (Merluccius merluccius) inferred by SNP variation and applications to traceability

The research presented in my PhD thesis is part of a wider European project, FishPopTrace, focused on traceability of fish populations and products. My work was aimed at developing and analyzing novel genetic tools for a widely distributed marine fish species, the European hake (Merluccius merluccius), in order to investigate population genetic structure and explore potential applications to traceability scenarios. A total of 395 SNPs (Single Nucleotide Polymorphisms) were discovered from a massive collection of Expressed Sequence Tags, obtained by high-throughput sequencing, and validated on 19 geographic samples from Atlantic and Mediterranean. Genome-scan approaches were applied to identify polymorphisms on genes potentially under divergent selection (outlier SNPs), showing higher genetic differentiation among populations respect to the average observed across loci. Comparative analysis on population structure were carried out on putative neutral and outlier loci at wide (Atlantic and Mediterranean samples) and regional (samples within each basin) spatial scales, to disentangle the effects of demographic and adaptive evolutionary forces on European hake populations genetic structure. Results demonstrated the potential of outlier loci to unveil fine scale genetic structure, possibly identifying locally adapted populations, despite the weak signal showed from putative neutral SNPs. The application of outlier SNPs within the framework of fishery resources management was also explored. A minimum panel of SNP markers showing maximum discriminatory power was selected and applied to a traceability scenario aiming at identifying the basin (and hence the stock) of origin, Atlantic or Mediterranean, of individual fish. This case study illustrates how molecular analytical technologies have operational potential in real-world contexts, and more specifically, potential to support fisheries control and enforcement and fish and fish product traceability.

Conservation genetics of European wildcat (Felis silvestris silvestris): a wide and integrating analysis protocol for admixture inferences and population structure

Introgression of domestic cat genes into European wildcat (Felis silvestris silvestris) populations and reduction of wildcats’ range in Europe, leaded by habitat loss and fragmentation, are considered two of the main conservation problems for this endangered feline. This thesis addressed the questions related with the artificial hybridization and populations’ fragmentation, using a conservation genetics perspective. We combined the use of highly polymorphic loci, Bayesian statistical inferences and landscape analyses tools to investigate the origin of the geographic-genetic substructure of European wildcats (Felis silvestris silvestris) in Italy and Europe. The genetic variability of microsatellites evidenced that European wildcat populations currently distributed in Italy differentiated in, and expanded from two distinct glacial refuges during the Last Glacial Maximum. The genetic and geographic substructure detected between the eastern and western sides of the Apennine ridge, resulted by adaptation to specific ecological conditions of the Mediterranean habitats. European wildcat populations in Europe are strongly structured into 5 geographic-genetic macro clusters corresponding to: the Italian peninsular & Sicily; Balkans & north-eastern Italy; Germany eastern; central Europe; and Iberian Peninsula. Central European population might have differentiated in the extra-Mediterranean Würm ice age refuge areas (Northern Alps, Carpathians, and the Bulgarian mountain systems), while the divergence among and within the southern European populations might have resulted by the Pleistocene bio geographical framework of Europe, with three southern refugia localized in the Balkans, Italian Peninsula and Iberia Peninsula. We further combined the use of most informative autosomal SNPs with uniparental markers (mtDNA and Y-linked) for accurately detecting parental genotypes and levels of introgressive hybridization between European wild and domestic cats. A total of 11 hybrids were identified. The presence of domestic mitochondrial haplotypes shared with some wild individuals led us to hypnotize the possibility that ancient introgressive events might have occurred and that further investigation should be recommended.

Conservation genetics of North Atlantic loggerhead sea turtles: analysis of nuclear and mitochondrial DNA

[EN] Complex population structure has been described for the loggerhead sea turtle (Caretta caretta), revealing lower levels of population genetic structure in nuclear compared to mitochondrial DNA assays. This may result from mating during spatially overlapping breeding migrations, or male-biased dispersal as previously found for the green turtle (Chelonia mydas). To further investigate these multiple possibilities, we carried out a comparative analysis from twelve newly developed microsatellite loci and the mitochondrial DNA control region (~804 bp) in adult females of the Cape Verde Islands (n=158), and Georgia, USA (n=17).