Limited recombination events in merozoite surface protein-1 alleles of Plasmodium falciparum on islands

Little Polymorphism at the K13 Propeller Locus in Worldwide Plasmodium falciparum Populations Prior to the Introduction of Artemisinin Combination Therapies

The emergence and spread of artemisinin-resistant Plasmodium falciparum is of huge concern for the global effort toward malaria control and elimination. Artemisinin resistance, defined as a delayed time to parasite clearance following administration of artemisinin, is associated with mutations in the Pfkelch13 gene of resistant parasites. To date, as many as 60 nonsynonymous mutations have been identified in this gene, but whether these mutations have been selected by artemisinin usage or merely reflect natural polymorphism independent of selection is currently unknown. To clarify this, we sequenced the Pfkelch13 propeller domain in 581 isolates collected before (420 isolates) and after (161 isolates) the implementation of artemisinin combination therapies (ACTs), from various regions of endemicity worldwide. Nonsynonymous mutations were observed in 1% of parasites isolated prior to the introduction of ACTs. Frequencies of mutant isolates, nucleotide diversity, and haplotype diversity were significantly higher in the parasites isolated from populations exposed to artemisinin than in those from populations that had not been exposed to the drug. In the artemisinin-exposed population, a significant excess of dN compared to dS was observed, suggesting the presence of positive selection. In contrast, pairwise comparison of dN and dS and the McDonald and Kreitman test indicate that purifying selection acts on the Pfkelch13 propeller domain in populations not exposed to ACTs. These population genetic analyses reveal a low baseline of Pfkelch13 polymorphism, probably due to purifying selection in the absence of artemisinin selection. In contrast, various Pfkelch13 mutations have been selected under artemisinin pressure.

Genetic diversity of the Plasmodium falciparum apical membrane antigen I gene in parasite population from the China-Myanmar border area

To investigate the genetic diversity of the Plasmodium falciparum apical membrane antigen 1 (PfAMA1) gene in Southeast Asia, we determined PfAMA1 sequences from 135 field isolates collected from the China-Myanmar border area and compared them with 956 publically available PfAMA1 sequences from seven global P. falciparum populations. This analysis revealed high genetic diversity of PfAMA1 in global P. falciparum populations with a total of 229 haplotypes identified. The genetic diversity of PfAMA1 gene from the China-Myanmar border is not evenly distributed in the different domains of this gene. Sequence diversity in PfAMA1 from the China-Myanmar border is lower than that observed in Thai, African and Oceanian populations, but higher than that in the South American population. This appeared to correlate well with the levels of endemicity of different malaria-endemic regions, where hyperendemic regions favor genetic cross of the parasite isolates and generation of higher genetic diversity. Neutrality tests show significant departure from neutrality in the entire ectodomain and Domain I of PfAMA1 in the China-Myanmar border parasite population. We found evidence supporting a substantial continent-wise genetic structure among P. falciparum populations, with the highest genetic differentiation detected between the China-Myanmar border and the South American populations. Whereas no alleles were unique to a specific region, there were considerable geographical differences in major alleles and their frequencies, highlighting further necessity to include more PfAMA1 alleles in vaccine designs.

Plasmodium vivax and Plasmodium falciparum at the crossroads of exchange among islands in Vanuatu: implications for malaria elimination strategies

Understanding the transmission and movement of Plasmodium parasites is crucial for malaria elimination and prevention of resurgence. Located at the limit of malaria transmission in the Pacific, Vanuatu is an ideal candidate for elimination programs due to low endemicity and the isolated nature of its island setting. We analyzed the variation in the merozoite surface protein 1 (msp1) and the circumsporozoite protein (csp) of P. falciparum and P. vivax populations to examine the patterns of gene flow and population structures among seven sites on five islands in Vanuatu. Genetic diversity was in general higher in P. vivax than P. falciparum from the same site. In P. vivax, high genetic diversity was likely maintained by greater extent of gene flow among sites and among islands. Consistent with the different patterns of gene flow, the proportion of genetic variance found among islands was substantially higher in P. falciparum (28.81–31.23%) than in P. vivax (-0.53–3.99%). Our data suggest that the current island-by-island malaria elimination strategy in Vanuatu, while adequate for P. falciparum elimination, might need to be complemented with more centrally integrated measures to control P. vivax movement across islands.

Patterns and dynamics of genetic diversity in Plasmodium falciparum: what past human migrations tell us about malaria

Plasmodium falciparum is the main agent of malaria, one of the major human infectious diseases affecting millions of people worldwide. The genetic diversity of P. falciparum populations is an essential factor in the parasite's ability to adapt to changes in its environment, enabling the development of drug resistance and the evasion from the host immune system through antigenic variation. Therefore, characterizing these patterns and understanding the main drivers of the pathogen's genetic diversity can provide useful inputs for informing control strategies. In this paper, we review the pioneering work led by Professor Kazuyuki Tanabe on the genetic diversity of P. falciparum populations. In a first part, we recall basic results from population genetics for quantifying within-population genetic diversity, and discuss the main mechanisms driving this diversity. Then, we show how these approaches have been used for reconstructing the historical spread of malaria worldwide, and how current patterns of genetic diversity suggest that the pathogen followed our ancestors in their journey out of Africa. Because these results are robust to different types of genetic markers, they provide a baseline for predicting the pathogen's diversity in unsampled populations, and some useful elements for predicting vaccine efficacy and informing malaria control strategies.

Characteristic age distribution of Plasmodium vivax infections after malaria elimination on Aneityum Island, Vanuatu

Resurgence is a major concern after malaria elimination. After the initiation of the elimination program on Aneityum Island in 1991, microscopy showed that Plasmodium falciparum disappeared immediately, whereas P. vivax disappeared from 1996 onward, until P. vivax cases were reported in January 2002. By conducting malariometric surveys of the entire population of Aneityum, we investigated the age distribution of individuals with parasites during this epidemic in the context of antimalarial antibody levels and parasite antigen diversity. In July 2002, P. vivax infections were detected by microscopy in 22/759 individuals: 20/298 born after the beginning of the elimination program in 1991, 2/126 born between 1982 and 1991, and none of 335 born before 1982. PCR increased the number of infections detected to 77, distributed among all age groups. Prevalences were 12.1%, 16.7%, and 6.0%, respectively (P < 0.001). In November, a similar age pattern was found, but with fewer infections: 6/746 and 39/741 individuals were found to be infected by microscopy and PCR, respectively. The frequencies of antibody responses to P. vivax were significantly higher in individuals born before 1991 than in younger age groups and were similar to those on Malakula Island, an area of endemicity. Remarkably low antigen diversity (h, 0.15) of P. vivax infections was observed on Aneityum compared with the other islands (h, 0.89 to 1.0). A P. vivax resurgence was observed among children and teenagers on Aneityum, an age distribution similar to those before elimination and on islands where P. vivax is endemic, suggesting that in the absence of significant exposure, immunity may persist, limiting infection levels in adults. The limited parasite gene pool on islands may contribute to this protection.

Plasmodium falciparum: genetic diversity and complexity of infections in an isolated village in western Thailand

Genetic diversity of Plasmodium falciparum is intimately associated with morbidity, mortality and malaria control strategies. It is therefore imperative to study genetic makeup and population structure of this parasite in endemic areas. In Kong Mong Tha, an isolated village in western Thailand, the majority of P. falciparum infections are asymptomatic. In this study we investigated complexity of infections and single nucleotide polymorphisms (SNPs) in the P. falciparum population of Kong Mong Tha, and compared results with those previously obtained from Mae Sod, in northwestern Thailand, where the majority of infections were symptomatic. Using PCR-based determination of the 5' merozoite surface protein 1 gene (msp1) recombinant types, we found that 39% of 59 P. falciparum isolates from Kong Mong Tha had multiple 5' recombinant types with a mean number of 1.54. These values were much lower than those obtained from Mae Sod: 96% for multiple infections and with a mean number of 3.61. Analysis of full-length sequences of two housekeeping genes, the P-type Ca(2+)-transporting ATPase gene (n=33) plus adenylosuccinate lyase gene (n=33), and three vaccine candidate antigen genes, msp1 (n=26), the circumsporozoite protein gene, csp (n=30) and the apical membrane antigen 1 gene, ama 1 (n=32), revealed that in all of these genes within-population SNP diversity was at similar levels between Kong Mong Tha and Mae Sod, suggesting that the extent of MOI and clinical manifestations of malaria are not strongly associated with genetic diversity. Additionally, we did not detect significant genetic differentiation between the two parasite populations, as estimated by the Wright's fixation index of inter-population variance in allele frequencies, suggesting that gene flow prevented the formation of population structuring. Thus, this study highlights unique features of P. falciparum populations in Thailand. The implications of these finding are discussed.

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.

Travellers as sentinels: Assaying the worldwide distribution of polymorphisms associated with artemisinin combination therapy resistance in Plasmodium falciparum using malaria cases imported into Scotland

There is growing evidence that Plasmodium falciparum parasites in southeastern Asia have developed resistance to artemisinin combination therapy. The resistance phenotype has recently been shown to be associated with four single nucleotide polymorphisms in the parasite's genome. We assessed the prevalence of two of these single nucleotide polymorphisms in P. falciparum parasites imported into Scotland between 2009 and 2012, and in additional field samples from six countries in southeastern Asia. We analysed 28 samples from 11 African countries, and 25 samples from nine countries in Asia/southeastern Asia/Oceania. Single nucleotide polymorphisms associated with artemisinin combination therapy resistance were not observed outside Thailand and Cambodia.

Large-scale survey for novel genotypes of Plasmodium falciparum chloroquine-resistance gene pfcrt

Background

In Plasmodium falciparum, resistance to chloroquine (CQ) is conferred by a K to T mutation at amino acid position 76 (K76T) in the P. falciparum CQ transporter (PfCRT). To date, at least 15 pfcrt genotypes, which are represented by combinations of five amino acids at positions 72-76, have been described in field isolates from various endemic regions. To identify novel mutant pfcrt genotypes and to reveal the genetic relatedness of pfcrt genotypes, a large-scale survey over a wide geographic area was performed.

Methods

Sequences for exon 2 in pfcrt, including known polymorphic sites at amino acid positions 72, 74, 75 and 76, were obtained from 256 P. falciparum isolates collected from eight endemic countries in Asia (Bangladesh, Cambodia, Lao P.D.R., the Philippines and Thailand), Melanesia (Papua New Guinea and Vanuatu) and Africa (Ghana). A haplotype network was constructed based on six microsatellite markers located -29 kb to 24 kb from pfcrt in order to examine the genetic relatedness among mutant pfcrt genotypes.

Results

In addition to wild type (CVMNK at positions 72-76), four mutant pfcrt were identified; CVIET, CVIDT, SVMNT and CVMNT (mutated amino acids underlined). Haplotype network revealed that there were only three mutant pfcrt lineages, originating in Indochina, Philippines and Melanesia. Importantly, the Indochina lineage contained two mutant pfcrt genotypes, CVIET (n = 95) and CVIDT (n = 14), indicating that CVIDT shares a common origin with CVIET. Similarly, one major haplotype in the Melanesian lineage contained two pfcrt genotypes; SVMNT (n = 71) and CVMNT (n = 3). In Africa, all mutant pfcrt genotypes were the CVIET of the Indochina lineage, probably resulting from the intercontinental migration of CQ resistance from Southeast Asia.

Conclusions

The number of CQ-mutant lineages observed in this study was identical to that found in previous studies. This supports the hypothesis that the emergence of novel CQ resistance is rare. However, in the mutant pfcrt genotypes, amino acid changes at positions 72, 74 and 75 appear to have recently been generated at least several times, producing distinct pfcrt mutant genotypes. The occurrence of new mutations flanking K76T may yield stronger resistance to CQ and/or a higher fitness than the original pfcrt mutant.

Plasmodium falciparum accompanied the human expansion out of Africa

Plasmodium falciparum is distributed throughout the tropics and is responsible for an estimated 230 million cases of malaria every year, with a further 1.4 billion people at risk of infection. Little is known about the genetic makeup of P. falciparum populations, despite variation in genetic diversity being a key factor in morbidity, mortality, and the success of malaria control initiatives. Here we analyze a worldwide sample of 519 P. falciparum isolates sequenced for two housekeeping genes (63 single nucleotide polymorphisms from around 5000 nucleotides per isolate). We observe a strong negative correlation between within-population genetic diversity and geographic distance from sub-Saharan Africa (R(2) = 0.95) over Africa, Asia, and Oceania. In contrast, regional variation in transmission intensity seems to have had a negligible impact on the distribution of genetic diversity. The striking geographic patterns of isolation by distance observed in P. falciparum mirror the ones previously documented in humans and point to a joint sub-Saharan African origin between the parasite and its host. Age estimates for the expansion of P. falciparum further support that anatomically modern humans were infected prior to their exit out of Africa and carried the parasite along during their colonization of the world.

Active case detection, treatment of falciparum malaria with combined chloroquine and sulphadoxine/pyrimethamine and vivax malaria with chloroquine and molecular markers of anti-malarial resistance in the Republic of Vanuatu

Background

Chloroquine-resistant Plasmodium falciparum was first described in the Republic of Vanuatu in the early 1980s. In 1991, the Vanuatu Ministry of Health instituted new treatment guidelines for uncomplicated P. falciparum infection consisting of chloroquine/sulphadoxine-pyrimethamine combination therapy. Chloroquine remains the recommended treatment for Plasmodium vivax.

Methods

In 2005, cross-sectional blood surveys at 45 sites on Malo Island were conducted and 4,060 adults and children screened for malaria. Of those screened, 203 volunteer study subjects without malaria at the time of screening were followed for 13 weeks to observe peak seasonal incidence of infection. Another 54 subjects with malaria were followed over a 28-day period to determine efficacy of anti-malarial therapy; chloroquine alone for P. vivax and chloroquine/sulphadoxine-pyrimethamine for P. falciparum infections.

Results

The overall prevalence of parasitaemia by mass blood screening was 6%, equally divided between P. falciparum and P. vivax. Twenty percent and 23% of participants with patent P. vivax and P. falciparum parasitaemia, respectively, were febrile at the time of screening. In the incidence study cohort, after 2,303 person-weeks of follow-up, the incidence density of malaria was 1.3 cases per person-year with P. vivax predominating. Among individuals participating in the clinical trial, the 28-day chloroquine P. vivax cure rate was 100%. The 28-day chloroquine/sulphadoxine-pyrimethamine P. falciparum cure rate was 97%. The single treatment failure, confirmed by merozoite surface protein-2 genotyping, was classified as a day 28 late parasitological treatment failure. All P. falciparum isolates carried the Thr-76 pfcrt mutant allele and the double Asn-108 + Arg-59 dhfr mutant alleles. Dhps mutant alleles were not detected in the study sample.

Conclusion

Peak seasonal malaria prevalence on Malo Island reached hypoendemic levels during the study observation period. The only in vivo malaria drug efficacy trial thus far published from the Republic of Vanuatu showed chloroquine/sulphadoxine-pyrimethamine combination therapy for P. falciparum and chloroquine alone for P. vivax to be highly efficacious. Although the chloroquine-resistant pfcrt allele was present in all P. falciparum isolates, mutant alleles in the dhfr and dhps genes do not yet occur to the extent required to confer sulphadoxine-pyrimethamine resistance in this population.

Allelic diversity in the merozoite surface protein 1 gene of Plasmodium falciparum on Palawan Island, the Philippines

Allelic diversity of the Plasmodium falciparum merozoite surface protein 1 gene (msp1) is mainly generated by meiotic recombination at the mosquito stage. We investigated recombination-based allelic diversity of msp1 in a P. falciparum population from Palawan Island, the Philippines, where malaria transmission is moderate. We identified the 5' recombinant types, 3' sequence types and msp1 haplotypes (unique combinations of 5' recombinant type and 3' sequence type), and compared them with those of P. falciparum from the Solomon Islands, where malaria transmission is high. The mean number of 5' recombinant types per patient in Palawan was 1.44, which is comparable to the number for the Solomon Islands (1.41). The Palawan parasite population had 15 msp1 haplotypes, whereas the Solomon Islands population had only 8 haplotypes. The Palawan population showed strong linkage disequilibrium between polymorphic blocks/sites within msp1, which is comparable to the results for the Solomon Islands. These findings support our hypothesis that the extent of allelic diversity of msp1 is determined not only by the transmission intensity but also by the number of msp1 alleles prevalent in the local parasite population and the extent of mixed-allele infections. Contribution of a high prevalence of the chloroquine (CQ)-sensitive allele of P. falciparum CQ resistance transporter (pfcrt) to the relatively high msp1 diversity in the Palawan population is discussed.

Genetic diversity and gene flow of humans, Plasmodium falciparum, and Anopheles farauti s.s. of Vanuatu: inferred malaria dispersal and implications for malaria control

A comparison of the patterns of gene flow within and between islands and the genetic diversities of the three species required for malaria transmission (humans, Plasmodium falciparum, and Anopheles farauti s.s.) within the model island system of Vanuatu, shows that the active dispersal of An. farauti s.s. is responsible for within island movement of parasites. In contrast, since both P. falciparum and An. farauti s.s. populations are largely restricted to islands, movement of parasites between islands is likely due to human transport. Thus, control of vectors is crucial for controlling malaria within islands, while control of human movement is essential to control malaria transmission across the archipelago.

Molecular evolution and intragenic recombination of the merozoite surface protein MSP-3alpha from the malaria parasite Plasmodium vivax in Thailand

The merozoite surface antigens of malaria parasites are prime anti-morbidity/mortality vaccine candidates. However, their highly polymorphic nature requires extensive surveys of parasite populations to validate vaccine designs. Previous studies have found 3 molecular types (A, B and C) of the Plasmodium vivax merozoite surface protein 3a (PvMSP-3alpha) among parasite field populations. Here we analysed complete PvMSP-3alpha sequences from 17 clinical P. vivax isolates from Thailand and found that the nucleotide diversity was as high as that from samples widely separated by time and space. The polymorphic sites were not randomly distributed but concentrated in the N-terminal Ala-rich domain (block 2A), which is partially deleted in type B and C sequences. The size variations among type A sequences were due to small indels occurring in block 2A, whereas type B and C sequences were uniform in length with each type having a different large deletion. Analysis of synonymous and non-synonymous substitutions suggested that different selection forces were operating on different regions of the molecule. The numerous recombination sites detected within the Ala-rich domain suggested that intragenic recombination was at least partially responsible for the observed genetic diversity of the PvMSP-3alpha gene. Phylogenetic analysis failed to link any alleles to a specific geographical origin, even when different domains of PvMSP-3alpha were used for analysis. The highly polymorphic nature and lack of geographical clustering of isolates suggest that more systematic investigations of the PvMSP-3alpha gene are needed to explore its evolution and vaccine potential.

Relative frequencies of polymorphisms of variation in Block 2 repeats and 5' recombinant types of Plasmodium falciparum msp1 alleles

The mechanisms producing the genetic polymorphism at Plasmodium falciparum merozoite surface antigen-1 locus (pfmsp1) include the insertion and deletion of the different type of dimorphic Block 2 9-nucleotide repeat units as well as the intragenic recombination. To study relative occurrence frequencies of these two distinct mechanisms, we have developed a sensitive PCR strategy to identify both 5' recombinant types and the number of Block 2 repeats from the same sample. This method can specifically detect the target 5' recombinant type (Blocks 2-6) at the sensitivity of 1-4 copies of the pfmsp1. Applying the new method to field isolates from the Solomon Islands enabled us to identify six different 5' recombinant types and variation in Block 2 repeat number in three of them, thus distinguishing 10 different alleles. Distribution of these alleles in local three villages in the study area suggests that frequencies of variation in the number of Block 2 9-bp repeats and recombination events within Blocks 2-6 are mutually independent and the frequency of repeat variation is relatively high as compared to that of recombination events at the pfmsp1 locus in P. falciparum populations from the Solomon Islands.

Malaria dispersal among islands: human mediated Plasmodium falciparum gene flow in Vanuatu, Melanesia

A comparison of human and Plasmodium falciparum gene flow patterns in the model island system of Vanuatu, the limit of malaria in the Pacific reveals that human movement is essential for long, but not short distance P. falciparum gene flow. This suggests that long distance movement of humans may accelerate the evolution and spread of drug resistance and therefore exacerbate the global malaria problem.

Mosaic organization and heterogeneity in frequency of allelic recombination of the Plasmodium vivax merozoite surface protein-1 locus

The organization and allelic recombination of the merozoite surface protein-1 gene of Plasmodium vivax (PvMsp-1), the most widely prevalent human malaria parasite, were evaluated in complete nucleotide sequences of 40 isolates from various geographic areas. Alignment of 31 distinct alleles revealed the mosaic organization of PvMsp-1, consisting of seven interallele conserved blocks flanked by six variable blocks. The variable blocks showed extensive variation in repeats and nonrepeat unique sequences. Numerous recombination sites were distributed throughout PvMsp-1, in both conserved blocks and variable block unique sequences, and the distribution was not uniform. Heterozygosity of PvMsp-1 alleles was higher in Asia (0.953 +/- 0.009) than in Brazil (0.813 +/- 0.047). No identical alleles were shared between Asia and Brazil, whereas all but one variable block nonrepeat sequence found in Brazil occurred in Asia. These observations suggest that P. vivax populations in Asia are ancestral to Brazilian populations, and that PvMsp-1 has heterogeneity in frequency of allelic recombination events. Recurrent origins of new PvMsp-1 alleles by repeated recombination events were supported by a rapid decline in linkage disequilibrium between pairs of synonymous sites with increasing nucleotide distance, with little linkage disequilibrium at a distance of over 3 kb in a P. vivax population from Thailand, evidence for an effectively high recombination rate of the parasite. Meanwhile, highly reduced nucleotide diversity was noted in a region encoding the 19-kDa C-terminal epidermal growth factor-like domain of merozoite surface protein-1, a vaccine candidate.