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.

Frequency-dependent sexual selection with respect to offspring fitness returns is consistent with predictions from rock-paper-scissors dynamics in the European common lizard

Genetic polymorphism can be maintained over time by negative frequency-dependent (FD) selection induced by Rock-paper-scissors (RPS) social systems. RPS games produce cyclic dynamics, and have been suggested to exist in lizards, insects, isopods, plants, and bacteria. Sexual selection is predicted to accentuate the survival of the future progeny during negative FD survival selection. More specifically, females are predicted to select mates that produce progeny genotypes that exhibit highest survival during survival selection imposed by adult males. However, no empirical evidence demonstrates the existence of FD sexual selection with respect to fitness payoffs of genetic polymorphisms. Here we tested this prediction using the common lizard Zootoca vivipara, a species with three male color morphs (orange, white, yellow) that exhibit morph frequency cycles. In a first step we tested the congruence of the morph frequency change with the predicted change in three independent populations, differing in male color morph frequency and state of the FD morph cycle. Thereafter we ran standardized sexual selection assays in which we excluded alternative mechanisms that potentially induce negative FD selection, and we quantified inter-sexual behavior. The patterns of sexual selection and the observed behavior were in line with context-dependent female mate choice and male behavior played a minor role. Moreover, the strength of the sexual selection was within the magnitude of selection required to produce the observed 3-4-year and 6-8 year morph frequency cycles at low and high altitudes, respectively. In summary, the study provides the first experimental evidence that underpins the crucial assumption of the RPS games suggested to exist in lizards, insects, isopods, and plants; namely, that sexual selection produces negative-FD selection. This indicates that sexual selection, in our study exert by females, might be a crucial driver of the maintenance of genetic polymorphisms.

Selection of Beauveria bassiana and Metarhizium anisopliae Isolates to Control Triatoma infestans

Twenty three isolates of Beauveria bassiana and 13 isolates of Metarhizium anisopliae were tested on third instar nymphs of Triatoma infestans, a serious vector of Chagas disease. Pathogenicity tests at saturated humidity showed that this insect is very susceptible to fungal infection. At lower relative humidity (50%), conditions expected in the vector microhabitat, virulence was significantly different among isolates. Cumulative mortality 15 days after treatment varied from 17.5 to 97.5%, and estimates of 50% survival time varied from 6 to 11 days. Maintaining lower relative humidity, four B. bassiana and two M. anisopliae isolates were selected for analysis of virulence at different conidial concentrations and temperatures. Lethal concentrations sufficient to kill 50% of insects (LC50) varied from 7.1x105 to 4.3x106 conidia/ml, for a B. bassiana isolate (CG 14) and a M. anisopliae isolate (CG 491) respectively. Most isolates, particularly B. bassiana isolates CG 24 and CG 306, proved to be more virulent at 25 and 30°C, compared to 15 and 20°C. The differential virulence at 50% humidity observed among some B. bassiana isolates was not correlated to phenetic groups in cluster analysis of RAPD markers. In fact, the B. bassiana isolates analyzed presented a high homogeneity (> 73% similarity).

Identifying genetic signatures of selection in a non-model species, alpine gentian (Gentiana nivalis L.), using a landscape genetic approach

It is generally accepted that most plant populations are locally adapted. Yet, understanding how environmental forces give rise to adaptive genetic variation is a challenge in conservation genetics and crucial to the preservation of species under rapidly changing climatic conditions. Environmental variation, phylogeographic history, and population demographic processes all contribute to spatially structured genetic variation, however few current models attempt to separate these confounding effects. To illustrate the benefits of using a spatially-explicit model for identifying potentially adaptive loci, we compared outlier locus detection methods with a recently-developed landscape genetic approach. We analyzed 157 loci from samples of the alpine herb Gentiana nivalis collected across the European Alps. Principle coordinates of neighbor matrices (PCNM), eigenvectors that quantify multi-scale spatial variation present in a data set, were incorporated into a landscape genetic approach relating AFLP frequencies with 23 environmental variables. Four major findings emerged. 1) Fifteen loci were significantly correlated with at least one predictor variable (R (adj) (2) > 0.5). 2) Models including PCNM variables identified eight more potentially adaptive loci than models run without spatial variables. 3) When compared to outlier detection methods, the landscape genetic approach detected four of the same loci plus 11 additional loci. 4) Temperature, precipitation, and solar radiation were the three major environmental factors driving potentially adaptive genetic variation in G. nivalis. Techniques presented in this paper offer an efficient method for identifying potentially adaptive genetic variation and associated environmental forces of selection, providing an important step forward for the conservation of non-model species under global change.

The Body Action and Posture Coding System (BAP): Development and Reliability

Several methods are available for coding body movement in nonverbal behavior research, but there is no consensus on a reliable coding system that can be used for the study of emotion expression. Adopting an integrative approach, we developed a new method, the Body Action and Posture (BAP) coding system, for the time-aligned micro description of body movement on an anatomical level (different articulations of body parts), a form level (direction and orientation of movement), and a functional level (communicative and self-regulatory functions). We applied the system to a new corpus of acted emotion portrayals, examined its comprehensiveness and demonstrated intercoder reliability at three levels: a) occurrence, b) temporal precision and c) segmentation. We discuss issues for further validation and propose some research applications.

Improving the performance of positive selection inference by filtering unreliable alignment regions.

Errors in the inferred multiple sequence alignment may lead to false prediction of positive selection. Recently, methods for detecting unreliable alignment regions were developed and were shown to accurately identify incorrectly aligned regions. While removing unreliable alignment regions is expected to increase the accuracy of positive selection inference, such filtering may also significantly decrease the power of the test, as positively selected regions are fast evolving, and those same regions are often those that are difficult to align. Here, we used realistic simulations that mimic sequence evolution of HIV-1 genes to test the hypothesis that the performance of positive selection inference using codon models can be improved by removing unreliable alignment regions. Our study shows that the benefit of removing unreliable regions exceeds the loss of power due to the removal of some of the true positively selected sites.

Selection of Trypanosoma cruzi clonal genotypes (clonet 20 and 39) isolated from Bolivian triatomines following subculture in liquid medium

Previous studies showed that two groups of Trypanosoma cruzi clonal genotypes named clonet 20 and clonet 39 were predominant in Triatoma infestans, the unique vector of Chagas disease in Bolivia. These groups of clones correspond to distinct genetic clusters. These clonets were detected in T. infestans and Rhodnius pictipes fecal samples before isolation and after culture by kDNA PCR (polymerase chain rreaction) and hybridization of the amplified products with clonet specific kDNA probes named 20 and 39 as previously reported. Forty eight T. infestans and three R. pictipes infected insects captured at random in different Bolivian departments were proceeded. As previously reported the direct identification of the two major clonets in fecal samples allowed the detection of abundant mixed infections: 41% in the original sample, however after culture, only 6% of mixed infections were detected. Among the 21 parasite stocks isolated from digestive tracts where mixed infections were initially detected (clonet 20 + 39) clonet 20 alone was detected in 81% of them. This result clearly showed that the culture step selected clonet 20 parasites over those belonging to clonet 39. The taxonomic status of the isolated stocks was also confirmed by isoenzyme typing, and correlation was observed between clustering topology and hybridization patterns with the probes 20 and 39.

Thymic selection generates a large T cell pool recognizing a self-peptide in humans.

The low frequency of self-peptide-specific T cells in the human preimmune repertoire has so far precluded their direct evaluation. Here, we report an unexpected high frequency of T cells specific for the self-antigen Melan-A/MART-1 in CD8 single-positive thymocytes from human histocompatibility leukocyte antigen-A2 healthy individuals, which is maintained in the peripheral blood of newborns and adults. Postthymic replicative history of Melan-A/MART-1-specific CD8 T cells was independently assessed by quantifying T cell receptor excision circles and telomere length ex vivo. We provide direct evidence that the large T cell pool specific for the self-antigen Melan-A/MART-1 is mostly generated by thymic output of a high number of precursors. This represents the only known naive self-peptide-specific T cell repertoire directly accessible in humans.

Mosquitocidal bacterial toxins: diversity, mode of action and resistance phenomena

Bacteria active against dipteran larvae (mosquitoes and black flies) include a wide variety of Bacillus thuringiensis and B. sphaericus strains, as well as isolates of Brevibacillus laterosporus and Clostridium bifermentans. All display different spectra and levels of activity correlated with the nature of the toxins, mainly produced during the sporulation process. This paper describes the structure and mode of action of the main mosquitocidal toxins, in relationship with their potential use in mosquito and/or black fly larvae control. Investigations with laboratory and field colonies of mosquitoes that have become highly resistant to the B. sphaericus Bin toxin have shown that several mechanisms of resistance are involved, some affecting the toxin/receptor binding step, others unknown.

Strategy selection during exploratory behavior: sex differences

This study was designed to assess sex-related differences in the selection of an appropriate strategy when facing novelty. A simple visuo-spatial task was used to investigate exploratory behavior as a specific response to novelty. The exploration task was followed by a visual discrimination task, and the responses were analyzed using signal detection theory. During exploration women selected a local searching strategy in which the metric distance between what is already known and what is unknown was reduced, whereas men adopted a global strategy based on an approximately uniform distribution of choices. Women's exploratory behavior gives rise to a notion of a secure base warranting a sense of safety while men's behavior does not appear to be influenced by risk. This sex-related difference was interpreted as a difference in beliefs concerning the likelihood of uncertain events influencing risk evaluation. Keywords: exploration, spontaneous strategies, sex differences, decision-making.

Public sector rationing and private sector selection [paper accepted and forthcoming in Journal of Public Economic Theory]

We study the interaction between nonprice public rationing and prices in the private market. Under a limited budget, the public supplier uses a rationing policy. A private firm may supply the good to those consumers who are rationed by the public system. Consumers have different amounts of wealth, and costs of providing the good to them vary. We consider two regimes. First, the public supplier observes consumers' wealth information; second, the public supplier observes both wealth and cost information. The public supplier chooses a rationing policy, and, simultaneously, the private firm, observing only cost but not wealth information, chooses a pricing policy. In the first regime, there is a continuum of equilibria. The Pareto dominant equilibrium is a means-test equilibrium: poor consumers are supplied while rich consumers are rationed. Prices in the private market increase with the budget. In the second regime, there is a unique equilibrium. This exhibits a cost-effectiveness rationing rule; consumers are supplied if and only if their costbenefit ratios are low. Prices in the private market do not change with the budget. Equilibrium consumer utility is higher in the cost-effectiveness equilibrium than the means-test equilibrium [Authors]

Vertebrate conserved non coding DNA regions have a high persistence length and a short persistence time.

BACKGROUND: The comparison of complete genomes has revealed surprisingly large numbers of conserved non-protein-coding (CNC) DNA regions. However, the biological function of CNC remains elusive. CNC differ in two aspects from conserved protein-coding regions. They are not conserved across phylum boundaries, and they do not contain readily detectable sub-domains. Here we characterize the persistence length and time of CNC and conserved protein-coding regions in the vertebrate and insect lineages. RESULTS: The persistence length is the length of a genome region over which a certain level of sequence identity is consistently maintained. The persistence time is the evolutionary period during which a conserved region evolves under the same selective constraints.Our main findings are: (i) Insect genomes contain 1.60 times less conserved information than vertebrates; (ii) Vertebrate CNC have a higher persistence length than conserved coding regions or insect CNC; (iii) CNC have shorter persistence times as compared to conserved coding regions in both lineages. CONCLUSION: Higher persistence length of vertebrate CNC indicates that the conserved information in vertebrates and insects is organized in functional elements of different lengths. These findings might be related to the higher morphological complexity of vertebrates and give clues about the structure of active CNC elements.Shorter persistence time might explain the previously puzzling observations of highly conserved CNC within each phylum, and of a lack of conservation between phyla. It suggests that CNC divergence might be a key factor in vertebrate evolution. Further evolutionary studies will help to relate individual CNC to specific developmental processes.

Cefotaxime acts synergistically with levofloxacin in experimental meningitis due to penicillin-resistant pneumococci and prevents selection of levofloxacin-resistant mutants in vitro.

Cefotaxime, given in two doses (each 100 mg/kg of body weight), produced a good bactericidal activity (-0.47 Deltalog(10) CFU/ml. h) which was comparable to that of levofloxacin (-0.49 Deltalog(10) CFU/ml. h) against a penicillin-resistant pneumococcal strain WB4 in experimental meningitis. Cefotaxime combined with levofloxacin acted synergistically (-1.04 Deltalog(10) CFU/ml. h). Synergy between cefotaxime and levofloxacin was also demonstrated in vitro in time killing assays and with the checkerboard method for two penicillin-resistant strains (WB4 and KR4). Using in vitro cycling experiments, the addition of cefotaxime in sub-MIC concentrations (one-eighth of the MIC) drastically reduced levofloxacin-induced resistance in the same two strains (64-fold increase of the MIC of levofloxacin after 12 cycles versus 2-fold increase of the MIC of levofloxacin combined with cefotaxime). Mutations detected in the genes encoding topoisomerase IV (parC and parE) and gyrase (gyrA and gyrB) confirmed the levofloxacin-induced resistance in both strains. Addition of cefotaxime in low doses was able to suppress levofloxacin-induced resistance.