Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence

Background

The malaria parasite Plasmodium falciparum exhibits abundant genetic diversity, and this diversity is key to its success as a pathogen. Previous efforts to study genetic diversity in P. falciparum have begun to elucidate the demographic history of the species, as well as patterns of population structure and patterns of linkage disequilibrium within its genome. Such studies will be greatly enhanced by new genomic tools and recent large-scale efforts to map genomic variation. To that end, we have developed a high throughput single nucleotide polymorphism (SNP) genotyping platform for P. falciparum.

Results

Using an Affymetrix 3,000 SNP assay array, we found roughly half the assays (1,638) yielded high quality, 100% accurate genotyping calls for both major and minor SNP alleles. Genotype data from 76 global isolates confirm significant genetic differentiation among continental populations and varying levels of SNP diversity and linkage disequilibrium according to geographic location and local epidemiological factors. We further discovered that nonsynonymous and silent (synonymous or noncoding) SNPs differ with respect to within-population diversity, inter-population differentiation, and the degree to which allele frequencies are correlated between populations.

Conclusions

The distinct population profile of nonsynonymous variants indicates that natural selection has a significant influence on genomic diversity in P. falciparum, and that many of these changes may reflect functional variants deserving of follow-up study. Our analysis demonstrates the potential for new high-throughput genotyping technologies to enhance studies of population structure, natural selection, and ultimately enable genome-wide association studies in P. falciparum to find genes underlying key phenotypic traits.     

Large-scale Genotyping and Genetic Mapping in Plasmodium Parasites

The completion of many malaria parasite genomes provides great opportunities for genomewide characterization of gene expression and high-throughput genotyping. Substantial progress in malaria genomics and genotyping has been made recently, particularly the development of various microarray platforms for large-scale characterization of the Plasmodium falciparum genome. Microarray has been used for gene expression analysis, detection of single nucleotide polymorphism (SNP) and copy number variation (CNV), characterization of chromatin modifications, and other applications. Here we discuss some recent advances in genetic mapping and genomic studies of malaria parasites, focusing on the use of high-throughput arrays for the detection of SNP and CNV in the P. falciparum genome. Strategies for genetic mapping of malaria traits are also discussed.     

Prospective Identification of Malaria Parasite Genes under Balancing Selection

Background

Endemic human pathogens are subject to strong immune selection, and interrogation of pathogen genome variation for signatures of balancing selection can identify important target antigens. Several major antigen genes in the malaria parasite Plasmodium falciparum have shown such signatures in polymorphism-versus-divergence indices (comparing with the chimpanzee parasite P. reichenowi), and in allele frequency based indices.

Methodology/Principal Findings

To compare methods for prospective identification of genes under balancing selection, 26 additional genes known or predicted to encode surface-exposed proteins of the invasive blood stage merozoite were first sequenced from a panel of 14 independent P. falciparum cultured lines and P. reichenowi. Six genes at the positive extremes of one or both of the Hudson-Kreitman-Aguade (HKA) and McDonald-Kreitman (MK) indices were identified. Allele frequency based analysis was then performed on a Gambian P. falciparum population sample for these six genes and three others as controls. Tajima''s D (TjD) index was most highly positive for the msp3/6-like PF10_0348 (TjD = 1.96) as well as the positive control ama1 antigen gene (TjD = 1.22). Across the genes there was a strong correlation between population TjD values and the relative HKA indices (whether derived from the population or the panel of cultured laboratory isolates), but no correlation with the MK indices.

Conclusions/Significance

Although few individual parasite genes show significant evidence of balancing selection, analysis of population genomic and comparative sequence data with the HKA and TjD indices should discriminate those that do, and thereby identify likely targets of immunity.     

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).     

Historical Shifts in Brazilian P. falciparum Population Structure and Drug Resistance Alleles

Previous work suggests that Brazilian Plasmodium falciparum has limited genetic diversity and a history of bottlenecks, multiple reintroductions due to human migration, and clonal expansions. We hypothesized that Brazilian P. falciparum would exhibit clonal structure. We examined isolates collected across two decades from Amapá, Rondônia, and Pará state (n = 190). By examining more microsatellites markers on more chromosomes than previous studies, we hoped to define the extent of low diversity, linkage disequilibrium, bottlenecks, population structure, and parasite migration within Brazil. We used retrospective genotyping of samples from the 1980s and 1990s to explore the population genetics of SP resistant dhfr and dhps alleles. We tested an existing hypothesis that the triple mutant dhfr mutations 50R/51I/108N and 51I/108N/164L developed in southern Amazon from a single origin of common or similar parasites. We found that Brazilian P. falciparum had limited genetic diversity and isolation by distance was rejected, which suggests it underwent bottlenecks followed by migration between sites. Unlike Peru, there appeared to be gene flow across the Brazilian Amazon basin. We were unable to divide parasite populations by clonal lineages and pairwise FST were common. Most parasite diversity was found within sites in the Brazilian Amazon, according to AMOVA. Our results challenge the hypothesis that triple mutant alleles arose from a single lineage in the Southern Amazon. SP resistance, at both the double and triple mutant stages, developed twice and potentially in different regions of the Brazilian Amazon. We would have required samples from before the 1980s to describe how SP resistance spread across the basin or describe the complex internal migration of Brazilian parasites after the colonization efforts of past decades. The Brazilian Amazon basin may have sufficient internal migration for drug resistance reported in any particular region to rapidly spread to other parts of basin under similar drug pressure.     

New insights into the evolutionary history of Plasmodium falciparum from mitochondrial genome sequence analyses of Indian isolates

Estimating genetic diversity and inferring the evolutionary history of Plasmodium falciparum could be helpful in understanding origin and spread of virulent and drug‐resistant forms of the malaria pathogen and therefore contribute to malaria control programme. Genetic diversity of the whole mitochondrial (mt) genome of P. falciparum sampled across the major distribution ranges had been reported, but no Indian P. falciparum isolate had been analysed so far, even though India is highly endemic to P. falciparum malaria. We have sequenced the whole mt genome of 44 Indian field isolates and utilized published data set of 96 genome sequences to present global genetic diversity and to revisit the evolutionary history of P. falciparum. Indian P. falciparum presents high genetic diversity with several characteristics of ancestral populations and shares many of the genetic features with African and to some extent Papua New Guinean (PNG) isolates. Similar to African isolates, Indian P. falciparum populations have maintained high effective population size and undergone rapid expansion in the past with oldest time to the most recent common ancestor (TMRCA). Interestingly, one of the four single nucleotide polymorphisms (SNPs) that differentiates P. falciparum from P. falciparum‐like isolates (infecting non‐human primates in Africa) was found to be segregating in five Indian P. falciparum isolates. This SNP was in tight linkage with other two novel SNPs that were found exclusively in these five Indian isolates. The results on the mt genome sequence analyses of Indian isolates on the whole add to the current understanding on the evolutionary history of P. falciparum.     

Plasmodium falciparum Mating Patterns and Mosquito Infectivity of Natural Isolates of Gametocytes

Plasmodium falciparum infections in malaria endemic areas often harbor multiple clones of parasites. However, the transmission success of the different genotypes within the mosquito vector has remained elusive so far. The genetic diversity of malaria parasites was measured by using microsatellite markers in gametocyte isolates from 125 asymptomatic carriers. For a subset of 49 carriers, the dynamics of co-infecting genotypes was followed until their development within salivary glands. Also, individual oocysts from midguts infected with blood from 9 donors were genotyped to assess mating patterns. Multiplicity of infection (MOI) was high both in gametocyte isolates and sporozoite populations, reaching up to 10 genotypes. Gametocyte isolates with multiple genotypes gave rise to lower infection prevalence and intensity. Fluctuations of genotype number occurred during the development within the mosquito and sub-patent genotypes, not detected in gametocyte isolates, were identified in the vector salivary glands. The inbreeding coefficient Fis was positively correlated to the oocyst loads, suggesting that P. falciparum parasites use different reproductive strategies according to the genotypes present in the gametocyte isolate. The number of parasite clones within an infection affects the transmission success and the mosquito has an important role in maintaining P. falciparum genetic diversity. Our results emphasize the crucial importance of discriminating between the different genotypes within an infection when studying the A. gambiae natural resistance to P. falciparum, and the need to monitor parasite diversity in areas where malaria control interventions are implemented.     

Plasmodium falciparum mitochondrial genetic diversity exhibits isolation-by-distance patterns supporting a sub-Saharan African origin

The geographical distribution of single nucleotide polymorphism (SNP) in the mitochondrial genome of the human malaria parasite Plasmodium falciparum was investigated. We identified 88 SNPs in 516 isolates from seven parasite populations in Africa, Southeast Asia and Oceania. Analysis of the SNPs postulated a sub-Saharan African origin and recovered a strong negative correlation between within-population SNP diversity and geographic distance from the putative African origin over Southeast Asia and Oceania. These results are consistent with those previously obtained for nuclear genome-encoded housekeeping genes, indicating that the pattern of inheritance does not substantially affect the geographical distribution of SNPs.     

ama1 genes of sympatric Plasmodium vivax and P. falciparum from Venezuela differ significantly in genetic diversity and recombination frequency

Background

We present the first population genetic analysis of homologous loci from two sympatric human malaria parasite populations sharing the same human hosts, using full-length sequences of ama1 genes from Plasmodium vivax and P. falciparum collected in the Venezuelan Amazon.

Methodology/Principal Findings

Significant differences between the two species were found in genetic diversity at the ama1 locus, with 18 distinct haplotypes identified among the 73 Pvama1 sequences obtained, compared to 6 unique haplotypes from 30 Pfama1 sequences, giving overall diversity estimates of h = 0.9091, and h = 0.538 respectively. Levels of recombination were also found to differ between the species, with P. falciparum exhibiting very little recombination across the 1.77kb sequence. In contrast, analysis of patterns of nucleotide substitutions provided evidence that polymorphisms in the ama1 gene of both species are maintained by balancing selection, particularly in domain I. The two distinct population structures observed are unlikely to result from different selective forces acting upon the two species, which share both human and mosquito hosts in this setting. Rather, the highly structured P. falciparum population appears to be the result of a population bottleneck, while the much less structured P. vivax population is likely to be derived from an ancient pool of diversity, as reflected in a larger estimate of effective population size for this species. Greatly reduced mosquito transmission in 1997, due to low rainfall prior to the second survey, was associated with far fewer P. falciparum infections, but an increase in P. vivax infections, probably due to hypnozoite activation.

Conclusions/Significance

The relevance of these findings to putative competitive interactions between these two important human pathogen species is discussed. These results highlight the need for future control interventions to employ strategies targeting each of the parasite species present in endemic areas.     

Whole Genome Sequencing of Field Isolates Provides Robust Characterization of Genetic Diversity in Plasmodium vivax

Background

An estimated 2.85 billion people live at risk of Plasmodium vivax transmission. In endemic countries vivax malaria causes significant morbidity and its mortality is becoming more widely appreciated, drug-resistant strains are increasing in prevalence, and an increasing number of reports indicate that P. vivax is capable of breaking through the Duffy-negative barrier long considered to confer resistance to blood stage infection. Absence of robust in vitro propagation limits our understanding of fundamental aspects of the parasite''s biology, including the determinants of its dormant hypnozoite phase, its virulence and drug susceptibility, and the molecular mechanisms underlying red blood cell invasion.

Methodology/Principal Findings

Here, we report results from whole genome sequencing of five P. vivax isolates obtained from Malagasy and Cambodian patients, and of the monkey-adapted Belem strain. We obtained an average 70–400 X coverage of each genome, resulting in more than 93% of the Sal I reference sequence covered by 20 reads or more. Our study identifies more than 80,000 SNPs distributed throughout the genome which will allow designing association studies and population surveys. Analysis of the genome-wide genetic diversity in P. vivax also reveals considerable allele sharing among isolates from different continents. This observation could be consistent with a high level of gene flow among parasite strains distributed throughout the world.

Conclusions

Our study shows that it is feasible to perform whole genome sequencing of P. vivax field isolates and rigorously characterize the genetic diversity of this parasite. The catalogue of polymorphisms generated here will enable large-scale genotyping studies and contribute to a better understanding of P. vivax traits such as drug resistance or erythrocyte invasion, partially circumventing the lack of laboratory culture that has hampered vivax research for years.     

Cloning of Plasmodium falciparum by single-cell sorting

Malaria parasite cloning is traditionally carried out mainly by using the limiting dilution method, which is laborious, imprecise, and unable to distinguish multiply-infected RBCs. In this study, we used a parasite engineered to express green fluorescent protein (GFP) to evaluate a single-cell sorting method for rapidly cloning Plasmodium falciparum. By dividing a two-dimensional scattergram from a cell sorter into 17 gates, we determined the parameters for isolating singly-infected erythrocytes and sorted them into individual cultures. Pre-gating of the engineered parasites for GFP allowed the isolation of almost 100% GFP-positive clones. Compared with the limiting dilution method, the number of parasite clones obtained by single-cell sorting was much higher. Molecular analyses showed that parasite isolates obtained by single-cell sorting were highly homogenous. This highly efficient single-cell sorting method should prove very useful for cloning both P. falciparum laboratory populations from genetic manipulation experiments and clinical samples.     

The Plasmodium vivax Merozoite Surface Protein 3β Sequence Reveals Contrasting Parasite Populations in Southern and Northwestern Thailand

Background

Malaria control efforts have a significant impact on the epidemiology and parasite population dynamics. In countries aiming for malaria elimination, malaria transmission may be restricted to limited transmission hot spots, where parasite populations may be isolated from each other and experience different selection forces. Here we aim to examine the Plasmodium vivax population divergence in geographically isolated transmission zones in Thailand.

Methodology

We employed the P. vivax merozoite surface protein 3β (PvMSP3β) as a molecular marker for characterizing P. vivax populations based on the extensive diversity of this gene in Southeast Asian parasite populations. To examine two parasite populations with different transmission levels in Thailand, we obtained 45 P. vivax isolates from Tak Province, northwestern Thailand, where the annual parasite incidence (API) was more than 2%, and 28 isolates from Yala and Narathiwat Provinces, southern Thailand, where the API was less than 0.02%. We sequenced the PvMSP3β gene and examined its genetic diversity and molecular evolution between the parasite populations.

Principal Findings

Of 58 isolates containing single PvMSP3β alleles, 31 sequence types were identified. The overall haplotype diversity was 0.77±0.06 and nucleotide diversity 0.0877±0.0054. The northwestern vivax malaria population exhibited extensive haplotype diversity (HD) of PvMSP3β (HD = 1.0). In contrast, the southern parasite population displayed a single PvMSP3β allele (HD = 0), suggesting a clonal population expansion. This result revealed that the extent of allelic diversity in P. vivax populations in Thailand varies among endemic areas.

Conclusion

Malaria parasite populations in a given region may vary significantly in genetic diversity, which may be the result of control and influenced by the magnitude of malaria transmission intensity. This is an issue that should be taken into account for the implementation of P. vivax control measures such as drug policy and vaccine development.     

Evolution of Genetic Polymorphisms of Plasmodium falciparum Merozoite Surface Protein (PfMSP) in Thailand

Plasmodium falciparum malaria is a major public health problem in Thailand due to the emergence of multidrug resistance. The understanding of genetic diversity of malaria parasites is essential for developing effective drugs and vaccines. The genetic diversity of the merozoite surface protein-1 (PfMSP-1) and merozoite surface protein-2 (PfMSP-2) genes was investigated in a total of 145 P. falciparum isolates collected from Mae Sot District, Tak Province, Thailand during 3 different periods (1997-1999, 2005-2007, and 2009-2010). Analysis of genetic polymorphisms was performed to track the evolution of genetic change of P. falciparum using PCR. Both individual genes and their combination patterns showed marked genetic diversity during the 3 study periods. The results strongly support that P. falciparum isolates in Thailand are markedly diverse and patterns changed with time. These 2 polymorphic genes could be used as molecular markers to detect multiple clone infections and differentiate recrudescence from reinfection in P. falciparum isolates in Thailand.     

Effective Preparation of Plasmodium vivax Field Isolates for High-Throughput Whole Genome Sequencing

Whole genome sequencing (WGS) of Plasmodium vivax is problematic due to the reliance on clinical isolates which are generally low in parasitaemia and sample volume. Furthermore, clinical isolates contain a significant contaminating background of host DNA which confounds efforts to map short read sequence of the target P. vivax DNA. Here, we discuss a methodology to significantly improve the success of P. vivax WGS on natural (non-adapted) patient isolates. Using 37 patient isolates from Indonesia, Thailand, and travellers, we assessed the application of CF11-based white blood cell filtration alone and in combination with short term ex vivo schizont maturation. Although CF11 filtration reduced human DNA contamination in 8 Indonesian isolates tested, additional short-term culture increased the P. vivax DNA yield from a median of 0.15 to 6.2 ng µl⁻¹ packed red blood cells (pRBCs) (p = 0.001) and reduced the human DNA percentage from a median of 33.9% to 6.22% (p = 0.008). Furthermore, post-CF11 and culture samples from Thailand gave a median P. vivax DNA yield of 2.34 ng µl⁻¹ pRBCs, and 2.65% human DNA. In 22 P. vivax patient isolates prepared with the 2-step method, we demonstrate high depth (median 654X coverage) and breadth (≥89%) of coverage on the Illumina GAII and HiSeq platforms. In contrast to the A+T-rich P. falciparum genome, negligible bias was observed in coverage depth between coding and non-coding regions of the P. vivax genome. This uniform coverage will greatly facilitate the detection of SNPs and copy number variants across the genome, enabling unbiased exploration of the natural diversity in P. vivax populations.     

Neutral Polymorphisms in Putative Housekeeping Genes and Tandem Repeats Unravels the Population Genetics and Evolutionary History of Plasmodium vivax in India

The evolutionary history and age of Plasmodium vivax has been inferred as both recent and ancient by several studies, mainly using mitochondrial genome diversity. Here we address the age of P. vivax on the Indian subcontinent using selectively neutral housekeeping genes and tandem repeat loci. Analysis of ten housekeeping genes revealed a substantial number of SNPs (n = 75) from 100 P. vivax isolates collected from five geographical regions of India. Neutrality tests showed a majority of the housekeeping genes were selectively neutral, confirming the suitability of housekeeping genes for inferring the evolutionary history of P. vivax. In addition, a genetic differentiation test using housekeeping gene polymorphism data showed a lack of geographical structuring between the five regions of India. The coalescence analysis of the time to the most recent common ancestor estimate yielded an ancient TMRCA (232,228 to 303,030 years) and long-term population history (79,235 to 104,008) of extant P. vivax on the Indian subcontinent. Analysis of 18 tandem repeat loci polymorphisms showed substantial allelic diversity and heterozygosity per locus, and analysis of potential bottlenecks revealed the signature of a stable P. vivax population, further corroborating our ancient age estimates. For the first time we report a comparable evolutionary history of P. vivax inferred by nuclear genetic markers (putative housekeeping genes) to that inferred from mitochondrial genome diversity.     

Practical PCR genotyping protocols for Plasmodium vivax using Pvcs and Pvmsp1

Background

Plasmodium vivax is the second most prevalent malaria parasite affecting more than 75 million people each year, mostly in South America and Asia. In addition to major morbidity this parasite is associated with relapses and a reduction in birthweight. The emergence and spread of drug resistance in Plasmodium falciparum is a major factor in the resurgence of this parasite. P. vivax resistance to drugs has more recently emerged and monitoring the situation would be helped, as for P. falciparum, by molecular methods that can be used to characterize parasites in field studies and drug efficacy trials.

Methods

Practical PCR genotyping protocols based on polymorphic loci present in two P. vivax genetic markers, Pvcs and Pvmsp1, were developed. The methodology was evaluated using 100 P. vivax isolates collected in Thailand.

Results and Discussion

Analysis revealed that P. vivax populations in Thailand are highly diverse genetically, with mixed genotype infections found in 26 % of the samples (average multiplicity of infection = 1.29). A large number of distinguishable alleles were found for the two markers, 23 for Pvcs and 36 for Pvmsp1. These were generally randomly distributed amongst the isolates. A total of 68 distinct genotypes could be enumerated in the 74 isolates with a multiplicity of infection of 1.

Conclusion

These results indicate that the genotyping protocols presented can be useful in the assessment of in vivo drug efficacy clinical trials conducted in endemic areas and for epidemiological studies of P. vivax infections.     

The population genetics of Pseudomonas aeruginosa isolates from different patient populations exhibits high-level host specificity

Objective

To determine whether highly prevalent P. aeruginosa sequence types (ST) in Dutch cystic fibrosis (CF) patients are specifically linked to CF patients we investigated the population structure of P. aeruginosa from different clinical backgrounds. We first selected the optimal genotyping method by comparing pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST) and multilocus variable number tandem-repeat analysis (MLVA).

Methods

Selected P. aeruginosa isolates (n = 60) were genotyped with PFGE, MLST and MLVA to determine the diversity index (DI) and congruence (adjusted Rand and Wallace coefficients). Subsequently, isolates from patients admitted to two different ICUs (n = 205), from CF patients (n = 100) and from non-ICU, non-CF patients (n = 58, of which 19 were community acquired) were genotyped with MLVA to determine distribution of genotypes and genetic diversity.

Results

Congruence between the typing methods was >79% and DIs were similar and all >0.963. Based on costs, ease, speed and possibilities to compare results between labs an adapted MLVA scheme called MLVA9-Utrecht was selected as the preferred typing method. In 363 clinical isolates 252 different MLVA types (MTs) were identified, indicating a highly diverse population (DI  = 0.995; CI  = 0.993–0.997). DI levels were similarly high in the diverse clinical sources (all >0.981) and only eight genotypes were shared. MTs were highly specific (>80%) for the different patient populations, even for similar patient groups (ICU patients) in two distinct geographic regions, with only three of 142 ICU genotypes detected in both ICUs. The two major CF clones were unique to CF patients.

Conclusion

The population structure of P. aeruginosa isolates is highly diverse and population specific without evidence for a core lineage in which major CF, hospital or community clones co-cluster. The two genotypes highly prevalent among Dutch CF patients appeared unique to CF patients, suggesting specific adaptation of these clones to the CF lung.     

Validation of Genotyping-By-Sequencing Analysis in Populations of Tetraploid Alfalfa by 454 Sequencing

Genotyping-by-sequencing (GBS) is a relatively low-cost high throughput genotyping technology based on next generation sequencing and is applicable to orphan species with no reference genome. A combination of genome complexity reduction and multiplexing with DNA barcoding provides a simple and affordable way to resolve allelic variation between plant samples or populations. GBS was performed on ApeKI libraries using DNA from 48 genotypes each of two heterogeneous populations of tetraploid alfalfa (Medicago sativa spp. sativa): the synthetic cultivar Apica (ATF0) and a derived population (ATF5) obtained after five cycles of recurrent selection for superior tolerance to freezing (TF). Nearly 400 million reads were obtained from two lanes of an Illumina HiSeq 2000 sequencer and analyzed with the Universal Network-Enabled Analysis Kit (UNEAK) pipeline designed for species with no reference genome. Following the application of whole dataset-level filters, 11,694 single nucleotide polymorphism (SNP) loci were obtained. About 60% had a significant match on the Medicago truncatula syntenic genome. The accuracy of allelic ratios and genotype calls based on GBS data was directly assessed using 454 sequencing on a subset of SNP loci scored in eight plant samples. Sequencing depth in this study was not sufficient for accurate tetraploid allelic dosage, but reliable genotype calls based on diploid allelic dosage were obtained when using additional quality filtering. Principal Component Analysis of SNP loci in plant samples revealed that a small proportion (<5%) of the genetic variability assessed by GBS is able to differentiate ATF0 and ATF5. Our results confirm that analysis of GBS data using UNEAK is a reliable approach for genome-wide discovery of SNP loci in outcrossed polyploids.     

Differing Patterns of Selection and Geospatial Genetic Diversity within Two Leading Plasmodium vivax Candidate Vaccine Antigens

Although Plasmodium vivax is a leading cause of malaria around the world, only a handful of vivax antigens are being studied for vaccine development. Here, we investigated genetic signatures of selection and geospatial genetic diversity of two leading vivax vaccine antigens – Plasmodium vivax merozoite surface protein 1 (pvmsp-1) and Plasmodium vivax circumsporozoite protein (pvcsp). Using scalable next-generation sequencing, we deep-sequenced amplicons of the 42 kDa region of pvmsp-1 (n = 44) and the complete gene of pvcsp (n = 47) from Cambodian isolates. These sequences were then compared with global parasite populations obtained from GenBank. Using a combination of statistical and phylogenetic methods to assess for selection and population structure, we found strong evidence of balancing selection in the 42 kDa region of pvmsp-1, which varied significantly over the length of the gene, consistent with immune-mediated selection. In pvcsp, the highly variable central repeat region also showed patterns consistent with immune selection, which were lacking outside the repeat. The patterns of selection seen in both genes differed from their P. falciparum orthologs. In addition, we found that, similar to merozoite antigens from P. falciparum malaria, genetic diversity of pvmsp-1 sequences showed no geographic clustering, while the non-merozoite antigen, pvcsp, showed strong geographic clustering. These findings suggest that while immune selection may act on both vivax vaccine candidate antigens, the geographic distribution of genetic variability differs greatly between these two genes. The selective forces driving this diversification could lead to antigen escape and vaccine failure. Better understanding the geographic distribution of genetic variability in vaccine candidate antigens will be key to designing and implementing efficacious vaccines.     

Polymorphisms in K13 and Falcipain-2 Associated with Artemisinin Resistance Are Not Prevalent in Plasmodium falciparum Isolated from Ugandan Children

The emergence of resistance to artemisinin derivatives in Southeast Asia, manifested as delayed clearance of Plasmodium falciparum following treatment with artemisinins, is a major concern. Recently, the artemisinin resistance phenotype was attributed to mutations in portions of a P. falciparum gene (PF3D7_1343700) encoding kelch (K13) propeller domains, providing a molecular marker to monitor the spread of resistance. The P. falciparum cysteine protease falcipain-2 (FP2; PF3D7_1115700) has been shown to contribute to artemisinin action, as hemoglobin degradation is required for potent drug activity, and a stop mutation in the FP2 gene was identified in parasites selected for artemisinin resistance. Although delayed parasite clearance after artemisinin-based combination therapy (ACT) has not yet been noted in Uganda and ACTs remain highly efficacious, characterizing the diversity of these genes is important to assess the potential for resistance selection and to provide a baseline for future surveillance. We therefore sequenced the K13-propeller domain and FP2 gene in P. falciparum isolates from children previously treated with ACT in Uganda, including samples from 2006–7 (n = 49) and from 2010–12 (n = 175). Using 3D7 as the reference genome, we identified 5 non-synonymous polymorphisms in the K13-propeller domain (133 isolates) and 35 in FP2 (160 isolates); these did not include the polymorphisms recently associated with resistance after in vitro selection or identified in isolates from Asia. The prevalence of K13-propeller and FP2 polymorphisms did not increase over time, and was not associated with either time since prior receipt of an ACT or the persistence of parasites ≥2 days following treatment with an ACT. Thus, the K13-propeller and FP2 polymorphisms associated with artemisinin resistance are not prevalent in Uganda, and we did not see evidence for selection of polymorphisms in these genes.