Population dynamics of genetically diverse Plasmodium falciparum lineages: community-based prospective study in rural Amazonia

Monitoring Plasmodium falciparum and Plasmodium vivax using microsatellite markers indicates limited changes in population structure after substantial transmission decline in Papua New Guinea

Monitoring the genetic structure of pathogen populations may be an economical and sensitive approach to quantify the impact of control on transmission dynamics, highlighting the need for a better understanding of changes in population genetic parameters as transmission declines. Here we describe the first population genetic analysis of two major human malaria parasites, Plasmodium falciparum (Pf) and Plasmodium vivax (Pv), following nationwide distribution of long-lasting insecticide-treated nets (LLINs) in Papua New Guinea (PNG). Parasite isolates from pre- (2005-2006) and post-LLIN (2010-2014) were genotyped using microsatellite markers. Despite parasite prevalence declining substantially (East Sepik Province: Pf = 54.9%-8.5%, Pv = 35.7%-5.6%, Madang Province: Pf = 38.0%-9.0%, Pv: 31.8%-19.7%), genetically diverse and intermixing parasite populations remained. Pf diversity declined modestly post-LLIN relative to pre-LLIN (East Sepik: R_s  = 7.1-6.4, H_E  = 0.77-0.71; Madang: R_s  = 8.2-6.1, H_E  = 0.79-0.71). Unexpectedly, population structure present in pre-LLIN populations was lost post-LLIN, suggesting that more frequent human movement between provinces may have contributed to higher gene flow. Pv prevalence initially declined but increased again in one province, yet diversity remained high throughout the study period (East Sepik: R_s  = 11.4-9.3, H_E  = 0.83-0.80; Madang: R_s  = 12.2-14.5, H_E  = 0.85-0.88). Although genetic differentiation values increased between provinces over time, no significant population structure was observed at any time point. For both species, a decline in multiple infections and increasing clonal transmission and significant multilocus linkage disequilibrium post-LLIN were positive indicators of impact on the parasite population using microsatellite markers. These parameters may be useful adjuncts to traditional epidemiological tools in the early stages of transmission reduction.

Estimating Relatedness Between Malaria Parasites

Understanding the relatedness of individuals within or between populations is a common goal in biology. Increasingly, relatedness features in genetic epidemiology studies of pathogens. These studies are relatively new compared to those in humans and other organisms, but are important for designing interventions and understanding pathogen transmission. Only recently have researchers begun to routinely apply relatedness to apicomplexan eukaryotic malaria parasites, and to date have used a range of different approaches on an ad hoc basis. Therefore, it remains unclear how to compare different studies and which measures to use. Here, we systematically compare measures based on identity-by-state (IBS) and identity-by-descent (IBD) using a globally diverse data set of malaria parasites, Plasmodium falciparum and P. vivax, and provide marker requirements for estimates based on IBD. We formally show that the informativeness of polyallelic markers for relatedness inference is maximized when alleles are equifrequent. Estimates based on IBS are sensitive to allele frequencies, which vary across populations and by experimental design. For portability across studies, we thus recommend estimates based on IBD. To generate estimates with errors below an arbitrary threshold of 0.1, we recommend ∼100 polyallelic or 200 biallelic markers. Marker requirements are immediately applicable to haploid malaria parasites and other haploid eukaryotes. C.I.s facilitate comparison when different marker sets are used. This is the first attempt to provide rigorous analysis of the reliability of, and requirements for, relatedness inference in malaria genetic epidemiology. We hope it will provide a basis for statistically informed prospective study design and surveillance strategies.

Quantifying connectivity between local Plasmodium falciparum malaria parasite populations using identity by descent

With the rapidly increasing abundance and accessibility of genomic data, there is a growing interest in using population genetic approaches to characterize fine-scale dispersal of organisms, providing insight into biological processes across a broad range of fields including ecology, evolution and epidemiology. For sexually recombining haploid organisms such as the human malaria parasite P. falciparum, however, there have been no systematic assessments of the type of data and methods required to resolve fine scale connectivity. This analytical gap hinders the use of genomics for understanding local transmission patterns, a crucial goal for policy makers charged with eliminating this important human pathogen. Here we use data collected from four clinics with a catchment area spanning approximately 120 km of the Thai-Myanmar border to compare the ability of divergence (FST) and relatedness based on identity by descent (IBD) to resolve spatial connectivity between malaria parasites collected from proximal clinics. We found no relationship between inter-clinic distance and FST, likely due to sampling of highly related parasites within clinics, but a significant decline in IBD-based relatedness with increasing inter-clinic distance. This association was contingent upon the data set type and size. We estimated that approximately 147 single-infection whole genome sequenced parasite samples or 222 single-infection parasite samples genotyped at 93 single nucleotide polymorphisms (SNPs) were sufficient to recover a robust spatial trend estimate at this scale. In summary, surveillance efforts cannot rely on classical measures of genetic divergence to measure P. falciparum transmission on a local scale. Given adequate sampling, IBD-based relatedness provides a useful alternative, and robust trends can be obtained from parasite samples genotyped at approximately 100 SNPs.

Genetic diversity of Plasmodium falciparum isolates from uncomplicated malaria cases in Ghana over a decade

Background

Genotyping malaria parasites to assess their diversity in different geographic settings have become necessary for the selection of antigenic epitopes for vaccine development and for antimalarial drug efficacy or resistance investigations. This study describes the genetic diversity of Plasmodium falciparum isolates from uncomplicated malaria cases over a ten year period (2003–2013) in Ghana using the polymorphic antigenic marker, merozoite surface protein 2 (msp2).

Methods

Archived filter paper blood blots from children aged nine years and below with uncomplicated malaria collected from nine sites in Ghana were typed for the presence of the markers. A total of 880 samples were genotyped for msp2 for the two major allelic families, FC27 and 3D7, using nested polymerase chain reaction (PCR). The allele frequencies and the multiplicity of infection were determined for the nine sites for five time points over a period of ten years, 2003–2004, 2005–2006, 2007–2008, 2010 and 2012–2013 malaria transmission seasons.

Results

The number of different alleles detected for the msp2 gene by resolving PCR products on agarose gels was 14. Both of the major allelic families, 3D7 and FC27 were common in all population samples. The highest multiplicity of infection (MOI) was observed in isolates from Begoro (forest zone, rural site): 3.31 for the time point 2007–2008. A significant variation was observed among the sites in the MOIs detected per infection (Fisher's exact test, P < 0.001) for the 2007 isolates and also at each of the three sites with data for three different years, Hohoe, P = 0.03; Navrongo, P < 0.001; Cape Coast, P < 0.001. Overall, there was no significant difference between the MOIs of the three ecological zones over the years (P = 0.37) and between the time points when data from all sites were pooled (P = 0.40).

Conclusions

The diversity and variation between isolates detected using the msp2 gene in Ghanaian isolates were observed to be profound; however, there was homogeneity throughout the three ecological zones studied. This is indicative of gene flow between the parasite populations across the country probably due to human population movements (HPM).

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1692-1) contains supplementary material, which is available to authorized users.

Comparison of allele frequencies of Plasmodium falciparum merozoite antigens in malaria infections sampled in different years in a Kenyan population

Background

Plasmodium falciparum merozoite antigens elicit antibody responses in malaria-endemic populations, some of which are clinically protective, which is one of the reasons why merozoite antigens are the focus of malaria vaccine development efforts. Polymorphisms in several merozoite antigen-encoding genes are thought to arise as a result of selection by the human immune system.

Methods

The allele frequency distribution of 15 merozoite antigens over a two-year period, 2007 and 2008, was examined in parasites obtained from children with uncomplicated malaria. In the same population, allele frequency changes pre- and post-anti-malarial treatment were also examined. Any gene which showed a significant shift in allele frequencies was also assessed longitudinally in asymptomatic and complicated malaria infections.

Results

Fluctuating allele frequencies were identified in codons 147 and 148 of reticulocyte-binding homologue (Rh) 5, with a shift from HD to YH haplotypes over the two-year period in uncomplicated malaria infections. However, in both the asymptomatic and complicated malaria infections YH was the dominant and stable haplotype over the two-year and ten-year periods, respectively. A logistic regression analysis of all three malaria infection populations between 2007 and 2009 revealed, that the chance of being infected with the HD haplotype decreased with time from 2007 to 2009 and increased in the uncomplicated and asymptomatic infections.

Conclusion

Rh5 codons 147 and 148 showed heterogeneity at both an individual and population level and may be under some degree of immune selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1304-8) contains supplementary material, which is available to authorized users.

Plasmodium vivax populations are more genetically diverse and less structured than sympatric Plasmodium falciparum populations

Introduction

The human malaria parasite, Plasmodium vivax, is proving more difficult to control and eliminate than Plasmodium falciparum in areas of co-transmission. Comparisons of the genetic structure of sympatric parasite populations may provide insight into the mechanisms underlying the resilience of P. vivax and can help guide malaria control programs.

Methodology/Principle findings

P. vivax isolates representing the parasite populations of four areas on the north coast of Papua New Guinea (PNG) were genotyped using microsatellite markers and compared with previously published microsatellite data from sympatric P. falciparum isolates. The genetic diversity of P. vivax (H_e = 0.83–0.85) was higher than that of P. falciparum (H_e = 0.64–0.77) in all four populations. Moderate levels of genetic differentiation were found between P. falciparum populations, even over relatively short distances (less than 50 km), with 21–28% private alleles and clear geospatial genetic clustering. Conversely, very low population differentiation was found between P. vivax catchments, with less than 5% private alleles and no genetic clustering observed. In addition, the effective population size of P. vivax (30353; 13043–69142) was larger than that of P. falciparum (18871; 8109–42986).

Conclusions/Significance

Despite comparably high prevalence, P. vivax had higher diversity and a panmictic population structure compared to sympatric P. falciparum populations, which were fragmented into subpopulations. The results suggest that in comparison to P. falciparum, P. vivax has had a long-term large effective population size, consistent with more intense and stable transmission, and limited impact of past control and elimination efforts. This underlines suggestions that more intensive and sustained interventions will be needed to control and eventually eliminate P. vivax. This research clearly demonstrates how population genetic analyses can reveal deeper insight into transmission patterns than traditional surveillance methods.

The neglected human malaria parasite Plasmodium vivax is responsible for a large proportion of the global malaria burden. Efforts to control malaria have revealed that P. vivax is more resilient than the other major human malaria parasite, Plasmodium falciparum. This study utilised population genetics to compare patterns of P. vivax and P. falciparum transmission in Papua New Guinea, a region where infection rates of the two species are similar. The results demonstrated that P. vivax populations are more genetically diverse than those of P. falciparum suggestive of a parasite population that is more resilient to environmental challenges, undergoing higher levels of interbreeding locally and between distant parasite populations. Unique characteristics of P. vivax such as relapse, which allows different strains from past infections to produce subsequent infections, may provide more opportunities for the exchange and dissemination of genetic material. The contrasting patterns observed for the two species may be the result of a differential impact of past elimination attempts and indicate that more rigorous interventions will be needed in efforts to control and eventually eliminate P. vivax.

Population structure of a microparasite infecting Daphnia: spatio-temporal dynamics

Background

Detailed knowledge of spatial and temporal variation in the genetic population structure of hosts and parasites is required for understanding of host − parasite coevolution. As hot-spots of contemporary coevolution in natural systems are difficult to detect and long-term studies are restricted to few systems, additional population genetic data from various host − parasite systems may provide important insights into the topic. This is particularly true for parasites, as these players have been under-investigated so far due to the lower availability of suitable molecular markers. Here, we traced genetic variation (based on sequence variants in the internal transcribed spacer region, ITS) among seven geographically isolated populations of the ichthyosporean Caullerya mesnili, a common microparasite of the cladoceran Daphnia (here, the D. longispina hybrid complex). At three sites, we also studied parasite genetic variation over time (three to four sampling points) and tested for associations between parasite genotypes and host species.

Results

Parasite (and host) populations were significantly structured across space, indicating limited dispersal. Moreover, the frequency of parasite genotypes varied significantly over time, suggesting rapid evolutionary change in Caullerya. However, the distribution of parasite genotypes was similar across different host species, which might in turn have important consequences for parasite epidemiology.

Conclusions

The approach proposed here can be applied to track spatial and temporal changes in the population structure of other microparasite species for which sequence variation in the ITS or other highly variable genome regions has been documented but other types of polymorphic markers are lacking. Screening of parasite sequence variants allows for reliable detection of cross-species infections and, using advanced sequencing techniques in the near future, for detailed studies of parasite evolution in natural host − parasite systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0247-3) contains supplementary material, which is available to authorized users.

Genetic structure of Plasmodium falciparum populations across the Honduras-Nicaragua border

Background

The Caribbean coast of Central America remains an area of malaria transmission caused by Plasmodium falciparum despite the fact that morbidity has been reduced in recent years. Parasite populations in that region show interesting characteristics such as chloroquine susceptibility and low mortality rates. Genetic structure and diversity of P. falciparum populations in the Honduras-Nicaragua border were analysed in this study.

Methods

Seven neutral microsatellite loci were analysed in 110 P. falciparum isolates from endemic areas of Honduras (n = 77) and Nicaragua (n = 33), mostly from the border region called the Moskitia. Several analyses concerning the genetic diversity, linkage disequilibrium, population structure, molecular variance, and haplotype clustering were conducted.

Results

There was a low level of genetic diversity in P. falciparum populations from Honduras and Nicaragua. Expected heterozigosity (H_e) results were similarly low for both populations. A moderate differentiation was revealed by the F_ST index between both populations, and two putative clusters were defined through a structure analysis. The main cluster grouped most of samples from Honduras and Nicaragua, while the second cluster was smaller and included all the samples from the Siuna community in Nicaragua. This result could partially explain the stronger linkage disequilibrium (LD) in the parasite population from that country. These findings are congruent with the decreasing rates of malaria endemicity in Central America.

Higher microsatellite diversity in Plasmodium vivax than in sympatric Plasmodium falciparum populations in Pursat, Western Cambodia

Previous microsatellite analyses of sympatric populations of Plasmodium vivax and Plasmodium falciparum in Brazil revealed higher diversity in the former species. However, it remains unclear whether regional species-specific differences in prevalence and transmission levels might account for these findings. Here, we examine sympatric populations of P. vivax (n=87) and P. falciparum (n=164) parasites from Pursat province, Western Cambodia, where both species are similarly prevalent. Using 10 genome-wide microsatellites for P. falciparum and 13 for P. vivax, we found that the P. vivax population was more diverse than the sympatric P. falciparum population (average virtual heterozygosity [HE], 0.87 vs. 0.66, P=0.003), with more multiple-clone infections (89.6% vs. 47.6%) and larger mean number of alleles per marker (16.2 vs. 11.1, P=0.07). Both populations showed significant multi-locus linkage disequilibrium suggestive of a predominantly clonal mode of parasite reproduction. The higher microsatellite diversity found in P. vivax isolates, compared to sympatric P. falciparum isolates, does not necessarily result from local differences in transmission level and may reflect differences in population history between species or increased mutation rates in P. vivax.

Long term persistence of clonal malaria parasite Plasmodium falciparum lineages in the Colombian Pacific region

Background

Resistance to chloroquine and antifolate drugs has evolved independently in South America, suggesting that genotype - phenotype studies aimed at understanding the genetic basis of resistance to these and other drugs should be conducted in this continent. This research was conducted to better understand the population structure of Colombian Plasmodium falciparum in preparation for such studies.

Results

A set of 384 SNPs were genotyped in blood spot DNA samples from 447 P. falciparum infected subjects collected over a ten year period from four provinces of the Colombian Pacific coast to evaluate clonality, population structure and linkage disequilibrium (LD). Most infections (81%) contained a single predominant clone. These clustered into 136 multilocus genotypes (MLGs), with 32% of MLGs recovered from multiple (2 – 28) independent subjects. We observed extremely low genotypic richness (R = 0.42) and long persistence of MLGs through time (median = 537 days, range = 1 – 2,997 days). There was a high probability (>5%) of sampling parasites from the same MLG in different subjects within 28 days, suggesting caution is needed when using genotyping methods to assess treatment success in clinical drug trials. Panmixia was rejected as four well differentiated subpopulations (F_ST = 0.084 - 0.279) were identified. These occurred sympatrically but varied in frequency within the four provinces. Linkage disequilibrium (LD) decayed more rapidly (r² = 0.17 for markers <10 kb apart) than observed previously in South American samples.

Conclusions

We conclude that Colombian populations have several advantages for association studies, because multiple clone infections are uncommon and LD decays over the scale of one or a few genes. However, the extensive population structure and low genotype richness will need to be accounted for when designing and analyzing association studies.

Microsatellite analysis of malaria parasites

Microsatellites have been increasingly used to investigate the population structure of malaria parasites, to map genetic loci contributing to phenotypes such as drug resistance and virulence in laboratory crosses and genome-wide association studies and to distinguish between treatment failures and new infections in clinical trials. Here, we provide optimized protocols for genotyping highly polymorphic microsatellites sampled from across the genomes of the human malaria parasites Plasmodium falciparum and P. vivax that have been extensively used in research laboratories worldwide.

Adaptive radiation within the vaccine target tetraspanin-23 across nine Schistosoma species from Africa

High levels of polymorphism in DNA sequences of tetraspanin-23 (TSP-23) were revealed within and between nine different species of Schistosoma from Africa including Schistosoma mansoni, Schistosoma rodhaini, Schistosoma margrebowiei, Schistosoma mattheei, Schistosoma intercalatum, Schistosoma haematobium, Schistosoma guineensis, Schistosoma curassoni and Schistosoma bovis. The greatest levels of diversity coincided with evidence of positive selection (d(N)/d(S)>1) within regions that code for extracellular loops of TSP-23 believed to interact with the host immune system. Kolaskar and Tongaonkar antigenicity predictions of protein sequences were compared across species and high levels of variation in antigenicity were also identified with each species which possessed their own unique antigenic profile. Phylogenetic analysis of TSP-23 proteins suggested evidence of convergent evolution in antigenic lineages as no true inter-species phylogenetic relationships were seen. This could be indicative of host-specific evolution of antigens in different species of schistosomes, a factor that should be considered carefully when developing vaccine targets.

Microsatellite loci over a thirty-three year period for a malaria parasite (Plasmodium mexicanum): bottleneck in effective population size and effect on allele frequencies

Changes in population allele frequencies may be driven by several forces, including selection and drift, and are revealed only by sampling over many generations. Such studies, however, are rare for protist parasites. Microsatellite allele frequencies for 4 loci were followed in a population of Plasmodium mexicanum, a malaria parasite of lizards in California USA at 1 site from 1978 to 2010. Rapid turnover of the lizards indicates the parasite was studied for a minimum of 33 transmission cycles and possibly twice that number. Sample sizes ranged from 841 to 956 scored parasite clones per locus. DNA was extracted from frozen dried blood and blood removed from stained blood smears from the earliest years, and a verification study demonstrated DNA from the blood smears provided valid genetic data. Parasite prevalence and effective population size (Ne) dropped after 2000, remaining lower for the next decade. For 2 loci, allele frequencies appeared stable for the first 2 decades of the study, but changed more rapidly after the decline in prevalence. Allele frequencies changed more gradually for the other 2 loci. Genetic drift could account for changes in allele frequencies, especially after the drop in prevalence and Ne, but the force of selection could also have driven the observed patterns.

Single-nucleotide polymorphism and copy number variation of the multidrug resistance-1 locus of Plasmodium vivax: local and global patterns

Emerging resistance to chloroquine (CQ) poses a major challenge for Plasmodium vivax malaria control, and nucleotide substitutions and copy number variation in the P. vivax multidrug resistance 1 (pvmdr-1) locus, which encodes a digestive vacuole membrane transporter, may modulate this phenotype. We describe patterns of genetic variation in pvmdr-1 alleles from Acre and Amazonas in northwestern Brazil, and compare then with those reported in other malaria-endemic regions. The pvmdr-1 mutation Y976F, which is associated with CQ resistance in Southeast Asia and Oceania, remains rare in northwestern Brazil (1.8%) and its prevalence mirrors that of CQ resistance worldwide. Gene amplification of pvmdr-1, which is associated with mefloquine resistance but increased susceptibility to CQ, remains relatively rare in northwestern Brazil (0.9%) and globally (< 4%), but became common (> 10%) in Tak Province, Thailand, possibly because of drug-mediated selection. The global database we have assembled provides a baseline for further studies of genetic variation in pvmdr-1 and drug resistance in P. vivax malaria.

Amazonian malaria: asymptomatic human reservoirs, diagnostic challenges, environmentally driven changes in mosquito vector populations, and the mandate for sustainable control strategies

Across the Americas and the Caribbean, nearly 561,000 slide-confirmed malaria infections were reported officially in 2008. The nine Amazonian countries accounted for 89% of these infections; Brazil and Peru alone contributed 56% and 7% of them, respectively. Local populations of the relatively neglected parasite Plasmodium vivax, which currently accounts for 77% of the regional malaria burden, are extremely diverse genetically and geographically structured. At a time when malaria elimination is placed on the public health agenda of several endemic countries, it remains unclear why malaria proved so difficult to control in areas of relatively low levels of transmission such as the Amazon Basin. We hypothesize that asymptomatic parasite carriage and massive environmental changes that affect vector abundance and behavior are major contributors to malaria transmission in epidemiologically diverse areas across the Amazon Basin. Here we review available data supporting this hypothesis and discuss their implications for current and future malaria intervention policies in the region. Given that locally generated scientific evidence is urgently required to support malaria control interventions in Amazonia, we briefly describe the aims of our current field-oriented malaria research in rural villages and gold-mining enclaves in Peru and a recently opened agricultural settlement in Brazil.

Quantification of Plasmodium falciparum malaria from complex infections in the Peruvian Amazon using quantitative PCR of the merozoite surface protein 1, block 2 (PfMSP1-B2): in vitro dynamics reveal density-dependent interactions

The majority of Plasmodium falciparum field isolates are defined as complex infections because they contain multiple genetically distinct clones. Studying interactions between clones in complex infections in vivo and in vitro could elucidate important phenomena in malaria infection, transmission and treatment. Using quantitative PCR (qPCR) of the P. falciparum merozoite surface protein 1, block 2 (PfMSP1-B2), we provide a sensitive and efficient genotyping method. This is important for epidemiological studies because it makes it possible to study genotype-specific growth dynamics. We compared 3 PfMSP1-B2 genotyping methods by analysing 79 field isolates from the Peruvian Amazon. In vivo observations from other studies using these techniques led to the hypothesis that clones within complex infections interact. By co-culturing clones with different PfMSP1-B2 genotypes, and measuring parasitaemia using qPCR, we found that suppression of clonal expansion was a factor of the collective density of all clones present in a culture. PfMSP1-B2 qPCR enabled us to find in vitro evidence for parasite-parasite interactions and could facilitate future investigations of growth trends in naturally occurring complex infections.

Plasmodium vivax populations revisited: mitochondrial genomes of temperate strains in Asia suggest ancient population expansion

BACKGROUND

Plasmodium vivax is the most widely distributed human malaria parasite outside of Africa, and its range extends well into the temperate zones. Previous studies provided evidence for vivax population differentiation, but temperate vivax parasites were not well represented in these analyses. Here we address this deficit by using complete mitochondrial (mt) genome sequences to elucidate the broad genetic diversity and population structure of P. vivax from temperate regions in East and Southeast Asia.

RESULTS

From the complete mtDNA sequences of 99 clinical samples collected in China, Myanmar and Korea, a total of 30 different haplotypes were identified from 26 polymorphic sites. Significant differentiation between different East and Southeast Asian parasite populations was observed except for the comparison between populations from Korea and southern China. Haplotype patterns and structure diversity analysis showed coexistence of two different groups in East Asia, which were genetically related to the Southeast Asian population and Myanmar population, respectively. The demographic history of P. vivax, examined using neutrality tests and mismatch distribution analyses, revealed population expansion events across the entire P. vivax range and the Myanmar population. Bayesian skyline analysis further supported the occurrence of ancient P. vivax population expansion.

CONCLUSIONS

This study provided further resolution of the population structure and evolution of P. vivax, especially in temperate/warm-temperate endemic areas of Asia. The results revealed divergence of the P. vivax populations in temperate regions of China and Korea from other populations. Multiple analyses confirmed ancient population expansion of this parasite. The extensive genetic diversity of the P. vivax populations is consistent with phenotypic plasticity of the parasites, which has implications for malaria control.

Molecular markers and genetic diversity of Plasmodium vivax

Enhanced understanding of the transmission dynamics and population genetics for Plasmodium vivax is crucial in predicting the emergence and spread of novel parasite phenotypes with major public health implications, such as new relapsing patterns, drug resistance and increased virulence. Suitable molecular markers are required for these population genetic studies. Here, we focus on two groups of molecular markers that are commonly used to analyse natural populations of P. vivax. We use markers under selective pressure, for instance, antigen-coding polymorphic genes, and markers that are not under strong natural selection, such as most minisatellite and microsatellite loci. First, we review data obtained using genes encoding for P. vivax antigens: circumsporozoite protein, merozoite surface proteins 1 and 3α, apical membrane antigen 1 and Duffy binding antigen. We next address neutral or nearly neutral molecular markers, especially microsatellite loci, providing a complete list of markers that have already been used in P. vivax populations studies. We also analyse the microsatellite loci identified in the P. vivax genome project. Finally, we discuss some practical uses for P. vivax genotyping, for example, detecting multiple-clone infections and tracking the geographic origin of isolates.

The evolution of drug resistance and the curious orthodoxy of aggressive chemotherapy

The evolution of drug-resistant pathogens is a major challenge for 21st century medicine. Drug use practices vigorously advocated as resistance management tools by professional bodies, public health agencies, and medical schools represent some of humankind's largest attempts to manage evolution. It is our contention that these practices have poor theoretical and empirical justification for a broad spectrum of diseases. For instance, rapid elimination of pathogens can reduce the probability that de novo resistance mutations occur. This idea often motivates the medical orthodoxy that patients should complete drug courses even when they no longer feel sick. Yet "radical pathogen cure" maximizes the evolutionary advantage of any resistant pathogens that are present. It could promote the very evolution it is intended to retard. The guiding principle should be to impose no more selection than is absolutely necessary. We illustrate these arguments in the context of malaria; they likely apply to a wide range of infections as well as cancer and public health insecticides. Intuition is unreliable even in simple evolutionary contexts; in a social milieu where in-host competition can radically alter the fitness costs and benefits of resistance, expert opinion will be insufficient. An evidence-based approach to resistance management is required.

Interrupting malaria transmission: quantifying the impact of interventions in regions of low to moderate transmission

Malaria has been eliminated from over 40 countries with an additional 39 currently planning for, or committed to, elimination. Information on the likely impact of available interventions, and the required time, is urgently needed to help plan resource allocation. Mathematical modelling has been used to investigate the impact of various interventions; the strength of the conclusions is boosted when several models with differing formulation produce similar data. Here we predict by using an individual-based stochastic simulation model of seasonal Plasmodium falciparum transmission that transmission can be interrupted and parasite reintroductions controlled in villages of 1,000 individuals where the entomological inoculation rate is <7 infectious bites per person per year using chemotherapy and bed net strategies. Above this transmission intensity bed nets and symptomatic treatment alone were not sufficient to interrupt transmission and control the importation of malaria for at least 150 days. Our model results suggest that 1) stochastic events impact the likelihood of successfully interrupting transmission with large variability in the times required, 2) the relative reduction in morbidity caused by the interventions were age-group specific, changing over time, and 3) the post-intervention changes in morbidity were larger than the corresponding impact on transmission. These results generally agree with the conclusions from previously published models. However the model also predicted changes in parasite population structure as a result of improved treatment of symptomatic individuals; the survival probability of introduced parasites reduced leading to an increase in the prevalence of sub-patent infections in semi-immune individuals. This novel finding requires further investigation in the field because, if confirmed, such a change would have a negative impact on attempts to eliminate the disease from areas of moderate transmission.

Plasmodium falciparum genetic diversity maintained and amplified over 5 years of a low transmission endemic in the Peruvian Amazon

Plasmodium falciparum entered into the Peruvian Amazon in 1994, sparking an epidemic between 1995 and 1998. Since 2000, there has been sustained low P. falciparum transmission. The Malaria Immunology and Genetics in the Amazon project has longitudinally followed members of the community of Zungarococha (N = 1,945, 4 villages) with active household and health center-based visits each year since 2003. We examined parasite population structure and traced the parasite genetic diversity temporally and spatially. We genotyped infections over 5 years (2003–2007) using 14 microsatellite (MS) markers scattered across ten different chromosomes. Despite low transmission, there was considerable genetic diversity, which we compared with other geographic regions. We detected 182 different haplotypes from 302 parasites in 217 infections. Structure v2.2 identified five clusters (subpopulations) of phylogenetically related clones. To consider genetic diversity on a more detailed level, we defined haplotype families (hapfams) by grouping haplotypes with three or less loci differences. We identified 34 different hapfams identified. The F_st statistic and heterozygosity analysis showed the five clusters were maintained in each village throughout this time. A minimum spanning network (MSN), stratified by the year of detection, showed that haplotypes within hapfams had allele differences and haplotypes within a cluster definition were more separated in the later years (2006–2007). We modeled hapfam detection and loss, accounting for sample size and stochastic fluctuations in frequencies overtime. Principle component analysis of genetic variation revealed patterns of genetic structure with time rather than village. The population structure, genetic diversity, appearance/disappearance of the different haplotypes from 2003 to 2007 provides a genome-wide “real-time” perspective of P. falciparum parasites in a low transmission region.

Single-nucleotide polymorphism, linkage disequilibrium and geographic structure in the malaria parasite Plasmodium vivax: prospects for genome-wide association studies

Background

The ideal malaria parasite populations for initial mapping of genomic regions contributing to phenotypes such as drug resistance and virulence, through genome-wide association studies, are those with high genetic diversity, allowing for numerous informative markers, and rare meiotic recombination, allowing for strong linkage disequilibrium (LD) between markers and phenotype-determining loci. However, levels of genetic diversity and LD in field populations of the major human malaria parasite P. vivax remain little characterized.

Results

We examined single-nucleotide polymorphisms (SNPs) and LD patterns across a 100-kb chromosome segment of P. vivax in 238 field isolates from areas of low to moderate malaria endemicity in South America and Asia, where LD tends to be more extensive than in holoendemic populations, and in two monkey-adapted strains (Salvador-I, from El Salvador, and Belem, from Brazil). We found varying levels of SNP diversity and LD across populations, with the highest diversity and strongest LD in the area of lowest malaria transmission. We found several clusters of contiguous markers with rare meiotic recombination and characterized a relatively conserved haplotype structure among populations, suggesting the existence of recombination hotspots in the genome region analyzed. Both silent and nonsynonymous SNPs revealed substantial between-population differentiation, which accounted for ~40% of the overall genetic diversity observed. Although parasites clustered according to their continental origin, we found evidence for substructure within the Brazilian population of P. vivax. We also explored between-population differentiation patterns revealed by loci putatively affected by natural selection and found marked geographic variation in frequencies of nucleotide substitutions at the pvmdr-1 locus, putatively associated with drug resistance.

Conclusion

These findings support the feasibility of genome-wide association studies in carefully selected populations of P. vivax, using relatively low densities of markers, but underscore the risk of false positives caused by population structure at both local and regional levels.

See commentary: http://www.biomedcentral.com/1741-7007/8/90

Multilocus genotyping reveals high heterogeneity and strong local population structure of the Plasmodium vivax population in the Peruvian Amazon

Background

Peru is one of the Latin American countries with the highest malaria burden, mainly due to Plasmodium vivax infections. However, little is known about P. vivax transmission dynamics in the Peruvian Amazon, where most malaria cases occur. The genetic diversity and population structure of P. vivax isolates collected in different communities around Iquitos city, the capital of the Peruvian Amazon, was determined.

Methods

Plasmodium vivax population structure was determined by multilocus genotyping with 16 microsatellites on 159 P. vivax infected blood samples (mono-infections) collected in four sites around Iquitos city. The population characteristics were assessed only in samples with monoclonal infections (n = 94), and the genetic diversity was determined by calculating the expected heterozygosity and allelic richness. Both linkage disequilibrium and the genetic differentiation (θ) were estimated.

Results

The proportion of polyclonal infections varied substantially by site (11% - 70%), with the expected heterozygosity ranging between 0.44 and 0.69; no haplotypes were shared between the different populations. Linkage disequilibrium was present in all populations (I_A^S 0.14 - 0.61) but was higher in those with fewer polyclonal infections, suggesting inbreeding and a clonal population structure. Strong population differentiation (θ = 0.45) was found and the Bayesian inference cluster analysis identified six clusters based on distinctive allele frequencies.

Conclusion

The P. vivax populations circulating in the Peruvian Amazon basin are genetically diverse, strongly differentiated and they have a low effective recombination rate. These results are in line with the low and clustered pattern of malaria transmission observed in the region around Iquitos city.