Effects of drug pressure and human migration on antimalarial resistance in circulating Plasmodium falciparum malaria parasites in Ecuador.

Antimalarial resistance in Plasmodium falciparum is a public health problem in the fight against malaria in Ecuador. Characterizing the molecular epidemiology of drug resistance genes helps to understand the emergence and spread of resistant parasites. In this study, the effects of drug pressure and human migration on antimalarial resistance in P. falciparum were evaluated. Sixty-seven samples from northwestern Ecuador from the 2019-2021 period were analyzed. SNPs in Pfcrt , Pfdhps , Pfdhfr , Pfmdr-1 , Pfk13 and Pfaat1 were identified by Sanger sequencing and whole-genome sequencing. A comparison of the frequencies of the haplotypes was made with data from the 2013-2015 period. Also, nucleotide and haplotype diversity were calculated. The frequencies of the mutant haplotypes, CVM ET in Pfcrt and C I C N I in Pfdhfr , increased. NED F S D F Y in Pfmdr-1 was detected for the first time. While the wild-type haplotypes, SAKAA in Pfdhps and MYRIC in Pfk13 , remained dominant. Interestingly, the A16 V mutation in Pfdhfr that gives resistance to proguanil is reported in Ecuador. In conclusion, parasites resistant to chloroquine ( Pfcrt ) and pyrimethamine ( Pfdhfr ) increased in recent years, while parasites sensitive to sulfadoxine ( Pfdhps ) and artemisinin ( Pfk13 ) prevail in Ecuador. Therefore, the current treatment is still useful against P. falciparum . The frequent human migration between Ecuador and Colombia has likely contributed to the spread of resistant parasites. Keys words : Plasmodium falciparum , resistance, antimalarial, selective pressure, human migration.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana.

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. In low transmission settings, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. We performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n = 1,409) through estimation of identity-by-descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 multi-isolate clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally observed 61 nonsynonymous substitutions that increased markedly in frequency over the study period as well as a novel pfk13 mutation (G718S). However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions, given that clones carrying drug resistance polymorphisms do not demonstrate enhanced persistence or higher abundance than clones carrying polymorphisms of comparable frequency that are unrelated to resistance. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana.

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. However, in low transmission settings where most mosquitoes become infected with only a single parasite clone, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. To investigate whether this clonality was potentially associated with the persistence and spatial spread of the mutation, we performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n=1,409) through estimation of identity by descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally report polymorphisms exhibiting evidence of selection for drug resistance or other phenotypes and reported a novel pfk13 mutation (G718S) as well as 61 nonsynonymous substitutions that increased markedly in frequency. However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Resolving drug selection and migration in an inbred South American Plasmodium falciparum population with identity-by-descent analysis.

The human malaria parasite Plasmodium falciparum is globally widespread, but its prevalence varies significantly between and even within countries. Most population genetic studies in P. falciparum focus on regions of high transmission where parasite populations are large and genetically diverse, such as sub-Saharan Africa. Understanding population dynamics in low transmission settings, however, is of particular importance as these are often where drug resistance first evolves. Here, we use the Pacific Coast of Colombia and Ecuador as a model for understanding the population structure and evolution of Plasmodium parasites in small populations harboring less genetic diversity. The combination of low transmission and a high proportion of monoclonal infections means there are few outcrossing events and clonal lineages persist for long periods of time. Yet despite this, the population is evolutionarily labile and has successfully adapted to changes in drug regime. Using newly sequenced whole genomes, we measure relatedness between 166 parasites, calculated as identity by descent (IBD), and find 17 distinct but highly related clonal lineages, six of which have persisted in the region for at least a decade. This inbred population structure is captured in more detail with IBD than with other common population structure analyses like PCA, ADMIXTURE, and distance-based trees. We additionally use patterns of intra-chromosomal IBD and an analysis of haplotypic variation to explore past selection events in the region. Two genes associated with chloroquine resistance, crt and aat1, show evidence of hard selective sweeps, while selection appears soft and/or incomplete at three other key resistance loci (dhps, mdr1, and dhfr). Overall, this work highlights the strength of IBD analyses for studying parasite population structure and resistance evolution in regions of low transmission, and emphasizes that drug resistance can evolve and spread in small populations, as will occur in any region nearing malaria elimination.

Plasmodium simium: Population Genomics Reveals the Origin of a Reverse Zoonosis.

BACKGROUND: The population history of Plasmodium simium, which causes malaria in sylvatic Neotropical monkeys and humans along the Atlantic Coast of Brazil, remains disputed. Genetically diverse P vivax populations from various sources, including the lineages that founded the species P simium, are thought to have arrived in the Americas in separate migratory waves. METHODS: We use population genomic approaches to investigate the origin and evolution of P simium. RESULTS: We find a minimal genome-level differentiation between P simium and present-day New World P vivax isolates, consistent with their common geographic origin and subsequent divergence on this continent. The meagre genetic diversity in P simium samples from humans and monkeys implies a recent transfer from humans to non-human primates - a unique example of malaria as a reverse zoonosis of public health significance. Likely genomic signatures of P simium adaptation to new hosts include the deletion of >40% of a key erythrocyte invasion ligand, PvRBP2a, which may have favored more efficient simian host cell infection. CONCLUSIONS: New World P vivax lineages that switched from humans to platyrrhine monkeys founded the P simium population that infects nonhuman primates and feeds sustained human malaria transmission in the outskirts of major cities.

Population genomics of Plasmodium vivax in Panama to assess the risk of case importation on malaria elimination.

Malaria incidence in Panama has plateaued in recent years in spite of elimination efforts, with almost all cases caused by Plasmodium vivax. Notwithstanding, overall malaria prevalence remains low (fewer than 1 case per 1000 persons). We used selective whole genome amplification to sequence 59 P. vivax samples from Panama. The P. vivax samples were collected from two periods (2007-2009 and 2017-2019) to study the population structure and transmission dynamics of the parasite. Imported cases resulting from increased levels of human migration could threaten malaria elimination prospects, and four of the samples evaluated came from individuals with travel history. We explored patterns of recent common ancestry among the samples and observed that a highly genetically related lineage (termed CL1) was dominant among the samples (47 out of 59 samples with good sequencing coverage), spanning the entire period of the collection (2007-2019) and all regions of the country. We also found a second, smaller clonal lineage (termed CL2) of four parasites collected between 2017 and 2019. To explore the regional context of Panamanian P. vivax we conducted principal components analysis and constructed a neighbor-joining tree using these samples and samples collected worldwide from a previous study. Three of the four samples with travel history clustered with samples collected from their suspected country of origin (consistent with importation), while one appears to have been a result of local transmission. The small number of Panamanian P. vivax samples not belonging to either CL1 or CL2 clustered with samples collected from Colombia, suggesting they represent the genetically similar ancestral P. vivax population in Panama or were recently imported from Colombia. The low diversity we observe in Panama indicates that this parasite population has been previously subject to a severe bottleneck and may be eligible for elimination. Additionally, while we confirmed that P. vivax is imported to Panama from diverse geographic locations, the lack of impact from imported cases on the overall parasite population genomic profile suggests that onward transmission from such cases is limited and that imported cases may not presently pose a major barrier to elimination.

Local emergence in Amazonia of Plasmodium falciparum k13 C580Y mutants associated with in vitro artemisinin resistance.

Antimalarial drug resistance has historically arisen through convergent de novo mutations in Plasmodium falciparum parasite populations in Southeast Asia and South America. For the past decade in Southeast Asia, artemisinins, the core component of first-line antimalarial therapies, have experienced delayed parasite clearance associated with several pfk13 mutations, primarily C580Y. We report that mutant pfk13 has emerged independently in Guyana, with genome analysis indicating an evolutionary origin distinct from Southeast Asia. Pfk13 C580Y parasites were observed in 1.6% (14/854) of samples collected in Guyana in 2016-2017. Introducing pfk13 C580Y or R539T mutations by gene editing into local parasites conferred high levels of in vitro artemisinin resistance. In vitro growth competition assays revealed a fitness cost associated with these pfk13 variants, potentially explaining why these resistance alleles have not increased in frequency more quickly in South America. These data place local malaria control efforts at risk in the Guiana Shield.

All recommended treatments against malaria include a drug called artemisinin or some of its derivatives. However, there are concerns that Plasmodium falciparum, the parasite that causes most cases of malaria, will eventually develop widespread resistance to the drug. A strain of P. falciparum partially resistant to artemisinin was seen in Cambodia in 2008, and it has since spread across Southeast Asia. The resistance appears to be frequently linked to a mutation known as pfk13 C580Y. Southeast Asia and Amazonia are considered to be hotspots for antimalarial drug resistance, and the pfk13 C580Y mutation was detected in the South American country of Guyana in 2010. To examine whether the mutation was still circulating in this part of the world, Mathieu et al. collected and analyzed 854 samples across Guyana between 2016 and 2017. Overall, 1.6% of the samples had the pfk13 C580Y mutation, but this number was as high as 8.8% in one region. Further analyses revealed that the mutation in Guyana had not spread from Southeast Asia, but that it had occurred in Amazonia independently. To better understand the impact of the pfk13 C580Y mutation, Mathieu et al. introduced this genetic change into non-resistant parasites from a country neighbouring Guyana. As expected, the mutation made P. falciparum highly resistant to artemisinin, but it also slowed the growth rate of the parasite. This disadvantage may explain why the mutation has not spread more rapidly through Guyana in recent years. Artemisinin and its derivatives are always associated with other antimalarial drugs to slow the development of resistance; there are concerns that reduced susceptibility to artemisinin leads to the parasites becoming resistant to the partner drugs. Further research is needed to evaluate how the pfk13 C580Y mutation affects the parasite's response to the typical combination of drugs that are given to patients.