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.

Association between circulating cytokine levels, diabetes and insulin resistance in a population-based sample (CoLaus study).

OBJECTIVE: The associations between inflammation, diabetes and insulin resistance remain controversial. Hence, we assessed the associations between diabetes, insulin resistance (using HOMA-IR) and metabolic syndrome with the inflammatory markers high-sensitive C-reactive protein (hs-CRP), interleukin-1 beta (IL-1β), interleukin-6 (IL-6) and tumour necrosis factor-α (TNF-α). DESIGN: Cross-sectional study. PARTICIPANTS: Two thousand eight hundred and eighty-four men and 3201 women, aged 35-75, participated in this study. METHODS: C-reactive protein was assessed by immunoassay and cytokines by multiplexed flow cytometric assay. In a subgroup of 532 participants, an oral glucose tolerance test (OGTT) was performed to screen for impaired glucose tolerance (IGT). RESULTS: IL-6, TNF-α and hs-CRP were significantly and positively correlated with fasting plasma glucose (FPG), insulin and HOMA-IR. Participants with diabetes had higher IL-6, TNF-α and hs-CRP levels than participants without diabetes; this difference persisted for hs-CRP after multivariate adjustment. Participants with metabolic syndrome had increased IL-6, TNF-α and hs-CRP levels; these differences persisted after multivariate adjustment. Participants in the highest quartile of HOMA-IR had increased IL-6, TNF-α and hs-CRP levels; these differences persisted for TNF-α and hs-CRP after multivariate adjustment. No association was found between IL-1β levels and all diabetes and insulin resistance markers studied. Finally, participants with IGT had higher hs-CRP levels than participants with a normal OGTT, but this difference disappeared after controlling for body mass index (BMI). CONCLUSION: We found that subjects with diabetes, metabolic syndrome and increased insulin resistance had increased levels of IL6, TNF-α and hs-CRP, while no association was found with IL-1β. The increased inflammatory state of subjects with IGT is partially explained by increased BMI.

Population pharmacokinetic modelling and evaluation of different dosage regimens for darunavir and ritonavir in HIV-infected individuals.

OBJECTIVES: Darunavir is a protease inhibitor that is administered with low-dose ritonavir to enhance its bioavailability. It is prescribed at standard dosage regimens of 600/100 mg twice daily in treatment-experienced patients and 800/100 mg once daily in naive patients. A population pharmacokinetic approach was used to characterize the pharmacokinetics of both drugs and their interaction in a cohort of unselected patients and to compare darunavir exposure expected under alternative dosage regimens. METHODS: The study population included 105 HIV-infected individuals who provided darunavir and ritonavir plasma concentrations. Firstly, a population pharmacokinetic analysis for darunavir and ritonavir was conducted, with inclusion of patients' demographic, clinical and genetic characteristics as potential covariates (NONMEM(®)). Then, the interaction between darunavir and ritonavir was studied while incorporating levels of both drugs into different inhibitory models. Finally, model-based simulations were performed to compare trough concentrations (Cmin) between the recommended dosage regimen and alternative combinations of darunavir and ritonavir. RESULTS: A one-compartment model with first-order absorption adequately characterized darunavir and ritonavir pharmacokinetics. The between-subject variability in both compounds was important [coefficient of variation (CV%) 34% and 47% for darunavir and ritonavir clearance, respectively]. Lopinavir and ritonavir exposure (AUC) affected darunavir clearance, while body weight and darunavir AUC influenced ritonavir elimination. None of the tested genetic variants showed any influence on darunavir or ritonavir pharmacokinetics. The simulations predicted darunavir Cmin much higher than the IC50 thresholds for wild-type and protease inhibitor-resistant HIV-1 strains (55 and 550 ng/mL, respectively) under standard dosing in >98% of experienced and naive patients. Alternative regimens of darunavir/ritonavir 1200/100 or 1200/200 mg once daily also had predicted adequate Cmin (>550 ng/mL) in 84% and 93% of patients, respectively. Reduction of darunavir/ritonavir dosage to 600/50 mg twice daily led to a 23% reduction in average Cmin, still with only 3.8% of patients having concentrations below the IC50 for resistant strains. CONCLUSIONS: The important variability in darunavir and ritonavir pharmacokinetics is poorly explained by clinical covariates and genetic influences. In experienced patients, treatment simplification strategies guided by drug level measurements and adherence monitoring could be proposed.

Dosing regimen of gentamicin in neonates: evaluation of current guidelines based on a large population pharmacokinetic analysis

Objectives: Gentamicin is among the most commonly prescribed antibiotics in newborns, but large interindividual variability in exposure levels exists. Based on a population pharmacokinetic analysis of a cohort of unselected neonates, we aimed to validate current dosing recommendations from a recent reference guideline (Neofax®). Methods: From 3039 concentrations collected in 994 preterm (median gestational age 32.3 weeks, range 24.2-36.5) and 455 term newborns, treated at the University Hospital of Lausanne between 2006 and 2011, a population pharmacokinetic analysis was performed with NONMEM®. Model-based simulations were used to assess the ability of dosing regimens to bring concentrations into targets: trough ≤ 1mg/L and peak ~ 8mg/L. Results: A two-compartment model best characterized gentamicin pharmacokinetics. Model parameters are presented in the table. Body weight, gestational age and postnatal age positively influence clearance, which decreases under dopamine administration. Body weight and gestational age influence the distribution volume. Model based simulations confirm that preterm infants need doses superior to 4 mg/kg, and extended dosage intervals, up to 48 hours for very preterm newborns, whereas most term newborns would achieve adequate exposure under 4 mg/kg q. 24 h. More than 90% of neonates would achieve trough concentrations below 2 mg/L and peaks above 6 mg/L following most recent guidelines. Conclusions: Simulated gentamicin exposure demonstrates good accordance with recent dosing recommendations for target concentration achievement.

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.

Population pharmacokinetics and clinical response for artemether-lumefantrine in pregnant and nonpregnant women with uncomplicated Plasmodium falciparum malaria in Tanzania.

Artemether-lumefantrine (AL) is the first-line treatment for uncomplicated malaria in the second and third trimesters of pregnancy. Its efficacy during pregnancy has recently been challenged due to altered pharmacokinetic (PK) properties in this vulnerable group. The aim of this study was to determine the PK profile of AL in pregnant and nonpregnant women and assess their therapeutic outcome. Thirty-three pregnant women and 22 nonpregnant women with malaria were treated with AL (80/480 mg) twice daily for 3 days. All patients provided five venous plasma samples for drug quantification at random times over 7 days. Inter- and intraindividual variability was assessed, and the effects of covariates were quantified using a nonlinear mixed-effects modeling approach (NONMEM). A one-compartment model with first-order absorption and elimination with linear metabolism from drug to metabolite fitted the data best for both arthemether (AM) and lumefantrine (LF) and their metabolites. Pregnancy status and diarrhea showed a significant influence on LF PK. The relative bioavailability of lumefantrine and its metabolism rate into desmethyl-lumefantrine were, respectively, 34% lower and 78% higher in pregnant women than in nonpregnant patients. The overall PCR-uncorrected treatment failure rates were 18% in pregnant women and 5% in nonpregnant women (odds ratio [OR] = 4.04; P value of 0.22). A high median day 7 lumefantrine concentration was significantly associated with adequate clinical and parasitological response (P = 0.03). The observed reduction in the relative bioavailability of lumefantrine in pregnant women may explain the higher treatment failure in this group, mostly due to lower posttreatment prophylaxis. Hence, a modified treatment regimen of malaria in pregnancy should be considered.

Alcohol drinking and cardiovascular risk in a population with high mean alcohol consumption.

Moderate alcohol consumption has been associated with lower coronary artery disease (CAD) risk. However, data on the CAD risk associated with high alcohol consumption are conflicting. The aim of this study was to examine the impact of heavier drinking on 10-year CAD risk in a population with high mean alcohol consumption. In a population-based study of 5,769 adults (aged 35 to 75 years) without cardiovascular disease in Switzerland, 1-week alcohol consumption was categorized as 0, 1 to 6, 7 to 13, 14 to 20, 21 to 27, 28 to 34, and > or =35 drinks/week or as nondrinkers (0 drinks/week), moderate (1 to 13 drinks/week), high (14 to 34 drinks/week), and very high (> or =35 drinks/week). Blood pressure and lipids were measured, and 10-year CAD risk was calculated according to the Framingham risk score. Seventy-three percent (n = 4,214) of the participants consumed alcohol; 16% (n = 909) were high drinkers and 2% (n = 119) very high drinkers. In multivariate analysis, increasing alcohol consumption was associated with higher high-density lipoprotein cholesterol (from a mean +/- SE of 1.57 +/- 0.01 mmol/L in nondrinkers to 1.88 +/- 0.03 mmol/L in very high drinkers); triglycerides (1.17 +/- 1.01 to 1.32 +/- 1.05 mmol/L), and systolic and diastolic blood pressure (127.4 +/- 0.4 to 132.2 +/- 1.4 mm Hg and 78.7 +/- 0.3 to 81.7 +/- 0.9 mm Hg, respectively) (all p values for trend <0.001). Ten-year CAD risk increased from 4.31 +/- 0.10% to 4.90 +/- 0.37% (p = 0.03) with alcohol use, with a J-shaped relation. Increasing wine consumption was more related to high-density lipoprotein cholesterol levels, whereas beer and spirits were related to increased triglyceride levels. In conclusion, as measured by 10-year CAD risk, the protective effect of alcohol consumption disappears in very high drinkers, because the beneficial increase in high-density lipoprotein cholesterol is offset by the increases in blood pressure levels.

Biocontrol strain Pseudomonas fluorescens CHA0 and its genetically modified derivative with enhanced biocontrol capability exert comparable effects on the structure of a Sinorhizobium meliloti population in a gnotobiotic system

The impact of biocontrol strain Pseudomonas fluorescens CHA0 and of its genetically modified, antibiotic-overproducing derivative CHA0/pME3424 on a reconstructed population of the plant-beneficial Sinorhizobium meliloti bacteria was assessed in gnotobiotic systems. In sterile soil, the final density of the reconstructed S. meliloti population decreased by more than one order of magnitude in the presence of either of the Pseudomonas strains when compared to a control without addition of P. fluorescens. Moreover, there was a change in the proportion of each individual S. meliloti strain within the population. Plant tests also revealed changes in the nodulating S. meliloti population in the presence of strains CHA0 or CHA0/pME3424. In both treatments one S. meliloti strain, f43, was significantly reduced in its root nodule occupancy. Analysis of alfalfa yields showed a slight but statistically significant increase in shoot dry weight when strain CHA0 was added to the reconstructed S. meliloti population whereas no such effect was observed with CHA0/pME3424.

Patterns of sex work contact among men in the general population of Switzerland, 1987-2000

OBJECTIVES: To analyse the prevalence of lifetime recourse to prostitution (LRP) among men in the general population of Switzerland from a trend and cohort perspective. METHODS: Using nine repeated representative cross-sectional surveys from 1987 to 2000, age-specific estimates of LRP were computed. Trends and period effect were analysed as the evolution of cross-sectional population estimates within age groups and overall. Cohort analysis relied on cohorts constructed from the 1989 survey and followed in subsequent waves. Age and cohort effects were modelled using logistic regression and non-parametric monotone regression. RESULTS: Whereas prevalence for the younger groups was found to be logically lower, there was no consistent increasing or decreasing trend over the years; there was no significant period effect. For the 17-30 year age group, the mean estimate over 1987-2000 was 11.5% (range 8.3 to 12.7%); for the 31-45 year group, the mean was 21.5% (range over 1989-2000 20.3 to 23.0%). Regarding cohort analysis, the prevalence of LRP was found to increase steeply in the youngest ages before reaching a plateau near the age of 40 years. At the age of 43 years, the prevalence was estimated to be 22.6% (95% CI 21.1% to 24.1%). CONCLUSIONS: The steep increase in the cohort-wise prevalence of LRP in younger ages calls for a concentration of prevention activities in young people. If the plateauing at approximately 40 years of age is not followed by a further increase later in life, which is not known, then consumers of paid sex would be repeat buyers only, a fact that should be taken into account by prevention.

Paper 3: EUROCAT data quality indicators for population-based registries of congenital anomalies.

The European Surveillance of Congenital Anomalies (EUROCAT) network of population-based congenital anomaly registries is an important source of epidemiologic information on congenital anomalies in Europe covering live births, fetal deaths from 20 weeks gestation, and terminations of pregnancy for fetal anomaly. EUROCAT's policy is to strive for high-quality data, while ensuring consistency and transparency across all member registries. A set of 30 data quality indicators (DQIs) was developed to assess five key elements of data quality: completeness of case ascertainment, accuracy of diagnosis, completeness of information on EUROCAT variables, timeliness of data transmission, and availability of population denominator information. This article describes each of the individual DQIs and presents the output for each registry as well as the EUROCAT (unweighted) average, for 29 full member registries for 2004-2008. This information is also available on the EUROCAT website for previous years. The EUROCAT DQIs allow registries to evaluate their performance in relation to other registries and allows appropriate interpretations to be made of the data collected. The DQIs provide direction for improving data collection and ascertainment, and they allow annual assessment for monitoring continuous improvement. The DQI are constantly reviewed and refined to best document registry procedures and processes regarding data collection, to ensure appropriateness of DQI, and to ensure transparency so that the data collected can make a substantial and useful contribution to epidemiologic research on congenital anomalies.

Anti-basal ganglia antibodies and Tourette's syndrome: a voxel-based morphometry and diffusion tensor imaging study in an adult population.

Anti-basal ganglia antibodies (ABGAs) have been suggested to be a hallmark of autoimmunity in Gilles de la Tourette's syndrome (GTS), possibly related to prior exposure to streptococcal infection. In order to detect whether the presence of ABGAs was associated with subtle structural changes in GTS, whole-brain analysis using independent sets of T(1) and diffusion tensor imaging MRI-based methods were performed on 22 adults with GTS with (n = 9) and without (n = 13) detectable ABGAs in the serum. Voxel-based morphometry analysis failed to detect any significant difference in grey matter density between ABGA-positive and ABGA-negative groups in caudate nuclei, putamina, thalami and frontal lobes. These results suggest that ABGA synthesis is not related to structural changes in grey and white matter (detectable with these methods) within frontostriatal circuits.