Landscape genetics of the Alpine newt (Mesotriton alpestris) inferred from a strip-based approach

Habitat destruction and fragmentation are known to strongly affect dispersal by altering the quality of the environment between populations. As a consequence, lower landscape connectivity is expected to enhance extinction risks through a decrease in gene flow and the resulting negative effects of genetic drift, accumulation of deleterious mutations and inbreeding depression. Such phenomena are particularly harmful for amphibian species, characterized by disjunct breeding habitats. The dispersal behaviour of amphibians being poorly understood, it is crucial to develop new tools, allowing us to determine the influence of landscape connectivity on the persistence of populations. In this study, we developed a new landscape genetics approach that aims at identifying land-uses affecting genetic differentiation, without a priori assumptions about associated ecological costs. We surveyed genetic variation at seven microsatellite loci for 19 Alpine newt (Mesotriton alpestris) populations in western Switzerland. Using strips of varying widths that define a dispersal corridor between pairs of populations, we were able to identify land-uses that act as dispersal barriers (i.e. urban areas) and corridors (i.e. forests). Our results suggest that habitat destruction and landscape fragmentation might in the near future affect common species such as M. alpestris. In addition, by identifying relevant landscape variables influencing population structure without unrealistic assumptions about dispersal, our method offers a simple and flexible tool of investigation as an alternative to least-cost models and other approaches.

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.

Genomics of the divergence continuum in an African plant biodiversity hotspot, I: drivers of population divergence in Restio capensis (Restionaceae).

Understanding the drivers of population divergence, speciation and species persistence is of great interest to molecular ecology, especially for species-rich radiations inhabiting the world's biodiversity hotspots. The toolbox of population genomics holds great promise for addressing these key issues, especially if genomic data are analysed within a spatially and ecologically explicit context. We have studied the earliest stages of the divergence continuum in the Restionaceae, a species-rich and ecologically important plant family of the Cape Floristic Region (CFR) of South Africa, using the widespread CFR endemic Restio capensis (L.) H.P. Linder & C.R. Hardy as an example. We studied diverging populations of this morphotaxon for plastid DNA sequences and >14 400 nuclear DNA polymorphisms from Restriction site Associated DNA (RAD) sequencing and analysed the results jointly with spatial, climatic and phytogeographic data, using a Bayesian generalized linear mixed modelling (GLMM) approach. The results indicate that population divergence across the extreme environmental mosaic of the CFR is mostly driven by isolation by environment (IBE) rather than isolation by distance (IBD) for both neutral and non-neutral markers, consistent with genome hitchhiking or coupling effects during early stages of divergence. Mixed modelling of plastid DNA and single divergent outlier loci from a Bayesian genome scan confirmed the predominant role of climate and pointed to additional drivers of divergence, such as drift and ecological agents of selection captured by phytogeographic zones. Our study demonstrates the usefulness of population genomics for disentangling the effects of IBD and IBE along the divergence continuum often found in species radiations across heterogeneous ecological landscapes.

Extreme Precipitation Modelling Using Geostatistics and Machine Learning Algorithms

The paper presents an approach for mapping of precipitation data. The main goal is to perform spatial predictions and simulations of precipitation fields using geostatistical methods (ordinary kriging, kriging with external drift) as well as machine learning algorithms (neural networks). More practically, the objective is to reproduce simultaneously both the spatial patterns and the extreme values. This objective is best reached by models integrating geostatistics and machine learning algorithms. To demonstrate how such models work, two case studies have been considered: first, a 2-day accumulation of heavy precipitation and second, a 6-day accumulation of extreme orographic precipitation. The first example is used to compare the performance of two optimization algorithms (conjugate gradients and Levenberg-Marquardt) of a neural network for the reproduction of extreme values. Hybrid models, which combine geostatistical and machine learning algorithms, are also treated in this context. The second dataset is used to analyze the contribution of radar Doppler imagery when used as external drift or as input in the models (kriging with external drift and neural networks). Model assessment is carried out by comparing independent validation errors as well as analyzing data patterns.

Break zones in the distributions of alleles and species in alpine plants

Aim  We test for the congruence between allele-based range boundaries (break zones) in silicicolous alpine plants and species-based break zones in the silicicolous flora of the European Alps. We also ask whether such break zones coincide with areas of large elevational variation.Location  The European Alps.Methods  On a regular grid laid across the entire Alps, we determined areas of allele- and species-based break zones using respective clustering algorithms, identifying discontinuities in cluster distributions (breaks), and quantifying integrated break densities (break zones). Discontinuities were identified based on the intra-specific genetic variation of 12 species and on the floristic distribution data from 239 species, respectively. Coincidence between the two types of break zones was tested using Spearman's correlation. Break zone densities were also regressed on topographical complexity to test for the effect of elevational variation.Results  We found that two main break zones in the distribution of alleles and species were significantly correlated. Furthermore, we show that these break zones are in topographically complex regions, characterized by massive elevational ranges owing to high mountains and deep glacial valleys. We detected a third break zone in the distribution of species in the eastern Alps, which is not correlated with topographic complexity, and which is also not evident from allelic distribution patterns. Species with the potential for long-distance dispersal tended to show larger distribution ranges than short-distance dispersers.Main conclusions  We suggest that the history of Pleistocene glaciations is the main driver of the congruence between allele-based and species-based distribution patterns, because occurrences of both species and alleles were subject to the same processes (such as extinction, migration and drift) that shaped the distributions of species and genetic lineages. Large elevational ranges have had a profound effect as a dispersal barrier for alleles during post-glacial immigration. Because plant species, unlike alleles, cannot spread via pollen but only via seed, and thus disperse less effectively, we conclude that species break zones are maintained over longer time spans and reflect more ancient patterns than allele break zones.Conny Thiel-Egenter and Nadir Alvarez contributed equally to this paper and are considered joint first authors.

Invasion and fixation of sex-reversal genes.

We simulated a meta-population with random dispersal among demes but local mating within demes to investigate conditions under which a dominant female-determining gene W, with no individual selection advantage, can invade and become fixed in females, changing the population from male to female heterogamety. Starting with one mutant W in a single deme, the interaction of sex ratio selection and random genetic drift causes W to be fixed among females more often than a comparable neutral mutation with no influence on sex determination, even when YY males have slightly reduced viability. Meta-population structure and interdeme selection can also favour the fixation of W. The reverse transition from female to male heterogamety can also occur with higher probability than for a comparable neutral mutation. These results help to explain the involvement of sex-determining genes in the evolution of sex chromosomes and in sexual selection and speciation.

An alternative plate tectonic model for the Palaeozoic-Early Mesozoic Palaeotethyan evolution of Southeast Asia (northern Thailand-Myanmar)

An alternative model for the geodynamic evolution of Southeast Asia is proposed and inserted in a modern plate tectonic model. The reconstruction methodology is based on dynamic plate boundaries, constrained by data such as spreading rates and subduction velocities; in this way it differs from classical continental drift models proposed so far. The different interpretations about the location of the Palaeotethys suture in Thailand are revised, the Tertiary Mae Yuam fault is seen as the emplacement of the suture. East of the suture we identify an Indochina derived terrane for which we keep the name Shan-Thai, formerly used to identify the Cimmerian block present in Southeast Asia, now called Sibumasu. This nomenclatural choice was made on the basis of the geographic location of the terrane (Eastern Shan States in Burma and Central Thailand) and in order not to introduce new confusing terminology. The closure of the Eastern Palaeotethys is related to a southward subduction of the ocean, that triggered the Eastern Neotethys to open as a back-arc, due to the presence of Late Carboniferous-Early Permian arc magmatism in Mergui (Burma) and in the Lhasa block (South Tibet), and to the absence of arc magmatism of the same age East of the suture. In order to explain the presence of Carboniferous-Early Permian and Permo-Triassic volcanic arcs in Cambodia, Upper Triassic magmatism in Eastern Vietnam and Lower Permian-Middle Permian arc volcanites in Western Sumatra, we introduce the Orang Laut terranes concept. These terranes were detached from Indochina and South China during back-arc opening of the Poko-Song Ma system, due to the westward subduction of the Palaeopacific. This also explains the location of the Cathaysian West Sumatra block to the West of the Cimmerian Sibumasu block.

Genetic diversity in Brazilian populations of Aedes albopictus

Random amplified polymorphic DNA (RAPD) analysis technique was undertaken in Aedes albopictus populations from three states in Brazil, Rio de Janeiro (RJ), Minas Gerais (MG) and Pernambuco (PE), to estimate the level of genetic variability and levels of genetic exchange between populations. Allele and genotype frequencies were measured on 47 RAPD loci. Average observed heterozigosity (Ho) ranged from 0.282 in MG to 0.355 in Casa Forte (PE) population. Genetic distances estimates indicated that RJ and MG were more genetically similar than populations from PE. Genetic variation observed in local Brazilian populations was attributed to genetic drift associated with restricted gene flow in recently established populations.

Superiority of prone position in free-breathing 3D coronary MRA in patients with coronary disease.

Navigator-gated and corrected 3D coronary MR angiography (MRA) allows submillimeter image acquisition during free breathing. However, cranial diaphragmatic drift and relative phase shifts of chest-wall motion are limiting factors for image quality and scanning duration. We hypothesized that image acquisition in the prone position would minimize artifacts related to chest-wall motion and suppress diaphragmatic drift. Twelve patients with radiographically-confirmed coronary artery disease and six healthy adult volunteers were studied in both the prone and the supine position during free-breathing navigator-gated and corrected 3D coronary MRA. Image quality and the diaphragmatic positions were objectively compared. In the prone position, there was a 36% improvement in signal-to-noise ratio (SNR; 15.5 +/- 2.7 vs. 11.4 +/- 2.6; P < 0.01) and a 34% improvement in CNR (12.5 +/- 3.3 vs. 9.3 +/- 2.5, P < 0.01). The prone position also resulted in a 17% improvement in coronary vessel definition (P < 0.01). Cranial end-expiratory diaphragmatic drift occurred less frequently in the prone position (23% +/- 17% vs. 40% +/- 26% supine; P <0.05), and navigator efficiency was higher. Prone coronary MRA results in improved SNR and CNR with enhanced coronary vessel definition. Cranial end-expiratory diaphragmatic drift also was reduced, and navigator efficiency was enhanced. When feasible, prone imaging is recommended for free-breathing coronary MRA.

Multivariate QST-FST comparisons: a neutrality test for the evolution of the g matrix in structured populations.

Neutrality tests in quantitative genetics provide a statistical framework for the detection of selection on polygenic traits in wild populations. However, the existing method based on comparisons of divergence at neutral markers and quantitative traits (Q(st)-F(st)) suffers from several limitations that hinder a clear interpretation of the results with typical empirical designs. In this article, we propose a multivariate extension of this neutrality test based on empirical estimates of the among-populations (D) and within-populations (G) covariance matrices by MANOVA. A simple pattern is expected under neutrality: D = 2F(st)/(1 - F(st))G, so that neutrality implies both proportionality of the two matrices and a specific value of the proportionality coefficient. This pattern is tested using Flury's framework for matrix comparison [common principal-component (CPC) analysis], a well-known tool in G matrix evolution studies. We show the importance of using a Bartlett adjustment of the test for the small sample sizes typically found in empirical studies. We propose a dual test: (i) that the proportionality coefficient is not different from its neutral expectation [2F(st)/(1 - F(st))] and (ii) that the MANOVA estimates of mean square matrices between and among populations are proportional. These two tests combined provide a more stringent test for neutrality than the classic Q(st)-F(st) comparison and avoid several statistical problems. Extensive simulations of realistic empirical designs suggest that these tests correctly detect the expected pattern under neutrality and have enough power to efficiently detect mild to strong selection (homogeneous, heterogeneous, or mixed) when it is occurring on a set of traits. This method also provides a rigorous and quantitative framework for disentangling the effects of different selection regimes and of drift on the evolution of the G matrix. We discuss practical requirements for the proper application of our test in empirical studies and potential extensions.

Biocontrol by phenazine-1-carboxamide-producing Pseudomonas chlororaphis PCL1391of tomato root rot caused by Fusarium oxysporum f. sp. radicis-lycopersici

Seventy bacterial isolates from the rhizosphere of tomato were screened for antagonistic activity against the tomato foot and root rot-causing fungal pathogen Fusarium oxysporum f. sp. radicis-lycopersici. One isolate, strain PCL1391, appeared to be an efficient colonizer of tomato roots and an excellent biocontrol strain in an F. oxysporum/tomato test system. Strain PCL1391 was identified as Pseudomonas chlororaphis and further characterization showed that it produces a broad spectrum of antifungal factors (AFFs), including a hydrophobic compound, hydrogen cyanide, chitinase(s), and protease(s). Through mass spectrometry and nuclear magnetic resonance, the hydrophobic compound was identified as phenazine-1-carboxamide (PCN). We have studied the production and action of this AFF both in vitro and in vivo. Using a PCL1391 transposon mutant, with a lux reporter gene inserted in the phenazine biosynthetic operon (phz), we showed that this phenazine biosynthetic mutant was substantially decreased in both in vitro antifungal activity and biocontrol activity. Moreover, with the same mutant it was shown that the phz biosynthetic operon is expressed in the tomato rhizosphere. Comparison of the biocontrol activity of the PCN-producing strain PCL1391 with those of phenazine-1-carboxylic acid (PCA)-producing strains P. fluorescens 2-79 and P. aureofaciens 30-84 showed that the PCN-producing strain is able to suppress disease in the tomato/F. oxysporum system, whereas the PCA-producing strains are not. Comparison of in vitro antifungal activity of PCN and PCA showed that the antifungal activity of PCN was at least 10 times higher at neutral pH, suggesting that this may contribute to the superior biocontrol performance of strain PCL1391 in the tomato/F. oxysporum system.

The spread of incompatibility-inducing parasites in sub-divided host populations.

BACKGROUND:Maternally transmitted symbionts have evolved a variety of ways to promote their spread through host populations. One strategy is to hamper the reproduction of uninfected females by a mechanism called cytoplasmic incompatibility (CI). CI occurs in crosses between infected males and uninfected females and leads to partial to near-complete infertility. CI-infections are under positive frequency-dependent selection and require genetic drift to overcome the range of low frequencies where they are counter-selected. Given the importance of drift, population sub-division would be expected to facilitate the spread of CI. Nevertheless, a previous model concluded that variance in infection between competing groups of breeding individuals impedes the spread of CI.RESULTS:In this paper we derive a model on the spread of CI-infections in populations composed of demes linked by restricted migration. Our model shows that population sub-division facilitates the invasion of CI. While host philopatry (low migration) favours the spread of infection, deme size has a non-monotonous effect, with CI-invasion being most likely at intermediate deme size. Individual-based simulations confirm these predictions and show that high levels of local drift speed up invasion but prevent high levels of prevalence across the entire population. Additional simulations with sex-specific migration rates further show that low migration rates of both sexes are required to facilitate the spread of CI.CONCLUSION:Our analyses show that population structure facilitates the invasion of CI-infections. Since some level of sub-division is likely to occur in most natural populations, our results help to explain the high incidence of CI-infections across species of arthropods. Furthermore, our work has important implications for the use of CI-systems in order to genetically modify natural populations of disease vectors.

Quantification of the genetic change in the transition of Rhodnius pallescens Barber, 1932 (Hemiptera: Reduviidae) from field to laboratory

Previous studies have reported genetic differences between wild-caught sylvatic, domestic and laboratory pop-ulations of several Triatominae species. The differences between sylvatic and laboratory colonies parallel are similar to the differences observed between sylvatic and domestic populations. Laboratory colonies are frequently used as references for field populations, but the consequences of founder events on the genetic makeup of laboratory or domestic populations are rarely quantified. Our goal was to quantify the genetic change in Rhodnius pallescens populations artificially submitted to founder effects via laboratory colonization. We compared the genetic makeup of two sylvatic populations and their laboratory descendants using a panel of 10 microsatellite markers. Both sylvatic populations were initially collected from palm trees, but the colonies differed in the number of founder insects and amount of time kept in the laboratory. We evaluated allelic polymorphism, differences between expected and observed heterozygosity, estimates of population differentiation (Fst) and inbreeding (Fis, Fit) and cluster analyses based on Nei's distances. We found a unique genetic structure for each sample population, with significant differentiation between the field insects and each of the laboratory generations. These analyses showed strong founder effects and showed that genetic drift had led to a genetic equilibrium over several generations of isolation. Our results suggest that laboratory colonies of R. pallescens have a different genetic structure than their wild relatives and similar processes likely affect other Triatominae laboratory stocks.

Evolutionary implications of a high selfing rate in the freshwater snail Lymnaea truncatula.

Self-compatible hermaphroditic organisms that mix self-fertilization and outcrossing are of great interest for investigating the evolution of mating systems. We investigate the evolution of selfing in Lymnaea truncatula, a self-compatible hermaphroditic freshwater snail. We first analyze the consequences of selfing in terms of genetic variability within and among populations and then investigate how these consequences along with the species ecology (harshness of the habitat and parasitism) might govern the evolution of selfing. Snails from 13 localities (classified as temporary or permanent depending on their water availability) were sampled in western Switzerland and genotyped for seven microsatellite loci. F(IS) (estimated on adults) and progeny array analyses (on hatchlings) provided similar selfing rate estimates of 80%. Populations presented a low polymorphism and were highly differentiated (F(ST) = 0.58). Although the reproductive assurance hypothesis would predict higher selfing rate in temporary populations, no difference in selfing level was observed between temporary and permanent populations. However, allelic richness and gene diversity declined in temporary habitats, presumably reflecting drift. Infection levels varied but were not simply related to either estimated population selfing rate or to differences in heterozygosity. These findings and the similar selfing rates estimated for hatchlings and adults suggest that within-population inbreeding depression is low in L. truncatula.

Vision-based localization of an underwater robot in a structured environment

This paper presents a vision-based localization approach for an underwater robot in a structured environment. The system is based on a coded pattern placed on the bottom of a water tank and an onboard down looking camera. Main features are, absolute and map-based localization, landmark detection and tracking, and real-time computation (12.5 Hz). The proposed system provides three-dimensional position and orientation of the vehicle along with its velocity. Accuracy of the drift-free estimates is very high, allowing them to be used as feedback measures of a velocity-based low-level controller. The paper details the localization algorithm, by showing some graphical results, and the accuracy of the system