Search

Riverine floodplains exhibit high floral and faunal diversity as a consequence of their biophysical complexity. Extension of such niche partitioning processes to microbial communities is far less resolved or supported. Here, we evaluated the responses of aquatic biofilms diversity to environmental gradients across ten riverine floodplains with differing degrees of flow alteration and habitat diversity to assess whether complex floodplains support biofilm communities with greater biodiversity and species interactions. No significant evidence was found to support a central role for habitat diversity in promoting microbial diversity across 116 samples derived from 62 aquatic habitats, as neither α (H’: 2.8–4.1) nor β (Sørensen: 0.3–0.39) diversity were positively related to floodplain complexity across the ten floodplains. In contrast, our results documented the sensitivity of biofilm communities to regional templates manifested as gradients of carbon, nitrogen, and phosphorous availability. Large-scale conditions reflecting nitrogen limitation increased the relative abundance of N-fixing cyanobacteria (up to 0.34 as fraction of total reads), constrained the total number of interactions among bacterial taxa, and reinforced negative over positive interactions, generating unique microbial communities and networks that reflect large-scale species sorting in response to regional geochemical gradients.

The sustainability of malaria control in Africa is threatened by the rise of insecticide resistance in Anopheles mosquitoes, which transmit the disease(1). To gain a deeper understanding of how mosquito populations are evolving, here we sequenced the genomes of 765 specimens of Anopheles gambiae and Anopheles coluzzii sampled from 15 locations across Africa, and identified over 50 million single nucleotide polymorphisms within the accessible genome. These data revealed complex population structure and patterns of gene flow, with evidence of ancient expansions, recent bottlenecks, and local variation in effective population size. Strong signals of recent selection were observed in insecticide-resistance genes, with several sweeps spreading over large geographical distances and between species. The design of new tools for mosquito control using gene-drive systems will need to take account of high levels of genetic diversity in natural mosquito populations.

Estimating the effective population size and effective number of breeders per year (Nb ) can facilitate early detection of population declines. We used computer simulations to quantify bias and precision of the one-sample LDNe estimator of Nb in age-structured populations using a range of published species life history types, sample sizes, and DNA markers. Nb estimates were biased by ~5%-10% when using SNPs or microsatellites in species ranging from fishes to mosquitoes, frogs, and seaweed. The bias (high or low) was similar for different life history types within a species suggesting that life history variation in populations will not influence Nb estimation. Precision was higher for 100 SNPs (H ≈ 0.30) than for 15 microsatellites (H ≈ 0.70). Confidence intervals (CIs) were occasionally too narrow, and biased high when Nb was small (Nb

Background: High altitude is characterised by hypoxia, arid climate, high solar radiation and other environmental stressors. Andeans have resided at high elevations for several millennia and have developed an unique array of physiological adaptations. However, the genetic changes that led to these adaptations remain largely elusive. Objective: The aim of this study is to identify genetic variants that harbour signatures of recent positive selection and may facilitate physiological adaptations to hypobaric hypoxia. Methods: We conducted whole genome sequencing and lung function tests in 19 Argentinean highlanders (>3500 m) comparing them to 16 lowlanders. We developed a new statistical procedure using a combination of population branch statistics and number of segregating sites by length to detect beneficial alleles that arose since the settlement of the Andes and are currently present in 15-50% of the population. Results: We identified two missense variants as significant targets of selection. One variant is located within the GPR126 gene and has been previously associated with the forced expiratory volume/forced vital capacity ratio. The derived allele of GPR126 is also associated with lung function in our sample of highlanders (p Conclusions: These variants may contribute to the physiological adaptations to hypobaric hypoxia, possibly by altering lung function and preventing excessive erythrocytosis and hypoxia induced pulmonary hypertension in Andean populations.

We present Bioconda (https://bioconda.github.io), a distribution of bioinformatics software for the lightweight, multiplatform and language-agnostic package manager Conda. Currently, Bioconda offers a collection of over 3000 software packages, which is continuously maintained, updated, and extended by a growing global community of more than 200 contributors. Bioconda improves analysis reproducibility by allowing users to define isolated environments with defined software versions, all of which are easily installed and managed without administrative privileges.

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.

Arctic populations live in an environment characterized by extreme cold and the absence of plant foods for much of the year and are likely to have undergone genetic adaptations to these environmental conditions in the time they have been living there. Genome-wide selection scans based on genotype data from native Siberians have previously highlighted a 3 Mb chromosome 11 region containing 79 protein-coding genes as the strongest candidates for positive selection in Northeast Siberians. However, it was not possible to determine which of the genes might be driving the selection signal. Here, using whole-genome high-coverage sequence data, we identified the most likely causative variant as a nonsynonymous G>A transition (rs80356779; c.1436C>T [p.Pro479Leu] on the reverse strand) in CPT1A, a key regulator of mitochondrial long-chain fatty-acid oxidation. Remarkably, the derived allele is associated with hypoketotic hypoglycemia and high infant mortality yet occurs at high frequency in Canadian and Greenland Inuits and was also found at 68% frequency in our Northeast Siberian sample. We provide evidence of one of the strongest selective sweeps reported in humans; this sweep has driven this variant to high frequency in circum-Arctic populations within the last 6–23 ka despite associated deleterious consequences, possibly as a result of the selective advantage it originally provided to either a high-fat diet or a cold environment.

In iteroparous species, it is easier to estimate Nb (effective number of breeders in one reproductive cycle) than Ne (effective population size per generation). Nb can be used as a proxy for Ne and also can provide crucial insights into eco-evolutionary processes that occur during reproduction. We used analytical and numerical methods to evaluate effects of intermittent breeding and litter/clutch size on inbreeding Nb and Ne . Fixed or random litter sizes ≥ 3 have little effect on either effective-size parameter; however, in species (e.g., many large mammals) in which females can produce only one offspring per cycle, female Nb = ∞ and overall Nb = 4Nb (male) . Intermittent breeding reduces the pool of female breeders, which reduces both female and overall Nb ; reductions are larger in high-fecundity species with high juvenile mortality and increase when multiple reproductive cycles are skipped. Simulated data for six model species showed that both intermittent breeding and litter-size constraints increase Ne , but only slightly. We show how to quantitatively account for these effects, which are important to consider when (1) using Nb to estimate Ne , or (2) drawing inferences about male reproductive success based on estimates of female Nb .

The aim of this study is to identify genetic variants that harbour signatures of recent positive selection and may facilitate physiological adaptations to hypobaric hypoxia. To achieve this, we conducted whole genome sequencing and lung function tests in 19 Argentinean highlanders (>3500 m) comparing them to 16 Native American lowlanders. We developed a new statistical procedure using a combination of population branch statistics (PBS) and number of segregating sites by length (nSL) to detect beneficial alleles that arose since the settlement of the Andes and are currently present in 15–50% of the population. We identified two missense variants as significant targets of selection. One of these variants, located within the GPR126 gene, has been previously associated with the forced expiratory volume/forced vital capacity ratio. The other novel missense variant mapped to the EPAS1 gene encoding the hypoxia inducible factor 2α. EPAS1 is known to be the major selection candidate gene in Tibetans. The derived allele of GPR126 is associated with lung function in our sample of highlanders (p

The sustainability of malaria control in Africa is threatened by rising levels of insecticide resistance, and new tools to prevent malaria transmission are urgently needed. To gain a better understanding of the mosquito populations that transmit malaria, we sequenced the genomes of 765 wild specimens ofAnopheles gambiaeandAnopheles coluzziisampled from 15 locations across Africa. The data reveal high levels of genetic diversity, with over 50 million single nucleotide polymorphisms across the 230 Mbp genome. We observe complex patterns of population structure and marked variations in local population size, some of which may be due at least in part to malaria control interventions. Insecticide resistance genes show strong signatures of recent selection associated with multiple independent mutations spreading over large geographical distances and between species. The genetic variability of natural populations substantially reduces the target space for novel gene-drive strategies for mosquito control. This large dataset provides a foundation for tracking the emergence and spread of insecticide resistance and developing new vector control tools.

Humans have more copies of amylase genes than other primates. It is still poorly understood, however, when the copy number expansion occurred and whether its spread was enhanced by selection. Here we assess amylase copy numbers in a global sample of 480 high coverage genomes and find that regions flanking the amylase locus show notable depression of genetic diversity both in African and non-African populations. Analysis of genetic variation in these regions supports the model of an early selective sweep in the human lineage after the split of humans from Neanderthals which led to the fixation of multiple copies of AMY1 in place of a single copy. We find evidence of multiple secondary losses of copy number with the highest frequency (52%) of a deletion of AMY2A and associated low copy number of AMY1 in Northeast Siberian populations whose diet has been low in starch content.

International audience; The history of human settlement in Southeast Asia has been complex and involved several distinct dispersal events. Here, we report the analyses of 1825 individuals from Southeast Asia including new genome-wide genotype data for 146 individuals from three Mainland Southeast Asian (Burmese, Malay and Vietnamese) and four Island Southeast Asian (Dusun, Filipino, Kankanaey and Murut) populations. While confirming the presence of previously recognised major ancestry components in the Southeast Asian population structure, we highlight the Kankanaey Igorots from the highlands of the Philippine Mountain Province as likely the closest living representatives of the source population that may have given rise to the Austronesian expansion. This conclusion rests on independent evidence from various analyses of autosomal data and uniparental markers. Given the extensive presence of trade goods, cultural and linguistic evidence of Indian influence in Southeast Asia starting from 2.5 kya, we also detect traces of a South Asian signature in different populations in the region dating to the last couple of thousand years.

High-coverage whole-genome sequence studies have so far focused on a limited number1 of geographically restricted populations2, 3, 4, 5, or been targeted at specific diseases, such as cancer6. Nevertheless, the availability of high- resolution genomic data has led to the development of new methodologies for inferring population history7, 8, 9 and refuelled the debate on the mutation rate in humans10. Here we present the Estonian Biocentre Human Genome Diversity Panel (EGDP), a dataset of 483 high-coverage human genomes from 148 populations worldwide, including 379 new genomes from 125 populations, which we group into diversity and selection sets. We analyse this dataset to refine estimates of continent-wide patterns of heterozygosity, long- and short-distance gene flow, archaic admixture, and changes in effective population size through time as well as for signals of positive or balancing selection. We find a genetic signature in present-day Papuans that suggests that at least 2% of their genome originates from an early and largely extinct expansion of anatomically modern humans (AMHs) out of Africa. Together with evidence from the western Asian fossil record11, and admixture between AMHs and Neanderthals predating the main Eurasian expansion12, our results contribute to the mounting evidence for the presence of AMHs out of Africa earlier than 75, 000 years ago.

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.

Early detection of population declines is essential to prevent extinctions and to ensure sustainable harvest. We evaluated the performance of two Ne estimators to detect population declines: the two-sample temporal method and a one-sample method based on linkage disequilibrium (LD). We used simulated data representing a wide range of population sizes, sample sizes and number of loci. Both methods usually detect a population decline only one generation after it occurs if Ne drops to less than approximately 100, and 40 microsatellite loci and 50 individuals are sampled. However, the LD method often out performed the temporal method by allowing earlier detection of less severe population declines (Ne approximately 200). Power for early detection increased more rapidly with the number of individuals sampled than with the number of loci genotyped, primarily for the LD method. The number of samples available is therefore an important criterion when choosing between the LD and temporal methods. We provide guidelines regarding design of studies targeted at monitoring for population declines. We also report that 40 single nucleotide polymorphism (SNP) markers give slightly lower precision than 10 microsatellite markers. Our results suggest that conservation management and monitoring strategies can reliably use genetic based methods for early detection of population declines.

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.

Retrieving and organizing data from complete genomes is a time-consuming task, even more so if the interest lies only in part of the genome (for nongenomic analysis). Furthermore, when comparing several genomes or genes, data retrieval has to be repeated multiple times. We present baca, a software for retrieving, organizing and visualizing multiple mitochondrial genomes. baca takes a GenBank query, retrieves all related genomes and generates multiple fasta files organized both by genomes and genes. A web-based user interface and an interactive graphical map of all genomes with all genes are also provided. The program is available from http://cibio.up.pt/software/baca.

Arsenic is a carcinogen associated with skin lesions and cardiovascular diseases. The Colla population from the Puna region in Northwest Argentinean is exposed to levels of arsenic in drinking water exceeding the recommended maximum by a factor of 20. Yet, they thrive in this challenging environment since thousands of years and therefore we hypothesize strong selection signatures in genes involved in arsenic metabolism. We analyzed genome-wide genotype data for 730,000 loci in 25 Collas, considering 24 individuals of the neighbouring Calchaquíes and 24 Wichí from the Gran Chaco region in the Argentine province of Salta as control groups. We identified a strong signal of positive selection in the main arsenic methyltransferase AS3MT gene, which has been previously associated with lower concentrations of the most toxic product of arsenic metabolism monomethylarsonic acid. This study confirms recent studies reporting selection signals in the AS3MT gene albeit using different samples, tests and control populations.

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 (F ST>0.2). Outside Central Asia F ST 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., Author Summary Plasmodium vivax is the predominant malaria parasite in Latin America, Asia and the South Pacific. Different factors are expected to shape diversity and population structure across continents, e.g. transmission intensity which is much lower in South America as compared to Southeast-Asia and the South Pacific, or geographical isolation of P. vivax populations in the South Pacific. We have compiled data from 841 isolates from South and Central America, Africa, Central Asia, Southeast-Asia and the South Pacific typed with a panel of 11 microsatellite markers. Diversity was highest in Southeast-Asia, where transmission is intermediate-high and migration of infected hosts is high, and lowest in South America and Central Asia where malaria transmission is low and focal. Reducing the panel of microsatellites showed that 2–6 markers are sufficient for genotyping for most drug trials and epidemiological studies, as these markers can identify >90% of all haplotypes. Parasites clustered according to continental origin, with high population differentiation between South American and Central Asian populations and the other populations, and lowest differences between Southeast-Asia and the South Pacific. Current attempts to reduce malaria transmission might change this pattern, but only after transmission is reduced for an extended period of time.

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: see text] was a linear (r(2) = 0.98) function of the harmonic mean of Ne and Nb. We provide a quantitative bias correction for raw [Formula: see text] 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).

This paper discusses the development of a debugging tool for parallel programs showing how the requirements posed by high-level tools for parallel program development have influenced the design of the debugging system since its early stages of development. We concentrate our attention upon the interfacing of the debugger with other tools of a parallel software engineering environment, namely a graphical programming language and a testing and debugging tool. This is illustrated with the results of our experimentation with the design and implementation of DDBG, a debugger for the PVM environment.

International audience; Linguistic and cultural evidence suggest that Madagascar was the final point of two major dispersals of Austronesian- and Bantu-speaking populations. Today, the Mikea are described as the last-known Malagasy population reported to be still practicing a hunter-gatherer lifestyle. It is unclear, however, whether the Mikea descend from a remnant population that existed before the arrival of Austronesian and Bantu agriculturalists or whether it is only their lifestyle that separates them from the other contemporary populations of South Madagascar. To address these questions we have performed a genome-wide analysis of >700,000 SNP markers on 21 Mikea, 24 Vezo, and 24 Temoro individuals, together with 50 individuals from Bajo and Lebbo populations from Indonesia. Our analyses of these data in the context of data available from other Southeast Asian and African populations reveal that all three Malagasy populations are derived from the same admixture event involving Austronesian and Bantu sources. In contrast to the fact that most of the vocabulary of the Malagasy speakers is derived from the Barito group of the Austronesian language family, we observe that only one-third of their genetic ancestry is related to the populations of the Java-Kalimantan-Sulawesi area. Because no additional ancestry components distinctive for the Mikea were found, it is likely that they have adopted their hunter-gatherer way of life through cultural reversion, and selection signals suggest a genetic adaptation to their new lifestyle.

During their migrations out of Africa, humans successfully colonised and adapted to a wide range of habitats, including extreme high altitude environments, where reduced atmospheric oxygen (hypoxia) imposes a number of physiological challenges. This study evaluates genetic and phenotypic variation in the Colla population living in the Argentinean Andes above 3500 m and compares it to the nearby lowland Wichí group in an attempt to pinpoint evolutionary mechanisms underlying adaptation to high altitude hypoxia. We genotyped 730,525 SNPs in 25 individuals from each population. In genome-wide scans of extended haplotype homozygosity Collas showed the strongest signal around VEGFB, which plays an essential role in the ischemic heart, and ELTD1, another gene crucial for heart development and prevention of cardiac hypertrophy. Moreover, pathway enrichment analysis showed an overrepresentation of pathways associated with cardiac morphology. Taken together, these findings suggest that Colla highlanders may have evolved a toolkit of adaptative mechanisms resulting in cardiac reinforcement, most likely to counteract the adverse effects of the permanently increased haematocrit and associated shear forces that characterise the Andean response to hypoxia. Regulation of cerebral vascular flow also appears to be part of the adaptive response in Collas. These findings are not only relevant to understand the evolution of hypoxia protection in high altitude populations but may also suggest new avenues for medical research into conditions where hypoxia constitutes a detrimental factor

In iteroparous species, it is easier to estimate Nb (effective number of breeders in one reproductive cycle) than Ne (effective population size per generation). Nb can be used as a proxy for Ne and also can provide crucial insights into eco-evolutionary processes that occur during reproduction. We used analytical and numerical methods to evaluate effects of intermittent breeding and litter/clutch size on inbreeding Nb and Ne . Fixed or random litter sizes ≥ 3 have little effect on either effective-size parameter; however, in species (e.g., many large mammals) in which females can produce only one offspring per cycle, female Nb = ∞ and overall Nb = 4Nb (male) . Intermittent breeding reduces the pool of female breeders, which reduces both female and overall Nb ; reductions are larger in high-fecundity species with high juvenile mortality and increase when multiple reproductive cycles are skipped. Simulated data for six model species showed that both intermittent breeding and litter-size constraints increase Ne , but only slightly. We show how to quantitatively account for these effects, which are important to consider when (1) using Nb to estimate Ne , or (2) drawing inferences about male reproductive success based on estimates of female Nb .

INTRODUCTION: The importance of Plasmodium vivax in malaria elimination is increasingly being recognized yet little is known about its population size and population genetic structure in the South Pacific an area that is the focus of intensified malaria control. METHODS: We have genotyped 13 microsatellite markers in 295 P. vivax isolates from four geographically distinct sites in Papua New Guinea (PNG) and one site from Solomon Islands representing different transmission intensities. RESULTS: Diversity was very high with expected heterozygosity values ranging from 0.62 to 0.98 for the different markers. Effective population size was high (12'872 to 19'533 per site). In PNG population structuring was limited with moderate levels of genetic differentiation. F ST values (adjusted for high diversity of markers) were 0.14 0.15. Slightly higher levels were observed between PNG populations and Solomon Islands (F ST?=?0.16). CONCLUSIONS: Low levels of population structure despite geographical barriers to transmission are in sharp contrast to results from regions of low P. vivax endemicity. Prior to intensification of malaria control programs in the study area parasite diversity and effective population size remained high.

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. Contact: tra@popgen.eu Supplementary information: Supplementary data is available at Bioinformatics online.

Plasmodium falciparum malaria is subject to artificial selection from antimalarial drugs that select for drug-resistant parasites. We describe and apply a flexible new approach to investigate how epistasis, inbreeding, selection heterogeneity and multiple simultaneous drug deployments interact to influence the spread of drug-resistant malaria. This framework recognizes that different human ‘environments’ within which treatment may occur (such as semi- and non-immune humans taking full or partial drug courses) influence the genetic interactions between parasite loci involved in resistance. Our model provides an explanation for how the rate of spread varies according to different malaria transmission intensities, why resistance might stabilize at intermediate frequencies and also identifies several factors that influence the decline of resistance after a drug is removed. Results suggest that studies based on clinical outcomes might overestimate the spread of resistant parasites, especially in high-transmission areas. We show that when transmission decreases, prevalence might decrease without a corresponding change in frequency of resistance and that this relationship is heavily influenced by the extent of linkage disequilibrium between loci. This has important consequences on the interpretation of data from areas where control is being successful and suggests that reducing transmission might have less impact on the spread of resistance than previously expected.

Early detection of population declines is essential to prevent extinctions and to ensure sustainable harvest. We evaluated the performance of two N e estimators to detect population declines: the two-sample temporal method and a one-sample method based on linkage disequilibrium (LD). We used simulated data representing a wide range of population sizes, sample sizes and number of loci. Both methods usually detect a population decline only one generation after it occurs if N e drops to less than approximately 100, and 40 microsatellite loci and 50 individuals are sampled. However, the LD method often out performed the temporal method by allowing earlier detection of less severe population declines (N e approximately 200). Power for early detection increased more rapidly with the number of individuals sampled than with the number of loci genotyped, primarily for the LD method. The number of samples available is therefore an important criterion when choosing between the LD and temporal methods. We provide guidelines regarding design of studies targeted at monitoring for population declines. We also report that 40 single nucleotide polymorphism (SNP) markers give slightly lower precision than 10 microsatellite markers. Our results suggest that conservation management and monitoring strategies can reliably use genetic based methods for early detection of population declines.

Summary: MODELER4SIMCOAL2 (M4S2) is an extensible graphical tool to model linked loci and population demographies. M4S2 is easy to use, allowing for the modeling of complicated scenarios, making coalescent simulation modeling accessible to biologists with limited computer skills. The software includes an extension system allowing for new models to be created, published and downloaded from the Internet.Availability: M4S2 is available from http://popgen.eu/soft/m4s2 under a GPL license. The web site also contains guides, screen shots and tutorials.Contact: tra@fc.up.pt

Huijben et al. [1] recently suggested that subcurative antimalarial drug treatment of clinical cases might slow the spread of drug resistance through increased competition between drug resistant and sensitive infections. This was further discussed in letters to Trends in Parasitology. Goncalves and Paul [2] raised important questions about the potential use of subcurative treatment: What is the epidemiological and clinical impact of leaving a patient with circulating parasites? How is the correct drug sub-clearance level determined? Hastings [3] added several ethical and operational arguments against such a proposal namely that patients treated with subcurative drug levels could have repeated episodes of recurrent malaria and/or succumb to a secondary infection after being weakened by the primary malaria episode.

Human pygmy populations inhabit different regions of the world, from Africa to Melanesia. In Asia, short-statured populations are often referred to as "negritos." Their short stature has been interpreted as a consequence of thermoregulatory, nutritional, and/or locomotory adaptations to life in tropical forests. A more recent hypothesis proposes that their stature is the outcome of a life history trade-off in high-mortality environments, where early reproduction is favored and, consequently, early sexual maturation and early growth cessation have coevolved. Some serological evidence of deficiencies in the growth hormone/insulin-like growth factor axis have been previously associated with pygmies' short stature. Using genome-wide single-nucleotide polymorphism genotype data, we first tested whether different negrito groups living in the Philippines and Papua New Guinea are closely related and then investigated genomic signals of recent positive selection in African, Asian, and Papuan pygmy populations. We found that negritos in the Philippines and Papua New Guinea are genetically more similar to their nonpygmy neighbors than to one another and have experienced positive selection at different genes. These results indicate that geographically distant pygmy groups are likely to have evolved their short stature independently. We also found that selection on common height variants is unlikely to explain their short stature and that different genes associated with growth, thyroid function, and sexual development are under selection in different pygmy groups.