Search

Background: Plasmodium vivax has been recently discovered as a significant cause of malaria in Mauritania, although very rare elsewhere in West Africa. It has not been known if this is a recently introduced or locally remnant parasite population, nor whether the genetic structure reflects epidemic or endemic transmission., Methodology/principal Findings: To investigate the P. vivax population genetic structure in Mauritania and compare with populations previously analysed elsewhere, multi-locus genotyping was undertaken on 100 clinical isolates, using a genome-wide panel of 38 single nucleotide polymorphisms (SNPs), plus seven SNPs in drug resistance genes. The Mauritanian P. vivax population is shown to be genetically diverse and divergent from populations elsewhere, indicated consistently by genetic distance matrix analysis, principal components analyses, and fixation indices. Only one isolate had a genotype clearly indicating recent importation, from a southeast Asian source. There was no linkage disequilibrium in the local parasite population, and only a small number of infections appeared to be closely genetically related, indicating that there is ongoing genetic recombination consistent with endemic transmission. The P. vivax diversity in a remote mining town was similar to that in the capital Nouakchott, with no indication of local substructure or of epidemic population structure. Drug resistance alleles were virtually absent in Mauritania, in contrast with P. vivax in other areas of the world., Conclusions/significance: The molecular epidemiology indicates that there is long-standing endemic transmission that will be very challenging to eliminate. The virtual absence of drug resistance alleles suggests that most infections have been untreated, and that this endemic infection has been more neglected in comparison to P. vivax elsewhere., Competing Interests: The authors have declared that no competing interests exist.

Plasmodium vivax is the most widespread and difficult to treat cause of human malaria. The development of vaccines against the blood stages of P. vivax remains a key objective for the control and elimination of vivax malaria. Erythrocyte binding-like (EBL) protein family members such as Duffy binding protein (PvDBP) are of critical importance to erythrocyte invasion and have been the major target for vivax malaria vaccine development. In this study, we focus on another member of EBL protein family, P. vivax erythrocyte binding protein (PvEBP). PvEBP was first identified in Cambodian (C127) field isolates and has subsequently been showed its preferences for binding reticulocytes which is directly inhibited by antibodies. We analysed PvEBP sequence from 316 vivax clinical isolates from eight countries including China (n = 4), Ethiopia (n = 24), Malaysia (n = 53), Myanmar (n = 10), Papua New Guinea (n = 16), Republic of Korea (n = 10), Thailand (n = 174), and Vietnam (n = 25). PvEBP gene exhibited four different phenotypic clusters based on the insertion/deletion (indels) variation. PvEBP-RII (179-479 aa.) showed highest polymorphism similar to other EBL family proteins in various Plasmodium species. Whereas even though PvEBP-RIII-V (480-690 aa.) was the most conserved domain, that showed strong neutral selection pressure for gene purifying with significant population expansion. Antigenicity of both of PvEBP-RII (16.1%) and PvEBP-RIII-V (21.5%) domains were comparatively lower than other P. vivax antigen which expected antigens associated with merozoite invasion. Total IgG recognition level of PvEBP-RII was stronger than PvEBP-RIII-V domain, whereas total IgG inducing level was stronger in PvEBP-RIII-V domain. These results suggest that PvEBP-RII is mainly recognized by natural IgG for innate protection, whereas PvEBP-RIII-V stimulates IgG production activity by B-cell for acquired immunity. Overall, the low antigenicity of both regions in patients with vivax malaria likely reflects genetic polymorphism for strong positive selection in PvEBP-RII and purifying selection in PvEBP-RIII-V domain. These observations pose challenging questions to the selection of EBP and point out the importance of immune pressure and polymorphism required for inclusion of PvEBP as a vaccine candidate., Competing Interests: The authors have declared that no competing interests exist.

The zoonotic Plasmodium knowlesi parasite is the most common cause of human malaria in Malaysia. Genetic analysis has shown that the parasites are divided into three subpopulations according to their geographic origin (Peninsular or Borneo) and, in Borneo, their macaque host (Macaca fascicularis or M. nemestrina). Whilst evidence suggests that genetic exchange events have occurred between the two Borneo subpopulations, the picture is unclear in less studied Peninsular strains. One difficulty is that P. knowlesi infected individuals tend to present with low parasitaemia leading to samples with insufficient DNA for whole genome sequencing. Here, using a parasite selective whole genome amplification approach on unprocessed blood samples, we were able to analyse recent genomes sourced from both Peninsular Malaysia and Borneo. The analysis provides evidence that recombination events are present in the Peninsular Malaysia parasite subpopulation, which have acquired fragments of the M. nemestrina associated subpopulation genotype, including the DBPβ and NBPXa erythrocyte invasion genes. The NBPXb invasion gene has also been exchanged within the macaque host-associated subpopulations of Malaysian Borneo. Our work provides strong evidence that exchange events are far more ubiquitous than expected and should be taken into consideration when studying the highly complex P. knowlesi population structure.

Background: Plasmodium vivax malaria remains a major public health burden in Myanmar. Resistance to chloroquine (CQ), the first-line treatment for P. vivax, has been reported in the country and has potential to undermine local control efforts., Methods: Patients over 6 years of age with uncomplicated P. vivax mono-infection were enrolled into clinical efficacy studies in Myawaddy in 2014 and Kawthoung in 2012. Study participants received a standard dose of CQ (25 mg/kg over 3 days) followed by weekly review until day 28. Pvmdr1 copy number (CN) and microsatellite diversity were assessed on samples from the patients enrolled in the clinical study and additional cross-sectional surveys undertaken in Myawaddy and Shwegyin in 2012., Results: A total of 85 patients were enrolled in the CQ clinical studies, 25 in Myawaddy and 60 in Kawthoung. One patient in Myawaddy (1.2%) had an early treatment failure and two patients (2.3%) in Kawthoung presented with late treatment failures on day 28. The day 28 efficacy was 92.0% (95% CI 71.6-97.9) in Myawaddy and 98.3% (95% CI 88.7-99.8) in Kawthoung. By day 2, 92.2% (23/25) in Myawaddy and 85.0% (51/60) in Kawthoung were aparasitaemic. Genotyping and pvmdr1 CN assessment was undertaken on 43, 52 and 46 clinical isolates from Myawaddy, Kawthoung and Shwegyin respectively. Pvmdr1 amplification was observed in 3.2% (1/31) of isolates in Myawaddy, 0% (0/49) in Kawthoung and 2.5% (1/40) in Shwegyin. Diversity was high in all sites (H E 0.855-0.876), with low inter-population differentiation (F ST 0.016-0.026, P < 0.05)., Conclusions: Treatment failures after chloroquine were observed following chloroquine monotherapy, with pvmdr1 amplification present in both Myawaddy and Shwegyin. The results emphasize the importance of ongoing P. vivax drug resistance surveillance in Myanmar, particularly given the potential connectivity between parasite population at different sites.

The widespread distribution and relapsing nature of Plasmodium vivax infection present major challenges for the elimination of malaria. To characterize the genetic diversity of this parasite in individual infections and across the population, we performed deep genome sequencing of >200 clinical samples collected across the Asia-Pacific region and analyzed data on >300,000 SNPs and nine regions of the genome with large copy number variations. Individual infections showed complex patterns of genetic structure, with variation not only in the number of dominant clones but also in their level of relatedness and inbreeding. At the population level, we observed strong signals of recent evolutionary selection both in known drug resistance genes and at new loci, and these varied markedly between geographical locations. These findings demonstrate a dynamic landscape of local evolutionary adaptation in the parasite population and provide a foundation for genomic surveillance to guide effective strategies for control and elimination of P. vivax., Competing Interests: The authors declare no competing financial interests.

Background: Vivax malaria was successfully eliminated from the Republic of Korea (ROK) in the late 1970s but re-emerged in 1993. Two decades later as the ROK enters the final stages of malaria elimination, dedicated surveillance of the local P. vivax population is critical. We apply a population genetic approach to gauge P. vivax transmission dynamics in the ROK between 2010 and 2012., Methodology/principal Findings: P. vivax positive blood samples from 98 autochthonous cases were collected from patients attending health centers in the ROK in 2010 (n = 27), 2011 (n = 48) and 2012 (n = 23). Parasite genotyping was undertaken at 9 tandem repeat markers. Although not reaching significance, a trend of increasing population diversity was observed from 2010 (HE = 0.50 ± 0.11) to 2011 (HE = 0.56 ± 0.08) and 2012 (HE = 0.60 ± 0.06). Conversely, linkage disequilibrium declined during the same period: IAS = 0.15 in 2010 (P = 0.010), 0.09 in 2011 (P = 0.010) and 0.05 in 2012 (P = 0.010). In combination with data from other ROK studies undertaken between 1994 and 2007, our results are consistent with increasing parasite divergence since re-emergence. Polyclonal infections were rare (3% infections) suggesting that local out-crossing alone was unlikely to explain the increased divergence. Cases introduced from an external reservoir may therefore have contributed to the increased diversity. Aside from one isolate, all infections carried a short MS20 allele (142 or 149 bp), not observed in other studies in tropical endemic countries despite high diversity, inferring that these regions are unlikely reservoirs., Conclusions: Whilst a number of factors may explain the observed population genetic trends, the available evidence suggests that an external geographic reservoir with moderate diversity sustains the majority of P. vivax infection in the ROK, with important implications for malaria elimination.

Background: P. vivax is an important public health burden in Ethiopia, accounting for almost half of all malaria cases. Owing to heterogeneous transmission across the country, a stronger evidence base on local transmission dynamics is needed to optimise allocation of resources and improve malaria interventions., Methodology and Principal Findings: In a pilot evaluation of local level P. vivax molecular surveillance in southern Ethiopia, the diversity and population structure of isolates collected between May and November 2013 were investigated. Blood samples were collected from microscopy positive P. vivax patients recruited to clinical and cross-sectional surveys from four sites: Arbaminch, Halaba, Badawacho and Hawassa. Parasite genotyping was undertaken at nine tandem repeat markers. Eight loci were successfully genotyped in 197 samples (between 36 and 59 per site). Heterogeneity was observed in parasite diversity and structure amongst the sites. Badawacho displayed evidence of unstable transmission, with clusters of identical clonal infections. Linkage disequilibrium in Badawacho was higher (IAS = 0.32, P = 0.010) than in the other populations (IAS range = 0.01-0.02) and declined markedly after adjusting for identical infections (IAS = 0.06, P = 0.010). Other than Badawacho (HE = 0.70), population diversity was equivalently high across the sites (HE = 0.83). Polyclonal infections were more frequent in Hawassa (67%) than the other populations (range: 8-44%). Despite the variable diversity, differentiation between the sites was low (FST range: 5 x 10-3-0.03)., Conclusions: Marked variation in parasite population structure likely reflects differing local transmission dynamics. Parasite genotyping in these heterogeneous settings has potential to provide important complementary information with which to optimise malaria control interventions.

Background: In Central China the declining incidence of Plasmodium vivax has been interrupted by epidemic expansions and imported cases. The impact of these changes on the local parasite population, and concurrent risks of future resurgence, was assessed., Methods: Plasmodium vivax isolates collected from Anhui and Jiangsu provinces, Central China between 2007 and 2010 were genotyped using capillary electrophoresis at seven polymorphic short tandem repeat markers. Spatial and temporal analyses of within-host and population diversity, population structure, and relatedness were conducted on these isolates., Results: Polyclonal infections were infrequent in the 94 isolates from Anhui (4%) and 25 from Jiangsu (12%), with a trend for increasing frequency from 2008 to 2010 (2 to 19%) when combined. Population diversity was high in both provinces and across the years tested (H(E) = 0.8 - 0.85). Differentiation between Anhui and Jiangsu was modest (F'(ST) = 0.1). Several clusters of isolates with identical multi-locus haplotypes were observed across both Anhui and Jiangsu. Linkage disequilibrium was strong in both populations and in each year tested (I(A)(S) = 0.2 - 0.4), but declined two- to four-fold when identical haplotypes were accounted for, indicative of occasional epidemic transmission dynamics. None of five imported isolates shared identical haplotypes to any of the central Chinese isolates., Conclusions: The population genetic structure of P. vivax in Central China highlights unstable transmission, with limited barriers to gene flow between the central provinces. Despite low endemicity, population diversity remained high, but the reservoirs sustaining this diversity remain unclear. The challenge of imported cases and risks of resurgence emphasize the need for continued surveillance to detect early warning signals. Although parasite genotyping has potential to inform the management of outbreaks, further studies are required to identify suitable marker panels for resolving local from imported P. vivax isolates.

Despite significant progress in the control of malaria in Malaysia, the complex transmission dynamics of P. vivax continue to challenge national efforts to achieve elimination. To assess the impact of ongoing interventions on P. vivax transmission dynamics in Sabah, we genotyped 9 short tandem repeat markers in a total of 97 isolates (8 recurrences) from across Sabah, with a focus on two districts, Kota Marudu (KM, n = 24) and Kota Kinabalu (KK, n = 21), over a 2 year period. STRUCTURE analysis on the Sabah-wide dataset demonstrated multiple sub-populations. Significant differentiation (F ST  = 0.243) was observed between KM and KK, located just 130 Km apart. Consistent with low endemic transmission, infection complexity was modest in both KM (mean MOI  = 1.38) and KK (mean MOI  = 1.19). However, population diversity remained moderate (H E  = 0.583 in KM and H E  = 0.667 in KK). Temporal trends revealed clonal expansions reflecting epidemic transmission dynamics. The haplotypes of these isolates declined in frequency over time, but persisted at low frequency throughout the study duration. A diverse array of low frequency isolates were detected in both KM and KK, some likely reflecting remnants of previous expansions. In accordance with clonal expansions, high levels of Linkage Disequilibrium (I A (S) >0.5 [P<0.0001] in KK and KM) declined sharply when identical haplotypes were represented once (I A (S)  = 0.07 [P = 0.0076] in KM, and I A (S) = -0.003 [P = 0.606] in KK). All 8 recurrences, likely to be relapses, were homologous to the prior infection. These recurrences may promote the persistence of parasite lineages, sustaining local diversity. In summary, Sabah's shrinking P. vivax population appears to have rendered this low endemic setting vulnerable to epidemic expansions. Migration may play an important role in the introduction of new parasite strains leading to epidemic expansions, with important implications for malaria elimination.

Summary: There is an immediate need for tools to both analyse and visualize in real-time single-nucleotide polymorphisms, insertions and deletions, and other structural variants from new sequence file formats. We have developed VarB software that can be used to visualize variant call format files in real time, as well as identify regions under balancing selection and informative markers to differentiate user-defined groups (e.g. populations). We demonstrate its utility using sequence data from 50 Plasmodium falciparum isolates comprising two different continents and confirm known signals from genomic regions that contain important antigenic and anti-malarial drug-resistance genes., Availability and Implementation: The C++-based software VarB and user manual are available from www.pathogenseq.org/varb., Contact: taane.clark@lshtm.ac.uk

The prevalence of CD36 deficiency in East Asian and African populations suggests that the causal variants are under selection by severe malaria. Previous analysis of data from the International HapMap Project indicated that a CD36 haplotype bearing a nonsense mutation (T1264G; rs3211938) had undergone recent positive selection in the Yoruba of Nigeria. To investigate the global distribution of this putative selection event, we genotyped T1264G in 3420 individuals from 66 populations. We confirmed the high frequency of 1264G in the Yoruba (26%). However, the 1264G allele is less common in other African populations and absent from all non-African populations without recent African admixture. Using long-range linkage disequilibrium, we studied two West African groups in depth. Evidence for recent positive selection at the locus was demonstrable in the Yoruba, although not in Gambians. We screened 70 variants from across CD36 for an association with severe malaria phenotypes, employing a case-control study of 1350 subjects and a family study of 1288 parent-offspring trios. No marker was significantly associated with severe malaria. We focused on T1264G, genotyping 10,922 samples from four African populations. The nonsense allele was not associated with severe malaria (pooled allelic odds ratio 1.0; 95% confidence interval 0.89-1.12; P = 0.98). These results suggest a range of possible explanations including the existence of alternative selection pressures on CD36, co-evolution between host and parasite or confounding caused by allelic heterogeneity of CD36 deficiency.

Evidence from autopsy and in vitro binding studies suggests that adhesion of erythrocytes infected with Plasmodium falciparum to the human host intercellular adhesion molecule (ICAM)-1 receptor is important in the pathogenesis of severe malaria. Previous association studies between polymorphisms in the ICAM1 gene and susceptibility to severe malarial phenotypes have been inconclusive and often contradictory. We performed genetic association studies with 15 single nucleotide polymorphisms (SNPs) around the ICAM1 locus. All SNPs were screened in a family study of 1071 trios from The Gambia, Malawi and Kenya. Two key non-synonymous SNPs with previously reported associations, rs5491 (K56M or 'ICAM-1(Kilifi)') and rs5498 (K469E), were tested in an additional 708 Gambian trios and a case-control study of 4058 individuals. None of the polymorphisms were associated with severe malaria phenotypes. Pooled results across our studies for ICAM-1(Kilifi) were, in severe malaria, odds ratio (OR) 1.02, 95% confidence interval (CI) 0.96-1.09, P=0.54, and cerebral malaria OR 1.07, CI 0.97-1.17, P=0.17. We assess the available epidemiological, population genetic and functional evidence that links ICAM-1(Kilifi) to severe malaria susceptibility.

Background: Localising regulatory variants that control gene expression is a challenge for genome research. Several studies have recently identified non-coding polymorphisms associated with inter-individual differences in gene expression. These approaches rely on the identification of signals of association against a background of variation due to other genetic and environmental factors. A complementary approach is to use an Allele-Specific Expression (ASE) assay, which is more robust to the effects of environmental variation and trans-acting genetic factors., Methodology/principal Findings: Here we apply an ASE method which utilises heterozygosity within an individual to compare expression of the two alleles of a gene in a single cell. We used individuals from three HapMap population groups and analysed the allelic expression of genes with cis-regulatory regions previously identified using total gene expression studies. We were able to replicate the results in five of the six genes tested, and refined the cis- associated regions to a small number of variants. We also showed that by using multi-populations it is possible to refine the associated cis-effect DNA regions., Conclusions/significance: We discuss the efficacy and drawbacks of both total gene expression and ASE approaches in the discovery of cis-acting variants. We show that the ASE approach has significant advantages as it is a cleaner representation of cis-acting effects. We also discuss the implication of using different populations to map cis-acting regions and the importance of finding regulatory variants which contribute to human phenotypic variation.