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.

Genetic Structure of Triatoma sordida (Hemiptera: Reduviidae) Domestic Populations from Bolivia: Application on Control Interventions

The genetic population of Triatoma sordida group 1, a secondary vector of Chagas disease in Bolivia, was studied by multi-locus enzyme electrophoresis. A total of 253 nymphal and adult specimens collected from seven neighbouring localities in the Velasco Province, Department of Santa Cruz, were processed. The relatively low genetic variability was confirmed for this species (rate of polymorphism: 0.20). The absence of genetic disequilibrium detected within the seven localities was demonstrated. A geographical structuration appears between localities with distances greater than 20 km apart. Although T. sordida presents a relatively reduced dispersive capacity, its panmictic unit is wider than compared with T. infestans. Genetic distances between T. sordida populations were correlated with geographic distance. Gene flow between geographic populations of T. sordida provides an efficient framework for effective vigilance and control protocols.

Penetrance of marked cognitive impairment in older male carriers of the FMR1 gene premutation.

BACKGROUND: Male carriers of the FMR1 premutation are at risk of developing the fragile X-associated tremor/ataxia syndrome (FXTAS), a newly recognised and largely under-diagnosed late onset neurodegenerative disorder. Patients affected with FXTAS primarily present with cerebellar ataxia and intention tremor. Cognitive decline has also been associated with the premutation, but the lack of data on its penetrance is a growing concern for clinicians who provide genetic counselling. METHODS: The Mattis Dementia Rating Scale (MDRS) was administered in a double blind fashion to 74 men aged 50 years or more recruited from fragile X families (35 premutation carriers and 39 intrafamilial controls) regardless of their clinical manifestation. Based on previous publications, marked cognitive impairment was defined by a score <or=123 on the MDRS. RESULTS: Both logistic and survival models confirmed that in addition to age and education level, premutation size plays a significant (p<0.01 and p<0.03 for logistic and survival model, respectively) role in cognitive impairment. The estimated penetrance of marked cognitive impairment in our sample (adjusted for the mean age 63.4 years and mean education level 9.7 years) for midsize/large (70-200 CGG) and small (55-69 CGG) premutation alleles was 33.3% (relative risk (RR) 6.5; p = 0.01) and 5.9% (RR 1.15; p = 0.9) respectively. Penetrance in the control group was 5.1%. CONCLUSIONS: Male carriers of midsize to large premutation alleles had a sixfold increased risk of developing cognitive decline and the risk increases with allele size. In addition, it was observed that cognitive impairment may precede motor symptoms. These data provide guidance for genetic counselling although larger samples are required to refine these estimates.

Sylvatic population of Triatoma infestans from the Bolivian Chaco: from field collection to characterization

A sylvatic Triatoma infestans DM (dark morph) population detected in the Bolivian Chaco was characterized and compared with various domestic ones. The degree of differentiation of DM was clearly within the T. infestans intra-specific level. Nevertheless marked chromatic and morphometric differences as well as differences in antennal pattern, chromosome banding and randomly amplified polymorphic DNA support the hypothesis of a distinct population. Continuous exchange of insects between wild and domestic habitats seems unlikely in the Chaco.

Non-nest mate discrimination and clonal colony structure in the parthenogenetic ant Cerapachys biroi

Understanding the interplay between cooperation and conflict in social groups is a major goal of biology. One important factor is genetic relatedness, and animal societies are usually composed of related but genetically different individuals, setting the stage for conflicts over reproductive allocation. Recently, however, it has been found that several ant species reproduce predominantly asexually. Although this can potentially give rise to clonal societies, in the few well-studied cases, colonies are often chimeric assemblies of different genotypes, due to worker drifting or colony fusion. In the ant Cerapachys biroi, queens are absent and all individuals reproduce via thelytokous parthenogenesis, making this species an ideal study system of asexual reproduction and its consequences for social dynamics. Here, we show that colonies in our study population on Okinawa, Japan, recognize and effectively discriminate against foreign workers, especially those from unrelated asexual lineages. In accord with this finding, colonies never contained more than a single asexual lineage and average pairwise genetic relatedness within colonies was extremely high (r = 0.99). This implies that the scope for social conflict in C. biroi is limited, with unusually high potential for cooperation and altruism.

Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study.

The identification of genetically homogeneous groups of individuals is a long standing issue in population genetics. A recent Bayesian algorithm implemented in the software STRUCTURE allows the identification of such groups. However, the ability of this algorithm to detect the true number of clusters (K) in a sample of individuals when patterns of dispersal among populations are not homogeneous has not been tested. The goal of this study is to carry out such tests, using various dispersal scenarios from data generated with an individual-based model. We found that in most cases the estimated 'log probability of data' does not provide a correct estimation of the number of clusters, K. However, using an ad hoc statistic DeltaK based on the rate of change in the log probability of data between successive K values, we found that STRUCTURE accurately detects the uppermost hierarchical level of structure for the scenarios we tested. As might be expected, the results are sensitive to the type of genetic marker used (AFLP vs. microsatellite), the number of loci scored, the number of populations sampled, and the number of individuals typed in each sample.

20 years of eosinophilic esophagitis in Olten county: population-based, prospective documentation of an increasing incidence

Background and Aims: Eosinophilic Esophagitis (EoE) is reported with increasing frequency over the last two decades. However, it is still unknown whether this reflects a true increase in incidence or just an increased awareness by gastroenterologists. Therefore, we evaluated the incidence and cumulative prevalence of EoE in Olten county over the last 20 years. Methods: Olten county is an area of approximately 91,000 inhabitants without pronounced demographic changes in the last two decades. EoE evaluation is based upon two gastroenterology centers and one pathology center. No public programs for increased EoE awareness were implemented in this region. All EoE patients diagnosed from 1989 to 2009 were entered prospectively into the Olten county database. Results: Fourty-six patients (76% males, mean age 41±16 yrs) were diagnosed with EoE from 1989 to 2009. Ninety-four percent presented with dysphagia. In 70% of the patients concomitant allergies were found. The number of upper endoscopies per year was stable during the entire observation period. An average annual incidence rate of 2/100,000 was found (range 0-8) with a marked increase in the period from 2001 to 2009. A current cumulative EoE prevalence of 43/100,000 inhabitants was calculated. The mean diagnostic delay (time from first symptoms to diagnosis) was 4.3 years from 1989 to 1998 and 4.8 years from 1999 to 2009. Conclusions: Over the last 20 years, a significant increase in EoE incidence was found in a stable indicator region of Switzerland. The constant rate of upper endoscopies, the constant diagnostic delay, as well as the lack of EoE awareness programs in Olten county indicate a true increase in EoE incidence.