Plasmodium vivax Diversity and Population Structure across Four Continents

Plasmodium vivax is the geographically most widespread human malaria parasite. To analyze patterns of microsatellite diversity and population structure across countries of different transmission intensity, genotyping data from 11 microsatellite markers was either generated or compiled from 841 isolates from four continents collected in 1999–2008. Diversity was highest in South-East Asia (mean allelic richness 10.0–12.8), intermediate in the South Pacific (8.1–9.9) Madagascar and Sudan (7.9–8.4), and lowest in South America and Central Asia (5.5–7.2). A reduced panel of only 3 markers was sufficient to identify approx. 90% of all haplotypes in South Pacific, African and SE-Asian populations, but only 60–80% in Latin American populations, suggesting that typing of 2–6 markers, depending on the level of endemicity, is sufficient for epidemiological studies. Clustering analysis showed distinct clusters in Peru and Brazil, but little sub-structuring was observed within Africa, SE-Asia or the South Pacific. Isolates from Uzbekistan were exceptional, as a near-clonal parasite population was observed that was clearly separated from all other populations (FST>0.2). Outside Central Asia FST values were highest (0.11–0.16) between South American and all other populations, and lowest (0.04–0.07) between populations from South-East Asia and the South Pacific. These comparisons between P. vivax populations from four continents indicated that not only transmission intensity, but also geographical isolation affect diversity and population structure. However, the high effective population size results in slow changes of these parameters. This persistency must be taken into account when assessing the impact of control programs on the genetic structure of parasite populations

Massive introgression drives species radiation at the range limit of Anopheles gambiae

Impacts of introgressive hybridisation may range from genomic erosion and species collapse to rapid adaptation and speciation but opportunities to study these dynamics are rare. We investigated the extent, causes and consequences of a hybrid zone between Anopheles coluzzii and Anopheles gambiae in Guinea-Bissau, where high hybridisation rates appear to be stable at least since the 1990s. Anopheles gambiae was genetically partitioned into inland and coastal subpopulations, separated by a central region dominated by A. coluzzii. Surprisingly, whole genome sequencing revealed that the coastal region harbours a hybrid form characterised by an A. gambiae-like sex chromosome and massive introgression of A. coluzzii autosomal alleles. Local selection on chromosomal inversions may play a role in this process, suggesting potential for spatiotemporal stability of the coastal hybrid form and providing resilience against introgression of medically-important loci and traits, found to be more prevalent in inland A. gambiae

Local selection in the presence of high levels of gene flow: Evidence of heterogeneous insecticide selection pressure across Ugandan Culex quinquefasciatus populations

Background: Culex quinquefasciatus collected in Uganda, where no vector control interventions directly targeting this species have been conducted, was used as a model to determine if it is possible to detect heterogeneities in selection pressure driven by insecticide application targeting other insect species. Methodology/Principal findings: Population genetic structure was assessed through microsatellite analysis, and the impact of insecticide pressure by genotyping two target-site mutations, Vgsc-1014F of the voltage-gated sodium channel target of pyrethroid and DDT insecticides, and Ace1-119S of the acetylcholinesterase gene, target of carbamate and organophosphate insecticides. No significant differences in genetic diversity were observed among populations by microsatellite markers with HE ranging from 0.597 to 0.612 and low, but significant, genetic differentiation among populations (FST = 0.019, P = 0.001). By contrast, the insecticide-resistance markers display heterogeneous allelic distributions with significant differences detected between Central Ugandan (urban) populations relative to Eastern and Southwestern (rural) populations. In the central region, a frequency of 62% for Vgsc-1014F, and 32% for the Ace1-119S resistant allele were observed. Conversely, in both Eastern and Southwestern regions the Vgsc-1014F alleles were close to fixation, whilst Ace1-119S allele frequency was 12% (although frequencies may be underestimated due to copy number variation at both loci). Conclusions/Significance: Taken together, the microsatellite and both insecticide resistance target-site markers provide evidence that in the face of intense gene flow among populations, disjunction in resistance frequencies arise due to intense local selection pressures despite an absence of insecticidal control interventions targeting Culex

ogaraK: a population genetics simulator for malaria.

Motivation: The evolution of resistance in Plasmodium falciparum malaria against most available treatments is a major global health threat. Population genetics approaches are commonly used to model the spread of drug resistance. Due to uncommon features in malaria biology, existing forward-time population genetics simulators cannot suitably model Plasmodium falciparum malaria. Results: Here we present ogaraK, a population genetics simulator for modelling the spread of drug-resistant malaria. OgaraK is designed to make malaria simulation computationally tractable as it models infections, not individual parasites. OgaraK is also able to model the life cycle of the parasite which includes both haploid and diploid phases and sexual and asexual reproduction. We also allow for the simulation of different inbreeding levels, an important difference between high and low transmission areas and a fundamental factor influencing the outcome of strategies to control or eliminate malaria. Availability: OgaraK is available as free software (GPL) from the address http://popgen.eu/soft/ogaraK

Mcheza: a workbench to detect selection using dominant markers.

Motivation: Dominant markers (DArTs and AFLPs) are commonly used for genetic analysis in the fields of evolutionary genetics, ecology and conservation of genetic resources. The recent prominence of these markers has coincided with renewed interest in detecting the effects of local selection and adaptation at the level of the genome. Results: We present Mcheza, an application for detecting loci under selection based on a well-evaluated FST-outlier method. The application allows robust estimates to be made of model parameters (e.g. genome-wide average, neutral FST), provides data import and export functions, iterative contour smoothing and generation of graphics in an easy to use graphical user interface with a computation engine that supports multicore processors for enhanced performance. Mcheza also provides functionality to mitigate common analytical errors when scanning for loci under selection. Availability: Mcheza is freely available under GPL version 3 from http://popgen.eu/soft/mcheza

Effects of Overlapping Generations on Linkage Disequilibrium Estimates of Effective Population Size

Use of single-sample genetic methods to estimate effective population size has skyrocketed in recent years. Although the underlying models assume discrete generations, they are widely applied to age-structured species. We simulated genetic data for 21 iteroparous animal and plant species to evaluate two untested hypotheses regarding performance of the single-sample method based on linkage disequilibrium (LD): (1) estimates based on single-cohort samples reflect the effective number of breeders in one reproductive cycle (Nb), and (2) mixed-age samples reflect the effective size per generation (Ne). We calculated true Ne and Nb, using the model species’ vital rates, and verified these with individual-based simulations. We show that single-cohort samples should be equally influenced by Nb and Ne and confirm this with simulated results: Formula was a linear (r2 = 0.98) function of the harmonic mean of Ne and Nb. We provide a quantitative bias correction for raw Formula based on the ratio Nb/Ne, which can be estimated from two or three simple life history traits. Bias-adjusted estimates were within 5% of true Nb for all 21 study species and proved robust when challenged with new data. Mixed-age adult samples produced downwardly biased estimates in all species, which we attribute to a two-locus Wahlund effect (mixture LD) caused by combining parents from different cohorts in a single sample. Results from this study will facilitate interpretation of rapidly accumulating genetic estimates in terms of both Ne (which influences long-term evolutionary processes) and Nb (which is more important for understanding eco-evolutionary dynamics and mating systems)

Genetic and phenotypic differentiation of an Andean intermediate altitude population

Highland populations living permanently under hypobaric hypoxia have been subject of extensive research because of the relevance of their physiological adaptations for the understanding of human health and disease. In this context, what is considered high altitude is a matter of interpretation and while the adaptive processes at high altitude (above 3000 m) are well documented, the effects of moderate altitude (below 3000 m) on the phenotype are less well established. In this study, we compare physiological and anthropometric characteristics as well as genetic variations in two Andean populations: the Calchaquíes (2300 m) and neighboring Collas (3500 m). We compare their phenotype and genotype to the sea‐level Wichí population. We measured physiological (heart rate, oxygen saturation, respiration rate, and lung function) as well as anthropometric traits (height, sitting height, weight, forearm, and tibia length). We conducted genome‐wide genotyping on a subset of the sample (n = 74) and performed various scans for positive selection. At the phenotypic level (n = 179), increased lung capacity stood out in both Andean groups, whereas a growth reduction in distal limbs was only observed at high altitude. At the genome level, Calchaquíes revealed strong signals around PRKG1, suggesting that the nitric oxide pathway may be a target of selection. PRKG1 was highlighted by one of four selection tests among the top five genes using the population branch statistic. Selection tests results of Collas were reported previously. Overall, our study shows that some phenotypic and genetic differentiation occurs at intermediate altitude in response to moderate lifelong selection pressures

Conserving genomic variability in large mammals: Effect of population fluctuations and variance in male reproductive success on variability in Yellowstone bison

Loss of genetic variation through genetic drift can reduce population viability. However, relatively little is known about loss of variation caused by the combination of fluctuating population size and variance in reproductive success in age structured populations. We built an individual-based computer simulation model to examine how actual culling and hunting strategies influence the effective population size (Ne) and allelic diversity in Yellowstone bison over 200 years (∼28 generations). The Ne for simulated populations ranged from 746 in stable populations of size 2000 up to 1165 in fluctuating populations whose census size fluctuates between 3000 and 3500 individuals. Simulations suggested that ∼93% of allelic diversity, for loci with five alleles will be maintained over 200 years if the population census size remains well above ∼2000 bison (and if variance in male reproductive success is high). However for loci with 20 alleles, only 83% of allelic diversity will be maintained over 200 years. Removal of only juveniles (calves and yearlings) resulted in longer generation intervals which led to higher maintenance of allelic diversity (96%) after 200 years compared to the culling of adults (94%) when the mean census size was 3250 (for loci with five alleles). These simulations suggest that fluctuations in population census size do not necessarily accelerate the loss of genetic variation, at least for the relatively large census size and growing populations such as in Yellowstone bison. They also suggest that the conservation of high allelic diversity (>95%) at loci with many alleles (e.g., ⩾5) will require maintenance of a populations size greater than approximately 3250 and removal of mainly or only juveniles

Data from: Plasmodium vivax diversity and population structure across four continents

Plasmodium vivax is the geographically most widespread human malaria parasite. To analyze patterns of microsatellite diversity and population structure across countries of different transmission intensity, genotyping data from 11 microsatellite markers was either generated or compiled from 841 isolates from four continents collected in 1999–2008. Diversity was highest in South-East Asia (mean allelic richness 10.0–12.8), intermediate in the South Pacific (8.1–9.9) Madagascar and Sudan (7.9–8.4), and lowest in South America and Central Asia (5.5–7.2). A reduced panel of only 3 markers was sufficient to identify approx. 90% of all haplotypes in South Pacific, African and SE-Asian populations, but only 60–80% in Latin American populations, suggesting that typing of 2–6 markers, depending on the level of endemicity, is sufficient for epidemiological studies. Clustering analysis showed distinct clusters in Peru and Brazil, but little sub-structuring was observed within Africa, SE-Asia or the South Pacific. Isolates from Uzbekistan were exceptional, as a near-clonal parasite population was observed that was clearly separated from all other populations (FST>0.2). Outside Central Asia FST values were highest (0.11–0.16) between South American and all other populations, and lowest (0.04–0.07) between populations from South-East Asia and the South Pacific. These comparisons between P. vivax populations from four continents indicated that not only transmission intensity, but also geographical isolation affect diversity and population structure. However, the high effective population size results in slow changes of these parameters. This persistency must be taken into account when assessing the impact of control programs on the genetic structure of parasite populations